杨 帆， 汪 洁， 林 海， 葛俊祥

(南京信息工程大学电子信息技术与装备研究院， 江苏南京 210044)

0 引言

目前，传统的针对阵列模型或导向矢量失配等问题提出的稳健算法为：空间谱估计算法、稀疏矩阵重构算法、对角加载(Diagnoal Loading，DL)算法、特征子空间算法(Eigenspace-based，ESB)和不确定集约束算法。对角加载波束形成算法通过在协方差矩阵上添加载平衡因子使得自适应和鲁棒性达到最佳平衡，但是因子大小的不易估计，改善效果不明显。特征子空间算法和投影方法本质上相同，已知期望和干扰信号导向矢量和特征向量，将所需信号投影到干扰互补子空间，从而得到最优的加权矢量，更好地消除了导向矢量失配的误差，但是在低信噪比时，信号功率和噪声功率大小非常接近，通过特征值大小分辨出信号和噪声子空间十分困难。稀疏矩阵算法是利用矩阵信号干扰空间的稀疏性，对信号的接收矩阵进行了重构，利用较低的快拍数恢复出原来的信号，以此降低算法计算量。不确定集约束的方法不同与以上对协方差矩阵的改进，可看成是对权值的优化约束问题。当采样协方差矩阵中包含了期望信号时，会出现信号“自消”现象，传统的稳健自适应算法在期望导向量严重失配或高信噪比时，性能会下降。

为了提高自适应波束算法的稳健性，文献[10]指出期望信号在空域中的功率分布中有残留噪声，消除残留噪声可以重构出更加精确的期望导向矢量和干扰加噪声矩阵协方差，但随着信噪比的提高，干扰抑制减弱，算法性能变差。文献[11-12]提出了对传统的波束形成算法进行加权，形成改进算法，虽然相比传统算法副瓣电平大大降低，小样本情况下方向图稳定性也得到了极大的改善。但是由于没有有效地提取出干扰噪声信号协方差矩阵，故容易产生“自消”。文献[13]在基于标准Capon波束形成器对期望矢量先进行波束预处理，降低了导向矢量失配的影响，但是对干扰移动和低样拍等情况下的算法性能没有改善。文献[14]在文献[13]的基础上展宽了干扰零陷，但是展宽干扰零陷引入的计算误差会影响算法的性能，并且算法需要对所有空域进行积分，运算量大。文献[15]基于Capon的空间谱估计算法，重构出精度较高的干扰噪声协方差矩阵和期望信号，降低了失配对算法的影响。为了降低小样本对波束性能的影响，文献[16]提出了一种基于稀疏干扰矩阵的稳健波束算法，该算法利用信号在空域的稀疏性，仅设置扰动参量对干扰方向的谱估计进行修正并重构，性能良好且计算复杂度低。文献[17-18]提出了缺少先验信息的协方差矩阵重构算法，在未知阵列结构的情况下性能良好，但是在估计期望信号导向矢量的时候不确定约束集的大小很难选择，增大了算法的复杂性。文献[19]使用去相关性算法得到新的干扰噪声协方差矩阵，有效地降低了干扰和期望信号的相关性，但该算法提高波束形成的性能效果有限。文献[20]针对干扰运动的抑制问题提出了一种零陷加宽的方法。

本文利用残留噪声功率消除重构期望导向矢量和干扰协方差矩阵减小算法计算量，再基于改进的MVDR算法有效地降低了副瓣和样本对算法性能的影响，最后在干扰方向展宽零陷解决了干扰波动的问题。

1 阵列信号模型

()=()+()=

()+()

(1)

式中，=[(),(),…,()]为×(+1)维阵列流型矩阵，[]为矩阵转置；()为采样信号波形，为快拍数；()为期望导向矢量，()为干扰导向矢量，()为独立同分布的高斯白噪声向量。

元均匀线阵列天线的接收信号协方差矩阵为

=[()()]=

+=

(2)

根据最大化输出信干噪比(MSINR)准则，即

(3)

标准Capon波束形成可以表示期望信号无失真响应下，阵列输出功率最小的最优化问题，即最小方差无失真响应：

(4)

由此可得最优权值为

(5)

(6)

2 本文提出的算法2.1 期望信号矢量估计

根据文献[6]，MUSIC谱估计算法功率谱表示为

(7)

(8)

