CN107167785B - 一种稳健的大阵列mimo雷达目标收发角联合估计方法 - Google Patents
一种稳健的大阵列mimo雷达目标收发角联合估计方法 Download PDFInfo
- Publication number
- CN107167785B CN107167785B CN201710563441.6A CN201710563441A CN107167785B CN 107167785 B CN107167785 B CN 107167785B CN 201710563441 A CN201710563441 A CN 201710563441A CN 107167785 B CN107167785 B CN 107167785B
- Authority
- CN
- China
- Prior art keywords
- matrix
- transmitting
- receiving
- angle
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 115
- 239000011159 matrix material Substances 0.000 claims abstract description 119
- 230000008602 contraction Effects 0.000 claims abstract description 14
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 12
- 238000005516 engineering process Methods 0.000 claims abstract description 10
- 238000001228 spectrum Methods 0.000 claims abstract description 10
- 239000013598 vector Substances 0.000 claims description 20
- 230000003595 spectral effect Effects 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 16
- 230000007812 deficiency Effects 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 229940050561 matrix product Drugs 0.000 claims description 6
- 238000005315 distribution function Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000003491 array Methods 0.000 claims description 2
- 238000006467 substitution reaction Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 240000001973 Ficus microcarpa Species 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000005658 nuclear physics 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
- 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/414—Discriminating targets with respect to background clutter
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Peptides Or Proteins (AREA)
Abstract
本发明提供一种稳健的大阵列MIMO雷达目标收发角联合估计方法,属于多输入多输出雷达目标参数估计领域。该方法基于线性收缩和随机矩阵理论,其实现步骤包括:利用回波信号构造观测数据的随机矩阵模型,在非高斯噪声背景下,通过定点迭代求得其协方差矩阵的极大似然估计,当收发阵元数大于快拍数时,利用线性收缩技术得到稳健的协方差矩阵估计,利用围线积分、Stieltjes变换以及留数定理构造稳健的代价函数,对代价函数进行谱峰搜索求解目标收发角。优点是所提出方法在非高斯噪声背景下具有稳健的估计性能,同时适用于收发阵元数大于快拍数的大阵列***。
Description
技术领域
本发明属于多输入多输出(MIMO)雷达目标参数估计领域;涉及收发阵元数与快拍数相接近的大阵列MIMO雷达***;涉及非高斯噪声背景下双基地MIMO雷达的目标收发角联合估计;涉及收发阵元数大于快拍数情况下的MIMO雷达目标收发角度估计。
背景技术
雷达(Radar)通过发射无线电波来检测目标并测定目标的空间位置,从而对目标进行定位。作为一种新体制雷达,多输入多输出(MIMO)雷达由多个发射天线同时发射多个信号探测目标,并使用多个接收天线接收目标回波信号,通过采用多种可区分的发射信号(本发明采用正交信号),为雷达引入参数上有差异的目标观测通道。由于各收发通道之间采用相干处理,MIMO雷达可充分利用收发阵元获得虚拟孔径,提高收发角估计的性能。
MIMO雷达的目标收发角联合估计问题,是指在噪声背景下对阵列接收的回波信号进行有效处理,进而对目标的发射角(DOD)和接收角(DOA)进行估计。根据发射接收天线的阵列布局不同,MIMO雷达可分为分布式MIMO雷达和集中式MIMO雷达两类。分布式MIMO雷达的发射天线之间距离较远以获得空间分集增益;集中式MIMO雷达的发射阵列和接收阵列的各个阵元间距通常在半波长以内,通过相干处理提高估计增益,比如单/双基地MIMO雷达。在双基地MIMO雷达中,利用收发角联合估计可实现目标定位。现有的收发角联合估计方法主要有2D-Capon方法、2D-MUSIC方法、分维ESPRIT方法等。这些方法虽然能有效估计出收发角度,但仍存在很多问题亟需解决,主要包括:(1)现有方法通常将噪声模型简化为理想的高斯噪声,其分布函数为正态分布;而在实际环境中,噪声不是简单的高斯噪声,而是表现为一些非高斯特性,比如服从重尾分布,若仍采用高斯噪声背景下的目标参数估计方法,将导致估计性能严重下降;(2)现有方法的前提条件都是假定收发阵元数远远小于快拍数,而当收发阵元数大于快拍数时,在求解采样协方差矩阵定点迭代的过程中会出现因秩缺失而无解的问题,此时无法对收发角度进行有效估计;(3)现有方法所考虑的MIMO雷达***均假设快拍数足够多,而实际中,当快拍数与收发阵元数相接近时,现有估计方法不再适用。在实际的雷达例如超视距MIMO雷达***中,往往会出现收发阵元数与快拍数相接近甚至收发阵元数大于快拍数的情况,即大阵列***。迄今为止,非高斯噪声背景下大阵列MIMO雷达***的目标收发角联合问题一直是一个难题。
近年来,随机矩阵理论在金融数学、核物理、数字通信等领域应用日益广泛,已成为学术界的一个新的研究热点。然而,随机矩阵理论在大阵列MIMO雷达目标参数估计领域的应用研究甚少,尤其是在非高斯噪声背景且收发阵元数大于快拍数情况下,针对大阵列双基地MIMO雷达的目标收发角联合估计方法,国内外尚未见相关成果的研究报道。
随机矩阵渐进谱理论是符合理想高斯噪声背景下大阵列MIMO雷达目标收发角联合估计的有效方法。然而,实际环境中,噪声往往不再服从正态分布,而是呈现重尾分布。针对非高斯噪声情况,由于采样协方差矩阵将不再是统计协方差矩阵的一致性估计,导致高斯噪声背景下提出的收发角联合估计方法不再适用。为解决这一问题,人们利用极大似然估计对采样协方差矩阵进行定点迭代,作为协方差矩阵的一致估计。然而,当阵元数大于快拍数时,定点迭代的过程中会出现因秩缺失而无解的情况,这一问题是目前亟待解决的难题。
近年来出现的基于线性收缩(linear Shrinkage)的协方差矩阵估计方法,已成为一种新的信号处理手段。收缩估计可以在矩阵维数大于快拍数的情况下,对大维采样协方差矩阵的特征值分布区间进行收缩,解决秩缺失的问题。目前,已经有文献研究了大维***中协方差矩阵的线性收缩估计问题以及一些改进的方法,广义收缩估计器已经应用于空时自适应处理(STAP)检测问题中。然而,在大阵列MIMO雷达收发角联合估计领域,国内外尚未见任何有关线性收缩应用的报道。
发明内容
本发明提供一种稳健的大阵列MIMO雷达目标收发角联合估计方法,以解决估计性能严重下降、无法对收发角度进行有效估计、非高斯噪声背景下大阵列MIMO雷达***的目标收发角联合估计困难的问题。
本发明采取的技术方案是:建立双基地MIMO雷达的***模型,M个发射阵元同时发射窄带信号照射K个远场目标,在接收端采用N个接收阵元,其中发射阵列和接收阵列均采用均匀线阵,阵元间距dt=dr=λ/2,发射角度为θ,接收角度为第k个目标的发射角度和接收角度分别为θk和L为快拍数,其特征在于包括下列步骤:
步骤一、MIMO雷达接收端N个接收阵元对多次快拍得到的回波数据信号进行收集和处理,第l个快拍期间回波数据可以写成一个N×Lt的矩阵X(l):
其中,表示M行1列的发射方向矢量,t表示发射端,表示N行1列的接收方向矢量,r表示接收端;S表示发射波形矩阵,这里sm(m=1,…,M)表示第m个发射阵元所发射的波形,其编码长度为Lt,且S为归
一化正交矩阵,SSH=IM;表示在第l个快拍期间第k个目标的散射系数;
表示一个N行Lt列的非高斯噪声矩阵,服从Student-T分布;
步骤二、首先对N个接收阵元所接收的信号X(l)进行匹配滤波处理,得到:
其中,为N行M列的匹配滤波后的信号;
步骤三、对匹配滤波后的结果进行矢量化处理,即y(l)=vec(X(l)SH),在第l个快拍期间,经过预处理后的观测数据可表示为:
其中,为目标散射系数矩阵;是在第l个快拍期间处理后的噪声矢量,维数为MN行1列;是联合导向矢量矩阵,维数为MN行K列at(θk)和分别为发射和接收方向矢量,是Khatri-Rao矩阵积,是Kronecker矩阵积;
步骤四、收集L个快拍期间的数据y(l)(l=1,…,L),可得到一个观测数据矩阵Y,它是一个维数为MN×L的随机矩阵,即:
其中,Y=[y(1),…,y(L)],B=[b(1),…,b(L)],N=[n(1),…,n(L)],这里表示包含收发角度的联合导向矩阵,B是服从均值为0,方差为1的目标散射系数矩阵,N为MN×L维的服从Student-T分布非高斯噪声矩阵,由高斯噪声矩阵的元素乘上服从伽马分布Γ(v,1/v)的因子构成,其中,v=0.02;
步骤五、接收信号的采样协方差矩阵可以表示为求解协方差矩阵的极大似然估计;
步骤六、当收发阵元数大于快拍数时,L<MN,上述极大似然估计方法因为秩缺失而导致无解,此时采用线性收缩方法使得定点迭代的过程能持续进行;
步骤七、在获得稳健协方差矩阵估计的基础上,利用2D-MW方法可以得到收发角联合估计的谱峰搜索函数;
步骤八、当收发阵元数与快拍数在同一数量级时,即MN→∞,L→∞,MN/L=c,c为一个常数,Y可以看作一个大维随机矩阵,传统方法性能急剧下降,不能正确估计出收发角度;因此,应用随机矩阵理论中的渐近谱理论和围线积分技术来产生一个稳健的谱峰搜索函数;
步骤九、对进行谱峰搜索可以获得发射角θ和接收角
本发明所述步骤五的具体实现方法是:对接收信号y(l)进行归一化yi的概率密度函数为其中U(MN/2)/2πMN/2为单位圆上的均匀密度分布函数,对p(yi:R)求导并令导数等于零可以得到经过定点迭代可以得到极大似然估计定点迭代的初始矩阵可以是任意正定矩阵,令为单位矩阵。
本发明所述步骤六的收缩方法如下, 其中0<ρ<1为收缩系数,将定点迭代完成之后最终收敛的协方差矩阵记为收缩系数ρo由最小均方误差准则得到通过化简可以解得收缩系数为:
由于R未知,用替代R可以得到:
本发明所述步骤七具体实现方法是:对采样协方差矩阵进行特征值分解得到
其中,(λmn>0),EEH=I。定义0<z<∞上的一个连续的函数f(z),其反函数为可以得到
收发角联合估计估计的谱峰搜索函数可以表示为
改变函数f(z)可以得到不同的谱峰搜索函数,令函数f(z)=ln(z),得到Midway(MW)方法的代价函数
本发明所述步骤八具体实现方法是:首先根据留数定理,可以表示为:
其中,ζ-是一条在上的顺时针方向曲线,并且该曲线包含R的所有特征值;
对曲线ζ-进行参数化处理,令其中,是一个包含所有特征值的矩形,表示该矩形的边界,wMN(z)与特征值有关,其关系式如下:
为了简化公式,定义两个变量bMN(z)和mMN(z)来分别描述特征值和特征向量:
得到参数化表示形式为
通过变量bMN(z),得到wMN(z)的导数对上式进一步化简得到:
当MN→∞,L→∞,MN/L=c时,用bMN(z)和mMN(z)来描述特征值和特征向量不再准确,此时需要研究特征值和特征向量的渐近表现,在随机矩阵理论中,采样协方差矩阵特征值的瞬时分布函数的Stieltjes变换定义如下:
仅仅体现了采样协方差矩阵的特征值渐近特性,为了同时体现特征值和特征向量的渐近特性,定义如下函数:
此时,我们用和替换bMN(z)和mMN(z),得到:
其中,
为了求解定义一个关于的函数:
其中,x∈[0,c];
在x∈(0,c)区间上对求导:
利用留数定理可以得到:
其中,是下面等式的解:
可以由柯西留数定理直接得到,再此基础上对求积分,可以得到如下:
令x=c,经过化简可以得到ξmn(c):
其中,是下面等式的解:
至此,基于线性收缩和随机矩阵理论的MIMO雷达收发角联合估计谱峰搜索函数可以写为:
其中,如上所述。
本发明在非高斯背景下,且以Student-T噪声为例,计算接收信号的采样协方差矩阵,并利用定点迭代技术产生其协方差矩阵的一致估计,即:其中,yi为归一化接收信号向量,j为迭代的次数。当收发阵元数的乘积MN大于快拍数L时,定点迭代的过程中会出现因秩缺失而无法求解采样协方差矩阵的问题,因此本发明采用线性收缩的方法解决这一问题,得到稳健的采样协方差矩阵估计:
其中,ρ为收缩系数,本发明中采用最小均方误差准则对其进行求解。
本发明将基于Pisarenko结构的Midway(MW)方法应用到MIMO雷达收发角联合估计当中,提出一种2D-MW方法,当收发阵元数和快拍数在同一数量级时,采用随机矩阵理论中的围线积分技术、Stieltjes变换以及留数定理,推导了大维条件下稳健的代价函数,提出一种基于随机矩阵理论的2D-RMW方法,之后通过谱峰搜索可以得到精确的收发角。在相同仿真条件下,将2D-RMW方法与2D-MW方法进行了对比,验证了方法的有效性,代价函数如下所示:
其中,ξmn(mn=1,…,MN)为基于随机矩阵理论推导出的谱峰搜索函数的修正值,为联合导向矢量,为特征向量。
本发明在非高斯噪声背景且收发阵元数大于快拍数的情况下,将线性收缩技术与随机矩阵理论结合,通过收缩得到采样协方差矩阵的一致估计,再利用高斯噪声背景下基于随机矩阵理论的2D-RMW方法,提出了一种新的收发角联合估计方法----2D-RLSMW方法,代价函数如下:
其中,为基于线性收缩和随机矩阵理论结合推导出的谱峰搜索函数的修正值,LS指的是线性收缩,为联合导向矢量,为特征向量。
本发明为了分析发明方法的收发角联合估计性能,在大阵列双基地MIMO雷达***中,通过Matlab仿真实验证明了高斯环境背景下2D-RMW方法的有效性,检验了Student-T分布的非高斯噪声背景下且收发阵元数大于快拍数时本发明所提出的2D-RLSMW方法的优越性。
本发明的优点在于:(1)在大阵列MIMO雷达收发角联合估计中,利用极大似然估计,通过定点迭代产生协方差矩阵的一致估计,解决非高斯噪声背景下收发角估计不准确的问题;(2)首次将线性收缩技术应用到大阵列MIMO雷达目标收发角联合估计中,解决了收发阵元数大于快拍数时,求解稳健协方差矩阵的过程中产生的秩缺失问题;(3)首次将用于瞬时功率估计中的MW方法应用到大阵列MIMO雷达收发角联合估计中,并采用随机矩阵的渐进特征值分布理论作为数学工具,以回波数据样本协方差矩阵的对数形式为基础,基于随机矩阵的围线积分技术和Stieltjes变换推导新的代价函数,在收发阵元数和快拍数在同一数量级时具有较好的性能。
本发明提出的方法得到国家自然科学基金项目“基于大维随机矩阵理论的MIMO雷达稳健目标检测与估计(No.61371158)”的资助。
附图说明
图1是本发明中所用到的双基地MIMO雷达***模型;
图2是在Student-T分布的非高斯噪声背景下基于线性收缩和随机矩阵理论的大阵列双基地MIMO雷达目标收发角联合估计方法流程图;
图3是为高斯噪声背景下2D-MW方法、2D-RMW方法和2D-RLSMW三种方法的均方根误差随信噪比变化的曲线,仿真参数设置如下:MN=36,L=50;
图4是在Student-T分布的非高斯噪声背景下2D-MW方法、2D-RMW方法和2D-RLSMW三种方法的均方根误差随信噪比变化的曲线,仿真参数设置如下:MN=36,L=50;
图5(a)是在Student-T分布的非高斯噪声背景下利用2D-RMW方法进行收发角联合估计的谱峰搜索图,仿真参数设置如下:MN=36,L=30;
图5(b)是图5(a)的等高线图;
图6(a)是在Student-T分布的非高斯噪声背景下利用2D-RLSMW方法进行收发角联合估计的谱峰搜索图,仿真参数设置如下:MN=36,L=30;
图6(b)是图6(a)的等高线图;
图7在Student-T分布的非高斯噪声背景下分别用2D-RMW方法和2D-RLSMW方法进行收发角估计的均方根误差随信噪比变化的仿真曲线,仿真参数设置如下:MN=36,L=30。
具体实施方式
图1是双基地MIMO雷达的***模型,M个发射阵元同时发射窄带信号照射K个远场目标,在接收端通过N个接收阵元对回波信号进行综合处理,本发明中发射阵列和接收阵列均采用均匀线阵,阵元间距dt=dr=λ/2,发射角度为θ,接收角度为第k个目标的发射角度和接收角度分别为θk和快拍数为L。图2是在Student-T分布的非高斯噪声背景下基于线性收缩和随机矩阵理论的大阵列双基地MIMO雷达目标收发角估计方法的流程图,其具体实施步骤如下:
步骤一:MIMO雷达接收端N个接收阵元对多次快拍得到的回波数据信号进行收集和处理。第l个快拍期间回波数据可以写成一个N×Lt的矩阵X(l):
其中,表示M行1列的发射方向矢量,t表示发射端;表示N行1列的接收方向矢量,r表示接收端;S表示发射波形矩阵,这里sm(m=1,…,M)表示第m个发射阵元所发射的波形,其编码长度为Lt,且S为归一化正交矩阵,SSH=IM;表示在第l个快拍期间第k个目标的散射系数;表示一个N行Lt列的非高斯噪声矩阵,服从Student-T分布;
步骤二:首先对N个接收阵元所接收的信号X(l)进行匹配滤波处理,得到:
其中,为N行M列的匹配滤波后的信号;
步骤三:对匹配滤波后的结果进行矢量化处理,即y(l)=vec(X(l)SH)。在第l个快拍期间,经过预处理后的观测数据可表示为
其中,为目标散射系数矩阵;是在第l个快拍期间处理后的噪声矢量,维数为MN行1列;是联合导向矢量矩阵,维数为MN行K列at(θk)和分别为发射和接收方向矢量,是Khatri-Rao矩阵积,是Kronecker矩阵积;
步骤四:收集L个快拍期间的数据y(l)(l=1,…,L),可得到一个观测数据矩阵Y,它是一个维数为MN×L的随机矩阵,即:
其中,Y=[y(1),…,y(L)],B=[b(1),…,b(L)],N=[n(1),…,n(L)],这里表示包含收发角度的联合导向矩阵,B是服从均值为0,方差为1的目标散射系数矩阵,N为MN×L维的服从Student-T分布非高斯噪声矩阵,由高斯噪声矩阵的元素乘上服从伽马分布Γ(v,1/v)的因子构成,其中,v=0.02;
步骤五:接收信号的采样协方差矩阵可以表示为求解协方差矩阵的极大似然估计,对接收信号y(l)进行归一化yi的概率密度函数为其中U(MN/2)/2πMN/2为单位圆上的均匀密度分布函数,对p(yi:R)求导并令导数等于零可以得到经过定点迭代可以得到极大似然估计定点迭代的初始矩阵可以是任意正定矩阵,本发明中令为单位矩阵。
步骤六:当收发阵元数大于快拍数时,L<MN,上述极大似然估计方法因为秩缺失而导致无解,此时采用线性收缩技术使得定点迭代的过程能持续进行;收缩过程如下:其中0<ρ<1为收缩系数,将定点迭代完成之后最终收敛的协方差矩阵记为收缩系数ρo由最小均方误差准则得到通过化简可以解得收缩系数为由于R未知,用替代R可以得到
步骤七:在获得稳健协方差矩阵估计的基础上,利用2D-MW方法可以得到收发角联合估计的谱峰搜索函数。MW方法是一种Pisarenko结构方法,用于瞬时功率估计,本发明将其应用到双基地MIMO雷达目标收发角联合估计中,提出2D-MW方法。Pisarenko结构过程如下,对采样协方差矩阵进行特征值分解得到:
其中,(λmn>0),EEH=I,定义0<z<∞上的一个连续的函数f(z),其反函数为可以得到
收发角联合估计估计的谱峰搜索函数可以表示为
改变函数f(z)可以得到不同的谱峰搜索函数,令函数f(z)=ln(z),得到Midway(MW)方法的代价函数
步骤八:当收发阵元数与快拍数在同一数量级时,即MN→∞,L→∞,MN/L=c,c为一个常数,Y可以看作一个大维随机矩阵,传统方法性能急剧下降,不能正确估计出收发角度,因此,本发明应用随机矩阵理论中的渐近谱理论和围线积分技术来产生一个稳健的谱峰搜索函数,方法如下:首先根据留数定理,可以表示为
其中,ζ-是一条在上的顺时针方向曲线,并且该曲线包含R的所有特征值。
对曲线ζ-进行参数化处理,令其中,是一个包含所有特征值的矩形,表示该矩形的边界,wMN(z)与特征值有关,其关系式如下:
为了简化公式,定义两个变量bMN(z)和mMN(z)来分别描述特征值和特征向量:
得到参数化表示形式为
通过变量bMN(z),得到wMN(z)的导数对上式进一步化简得到:
当MN→∞,L→∞,MN/L=c时,用bMN(z)和mMN(z)来描述特征值和特征向量不再准确,此时需要研究特征值和特征向量的渐近表现。在随机矩阵理论中,采样协方差矩阵特征值的瞬时分布函数的Stieltjes变换定义如下:
仅仅体现了采样协方差矩阵的特征值渐近特性,为了同时体现特征值和特征向量的渐近特性,定义如下函数:
此时,我们用和替换bMN(z)和mMN(z),得到:
其中,
为了求解定义一个关于的函数:
其中,x∈[0,c];
在x∈(0,c)区间上对求导:
利用留数定理可以得到:
其中,是下面等式的解:
可以由柯西留数定理直接得到,再此基础上对求积分,可以得到如下:
令x=c,经过化简可以得到ξmn(c):
其中,是下面等式的解:
至此,基于线性收缩和随机矩阵理论的MIMO雷达收发角联合估计谱峰搜索函数可以写为:
其中,如上所述;
步骤九:对进行谱峰搜索可以获得发射角θ和接收角
下面通过Matlab仿真实验来进一步说明本发明的优点。分别在高斯噪声背景下和Student-T分布的非高斯噪声背景下比较了2D-MW、基于随机矩阵理论的2D-RMW、基于线性收缩和随机矩阵理论相结合的2D-RLSMW的性能,验证了本发明所提出方法的有效性。考虑发射阵元数为M=6、接收阵元数为N=6,收发阵列为均匀阵列,阵元间距dt=dr=λ/2,K=2个远场目标,目标位置选取快拍数分别取L=50和L=30。
图3和图4分别为高斯噪声背景下和Student-T分布的非高斯噪声背景下,当收发阵元数和快拍数在同一数量级时,2D-MW方法、2D-RMW方法和2D-RLSMW三种方法的均方根误差随信噪比变化的曲线。仿真条件如下:噪声矩阵N为非高斯噪声,服从Student-T分布,矩阵由高斯噪声矩阵的元素乘上服从伽马分布Γ(v,1/v)的因子构成,其中,v=0.02,快拍数L=50,两个目标的坐标为进行40次蒙特卡罗仿真实验。从图3可以看出,2D-RMW方法和2D-RLSMW方法均方根误差曲线相近,说明在高斯噪声环境下两种方法的估计性能相近;从图4可以看出,相同信噪比下,2D-RLSMW方法的均方根误差要小于2D-RMW方法,说明在Student-T分布的非高斯噪声背景下,基于线性收缩和随机矩阵理论的2D-RLSMW方法具有更好的优越性。
图5(a)、图5(b)和图6(a)、图6(b)分别为Student-T分布的非高斯噪声背景下收发阵元数大于快拍数时2D-RMW方法和2D-RLSMW方法的谱峰搜索示意图。将快拍数设为L=30,即收发阵元数大于快拍数的情况,目标位置信噪比为SNR=-5dB。从图5(a)和图5(b)可以看出,2D-RMW方法谱峰搜索具有较大的伪峰,无法正确的估计出收发角度;从图6(a)和图6(b)可以看出,2D-RLSMW方法通过谱峰搜索完全可以正确地估计出收发角,充分验证了2D-RLSMW方法的有效性。
图7为相同条件下进行40次蒙特卡罗仿真实验,利用2D-RMW和2D-RLSMW两种方法进行收发角估计的均方根误差(RMSE)随信噪比变化的曲线。从图7可以看出,随着信噪比的增加,2D-RMW方法的均方根误差在10°和100°之间浮动,该方法完全失效,无法估计出目标收发角,所估计的收发角度发生随机性变化;2D-RLSMW方法的均方根误差曲线随信噪比的增加呈下降趋势,而且在信噪比为SNR=5dB时,角度估计误差小于0.1°,说明在Student-T分布的非高斯噪声环境下,当收发阵元数大于快拍数时,基于线性收缩和随机矩阵理论的2D-RLSMW方法具有较大的优势。因此,本发明所提出的方法显著提高了大阵列MIMO雷达目标收发角联合估计的稳健性。
Claims (2)
1.一种稳健的大阵列MIMO雷达目标收发角联合估计方法,建立双基地MIMO雷达的***模型,M个发射阵元同时发射窄带信号照射K个远场目标,在接收端采用N个接收阵元,其中发射阵列和接收阵列均采用均匀线阵,阵元间距dt=dr=λ/2,发射角度为θ,接收角度为第k个目标的发射角度和接收角度分别为θk和k=1,…,K,L为快拍数,其特征在于包括下列步骤:
步骤一、MIMO雷达接收端N个接收阵元对多次快拍得到的回波数据信号进行收集和处理,第l个快拍期间回波数据可以写成一个N×Lt的矩阵X(l):
其中,表示M行1列的发射方向矢量,t表示发射端,表示N行1列的接收方向矢量,r表示接收端;S表示发射波形矩阵,这里sm(m=1,…,M)表示第m个发射阵元所发射的波形,其编码长度为Lt,且S为归一化正交矩阵,SSH=IM;表示在第l个快拍期间第k个目标的散射系数;
表示一个N行Lt列的非高斯噪声矩阵,服从Student-T分布;
步骤二、首先对N个接收阵元所接收的信号X(l)进行匹配滤波处理,得到:
其中,为N行M列的匹配滤波后的信号;
步骤三、对匹配滤波后的结果进行矢量化处理,即y(l)=vec(X(l)SH),在第l个快拍期间,经过预处理后的观测数据可表示为:
其中,为目标散射系数矩阵;是在第l个快拍期间处理后的噪声矢量,维数为MN行1列;是联合导向矢量矩阵,维数为MN行K列at(θk)和分别为发射和接收方向矢量,⊙是Khatri-Rao矩阵积,是Kronecker矩阵积;
步骤四、收集L个快拍期间的数据y(l)(l=1,…,L),可得到一个观测数据矩阵Y,它是一个维数为MN×L的随机矩阵,即:
其中,Y=[y(1),…,y(L)],B=[b(1),…,b(L)],N=[n(1),…,n(L)],这里表示包含收发角度的联合导向矩阵,B是服从均值为0,方差为1的目标散射系数矩阵,N为MN×L维的服从Student-T分布非高斯噪声矩阵,由高斯噪声矩阵的元素乘上服从伽马分布Γ(v,1/v)的因子构成,其中,v=0.02;
步骤五、接收信号的采样协方差矩阵可以表示为求解协方差矩阵的极大似然估计;
步骤六、当收发阵元数大于快拍数时,L<MN,上述极大似然估计方法因为秩缺失而导致无解,此时采用线性收缩方法使得定点迭代的过程能持续进行;收缩过程如下:
其中0<ρ<1为收缩系数,将定点迭代完成之后最终收敛的协方差矩阵记为收缩系数ρo由最小均方误差准则得到通过化简可以解得收缩系数为:
由于R未知,用替代R可以得到:
步骤七、在获得稳健协方差矩阵估计的基础上,利用2D-MW方法可以得到收发角联合估计的谱峰搜索函数;对采样协方差矩阵进行特征值分解得到:
其中,EEH=I,定义0<z<∞上的一个连续的函数f(z),其反函数为可以得到:
收发角联合估计的谱峰搜索函数可以表示为:
改变函数f(z)可以得到不同的谱峰搜索函数,令函数f(z)=ln(z),得到Midway(MW)方法的代价函数
步骤八、当收发阵元数与快拍数在同一数量级时,即MN→∞,L→∞,MN/L=c,c为一个常数,Y可以看作一个大维随机矩阵,传统方法性能急剧下降,不能正确估计出收发角度;因此,应用随机矩阵理论中的渐近谱理论和围线积分技术来产生一个稳健的谱峰搜索函数;方法如下:首先根据留数定理,可以表示为:
其中,ζ-是一条在上的顺时针方向曲线,并且该曲线包含R的所有特征值;
对曲线ζ-进行参数化处理,令其中,是一个包含所有特征值的矩形,表示该矩形的边界,wMN(z)与特征值有关,其关系式如下:
为了简化公式,定义两个变量bMN(z)和mMN(z)来分别描述特征值和特征向量:
得到参数化表示形式为:
通过变量bMN(z),得到wMN(z)的导数对上式进一步化简得到:
当MN→∞,L→∞,MN/L=c时,用bMN(z)和mMN(z)来描述特征值和特征向量不再准确,此时需要研究特征值和特征向量的渐近表现,在随机矩阵理论中,采样协方差矩阵特征值的瞬时分布函数的Stieltjes变换定义如下:
仅仅体现了采样协方差矩阵的特征值渐近特性,为了同时体现特征值和特征向量的渐近特性,定义如下函数:
此时,我们用和替换bMN(z)和mMN(z),得到:
其中,
为了求解定义一个关于的函数:
其中,x∈[0,c];
在x∈(0,c)区间上对求导:
利用留数定理可以得到:
其中,是下面等式的解:
可以由柯西留数定理直接得到,再此基础上对求积分,可以得到如下:
令x=c,经过化简可以得到ξmn(c):
其中,是下面等式的解:
至此,基于线性收缩和随机矩阵理论的MIMO雷达收发角联合估计谱峰搜索函数可以写为:
其中,如上所述;
步骤九、对进行谱峰搜索可以获得发射角θ和接收角
2.根据权利要求1所述的一种稳健的大阵列MIMO雷达目标收发角联合估计方法,其特征在于:所述步骤五的具体实现方法是:对接收信号y(l)进行归一化yi的概率密度函数为其中U(MN/2)/2πMN/2为单位圆上的均匀密度分布函数,对p(yi:R)求导并令导数等于零可以得到经过定点迭代可以得到极大似然估计定点迭代的初始矩阵可以是任意正定矩阵,令为单位矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710563441.6A CN107167785B (zh) | 2017-07-11 | 2017-07-11 | 一种稳健的大阵列mimo雷达目标收发角联合估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710563441.6A CN107167785B (zh) | 2017-07-11 | 2017-07-11 | 一种稳健的大阵列mimo雷达目标收发角联合估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107167785A CN107167785A (zh) | 2017-09-15 |
CN107167785B true CN107167785B (zh) | 2019-10-01 |
Family
ID=59823967
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710563441.6A Active CN107167785B (zh) | 2017-07-11 | 2017-07-11 | 一种稳健的大阵列mimo雷达目标收发角联合估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107167785B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109001690A (zh) * | 2018-07-05 | 2018-12-14 | 电子科技大学 | 基于馈电网络的时域空间域结合的雷达目标检测方法 |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108594194A (zh) * | 2018-03-19 | 2018-09-28 | 长江大学 | 基于四线性分解的双基地mimo雷达角度估算方法 |
CN108957387B (zh) * | 2018-05-21 | 2022-06-21 | 西安电子科技大学 | 一种卫星信号二维到达角估计方法及*** |
JP7166101B2 (ja) * | 2018-08-13 | 2022-11-07 | 株式会社東芝 | レーダシステム及びそのレーダ信号処理方法 |
CN111198351A (zh) * | 2018-11-19 | 2020-05-26 | 中移(杭州)信息技术有限公司 | 一种基于doa的定位方法、装置、设备及存储介质 |
JP7186596B2 (ja) * | 2018-12-07 | 2022-12-09 | 株式会社東芝 | アンテナシステム、レーダシステム及び信号処理方法 |
CN109870676B (zh) * | 2019-03-13 | 2021-01-12 | 电子科技大学 | 基于定位估计率的雷达通信一体化***性能界限测量方法 |
CN110146854B (zh) * | 2019-05-30 | 2022-12-23 | 西安电子科技大学 | 一种fda-mimo雷达稳健抗干扰方法 |
CN110596687B (zh) * | 2019-09-19 | 2022-09-27 | 吉林大学 | 一种基于黎曼流形的单基地mimo雷达目标检测方法 |
CN111427022A (zh) * | 2020-05-08 | 2020-07-17 | 北京理工大学重庆创新中心 | 一种基于极大似然估计的阵列雷达测角方法 |
CN111651719B (zh) * | 2020-06-04 | 2023-04-25 | 湖北工业大学 | 基于收缩估计协方差矩阵重构稳健自适应波束形成方法 |
CN112990256A (zh) * | 2020-12-29 | 2021-06-18 | 北京电子工程总体研究所 | 多传感器噪声计算方法及信息融合方法、存储介质及设备 |
CN112800599B (zh) * | 2021-01-15 | 2022-04-15 | 吉林大学 | 一种阵元失配情况下基于admm的无网格doa估计方法 |
CN113985394B (zh) * | 2021-10-05 | 2024-04-05 | 吉林大学 | 基于nomp的单基地mimo雷达远近场混合目标定位方法 |
CN114696876B (zh) * | 2022-03-30 | 2023-04-28 | 电子科技大学 | 一种针对MIMO Saleh-Valenzuela信道的收发角度联合估计方法 |
CN115047456B (zh) * | 2022-08-16 | 2022-10-25 | 中国人民解放军国防科技大学 | 无人驾驶车三维全景感知方法和雷达*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104360334A (zh) * | 2014-11-21 | 2015-02-18 | 吉林大学 | 相关噪声背景下基于特征值mimo雷达目标盲检测方法 |
CN106154242A (zh) * | 2016-06-21 | 2016-11-23 | 大连大学 | 基于分数低阶类相关熵的目标参数联合估计新算法 |
CN106199550A (zh) * | 2016-07-12 | 2016-12-07 | 北京理工大学 | 一种波离角和二维波达角自动配对联合估计方法 |
CN106680797A (zh) * | 2016-06-21 | 2017-05-17 | 大连大学 | 基于宽带模糊函数的目标参数估计新方法 |
-
2017
- 2017-07-11 CN CN201710563441.6A patent/CN107167785B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104360334A (zh) * | 2014-11-21 | 2015-02-18 | 吉林大学 | 相关噪声背景下基于特征值mimo雷达目标盲检测方法 |
CN106154242A (zh) * | 2016-06-21 | 2016-11-23 | 大连大学 | 基于分数低阶类相关熵的目标参数联合估计新算法 |
CN106680797A (zh) * | 2016-06-21 | 2017-05-17 | 大连大学 | 基于宽带模糊函数的目标参数估计新方法 |
CN106199550A (zh) * | 2016-07-12 | 2016-12-07 | 北京理工大学 | 一种波离角和二维波达角自动配对联合估计方法 |
Non-Patent Citations (1)
Title |
---|
《RANDOM MATRIX METHOD FOR JOINT DOD ESTIMATION FOR LARGE SCALE MIMO RADAR IN NON-GAUSSIAN NOISE》;Hong Jiang et.al;《College of Communication Engineering》;20161231;第3031-3035页 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109001690A (zh) * | 2018-07-05 | 2018-12-14 | 电子科技大学 | 基于馈电网络的时域空间域结合的雷达目标检测方法 |
CN109001690B (zh) * | 2018-07-05 | 2020-05-12 | 电子科技大学 | 基于馈电网络的时域空间域结合的雷达目标检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107167785A (zh) | 2017-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107167785B (zh) | 一种稳健的大阵列mimo雷达目标收发角联合估计方法 | |
CN108732549B (zh) | 一种基于协方差矩阵重构的阵元缺损mimo雷达doa估计方法 | |
CN109471082B (zh) | 基于信号子空间重构的阵元缺损mimo雷达角度估计方法 | |
CN105699945B (zh) | 频控阵mimo雷达***中的波形优化设计方法 | |
CN106707257B (zh) | 基于嵌套阵列的mimo雷达波达方向估计方法 | |
CN107290730B (zh) | 互耦条件下双基地mimo雷达角度估算方法 | |
CN103076596B (zh) | 基于先验信息的mimo雷达发射方向图设计方法 | |
CN106646387B (zh) | 基于发射波束域的mimo雷达抗有源干扰方法 | |
CN107064892B (zh) | 基于张量子空间和旋转不变的mimo雷达角度估计算法 | |
CN104251989B (zh) | 基于压缩空间谱的单基地mimo雷达目标波达方向估计方法 | |
CN107576940B (zh) | 一种低复杂度单基地mimo雷达非圆信号角度估计方法 | |
CN108717189B (zh) | 基于压缩感知理论的双基地mimo雷达成像方法 | |
CN104020469B (zh) | 一种mimo雷达距离-角度二维超分辨率成像算法 | |
CN108303683B (zh) | 单基地mimo雷达实值esprit非圆信号角度估计方法 | |
CN104360334B (zh) | 相关噪声背景下基于特征值mimo雷达目标盲检测方法 | |
CN106909779A (zh) | 基于分布式处理的mimo雷达克拉美罗界计算方法 | |
CN106646388B (zh) | 基于嵌套阵列的mimo雷达抗干扰方法 | |
CN111651719A (zh) | 基于收缩估计协方差矩阵重构稳健自适应波束形成方法 | |
CN108469601B (zh) | 改善mimo-stap最差检测性能的收发联合稳健优化方法 | |
CN109188386B (zh) | 基于改进二维esprit算法的mimo雷达高分辨率参数估计方法 | |
CN109245814B (zh) | 基于极大似然重采样的自适应波束形成方法 | |
CN103777197A (zh) | 单基地mimo雷达中降维传播算子的方位估计方法 | |
Feng et al. | Jointly iterative adaptive approach based space time adaptive processing using MIMO radar | |
CN109901103B (zh) | 基于非正交波形的mimo雷达doa估算方法及设备 | |
CN106125039B (zh) | 基于局域联合处理的改进空时自适应单脉冲测角方法 |
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 |