本发明涉及一种雷达信号处理技术,具体的说涉及一种基于music算法的双基地mimo雷达角度估算方法。
背景技术:
多输入多输出(multiple-inputmultiple-output,mimo)雷达是未来雷达发展的主流趋势,从其概念被提出之日便引起国内外学者的广泛关注。与现有相控阵雷达系统不同,mimo雷达在发射端发射多个相互正交的波形,在接收端利用匹配滤波器分离出各个通道的信息,利用分集增益大大改善目标探测性能。相比相控阵雷达系统,mimo雷达在分辨率、抗衰落性、可辨识性以及抑制噪声等方面具有潜在的优势。
角度估计是mimo雷达目标定位的关键一环,至今已涌现出大量优秀的角度估计算法。典型的算法有capon算法、多重谱峰分类(multiplesignalclassification,music)算法、基于旋转不变技术的参数估计(estimationmethodofsignalparametersviarotational,esprit)算法、传播算子(propagatormethod,pm)算法、高阶子空间分解(higherordersingularvaluedecomposition,hosvd)算法、平行因子(parallelfactor,parafac)算法、基于稀疏表示的估计算法等。然而上述算法的优异性能均是在假设接收噪声为均匀高斯白噪声的条件下获得的。由于实际工程中探测背景的复杂性,噪声往往呈现非高斯特性。空域色噪声是阵列mimo雷达中最常见的一类非高斯噪声,在存在空域色噪声的条件下,信号和噪声(或者信号子空间和噪声子空间)无法有效分离,从而引起角度估计性能的下降。特别是在低信噪比条件下,现有算法可能会至完全失效。mimo雷达中的色噪声问题已经引起部分学者的注意,目前已经有部分学者提出了一些色噪声抑制算法。按照抑噪机理的不同,可以将现有抑噪算法划分为两大类——空域协方差算法和时域协方差算法。前者主要是利用不同匹配滤波器输出的色噪声不相关的性质,后者主要利用了不同快拍的空域色噪声不相关特性。前者的典型代表有jin等人的算法(jinm,liaog,lij.jointdodanddoaestimationforbistaticmimoradar[j].signalprocessing,2009,89(2):244-251.)chen等人的算法(anewmethodforjointdodanddoaestimationinbistaticmimoradar)、jiang等人的算法(jiangh,zhangjk,wongkm.jointdodanddoaestimationforbistaticmimoradarinunknowncorrelatednoise[j].ieeetransactionsonvehiculartechnology,2015,64(11):5113-5125.)、wang等人的算法(wangx,wangw,lix,etal.atensor-basedsubspaceapproachforbistaticmimoradarinspatialcolorednoise[j].sensors,2014,14(3):3897-3907.)。其中jin等人的算法仅适用于三个发射阵元的mimo雷达,chen等人的算法是对jin等人方法的推广,该算法适用于多于三个阵元的mimo雷达配置。其将发射阵元分为两个子阵,利用两个子阵数据的互协方差抑制色噪声;jiang等人将收发阵元均分成两个子阵,该算法需要估计四维角度;wang等人的算法能在抑制色噪声的同时利用阵列数据的多维结构,从而提升角度估计精度。空域协方差算法的最大缺陷是存在阵列虚拟孔径的损失,从而降低了参数估计的精度。时域协方差法能有效的避免该缺陷,这类算法主要有符渭波等人的算法(符渭波,苏涛,赵永波,等.空间色噪声环境下基于时空结构的双基地mimo雷达角度和多普勒频率联合估计方法[j].电子与信息学报,2011,33(7):1649-1654.)和wen等人的算法(wenfq,xiongxd,suj,etal.angleestimationforbistaticmimoradarinthepresenceofspatialcolorednoise[j].signalprocessing,2017,134,261-267.)。其中,wen等人的算法能有效利用阵列数据的多维结构,获得比符渭波等人算法更好的角度数估计的精度。然而,现有算法并不适用于非均匀阵列配置,其适用范围存在一定的限制。此外,现有抑噪算法均忽略了噪声结构的先验信息,因而角度估计的精度有待于进一步提高。
技术实现要素:
鉴于以上原因,有必要提供一种能够对协方差矩阵具有厄米特对称toeplitz类型的色噪声的具有明显的抑制效果,提高双基地mimo雷达角度估算的精度;且适用于均匀阵列和非均匀阵列配置的基于music算法的双基地mimo雷达角度估算方法。
本发明提供一种基于music算法的双基地mimo雷达角度估算方法,所述基于music算法的双基地mimo雷达角度估算方法包括如下步骤:
s1、构建接收阵列数据在匹配滤波后的差分协方差矩阵;
s2、利用特征值分解的方法对协方差矩阵进行特征分解,获得噪声子空间的估计;
s3、利用谱峰搜索中的极小值获得目标角度的估计;
s4、通过配对计算确定真实dod值及所估计角度的配对。
本发明所述基于music算法的双基地mimo雷达角度估算方法,其通过构建接收阵列数据在匹配滤波后的差分协方差矩阵,利用特征值分解的方法对协方差矩阵进行特征分解,获得噪声子空间的估计;并利用谱峰搜索中的极小值即可获得目标角度的估计;最后通过配对计算确定真实dod值及所估计角度的配对。本发明所述估算方法不仅对协方差矩阵具有厄米特对称toeplitz类型的色噪声的具有明显的抑制效果,提高了双基地mimo雷达角度估算的精度;还对mimo雷达的阵列流形要求低,适用于均匀阵列和非均匀阵列配置;且相对二维music算法,本发明所述基于music算法的双基地mimo雷达角度估算方法的计算复杂度低。
附图说明
图1是双基地mimo雷达角度估计示意图;
图2是本发明所述基于music算法的双基地mimo雷达角度估算方法在空域色噪声背景下对doa进行谱峰搜索的结果;
图3是本发明所述基于music算法的双基地mimo雷达角度估算方法在空域色噪声背景下对dod进行谱峰搜索的结果;
图4是本发明所述基于music算法的双基地mimo雷达角度估算方法与现有其他算法在空域色噪声背景下的rmse比较。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明,应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明提供一种基于music算法的双基地mimo雷达角度估算方法,所述基于music算法的双基地mimo雷达角度估算方法包括如下步骤:
s1、构建接收阵列数据在匹配滤波后的差分协方差矩阵;
s2、利用子特征值分解的方法对协方差矩阵进行特征分解,获得噪声子空间的估计;
s3、利用谱峰搜索中的极小值即可获得目标角度的估计;
s4、通过配对计算确定真实dod值及所估计角度的配对。
具体的,如图1所示,假设天线系统由m个发射阵元和n个接收阵元构成,二者均为线性阵列,且收发阵元间的间距均为λ/2,λ为发射信号波长。假设第m(m=1,…,m)个发射阵元在一个脉冲重复周期内发射码长为q的归一化基带编码信号
上式中,
容易得知,如果噪声是平稳的,则c具有toeplitz的结构先验,对每个接收阵元均用sm/q进行匹配滤波处理,并将阵列数据表述成矩阵形式。则匹配滤波输出结果可以被表述成
y=[at⊙ar]b+n=ab+n表达式2
上式中,
在传统子空间算法中,需要对接收阵列信号的协方差矩阵ry进行分解。在本发明的信号模型中,阵列信号协方差矩阵为
ry=e{yyh}=arbah+rn表达式3
其中,rb=e{bbh}为目标特性协方差矩阵,在不考虑多径效应的条件下rb近似为一个对角矩阵rn=e{nnh}为噪声协方差矩阵。在样本数量有限的条件下阵列的协方差矩阵ry的估计为
在传统高斯白噪声假设下,rn=σ2imn,σ2为噪声功率,从而对
上述公式的第一个等号处用到矢量化的性质vec(abc)=(cta)vec(b),第三个等号处用到了克罗内科积的性质
因此,定义如下差分协方差矩阵
根据表达式5,可知δry中的色噪声已经被有效的抑制。由于
上式中
其中
显然,上式中
根据对表达式9的分析,可知δry具有共轭对称特性。因此δry可以表示为如下特征分解的形式
其中,
通过谱峰搜索上式中的极小值即可获得目标角度的估计。然而,上述搜索过程需要二维联合谱峰搜索,计算复杂度非常高。为降低谱峰搜索的复杂度,可利用克罗内克积的性质对上式进行如下简化
由于
由于doa与dod是分别获取的,因此该算法还需要额外的配对计算。通过组合计算表达式11中的值,前2k个最小值的位置即可配对所估计的角度。此外,由于dod的峰值是在
如果上述猜测是正确的,则g的上半对角线元素全为正,下半对角线元素全为负。如果真实的角度不是上述组合,则上述现象不存在,但通过对g的分析仍然可以正确的确定正确角度的组合,值得注意的是,本发明所述估算方法对线性阵列并无进一步的位置要求,因此本发明所述估算方法同样适合非均匀阵列。
进一步的,技术人员针对本发明提出的空域色噪声背景下双基地mimo雷达角度估计方法进行了大量的仿真实验。假设k=3个目标处于远场,其doa和dod分别为
附图2和附图3分别为本发明所提算在snr=-15db时,角度估计效果与直接对接收信号进行子空间分解后再利用降维music(rd-music)算法进行角度估计的效果图比较。由仿真结果可以看出,受色噪声的影响,rd-music算法此时不能有效的进行角度估计,而本发明所提的算法在dod估计时尽管存在角度模糊,但其能有效的估计出目标角度。
为比较本发明所述估算方法同现有算法估计精度的比较,对算法进行500次蒙特卡洛仿真,角度估计的精度用均方根误差(rootmeansquarederror,rmse),其中rmse定义为
式中
本发明所述基于music算法的双基地mimo雷达角度估算方法,其通过构建接收阵列数据在匹配滤波后的差分协方差矩阵;利用子特征值分解的方法对协方差矩阵进行特征分解,获得噪声子空间的估计;并利用谱峰搜索中的极小值即可获得目标角度的估计;最后通过配对计算确定真实dod值及所估计角度的配对。本发明所述估算方法不仅对协方差矩阵具有厄米特对称toeplitz类型的色噪声的具有明显的抑制效果,提高了双基地mimo雷达角度估算的精度;还对mimo雷达的阵列流形要求低,适用于均匀阵列和非均匀阵列配置;且相对二维music算法,本发明所述基于music算法的双基地mimo雷达角度估算方法的计算复杂度低。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。