基于小波分析的改进旁瓣相消宽带信号波束形成算法

刘 影，谢 驰

(1. 电子科技大学能源科学与工程学院 成都 611731; 2. 四川大学电子信息学院 成都 610065)

目前自适应阵列信号处理技术在雷达，声纳和移动通信系统中有着越来越多的重要应用[1]。宽带信号波束形成算法可以采用时域和频域分析方法。频域方法是对阵列接收数据按照频率进行分块，并通过DFT变换将时域数据转换成频域数据，由于在子块中信号频宽较窄，可以通过独立的窄带波束形成器实现波束成形，最后在通过合成实现宽带信号波束成形。频域方法复杂度较小，但不能得到连续的时域波束输出[2]。文献[3]和文献[4]分别采用频率矩和滑动窗函数的方法，在有效地进行宽带波束形成的同时也抑制了宽带干扰。文献[5]将分数阶傅里叶变换应用于宽带信号频域多波束形成中，从而提高了运算的速度。文献[6]采用基于均匀一致性旁瓣级的设计方法，将宽带信号用一组FIR(finite impulse response)滤波器拟合，从而实现恒定束宽的波束形成，但该方法存在着迭代收敛的步长难以确定的问题，因而算法的精度不高，收敛速度慢。文献[7-9]基于高斯模型的超宽带信号，提出数字延迟线和分数时延滤波器相结合的方法，但由于在时域处理方法中往往需要较多的延时单元，并且涉及到高阶协方差矩阵的求逆问题，因此对于宽带或超宽带信号处理，该方法采样频率高，数据量大，难于在工程中实现。在宽带信号处理中，有限脉冲响应(FIR)抽头系数较多，为了减轻矩阵运算的高复杂性可以采用无限脉冲响应(IIR)，然而IIR算法的不稳定性也使得收敛速度变慢[10]。

对于阵列信号处理技术中波束形成目前最广泛采用的是最小方差无失真响应(minimum variance distortionless response, MVDR)波束合成、对角加载方法、旁瓣对消器(generalized sidelobe canceller,GSC)[11-14]。MVDR是通过调节阵列的输出功率达到最佳的阵列权重向量。对角加载算法是通过对投影变换技术对阵列接收数据协方差矩阵进行预处理，在约束区域形成稳定的响应幅度。GSC是一个自适应波束形成器，通过求解线性约束得最小方差合成波束实现。最小方差的求解可以利用众所周知的自适应滤波算法，如最小均方(least mean square,LMS)、递归最小二乘(recursive least square, RLS)进行无约束的自适应优化[15-16]。然而，这些优化方法都涉及对矩阵求逆的计算，对于宽带接受信号，矩阵计算量非常大。因此，本文提出了一种基于小波滤波器组的改进GSC自适应波束形成方法，该方法在降低计算低复杂度和提高收敛速度方面具有一定优势。

1 旁瓣相消波束形成原理

图1为GSC的结构示意图。GSC自适应的原理是：将接收的已知期望信号方向信息的阵列信号x(n)变换为上下两个支路信号。上支路信号通过归一化期望信号的导向矢量h0变换为参考信号d0(n)，d0(n)同时含期望信号和干扰信号；下支路信号通过正交阻塞矩阵B0阻塞掉期望信号，这样x0(n)中就只含干扰，再将变换后的上下支路信号进行维纳滤波wx0，从而自适应的抵消干扰。

图1 旁瓣相消器结构

GSC波束形成结构中阵列相关输出量为：

式中，h0为归一的期望信号矢量，这里取归一化为：

则T为非奇异阵，得：

由GSC波束形成中的自适应对消干扰原理，得参考信号表达式为：

由式(8)可知：GSC波束形成结构中上下支路的变换是把阵列数据向量向两个正交的子空间投影，h0将x(n)投影到期望信号对应的一维子空间中，而将x(n)投影到与期望信号子空间正交的N维子空间，输出信号x(n)经正交投影变换后进行维纳滤波器实现干扰对消。

2 改进的GSC波束形成算法

GSC波束结构中存在构造满秩阻塞矩阵B0的问题，在实际系统中由于各种误差的存在，如阵元位置误差、阵元幅相误差、阵元之间的互耦和通道频率特性失配等，使得构造的B0矩阵难以满足条件，从而存在导向矢量失配的问题，导致阻塞矩阵B0不能完全阻塞掉期望信号，造成部分有用信号互相抵消，输出信干噪比(signal interference noise ration,SINR)下降，波束图中的旁瓣电平被升高的现象。

改进的GSC波束形成算法采用小波滤波器组在频域中阻塞接收信号的有用分量，然后对输出信号进行重构并求出协方差矩阵，再用重构的协方差矩阵形成自适应波束的权向量，其结构图如图2所示。

图2 改进GSC波束形成算法的结构图

改进的GSC波束形成算法结构中，B0为行滿秩正交阻塞矩阵，其作用是阻塞接收信号中的期望信号，但在实际构造B0阵列中如果存在误差，则将有一部分期望信号泄露到辅助支路中，从而引起期望信号相消的现象。对于阵列接收信号x(t)的协方差矩阵为：

