[发明专利]基于最大似然的嵌套阵稀疏表示波达方向估计方法

说明书

本发明公开了一种基于最大似然的嵌套阵稀疏表示波达方向估计方法，首先根据目标到嵌套阵的波达方向以及嵌套阵的布阵结构，计算接收阵列的回波信号、进而计算其协方差矩阵；根据稀疏表示理论，在角度维划分整个空域，得到角度集合；矢量化协方差矩阵，在角度集合上稀疏展开，得到嵌套阵波达方向估计稀疏模型；构造块对角矩阵，剔除稀疏模型中的噪声项，得到去噪后嵌套阵波达方向估计稀疏模型；再结合协方差矩阵计算噪声白化矩阵，结合去噪后嵌套阵波达方向估计稀疏模型，计算噪声白化后的数据模型，进而计算栅格最大似然估计，最后得到波达方向最大似然估计值。本发明以实现嵌套阵在低信噪比和少快拍下的波达方向估计性能。

技术领域

本发明属于雷达技术领域，具体涉及一种基于最大似然的嵌套阵稀疏表示波达方向估计方法。

背景技术

波达方向估计技术是实现目标定位不可或缺的技术手段，通过波达方向估计结果可得目标在空域的方位信息。传统典型的方法有多重信号分类方法，旋转不变子空间方法，这两种方法突破瑞利限，可实现目标的超分辨，但在低信噪比或者少快拍环境下性能严重下降。

最大似然估计假设信号源为一具有已知分布的随机过程，利用已知的样本结果信息，反推最具有可能导致这些样本结果出现的待估计参数值。在低信噪比或者相干信号源情况下仍然可以获得较好的估计性能，但是求解比较困难，需要在多维参数空间上寻优，这将导致算法的运算量很大，并且算法的性能依赖于初始值的选取。

嵌套阵是由多个均匀线阵嵌套而成，通过非均匀布阵减少空间采样的冗余性，可提高线性阵列的自由度，并且阵元位置的闭式解和嵌套阵所能提高的自由度均可通过阵元总数计算。

稀疏表示理论的基本思想是用冗余函数集合组成的原子字典代替基函数集合，进而可以将信号表示为原子字典中少数原子列向量的线性组合。当空域分布少量点目标时，对整个空域角度而言目标具有稀疏性，因此可以将稀疏表示理论应用到波达方向估计中。但需假设目标正好落在划分的角度栅格上，这将不可避免的存在模型失配问题。

针对稀疏表示波达方向估计方法中模型失配问题，现有方法主要分为两大类，一类是多次重构法，算法的计算量随着栅格个数增加按指数级增长，另一类是贝叶斯估计法，虽然降低了字典矩阵栅格失配问题，但仍然存在失配量。

发明内容

本发明的目的是提供一种基于最大似然的嵌套阵稀疏表示波达方向估计方法，以实现嵌套阵在低信噪比和少快拍下的波达方向估计性能。

本发明所采用的技术方案是，基于最大似然的嵌套阵稀疏表示波达方向估计方法，具体按照以下步骤实施：

步骤1、根据目标到嵌套阵雷达系统接收阵列的波达方向以及嵌套阵的布阵结构，计算接收阵列的回波信号；

步骤2、根据步骤1得到的接收阵列的回波信号，计算嵌套阵接收数据的协方差矩阵；

步骤3、根据稀疏表示理论，在角度维划分整个空域，得到角度集合；矢量化步骤2得到的嵌套阵接收数据协方差矩阵，并在角度集合上稀疏展开，得到嵌套阵波达方向估计稀疏模型；

步骤4、构造块对角矩阵，剔除步骤3稀疏模型中的噪声项，得到去噪后嵌套阵波达方向估计稀疏模型；再结合步骤2得到的协方差矩阵计算噪声白化矩阵，结合得到的去噪后嵌套阵波达方向估计稀疏模型，计算噪声白化后的稀疏模型；

步骤5、根据步骤4建立的噪声白化后的稀疏模型，计算目标波达方向。

本发明的特点还在于，

步骤1中接收阵列的回波信号y(t)计算如下：

y(t)＝A(θ)s(t)+n(t)，