一种抗阵列流型误差的波束形成零陷展宽方法

一种抗阵列流型误差的波束形成零陷展宽方法
赵宇;李文兴;毛晓军;张宁
【摘要】针对现有波束形成零陷展宽方法在阵列流型误差情况下性能严重下降的问题,提出了一种抗阵列流型误差的零陷展宽方法.该方法首先利用不确定集优化方法,对预定区域内的干扰导向矢量进行校正,再结合协方差矩阵锐化技术构造一个零陷展宽投影矩阵,对采样协方差矩阵做投影处理,最后用投影处理后的协方差矩阵进行波束形成.仿真结果表明:所提方法在阵列流型误差的情况下,与现有波束形成零陷展宽方法相比具有更好的零陷展宽效果和输出性能,且计算复杂度很低,是一种先进的稳健波束形成零陷展宽方法.%To address the serious degradation of the performance of existing null broadening beamforming when array calibration errors exist, a null broadening beamforming method against array calibration errors was proposed.In the method, uncertain quadratic optimization was utilized to correct the steering vectors of interferences in the predeter-mined area.Covariance matrix taper technology was then combined to construct a null broadening projection ma-trix, which was used to implement projection processing for the sampling covariance.Finally, the covariance matrix that experienced projection processing was used for beamforming.In case of array calibration errors, the proposed method presented better effect and output performances than the present null broadening beamforming method.Mo-reover, the computational complexity of the proposed method is low.The proposed technique is an advanced and robust null broadening beamforming method.
【期刊名称】《哈尔滨工程大学学报》
【年(卷),期】2018(039)001
【总页数】6页(P163-168)
【关键词】自适应波束形成;零陷展宽;阵列流型误差;投影变换;协方差矩阵锐化
【作者】赵宇;李文兴;毛晓军;张宁
【作者单位】哈尔滨工程大学信息与通信工程学院,黑龙江哈尔滨150001;哈尔滨
工程大学信息与通信工程学院,黑龙江哈尔滨150001;哈尔滨工程大学信息与通信
工程学院,黑龙江哈尔滨150001;北京机电工程研究所,北京100074
【正文语种】中文
【中图分类】TN911.7
自适应波束形成技术通过调整加权矢量，能够在干扰方向形成零陷，抑制干扰和噪声，被广泛应用于雷达、声呐、射电天文学、移动通信等领域[1-3]。

当干扰快速
运动或者天线平台振动时，自适应加权矢量的训练数据与应用数据可能会失配，导致干扰移出零陷位置，造成自适应波束形成器的性能严重下降，在干扰方向附近形成具有一定宽度的零陷是解决该问题的一种有效方法[4-5]。

协方差矩阵锐化(covariance matrix taper, CMT)[6-8]是一种经典的零陷展宽方法，其计算复杂度很低，但是零陷深度较浅。

武思军等[9]提出了干扰方向无关的自适
应零陷展宽技术，仍属于CMT方法。

对预定区域幅度响应进行二次约束(quadratic constraint sector suppressed, QCSS)[10- 11]也能够展宽零陷，QCSS方法与CMT方法相比，能够获得更深的零陷，但是二次约束的求解过程较
为复杂。

Amar等[12]在QCSS方法的基础上，提出了线性约束零陷展宽(linear constraint sector suppressed, LCSS)方法，将非线性的二次约束转化为线性约束，简化了求解过程。

李文兴等[13-14]提出了基于投影变换和对角加载的零陷展宽(projection and diagonal loading null broadening beamforming, PDNBB)方法，用投影矩阵修正采样协方差矩阵。

PDNBB和LCSS方法都通过构造导向矢量
相关矩阵提高了零陷展宽性能，因此性能都十分依赖于正确的阵列流型。

阵列天线在实际应用中，很可能会存在阵元位置误差、阵元互耦等非理想因素，导致假想的阵列流型与实际出现阵列流型误差。

在这种情况下，现有零陷展宽方法的性能会严重下降，甚至完全失效。

