CN110907923A - 基于平行因子算法的双基地emvs-mimo雷达角度估计算法及装置 - Google Patents

基于平行因子算法的双基地emvs-mimo雷达角度估计算法及装置 Download PDF

Info

Publication number
CN110907923A
CN110907923A CN201911037437.1A CN201911037437A CN110907923A CN 110907923 A CN110907923 A CN 110907923A CN 201911037437 A CN201911037437 A CN 201911037437A CN 110907923 A CN110907923 A CN 110907923A
Authority
CN
China
Prior art keywords
estimation
algorithm
matrix
factor
angle
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911037437.1A
Other languages
English (en)
Other versions
CN110907923B (zh
Inventor
黄江伟
文方青
王可
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yangtze University
Original Assignee
Yangtze University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Yangtze University filed Critical Yangtze University
Priority to CN201911037437.1A priority Critical patent/CN110907923B/zh
Publication of CN110907923A publication Critical patent/CN110907923A/zh
Application granted granted Critical
Publication of CN110907923B publication Critical patent/CN110907923B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/42—Diversity systems specially adapted for radar

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开一种基于平行因子算法的双基地EMVS‑MIMO雷达角度估计算法及装置，所述算法包括：根据匹配滤波后的阵列信号模型构建三阶平行因子张量模型；对所述三阶平行因子张量模型进行平行因子分解，得到因子矩阵的估计值；根据因子矩阵的估计值，采用ESPRIT技术获得对方位角的估计，利用矢量叉乘的性质获得俯仰角的估计；利用最小二乘技术获得二维极化角的估计。本发明所提算法可以获得低计算复杂度、高精度、自动配对的二维参数估计。

Description

基于平行因子算法的双基地EMVS-MIMO雷达角度估计算法及 装置
技术领域
本发明属于雷达定位技术领域，具体涉及一种基于平行因子分解的双基地MIMO雷达中的二维参数估计算法及装置。
背景技术
多输入多输出(Multiple-input Multiple-output，MIMO)技术是下一代移动通信的而核心技术，同时它也是未来雷达发展的主流方向。MIMO雷达采用多根发射天线发射正交波形，在接收端采用匹配滤波器分离出多根接收天线接收的信号，从而在发射-接收端之间形成多个虚拟通道。利用分集的思想，MIMO雷达可大大改善目标探测性能。其在在分辨率、抗衰落性、可辨识性以及抑制噪声等方面具有潜在的优势，因而引起国内外学术界和工程界的广泛关注。按照MIMO雷达收发阵元的分布不同可以将其分为分布式MIMO雷达和共址MIMO雷达两大类。分布式MIMO雷达采用广泛分布的收发阵元配置形式，这种体制的MIMO雷达能够克服RCS闪烁；共址MIMO雷达其收发阵元间距均较近，其可获得目标高分辨率的角度估计。本发明重点关注双基地MIMO雷达，其属于共址MIMO雷达的重要类型之一。
联合波离角(Direction-of-Departure，DOD)和波达角(Direction-of-Arrival，DOA)估计是双基地MIMO雷达目标定位的重要任务之一，也是目前MIMO雷达研究的热点问题之一。至今已涌现大量优秀的估计算法，典型的代表主要有Capon算法、多重谱峰分类(Multiple Signal Classification，MUSIC)算法、基于旋转不变技术的参数估计(Estimation Method of Signal Parameters via Rotational，ESPRIT)算法、传播算子(Propagator Method，PM)算法、最大似然法、矩阵束方法、高阶子空间分解(Higher OrderSingular Value Decomposition，HOSVD)算法、平行因子(Parallel Factor，PARAFAC)算法、基于稀疏表示的估计算法等。
然而，上述估计方法只能有效地获得一维波离角(Direction-of-Departure，DOD)和一维波达角(Direction-of-Arrival，DOA)。只有少数文献研究了2D-DOD和2D-DOA估计问题，比如Chen和Zhang的算法(C.Chen and X.Zhang,“A low-complexity joint 2D-DODand 2D-DOA estimation algorithm for MIMO radar with arbitrary arrays,”IntJ.Electron.,vol.100,no.10,pp.1455–1469,Jan.2013.)、Li和Zhang的算法(J.Li andX.Zhang,“Closed-form blind 2D-DOD and 2D-DOA estimation for MIMO radar witharbitrary arrays,”Wireless Per.Commun.,vol.69,no.1,pp.175–186,Mar.2013.)、Xia的算法(T.Xia,“Joint diagonalization based 2D-DOD and 2D-DOA estimation forbistatic MIMO radar,”Signal Process.,vol.116,pp.7–12,Nov.2015.)。上述估计方法的一个共同特点是使用非线性标量传感器阵列，例如：L形阵列、矩形阵列和任意阵列。与标量传感器不同，电磁矢量传感器(Electromagnetic VectorSensor,EMVS)可以测量二维(2D)方位角、俯仰角。此外，它还有几个固有的优势，例如：具有更好的识别能力、可以实现二维角度的自动配对、可提供信号源的附加极化信息(K.T.WongandX.Yuan,“vectorcross-productdirection-finding with an electromagnetic vector-sensor of sixorthogonally oriented but spatially noncollocating dipoles/loops,”IEEETrans.Signal Process.,vol.59,no.1,pp.160–171,Jan.2011.)。因此，Chintagunta等人提出了一种双基地电磁矢量传感器—多输入多输出(EMVS-MIMO)雷达系统，它在发射端和接收端均使用电磁矢量传感器阵列(S.Chintagunta and P.Palanisamy,“2D-DOD and 2D-DOA estimation using the electromagnetic vector sensors,”Signal Process.,vol.147,pp.163–172,Jun.2018.)。同时，该文献提出了一种改进的基于旋转不变技术的参数估计(Estimation method of signal parameters via rotational invarianttechnique,ESPRIT)算法，记为ESPRIT算法，用于测量2D-DOD、2D-DOA、二维传输极化角(Transmit-Polarization-Arrival,TPA)和二维接收极化角(Reception-Polarization-Arrival,RPA)。然而，该算法存在一系列缺陷。首先，该算法涉及到特征分解，计算复杂度很高。其次，多维样本间的张量结构被忽略了，算法的精度有待提高。此外，该方法需要额外配对所估计的2D-DOD和2D-DOA。
发明内容
针对上述缺陷，本申请提出了一种基于平行因子算法的双基地EMVS-MIMO雷达角度估计算法及装置，基于平行因子分解和张量结构解决了双基地MIMO雷达中的二维参数估计问题。
本发明第一方面，提出一种基于平行因子算法的双基地EMVS-MIMO雷达角度估计算法，包括：
S1、根据匹配滤波后的阵列信号模型构建三阶平行因子张量模型；
S2、对所述三阶平行因子张量模型进行平行因子分解，得到因子矩阵的估计值；
S3、根据因子矩阵的估计值，采用ESPRIT算法获得对方位角的估计，利用矢量叉乘的性质获得俯仰角的估计；
S4、利用最小二乘法获得二维极化角的估计。
优选的，所述步骤S1中，所述三阶平行因子张量模型具体的张量形式如下：
Figure BDA0002251916010000031
其中，
Figure BDA0002251916010000032
式中
Figure BDA0002251916010000033
Figure BDA0002251916010000034
bt,k是第k个目标的发射方向矢量；at,k是相应的发射极化响应矢量，br,k是第k个目标的接收方向矢量；ar,k是相应的接收极化响应矢量，F是目标特征矩阵，L为接收快拍数，
Figure BDA0002251916010000035
是相对应的噪声张量，下标×n表述张量的模n乘积。
优选的，所述步骤S2具体为：
S21、将对因子矩阵Dt，Dr，F的估计问题转化为求解如下优化问题：
Figure BDA0002251916010000041
S22、根据张量展开的定义，将
Figure BDA0002251916010000042
以矩阵的形式展开：
Figure BDA0002251916010000043
Figure BDA0002251916010000044
Figure BDA0002251916010000045
其中，
Figure BDA0002251916010000046
表示
Figure BDA0002251916010000047
的模n展开；将对因子矩阵Dt，Dr，F的估计问题转化为如下优化问题：
Figure BDA0002251916010000048
Figure BDA0002251916010000049
Figure BDA00022519160100000410
S23、通过三线性交替最小二乘法来求解步骤S22所述优化问题，其最优解分别为：
Figure BDA00022519160100000411
Figure BDA00022519160100000412
Figure BDA00022519160100000413
根据所述最优解得到因子矩阵的估计值
Figure BDA00022519160100000414
和
Figure BDA00022519160100000415
优选的，所述步骤S3中，2D-DOD估计过程具体为：
S301、计算
Figure BDA00022519160100000416
Figure BDA00022519160100000417
其中，
Figure BDA00022519160100000418
定义
Figure BDA00022519160100000419
Figure BDA00022519160100000420
表示置换矩阵，Δ1是K×K实对角矩阵；
S302、对
Figure BDA00022519160100000421
进行特征值分解，得到特征值
Figure BDA00022519160100000422
和相应的特征向量
Figure BDA00022519160100000423
θt,k的估计为：
Figure BDA0002251916010000051
S303、计算At的估计：
Figure BDA0002251916010000052
Figure BDA0002251916010000053
round{·}表示取近似值；令
Figure BDA0002251916010000054
和
Figure BDA00022519160100000518
分别代表
Figure BDA00022519160100000519
第k列的前三个和最后三个元素，则矢量叉积
Figure BDA0002251916010000057
φt,k的估计为：
Figure BDA0002251916010000058
其中，pt,k(1)和pt,k(2)分别表示pt,k的第一和第二个元素。
优选的，所述步骤S4中，2D-TPA估计为：
令
Figure BDA0002251916010000059
根据
Figure BDA00022519160100000510
构建方向矩阵Vt,k，则极化向量可以通过下式估计：
Figure BDA00022519160100000511
二维发射极化角的估计为：
Figure BDA00022519160100000512
Figure BDA00022519160100000513
其中gt,k(1)和gt,k(2)分别表示gt,k的第一和第二个元素。
优选的，所述步骤S3中，2D-DOA估计为：
S311、计算
Figure BDA00022519160100000514
Figure BDA00022519160100000515
其中
Figure BDA00022519160100000516
Figure BDA00022519160100000517
Figure BDA0002251916010000061
Δ2是K×K实对角矩阵；
S312、对
Figure BDA0002251916010000062
进行特征分解，λr,k为
Figure BDA0002251916010000063
的第k个对角元素，则θr,k的估计为：
Figure BDA0002251916010000064
S313、计算At的估计：
Figure BDA0002251916010000065
令
Figure BDA0002251916010000066
和
Figure BDA0002251916010000067
分别代表
Figure BDA0002251916010000068
第k列的前三个和最后三个元素，则矢量叉积
Figure BDA0002251916010000069
φt,k的估计为：
Figure BDA00022519160100000610
其中，pr,k(1)和pr,k(2)分别表示pr,k的第一和第二个元素。
优选的，所述步骤S4中，2D-RPA估计为：
令
Figure BDA00022519160100000611
根据
Figure BDA00022519160100000612
构建方向矩阵Vr,k，则极化向量可以通过下式估计：
Figure BDA00022519160100000613
二维发射极化角可以通过下式来计算：
Figure BDA00022519160100000614
Figure BDA00022519160100000615
其中gr,k(1)和gr,k(2)分别表示gr,k的第一和第二个元素。
本发明第二方面，提出一种基于平行因子算法的双基地EMVS-MIMO雷达角度估计装置，所述装置包括：
模型构建模块：根据匹配滤波后的阵列信号模型构建三阶平行因子张量模型；
因子分解模块：对所述三阶平行因子张量模型进行平行因子分解，得到因子矩阵的估计值；
角度估计模块：根据因子矩阵的估计值，采用ESPRIT技术获得对方位角的估计，利用矢量叉乘的性质获得俯仰角的估计；利用最小二乘技术获得二维极化角的估计。
本发明首先将双基地EMVS-MIMO雷达匹配滤波后的阵列信号表示成一个三阶平行因子张量模型，然后利用平行因子分解获得对因子矩阵的估计，再利用ESPRIT技术获得对方位角的估计，利用矢量叉乘的性质获得俯仰角的估计；最后利用最小二乘技术获得二维极化角的估计。本发明提出的算法及装置可以获得计算复杂度低、精度高、自动配对的二维参数估计。
附图说明
为了更清楚地说明本发明的技术方案，下面将对本发明技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本发明提供的基于平行因子算法的双基地EMVS-MIMO雷达角度估计算法流程示意图；
图2是本发明所提算法与对比算法在不同SNR下的RMSE比较；
图3是本发明所提算法与对比算法在不同N下的RMSE比较；
图4是本发明所提算法与对比算法在不同N下的平均运算时间的比较。
具体实施方式
为使得本发明的发明目的、特征、优点能够更加的明显和易懂，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，下面所描述的实施例仅仅是本发明一部分实施例，而非全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
请参阅图1，本发明提供的基于平行因子算法的双基地EMVS-MIMO雷达角度估计算法流程示意图。所述算法包括：
S1、根据匹配滤波后的阵列信号模型构建三阶平行因子张量模型；
S2、对所述三阶平行因子张量模型进行平行因子分解，得到因子矩阵的估计值；
S3、根据因子矩阵的估计值，采用ESPRIT算法获得对方位角的估计，利用矢量叉乘的性质获得俯仰角的估计；
S4、利用最小二乘法获得二维极化角的估计。
具体的，考虑一个双基地EMVS-MIMO雷达，其由M个EMVS发射阵元和N个EMVS接收真有组成，发射阵列和接收阵列都是具有λ/2间距的均匀线性阵列(ULA)，其中λ是发射波形的载波波长。假设在同一个距离单元中出现了K个远场目标。其匹配滤波后的阵列信号模型可以表示为：
Figure BDA0002251916010000081
其中，
Figure BDA0002251916010000082
是对应于第k个目标的发射方向矢量；
Figure BDA0002251916010000083
是相应的发射极化响应矢量，(θt,k,φt,k)和(γt,k,ηt,k)分别对应于第k个目标2D-DOD和的2D-TPA；
Figure BDA0002251916010000084
是对应于第k个目标的接收方向矢量；
Figure BDA0002251916010000085
是相应的接收极化响应矢量，(θr,k,φr,k)和(γr,k,ηr,k)分别对应于第k个目标2D-DOA和的2D-RPA；
Figure BDA0002251916010000086
是噪声样本；表达式1中
Figure BDA0002251916010000087
Figure BDA0002251916010000088
Bt和At分别表示发射阵列对应的发射方向矩阵和发射极化响应矩阵，Br和Ar分别表示接收阵列对应的接收方向矩阵和接收极化响应矩阵；F是目标特征矩阵，L为接收快拍数，
Figure BDA0002251916010000091
C为虚拟方向矩阵。
首先，我们构造一个三阶张量
Figure BDA0002251916010000092
其具体的张量形式如下
Figure BDA0002251916010000093
其中
Figure BDA00022519160100000917
是相对应的噪声张量，下标×n表述张量的模n乘积。因子矩阵Dt，Dr，F可以通过求解如下优化问题来估计
Figure BDA0002251916010000094
这里可以通过三线性交替最小二乘(TALS)法来解决。根据张量展开的定义，
Figure BDA0002251916010000095
可以以矩阵的形式展开为：
Figure BDA0002251916010000096
Figure BDA0002251916010000097
Figure BDA0002251916010000098
其中，
Figure BDA0002251916010000099
表示
Figure BDA00022519160100000910
的模n展开。因此，表达式3同时转化为如下优化问题：
Figure BDA00022519160100000911
Figure BDA00022519160100000912
Figure BDA00022519160100000913
上述问题转化成最小二乘法问题，可以用三线性交替最小二乘法来求解，其最优解分别为
Figure BDA00022519160100000914
Figure BDA00022519160100000915
Figure BDA00022519160100000916
在这里，我们使用COMFAC算法来加速迭代收敛，在COMFAC算法中，高维平行因子模型首先被压缩成较低维的模型，之后在压缩空间使用TALS算法，然后就再将解恢复到原始空间。
然后根据因子矩阵Dt，Dr，F的估计值分别进行方位角、俯仰角、极化角的估计：
1)2D-DOD和2D-TPA估计
众所周知，矩阵分解不是唯一的，除非矩阵满足一些严格的条件(例如矩阵正交)。然而，张量秩的分解往往是唯一的。如下定理(T.G.Kolda and B.W.Bader,“Tensordecompositions and applications,”SIAM Rev.,vol.51,no.3,pp.455–500,Aug.2009.)提供了一个平行因子分解唯一性的充分条件。根据该定理，如果
kDt+kDr+kF≥2K+3 表达式7
那么因子矩阵Dt，Dr，F的估计值对于列的置换和缩放是唯一的。置换不确定性意味着系数矩阵的列被重新排序，缩放指的是列通过乘以常数从而缩放。具体地，因子矩阵的估计值
Figure BDA0002251916010000101
和
Figure BDA0002251916010000102
可表示为：
Figure BDA0002251916010000103
Figure BDA0002251916010000104
Figure BDA0002251916010000105
其中
Figure BDA0002251916010000106
表示置换矩阵，Δ1、Δ2、Δ3是K×K实对角矩阵，其对角线元素为缩放因子，且Δ1Δ2Δ3＝IK。N1、N2和N3表示拟合误差。
定义
Figure BDA0002251916010000107
容易得到：
JM,2Bt＝JM,1Btψt 表达式9
在这里
Figure BDA0002251916010000108
令
Figure BDA0002251916010000109
我们可以得到：
Figure BDA00022519160100001011
忽略表达式8a中的噪声，可得：
Figure BDA00022519160100001010
对于置换矩阵∏和尺度矩阵Δ，易证：ΠΔ1＝Δ1Π-1，将其带入表达式11中得到：
Figure BDA0002251916010000111
对
Figure BDA0002251916010000112
使用特征值分解，可以得到特征值
Figure BDA0002251916010000113
和相应的特征向量
Figure BDA0002251916010000114
那么θt,k可以通过下式估计
Figure BDA0002251916010000115
Π可以通过下式近似
Figure BDA0002251916010000116
其中round{·}表示取近似值。因此，At可以通过下式估计
Figure BDA0002251916010000117
让
Figure BDA0002251916010000118
和
Figure BDA0002251916010000119
分别代表
Figure BDA00022519160100001110
第k列的前三个和最后三个元素。然后，可以计算下面的矢量叉积
Figure BDA00022519160100001111
φt,k可以通过下式获得
Figure BDA00022519160100001112
pt,k(1)和pt,k(2)分别表示pt,k的第一和第二个元素。
获得2D-DOD估计后，进行2D-TPA估计，令
Figure BDA00022519160100001113
一旦得到
Figure BDA00022519160100001114
我们可以构建方向矩阵Vt,k，则极化向量可以通过下式估计
Figure BDA00022519160100001115
二维发射极化角可以通过下式来计算
Figure BDA00022519160100001116
Figure BDA00022519160100001117
显然
Figure BDA0002251916010000121
可以自动匹配。
2)2D-DOA和2D-RPA估计
根据表达式8,
Figure BDA0002251916010000122
和
Figure BDA0002251916010000123
共享相同的置换矩阵，因此我们计算
Figure BDA0002251916010000124
类比于表达式10，我们可以得到：
Figure BDA00022519160100001220
其中，
Figure BDA0002251916010000125
Figure BDA0002251916010000126
之后我们可以得到
Figure BDA0002251916010000127
显然
Figure BDA0002251916010000128
让λr,k成为
Figure BDA0002251916010000129
的第k个对角元素，则θr,k可以通过：
Figure BDA00022519160100001221
来估算。因此，Ar可以通过下式估计
Figure BDA00022519160100001210
来估算，其中，
Figure BDA00022519160100001211
表示估计的ψr。类似于表达式15到表达式18，我们可以相应的得到
Figure BDA00022519160100001212
其具体过程为，令
Figure BDA00022519160100001213
和
Figure BDA00022519160100001214
分别代表
Figure BDA00022519160100001215
第k列的前三个和最后三个元素，则矢量叉积
Figure BDA00022519160100001216
φt,k的估计为：
Figure BDA00022519160100001217
其中，pr,k(1)和pr,k(2)分别表示pr,k的第一和第二个元素。
获得2D-DOA的估计后，进行2D-RPA估计，令
Figure BDA00022519160100001218
根据
Figure BDA00022519160100001219
构建方向矩阵Vr,k，则极化向量可以通过下式估计：
Figure BDA0002251916010000131
二维发射极化角可以通过下式来计算：
Figure BDA0002251916010000132
Figure BDA0002251916010000133
其中gr,k(1)和gr,k(2)分别表示gr,k的第一和第二个元素。
与上述算法相对应，本发明还提供一种基于平行因子算法的双基地EMVS-MIMO雷达角度估计装置，所述装置包括：
模型构建模块：根据匹配滤波后的阵列信号模型构建三阶平行因子张量模型；
因子分解模块：对所述三阶平行因子张量模型进行平行因子分解，得到因子矩阵的估计值；
角度估计模块：根据因子矩阵的估计值，采用ESPRIT技术获得对方位角的估计，利用矢量叉乘的性质获得俯仰角的估计；利用最小二乘技术获得二维极化角的估计。
以上装置实施例与算法实施例是一一对应的，装置实施例简略之处，参见算法实施例即可。
针对本发明提出的双基地EMVS-MIMO雷达角度估计算法，进行了大量的仿真实验。仿真中假设K＝3个目标处于远场，其2D-DOA、2D-DOD、2D-TPA和2D-RPA分别为θt＝(40°,20°),，30φt＝(15°,25°,35°)，γt＝(10°,22°,35°)，ηt＝(36°,48°,56°)，θr＝(24°,38°,16°)，φt＝(21°,32°,55°)，γt＝(42°,33°,60°)和ηt＝(17°,27°,39°)。仿真实验中发射EMVS的个数M＝6，接收EMVS的个数N，二者均为均匀线性阵列，阵元间距均为λ/2。仿真中接收快拍数设置为L＝200，信噪比(signal-to-noise ratio，SNR)定义为表达式(1)中信号与噪声功率的比值。为比较本发明所提算法同现有ESPRIT-Like算法估计精度的比较，对算法进行500次蒙特卡洛仿真，角度估计的精度用均方根误差(Root Mean Squared Error，RMSE)来评价，为简化显示的结果，我们只显示了方位角的RMSE均值和极化角的RMSE均值，分别在图中用后缀‘-d’和后缀‘-p’标注。本发明所提算法用‘Proposed’标注。
图2给出了本发明所提算法和ESPRIT算法及参数估计的克拉美-罗界(CRB)在不同SNR下性能比较的结果，其中N＝8。由仿真结果可知，所有算法的估计精度都会随着SNR的增大而改善。由于本发明利用了多维数据的结构信息，所提算法的参数估计精度明显优于现有ESPRIT算法。
图3与图4分别给出了在不同接收EMVS数目N下角度估计的精度与算法平均运算时间(Average running time，ART)的比较，其中SNR设置为0dB。可以明显的看到，所有算法随着N的增大，精度会有缓慢的改善，但是这种改善并不明显。同时，本发明所提算法的精度明显优于相对比的ESPRIT算法。此外，本发明所提的算法的运算复杂度明显低于ESPRIT算法，这种优势在N较大时尤其明显。
本发明提出一种基于平行因子分解的双基地MIMO雷达中的二维参数估计算法和装置。该算法首先将双基地MIMO雷达匹配滤波后的阵列信号表示成一个三阶平行因子张量模型，然后利用平行因子分解获得对因子矩阵的估计，再利用ESPRIT技术获得对方位角的估计。接着，利用矢量叉乘的性质获得俯仰角的估计。最后，利用最小二乘技术获得二维极化角的估计。本发明所提算法可以获得低计算复杂度、高精度、自动配对的二维参数估计。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述，在此不再赘述。
以上所述，以上实施例仅用以说明本发明的技术方案，而非对其限制，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (8)

