一种基于分子阵预处理的MUSIC波束形成方法与流程

本发明涉及声纳信号处理领域，特别涉及一种基于分子阵预处理的music波束形成方法。

背景技术：

水下目标检测与估计是阵列信号处理的一个重要分支。波束形成作为阵列信号处理中核心算法，输出波束中的背景噪声和旁瓣级一直是其设计中需要考虑的问题。低背景噪声和旁瓣级可以有效降低对位于强目标旁瓣区域内弱目标检测的漏报概率。

为了对波束形成输出背景噪声和旁瓣级实现控制，许多学者从不同方法对降低波束形成旁瓣级进行了深入研究，并取得了一定的研究成果，提出了很多方法，主要为chebyshev滤波方法，“凹槽噪声场”方法，静态波束图数字综合方法，反复迭代方法，多线性约束方法，非线性优化方法，凸优化(convexoptimization)方法，半无限二次规划(semi-infinitequadraticprogramming)方法，二阶锥(second-ordercone)约束方法，中心矩方法，虚拟干扰源构造能量聚焦矩阵方法，稀疏约束方法。在以上方法中，chebyshev滤波方法以其简单方便常被应用实际工程中，但存在旁瓣级设置和主瓣宽度控制折中选择问题。

目前的方法都未能解决针对music波束形成中背景噪声和旁瓣级对检测弱目标性能影响的问题。

技术实现要素：

本发明的目的在于针对music波束形成中背景噪声和旁瓣级对检测弱目标性能影响的问题，根据线列阵接收数据中信号和噪声相关性的差异，提出一种基于分子阵预处理的music波束形成方法。该方法首先需要将线列阵中2n-1个阵元接收数据通过分子阵预处理转变为n个阵元数据；然后再采用music波束形成思想对该n个阵元数据进行处理，可得到各方位处空间谱。由于该方法对线列阵接收数据进行了分子阵预处理，有效提高线列阵接收数据协方差矩阵中信号含有量和信噪比，降低了背景噪声和旁瓣级对music波束形成检测弱目标带来的影响，提高了music波束形成对弱目标检测的检测效果。

为了实现上述目的，本发明提出了一种基于分子阵预处理的music波束形成方法，所述方法包括：

根据线列阵接收数据中信号和噪声相关性的差异，将线列阵中2n-1个阵元接收数据通过分子阵预处理转变为n个阵元数据；然后再对n个阵元数据进行处理，得到music波束输出。

作为上述方法的一种改进，所述方法具体包括：

步骤1)对线列阵2n-1个阵元接收数据按下式进行分组处理：

式中，xn(fl)表示第n个阵元拾取的频率fl数据，表示为：

式中：s(fl)为目标辐射信号，nn(fl)为第n个阵元拾取的加性高斯白噪声数据，λ＝fl/c为波长，d为线列阵相邻阵元间距，θ0为目标相对线列阵入射角度，c为声速入射；

步骤2)在搜索角度θ,θ＝1,2,…180处，按下式对各组数据进行相移预处理：

步骤3)对各组数据预处理结果进行相加，得到一组新数据y(fl)为：

步骤4)求取y(fl)的协方差矩阵ry(fl)＝e[y(fl)^hy(fl)]并进行特征分解：

式中，e[·]为期望函数，λ＝diag(λ1,λ2,…,λn)为新数据协方差矩阵ry(fl)的特征值，u＝[u1,u2,…,un]为新数据的协方差矩阵ry(fl)的特征向量，uaps为新数据y(fl)的信号特征向量，uapv为新数据y(fl)的噪声特征向量；

步骤5)根据噪声特征向量和目标波达方向导向权向量的正交性，获得基于分子阵预处理的music波束形成输出的来波方向波束为：

式中，为导向权向量，τn＝(n-1)dcos(θ)/c,1≤n≤n。

作为上述方法的一种改进，所述方法还包括：

步骤6)按下式求取波束输出的宽带空间谱：

其中，l为频带数。

本发明的优势在于：

本发明方法对线列阵接收数据进行了分子阵预处理，有效提高线列阵接收数据协方差矩阵中信号含有量和信噪比，降低了背景噪声和旁瓣级对music波束形成检测弱目标带来的影响，提高了music波束形成对弱目标检测的检测效果。

附图说明

图1为本发明的拖线阵声纳结构示意图；

图2为本发明的方法与现有方法对31元线列阵波束形成的比较结果图；

图3为本发明的方法与现有方法对63元线列阵波束形成的比较结果图；

图4为本发明的方法与现有方法对31元线列阵波束形成的比较结果图(强、弱目标辐射信号谱级比为30db)；

图5为本发明的方法与现有方法对63元线列阵波束形成的比较结果图(强、弱目标辐射信号谱级比为30db)；

图6为本发明的方法与现有方法对63元线列阵波束形成的比较结果图(强、弱目标辐射信号谱级比为80db)。

具体实施方式

现结合附图对本发明作进一步的描述。