针对这一问题，本文提出一种抗阵列流型误差的波束形成零陷展宽(null broadening beamforming against array calibration errors, NBACE)方法，NBACE方法先利用不确定集优化方法对预定区域内的干扰导向矢量进行校正，再结合CMT方法构造一个零陷展宽投影矩阵，并对采样协方差矩阵做投影处理，最后用投影处理后的协方差矩阵进行波束形成，仿真结果验证了方法的有效性。

1 现有波束形成零陷展宽方法
1.1 阵列信号模型
考虑远场情况下的非相干窄带信号，天线为N元均匀直线阵，阵元间距为半波长，阵列天线接收数据X(t)为
X(t)=AS(t)+N(t)
(1)
式中：A=[a1(θ),a2(θ),…,aL(θ)]为阵列导向矢量矩阵，L为信源数，S(t)为信号复
包络，N(t)为自适应天线的噪声。

阵列天线接收数据的协方差矩阵为
R=E{X(t)XH(t)}=ARsAH+I
(2)
式中：Rs表示信号协方差矩阵是噪声功率。

实际计算中，阵列天线接收数据的协方差矩阵常用有限次快拍数据估计值代替，即
(3)
式中 K为快拍数。

标准MVDR波束形成器的自适应权矢量为
Wopt=
(4)
式中a(θd)为期望信号导向矢量。

对阵列天线波束形成器输出性能的衡量指标采用输出信干噪比(signal to interference plus noise ratio, SINR)，即阵列天线输出的期望信号功率与输出的干扰加噪声功率和的比值：
SINR=
(5)
式中表示期望信号输入功率，Ri+n表示干扰加噪声协方差矩阵。

为了进一步提高波束形成器的输出性能，可以采用对角加载技术，以下方法均默认采用对角加载技术进行优化。

1.2 CMT方法
CMT方法[6]的基本思想是用一束非相关的虚拟干扰源替代原有单个干扰源，用以在干扰信号方位附近形成具有一定宽度的零陷，具体实现方法是通过一个矩阵Tmn锐化采样协方差矩阵即
RCMT=RTCMT
(6)
式中：RCMT表示锐化后的协方差矩阵，TCMT第m行n列的元素为
Tmn==sin c((m-n)Δ/π)
(7)
式中：sin c(x)=sin(πx)/πx，Δ决定了零陷的宽度。

CMT方法能够自适应的在干扰方向附近形成具有一定宽度的零陷，但由于CMT 方法用一束虚拟干扰替代原有单个干扰，使得每个虚拟干扰的功率变小，导致零陷深度变浅，旁瓣上升。

1.3 LCSS方法
在零陷展宽的区域Θ内形成具有一定深度的零陷，即区域Θ的平均输出功率小于预设值η，可表示为
WHQW≤η
(8)
式中：W为权矢量，Q为导向矢量相关矩阵，即
(9)
对Θ区域平均输出功率抑制的最理想情况就是输出功率为0，即WHQW=0。

Q 可以被分解为Q=UΣUH，其中U是Q的特征向量矩阵，Σ是Q特征值按降序排列组成的对角阵。

假设Q的秩为r，Ur=[u1,u2,…,ur]是r个大特征值对应的特征向量组成的矩阵，可知当WHUr=0时，WHQW=0也能实现，即通过线性约束实现零陷展宽，最优权求解的问题转化为
(10)
式中：C=[a(θd),Ur]，e1=[1,0T]T。

再由拉格朗日乘子法可知，式(10)的解为
WLCSS(r)=-1C(CH-1C)-1e1
式中 r为Θ区域输出功率最小时的线性约束数。

1.4 PDNBB方法
PDNBB方法[13]是一种先进的零陷展宽方法，该方法先由式(9)构建导向矢量相关矩阵Q，对其特征分解，提取出较大特征值对应的特征向量矩阵Ur，将其作为投影子空间，构造投影矩阵T，再通过投影矩阵修正采样协方差矩阵，即
(12)
最后用RT进行波束形成。

PDNBB方法构建导向矢量相关矩阵Q，以及提取特征子空间Ur的过程与LCSS方法一致；但LCSS方法利用特征子空间构造线性约束实现零陷展宽，而PDNBB方法利用特征子空间进行投影变换实现零陷展宽。

