[发明专利]增广互质阵列脉冲噪声环境下相干信源的DOA估计方法

说明书

本发明公开了一种增广互质阵列高脉冲环境下相干信源的波达方向估计方法，使用增广互质阵列天线接收信号，得到量测信息；根据量测信息，计算相应的分数低阶矩估计矩阵和相位分数低阶矩估计协方差矩阵；对所得到的估计协方差矩阵进行向量化处理，删除冗余行并截取连续阵元部分，得到阵元间距为λ/2的虚拟均匀线阵接收信号信息；对虚拟接收信号信息进行Toeplitz矩阵重构，得到了重构后的协方差矩阵；用SIC DOA估计算法计算重后的协方差矩阵DOA的精确估计。本发明将FLOM，PFLOM估计矩阵与脉冲噪声环境下的互质阵列相干信源DOA估计问题相结合，脉冲噪声低信噪比环境下，本发明能够获得更好的DOA估计性能。

技术领域

本发明属于波达方向(direction of arrival,DOA)估计，空中监视，雷达声呐定位等技术领域。

背景技术

传统的DOA估计方法，如MUSIC和ESPRIT，使用N个阵元估计最多N-1个信源，空间自由度较小。同时，为了避免角度模糊问题，传统阵列间距通常需要小于接收信号波长的一半，即d＜λ/2。但阵元距离过近会带来很强的互耦合效应，从而降低估计精度。因此，此类阵列的优化设计和性能分析一般都不容易。

为了解决上述问题，最近，利用两个平行均匀线阵，且阵元数互质，第一阵元重叠的平行均匀线阵组合得到的增广互质阵列得到了关注，如果两个均匀线阵的阵元个数分别为2M和N，其中M和N互质，则2M+N-1阵元的增广互质阵列可以获得2M(N+1)-1的虚拟线阵连续自由度(degree of freedom,DOF)，而阵元个数相同的均匀线阵(uniform lineararray,ULA)只能获得2M+N的DOF。因此，增广互质阵列结构大大提高了阵列可探测源的数量。该阵列结构突破了传统天线阵阵元间距半波长的限制，使得天线孔径得到极大的扩展，能够获得角度估计性能的提升。同时，子阵单元间的间距为Mλ/2和Nλ/2，远大于半波长，有效地减弱了单元间的相互耦合效应。

此外，互质阵列中的大多数DOA估计方法均假设环境噪声为高斯分布。然而，实际中的噪声往往表现出非高斯特性，有时还伴有较强的脉冲。例如，汽车点火装置、微波炉和其他类型的自然或人为信号源可能会表现出高脉冲特性。近来，研究表明α稳定分布是一个合适的噪声模型来描述这种类型的噪声，它可以用S(α,β,γ,θ)表示，包括一系列分布，如高斯分布，柯西分布，莱维分布。参数的定义如下：α∈(0,2]是决定分布形状的特征指数；β∈[-1,1]是对称参数，指定分布是否为右偏β＞0或左偏β＜0；γ∈(0,+∞)是分散参数，类似于正态分布中的方差；δ∈R是位置参数作为正态分布的均值。通常选择β＝δ＝0来生成对称α稳定(SαS)过程。如果特征指数α＝2，稳定过程S(α,β,γ,θ)将是一个高斯过程，如果α＝1则是一个柯西过程。当选择不同的特征指数参数时，可以认为它是代表各种脉冲噪声的最大电位分布。

另外，由于传播环境的复杂性，入射到阵列的信号中是有相干信号源存在的。由于信号阵列会接收到不同方向上的相干信号，而相干信号会导致信源协方差矩阵的秩亏缺，从而使得信号特征向量发散到噪声子空间去。

发明内容

发明目的：为了解决现有技术存在的问题，本发明提供了一种增广互质阵列脉冲噪声环境下相干信源的DOA估计方法。

技术方案：本发明提供了一种增广互质阵列脉冲噪声环境下相干信源的DOA估计方法，包括如下步骤：

S1：使用增广互质阵列结构的阵列天线接收信号，得到量测信息Z；

S2：根据所述量测信息，计算分数低阶矩估计协方差矩阵R_FLOM和相位分数低阶矩估计协方差矩阵R_PFLOM；

S3：将R_FLOM和R_PFLOM分别进行向量化处理，并去掉获得的向量中的冗余行，得到虚拟阵列接收信号和