基于水下高速小目标末弹道测量的实孔径成像声呐波束形成研究

张思宇，何心怡，陈 菁，王 磊

(1.海军研究院，北京 100161；2.中国科学院声学研究所，北京 100097)

0 引 言

在声呐系统的设计过程中，波束形成是一个核心基础环节，利用一定几何形状的基阵和信号处理技术，使发射和接收的信号形成指向性波束。在发射信号时，通过指向性波束提高声呐的作用距离；在接收信号时，在空间域中抑制环境噪声和混响干扰，提高声呐接收机输入端的信噪比，获得更好的声呐图像和目标的分辨率，其直接表征为旁瓣级，旁瓣级越低，系统抑制旁瓣区域的干扰能力越强。

在没有水下弹道测量系统的条件下，设计测量水下高速小目标末弹道的声呐成像系统的主要考虑：1）目标接近成像系统的方向未知，可能从任一方向接近；2）目标运动速度快，回波多普勒偏移量大，可能达到宽带信号的标准；3）目标距离近，成像系统受到混响的干扰强等因素的限制[1]。在声呐系统的阵型设计中，只有使用均匀排布的圆形阵列，才具有 0°～360°的方位角范围内，形成均匀指向性波束的能力，避免目标的方位模糊，圆阵波束形成数学模型简单，设计的工程实现性强。但是圆阵的常规波束形成的旁瓣级高，图像的目标分辨率低，弹道测量误差大，本文以一种对均匀圆阵进行Dolph-Chebyshev 加权的波束形成方法，并给出在圆阵利用此法波束形成的约束条件，并通过仿真实验验证方法的可行性。

1 常规波束形成模型

由于声信号到达基阵各个阵元的距离不同，导致各个阵元输出的信号也不同，需要通过延时等信号处理方法补偿差异，这是常规波束形成的基本思想[2-3]。假设入射信号是平面波，基阵是N 个阵元组成。以其中某个阵元为参考点， 设第 i个阵元接收的信号为s[t +τi(θ0)]， 其中 θ0为预设的信号入射角，对各通道信号延时补偿 τi(θ0) ， 使各路输出信号均为 s(t)， 再对N 路信号求和、平方和积分就得到，其中为信号功率。此时，信号如从 θ角入射，则第 i个阵元输出信号为 s[t +τi(θ)−τi(θ0)]，阵列输出为：

图 1 常规波束形成模型Fig.1 Conventional beamforming model

2 圆阵波束形成的信号模型

目前，阵元发射信号受到换能器谐振频点的限制，一般只能发射窄带信号。接收信号则根据信号的不同分为窄带信号模型和宽带信号模型，通常当目标回波信号带宽大于信号中心频率的1/10 时，被认为是宽带信号，由于发射信号和接收信号的波束形成原理相同[3]，文中以接收的信号模型为例说明。

被测目标到成像系统的距离一般是远大于成像系统的圆阵直径，阵元接收的信号属于远场窄带信号。假设回波信号为 s(t) ， 圆阵的阵元数为 N，圆阵接收的信号可以表示为：

式中， f0为 信号的中心频率， a (t) 为信号幅度， φ (t)为信号相位，其解析信号可以表示为：

假设接收信号的时延量为τ,圆阵接收信号为：

但是，窄带信号的幅度是慢变化的，则a (t−τ)=a(t)。可见，窄带信号时延只会引起信号复包络的相移变化，对信号的幅值的影响基本可以忽略。假设圆阵的参考阵元的接收信号为第i个阵元接收信号相对参考阵元接收信号的时延量为 τi，其信号可表示为：

圆阵在第n 次快拍接收到的信号模型可用矩阵表示为：

3 圆阵波束形成的数理模型

3.1 均匀加权波束形成

假设圆阵半径为 R， 在其圆周上等间距排布 N 个阵元，以原点作为参考点，依次确定信号到达各阵元的时延量，如图2 所示。

图 2 均匀圆阵阵元分布图Fig.2 Distribution of uniform circular array elements

其中， ∆R 为参考原点和第i 个阵元的波程差。则第i 个阵元和参考原点之间的相位差为：

整个均匀圆阵的波束响应可以表示为：

式中， si(θ,φ)为 第 i 个 阵元的接收信号， ωi为 第 i 个阵元接收信号的加权值。

在均匀圆阵上，采用均匀加权法进行波束形成，则有：

整个均匀圆阵的波束响应可以表示为：

3.2 旁瓣抑制加权波束形成

通过Dolph-Chebyshev 加权法抑制旁瓣高度，改善圆阵波束形成的性能，但是以上加权法只适用于线性阵列，故需要通过相位模式激励法将圆阵等效成线列阵，再进行后续的波束形成[4-5]。为了方便讨论波束性能，令则方位角为 φ，此时参考阵元 i的接收信号响应为：

定义矩阵 J 和 F分别为：

图 3 圆阵响应模式转换模型图Fig.3 Model diagram of response mode transformation of circular array

经过转化后，阵元信号响应矩阵为：

形成了线列阵的响应模式，具有范德蒙行列式的结构形式，在此基础上，可以按照Dolph-Chebyshev 加权法对输出信号 yi(n),i = −h,···,h进行波束形成，假设波束指向 φ1方 向，设置权向量 W：

则整个整理的波束响应为：

B(φ)

即，波束响应 可以转化为：

令 β为主瓣幅值和旁瓣幅值的比值，将式（22）表示成根值 z ∈[−1,1]且 还有 2h项的Chebyshev 多项式：

式 中， z = z0cos(u(φ)),z ∈[−z0,z0] 。 其 中， z0有T2h(z0)= β ＞1决定，可以进一步得到：

圆阵响应模式可表示为：

此时，圆阵波束响应模式已经等效为均匀线列阵Dolph-Chebyshev 多项式的响应模式。因为 u(φ) 和 φ1成线性关系，对于任一取值的 φ1， 圆阵波束响应 B(φ)主瓣宽度都相同，即任一波束指向方向的圆阵波束响应的主瓣宽度一致。对于任一方向接近的目标，其波束响应的强度相同，有利于均衡图像的背景，提高目标的图像分辨率。

值得注意的是上述阵列模式等效转换需要满足矩阵 BJF可逆的条件，即圆阵的阵元数圆阵的阵元间距必须小于回波信号的半波长。下面通过仿真实验经过Dolph-Chebyshev 加权后的旁瓣抑制效果。

4 圆阵波束形成仿真实验及性能分析

仿真实验1：假设均匀圆阵共有33 个各向同性的阵元，圆阵半径为 R， 接收信号的波长为 λ。要求波束指向 (0°,0°)方向，经过Matlab 仿真得到圆阵的波束指向性图，如图4 所示。

由图4 可知，圆阵三维波束形成图可以获得均匀加权圆阵的俯仰角和方位角信息，反映了均匀加权波束形成性能。但是，实孔径声呐成像系统和被测目标基本在同一深度，故方位角信息更值得关注。为了方便分析圆阵波束形成水平方向的指向性特征及讨论不同阵元数和半径的波束指向性能。

图 4 均匀加权圆阵三维波束形成图Fig.4 Uniform weighted circular array three-dimensional beamforming

图 5 均匀加权圆阵二维波束形成图Fig.5 Two-dimensional beamforming of uniformly weighted circular array

由图5 可知，均匀加权圆阵水平方向的半功率波束宽度1 1.5°， 最大旁瓣高度约为 − 8 dB，远离预成方向的旁瓣呈现递减趋势。

仿真实验3：其他实验假设条件不变，保持 R/λ=0.66 不变，改变阵元数 N，得到图6。

图 6 信号频率一定不同阵元数圆阵二维波束形成图Fig.6 Two-dimensional beamforming of circular arrays with different signal frequencies

由图6 可知，由于圆阵的非线性特点，在调整圆阵半径和接收信号波长比值不变的情况下，无论怎样改变圆阵阵元的数量，均匀加权圆阵的第一旁瓣高度总是约为-8 dB，明显高于线列阵的第一旁瓣-13.5 dB的高度。高旁瓣会引起高的虚警概率，造成大的弹道测量误差，应当在设计中采用适当的方法抑制旁瓣，保证测量精度。

仿真实验4：实验1 的其他假设条件不变，改变圆阵接收信号频率，即改变圆阵的半径和信号波长比，得到图7。

图 7 圆阵半径一定不同频率信号二维波束形成图Fig.7 Two-dimensional beamforming of signals with certain circular array radius and different frequencies

由图7 所示，在圆阵阵形固定的条件下，提高发射信号频率可以减小波束形成的主瓣宽度，改善波束指向性能，但是最大旁瓣高度仍为-8 dB 左右，虚警概率和轨迹测量依然较大，且提高发射信号频率会导致信号吸收损失的增强，限制了声呐系统的作用距离，需要慎重选择发射频率。

仿真实验5：设圆阵具有32 个均匀排布的阵元且阵元各向同性，半径波长比为0.75，主瓣幅值和旁瓣幅值比为 β =10（ − 20 dB ），令波束指向 ( 0°,0°)方向，俯仰角，最大相位模式 h = 5，则圆阵chebyshev 加权部分可以等效为11 个阵元均匀排布的线列阵，利用Matlab 仿真，分别得到均匀加权圆阵波束响应图和Dolph-Chebyshev 加权圆阵波束响应图，如图8 和图9 所示。

5 结 语

图 8 均匀加权圆阵二维波束形成图Fig.8 Two-dimensional beamforming diagram of uniform weighted circular array

图 9 Dolph-Chebyshev 加权圆阵二维波束形成图Fig.9 Dolph-Chebyshev weighted circular array 2-D beamforming

设计基于水下高速小目标末弹道测量的实孔径成像声呐应该优先考虑均匀排布的圆阵，保证没有观测盲区，提高观测效率。在圆阵做波束形成时，Dolph-Chebyshev 加权法能够有效抑制旁瓣高度，克服均匀加权圆阵波束形成旁瓣高的缺陷，同时能够保证在各个波束方向上，波束响应的幅值相同，为后续图像信号处理提供高质量的输入信号。但是，圆阵阵元数必须满足才能等效成线列阵的响应模式进行Dolph-Chebyshev 加权，否则常规波束形成的方法旁瓣抑制效果均不理想，故为了获得良好的目标回波图像采用Dolph-Chebyshev 加权时，必须合理设计圆阵半径、发射信号频率和阵元数量的关系以满足 BJF矩阵可逆。