一种基于指向性模型数值求解的远场窄带信号来波方向估计方法

摘要

本发明公开了一种基于指向性模型数值求解的远场窄带信号来波方向估计方法，该方法包括：第一步：获取M通道阵元数据x

说明书

技术领域

本发明属于信号处理领域，尤其涉及一种基于指向性模型数值求解的远场窄带信号来波方向估计方法。

背景技术

目标来波方向估计是雷达、声纳、声学、语音和无线通信等阵列信号处理应用中的一个基本问题，尤其在水声和电子侦察处理中扮演了极其重要的角色。此外，在现代移动通信系统中，来波方向估计在提高无线系统性方面有着重要应用。

目前国内外学者提出了很多目标来波方向估计方法，如最小方差无失真响应、基于特征向量的来波方向估计算法(如MUSIC、Root-MUSIC、RV-Root-MUSIC)。这些算法通常能达到很高的精度。然而，这些基于峰值搜索的算法需要很大的计算复杂度。此外还存在着估计性能低、复杂度高和相位模糊等问题。同时，这些算法对信号失配非常敏感。

在波束形成方法中，常规波束形成是最常用的一种非自适应波束形成方法。常规波束形成使用一组固定的加权和时间延迟对来自阵列中传感器的信号进行组合叠加，这个叠加过程主要是使用关于传感器在空间中的位置和感兴趣的信号方向的信息。常规波束形成的优点是它对假定的和实际的信号波前之间的信号失配具有较强的鲁棒性，并且只需要少量样本即可稳健地估计目标来波方向，易于实现，因此在单目标场景中有着较好的性能。目前基于常规波束形成的来波方向估计方法主要有：峰值搜索算法和抛物线插值算法。

谱峰搜索法利用常规波束形成结果使角谱最大化，并返回M个最大值的位置作为M个目标的来波方向估计。谱峰搜索法对于白噪声中单个远场、窄带和电源的来波方向估计而言，是已知的最佳地最大似然估计方法。谱峰搜索法的来波方向估计精度取决于预引导波束扫描角的扫描间隔。为了获得更高的估计精度，需要更多的预引导波束扫描角。然而，更多的预引导波束扫描角将带来更多的计算负担。此外，受方位估计瑞利极限的限制，仅通过增加预引导波束扫描角的个数是无法无限提高来波方向估计精度的。

抛物线插值算法是基于对主瓣角谱的抛物线拟合，易于实现。对于大多数目标的来波方向估计，估计精度较高。然而，对于某些来波方向，角谱的主瓣与抛物线模型不匹配，在这种情况下，来波方向估计精度下降。

发明内容

技术问题：本发明提出了一种利用给定均匀线阵的理论波束图的模型求解算法。该算法是基于二分法的数值解，采用二分法数值求解信号来波方向相对偏差从而精确地估计远场窄带信号的来波方向。由于采用的是理论波束模式模型，因此不存在模型失配。与抛物线插值算法相比，该算法对信号来波方向相对偏差具有更强的鲁棒性，并且具有更好的估计性能。

技术方案：本发明的一种基于指向性模型数值求解远场窄带信号来波方向估计方法，包括以下步骤：

第一步：获取待处理的均匀线阵列M通道阵元远场窄带信号采样数据序列x

第二步：根据给定的阵列参数自动生成等余弦间隔的预引导波束扫描角θ

首先应根据给定的阵列参数计算总的预引导波束个数K，即

其中M为均匀线列阵的阵元个数，

然后，依据总的预引导波束个数K自动计算生成预引导波束扫描角θ

其中arccos[]代表取反余弦运算，k为预引导波束号。

第三步：对已知频率f

首先，计算待处理的均匀线列阵的M通道的远场窄带信号采样数据序列x

其中m为阵元通道号，m＝1,2,…,M，l为离散频率索引，j表示虚数单位，即

然后，依据各预引导波束扫描角θ

其中m为阵元通道号，m＝1,2,…,M，l为离散频率索引，k为预引导波束号，k＝1,2,…,K，d为均匀线列阵的间距，f

最后，计算已知频率f

其中||代表取模值运算，M为均匀线列阵通道数，l

l

其中round()代表四舍五入运算，f

