基于时延估计的分裂阵时域波束形成技术

李兰瑞，章新华,，李 鹏，刘心语

(1.海军大连舰艇学院 水武防化系，辽宁 大连 116018；2.哈尔滨工程大学 水声工程学院，黑龙江 哈尔滨 1500011；3.海军大连舰艇学院，辽宁 大连 116018)

0 引 言

拖曳线列阵相对于安装在平台外壳上的舰壳式声呐而言，具有不受载体平台安装尺寸限制、受本平台噪声影响相对较小、具有较大的阵列孔径和较低的工作频段等优势。对于大孔径拖曳线阵，在实际使用中受舰艇横向机动、洋流影响和水动力影响会产生一定的形变[1]。由于阵列形变，对于某θ方向上的来波信号，会导致阵元间的理论声程差和实际声程差不一致。这时，依据理论阵列流形计算出的时延差和真实时延差也会出现差距，不能实现信号的同相叠加,目标信号无法获得最大增益，减弱了波束分辨能力。针对这一问题目前常采用阵形估计算法估计各阵元的实时位置，柔性长拖曳线阵阵形估计主要有两大类：一类是安装传感器的借助硬件估计阵形方法；一类是采用水听器信号进行阵形估计的方法[2]。阵形估计算法在一定程度上还原了阵元的实时位置，对波束形成算法中计算导向矢量具有重要意义。但是阵形估计算法实时性不强、计算量大，并且阵形估计不准确时会引入新的误差。

常规时域波束形成算法存在增益低和波束分辨能力差的问题，Macdonald和和Schultheiss的研究表明分裂阵波束形成方法具有接近最优的目标估计精度，能够有效提高对目标的分辨能力[3-4]。分裂阵波束形成方法在线列阵精确测向中已得到了广泛应用，其中基于左右波束相位差测向的互谱法[5]和对左右波束采取和差运算的超波束形成算法[6]是2种高分辨率的分裂阵波束形成方法，但这2种算法均是在频域而非时域进行处理，无法输出听音波束信号。

针对以上问题，在无法进行阵形估计时，本文将大孔径拖曳线阵分为左右双子阵分别做波束形成，通过最大似然时延估计算法估算对应波束的时延差，再依据估算时延差对左右波束进行延时求和得到最终的波束信号。仿真和海试数据证明，相对于全阵直接做波束形成的方法，基于双子阵的时域波束形成技术有效提高了信号增益和波束分辨率。

1 阵形畸变对阵增益的影响

1.1 阵增益

设均匀直线阵总阵元数为2N，阵元无指向性，阵元间距为d，目标入射角相对阵列方向夹角为θ，理想条件下，第i个阵元接收信号记为：为第i路阵元信号相对第1号参考阵元的时延差，ni（t）为各自独立的噪声。理想情况下波束形成之后系统输出为：

当θ=θ0时,τi（θ）=τi（θ0），D理论（θ）=（2N）2+2N，其中为信号均方差、为噪声均方差。可以看出信号增强了（2N）2倍，而噪声增强了2N倍。波束形成所带来的增益Gs理论为[7]：

通常将空间增益取对数转换为功率形式：

1.2 阵形畸变对阵增益的影响

由于受阵形畸变影响，阵元位置发生偏移进而引入新的随机时延差τj（θ），第i号阵元实际接收到的信号为：

此时波束形成输出D实际（θ）为：

τj（θ）≠0 时τi（θ）+τj（θ）≠τi（θ0），波形失配使得此时的增益Gs实际＜Gs理论。阵形畸变带来的时延误差τj（θ）公式为：

式中α（j）为取值大于零的随机畸变系数。

对于全阵，在波束形成时阵形畸变带来的误差值累积量μ为:

s.e. 0.0167 0.0080 0.0215 0.0347 0.0457 0.0562 0.0568 0.0442

由误差累积公式可以看出，在阵元间距一定的条件下，阵形畸变带来的误差累积量同阵元数成正比例关系。

2 基于时延估计的分裂阵时域波束形成

长拖曳线阵畸变会引入畸变误差，降低增益和波束分辨能力,并且在阵元间距一定的条件下，阵形畸变累积值μ将随着阵元数目的增多而非线性增大。另外，阵增益同阵元数成正比例关系，而阵元数又同阵形畸变误差累积值成反比例关系。综上，在无法进行阵形估计时，为保证阵增益的同时降低阵形畸变带来的影响，本文提出一种基于广义互相关估计时延的分裂阵时域波束形成方法。将长阵分为左右2个子阵，左右子阵分别按照各自的参考阵元进行波束形成，之后再将两左右波束合成为一个波束输出。左右子阵分别做波束形成，阵元数降为一半，畸变误差累积值被降低。对左右波束重新估计时延之后再延时累加，又利用了全阵信息，此时的阵增益理论上同全阵直接做波束形成时相同。

