2. 中国人民解放军 91001部队,北京 100000
2. No. 91001 Unit of PLA, Beijing 100000, China
波束形成作为声呐在水下探测中用以提高阵列输出信号增益的基本方法,一直是研究的重点方向。传统波束形成算法归纳起来主要包括经典常规波束形成方法、分裂阵波束形成方法、自适应波束形成方法三类[1]。常规波束形成基于全阵列,优点是稳健性好、计算量小,缺点是其空间分辨力受限于瑞利限,且动态响应受旁瓣影响大[2];分裂阵波束形成是将线列阵分为若干个子阵,通过子阵间互相关等处理,实现测向[3];自适应波束形成是近年来的研究热点,典型的自适应波束形成算法有MVDR算法[4-5]、MUSIC算法[6-7]、ESPRIT算法[8]等,但这些算法对接收信号信噪比和阵列形态要求较高。综上,现有的波束形成算法的稳健性、高分辨力往往不可同时兼顾。
本文针对传统波束形成算法中旁瓣抑制问题,对均匀线阵的波束指向性特点进行分析,将常规波束形成作为一个线性系统,以常规波束形成输出结果作为反卷积算法的输入,利用Richardson-Lucy迭代[9]求解点声源散射函数,对常规波束形成输出结果与声源点散射函数进行卷积运算,实现水下目标方位估计。仿真对比常规波束形成和反卷积波束形成目标角度分辨力和目标噪声包络谱输出信噪比,并对湖试数据进行处理对比,验证了反卷积波束形成算法的弱目标检测和多目标分辨能力。
1 理论基础 1.1 均匀线阵指向性均匀线阵自然指向性可以用公式表示为[10]:
$ R(\theta ) = \left| {\frac{{\sin \left( {\displaystyle\frac{{N{\text π} d}}{\lambda }\sin \theta } \right)}}{{N\sin \left( {\displaystyle\frac{{{\text π} d}}{\lambda }\sin \theta } \right)}}} \right| 。$ | (1) |
其中:
$ \sin \theta = \pm m\frac{\lambda }{d} ,$ | (2) |
自然指向性函数的零点为:
$ \sin (\displaystyle\frac{{N{\text π} d\sin \theta }}{\lambda }) = 0 \text{,} \sin \theta = \pm m\frac{\lambda }{{Nd}},m = 1,2,... ,$ | (3) |
则第一个零点的位置为
$ \Delta \sin \theta = 2\frac{\lambda }{{Nd}}{\kern 1pt} 。$ | (4) |
当信号的入射角度不为
$ R(\theta ) = \left| {\frac{{\sin \left( {\displaystyle\frac{{N(\varphi - \beta )}}{2}} \right)}}{{N\sin \left( {\displaystyle\frac{{(\varphi - \beta )}}{2}} \right)}}} \right| ,$ | (5) |
指向性函数的零点变为:
$ \frac{N}{2}(\varphi - \beta ) = \pm m{\text π} ,m = 1,2,... ,$ | (6) |
即
$\frac{N}{2}\left( {\frac{{2{\rm{ }}{\text π} d}}{\lambda }\sin \theta - \beta } \right) = \pm m{\text π} ,m = 1,2,... ,$ | (7) |
此时的
$ \sin {\theta _0} = \beta \frac{\lambda }{{2{\text π} d}} \pm m\frac{\lambda }{{Nd}} ,m = 1,2, \cdots 。$ | (8) |
从式(3)和式(8)的对比可知,当相邻阵元间存在相移或时延时,线阵指向性函数的零点与主极大位置在
在常规波束形成为一个线性系统的设定条件下,其输出结果可以表示为声源的点散射函数与阵列指向性的线性卷积运算,如下式:
$ r(\sin \theta ) = s(\sin \theta ) * R(\sin \theta )。$ | (9) |
式中:
采用Richardson-Lucy公式对声源点散射函数进行求解可得:
$ \begin{aligned} & {s^{(n + 1)}}(\sin \theta ) = {s^n}(\sin \theta )\int_{ - \infty }^\infty \times\\ & {\frac{{R(\sin \alpha - \sin \theta )}}{{\displaystyle\int_{ - \infty }^\infty {R(\sin \alpha - \sin \theta ) \cdot } {s^n}(\sin \theta ){\rm{d}}\sin \theta }}} r(\sin \alpha ){\rm{d}}\sin \alpha 。\end{aligned} $ | (10) |
式中:
均匀直线阵中第
$ {x_n}(t) = \sum\limits_{i = 1}^K {{s_i}(t)} {e^{ - j\omega {\tau _{ni}}}} + {\rm nois}{e_n}(t), i = 1,2,...,K,$ | (11) |
式中,
$ {\tau _{ni}} = \frac{{d(n - 1)\sin {\theta _i}}}{\lambda } 。$ | (12) |
式中,
由式(12)可以看出,均匀线阵各阵元接收到的信号时延随信号入射角度的变化而变化。
当空间中存在多个目标时,利用反卷积波束形成的高分辨特性,估计各个目标的波达角度,然后根据声波入射角度与时延的对应关系,对各个阵元接收的舰船辐射噪声进行时延补偿,补偿后的信号滤除了背景噪声和其他干扰信号的影响,各阵元时延补偿后的信号相位对其、累加,即可得到纯净的目标辐射信号。最后计算该信号的调制谱,实现多目标的调制谱分离,提高目标检测性能,具体算法流程如图1所示。
![]() |
图 1 调制谱反卷积波束形成分离算法过程 Fig. 1 Algorithm flow of modulation spectrum separation |
仿真条件:采用48阵元的均匀直线水听器阵,相邻阵元间距0.5 m,入射信号为舰船辐射噪声,假定信噪比为10 dB,反卷积瑞利迭代次数为200次。图2给出了2种方法角度分辨力的比较图,其中横轴为角度,纵轴为归一化波束形成输出幅度。
![]() |
图 2 角度分辨力比较 Fig. 2 Angle resolution comparison |
从图2可知,反卷积波束形成的主瓣宽度比常规波束形成更窄,当阵元个数相同时,反卷积波束形成相比常规波束形成拥有更优的角度分辨能力。
图3和图4仿真分析了反卷积波束形成和常规波束形成算法的多目标角度分辨力。0~76 s时,2个目标相距很近,无论是反卷积还是常规算法都无法将二者分开,随着目标的运动,角度差别越来越大,反卷积方法的分辨能力增强,从图3清晰可看到2个目标的运动角度,而图4中,2个波束对应的主瓣较宽且叠加在一起,无法准确分辨2个目标的到达角度。反卷积波束形成角度分辨力虽然具有一定局限性,但是整体计算效果优于常规波束形成。
![]() |
图 3 反卷积波束形成双目标角度分辨力 Fig. 3 Angle resolution of deconvolution beamforming |
![]() |
图 4 常规波束形成双目标角度分辨力 Fig. 4 Angle resolution of conventional beamforming |
考虑2个等强度入射的舰船辐射噪声信号,信号1的波达角度为2.3°,基频为17 Hz,信号2的波达角度为4.6°,基频为47 Hz,阵元个数为48,间距为0.5 m,分析的信号带宽为1~2 kHz,信噪比为0 dB,进行200次反卷积波束形成瑞利迭代。
利用补偿后的信号提取得到的信号1和信号2的DEMON谱如图5和图6所示。从图中可以清晰地看到,阵元接收到的混合信号的包络谱被准确地分离开来,而且由于阵列处理的优势,分离后的包络谱的基频幅度增加,包络谱的输出信噪比较混合信号获得了约14 dB的增益。
![]() |
图 5 信号1包络谱分离 Fig. 5 Envelope spectrum separation of signal 1 |
![]() |
图 6 信号2包络谱分离 Fig. 6 Envelope spectrum separation of signal 2 |
利用反卷积波束形成和常规波束形成处理湖试数据,对比结果如图7和图8所示。从图7可以看到,3条明显的亮线,由左至右分别代表船辐射噪声的海底反射波、直达波和海面反射波。图8中虽然可以看到3条路径的时间历程,但明显信噪比不高,主瓣较宽。对比图7和图8可以看出,反卷积波束形成得到的船辐射噪声多途反射波具有更窄的主瓣宽度,多途角度更容易分辨。
![]() |
图 7 反卷积波束形成湖试结果 Fig. 7 Deconvolution beamforming data results |
![]() |
图 8 常规波束形成湖试结果 Fig. 8 Conventional beamforming data results |
针对波束形成中旁瓣抑制问题,提出一种反卷积波束形成高分辨算法。将常规波束形成作为一个线性系统,由均匀线阵的波束指向性出发,利用Richardson-Lucy对点声源散射函数进行迭代求解,以常规波束形成输出结果作为输入,与声源点散射函数相卷积,抑制阵列指向性函数的影响,获取更窄的主瓣和更低的旁瓣,精确估计多目标波达角度,延时补偿各阵元收到的舰船辐射噪声,补偿后的信号相位对齐并累加,即可得到纯净的目标辐射噪声信号。仿真对比了常规波束形成和反卷积波束形成单目标角度分辨力,湖试数据处理对比了常规波束形成和反卷积波束形成多目标角度分辨力,验证了算法的高分辨特性;仿真对比了2个等强度目标不同波达角下,利用反卷积波束形成与常规波束形成获取的目标噪声包络谱输出信噪比,验证了反卷积波束形成算法的性能。结果表明,反卷积波束形成算法具有更优的多目标角度分辨力,能够实现多目标的调制谱分离,有效提高目标检测性能。
[1] |
杜选民, 周胜增, 高源. 声纳阵列信号处理技术[M]. 北京: 电子工业出版社, 2018.
|
[2] |
季浩然, 马晓川, 张舒皓, 等. 功率受限的常规波束形成后处理拟合方法[J]. 声学学报, 2020, 45(1): 1-14. |
[3] |
李兰瑞, 章新华, 李鹏, 等. 基于时延估计的分裂阵时域波束形成技术[J]. 舰船科学技术, 2019, 41(11): 155-156. DOI:10.3404/j.issn.1672-7649.2019.11.032 |
[4] |
张思宇, 何心怡, 程善政. 基于少快拍条件均匀圆阵波束域MVDR方法[J]. 舰船科学技术, 2019, 41(11): 148-151. DOI:10.3404/j.issn.1672-7649.2019.11.030 |
[5] |
董晋, 胡鹏, 冯金鹿. 双约束条件的时域MVDR自适应波束形成方法[J]. 舰船科学技术, 2021, 43(7): 127-130. |
[6] |
姚昕彤, 王玉文, 刘奇. 基于MUSIC及其改进算法的DOA估计研究[J]. 通信技术, 2021, 54(6): 1363-1369. DOI:10.3969/j.issn.1002-0802.2021.06.012 |
[7] |
智婉君, 严胜刚, 李志舜. 宽带方位估计的波束域Root-MUSIC算法[J]. 电子与信息学报, 2002, 24(5): 644-648. |
[8] |
刘松, 庞育才. 基于扩展ESPRIT的随机阵列高效DoA估计算法[J]. 电子与信息学报, 2019, 41(6): 6: 8.
|
[9] |
YANG T C. Performance analysis of superdirectivity of circular arrays and implications for sonar systems[J]. IEEE Journal of Oceanic Engineering, 2019, 44(1): 156-166. DOI:10.1109/JOE.2018.2801144 |
[10] |
Michael A. Ainslie. 声呐性能建模原理[M]. 北京: 国防工业出版社, 2015.
|