1.一种基于平行因子算法的双基地EMVS-MIMO雷达角度估计算法，其特征在于，所述算法包括：
S1、根据匹配滤波后的阵列信号模型构建三阶平行因子张量模型；
S2、对所述三阶平行因子张量模型进行平行因子分解，得到因子矩阵的估计值；
S3、根据因子矩阵的估计值，采用ESPRIT算法获得对方位角的估计，利用矢量叉乘的性质获得俯仰角的估计；
S4、利用最小二乘法获得二维极化角的估计。
2.根据权利要求1所述基于平行因子算法的双基地EMVS-MIMO雷达角度估计算法，其特征在于，所述步骤S1中，所述三阶平行因子张量模型具体的张量形式如下：
Figure FDA0002251914000000011
其中，
Figure FDA0002251914000000012
式中
Figure FDA0002251914000000013
Figure FDA0002251914000000014
Bt和At分别表示发射阵列对应的发射方向矩阵和发射极化响应矩阵，Br和Ar分别表示接收阵列对应的接收方向矩阵和接收极化响应矩阵，⊙表示按列克罗内克积，F是目标特征矩阵，
Figure FDA0002251914000000015
是相对应的噪声张量，下标×n表述张量的模n乘积。
3.根据权利要求1所述基于平行因子算法的双基地EMVS-MIMO雷达角度估计算法，其特征在于，所述步骤S2具体为：
S21、将对因子矩阵Dt，Dr，F的估计问题转化为求解如下优化问题：
Figure FDA0002251914000000016
S22、根据张量展开的定义，将
Figure FDA0002251914000000017
以矩阵的形式展开：
Figure FDA0002251914000000018
Figure FDA0002251914000000019
Figure FDA0002251914000000021
其中，
Figure FDA0002251914000000022
表示
Figure FDA0002251914000000023
的模n展开；将对因子矩阵Dt，Dr，F的估计问题转化为如下优化问题：
Figure FDA0002251914000000024
Figure FDA0002251914000000025
Figure FDA0002251914000000026
S23、通过三线性交替最小二乘法来求解步骤S22所述优化问题，其最优解分别为：
Figure FDA0002251914000000027
Figure FDA0002251914000000028
Figure FDA0002251914000000029
根据所述最优解得到因子矩阵的估计值
Figure FDA00022519140000000210
和
Figure FDA00022519140000000211
4.根据权利要求1所述基于平行因子算法的双基地EMVS-MIMO雷达角度估计算法，其特征在于，所述步骤S3中，2D-DOD估计过程具体为：
S301、计算
Figure FDA00022519140000000212
Figure FDA00022519140000000213
其中，
Figure FDA00022519140000000214
定义
Figure FDA00022519140000000215
Figure FDA00022519140000000216
Figure FDA00022519140000000217
表示置换矩阵，Δ1是K×K实对角矩阵；
S302、对
Figure FDA00022519140000000218
进行特征值分解，得到特征值
Figure FDA00022519140000000219
和相应的特征向量
Figure FDA00022519140000000220
θt,k的估计为：
Figure FDA00022519140000000221
S303、计算At的估计：
Figure FDA00022519140000000222
Figure FDA0002251914000000031
round{·}表示取近似值；令
Figure FDA0002251914000000032
和
Figure FDA0002251914000000033
分别代表
Figure FDA0002251914000000034
第k列的前三个和最后三个元素，则矢量叉积
Figure FDA0002251914000000035
φt,k的估计为：
Figure FDA0002251914000000036
其中，pt,k(1)和pt,k(2)分别表示pt,k的第一和第二个元素。
5.根据权利要求4所述基于平行因子算法的双基地EMVS-MIMO雷达角度估计算法，其特征在于，所述步骤S4中，2D-TPA估计为：
令
Figure FDA0002251914000000037
根据
Figure FDA0002251914000000038
构建方向矩阵Vt,k，则极化向量可以通过下式估计：
Figure FDA0002251914000000039
二维发射极化角的估计为：
Figure FDA00022519140000000310
Figure FDA00022519140000000311
其中gt,k(1)和gt,k(2)分别表示gt,k的第一和第二个元素。
6.根据权利要求1所述基于平行因子算法的双基地EMVS-MIMO雷达角度估计算法，其特征在于，所述步骤S3中，2D-DOA估计为：
S311、计算
Figure FDA00022519140000000312
Figure FDA00022519140000000313
其中
Figure FDA00022519140000000314
Figure FDA00022519140000000315
Δ2是K×K实对角矩阵；
S312、对
Figure FDA00022519140000000316
进行特征分解，λr,k为
Figure FDA00022519140000000317
的第k个对角元素，则θr,k的估计为：
Figure FDA0002251914000000041
S313、计算At的估计：
Figure FDA0002251914000000042
令
Figure FDA0002251914000000043
和
Figure FDA0002251914000000044
分别代表
Figure FDA0002251914000000045
第k列的前三个和最后三个元素，则矢量叉积
Figure FDA0002251914000000046
φt,k的估计为：
Figure FDA0002251914000000047
其中，pr,k(1)和pr,k(2)分别表示pr,k的第一和第二个元素。
7.根据权利要求6所述基于平行因子算法的双基地EMVS-MIMO雷达角度估计算法，其特征在于，所述步骤S4中，2D-RPA估计为：
令
Figure FDA0002251914000000048
根据
Figure FDA0002251914000000049
构建方向矩阵Vr,k，则极化向量可以通过下式估计：
Figure FDA00022519140000000410
二维发射极化角可以通过下式来计算：
Figure FDA00022519140000000411
Figure FDA00022519140000000412
其中gr,k(1)和gr,k(2)分别表示gr,k的第一和第二个元素。
8.一种基于平行因子算法的双基地EMVS-MIMO雷达角度估计装置，其特征在于，所述装置包括：
模型构建模块：根据匹配滤波后的阵列信号模型构建三阶平行因子张量模型；
因子分解模块：对所述三阶平行因子张量模型进行平行因子分解，得到因子矩阵的估计值；
角度估计模块：根据因子矩阵的估计值，采用ESPRIT技术获得对方位角的估计，利用矢量叉乘的性质获得俯仰角的估计；利用最小二乘技术获得二维极化角的估计。
CN201911037437.1A 2019-10-29 2019-10-29 基于平行因子算法的双基地emvs-mimo雷达角度估计算法及装置 Active CN110907923B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911037437.1A CN110907923B (zh) 2019-10-29 2019-10-29 基于平行因子算法的双基地emvs-mimo雷达角度估计算法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911037437.1A CN110907923B (zh) 2019-10-29 2019-10-29 基于平行因子算法的双基地emvs-mimo雷达角度估计算法及装置