对式(9)进行矩阵求逆，得：

GSC波束形成器的权矢量为：

通过式(11)可知：

故式(12)可以改写为：

小波滤波器组能够有效地对期望信号进行检测，因此小波滤波器对阵列信号接收的数据进行频域滤波后，可以有效阻塞期望信号。相应的输出干扰和噪声的协方差矩阵Ri+n特征值分解为：

式中，ui+n为干扰和噪声的值；Λi=diag{λ1,λ2,…,λP-1}包含P-1个干扰信号的特征值；Λn=diag{λP,λP+1,…,λM}为噪声的特征值；Ei为干扰子空间；En为噪声子空间。对Ei干扰子空间重构得到干扰空间的协方差矩阵为：

当对阵列天线噪声估计不足时会造成协方差矩阵特征值分散，使得自适应波束的旁瓣升高。可以采用用噪声的平均功率进行对角加载，加载量的运用对协方差矩阵中的大特征值不产生影响，但使噪声对应的小特征值加大并近似等于加载量，从而抑制或降低了小特征值的扰动，得到较好的旁瓣压低性能。有：

转换到阵元空间为：

式中，G为小波综合滤波器的矩阵形式。根据式(13)在期望信号方向矢量和协方差矩阵准确已知的条件下，干扰噪声的协方差矩阵可以代替阵列的协方差矩阵，权矢量保持不变，对干扰子空间重构后的自适应波束的权向量表示为：

由于阵列信号经小波滤波器组后协方差矩阵中不包含期望信号分量，从而克服了存在系统误差时自适应波束形成中期望信号相消的问题，解决了在实际应用中GSC波束形成算法需要构造与期望信号的导向矢量完全正交的阻塞矩阵的难点问题。

3 仿真分析

1)考察存在输入信号方位估计误差和输入信噪比不同的情况时，改进算法的稳健性能。

仿真1：设实际的期望信号方向为13°，估计的期望信号方向为18°，存在Δ=5°输入信号方位估计误差，两个干扰信号分别来自-30°和60°，信噪比为20 dB。图3给出了输出信干噪比随输入信噪比变换的曲线。

图3中，输入信噪比高于5 dB时采用对角加载效果较好，而改进的算法对阵列信号的方位估计误差变化敏感程度低于基于MVDR方法的对角加载波束形成方法，并且随着输入信干噪比的提高，改进算法由于不需要构造正交阻塞矩阵滤除干扰和噪声，因此在阵列信号的方位估计存在误差时，其能够通过小波滤波器组更好地滤除噪声和干扰。通过式(18)可以看出改进算法所构造的协方差矩阵优于对角加载方法，因而其鲁棒性优于MVDR算法和对角加载算法。

图3 存在输入信号方位估计误差时性能比较

2)考察存在信号输入方位估计误差和阵元位置误差时算法的稳健性。

仿真2：阵元的位置误差为高斯分布，其方差为ε2，信噪比SNR=5 dB，其他参数与仿真1相同。图4和图5分别给出了阵列输出信干噪比随输入信号方位估计误差的变化曲线和阵列输出信干噪比随阵元位置误差变化的性能曲线。图6为系统在阵元位置误差为0.2λ时的波束图。

图4 输入信号方位估计误差变化性能比较

从图4～图6可以看出，改进的算法对阵元方位估计指向误差和阵元位置误差同样有较好的稳健性和波束保形性，其性能优于传统的MVDR方法和基于MVDR的对角加载方法。对于MVDR波束形成在较低的信干噪比情况下，由于协方差矩阵的特征值会发生抖动，因此在具有较大的阵元方位估计误差时，MVDR波束形成算法会将期望信号当成干扰信号，产生期望信号之间对消现象，使得输出信干噪比的性能损失随着信噪比的减小而变得越发严重。基于MVDR方法的对角加载波束形成在采用对角加载后使得协方差矩阵小特征值散布变小，从而改进算法相比于MVDR方法能够得到稳定的方向图和较低的副瓣，其性能优于MVDR方法。

图5 阵元位置误差变化性能比较

图6 阵元位置误差为0.2λ时的波束图

4 结 束 语

本文提出了基于小波分析的改进GSC波束形成算法，通过对GSC波束形成前构造小波滤波器实现对接收阵列信号的预处理，将信号空间与噪声空间分离，并滤出噪声信号。经算例研究、系统研制和仿真运行，结论如下：

1)解决由于阵列信号经小波滤波器组后其协方差矩阵中不包含期望信号分量，从而克服了存在系统误差时自适应波束形成中期望信号相消的问题。

2)解决在实际应用中GSC波束形成算法中需要构造与期望信号的导向矢量完全正交的阻塞矩阵的难点问题，提高GSC波束形成算法的收敛速度。

3)仿真研究在阵元存在误差下，对比了改进算法、MVDR算法、对角加载算法，验证了本文所提出算法的有效性。