式(8)可表明，每个回波信号辐射区域中的功率是有用信号功率与空域噪声功率之和，并且回波信号中包含的噪声功率为自然噪声功率的1。所以可推出空域噪声功率可以利用功率谱在非信号区间进行估计：

(9)

式中，为非期望信号角度区域，为在的采样点，为采样点个数。根据式(8)、式(9)，可获得实际的噪声功率和噪声协方差矩阵：

(10)

为了降低计算量，可先结合文献[13]对波束预处理，然后消除空域噪声成分后重构期望信号，其最大特征值对应的特征矢量等于期望信号导向矢量。

2.2 干扰噪声协方差矩阵估计

(()+())≈()

(11)

由式(11)可得

(12)

(13)

(14)

2.3 低副瓣加权和特征干扰相消

虽然重构了导向矢量和干扰噪声协方差矩阵，传统的RCB算法在导向矢量失配的情况下具有稳健性，但是当小样本和副瓣较高时，接收矩阵出现低秩和奇异解，从而导致算法产生较大的误差。本文利用类似输出方向图加窗的方法，对最优权值进行加窗，使其随单元分布成锥形减弱，从而达到副瓣电平降低的目的。现在采用Dolph-Chebyshev加权函数对加窗，加窗后的期望信号方向矢量为

()=·()

(15)

另外，通过文献[12]提出的特征干扰相消器，本文对干扰信号矢量建立新的线性约束，这样提高干扰信号抑制和噪声抑制自由度，优点是在小样本数下，减小噪声的不稳定对方向图主瓣形成的影响。对应线性最小方差约束(LCMV)算法，特征干扰相消器可称为线性特征干扰相消器(LCEC)。

LCEC的准则为

(16)

由拉格朗日乘子法求得式(16)的最佳权值为

=()(()())

(17)

2.4 干扰零陷展宽

本文结合了文献[12]算法对权值进行修正并对干扰输出功率进行参数约束，实现干扰区域零陷加宽。

要使上一节得到的权值更加接近最终权值，对权值进行二次约束，可得约束方程如下：

(18)

式中，为极小的自定义约束且保证>0。利用Lagrange乘子法求解式(18)得目标函数为

(19)

式中，为拉格朗日乘数。

对式(19)关于求偏导并使偏导等于零解得

(20)

最优拉格朗日乘子可以通过求解约束方程获得，将式(20)代入约束条件式(19)中，可得关于的方程：

(21)

从式(14)可知理想干扰协方差矩阵表达式为

(22)

(23)

(24)

可得的取值范围为

(25)

(26)

3 仿真实验分析

实验1： 不同约束参数正负拉格朗日乘数的算法分析

当输入=10 dB，干噪比均为30 dB，采样拍数=50。图1、图2的约束参数分别取较大= 10和较小= 10时波束对比。

图1 ε=10-2的方向图差异

图2 ε=10-12的方向图差异

从图1可以看出，当取值较大时，正负乘数都能实现期望导向失配修正和零陷变宽，且主瓣和副瓣吻合比较一致，波动较小。随着取值越来越小，即约束参数精度越高时，如图2所示，干扰零陷的深度会加深，功率转移到主波束方向，导致主瓣相应变宽和副瓣上升，算法性能略有下降。综上可知，约束参数变化对算法性能虽然有影响，但该算法仍然拥有良好的抗干扰运动和抗系统误差的鲁棒性。

实验2： 导向矢量失配时不同算法的波束图分析

图3 导向矢量失配下的归一波束图

图4 高信噪比下的归一波束图

图5 低样拍数下的归一波束图

从图3可以看出，在期望矢量失配时，几种稳健算法都能够满足失配下的校准和干扰抑制，深度有所差别。文献[14]和本文的算法不仅在干扰位置产生较深的零陷，并且对零陷有效的展宽，有利于抗干扰移动。但是从计算量来说，前者进行空间谱估计时使用了大量积分运算，而后者的计算主要是矩阵分解，计算量大大减小，所以本文算法更有利于实际应用。并且，由于本文使用了副瓣控制加权，在主瓣宽未发生大的扩展的前提下副瓣高度明显比其他算法低，大大地提高了算法的性能。随着信噪比的提高，如图4可知文献[10]算法由于信噪比过大，使()的近似估计误差会增大，干扰信号无法正常抑制，性能急剧下降。在低样拍数的情况下，如图5所示，文献[14]算法的样拍数低于阵元数后性能急剧的恶化，导向矢量发生偏移，这是因为在低样拍数下接收信号协方差矩阵发生误差，无法再近似为实际信号计算。综上分析，本文算法解决了在高信噪比和低样拍数的情况下性能恶化的问题。