Publications (2)

Publication Number Publication Date
CN110907923A true CN110907923A (zh) 2020-03-24
CN110907923B CN110907923B (zh) 2022-07-05

Family

ID=69816133

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911037437.1A Active CN110907923B (zh) 2019-10-29 2019-10-29 基于平行因子算法的双基地emvs-mimo雷达角度估计算法及装置

Country Status (1)

Country Link
CN (1) CN110907923B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112269172A (zh) * 2020-10-14 2021-01-26 中国人民解放军空军工程大学 一种基于张量结构的嵌套mimo雷达角度估计方法和装置
CN112904296A (zh) * 2021-01-18 2021-06-04 西安电子科技大学 基于实值处理cs-parafac算法的双基地mimo雷达角度估计方法
US11860294B2 (en) 2020-08-24 2024-01-02 Google Llc Electromagnetic vector sensors for a smart-device-based radar system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106610483A (zh) * 2016-11-11 2017-05-03 长江大学 基于张量子空间和谱峰搜索的mimo雷达角度估计算法
CN107064892A (zh) * 2016-11-11 2017-08-18 长江大学 基于张量子空间和旋转不变的mimo雷达角度估计算法
CN108594194A (zh) * 2018-03-19 2018-09-28 长江大学 基于四线性分解的双基地mimo雷达角度估算方法
CN108919231A (zh) * 2018-05-25 2018-11-30 长江大学 联合误差条件下双基地mimo雷达角度估计方法
CN108983143A (zh) * 2018-05-23 2018-12-11 长江大学 有色噪声背景下双基地mimo雷达角度估计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106610483A (zh) * 2016-11-11 2017-05-03 长江大学 基于张量子空间和谱峰搜索的mimo雷达角度估计算法
CN107064892A (zh) * 2016-11-11 2017-08-18 长江大学 基于张量子空间和旋转不变的mimo雷达角度估计算法
CN108594194A (zh) * 2018-03-19 2018-09-28 长江大学 基于四线性分解的双基地mimo雷达角度估算方法
CN108983143A (zh) * 2018-05-23 2018-12-11 长江大学 有色噪声背景下双基地mimo雷达角度估计方法
CN108919231A (zh) * 2018-05-25 2018-11-30 长江大学 联合误差条件下双基地mimo雷达角度估计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
BAOQING XU ET AL.: "A novel unitary PARAFAC method for DOD and DOA estimation in bistatic MIMO radar", 《SIGNAL PROCESSING》 *
张小飞等: "电磁矢量阵中基于平行因子压缩感知的角度估计算法", 《数据采集与处理》 *
樊劲宇等: "偶极子分离的矢量阵MIMO雷达多维角度估计算法", 《电子与信息学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11860294B2 (en) 2020-08-24 2024-01-02 Google Llc Electromagnetic vector sensors for a smart-device-based radar system
CN112269172A (zh) * 2020-10-14 2021-01-26 中国人民解放军空军工程大学 一种基于张量结构的嵌套mimo雷达角度估计方法和装置
CN112269172B (zh) * 2020-10-14 2024-03-08 中国人民解放军空军工程大学 一种基于张量结构的嵌套mimo雷达角度估计方法和装置
CN112904296A (zh) * 2021-01-18 2021-06-04 西安电子科技大学 基于实值处理cs-parafac算法的双基地mimo雷达角度估计方法
CN112904296B (zh) * 2021-01-18 2023-12-26 西安电子科技大学 基于实值处理cs-parafac算法的双基地mimo雷达角度估计方法