由此可知，在理想情况下，PDNBB方法和LCSS方法的性能接近；而且两种方法都构造了导向矢量相关矩阵，因此性能都十分依赖于正确的阵列流型信息。

2 抗阵列流型误差的零陷展宽方法
2.1 NBACE方法
为了提高零陷展宽方法在阵列流型误差情况下的性能，本文在CMT和PDNBB方法的基础上，结合不确定集优化方法[15]，提出了NBACE方法。

首先在零陷展宽区域Θ内任意取一个点，得到理想情况下的干扰导向矢量利用不确定集优化方法对其进行校正，即
(13)
式中表示校正后的导向矢量表示约束参数。

为减小计算量，可由Lagrange乘子法得到
式中为Largrange乘子因子可通过求解约束方程得到，即
(15)
由此得到校正后的并利用构造导向矢量相关矩阵:
(16)
然后结合CMT方法，用协方差锐化矩阵对进行锐化，即
CMT=CTCMT
(17)
将CMT特征分解，得到
(18)
式中：λθk表示特征值，θk表示与特征值对应的特征向量。

在CMT中选取M个大特征值，提取出大特征值对应的特征向量，构成一个投影子空间为
(19)
再由投影子空间得到零陷展宽投影矩阵即
(20)
用零陷展宽投影矩阵修正采样协方差矩阵，得到投影变换后的协方差矩阵为
(21)
最后，用进行波束形成。

值得注意的是:式(16)对应于一个零陷展宽区域的情况，
增加零陷展宽区域需要增加取点数，每个零陷展宽区域只需取一个点即可；因此，NBACE方法与PDNBB方法以及LCSS方法相比，不需要在零陷展宽区域内取出
J(J≫N)个间隔相等的采样点构建导向矢量相关矩阵，运算量大幅度减小。

2.2 计算复杂度分析
NBACE方法中不确定集优化的计算复杂度为O(N3)，特征分解的计算复杂度为
O(N3)，协方差矩阵求逆的计算复杂度为O(N3)，因此NBACE方法的计算复杂度为O(N3)；CMT方法的计算复杂度也为O(N3)；LCSS方法和PDNBB方法的计
算复杂度相等，其中构造导向矢量相关矩阵的计算复杂度为O(JN2)，特征分解的计算复杂度为O(N3)，协方差矩阵求逆的计算复杂度为O(N3)，因此，LCSS方法和PDNBB方法的计算复杂度都是O(JN2)。

由此可知，NBACE方法的计算复杂
度与CMT方法相当，远低于LCSS方法和PDNBB方法。

3 仿真实验与性能分析
考虑远场情况下的非相干窄带信号，一个阵元数N为10的均匀直线阵，阵元间距为半波长，干扰噪声比(interferences to noise ratio, INR)为30 dB，噪声为高斯白噪声，期望信号来波方向为3°，干扰位于-30°、40°方向，零陷宽度设为10°，快拍数为100，约束参数为2。

下面通过四个例子将NBACE方法与CMT方法[6]、LCSS方法[12]以及PDNBB方法[13]进行比较，验证其零陷展宽效果和输出性能。

3.1 理想情况
考察理想情况下，NBACE方法的性能。

图1为四种方法在输入信噪比(signal to noise ratio, SNR)为0dB时的波束图；图2为四种方法SINR随输入SNR的变化曲线。

图1 理想情况下四种方法的波束图Fig.1 Beam patterns of the four methods in the case without array calibration error
由图1可知，在理想情况下，四种方法都能够展宽零陷，其中LCSS方法和
PDNBB方法的零陷展宽效果相近，都能得到很深的零陷；NBACE方法和CMT
方法的零陷展宽效果相近，且NBACE方法的零陷略深于CMT方法。

由图2可知，在理想情况下，PDNBB方法和LCSS方法的输出SINR都接近最优；当SNR低于0 dB时，NBACE方法和CMT方法的输出SINR接近最优，当SNR高于0 dB时，NBACE方法的输出SINR高于CMT方法。

