在水下作战时,随着安静型/低噪声潜艇、声隐身水面舰等出现,降低了目标的辐射/反射噪声强度,使声呐系统探测能力快速下降。弱目标探测的难点在于干扰功率和背景噪声强度远大于弱目标噪声,干扰数的数量远大于弱目标数量。因此在被动声呐时间-方位历程(Bearing Time Recording,BTR)显示中,弱目标显示往往受到干扰目标和海洋环境背景噪声的影响,从而导致目标(特别是弱目标)的无法被检测和跟踪[1-2]。
为了提高弱目标的检测能力和显示效果,国内外已经对其做了大量的研究。M.Bono[3-4]提出了子频段峰值能量检测(Sub-band Peak Energy Detection,SPED)算法,该算法根据频率方位图信息,利用信号峰值的“空间一致性”,对每个子频段各个方位上的峰值能量进行求和作为波束输出,以改善声呐显示效果。文献[5-6]提出了一种波束域宽带峰值能量检测方法(Beam-domain Broadband Peak Energy Detection),该方法增强目标出现的波束显示强度,但是对于目标不存在的波束降低了显示强度,这有可能降低非目标所在波束的干扰,达到改善目标显示效果的目的。文献[7]提出一种利用Eckart后置处理滤波器来提高被动声呐宽带信号检测的方法。
本文综合利用文献[3, 7]的思想,利用子频段峰值检测和Eckart后置处理滤波器相结合,提出一种基于子频段的弱目标检测方法。仿真结果和海试数据分析结果表明,该方法能够有效提高弱目标的检测效果。
1 子频段峰值能量检测假设声呐阵元接收到的时间快拍为
$P(f,{\theta _i}) = w(f,{\theta _i}) \cdot X(f),$ | (1) |
其中
假设整个处理频带为
$P({f_l},{\theta _i}) = \sum\limits_{f = {f_{l,\min }}}^{{f_{l,\max }}} {P(f,{\theta _i})}, $ | (2) |
对每个
$\begin{split}& {P}_{pks}({f}_{l},{\theta }_{i})=\\ &\left\{ \begin{array}{*{20}{c}}P{\rm{'}}({f}_{l},{\theta }_{i}) &P{\rm{'}}({f}_{l},{\theta }_{i})>P{\rm{'}}({f}_{l},{\theta }_{i-1})\&\&P{\rm{'}}({f}_{l},{\theta }_{i})>P{\rm{'}}({f}_{l},{\theta }_{i+1}),\\ 0& \text{otherwise}\text{。}\end{array} \right.\end{split} $ | (3) |
最后对同一波束上的
${P_{subband}}({\theta _i}) = \sum\limits_{l = 1}^L {{P_{pks}}({f_l},{\theta _i})}\text{。}$ | (4) |
从式(2)和式(3)可以看出,如果在某一方位上存在宽带信号,则在该方位上对应的几个子频段中都会出现局部峰值,累加后信号得到加强。而噪声在各个子频段内随机出现峰值,累加后其能量被平滑。
2 基于子频段的弱目标检测方法基于子频段的弱目标检测方法把子频段峰值能量检测和Eckart加权方法相结合。Eckart[8]滤波器基于偏差最大准则,它使有信号存在时均值相关器输出的变化与无信号存在时相关器输出的标准偏差之比最大化。这种滤波器的特点是可以抑制噪声较强频段的响应,改善时延估计的精度。其加权函数的定义如下
$W(f) = {G_{ss}}(f)/\left[ {{G_{{n_1}{n_1}}}(f){G_{{n_2}{n_2}}}(f)} \right]\text{。}$ | (5) |
其中:
从定义中可以看出,Eckart滤波器需要了解信号和噪声的功率谱。在被动声呐目标检测中,利用序截断平均(OTA)[9]算法来估计背景噪声。该算法是Struzinski在讨论频率域历程检测时提出来的。
每个子频段进行波束形成后得到子频段谱矩阵为
$\begin{split}&P'({f_l},{\theta _{K + 1}}),\cdots,P'({f_l},{\theta _2}),P'({f_l},{\theta _1}),\cdots,\\ &P'({f_l},{\theta _I}),P'({f_l},{\theta _{I - 1}}),\cdots,P'({f_l},{\theta _{I - K}}),\end{split}$ | (6) |
处理后,左右端点在扩展后的数据中都处于窗口中,避免了边缘效应。将式(6)重新表示为:
$Y(1),Y(2),\cdots,Y(M + 2K),$ | (7) |
OTA算法只针对
$y(1),y(2),\cdots,y(2K + 1),$ | (8) |
式(5)序列中的中位数为
$\bar y = \frac{1}{{K + 1}}\sum\limits_{m = 1}^{K + 1} {y(K + m)}, $ | (9) |
根据Eckart滤波器定义,得到子频段加权系数矩阵,即
${{\rm{W}}}({f_l},{\theta _i}) = \left[ {\begin{array}{*{20}{c}} {P'({f_1},{\theta _1})/{N^2}({f_1},{\theta _1})}& \cdots &{P'({f_1},{\theta _I})/{N^2}({f_1},{\theta _I})} \\ \cdots & \ddots & \cdots \\ {P'({f_L},{\theta _1})/{N^2}({f_L},{\theta _1})}& \cdots &{P'({f_L},{\theta _I})/{N^2}({f_L},{\theta _I})} \end{array}} \right],$ | (10) |
对子频段进行加权后,得到子频段最终输出结果
${P_{subband}}({\theta _i}) = \sum\limits_{l = 1}^L {{P_{pks}}({f_l},{\theta _i})W({f_l},{\theta _i})}\text{。} $ | (11) |
为了验证算法的有效性,对其进行多目标仿真。仿真基阵为128元线列阵,阵元间距为0.1 m,其工作频段范围为1000~3500 Hz,假设目标的辐射噪声为带限信号,其速度大小为相对于平台的速度,目标的参数如表1所示。数据处理采样频率为15 kHz,采样点数为2048,OTA算法中K=11,子频段数目为100。
某一时刻的波束形成结果如图1所示。
从图2和图3可以看出,常规算法中目标1轨迹基本不可见,目标3轨迹隐约可见,目标4轨迹断断续续。从图4可以看出,本文所提方法4个目标的轨迹都连续可见。仿真试验表明,本文所提的方法能够提高弱目标的检测能力。
为了进一步验证所提算法的有效性,给出一组海试数据结果。数据为116元线列阵,海试数据的波束形成结果如图5所示,方位历程结果如图6~图7所示。
从图5和图6可以看出,目标1的轨迹几乎不可见,目标难以被发现;从图7中可以看出,目标1的轨迹较为清晰。图8提供了300 s时刻的波束形成归一化幅度图,CBF和CBF-SPED方法能够检测目标1,但其相对幅值较小,本文所提算法检测结果优于CBF和CBF-SPED。
本文结合了子频段峰值能量检测方法与Eckart后置处理滤波方法,提出一种基于子频段的弱目标检测方法。它通过对子频段谱进行Eckart加权,提高了弱目标在环境背景噪声和强干扰的检测能力,仿真试验和海试数据分析验证了该方法的有效性。在计算Eckart权值时,所提方法无法准确估计背景噪声功率谱和信号的功率谱,从而会影响弱目标的检测,在后续的工作中将继续研究改进。
[1] |
ZHANG Wei, CONG Mingyu, WANG Liping. Algorithms for optical weak small targets detection and tracking: review[c]// IEEE int. Conf. Neural Networks & Signal Processing, pp.643−647, 2003.
|
[2] |
BRIAN H. Maranda, JOHN A. Fawcett. Detection and localization of weak targets by space-time integration[J]. IEEE Journal of Oceanic Engineering, Vol. 16, No.2, pp.189−194, 1991.
|
[3] |
BONO M, SHAPO B. Sub-band energy detection in passive array processing[R].ADA405484, 2000.
|
[4] |
M.Bono. SEED ED vs. SPED ED Evaluation. Signal processing Working Group, Applied Research Laboratories: The University of Texas at Austin, ARL-TL-EV-00-61.
|
[5] |
杨晨辉, 马远良, 杨益新. 峰值能量检测及其在被动声纳显示中的应用[J]. 应用声学, 2003, 22(5): 31-35. DOI:10.3969/j.issn.1000-310X.2003.05.008 |
[6] |
杨晨辉. 被动声纳宽带信号检测与现实方法研究[D]. 西安: 西北工业大学, 2002.
|
[7] |
S.k.MEHTA, J.FAY, P.MACIEJEWSKI. A modified Eckart post-beamformer filter for improved detection using broadband features[c]//Proc.ICASSP 96, 3045−3048.
|
[8] |
邱天爽, 王宏禹. Eckart加权的自适应时间延迟估计[J]. 大连理工学报, 1996, 36(2): 230-235. |
[9] |
李启虎, 潘学宝, 尹力. 数字式声呐中一种新的背景均衡算法[J]. 声学学报, 2000, 25(1): 5-9. DOI:10.3321/j.issn:0371-0025.2000.01.002 |