零陷展宽对角载入算法

零陷展宽对角载入算法
Null broadening and diagonal loading algorithm
CHEN Tao，SUI Li?li ，WU Huan?xin
(College of Information and Communication Engineering ，Harbin Engineering University ，Harbin 150001，China )
：The null broadening and diagonal loading algorithm is proposed ，which not only addresses the problem of narrow nulls as the interference in the fast movement condition ，but also solves the problem of the bad algorithm stability as the error covariance matrix and steering vector errors existing. Meanwhile ，diagonal loading factor values of the algorithm are determined based on the relationship between sample covariance matrix and diagonal loading factor. The
computer simulation results show that the algorithm has good robustness ，as well as broadening nulls.
Keywords：MVDR；steering vector ；sample covariance matrix ；null broadening ；diagonal loading loading factor
0 引言波束形成技术被广泛地应用在雷达、无线通信领域[1?4] 。

其中，MVDR自适应波束形成算法，因具有良好的分辨率和干扰
抑制能力，被广泛的使用。

MVD算法一般假定期望信号的导向矢量是精确已知的，但在实际中，由于指向误差、阵元位置误差以及各阵元特性不一致等因素影响，假设的期望信号导向矢量可能不完全准确，与真实值之间有失配误差，从而导致MVDRI法
性能急剧下降⑸；同时在窄带条件下，MVDf算法在干扰方向形成的零陷较窄，当出现快速移动的干扰时，由于自适应权值的收敛速度比不上干扰的移动速度，在干扰方向不能形成有效的零陷，从而使MVD算法失效。

采用零陷展宽算法一一协方差锥化（CMT）[6?7] 法在干扰方向形成较宽的零陷，能有效地抑制快速移动的干扰。

但是，由于实际环境中存在各种误差，协方差锥化法在实际使用中往往不能得到预期的结果。

因此，寻找出一种既对各种误差有较强稳健性又对快速运动的干扰有抑制性的自适应波束形成算法是非常紧迫的。

近几十年来，为了增强在各种误差下自适应波束形成算法的
稳健性[8] ，学者们进行了大量的研究。

其中，特征空间（ESB）算法、线性约束最小方差（LCMV算法以及对角载入（LSMI）法是最具有代表性的算法。

ESB算法[9]收敛速度比较快，但它需要准确估计出信号子空间的维数；LCMV算法通过适当的约束条件，控制某些方向波束的增益，使自适应波束满足一定的稳健条件，但它只适合用在观察方向出现失配时的情况；LSMI算法是
一种简单有效的方法，但是载入因子难以确定[10?11] 。

本文针对MVD算法在各种误差下的性能下降及干扰快速变化时算法失效的问题，将对角载入算法与零陷展宽算法结合在一起，提出了零陷展宽对角载入算法。

同时根据载入因子与采样协方差矩阵间的关系来确定载入因子。

零陷展宽对角载入算法利用对角载入算法提高波束形成对系统误差的稳健性，同时利用零陷展宽算法能在干扰方向形成较宽的零陷，解决了干扰快速移动时，算法失效的问题。

1 传统算法描述
1.1 阵列信号模型
在窄带条件下，考虑［M］元均匀线阵，［D］个互不相关的信
号。

其中，期望信号波达方向是［e 0, ］［D-1］个干扰信号的波达方向分别是［{ e 1, e 2,…，e D-1 },］有［M>D］。

