基于矩阵重构的Coprime阵稳健自适应波束形成方法

摘要

本发明提供一种基于矩阵重构的Coprime阵稳健自适应波束形成方法，充分利用了Coprime阵稀疏阵列在虚拟孔径大的特性，在虚拟域估计出了信号和干扰的角度，并利用噪声功率和协方差矩阵的对应关系估计出了信号和干扰的功率，实现了INCM的重构。当Coprime阵配置为：两个子阵阵元数分别为2M1、M2，阵元间距分别为M2d、M1d、d为最小阵元间距，应用SERB算法时，其抑制干扰的自由度为min{2M1‑1,M2‑1}，而所提算法在Coprime阵中抑制干扰的自由度为2M1+M2‑2，明显高于SERB算法抑制干扰的自由度。

说明书

技术领域

本发明涉及阵列雷达技术领域，具体涉及基于矩阵重构的Coprime阵稳健自适应波束形 成方法。

背景技术

Coprime阵相较于传统均匀阵列采用了稀疏布阵的方式，虽然具有大孔径，高自由度的特 性，但其阵列结构的变化使得传统信号处理算法不能发挥出Coprime阵的优势，因此，诸多基 于Coprime阵的阵列信号处理算法相继被提出。相较于Coprime阵中波达方向估计的研究， Coprime阵中的波束形成问题研究并不深入。

为了提高Coprime阵波束形成的性能，有人提出了一种提出了基于子阵估计和干扰噪声协 方差矩阵重构的波束形成SERB算法。SERB算法利用Coprime阵的双子阵无模糊地估计出了入射 信号和干扰的功率，精确重构出了干扰噪声协方差矩阵，提高了自适应波束形成的稳健性。

由于SERB算法对信号和干扰角度及功率的估计是采用双子阵分别估计后得到的，虽然能 够估计出的干扰个数受到子阵阵元数的限制，且明显低于整个阵列的阵元数，限制了算法抑 制多干扰的能力。

发明内容

本发明的目的在于提供一种基于矩阵重构的Coprime阵稳健自适应波束形成方法，解决 以往算法中使用双子阵来估计信号功率，使干扰估计的自由度和抑制多干扰的能力均不足的 问题。

为解决上述的技术问题，本发明采用以下技术方案：

一种基于矩阵重构的Coprime阵稳健自适应波束形成方法，包括步骤：

第一步，对Coprime阵的阵元位置集合进行定义，将远场窄带信号入射到定义后的Coprime阵上，并根据Coprime阵接收的数据得出协方差矩阵；

第二步，接收数据的协方差矩阵等效构建一个虚拟阵列，在Coprime阵构成的虚拟阵列 中部分阵元构成虚拟ULA；

第三步，根据空间平滑原理，对虚拟ULA划分成U个阵元数为U的ULA，并计算得到空间平滑矩阵R

第四步，构建并定义Toeplitz矩阵R

第五步，将获得的矩阵R

第六步，根据重构实阵列的INCM得出Coprime阵的波束形成权矢量，并形成稳健自适应 波束。

进一步的技术方案是，所述Coprime阵的阵元位置集合定义为：

当L个远场窄带信号入射到Coprime阵上时，接收到的信号表示为：

此时接收数据的协方差矩阵表示为：

当

其中，[·]

更进一步的技术方案是，所述第m个阵元和第n个阵元之间的互相关统计量等效为物理 阵元位置0处和物理阵元位置m-n处形成的相位差，即接收数据的协方差矩阵等效为一个虚 拟阵列，且其上虚拟阵元的位置满足：

其中，u

虚拟阵列中部分阵元构成了ULA，且构成的虚拟ULA中虚拟阵元的分布可以表示为：

定义

其中，

根据虚拟ULA中每个阵元的虚拟接收数据由虚拟阵列接收数据z中对应虚拟ULA阵元 位置的元素获得，可表示为：

其中，

其中，U(u)定义为所有使得虚拟阵元位置为u的实阵元位置集合，表示为：