由此可知，在理想情况下，四种零陷
展宽方法都能有效展宽零陷，其中PDNBB方法和LCSS方法性能最好，NBACE
方法性能优于CMT方法。

图2 理想情况下四种方法输出SINR随输入SNR变化曲线Fig.2 Output SINR versus input SNR of the four methods in the case without array calibration error
3.2 阵元位置误差情况
考察出现阵元位置误差时，NBACE方法的性能。

假设真实阵元位置与理想阵元位置的误差服从[-0.075 λ,0.075 λ]上的随机分布，其中λ为信号波长。

图3为四种
方法在SNR为0 dB时的波束图；图4为四种方法SINR随输入SNR的变化曲线。

图3 阵元位置误差情况下四种方法的波束图Fig.3 Beam patterns of the four methods in the case without array location error
由图3可知，存在阵元位置误差的情况下，LCSS方法和PDNBB方法都已经不能展宽零陷，因为LCSS方法和PDNBB方法都构造了导向矢量相关矩阵，所以性能十分依赖于正确的阵列流型信息；NBACE方法和CMT方法仍能展宽零陷，NBACE方法的零陷略深于CMT方法，且旁瓣低于CMT方法。

由图4可知，存
在阵元位置误差的情况下，PDNBB方法和LCSS方法的输出性能很差，NBACE
方法的输出性能最好。

由此可知，存在阵元位置误差的情况下，NBACE方法能够有效展宽零陷，且其零陷展宽效果和输出性能都优于其他方法。

图4 阵元位置误差情况下四种方法输出SINR随输入SNR变化曲线Fig.4 Output
SINR versus input SNR of the four methods in the case of array location error
3.3 阵元互耦情况
考察阵元互耦时，NBACE方法的性能。

假设仅相邻阵元间存在互耦，且相邻阵元间的互耦系数为0.75e-jπ/3。

图5为四种方法在在SNR为0 dB时的波束图；图6为四种方法SINR随输入SNR的变化曲线。

图5 阵元互耦情况下四种方法的波束图Fig.5 Beam patterns of the four methods in the case of mutual coupling
由图5可知，存在阵元互耦的情况下，PDNBB方法和LCSS方法已经不能展宽零陷；NBACE方法和CMT方法仍能展宽零陷，其中NBACE方法的零陷深于CMT 方法，且旁瓣低于CMT方法，NBACE方法的零陷展宽效果最好。

由图6可知，存在阵元互耦的情况下，PDNBB方法和LCSS方法的输出性能很差，CMT方法在SNR高于0 dB后输出性能下降，NBACE方法的输出性能始终最好。

由此可知，存在阵元互耦的情况下，NBACE方法能够展宽零陷，且其零陷展宽效果和输出性能都优于其他方法。

图6 阵元互耦情况下四种方法输出SINR随输入SNR变化曲线Fig.6 Output SINR versus input SNR of the four methods in the case of mutual coupling 3.4 运动干扰情况
考察存在阵元位置误差时，NBACE方法抗运动干扰的性能。

设定一个干扰信号从-40°向0°方向运动，运动速度为0.05°/快拍，其余仿真条件同3.2。

图7为四种方法在SNR为0 dB时的波束图；图8为四种方法输出SINR随输入SNR的变化曲线。

图7 运动干扰情况下四种方法的波束图Fig.7 Beam patterns of the four methods in the case of jamman motion
由图7可知，同时存在阵元位置误差和运动干扰的情况下，LCSS方法和PDNBB
方法都不能展宽零陷；NBACE方法和CMT方法能够展宽零陷，其中NBACE方
法的零陷深于CMT方法，旁瓣低于CMT方法，NBACE方法的零陷展宽效果最好。

由图8可知，同时存在阵元位置误差和运动干扰的情况下，PDNBB方法和LCSS方法的输出性能很差，NBACE方法的输出性能最好。

由此可知，同时存在
阵元位置误差和运动干扰的情况下，NBACE方法的零陷展宽效果和输出性能都优于其他方法，NBACE方法能够有效的抑制运动干扰。