在对本发明的方法做详细说明前，首先对本发明的方法所适用接收阵加以描述。图1为一拖线阵声纳结构示意图，该拖线阵声纳包括6个部分，显控与信号处理机1、甲板缆2、绞车3、导缆架4、拖缆5、接收线阵6。其中的接收线阵6通过拖缆5与位于绞车3上的甲板缆2连接，所述拖缆5还安装在导缆架4上；接收线阵6所接收的信号传输给显控与信号处理机1。

下面对本发明方法做进一步的说明。

music波束形成数学表达式

对于间距为d的2n-1元等间隔水平线阵，有1个目标从θ0入射，则第n个阵元拾取的频率fl数据xn(fl)可表示为：

式中：s(fl)为目标辐射信号，nn(fl)为第n个阵元拾取的加性高斯白噪声数据，c为声速，λ＝fl/c为波长。

对线列阵各阵元接收数据构造数据矩阵，可表示为

x(fl)＝[x1(fl),x2(fl),…,x2n-1(fl)]^t(2)

在噪声与信号不相关的假设下，对协方差矩阵rx(fl)＝e[x(fl)x(fl)^h]作特征分解得到信号特征向量us和噪声特征向量uv，依据噪声向量与导向权向量的正交性，可获得来波方向输出波束为

式中，为导向权向量，τn＝(n-1)dcos(θ)/c，θ为搜索角度，c为声速。

基于分子阵预处理的music波束形成方法

为了进一步降低在非目标波达方向上，music波束形成输出波束在不同搜索角度θ形成的极大值，减小其对弱目标检测形成影响。本发明将根据线列阵接收数据协方差矩阵形成过程中信号和噪声相关性的差异，通过对线列阵接收数据进行分组预处理，以便得到高信噪比的协方差矩阵，进一步降低其在非目标方向上输出值。

以上面所示基本数据模型为基础，首先对线列阵2n-1个阵元接收数据按式(4)进行分组处理

然后，按在搜索角度θ处，按式(5)对各组数据进行相移预处理，可得

对各组数据预处理结果进行相加，得到一组新数据为

最后，求取y(fl)协方差矩阵ry(fl)＝e[y(fl)^hy(fl)]，并进行特征分解，可得

式中，e[·]为期望函数，λ＝diag(λ1,λ2,…,λn)为新数据协方差矩阵y(fl)的特征值，u＝[u1,u2,…,un]为新数据协方差矩阵y(fl)的特征向量，uaps为新数据的信号特征向量，uapv为新数据的噪声特征向量。

根据噪声特征向量和目标波达方向导向权向量的正交性，可获得基于分子阵预处理的music波束形成输出的来波方向波束为

式中，为导向权向量。

根据上面所述数据处理过程，可将本发明方法实现过程分为如下步骤：

步骤1)按式(4)所示，首先对线列阵2n-1个阵元接收数据进行分组处理，得到n组数据；

步骤2)按式(5)所示，在搜索角度θ处，对各组数据进行相移预处理，可得n组经过相移处理后的数据；

步骤3)按式(6)所示，对n组数据预处理结果进行相加，得到一组新数据y(fl)；

步骤4)求取y(fl)协方差矩阵ry(fl)＝e[y(fl)y(fl)^h]，并进行特征分解，然后根据噪声特征向量和目标波达方向导向权向量的正交性可得到该搜索角度对应波束值papmusic(fl,θ)；

步骤5)按下式求取本发明方法的宽带空间谱：

本发明的方法(简称apmusic)与现有技术中的方法(简称cbf、music)相比具有明显的优点。

下面结合实例，对本发明方法与现有技术中相关方法的效果进行比较。

为了验证本发明方法可以很好地降低背景噪声和旁瓣级在music波束形成输出空间谱的占有量。下面给出如下数值仿真结果，数值仿真中采用分别采用31、63元均匀线列阵作为接收阵，相邻阵元间距为2m；目标辐射信号频率为375hz，目标相对线列阵波达方向为90°，信号与背景噪声谱级比为0db。

从图2和图3显示结果可知，在非目标方向上，相比music波束形成，本发明方法输出波束中的背景噪声和旁瓣级得到有效降低，数值仿真结果与理论分析相一致。

同时，为了进一步验证本发明方法可以降低背景噪声和旁瓣级对弱目标检测的影响。下面给出如下数值仿真，数值仿真中采用分别采用31、63元均匀线列阵作为接收阵，相邻阵元间距为2m；强、弱目标辐射信号频率均为375hz，强、弱目标相对线列阵波达方向分别为90°和60°，强、弱目标辐射信号谱级比为30db，弱目标与背景噪声谱级比为0db。

从图4和图5显示结果可知，由于music波束形成输出空间谱中的背景噪声和旁瓣级较高，在该仿真条件下，60°方位处的弱目标已经不能很好地在music波束形成输出的空间谱中显示出来，而本发明方法所得空间谱可以很好的显示出60°方位处的弱目标，降低了背景噪声和旁瓣级对弱目标检测的影响。

图6为63元线列阵对强、弱目标辐射信号谱级比为80db时的波束形成结果。对比图5和图6可知，相比music波束形成，本发明方法对60°方位处的弱目标检测能力提高了50db以上，提高music波束形成在实际应用的中普适性。

最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。