第［i］个阵元端接收的信号为：
［xl (t) =i=0D-1si (t) e-j2 n 入(l-1 ) dsin e i+nl (t)］
(1)
式中：［si (t)］为信号的复包络；［入］为信号的波长；［d］
为阵元间距；［nl (t)］为第［l］个阵元上均值为0、方差［(T 2n］
为1 的白噪声。

阵列接收的信号表示成向量的形式为：
［X( k =AS( k +n( k =s0( k a(e 0 +i ( k +n( k ］
(2
式中：［X(k =［X1 (k , X2( k ,…, XM(k ］T］为接收向量；
［A=］［［a (e 0 , a( e 1 ,…, a(e D-1 ］］为阵
列
流型；［a (B i )］为信号［i］的导向矢量；［S=［sO , si,…，sD-1］T］为信号源；［a ( 6 0)］为期望信号的导向矢量；［s0 (k)］为期望信号；［i (k)］为干扰信号向量；［n (k)］为噪声向量。

阵
列输出为：
［y ( k)=wHX( k)］(3)
式中：［w=［w1 , w2,…，wM］T］为权重向量，［(?) T］表示
矩阵的转置，［(?) H］表示矩阵的共轭转置。

1.2 MVDR波束形成算法
MVD波束形成算法使干扰和噪声受到抑制而在阵列输出中的功率最小，又能使期望方向上的信号功率保持不变。

其代价函数为：［minwHRxw subject to wHa (6 0) =1］(4)
式中［Rx=EX( k) X (k) H］为接收信号的协方差矩阵。

利用Lagrange 乘子算法求出最优权重向量为：
［wopt=R-1xa ( 6 0)a( 6 0)HR-1xa(6 0)］( 5)
在实际应用中，［Rx］一般难以得到，往往用采样协方差矩阵［Rx］代替［Rx ,］假设快拍数为［K ,］则：
［Rx=1Ki=1KX(i )XH(i )］(6)
此时MVDF波束形成算法又称采样协方差矩阵求逆(SMI)算
法，则SMI算法的权重向量为：
［wSMI=R-1xa(6 0)a(6 0)HR-1xa(6 0)］(7)
从式(5)可看出在MVD算法中，需要精确知道期望信号的
导向矢量以及接收信号的协方差矩阵。

然而，在实际中，由于指向误差、阵元位置误差以及各阵元特性不一致等因素影响，使得导向矢量存在误差，同样，在实际应用中，接收数据是有限长的，因此采样协方差矩阵[Rx] 代替接收信号的协方差矩阵[Rx] 也会产生一定的误差。

导向矢量和协方差矩阵误差使得MVD算法的
性能下降，尤其在接收数据中包含期望信号时，算法的稳健性更差。

除此之外，在窄带条件下，MVD算法在干扰方向形成的零
陷较窄，当出现快速移动的干扰时，由于自适应权值的收敛速度比不上干扰的移动速度，在干扰方向不能形成有效地零陷，从而使MVDf算法失效。

针对MVDf算法在存在导向矢量误差、协方差矩阵误差，以及存在快速运动的干扰时，算法稳健性变差和干扰零陷过窄的情况，本文提出了一种稳健的MVDf算法一一零陷展宽对角加载算法。

2稳健的MVD算法
2.1 零陷展宽对角载入算法零陷展宽对角载入算法是把零陷展宽
算法和对角载入
（LSMI）算法相结合形成的一种稳健的自适应波束形成算法。

该
算法解决了干扰在快速运动时，干扰零陷过窄的问题，又解决了协方差矩阵误差和导向矢量误差存在时，算法稳定性变差的问题。

因为该算法具有零陷展宽算法能在干扰方向形成较宽零陷的特点，又具有LSMI 算法具有较高的稳健性的特点。

LSMI算法是用对角载入的协方差矩阵代替常规的采样协方
差矩阵，即：
[Rd匸Rx+ E I] (8)
式中[E ]为对角载入因子。

CMT零陷展宽算法是通过一个锥化矩阵[TMZ]对采样协方差
矩阵[Rx] 进行扩展。

扩展后的采样协方差矩阵为：
[RMZ=Rx?TMZ]( 9)
式中：“⑺”为Hadamard积，[TMZ]第[m]行[n]列的元素
可表示为：
[[TMZ] mn=sin (( m-n) △)( m-n) △ ] (10)
式中[△]为零限展宽参数。

通过式( 8)和式( 9)可得到零陷展宽对角载入算法的协方
差矩阵[RC-T，] 即：
[RC-T=RMZ+E I] (11)
零陷展宽对角载入算法权向量为：
[wCMT-LSMI=R-1C-Ta( 6 0) a ( 6 0) HR-1C-Ta ( 6 0)= (RMZ E I ) -1a ( 6 0) a ( 6 0) H(RMZ E I ) -1a ( 6 0) ] (12)
从式(12)中可以知道零陷展宽对角载入算法需要知道期望
信号的导向矢量、扩展后的采样协方差矩阵[RMZ]以及对角载入因子[E。

]但是对角载入因子[E ] 一般很难确定，而且对角载入因子的大小会影响自适应波束的效果。

当载入因子过小时，波束旁瓣的高度不能有效的控制；当载入因子过大时，对干扰抑制效果的灵敏度会有所影响。

针对以上情况，本文提出一种动态的确定对角载入的方法。

2.2 确定对角载入因子
首先，对协方差矩阵[Rx] 进行特征分解，可写成如下形式：
[Rx=i=1D(T 2iuiuHi+ 彷2nl=A A AH+(r 2nI] (13)
式中：[彷2i]代表第[i]个信号的功率，[ui]是[Rx]特征值对应的特征向量；[A =diag[彷21,彷22，…，彷2D]]为对应信号的功率构成的对角矩阵。

采样协方差矩阵[Rx] 和真实的协方差矩阵[Rx] 之间的关系：
[R X=R X+E E] (14)
式中：[E]是均值为0，方差为1的矩阵；[£ ]是个大于零的常数，表示采样协方差矩阵[Rx] 的估计误差。

对角载入后的协方差矩阵为：
[Rdl=Rx+ £E+E I] (15)
则对角载入后协方差矩阵的逆矩阵为：
[R-1dl= (Rx+£E+E I ) -1=(Rx+E I ) -1[I+ £E(Rx+E I)
-1]-1] (16)
又假设[ £ E 式(17) 中大括号部分导致了自适应波
束性能降低。

为了得到最优权，希望有：
[£ E +^ 2n0]，根据式(18)和(19)可得：
[£ w E
图1理想条件和实际条件下不同载入值LSMI算法
图 2 不同算法阵列方向图比较
仿真三：不同快拍［K］及输入信噪比时阵列的输出信干噪比。

考虑训练数据中包含期望信号时，信号到达角无偏差，和偏差为1°时阵列的输出信干噪比SINR (［6 s=1?］代表偏差为1 °)。

仿真结果如图 3 所示。

图3不同快拍［K］和不同信噪比SNR寸阵列输出SINR
从图3 (a)中可以看出，当没有误差时，本文算法和MVDR 算法的变化趋势相同，即快拍数越少输出的信干噪比SINR与理
想SINR之间的差距越大，越大越接近于理想值，当两方法快拍数［K］相同时，本文算法的输出SINR高于MVDRT法。

当存在误差时，MVDRT法的输出的SINR变小，而对本文算法的输出SINR 影响较小，可见对角载入算法提高了MVDR^法的稳健性，尤其
在快拍数较少及期望信号存在误差时，性能改善非常明显。

图3 ( b)中期望信号存在偏差时，随着信噪比SNR的增加，MVD算法的输出信干噪比性能下降越明显，而本文算法对信号方向偏差的敏感度较低，输出的信干噪比更接近理想值。

可见本文算法不仅能够提高MVD波束形成器在小快拍、高信噪比时的稳健性，而且能够提高它对导向误差的稳健性。

4 结论
针对MVDF波束形成算法在实际应用环境中性能下降的问题，本文提出了零陷展宽对角载入算法，该算法把零陷展宽算法和LSMI 算法结合在一起，然后根据协方差矩阵的估计误差来自适应地确定对角载入因子。

仿真结果表明：该算法在低快拍数，期望信号存在偏差即导向矢量存在误差的情况下仍具有较好的波束形成性能，且在干扰方向上形成较宽的零陷。

它能改善干扰的入射方向变化过快，可能使得干扰移出天线阵方向图零陷位置而无法得到有效的抑制的情况。

总之，该算法是一种稳健的且性能优越的波束形成算法。