肖 军,苏 航,曹宦植
(中国人民解放军 91388部队,广东 湛江 524022)
现代鱼雷自导信号多采用双脉冲形式,且脉冲间隔不固定,使用存储转发方式易导致声诱饵应答不完全,很容易就被识别为假目标,失去声诱饵的干扰作用。边发边收式声诱饵可以在发射应答信号的同时不关闭接收机,保证应答完全的同时大大减小了转发信号的延时,提高了声诱饵模拟的逼真度。
边收边发过程中发射和接收同时工作,转发信号在发射换能器端面方向存在声串漏,容易给接收端造成干扰。所以实现边发边收的关键是如何克服发射换能器端面方向信号及其混响对声诱饵接收机的干扰。用声隔离来定义声诱饵接收机对声诱饵发射信号的抑制能力。在以往的使用经验中,实现有效声隔离的手段主要有3种:1)设计收/发换能器的指向性[1-3];2)在收/发换能器之间放置隔声障板[1-4];3)自适应干扰抵消[3,5]。
本文提出的基于稳健自适应波束形成的声隔离技术是使用空域滤波对干扰进行抑制。首先建立边发边收声诱饵模型,针对 STMV波束形成器在基阵模型失配情况下容易出现信号自消的现象,研究解决波束形成稳健性的新方法即RCB算法,实现有效的干扰抑制和信号提取。仿真和试验数据表明,该方法对端射方向干扰具有良好的抑制效果,且能无失真提取出目标信号,为稳健自适应波束形成在声诱饵边发边收中的应用提供了参考。
图1为边发边收声诱饵的组成示意图。诱饵发射后,悬浮于水中,可近似认为收、发换能器中心在同一垂线上。如图所示,声诱饵接收端收到的干扰来自于该垂线方位附近,信号方位与干扰方位不重叠是实现空域滤波的基础。
图1 边发边收声诱饵的组成示意图Fig.1 Compositional diagram of a receivingwhile-emitting acoustic decoy
使用基阵接收完成声诱饵边发边收的信号处理流程如图2所示。波束形成一般在频域实现,所以首先将阵列输出向量通过傅里叶变换从时域转换到频域。鱼雷寻的信号为窄带信号,该窄带信号的频率在一定带宽范围内随机出现。检测背景噪声为宽带,当信噪比足够时,在进行波束形成前须进行预检测估计,找到信号的频带范围后,进行窄带波束形成,既可以提高信噪比,也可以减小计算量。根据波束形成结果完成信号有无判断,并提取峰值方位波束信号用于回波生成。
图2 基阵接收声诱饵信号处理流程图Fig.2 Signal processing flow chart of array receiving acoustic decoys
在窄带波束形成器中,通过设计合适的加权值,可以有选择地增强某一指定方向的信号,抑制其他方向到达的信号(成为干扰与噪声),提高输出信噪比。
常规波束形成器(CBF)的加权向量与输入信号无关,每个阵元、每个频点采用相同的值,具有最好的稳健性,但空间增益有限,且旁瓣过高,不能有效抑制来自旁瓣区域的干扰。STMV是一种基于导向协方差矩阵的时空统计的最优波束形成方法,其在保证感兴趣方位信号无失真输出的条件下,使输出功率最小,最大限度地提高输出信噪比,该波束形成器具有较强的干扰抑制能力。
通过对不同频点的互谱密度矩阵进行相干积累,计算一定频带的空间谱。STMV可以获得带宽增益并减少收敛时间[7]。导向协方差矩阵可表示为
由于实际中精确的互谱密度矩阵和导向向量很难获得,当快拍数比较少或是存在导向向量误差时,就会使STMV波束形成器性能下降,甚至出现“自消”现象。
为克服导向向量误差引起的方向图性能下降的问题,Li[8]等提出了一种求解稳健性波束形成问题的新方法即RCB算法。
RCB稳健性波束形成的优化问题可以表述为
式中,ε0为设定的导向向量误差范围上界,理想情况下ε0=ε。
从式(3)可以看出,RCB方法的主要目的是从功率估计角度出发,在一定的导向矢量误差范围下,求解最优as,使其对应的信号功率最大,这是RCB算法的物理意义。
其方法的具体设计步骤如下:
2)运用牛顿迭代法解式:
式中γ1和γM分别是最大和最小的特征值。
加权向量为
对角加载方法对导向向量误差与协方差矩阵误差都具有稳健性,RCB方法的信号功率估计能力较好,几乎能正确指示信号功率。ε0过大或过小都会引起性能下降,ε0越大,计算得到的值越小,对角加载量也就越大,RCB波束稳健性越强。同时,为了保证波束形成器对干扰的抑制能力,ε0应该选择尽可能小。
综上所述,RCB算法是基于导向矢量椭圆不确定性的约束条件,从功率估计的角度出发,根据协方差矩阵拟合原理提出的稳健性波束形成算法。该波束形成算法比 STMV有更好的稳健性。
选取阵元数为8的均匀线阵,按照30 kHz半波长布阵,采样频率160 kHz,单拍处理设置为1 024点。设定目标入射方位为0°,30 kHz单频信号,背景噪声为海洋环境噪声,信号功率为149 dB,信噪比-15 dB。仿真时不存在导向向量误差,RCB算法设定导向向量误差上界为0.001。分别进行常规波束形成、STMV波束形成、RCB波束形成,对这 3种方法形成的方位谱的比较,如图3所示。
图3 方位谱Fig.3 Azimuth spectrum
从图3中可以看到,在没有导向向量误差的情况下,3种方法均能正确指示信号功率,STMV和RCB波束宽度更窄,旁瓣干扰抑制效果更好,且RCB旁瓣干扰抑制效果略优于STMV。
湖试试验时接收阵采用8元线列阵,按照30 kHz半波长布阵,阵元间距为0.025 m,采样频率160 kHz,单拍处理设置为1 024点。信号频率为30 kHz,接收端标定为148.6 dB,分别从0°,10°,20°,30°,56°,80°,90°方向入射。干扰频率为20 kHz,接收端标定为155.4 dB,入射角为90°。
考虑RCB对导向向量误差上界ε0的选取,ε0过大或过小都会引起性能下降,其选取准则是尽可能接近真实导向向量误差。因此,首先对数据进行预处理,得到0°情况下,不同参数选取对信号功率的影响,如图4所示。可以看到当参数ε0大于0.3,即可对信号功率进行有效估计,所以0°方向的最佳参数选取值为0.3。运用同样的方法,分别对10°,20°,30°,56°,80°,90°入射角方向进行预处理,得到各个角度应该选取的最佳参数值,如图5所示。
图4 RCB在0°下参数选取与估计功率关系Fig.4 Relationship between parameter selection and power estimation at 0° of RCB
图5 RCB在不同角度下最佳参数选取Fig.5 Selection of optimal parameters from different angles of RCB
通过图 5可以看到,随着入射角度的增大,最佳参数选取值也逐渐增大。因此,参数选取应该随着从0°~90°逐渐增大。从图1可以看到,干扰来自于 90°附近,将[80° 90°]范围内的参数设定小于最佳参数,抑制干扰。将0°信号与90°干扰叠加后,进行常规波束形成、STMV波束形成、RCB波束形成,对这3种方法形成的方位谱进行比较,如图6所示。
图6 方位谱Fig.6 Azimuth spectrum
通过图6可以看到CBF和RCB波束形成都可以对信号功率进行准确估计,估计功率为148.6 dB,但CBF主波束宽,且无法消除干扰信号的影响,方位谱的峰值方位在干扰方位处;STMV无法准确估计信号功率,发生了自消,对信号功率的估计仅为117.5 dB,验证了实际存在导向向量误差的情况下,STMV误将信号视作了干扰进行抑制,产生自消;RCB相比于 CBF,能有效抑制干扰信号的影响,相比于 STMV可以准确指示信号功率,且在根据先验信息对参数进行合适选取后,可获得较窄主瓣和低旁瓣。
分别提取CBF和RCB的0°方位波束谱,如图7所示。可以看到RCB的0°方位波束中20 kHz干扰完全被抑制,产生了较深的零陷;CBF的0°方位波束中干扰仍然存在,相较初始输入干扰被抑制16 dB。
图7 方位谱0°方位波束谱Fig.7 Beam spectrum at 0° of azimuth spectrum
对每拍数据提取的CBF和RCB的0°方位波束谱作傅立叶逆变化,用于回波生成,如图8所示。通过谱分析,可以看到RCB可得到较纯净的目标信号。
图8 信号提取及分析Fig.8 Signal extraction and analysis
本文针对边发边收带来的发射声串漏问题和波束形成中导向向量失配问题,研究了一种稳健性自适应波束形成方法。仿真表明该方法可以正确指示信号功率,试验数据处理结果显示,其性能不受导向向量误差的影响,具有较好的稳健性,且可通过对试验数据预处理制定合理参数选取,抑制干扰对目标信号的影响。本文所研究的方法稳定,计算量小,可为稳健自适应波束形成在声诱饵边发边收中的应用提供参考。