Also Published As

Publication number Publication date
CN110907923B (zh) 2022-07-05

Similar Documents

Publication Publication Date Title
CN107290730B (zh) 互耦条件下双基地mimo雷达角度估算方法
CN107064892B (zh) 基于张量子空间和旋转不变的mimo雷达角度估计算法
CN110907923B (zh) 基于平行因子算法的双基地emvs-mimo雷达角度估计算法及装置
CN110412537A (zh) 一种双基地mimo雷达角度估计方法
CN104977558A (zh) 一种基于贝叶斯压缩感知的分布源中心波达方向估计方法
Sun et al. Target estimation in colocated MIMO radar via matrix completion
CN107064926B (zh) 空域色噪声背景下的双基地mimo雷达角度估算方法
Porozantzidou et al. Azimuth and elevation angles estimation using 2-D MUSIC algorithm with an L-shape antenna
CN109655799A (zh) 基于iaa的协方差矩阵向量化的非均匀稀疏阵列测向方法
CN110927711A (zh) 一种双基地emvs-mimo雷达的高精度定位算法及装置
Wu et al. Direction finding in bistatic MIMO radar with direction-dependent mutual coupling
CN108363048B (zh) 一种基于块稀疏的极化mimo雷达的角度估计方法
Mao et al. Target localization in bistatic EMVS-MIMO radar using tensor subspace method
CN108872930A (zh) 扩展孔径二维联合对角化doa估计方法
Grundmann et al. Investigation of direction of arrival estimation using characteristic modes
Li et al. Mutual coupling self-calibration algorithm for uniform linear array based on ESPRIT
Xu et al. Tensor‐Based Angle and Range Estimation Method in Monostatic FDA‐MIMO Radar
Tian et al. 2D-DOA estimation in arc-array with a DNN based covariance matrix completion strategy
Zou et al. High accuracy frequency and 2D-DOAs estimation of conformal array based on PARAFAC
CN110907926A (zh) 基于传播算子的双基地emvs-mimo雷达快速目标定位算法及装置
Shen et al. An improved propagator estimator for DOA estimation in monostatic MIMO radar in the presence of imperfect waveforms
Ma et al. Computation-efficient 2-D DOA estimation algorithm with array motion strategy
Gong et al. Three‐dimensional coprime array for massive MIMO: Array configuration design and 2D DOA estimation
Zheng et al. Joint DOD and DOA estimation for bistatic polarimetric MIMO radar
Gao et al. 2D-DOA estimation based on higher-order SVD-based using EMVS sparse array

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant