双约束条件的时域MVDR 自适应波束形成方法

董 晋，胡 鹏，冯金鹿

(中国船舶集团有限公司第七一五研究所，浙江 杭州 310023)

0 引 言

波束形成技术是阵列信号处理的关键技术之一，在雷达、声呐信号处理领域有着广泛的应用。在主动声呐阵列信号处理中，常规波束形成方法（CBF）主瓣较宽、旁瓣较高，空间分辨力的不足导致CBF 对多目标、弱目标的检测性能不足。自适应波束形成方法拥有较高的空间分辨力，在多目标和弱目标检测上体现了良好的性能，因此研究自适应检测方法具有重要的意义[1]。

MVDR[2]是已经较为成熟的一种自适应波束形成方法，传统MVDR 方法为频域MVDR，其缺陷是需要较长的数据长度进行累积来保证最优权向量[2]的稳定估计。文献[3]提到时域MVDR 方法TMVDR，该方法在构造时域解析信号的基础上，对每一路时域解析信号引入复权，提高了MVDR 最优权向量估计的稳定性。但在大孔径主动声呐信号处理的应用环境中，回波信号脉宽较短、有效快拍次数较少以及阵列可能存在的阵型失配等问题，必然会使MVDR 滤波器的最优权向量估计存在误差，最终导致其输出功率失真，降低了MVDR 的工程实用价值。本文在TMVDR 的基础上加入了误差恒定和功率恒定2 个约束条件[4–5]，实现了一种双约束条件的时域MVDR 自适应波束形成方法DCTMVDR（Double constraint Time-domain MVDR），本方法对原本MVDR 无约束的最优权输出进行了误差恒定和功率恒定的约束补偿，在保证较高空间分辨率的同时有效解决了MVDR 在恶劣环境下性能下降的问题，提高了算法鲁棒性。同时本方法的协方差矩阵估计采用对角加载法[6]，增加了协方差矩阵的稳定性，进一步减少误差。本方法可以在单个快拍下实现最优权向量的稳定估计，在实际工程中也能体现出较好的性能。

1 TMVDR 原理及主动处理的应用

假设主动声呐M个阵元接受到的时域数据为：

其中，n=1,2,···,N表示时间采样点，N表示一个处理拍快的数据长度。

m=1,2,···,M

假设为阵元序号，第m路阵元时域数据可以表示为：

其中，A表示信号幅值，n0表示回波信号到达时间点，Na表示发射信号脉宽长度，a表示发射信号：

其中，*表示共轭，⊙表示Hadamard 积。

则可得到阵列相关检测输出：

计算本次快拍下相关检测输出的协方差矩阵：

其中，R为M×M矩阵。

假设TMVDR 滤波器的权矢量为w，则其输出可表示为：

其中，H 表示共轭转置。

TMVDR 的约束问题表达式[3]如下：

其中，I表示M×1维单位向量。

用拉格朗日乘数法求解式（11）可解得TMVDR 滤波器的最优权矢量：

代入式（10）则可得到TMVDR 滤波器的最优权向量：

需要注意的是由于协方差矩阵的估计很难得到满秩且稳定的M×M维矩阵，采用对角加载方法可以缓解这一问题。

2 DCTMVDR 原理

假设式（12）中的期望值为单位向量I，当协防差矩阵R的估计产生误差时则必然使得输出的权矢量不是最优值，导致性能下降。

假设存在一个补偿权向量b，对R进行误差补偿后可以使wopt的输出达到最优。以b¯=I作为期望值，则以‖b−I‖2≤ε视为误差上限为 ε的误差恒定约束条件，以‖b‖2=M视为能量恒定约束条件，因此DCTMVDR的优化问题可归纳为如下表达式：

其中：ε代表b与I之间误差范围。我们将给定误差范围的b的值域视为一个以I为圆心、半径为的球体，DCTMVDR 的原理可以表示为搜索该球体以寻找一个同时满足2 个约束条件的理想值bˆ使得算法输出权向量wopt的值为最优。

在不考虑第二约束条件‖b‖2=M的前提下，利用拉格朗日算法求解式（14），可以得到单约束条件下的解：

其中：IM×M为单位对角阵。则结合约束条件一‖b−I‖2≤ε，可得到λ的函数g(λ)：

对R进行特征值分解，可以得到：

其中 Λ表示由特征值构建的对角矩阵，其中特征值按从大到小排序，即γ1>γ2···>γM，U表示对应各个特征值的特征向量矩阵。令z=UHI，则式（16）可化简为：

利用式（18）可求得 λ。需要注意的是，因为‖I‖2>ε，所以由式（16）可得到g(λ)是λ的单调递减函数，且g(0)>ε，故λ≠0。由式（18）得到0<ε，可知对于式（18）存在 λ的唯一解（λ>0）。由排列方式可知 γM，γ1分别为特征值 γ的最小值和最大值，由式（18）可解得 λ取值范围的上下界。去掉式（18）分母中的1，可以得到 λ的另一个上界，故λ取值范围为：

利用拉格朗日算法求解式（14），可以得到双约束条件下的最优解：

将式（18）得到的 λ代入得到DCTMVDR 的理想补偿值。

则DCTMVDR 滤波器的最优权矢量可以表示为：

DCTMVDR 滤波器的最优输出可以表示为：

3 数值仿真验证

通过数值仿真进行阵列输出数据模拟，分别在理想条件和阵型失配条件下，对CBF，TMVDR 和DCTMVDR 三种方法进行比较，验证DCTMVDR 的有效性。

1）理想条件下的波束响应

仿真条件：进行32 元等间距线阵数据模拟，发射信号模拟频率为1.2 kHz 的调频信号，带宽为500 Hz，脉宽为0.2 s，采样率为15 kHz，仿真单个目标回波，目标方位81°，输入信噪比30 dB。

图1 仿真结果表明，在理想条件下，相比较与CBF，TMVDR 和DCTMVDR 拥有更好的空间分辨力和较低的旁瓣级。而TMVDR 和DCTMVDR 体现的性能基本相同，其中DCTMVDR 的旁瓣级略低于TMVDR。

图1 理想条件下CBF，TMVDR，DCTMVDR 的波束输出Fig.1 Beam output of CBF，TMVDR and DCTMVDR under ideal condition

2）阵型失配条件下的波束响应

仿真条件：进行80 元等间距线阵数据模拟，发射信号模拟频率为2.5 kHz 的调频信号，带宽为500 Hz，脉宽为0.2 s，采样率为15 kHz，仿真2 个目标回波，目标1 方位81°，输入信噪比10 dB，目标2 方位86°，输入信噪比15 dB。进行阵型失配模拟，在每个通道加入随机的相位扰动。

图2 仿真结果表明，在条件恶劣的情况下，TMVDR 会发生输出功率失真的现象，相比之下DCTMVDR 输出功率稳定，在保证输出能量不失真的前提下提高了时间分辨率，体现了最优的性能。

图2 阵型失配条件下TMVDR，DCTMVDR 的波束输出Fig.2 Beam output of TMVDR and DCTMVDR under the condition of array mismatch

4 实验数据分析

选取某型号大型拖线阵声呐实验数据中的1 个数据段作为数据分析样本。通过对CBF，TMVDR 和DCTMVDR 三种方法进行性能比较，验证DCTMVDR 的有效性。

实验数据说明，选用实验数据中某一周期数据，该段数据处于阵型轻微失配环境，截取包括宽带信号目标回波在内的数据段，数据段长度为2 s，发射信号为宽带调频信号。

图3～图5 为3 种方法的方位历程图。由图可见，回波目标位于方位39°、距离0.85 s 的位置（图中已标注），位于80°左右有1 个强干扰目标。CBF 检测目标较为模糊，TMVDR 检测目标失败，DCTMVDR 检测目标清晰可见。

图3 CBF 方位历程图Fig.3 Azimuth history chart of CBF

图4 TMVDR 方位历程图Fig.4 Azimuth history chart of TMVDR

图5 DCTMVDR 方位历程图Fig.5 Azimuth history chart of DCTMVDR

分析结果可以表明，CBF 的空间分辨率较低，受到强干扰影响，检测出的目标较为模糊。由于环境恶劣，TMVDR 性能下降，没有能够检测出目标。相比较之下，DCTMVDR 不仅拥有更好的空间分辨率，有效地抑制了强干扰的旁瓣，也体现了稳健良好的目标检测性能。

5 结 语

TMVDR 方法可以有效地提高声呐波束输出的空间分辨率、降低旁瓣级，在大孔径主动声呐的应用环境下，由于单次快拍协方差矩阵估计的不稳定性以及可能存在的阵型失配等不确定因素，其存在输出功率失真的问题，导致其性能下降。相对于TMVDR，DCTMVDR 在一定程度上能够有效地避免由于最优权向量估计误差导致的性能下降问题，在估计最优输出的同时拥有更好的鲁棒性。数值仿真和实验数据验证了DCTMVDR 方法的有效性。