4)该算法依赖于小波算法对期望信号分析，因此小波母函数的选取和信号分解层次决定了在原始信号中能否有效重构期望信号，增加了算法全局的复杂度。

[1]BIN Cao, JIA Yu, YE Wang, et al. Enabling polarisation filtering in wireless communications： models, algorithms and characteristics[J]. IET Communications, 2013, 7(3)：247-254.

[2]李根, 马彦恒, 董建. 宽带阵列雷达发射波束形成方法[J].电光与控制, 2016, 23(6)： 100-104.LI Gen, MA Yan-heng, DONG Jian. Transmitting beam forming method for wideband array radar[J]. Elecronics Optics & Control, 2016, 23(6)： 100-104.

[3]陈希信, 韩彦明, 于景兰. 高频雷达自适应波束形成抗干扰研究[J]. 电波科学学报, 2010, 25(6)： 1169-1174.CHEN Xi-xin, HAN Yan-ming, YU Jing-lan. Research for interference suppression in HF radar based on ADBF[J].Chinese Journal of Radio Science, 2010, 25(6)： 1169-1174.

[4]LI Y, JANDHYALA V. Design of retrodirective antenna arrays for short-range wireless power transmission[J]. IEEE Trans on Antennas and Propagation, 2012, 60(1)： 206-211.

[5]陈小龙, 刘宁波, 王国庆. 基于Radon-分数阶傅里叶变换的雷达动目标检测方法[J]. 电子学报, 2014, 42(6)：1074-1080.CHEN Xiao-long, LIU Ning-bo, WANG Guo-qing. Radar detection method for moving target based on radon-fourier fractional fourier transform[J]. Acta Electronica sinica, 2014,42(6)： 1074-1080.

[6]范展, 梁国龙. 基于凸优化的最小旁瓣恒定束宽时域宽带波束形成[J]. 电子学报, 2013, 41(5)： 943-948.FAN Zhang, LIANG Guo-long. Broadband beamforming with minimum sidelobe constant beamwidth based on convex optimization[J]. Acta Electonica Sinica, 2013, 41(5)：943-948.

[7]ROHDE G K, BUCHOLTZ F, NICHOLS J M. Maximum empirical likelihood estimation of time delay in independently and identically distributed noise[J]. IET Signal Processing, 2014, 8(7)： 720-728.

[8]周鲁苹, 王福才. 宽带主动目标信号的频率时延联合参量估计[J]. 控制工程, 2015, 22(6)： 1053-1056.ZHOU Lu-ping, WANG Fu-cai. Joint parameters estimation of frequency and time delay for broadband active target signal[J]. Control Engineering of China, 2015, 22(6)：1053-1056.

[9]李萌, 黄忠华. 超宽带引信取样脉冲宽度与相关接收输出信号幅度关系研究[J]. 兵工学报, 2016, 37(11)：1989-1994.LI Meng, HUANG Zhong-hua. Research on relation between sampling pulse width and output signal amplitude for ultra-wideband radio fuze[J]. 2016, 37(11)： 1894-1994.

[10]ACOSTA R, HARIHARAN S I. 2-D-IIR time-delay-sum linear aperture arrays[J]. IEEE Antennas and Wireless Propagation Letters, 2014, 13： 591-594.

[11]YAN She-feng, SUN Hao-hai, SVENSSON U P, et al.Optimal modal beamforming for spherical microphone arrays[J]. IEEE Trans on Audio, Speech, and Language Processing, 2011, 19(2)： 361-371.

[12]李文兴, 毛晓军, 翟助群. 抗导向矢量失配的零陷展宽波束形成算法[J]. 哈尔滨工业大学学报, 2016, 48(11)：116-122.LI Wen-xing, MAO Xiao-jun, ZHAI Zhu-qun. Null broadening beamforming against steering vector mismatch[J]. Journal of Harbin Institute of Technology,2016, 48(11)： 116-122.

[13]SIM S L, ER M H. An effective quiescent pattern control strategy for GSC structure[J]. IEEE Signal Processing Letters, 1996, 3(8)： 236-238.

[14]孔永飞, 吴海洲. 一种基于广义旁瓣相消的稳健降维方法[J]. 无线电工程, 2016, 46(9)： 41-44.KONG Yong-fei, WU Hai-zhou. A robust dimension reduction method based on GSC[J]. Radio Engineering,2016, 46(9)： 41-44.

[15]刘成城, 刘亚奇, 赵拥军, 等. 基于广义旁瓣对消器的Laguerre宽带波束形成[J]. 雷达学报, 2015, 4(3)：295-300.LIU Cheng-cheng, LIU Ya-qi, ZHAO Yong-jun, et al.Laguerre wideband beamforming algorithm based on generalized sidelobe canceller[J]. Journal of Radars, 2015,4(3)： 295-300.

[16]ERNER S, APOLINARIO J A, de CAMPOS M L R. On the equivalence of RLS implementations of LCMV and GSC processors[J]. IEEE Signal Processing Letters, 2003,10(12)： 356-359.