CN107064926A - 空域色噪声背景下的双基地mimo雷达角度估算方法 - Google Patents
空域色噪声背景下的双基地mimo雷达角度估算方法 Download PDFInfo
- Publication number
- CN107064926A CN107064926A CN201710033413.3A CN201710033413A CN107064926A CN 107064926 A CN107064926 A CN 107064926A CN 201710033413 A CN201710033413 A CN 201710033413A CN 107064926 A CN107064926 A CN 107064926A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- propagation operator
- mimo radar
- array
- 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
Links
- 238000011156 evaluation Methods 0.000 title abstract 2
- 239000011159 matrix material Substances 0.000 claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 40
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 230000014509 gene expression Effects 0.000 claims description 12
- 238000000638 solvent extraction Methods 0.000 claims description 2
- 230000000903 blocking effect Effects 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 81
- 238000004088 simulation Methods 0.000 description 10
- 239000013598 vector Substances 0.000 description 9
- 230000001629 suppression Effects 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 7
- 238000003491 array Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 6
- 238000000354 decomposition reaction Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000000342 Monte Carlo simulation Methods 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- IBSREHMXUMOFBB-JFUDTMANSA-N 5u8924t11h Chemical compound O1[C@@H](C)[C@H](O)[C@@H](OC)C[C@@H]1O[C@@H]1[C@@H](OC)C[C@H](O[C@@H]2C(=C/C[C@@H]3C[C@@H](C[C@@]4(O3)C=C[C@H](C)[C@@H](C(C)C)O4)OC(=O)[C@@H]3C=C(C)[C@@H](O)[C@H]4OC\C([C@@]34O)=C/C=C/[C@@H]2C)/C)O[C@H]1C.C1=C[C@H](C)[C@@H]([C@@H](C)CC)O[C@]11O[C@H](C\C=C(C)\[C@@H](O[C@@H]2O[C@@H](C)[C@H](O[C@@H]3O[C@@H](C)[C@H](O)[C@@H](OC)C3)[C@@H](OC)C2)[C@@H](C)\C=C\C=C/2[C@]3([C@H](C(=O)O4)C=C(C)[C@@H](O)[C@H]3OC\2)O)C[C@H]4C1 IBSREHMXUMOFBB-JFUDTMANSA-N 0.000 description 1
- 208000011231 Crohn disease Diseases 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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/88—Radar or analogous systems specially adapted for specific applications
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种空域色噪声背景下的双基地MIMO雷达角度估算方法,其通过构建接收阵列数据在匹配滤波后的协方差矩阵;并通过协方差矩阵构建传播算子,得到传播算子估计的最优解;根据传播算子估计的最优解获得阵列信号子空间的估计;最后通过子空间的旋转不变特性通过最小二乘的方法获取配对的角度。本发明所提算法无需高复杂度的SVD运算或者HOSVD运算,计算复杂度低;且能够获得自动配对的DOD与DOA,无需进一步进行配对计算,适用于任何协方差矩阵位置的空域色噪声背景,算法鲁棒性高。
Description
技术领域
本发明涉及一种雷达信号处理技术,具体的说涉及一种空域色噪声背景下的双基地MIMO雷达角度估算方法。
背景技术
多输入多输出(Multiple-input Multiple-output,MIMO)系统是未来通信、探测等领域的核心技术之一,也是近年来研究的热点技术之一。在雷达探测领域,MIMO雷达将是最具前景的下一代雷达系统。与现有相控阵雷达系统不同,MIMO雷达采用多发多收的信号处理模式,利用分集增益大大改善目标探测性能。相比传统相控阵雷达系统,MIMO雷达在分辨率、抗衰落性、可辨识性以及抑制噪声等方面具有潜在的优势。根据MIMO雷达收发阵元配置的不同,可以将MIMO雷达划分为两类:统计MIMO雷达和共址MIMO雷达。其中,统计MIMO雷达采样分布式收发阵元配置,即收发阵元广泛的分布在空间中,这种阵元分布模式可以有效的抑制目标的雷达截面系数闪烁的问题;共址MIMO雷达中的发射阵元和接收阵元一般采用阵列的形式分布,阵元间距往往较近,这种雷达可以获得高精度的目标方位信息。本发明主要关注双基地MIMO雷达,其为共址雷达的重要一类。
角度估计是共址MIMO雷达目标定位的关键任务之一,目前关于共址MIMO雷达角度估计已有很多优秀的算法。典型的算法有Capon算法、多重谱峰分类(Multiple SignalClassification,MUSIC)算法、基于旋转不变技术的参数估计(Estimation Method ofSignal Parameters via Rotational,ESPRIT)算法、传播算子(Propagator Method,PM)算法、高阶子空间分解(Higher Order Singular Value Decomposition,HOSVD)算法、平行因子(Parallel Factor,PARAFAC)算法、基于稀疏表示的估计算法等等。然而上述算法均假设接收噪声为理想高斯白噪声的背景。在实际工程中,由于阵元在空间分布的间距过小,很容易使阵元相互影响从而产生空域色噪声。在存在空域色噪声的条件下,上述角度估计算法性能会下降,特别是在低信噪比条件下,算法性能会严重下降甚至完全失效。双基地MIMO雷达中的色噪声问题已经引起部分学者的注意,目前已经有部分学者提出利用MIMO雷达的阵列数据互协方差抑制色噪声。按照抑制色噪声的原理的不同,可以将现有抑制色噪声的算法划分为两大类——空域协方差算法和时域协方差算法。前者主要是利用不同匹配滤波器输出的色噪声不相关的性质,后者主要利用了不同快拍的空域色噪声不相关特性。前者的典型代表有Jin等人的算法(Jin M,Liao G,Li J.Joint DOD and DOA estimation forbistatic MIMO radar[J].Signal Processing,2009,89(2):244-251.)Chen等人的算法(ANew Method for Joint DOD and DOA Estimation in bistatic MIMO Radar)、Jiang等人的算法(Jiang H,Zhang J K,Wong K M.Joint DOD and DOA Estimation for BistaticMIMO Radar in Unknown Correlated Noise[J].IEEE Transactions on VehicularTechnology,2015,64(11):5113-5125.)、Wang等人的算法(Wang X,Wang W,Li X,et al.Atensor-based subspace approach for bistatic MIMO radar in spatial colorednoise[J].Sensors,2014,14(3):3897-3907.)。其中Jin等人的算法虽然能有效抑制空域色噪声,但其仅适用于三个发射阵元的MIMO雷达;Chen等人将发射阵元分为两个子阵,利用两个子阵数据的互协方差抑制色噪声,该算法适用于多于三个阵元的MIMO雷达配置;Jiang等人将收发阵元均分成两个子阵,该算法需要估计四维角度;Wang等人的算法能在抑制色噪声的同时利用阵列数据的多维结构,从而提升角度估计精度。空域协方差算法的最大缺陷是阵列虚拟孔径的损失,从而降低了参数估计的精度,基于时域协方差的色噪声抑制算法能有效的避免该缺陷。这类算法主要有符渭波等人的算法(符渭波,苏涛,赵永波,等.空间色噪声环 境下基于时空结构的双基地MIMO雷达角度和多普勒频率联合估计方法[J].电子与信息学报,2011,33(7):1649-1654.)和Wen等人的算法(Wen F Q,Xiong X D,Su J,etal.Angle Estimation for bistatic MIMO Radar in The Presence of SpatialColored Noise,Signal Processing,accepted,2017,doi:10.1016/j.sigpro.2016.12.017.)。其中,Wen等人的算法能有效利用阵列数据的多维结构,获得比符渭波等人算法更好的角度数估计的精度。然而上述时域协方差算法均存在计算复杂度大的问题,如符渭波等人的算法需要对阵列协方差数据进行奇异值分解(Singular ValueDecomposition,SVD),Wen等人的算法需要构造阵列互协方差张量,并对张量数据进行HOSVD,计算复杂度相比SVD更大。在MIMO雷达探测中,收发阵元往往有成千上万个,这些阵元将产生海量的数据。无论是SVD与HOSVD均无法胜任实际雷达系统的计算要求,因而上述算法并不适用实际MIMO雷达系统,寻求一种高效的色噪声背景下参数估计算法对未来MIMO雷达工程化具有非常急迫的现实意义。
发明内容
鉴于以上原因,有必要提供一种无需高复杂度的SVD运算或者HOSVD运算,计算复杂度低;且算法鲁棒性高的空域色噪声背景下的双基地MIMO雷达角度估算方法。
本发明提供一种空域色噪声背景下的双基地MIMO雷达角度估算方法,所述空域色噪声背景下的双基地MIMO雷达角度估算方法包括如下步骤:
S1、构建接收阵列数据在匹配滤波后的协方差矩阵;
S2、通过协方差矩阵构建传播算子,并得到传播算子估计的最优解;
S3、根据传播算子估计的最优解获得阵列信号子空间的估计;
S4、通过子空间的旋转不变特性通过最小二乘的方法获取配对的角度。
本发明所述空域色噪声背景下的双基地MIMO雷达角度估算方法,其通过构建接收阵列数据在匹配滤波后的协方差矩阵;并通过协方差矩阵构建传播算子, 得到传播算子估计的最优解;根据传播算子估计的最优解获得阵列信号子空间的估计;最后通过子空间的旋转不变特性通过最小二乘的方法获取配对的角度。本发明所提算法无需高复杂度的SVD运算或者HOSVD运算,计算复杂度低;且能够获得自动配对的DOD与DOA,无需进一步进行配对计算,适用于任何协方差矩阵位置的空域色噪声背景,算法鲁棒性高。
附图说明
图1是双基地MIMO雷达角度估计示意图;
图2是本发明所述空域色噪声背景下的双基地MIMO雷达角度估算方法在一阶AR空域色噪声背景下的估计散点图;
图3是本发明所述空域色噪声背景下的双基地MIMO雷达角度估算方法在二阶AR空域色噪声背景下的估计散点图;
图4是本发明所述空域色噪声背景下的双基地MIMO雷达角度估算方法与其他算法在一阶AR空域色噪声背景下的RMSE比较;
图5是本发明所述空域色噪声背景下的双基地MIMO雷达角度估算方法与其他算法在一阶AR空域色噪声背景下的PSD对比;
图6是本发明所述空域色噪声背景下的双基地MIMO雷达角度估算方法与其他算法在二阶AR空域色噪声背景下的RMSE比较;
图7是本发明所述空域色噪声背景下的双基地MIMO雷达角度估算方法与其他算法在二阶AR空域色噪声背景下的PSD对比。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明,应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明提供一种空域色噪声背景下的双基地MIMO雷达角度估算方法,所述 空域色噪声背景下的双基地MIMO雷达角度估算方法包括如下步骤:
S1、构建接收阵列数据在匹配滤波后的协方差矩阵。
具体的,本发明所述空域色噪声背景下的双基地MIMO雷达角度估算方法的模型如图1所示。假设天线系统由M个发射阵元和N个接收阵元构成,二者均为线性阵列,且收发阵元间的间距均为λ2,λ为发射信号波长。假设第m(m=1,…,M)个发射阵元在一个脉冲重复周期内发射码长为Q的归一化基带编码信号且不同发射阵元的发射信号相互正交,若K个非相干点目标位于雷达阵列远场位置,且第k个(1≤k≤K)点目标的方位是 其中为目标相对发射天线阵列的波离角(Direction-of-Departure,DOD),θk为目标相对于接收阵列的波达角(Direction-of-Arrival,DOA)。考虑MIMO雷达的一个相干处理时间(coherent processing interval)内共包含L个脉冲,则第l(l=1,2,…,L)个脉冲的阵列接收信号为
式中,为接收方向矩阵,其第k(k=1,2,…,K)个接收导引矢量为导引矢量中的第n(n=1,2,…,N)个元素为 为发射方向矩阵,第k个发射导引矢量为m(m=1,2,…,M);为第l个快拍目标回波特性矢量,其包含了目标的雷达散射截面(RadarCross Section,RCS)、多普勒频率和脉冲重复频率等信息,并假设所有目标的RCS在l个接收快拍内满足Swerling-I(慢起伏)模型;为发射信号矩阵,为接收阵列天线接收的第l个快拍的噪声矩阵,并假设其满足空域色高斯模型,即Wl的列向量为独立同分布的 高斯随机变量,其均值为0,协方差为C,即其中vec(·)为矢量化操作算子,IQ表述维数为Q的单位矩阵,表示克罗内克积。对每个接收阵元均用sm/Q进行匹配滤波处理,并将阵列数据表述成矩阵形式。则匹配滤波输出结果可以被表述成
Y=[AT⊙AR]BT+N=ABT+N 表达式2
式中,可以被视为维数为MN×K的虚拟方向矩阵,其中⊙为Khatri-Rao积(按列克罗内克积),可被视为虚拟的导向矢量。B=[b1,b2,…,bL]T为目标特性矩阵,为噪声矩阵匹配滤波的结果,其第l个列向量可以表示为N(l)=vec(WSH)。
在传统子空间算法中,需要对接收阵列信号的协方差矩阵RY进行分解。在样本数量有限的条件下阵列的协方差矩阵RY的估计为
其中,RB=BTB*,RN=NNH/L。在不考虑目标具有相同的多普勒频率的条件下RB近似为对角矩阵。在传统高斯白噪声假设下,RN=σ2IMN,σ2为噪声功率。从而对进行子空间分解即可获得虚拟方向矩阵A的估计,进而获得目标角度估计。但是在高斯色噪声的影响下,RN不再为对角矩阵。此时,噪声子空间会与信号子空间相互串扰,从而导致子空间估计误差变大,空域色噪声对信号子空间的这种影响在低信噪比条件下表现尤为突出。在空域协方差抑噪方法中,利用了不同匹配滤波器阵列输出的噪声的不相关特性,将发射阵列分为两个子阵。这种抑制方法获得的阵列协方差矩阵的维数会降低,从而引起阵列孔径损失。为避免这个问题,本发明利用不同快拍输出噪声的非相关特性抑制空域色噪声。其抑制色噪声的原理以表述为如下:
上述公式的第一个等号处用到矢量化的性质vec(ABC)=(CTA)vec(B),第三个等号处用到了克罗内科积的性质根据表达式4可以看出,不同快拍所对应的阵列噪声是不相关的,为利用空域色噪声的这个特性,令Y1和Y2分别由Y的前L-1列和后Y的后L-1列构成。类似表达式3,阵列数据的协方差矩阵RZ的估计可以表示成
式中,B1、B1分别为B的前L-1行与后L-1行。可以看到,在中,色噪声已被很好的抑制。
S2、通过协方差矩阵构建传播算子,并得到传播算子估计的最优解。
所述步骤S2还包括以下分步骤:
S21、将阵列虚拟方向矩阵及阵列协方差矩阵分块。
为减少方向矩阵估计的计算复杂度,本发明采用PM思想来计算,将阵列虚拟方向矩阵及阵列协方差矩阵分块:
其中A1、A2分别为A的前K行与后MN-K行构成,G和H分别为维数为MN×K及MN×(MN-K)的数据矩阵。
S22、通过协方差矩阵构建传播算子,并得到传播算子的最小化代价函数。
由于A1是非奇异的,在无噪声条件下,存在唯一线性变换P使得
在存在噪声时,尽管数据仍然可以分块,但是表达式7已经不再满足。此时传播算子P可以通过最小化如下代价函数获得
其中,||·||F表示Frobenius范数。
S23、传播算子的最小化代价函数为传播算子的估计的二次凸函数,根据传播算子的最小化代价函数得到的最优解。
表达式8是关于的二次凸函数,因此,其最优解表示为
式中,表示伪逆运算。
S3、根据传播算子估计的最优解获得阵列信号子空间的估计。
具体的,令
根据表达式7可得
也就是说,由Q的列向量张成的子空间span{Q}与span{A}是一致的。因此,可将Q视为信号子空间的估计,即
Es=Q=[AT⊙AR]T 表达式12
式中,T是维数为K×K的非奇异矩阵。
S4、通过子空间的旋转不变特性通过最小二乘的方法获取配对的角度。
具体的,当获得信号子空间后,便可以利用ESPRIT思想获取目标的DOD与DOA估计。对于均匀阵列,其方向矩阵均是范德蒙的,因此有
上式中,JM1=[IM-1,0(M-1)×1],JM2=[0(M-1)×1,IM-1];JN1=[IN-1,0(N-1)×1],JN2=[0(N-1)×1,IN-1]。根据表达式12和表达式13可知,可以使用最小二乘方法获取ΦT和ΦR的估计
所估计的与ΦT、与ΦR之间满足关系故和有相同的特征矢量,因此可以通过如下方法进行DOD和DOA的自动配对。对进行特征值分解,即
其中,VT分别为的特征向量与特征值矩阵。则可通过下式获得的特征值矩阵
最后,令VT和VR的第k(k=1,2,…,K)个对角线上的元素分别为和 则第k个目标的DOD与DOA分别为
式中angle{·}表示取相位操作。
本发明所提算法的计算复杂度统计如下:
计算表达式5中的需要复杂度为M2N2(L-1)次复乘运算;
获取表达式9中的需要复杂度为2MNK2+Ο(K3)+MNK(MN-K)次复乘运算;
计算和共需要复杂度为2Ο(K3)+2(M-1)NK2+2(N-1)MK2次复乘运算;
估计DOD与DOA共需复杂度为Ο(K3)次复乘运算。
故本发明的计算复杂度为:
M2N2(L-1)+2MNK2+MNK(MN-K)+2(M-1)NK2+2(N-1)MK2+5Ο(K3)次复乘运算。
如表1所示,表1给出了本发明所提算法的计算复杂度和现有其他算法的计算复杂度:
表1
由表1可知,所提算法的计算复杂度远低于其他算法。
进一步的,技术人员针对本发明提出的空域色噪声背景下双基地MIMO雷达角度估计方法进行了大量的仿真实验。
在仿真实验中,假设K=3个目标处于远场,其DOA和DOD分别为 和三个目标的RCS系数满足系数均为1的Swerling I模型,多普勒频率分别为f1=200Hz、f2=400Hz和f3=800Hz。仿真实验中发射阵元的个数M=10,接收阵元数的个数N=6,二者均为均匀线性阵列,阵元间距均为λ/2。发射的基带编码波形矩阵为HM为维数是Q×Q的哈达码矩阵的前M行。编码长度Q、脉冲重复频率fs、快拍数分 别设置为Q=256,fs=20KHz、L=100。为验证本发明所提算法对色噪声抑制的有效性,仿真中设置了两种色噪声背景:场景A)、一阶自回归(Autoregressive,AR)模型,自相关系数为0.9;场景B)、二阶AR模型,模型系数为z=[1,-1,0.8];仿真中的信噪比(signal-to-noise ratio,SNR)定义为:
SNR=10log10(||Xl-Wl||2/||Wl||2)[dB]。
附图2和附图3分别为本发明所提算法在场景A和场景B背景下,500次蒙特卡洛仿真所得到角度估计结果,仿真中SNR=-15dB。由仿真结果可以看出,在不同的色噪声背景下,本发明所提的算法均能有效估计并正确配对所估计的DOD与DOA。
为比较本发明所提算法同现有算法估计精度的比较,对算法进行500次蒙特卡洛仿真,角度估计的精度用均方根误差(Root Mean Squared Error,RMSE)和成功检测概率(Probability of Successful Detection,PSD)来评价,其中RMSE定义为:
式中和分别为第i次蒙特卡洛仿真中获得的对θk与的估计;RE为正确检测的次数,若一次蒙特卡洛仿真中每个目标的DOD及DOA与真实的DOD及DOA之差的绝对值都小于0.3°,则记录该次仿真成功检测。
附图4与附图5分别给出了本发明所提算法和ESPRIT算法(标记为ESPRIT)、Chen等人的算法(标记为Chen)、符渭波等人的算法(标记为Fu)在一阶AR空域色噪声的背景下性能比较的结果。附图4为RMSE性能的比较,附图5给PSD性能的比较。由仿真结果可知,由于色噪声的影响,在信噪比较低时,ESPRIT算法的性能会有所下降,使用了抑噪方法的算法性能都会有所改善。但是应该注意到,Chen等人的空域抑噪方法会带来阵列虚拟孔径的损失,从而在SNR≥-5dB 后,该算法的RMSE性能会弱于ESPRIT算法。相比之下,时域抑噪算法在低信噪比条件下的RMSE低于ESPRIT,即估计精度高于ESPRIT。即使信噪比较高,所提算法的估计精度仍然很接近ESPRIT算法。此外,可以看出,所提算法的估计性能非常接近符渭波等人的算法,但其计算复杂度远低于他们的算法,故本发明所提算法的效率更高。
附图6与附图7分别给出了本发明所提算法在二阶AR空域色噪声背景下角度估计的性能同其他算法的比较。此时的仿真结果与一阶AR空域色噪声背景下的参数估计性能类似,所提算法的RMSE和PSD性能非常接近符渭波等人的算法,且在低信噪比条件下性能优于ESPRIT,在高信噪比条件下RMSE性能优于Chen等人的算法。因此,本发明所提的算法能非常有效的对抗不同类型的色噪声。
本发明所述空域色噪声背景下的双基地MIMO雷达角度估算方法,其通过构建接收阵列数据在匹配滤波后的协方差矩阵;并通过协方差矩阵构建传播算子,得到传播算子估计的最优解;根据传播算子估计的最优解获得阵列信号子空间的估计;最后通过子空间的旋转不变特性通过最小二乘的方法获取配对的角度。本发明所提算法无需高复杂度的SVD运算或者HOSVD运算,计算复杂度低;且能够获得自动配对的DOD与DOA,无需进一步进行配对计算,适用于任何协方差矩阵位置的空域色噪声背景,算法鲁棒性高。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种空域色噪声背景下的双基地MIMO雷达角度估算方法,其特征在于,所述空域色噪声背景下的双基地MIMO雷达角度估算方法包括如下步骤:
S1、构建接收阵列数据在匹配滤波后的协方差矩阵;
S2、通过协方差矩阵构建传播算子,并得到传播算子估计的最优解;
S3、根据传播算子估计的最优解获得阵列信号子空间的估计;
S4、通过子空间的旋转不变特性通过最小二乘的方法获取配对的角度。
2.根据权利要求1所述空域色噪声背景下的双基地MIMO雷达角度估算方法,其特征在于,所述接收阵列数据的协方差矩阵RZ的估计如下:
<mrow>
<msub>
<mi>R</mi>
<msup>
<mi>B</mi>
<mo>&prime;</mo>
</msup>
</msub>
<mo>=</mo>
<msubsup>
<mi>B</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<msubsup>
<mi>B</mi>
<mn>2</mn>
<mo>*</mo>
</msubsup>
</mrow>
<mrow>
<msub>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mi>Z</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<msub>
<mi>Y</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>Y</mi>
<mn>2</mn>
<mi>H</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>AR</mi>
<msup>
<mi>B</mi>
<mo>&prime;</mo>
</msup>
</msub>
<msup>
<mi>A</mi>
<mi>H</mi>
</msup>
</mrow>
式中,为接收方向矩阵,为发射方向矩阵,Y1和Y2分别为Y的前L-1列和Y的后L-1列,B1、B1分别为B的前L-1行与后L-1行。
3.根据权利要求2所述空域色噪声背景下的双基地MIMO雷达角度估算方法,其特征在于,步骤S2还包括以下分步骤:
S21、将阵列虚拟方向矩阵及阵列协方差矩阵分块;
S22、通过协方差矩阵构建传播算子,并得到传播算子的最小化代价函数;
S23、传播算子的最小化代价函数为传播算子的估计的二次凸函数,根据传播算子的最小化代价函数得到的最优解。
4.根据权利要求3所述空域色噪声背景下的双基地MIMO雷达角度估算方法,其特征在于,将阵列虚拟方向矩阵及阵列协方差矩阵分块的表达式为:
<mrow>
<mi>A</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mtable>
<mtr>
<mtd>
<mrow>
<mo>}</mo>
<mo>,</mo>
<mi>K</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>}</mo>
<mo>,</mo>
<mi>M</mi>
<mi>N</mi>
<mo>-</mo>
<mi>K</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mrow>
<mrow>
<msub>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mi>Z</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<mi>G</mi>
<mo>,</mo>
<mi>H</mi>
<mo>&rsqb;</mo>
</mrow>
其中A1、A2分别为A的前K行与后MN-K行,G和H分别为维数为MN×K及MN×(MN-K)的数据矩阵。
5.根据权利要求4所述空域色噪声背景下的双基地MIMO雷达角度估算方法,其特征在于,传播算子的估计的最优解表示为:
式中,P为传播算子,为传播算子P的估计,表示伪逆运算。
6.根据权利要求5所述空域色噪声背景下的双基地MIMO雷达角度估算方法,其特征在于,信号子空间的估计Q的表达式为:
式中,Q为信号子空间的估计,T是维数为K×K的非奇异矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710033413.3A CN107064926B (zh) | 2017-01-18 | 2017-01-18 | 空域色噪声背景下的双基地mimo雷达角度估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710033413.3A CN107064926B (zh) | 2017-01-18 | 2017-01-18 | 空域色噪声背景下的双基地mimo雷达角度估算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107064926A true CN107064926A (zh) | 2017-08-18 |
CN107064926B CN107064926B (zh) | 2020-05-19 |
Family
ID=59598891
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710033413.3A Expired - Fee Related CN107064926B (zh) | 2017-01-18 | 2017-01-18 | 空域色噪声背景下的双基地mimo雷达角度估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107064926B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108562866A (zh) * | 2018-04-04 | 2018-09-21 | 长江大学 | 基于矩阵填充的双基地mimo雷达角度估算方法 |
CN108594194A (zh) * | 2018-03-19 | 2018-09-28 | 长江大学 | 基于四线性分解的双基地mimo雷达角度估算方法 |
CN108983143A (zh) * | 2018-05-23 | 2018-12-11 | 长江大学 | 有色噪声背景下双基地mimo雷达角度估计方法 |
CN109116337A (zh) * | 2018-07-30 | 2019-01-01 | 西北工业大学 | 一种基于矩阵滤波的稀疏近似最小方差方位估计方法 |
CN109738877A (zh) * | 2019-01-30 | 2019-05-10 | 大连大学 | 冲击噪声环境下目标参数联合估计算法 |
CN110907926A (zh) * | 2019-10-29 | 2020-03-24 | 长江大学 | 基于传播算子的双基地emvs-mimo雷达快速目标定位算法及装置 |
CN111208502A (zh) * | 2020-01-13 | 2020-05-29 | 长江大学 | 一种无人驾驶物流车辆的定位方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217671A (zh) * | 2013-04-10 | 2013-07-24 | 哈尔滨工程大学 | 色噪声环境下的多输入多输出雷达收发角度快速估计方法 |
CN103760547A (zh) * | 2014-01-24 | 2014-04-30 | 西安电子科技大学 | 基于互相关矩阵的双基mimo雷达角度估计方法 |
CN103777197A (zh) * | 2013-12-24 | 2014-05-07 | 南京航空航天大学 | 单基地mimo雷达中降维传播算子的方位估计方法 |
CN103777190A (zh) * | 2014-02-26 | 2014-05-07 | 南京信息工程大学 | 一种双基地mimo雷达高速高机动目标的角度估计方法 |
CN103926573A (zh) * | 2014-04-17 | 2014-07-16 | 哈尔滨工程大学 | 基于四阶累积量的单基地mimo雷达分布式目标角度估计方法 |
CN104215947A (zh) * | 2014-08-18 | 2014-12-17 | 哈尔滨工业大学深圳研究生院 | 一种双基地mimo雷达角度的估计方法 |
CN105259550A (zh) * | 2015-11-17 | 2016-01-20 | 西安电子科技大学 | 基于压缩感知的多输入多输出雷达二维角度估计方法 |
-
2017
- 2017-01-18 CN CN201710033413.3A patent/CN107064926B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217671A (zh) * | 2013-04-10 | 2013-07-24 | 哈尔滨工程大学 | 色噪声环境下的多输入多输出雷达收发角度快速估计方法 |
CN103777197A (zh) * | 2013-12-24 | 2014-05-07 | 南京航空航天大学 | 单基地mimo雷达中降维传播算子的方位估计方法 |
CN103760547A (zh) * | 2014-01-24 | 2014-04-30 | 西安电子科技大学 | 基于互相关矩阵的双基mimo雷达角度估计方法 |
CN103777190A (zh) * | 2014-02-26 | 2014-05-07 | 南京信息工程大学 | 一种双基地mimo雷达高速高机动目标的角度估计方法 |
CN103926573A (zh) * | 2014-04-17 | 2014-07-16 | 哈尔滨工程大学 | 基于四阶累积量的单基地mimo雷达分布式目标角度估计方法 |
CN104215947A (zh) * | 2014-08-18 | 2014-12-17 | 哈尔滨工业大学深圳研究生院 | 一种双基地mimo雷达角度的估计方法 |
CN105259550A (zh) * | 2015-11-17 | 2016-01-20 | 西安电子科技大学 | 基于压缩感知的多输入多输出雷达二维角度估计方法 |
Non-Patent Citations (6)
Title |
---|
BIN YANG等: "Angle estimation for bistatic MIMO radar based on core tensor PARAFAC decomposition", 《2016 IEEE INTERNATIONAL CONFERENCE ON UBIQUITOUS WIRELESS BROADBAND (ICUWB)》 * |
FANG SHENGLIANG等: "A New Algorithm for Revising Noise Covariance Matrix Disturbance in the Presence of Coherent Sources", 《 2009 2ND INTERNATIONAL CONGRESS ON IMAGE AND SIGNAL PROCESSING》 * |
刘兆霆等: "色噪声背景下基于声矢量阵列孔径扩展的相干目标波达方向估计", 《兵工学报》 * |
孙中伟等: "双基地MIMO雷达中基于传播算子的DOD和DOA估计算法", 《数据采集与处理》 * |
符渭波等: "空间色噪声环境下双基地MIMO雷达角度和多普勒频率联合估计方法", 《电子与信息学报》 * |
钟原: "双基地MIMO雷达的角度估计方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108594194A (zh) * | 2018-03-19 | 2018-09-28 | 长江大学 | 基于四线性分解的双基地mimo雷达角度估算方法 |
CN108562866A (zh) * | 2018-04-04 | 2018-09-21 | 长江大学 | 基于矩阵填充的双基地mimo雷达角度估算方法 |
CN108562866B (zh) * | 2018-04-04 | 2022-06-14 | 长江大学 | 基于矩阵填充的双基地mimo雷达角度估算方法 |
CN108983143A (zh) * | 2018-05-23 | 2018-12-11 | 长江大学 | 有色噪声背景下双基地mimo雷达角度估计方法 |
CN109116337A (zh) * | 2018-07-30 | 2019-01-01 | 西北工业大学 | 一种基于矩阵滤波的稀疏近似最小方差方位估计方法 |
CN109116337B (zh) * | 2018-07-30 | 2023-01-24 | 西北工业大学 | 一种基于矩阵滤波的稀疏近似最小方差方位估计方法 |
CN109738877A (zh) * | 2019-01-30 | 2019-05-10 | 大连大学 | 冲击噪声环境下目标参数联合估计算法 |
CN109738877B (zh) * | 2019-01-30 | 2023-03-31 | 大连大学 | 冲击噪声环境下目标参数联合估计方法 |
CN110907926A (zh) * | 2019-10-29 | 2020-03-24 | 长江大学 | 基于传播算子的双基地emvs-mimo雷达快速目标定位算法及装置 |
CN111208502A (zh) * | 2020-01-13 | 2020-05-29 | 长江大学 | 一种无人驾驶物流车辆的定位方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107064926B (zh) | 2020-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wen et al. | Joint 2D-DOD, 2D-DOA, and polarization angles estimation for bistatic EMVS-MIMO radar via PARAFAC analysis | |
CN108562866B (zh) | 基于矩阵填充的双基地mimo雷达角度估算方法 | |
CN107064926B (zh) | 空域色噪声背景下的双基地mimo雷达角度估算方法 | |
CN107064892B (zh) | 基于张量子空间和旋转不变的mimo雷达角度估计算法 | |
CN106610483B (zh) | 基于张量子空间和谱峰搜索的mimo雷达角度估计算法 | |
CN107015213B (zh) | 基于music算法的双基地mimo雷达角度估算方法 | |
Zhang et al. | A method for finding best channels in beam-space post-Doppler reduced-dimension STAP | |
CN108983143A (zh) | 有色噪声背景下双基地mimo雷达角度估计方法 | |
Xudong et al. | Improved ESPRIT Method for Joint Direction‐of‐Arrival and Frequency Estimation Using Multiple‐Delay Output | |
CN105319545A (zh) | 提高stap检测性能的mimo-ofdm雷达波形设计方法 | |
CN112379327A (zh) | 一种基于秩损估计的二维doa估计与互耦校正方法 | |
CN108594194A (zh) | 基于四线性分解的双基地mimo雷达角度估算方法 | |
Li et al. | Parameter estimation based on fractional power spectrum density in bistatic MIMO radar system under impulsive noise environment | |
CN107132514A (zh) | 改进esprit的双基地mimo雷达角度估算方法 | |
Gong et al. | Joint DOD and DOA estimation for bistatic multiple‐input multiple‐output radar target discrimination based on improved unitary ESPRIT method | |
WO2018014197A1 (zh) | 一种信道估计方法及装置 | |
CN114415137A (zh) | 一种基于矩阵束原理的频率分集mimo雷达参数估计方法 | |
Cheng et al. | Joint 4-D angle and Doppler shift estimation via tensor decomposition for MIMO array | |
Du et al. | Bayesian robust tensor factorization for angle estimation in bistatic MIMO radar with unknown spatially colored noise | |
CN113075649A (zh) | 一种适用于分布式网络化雷达的信号级直接定位方法 | |
Xu et al. | Tensor‐Based Angle and Range Estimation Method in Monostatic FDA‐MIMO Radar | |
Li et al. | DOD and DOA estimation for MIMO radar based on combined MUSIC and sparse Bayesian learning | |
Lonkeng et al. | Two‐Dimensional DOA Estimation Using Arbitrary Arrays for Massive MIMO Systems | |
Li | Joint parameter estimation and target localization for bistatic MIMO radar system in impulsive noise | |
Yu et al. | Multi‐mode propagation mode localisation and spread‐Doppler clutter suppression method for multiple‐input multiple‐output over‐the‐horizon radar |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200519 |
|
CF01 | Termination of patent right due to non-payment of annual fee |