第四步：搜索波束功率谱P(θ

根据下式搜索波束功率谱P(θ

其中

第五步：利用二分数值求解法求解信号来波方向相对偏差

首先，计算波束功率谱P(θ

然后，进行参数初始化：设置二分数值求解法的搜索上限a

这里

最后，利用二分数值求解法对来波方向相对偏差

(1)令第q次迭代得到的来波方向相对偏差值为

(2)计算理论波束功率谱的主瓣峰值和左右相邻的预引导波束所对应的波束功率谱值中较大的一个的比值，计算如下

其中,d为均匀线列阵的间距，f

(3)重置二分求解法的上下限，即：

如果

(4)判断是否满足q＝N

(5)计算求得来波方向相对偏差

第六步：利用

如果1＜k

如果k

如果k

其中

进一步的，本发明方法的第三步中，第m个通道的远场窄带信号采样数据序列x

进一步的，本发明方法的第五步中，通过利用二分数值求解法对来波方向相对偏差

本发明方法首先获取M通道阵元数据，然后根据给定的阵列参数生成等余弦间隔的预引导波束扫描角和总的预引导波束个数，接着对已知频率的远场窄带信号进行频域波束形成，并求得各预引导波束扫描角所对应的波束功率谱，搜索最大值对应的预引导波束号，然后利用二分数值求解法求解信号来波方向的相对偏差，最后利用该偏差估计远场窄带信号的来波方向。

有益效果：与现有技术相比，本发明的技术方案具有以下有益技术效果：

1.本发明的估计方法采用理论波束模式模型，因此不存在模型失配的问题。

2.本发明采用二分法数值求解信号来波方向相对偏差从而精确地估计远场窄带信号的来波方向，与抛物线插值算法相比，该算法对相对来波方向偏差具有更强的鲁棒性，并且具有更好的估计性能。

3.本发明仅利用角谱峰值附近的三个样本进行来波方向估计，因而具有与抛物线插值算法相当的计算量。

4.本发明可以推广应用于任意几何阵列的来波方向估计。

附图说明

图1所示为本发明的流程图。

图2所示为M通道阵元接收远场窄带信号示意图。

图3所示为实施例1中仿真M通道阵元的预引导波束扫描角所对应的波束功率谱。

图4所示为实施例2中仿真M通道阵元的预引导波束扫描角所对应的波束功率谱。

具体实施方式

下面结合实施例和说明书附图对本发明作进一步的说明。

如图1所示，一种基于指向性模型数值求解远场窄带信号来波方向估计方法，包括以下步骤：

第一步：获取待处理的均匀线阵列M通道阵元远场窄带信号采样数据序列x

第二步：根据给定的阵列参数自动生成等余弦间隔的预引导波束扫描角θ

首先应根据给定的阵列参数计算总的预引导波束个数K，即

其中M为均匀线列阵的阵元个数，

然后，依据总的预引导波束个数K自动计算生成预引导波束扫描角θ

其中arccos[]代表取反余弦运算，k为预引导波束序号。

第三步：对已知频率f

首先，计算待处理的均匀线列阵的M通道的远场窄带信号采样数据序列x

其中m为阵元通道号，m＝1,2,…,M，l为离散频率索引，j表示虚数单位，即

然后，依据各预引导波束扫描角θ

其中m为阵元通道号，m＝1,2,…,M，l为离散频率索引，k为预引导波束号，k＝1,2,…,K，d为均匀线列阵的间距，f

最后，计算已知频率f

其中||代表取模值运算，M为均匀线列阵通道数，l

l

其中round()代表四舍五入运算，f

第四步：搜索波束功率谱P(θ

根据下式搜索波束功率谱P(θ

其中

第五步：利用二分数值求解法求解信号来波方向相对偏差

首先，计算波束功率谱P(θ

然后，进行参数初始化：设置二分数值求解法的搜索上限a

这里

最后，利用二分数值求解法对来波方向相对偏差

(1)令第q次迭代得到的来波方向相对偏差值为

(2)计算理论波束功率谱的主瓣峰值和左右相邻的预引导波束所对应的波束功率谱值中较大的一个的比值，计算如下

其中d为均匀线列阵的间距，f

(3)重置二分求解法的上下限，即：

如果

(4)判断是否满足q＝N

(5)计算求得来波方向相对偏差

第六步：利用

如果1＜k

如果k

如果k

其中

本发明方法首先获取M通道阵元数据，然后根据给定的阵列参数生成等余弦间隔的预引导波束扫描角和总的预引导波束个数，接着对已知频率的远场窄带信号进行频域波束形成，并求得各预引导波束扫描角所对应的波束功率谱，搜索最大值对应的预引导波束号，然后利用二分数值求解法求解信号来波方向的相对偏差，最后利用该偏差估计远场窄带信号的来波方向。

本发明的实施例中，阵元接收远场窄带信号示意图如图2所示，仿真阵元接收远场窄带信号模型为：

其中A为均匀线性阵列M通道阵元接收的远场窄带信号的幅度，

实施例1

仿真信号参数分别设置为：信号幅度A＝1，初始相位

首先根据给定的阵列参数自动生成总的预引导波束个数K＝25，由此我们可以得到预引导波束扫描角θ

然后对已知频率f

接着通过搜索波束功率谱P(θ

然后设置二分数值求解法的搜索上限a

接着利用二分数值求解法对来波方向相对偏差

由于

估计来波方向相对误差为

实施例2

仿真信号参数分别设置为：信号幅度A＝1，初始相位

首先根据给定的阵列参数自动生成总的预引导波束个数K＝18，由此我们可以得到预引导波束扫描角θ

然后对已知频率f

接着通过搜索波束功率谱P(θ

然后设置二分数值求解法的搜索上限a

接着利用二分数值求解法对来波方向相对偏差

估计来波方向相对误差为

从实施例1和例2的结果可以看出，本发明估计方法可以获得良好的估计精度，而且计算简单，计算量小，适用于高精度快速估计多阵元远场窄带信号来波方向的场合。

上述实施例仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和等同替换，这些对本发明权利要求进行改进和等同替换后的技术方案，均落入本发明的保护范围。