图8 运动干扰情况下四种方法输出SINR随输入SNR变化曲线Fig.8 Output SINR versus input SNR of the four methods in the case of jammer motion
4 结论
1)不考虑阵列流型误差的情况下，NBACE方法的波束形成零陷展宽效果和输出性能优于CMT方法。

2)存在阵列流型误差和运动的干扰情况下， NBACE方法与现有的波束形成零陷展宽方法相比，具有更好的零陷展宽效果和输出性能，稳健性明显提高。

3)通过计算复杂度分析可知，NBACE方法的计算复杂度与CMT方法相当，远低
于LCSS和PDNBB方法。

本文考虑到天线阵在应用中很可能存在阵列流型误差这一实际问题，提出了NBACE方法，提高了波束形成零陷展宽方法在实际应用中的稳健性。

在以后的研究中，可以考虑导向矢量失配对NBACE方法的影响，进一步提高该方法的稳定性。

【相关文献】
[1] HUANG Lei, ZHANG Jing, XU Xu, et al. Robust adaptive beamforming with a novel interference-plus-noise covariance matrix reconstruction method[J]. IEEE transactions on
signal processing, 2015, 63(7): 1643-1650.
[2] LI Wenxing, MAO Xiaojun, YU Wenhua, et al. Robust adaptive array beamforming based on modified norm constraint algorithm[J]. Applied computational electromagnetic society, 2014, 29(12): 1060-1066.
[3] LI Hongtao, WANG Ke, WANG Chaoyu, et al. Robust adaptive beamforming based on worst-case and norm constraint[J]. International journal of antennas and propagation, 2015, 2015: 1-7.
[4] HOSSAIN M S, MILFORD G N, REED M C, et al. Robust efficient broadband antenna array pattern synthesis techniques[J]. IEEE transactions on antennas and propagation, 2014, 62(9): 4537-4546.
[5] GERSHMAN A B, NICKEL, BOHME J F. Adaptive beamforming algorithms with robustness against jammer motion[J]. IEEE transactions on signal processing, 1997, 45(7): 1878-1885.
[6] MAILLOUX R J. Covariance matrix augmentation to produce adaptive pattern troughs[J]. Electronics letters, 1995, 31(25): 771-772.
[7] ZATMAN M. Production of adaptive array troughs by dispersion synthesis[J]. Electronics letters, 1995, 31(25): 2141-2142.
[8] GUERCI J G. Theory and application of covariance matrix tapers for robust adaptive beamforming[J]. IEEE transactions on signal processing, 1999, 47(4): 977-985.
[9] 武思军, 张锦中, 张署. 阵列波束形成的零陷加宽算法研究[J]. 哈尔滨工程大学学报, 2004, 10(5): 658-661.
WU Sijun, ZHANG Jinzhong, ZHANG Shu. Research on beamforming of wide nulling algorithm[J]. Journal of Harbin Engineering University, 2004, 10(5): 658-661.
[10] ER M H. Technique for antenna array pattern synthesis with controlled broad nulls[J]. IEE proceedings H, 1988, 135(6): 375-380.
[11] GODARA L C, SAYYAH JAHROMI M R. An optimized sector nulling technique for broadband antenna array[C]// International Conference on Signal Processing and Communication System. 2008: 1-5.
[12] AMAR A, DORON M A. A linearly constrained minimum variance beamformer with a prespecified suppression level over a pre-defined broad null sector[J]. Signal processing, 2015, 109: 165-171.
[13] 李文兴, 毛晓军, 孙亚秀. 一种新的波束形成零陷展宽方法[J]. 电子与信息学报, 2014, 36(12): 2882-2888.
LI Wenxing, MAO Xiaojun, SUN Yaxiu. A new algorithm for null broadening[J]. Journal of electronics and information technology, 2014, 36(12): 2882-2888.
[14] MAO Xiaojun, LI Wenxing, LI Yingsong, et al. Robust adaptive beamforming against signal steering vector mismatch and jammer motion[J]. International journal of antennas
and propagations, 2015, 2015: 1-12.
[15] LI J, STOICA P, WANG Z. On robust capon beamforming and diagonal loading[J]. IEEE transactions on signal processing, 2003, 51(7): 1702-1714.。