2.1 算法原理与技术流程

设总阵元数为2N，阵元间距为d，目标入射角为θ，将均匀线阵等分为左右2个子阵，每个子阵的阵元数为N，其中1号、N/2+1号阵元分别为左右子阵的参考阵元，阵形如图1所示。

图1 分裂线阵示意图Fig.1 Schematic diagram of split line array

基于分裂阵时延估计时域波束形成算法计算流程如下：

1）将均匀线阵等分为左右2个子阵，采用辛格函数法进行时延的常规时域波束形成方法[8-9]对阵列信号分别进行时域波束形成，得到左右波束信号，其中M为波束数。由于左右阵元数目一致，左右波束的理论信信比为0 dB。

2）左右阵元各自做时域波束形成之后，两子阵可以各自声中心等效为2个阵元，对于θ方位入射的信号，等效阵元的理论时延差为：

由于阵形畸变，τ理论（θ）不等于等效阵元时延差的真实值，所以要对左右波束信号进行时延估计。另外由于阵形畸变，左右子阵对同一目标估计出的方位也会出现偏差。所以，估计左右波束时延差时，将分别做广义互相关计算。如果是的时延信号，且当与它们两者相关的时延相等时，相关函数出现最大值：

最大似然加权时延估计流程图如图2所示。

图2 最大似然加权时延估计流程图Fig.2 Maximum likelihood weighted time delay estimation flow chart

通过分子阵做波束形成，再用时延估计方法而非理论公式计算左右波束的时延差，减小了阵形畸变带来的影响。

4）对M路预成波束信号分别累加求能量得到空间谱。对空间谱进行时间累积得到方位历程图，算法流程图如图3所示。

图3 分裂阵算法流程图Fig.3 Flowchart of split matrix algorithm

3 算法验证与分析

相对于全阵直接做时域波束形成方法，基于时延估计的分裂阵时域波束形成方法减小了阵形畸变带来的误差，延时后的各路信号更趋近于同相叠加。为验证本文算法的科学性和有效性，利用仿真和海试数据从对算法进行验证。

3.1 仿真验证与分析

仿真阵列为均匀线阵，阵元数目512个，阵间距2.4 m，采样频率为5 000 Hz，目标源信号A和B均为海上实录的商船信号，其中目标A设置在97°方位，目标B设置在100°方位。线阵正常工作状态下应为直线阵，加入扰动后阵形畸变，扰动前后阵形示意图如图4所示。

图4 阵列扰动前后示意图Fig.4 Schematic diagram before and after array disturbance

对比加入扰动前后，全阵直接做时域波束形成与分子阵之后时延估计做波束形成得到的空间谱，如图5和图6所示。

图5 扰动前空间谱对比图Fig.5 Comparison of spatial spectrum before disturbance

图6 扰动后空间谱对比图Fig.6 Comparison of spatial spectrum after disturbance

图5和图6仿真结果表明，均匀线阵加入扰动前，常规全阵做波束形成同分子阵估计时延再做波束形成均可以分开目标A，B，目标检测信噪比基本一致。加入扰动后，常规全阵波束形成算法不能有效分开目标A，B，而分子阵估计时延的方法依然可以有效检测到目标A，B，且检测信噪比相对常规全阵做法更高。

3.2 海试数据验证

仿真信号验证了分裂阵方法的有效性，采用某海试数据进行验证，海试数据为某型柔性长拖曳声呐的海上实录数据。对比全阵直接做时域波束形成和分子阵波束形成方法输出的空间谱和方位历程图，如图7和图8所示。

对比图7和图8可以发现，图中标注的A，B，C三处，常规全阵时域波束形成算法无法有效分开临近双目标，而分裂阵算法提高了长拖曳线阵分辨率，可以有效检测到临近目标，检测信噪比更高。

图7 空间谱对比图Fig.7 Comparison Chart of spatial spectrum

图8 方位历程对比图Fig.8 Comparison of azimuth process

4 结 语

本文首先分析了柔性长拖曳线阵阵形畸变对波束形成带来的影响，之后结合频域分裂阵处理思想提出了基于时延估计的分裂阵时域波束形成算法。仿真和海试数据证明，分裂阵时域波束形成算法有效降低了阵形畸变带来的影响，提高了波束分辨率和信号增益。由于提高了信号增益，时域波束形成后输出的时域声音信号更有利于声呐兵的听音识别。在无法进行阵形估计时，可有效降低阵形畸变带来的影响。