U(u)＝{(m

更进一步的技术方案是，所述虚拟ULA中第t(t＝1,2,L,U)个虚拟子ULA的阵元分布

示为：

根据

建立模型

更进一步的技术方案是，所述构建的Toeplitz矩阵为：

并定义Toeplitz矩阵为：

更进一步的技术方案是，所述重构实阵列的INCM为：

更进一步的技术方案是，所述Coprime阵S的波束形成权矢量为

与现有技术相比，本发明的有益效果是：

1、利用构造的Toeplitz矩阵、空间平滑矩阵与协方差矩阵的对应关系，利用Toeplitz 矩阵的特征值直接估计出了干扰功率，简化了估计过程，提高了估计的自由度。

2、充分利用了Coprime阵稀疏阵列在虚拟孔径大的特性，在虚拟域估计出了信号和干扰 的角度，并利用噪声功率和协方差矩阵的对应关系估计出了信号和干扰的功率，实现了INCM 的重构，提高了波束形成的稳健性。

附图说明

图1为本发明基于矩阵重构的Coprime阵稳健自适应波束形成方法的示意图。

图2为输出SINR和SNR的关系，INR＝30dB，K＝200时，高自由度稀疏阵列自适应 波束形成算法在信号角度已知时的性能示意图。

图3为输出SINR和快拍数的关系，SNR＝10dB，INR＝30dB时，高自由度稀疏阵列 自适应波束形成算法在信号角度已知时的性能示意图。

图4为输出SINR和SNR的关系，INR＝30dB，K＝200时，高自由度稀疏阵列自适应 波束形成算法在信号角度存在随机误差时的性能示意图。

图5为输出SINR和快拍数的关系，SNR＝10dB，INR＝30dB时，高自由度稀疏阵列 自适应波束形成算法在信号角度存在随机误差时的性能示意图。

图6为输出SINR和SNR的关系，INR＝30dB，K＝200时，高自由度稀疏阵列自适应 波束形成算法在存在不相干近场散射时的性能示意图。

图7为输出SINR和快拍数的关系，SNR＝10dB，INR＝30dB时，高自由度稀疏阵列 自适应波束形成算法在存在不相干近场散射时的性能示意图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发 明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用 于限定本发明。

实施例：

图1-7示出了本发明基于矩阵重构的Coprime阵稳健自适应波束形成方法的一个较佳实 施方式，本实施例中的Coprime阵稳健自适应波束形成方法具体包括步骤：

第一步，对Coprime阵的阵元位置集合进行定义，将远场窄带信号入射到定义后的Coprime阵上，并根据Coprime阵接收的数据得出协方差矩阵；

第二步，接收数据的协方差矩阵等效构建一个虚拟阵列，在Coprime阵构成的虚拟阵列 中部分阵元构成虚拟ULA；

第三步，根据空间平滑原理，对虚拟ULA划分成U个阵元数为U的ULA，并计算得到空间平滑矩阵R

第四步，构建并定义Toeplitz矩阵R

第五步，将获得的矩阵R

第六步，根据重构实阵列的INCM得出Coprime阵的波束形成权矢量，并形成稳健自适 应波束。

在本发明中利用构造的Toeplitz矩阵、空间平滑矩阵与协方差矩阵的对应关系，利用 Toeplitz矩阵的特征值直接估计出了干扰功率，简化了估计过程，提高了估计的自由度，同 时充分利用了Coprime阵稀疏阵列在虚拟孔径大的特性，在虚拟域估计出了信号和干扰的角 度，并利用噪声功率和协方差矩阵的对应关系估计出了信号和干扰的功率，实现了INCM的重 构，提高了波束形成的稳健性。

更进一步的技术方案是：Coprime阵是由两个阵元数互质的均匀线阵。假设子阵1和子阵 2的阵元数分别是2M

当L个远场窄带信号入射到Coprime阵上时，接收到的信号表示为；

此时接收数据的协方差矩阵表示为：

对于

其中，[·]

第m个阵元和第n个阵元之间的互相关统计量可以等效为物理阵元位置0处和物理阵元 位置m-n处形成的相位差，换句话说，接收数据的协方差矩阵可以等效出一个虚拟阵列，且 其上虚拟阵元的位置满足：

其中，u

更进一步的技术方案是：在Coprime阵构成的虚拟阵列中部分阵元构成了ULA，且构成 的虚拟ULA中虚拟阵元的分布可以表示为：

其中，为了表示方便，u中忽略了最小阵元间距。此时定义

根据获得差分阵的过程可知，虚拟ULA中每个阵元的虚拟接收数据可由虚拟阵列接收数 据z中对应虚拟ULA阵元位置的元素获得，可表示为：

其中，

其中，U(u)定义为所有使得虚拟阵元位置为u的实阵元位置集合，定义为：

U(u)＝{(m

为了更方便理解，这里给出一个实例：当阵元位置分布

虚拟ULA上的接收数据建立模型：

其中，

据空间平滑原理，可对虚拟ULA划分成U个阵元数为U的ULA，其中第t(t＝1,2,L,U)个 虚拟子ULA的阵元分布可表示为：

据此，类比

根据

其中，

据此，可得到空间平滑矩阵R

其中，

由式

根据上面的分析可以发现计算空间平滑矩阵R

并定义Toeplitz矩阵：

据此，矩阵R

此时，有

发现构造的矩阵R

而对于矩阵R

其中：

据此可得：

其中，Z＝[z

由于干扰功率远大于噪声功率，据此，R

由上可知空间平滑矩阵中对应的干扰功率以及噪声功率近似是实际干扰功率和噪声功率 的平方。矩阵

其中，

直接将获得的矩阵R

其中，

根据重构实阵列的INCM得出Coprime阵的波束形成权矢量，并形成稳健自适应波束。

尽管这里参照本发明的多个解释性实施例对本发明进行了描述，但是，应该理解，本领 域技术人员可以设计出很多其他的修改和实施方式，这些修改和实施方式将落在本申请公开 的原则范围和精神之内。更具体地说，在本申请公开、附图和权利要求的范围内，可以对主 题组合布局的组成部件和/或布局进行多种变型和改进。除了对组成部件和/或布局进行的变 形和改进外，对于本领域技术人员来说，其他的用途也将是明显的。