实验3： 不同输入SNR下的输出SINR分析

图6 不同输入信噪比下的输出信干噪比

从图6可以看出，本文算法使用的展宽零陷算法带来了约束计算的误差，会导致性能的下降，但是用副瓣加权的算法可以优化旁瓣，从而弥补性能的下降，使其远远优于其他算法。并且，本文算法利用特征干扰相消算法重构干扰矩阵从根本上解决了在高信噪比下干扰抑制下降的问题，使输出SINR与理论情况相差很小。

实验4： 不同采样拍数下的输出SINR分析

图7 不同快拍数下的输出信干噪比

从图7可知，本文算法明显降低了快拍数对算法稳健性的影响，在极低的样本数下就能快速达到收敛，并且输出SINR接近理论值。相比本文算法，文献[14]达到性能收敛时样本数至少需要20。文献[12]和文献[20]都只考虑了单方面的因素，性能一般。理论波束算法在低样拍数时无法得到完整的样本协方差矩阵，所以算法性能出现急剧变化，直到快拍数为50才趋于收敛。综上可得，本文算法在不降低算法性能的前提下，大大提高了算法收敛的速度。

实验5： 随失配角度变化的输出SINR分析

当快拍数=50，输入=10 dB，=30 dB固定不变。图8表示期望信号的失配角度在[-12:1:12]范围内，输出SINR的曲线对比。

图8 不同失配角度下的输出信干噪比

由图8可得，本文的算法失配角度在[-8:1:8]有着良好的稳健性，使用了矩阵相消算法的文献[12]由于未对期望信号矢量进行估计，所以在角度失配方面稳健性差。文献[20]在期望信号失配时输出SINR发生了剧烈起伏，并且远离性能最优值。本文算法在导向矢量失配的情况下稳健性与空间谱估计算法近似，且宽度略宽。

实验6： 不同算法的复杂度对比分析

算法的运算速度是数字波束形成器实际应用的关键，为了减少FPGA和AD采样的压力，快速的自适应算法是对雷达工作性能的保证，所有算法在MATLAB上进行仿真，性能配备主频为3.79 GHz-CPU。在其他仿真条件不变的情况下，快拍数=100，输入=10 dB，干扰噪声比为=30 dB。当将期望回波信号和强干扰信号预估范围设置为10°，为信号空域内的总采样点数，为阵元个数，为空域回波信号总数。算法总运行的时间由100次独立的蒙特卡洛实验累计得到。

由表1可知，本文算法和文献[10]算法主要的运算是矩阵特征分解，因此处于同一个量级，文献[12]主要是矩阵求逆运算，而文献[14]算法在重

表1 不同算法的复杂度对比

构干扰噪声矩阵时比本文算法增多了特征分解和空间积分运算，由于≫，所以当阵列数较小时，L的影响远远大于N，故大大地增加了算法的运算时间。

综上所述，本文算法在阵元数较大的时候，保证了算法的复杂度，在阵元数较小的情况下，减少算法的计算量。

4 结束语

本文算法首先基于残留噪声消除算法对期望信号矢量进行了估计，然后利用干扰功率估计算法构造了新的干扰噪声协方差矩阵，此算法较空间积分算法而言，性能更优，且计算量更小；接着为了降低零陷展宽算法误差对波束性能的影响，对最优权值进行副瓣加权降低波束的旁瓣；并结合特征干扰相消，消除快拍数对波束形成算法的影响；最后在构造干扰矩阵零陷之后，使用二次约束准则实现零陷加宽，虽然干扰区域展宽会使得零陷深度降低，但是提高了算法的鲁棒性，更有利于实际应用。仿真结果证明了该方法在降低协方差矩阵构造计算量的同时，克服了比较算法在高信噪比和低样拍数的情况中性能下降的问题，提升了导向矢量失配和抗运动干扰的能力。本文算法在高信噪比和强干扰等恶劣环境中都拥有良好的鲁棒性，通过降低快拍数的影响和减少计算量，使算法在雷达休止期拥有更快的波束校准和波束形成。