(19)中华人民共和国国家知识产权局
(12)发明专利申请
(10)申请公布号 CN 110018465 A(43)申请公布日 2019.07.16
(21)申请号 201810018101.X(22)申请日 2018.01.09
(71)申请人 中国科学院声学研究所
地址 100190 北京市海淀区北四环西路21
号(72)发明人 郑恩明 陈新华 李媛 宋春楠 (74)专利代理机构 北京方安思达知识产权代理
有限公司 11472
代理人 陈琳琳 张红生(51)Int.Cl.
G01S 7/537(2006.01)G01S 15/88(2006.01)
权利要求书2页 说明书6页 附图3页
(54)发明名称
一种基于全相位预处理的MVDR波束形成方法
(57)摘要
本发明公开了一种基于全相位预处理的MVDR波束形成方法，包括：根据线列阵接收数据中信号和噪声相关性的差异，将线列阵中2N-1个阵元接收数据通过全相位预处理转变为N个阵元数据；然后对N个阵元数据进行处理，得到基于全相位预处理的MVDR波束输出。本发明方法对线列阵接收数据进行了全相位预处理，有效提高线列阵接收数据协方差矩阵中信号含有量和信噪比，降低了背景噪声和旁瓣级对MVDR波束形成检测弱目标带来的影响，提高了MVDR波束形成对弱目标检测的检测效果。CN 110018465 ACN 110018465 A
权 利 要 求 书
1/2页
1.一种基于全相位预处理的MVDR波束形成方法，包括：根据线列阵接收数据中信号和噪声相关性的差异，将线列阵中2N-1个阵元接收数据通过全相位预处理转变为N个阵元数据；然后对N个阵元数据进行处理，得到基于全相位预处理的MVDR波束输出。
2.根据权利要求1所述的基于全相位预处理的MVDR波束形成方法，其特征在于，所述方法具体包括：
步骤1)对线列阵2N-1个阵元接收数据按下式进行分组处理：
式中，Xn(fl)表示第n个阵元拾取的频率fl数据，可表示为：
式中：S(fl)为目标辐射信号，Nn(fl)为第n个阵元拾取的加性高斯白噪声数据，λ＝fl/c为波长，d为线列阵相邻阵元间距，θc为声速入射；0为目标相对线列阵入射角度，
步骤2)在搜索角度θ,θ＝1,2,…180处，按下式对各组数据进行相移预处理：
步骤3)对各组数据预处理结果进行相加，得到一组新数据为：
2
CN 110018465 A
权 利 要 求 书
2/2页
步骤4)求取Y(fl)协方差矩阵RY(fl)＝E[Y(fl)HY(fl)]，获得基于全相位预处理的MVDR波束形成输出的来波方向波束为：
式中，E[·]为期望函数，
为导向权向量，τn＝(n-1)dcos
(θ)/c,1≤n≤N。
3.根据权利要求2所述的基于全相位预处理的MVDR波束形成方法，其特征在于，所述方法还包括：
步骤5)按下式求取波束输出的宽带空间谱：
其中，L为频带数。
3
CN 110018465 A
说 明 书
一种基于全相位预处理的MVDR波束形成方法
1/6页
技术领域
[0001]本发明涉及声纳信号处理领域，特别涉及一种基于全相位预处理的MVDR波束形成方法。
背景技术
[0002]水下目标检测与估计是阵列信号处理的一个重要分支。波束形成作为阵列信号处理中核心算法，输出波束中的背景噪声和旁瓣级一直是其设计中需要考虑的问题。低背景噪声和旁瓣级可以有效降低对位于强目标旁瓣区域内弱目标检测的漏报概率。[0003]为了对波束形成输出背景噪声和旁瓣级实现控制，许多学者从不同方法对降低波束形成旁瓣级进行了深入研究，并取得了一定的研究成果，提出了很多方法，主要为Chebyshev滤波方法，“凹槽噪声场”方法，静态波束图数字综合方法，反复迭代方法，多线性约束方法，非线性优化方法，凸优化(Convex Optimization)方法，半无限二次规划(Semi-Infinite Quadratic Programming)方法，二阶锥(Second-Order Cone)约束方法，中心矩方法[，虚拟干扰源构造能量聚焦矩阵方法，稀疏约束方法。在以上方法中，Chebyshev滤波方法以其简单方便常被应用实际工程中，但存在旁瓣级设置和主瓣宽度控制折中选择问题。
[0004]目前的方法都未能解决针对MVDR波束形成中背景噪声和旁瓣级对检测弱目标性能影响的问题。
发明内容
[0005]本发明的目的在于针对MVDR波束形成中背景噪声和旁瓣级对检测弱目标性能影响的问题，根据线列阵接收数据中信号和噪声相关性的差异，提出一种基于全相位预处理的MVDR波束形成方法。该方法首先需要将线列阵中2N-1个阵元接收数据通过全相位预处理转变为N个阵元数据；然后再采用MVDR波束形成思想对该N个阵元数据进行处理，可得到各方位处空间谱。由于该方法对线列阵接收数据进行了全相位预处理，有效提高线列阵接收数据协方差矩阵中信号含有量和信噪比，降低了背景噪声和旁瓣级对MVDR波束形成检测弱目标带来的影响，提高了MVDR波束形成对弱目标检测的检测效果。[0006]为了实现上述目的，本发明一种基于全相位预处理的MVDR波束形成方法，所述方法包括：
[0007]根据线列阵接收数据中信号和噪声相关性的差异，将线列阵中2N-1个阵元接收数据通过全相位预处理转变为N个阵元数据；然后对N个阵元数据进行处理，得到基于全相位预处理的MVDR波束输出。
[0008]作为上述方法的一种改进，所述方法具体包括：
[0009]步骤1)对线列阵2N-1个阵元接收数据按下式进行分组处理：
4
CN 110018465 A
说 明 书
2/6页
[0010]
[0011][0012]
式中，Xn(fl)表示第n个阵元拾取的频率fl数据，可表示为：
式中：S(fl)为目标辐射信号，Nn(fl)为第n个阵元拾取的加性高斯白噪声数据，λ＝
fl/c为波长，d为线列阵相邻阵元间距，θc为声速入射；0为目标相对线列阵入射角度，[0014]步骤2)在搜索角度θ,θ＝1,2,…180处，按下式对各组数据进行相移预处理：
[0013]
[0015]
[0016]
步骤3)对各组数据预处理结果进行相加，得到一组新数据为：
[0017]
[0018]
步骤4)求取Y(fl)协方差矩阵RY(fl)＝E[Y(fl)HY(fl)]，获得基于全相位预处理的MVDR波束形成输出的来波方向波束为：
5
CN 110018465 A[0019][0020]
说 明 书
3/6页
式中，E[·]为期望函数，
为导向权向量，τn＝(n-1)
dcos(θ)/c,1≤n≤N。
[0021]作为上述方法的一种改进，所述方法还包括：[0022]步骤5)按下式求取波束输出的宽带空间谱：
[0023][0024]
其中，L为频带数。
[0025]本发明的优点在于：
[0026]本发明方法对线列阵接收数据进行了全相位预处理，有效提高线列阵接收数据协方差矩阵中信号含有量和信噪比，降低了背景噪声和旁瓣级对MVDR波束形成检测弱目标带来的影响，提高了MVDR波束形成对弱目标检测的检测效果。附图说明
[0027]图1是本发明所述拖线阵声纳结构示意图；
[0028]图2为本发明的方法的协方差矩阵信号含有增加量的示意图；[0029]图3为本发明的方法的协方差矩阵对角线能量增加量的示意图；
[0030]图4为本发明的方法与现有方法对31元线列阵波束形成的比较结果图；[0031]图5为本发明的方法与现有方法对63元线列阵波束形成的比较结果图；[0032]图6为本发明的方法与现有方法对63元线列阵波束形成的比较结果图(强、弱目标辐射信号谱级比为30dB)；
[0033]图7为本发明的方法与现有方法对63元线列阵波束形成的比较结果图(强、弱目标辐射信号谱级比为40dB)。
具体实施方式
[0034]现结合附图对本发明作进一步的描述。[0035]在对本发明的方法做详细说明前，首先对本发明的方法所适用接收阵加以描述。图1为一拖线阵声纳结构示意图，该拖线阵声纳包括6个部分，显控与信号处理机1、甲板缆2、绞车3、导缆架4、拖缆5、接收线阵6。其中的接收线阵6通过拖缆5与位于绞车3上的甲板缆2连接，所述拖缆5还安装在导缆架4上；接收线阵6所接收的信号传输给显控与信号处理机1。
[0036]下面对本发明方法做进一步的说明。[0037]MVDR波束形成数学表达式
[0038]对于间距为d的2N-1元等间隔水平线阵，有1个目标从θ则第n个阵元拾取的0入射，频率fl数据Xn(fl)可表示为：
[0039][0040]
式中：S(fl)为目标辐射信号，Nn(fl)为第n个阵元拾取的加性高斯白噪声数据，c为
6
CN 110018465 A
说 明 书
4/6页
声速，λ＝fl/c为波长。
[0041]对线列阵各阵元接收数据构造数据矩阵，可表示为[0042]X(fl)＝[X1(fl),X2(fl),…,X2N-1(fl)]T (2)[0043]然后，求取线列阵协方差矩阵RX(fl)＝E[X(fl)X(fl)H]，可获得来波方向输出波束为
[0044][0045]
式中，为导向权向量，τθ)/c，θ为搜索角n＝(n-1)dcos(
度，c为声速。
[0046]基于全相位预处理的MVDR波束形成方法[0047]数学模型
[0048]为了进一步降低在非目标波达方向上，MVDR波束形成输出波束在不同搜索角度θ形成的极大值，减小其对弱目标检测形成影响。本发明将根据线列阵接收数据协方差矩阵形成过程中信号和噪声相关性的差异，通过对线列阵接收数据进行分组预处理，以便得到高信噪比的协方差矩阵，进一步降低其在非目标方向上输出值。[0049]以上面所示基本数据模型为基础，首先对线列阵2N-1个阵元接收数据按式(4)进行分组处理
[0050]
[0051]
然后，按在搜索角度θ处，按式(5)对各组数据进行相移预处理，可得
[0052]
[0053]
对各组数据预处理结果进行相加，得到一组新数据为
7
CN 110018465 A
说 明 书
5/6页
[0054]
[0055]
最后，求取Y(fl)协方差矩阵RY(fl)＝E[Y(fl)HY(fl)]，可获得基于全相位预处理的
MVDR波束形成输出的来波方向波束为
[0056][0057][0058][0059]
式中，E[·]为期望函数，
为导向权向量。
根据上面所述数据处理过程，可将本发明方法实现过程分为如下步骤：步骤1)按式(4)所示，首先对线列阵2N-1个阵元接收数据进行分组处理，得到N组
数据；
步骤2)按式(5)所示，在搜索角度θ处，对各组数据进行相移预处理，可得N组经过
相移处理后的数据；
[0061]步骤3)按式(6)所示，对N组数据预处理结果进行相加，得到一组新数据Y(fl)；[0062]步骤4)求取Y(fl)协方差矩阵RY(fl)＝E[Y(fl)Y(fl)H]，并进行矩阵求逆，然后根式(7)可得到该搜索角度对应波束值PAPMVDR(fl,θ)；
[0063]步骤5)按下式求取本发明方法的宽带空间谱：
[00][0065][0060]
其中，L为频带数。
[0066]为了进一步验证本文方法中协方差矩阵信号增加量，进行如下数值仿真，仿真中采用8:8:128元均匀线列阵作为接收阵，接收数据信号和背景噪声谱级比为0dB，数值仿真结果如图2和图3所示，每一种阵元所得结果均是由100次统计所得。[0067]图2和图3数值仿真结果进一步验证了，经全相位预处理后的新数据协方差矩阵信号含有量增加了10lg(N2/(2N-1))dB，协方差矩阵主对角线元素能量由原来的变为线阵的
[0068]
的正确性。。
本发明的方法(简称APMVDR)与现有技术中的方法(简称CBF、MVDR)相比具有明显的优点。
[0069]下面结合实例，对本发明方法与现有技术中相关方法的效果进行比较。
8
CN 110018465 A[0070]
说 明 书
6/6页
为了验证本发明方法可以很好地降低背景噪声和旁瓣级在MVDR波束形成输出空
间谱的占有量。下面给出如下数值仿真结果，数值仿真中采用分别采用31、63元均匀线列阵作为接收阵，相邻阵元间距为2m；目标辐射信号频率为375Hz，目标相对线列阵波达方向为90°，信号与背景噪声谱级比为0dB。[0071]从图4和图5显示结果可知，在非目标方向上，相比MVDR波束形成，本发明方法输出波束中的背景噪声和旁瓣级得到有效降低，数值仿真结果与理论分析相一致。[0072]同时，为了进一步验证本发明方法可以降低背景噪声和旁瓣级对弱目标检测的影响。下面给出如下数值仿真，数值仿真中采用分别采用63元均匀线列阵作为接收阵，相邻阵元间距为2m；强、弱目标辐射信号频率均为375Hz，强、弱目标相对线列阵波达方向分别为90°和60°，强、弱目标辐射信号谱级比为30dB，弱目标与背景噪声谱级比为0dB。[0073]从图6显示结果可知，由于MVDR波束形成输出空间谱中的背景噪声和旁瓣级较高，在该仿真条件下，60°方位处的弱目标已经不能很好地在MVDR波束形成输出的空间谱中显示出来，而本发明方法所得空间谱可以很好的显示出60°方位处的弱目标，降低了背景噪声和旁瓣级对弱目标检测的影响。[0074]图7为63元线列阵对强、弱目标辐射信号谱级比为40dB时的波束形成结果。对比图6和图7可知，相比MVDR波束形成，本发明方法对60°方位处的弱目标检测能力提高了10dB以上，提高MVDR波束形成在实际应用的中普适性。[0075]最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。
9
CN 110018465 A
说 明 书 附 图
1/3页
图1
图2
图3
10
CN 110018465 A
说 明 书 附 图
2/3页
图4
图5
11
CN 110018465 A
说 明 书 附 图
3/3页
图6
图7
12
