一种基于导向矢量约束的恒模盲波束形成算法

刘　可， 钱华明， 马俊达

(哈尔滨工程大学 自动化学院， 哈尔滨 150001)



一种基于导向矢量约束的恒模盲波束形成算法

刘可， 钱华明， 马俊达

(哈尔滨工程大学 自动化学院， 哈尔滨 150001)

针对阵列信号波达角(direction of arrival，DOA)先验信息已知的情况，利用信号的恒模特性，在卡尔曼滤波(Kalman filter，KF)结构下，提出一种附加阵列导向矢量约束的自适应波束形成算法.对约束情况下的卡尔曼滤波目标函数运用拉格朗日乘子法，求得约束条件下的最优估计表达式，并将其推广到无迹卡尔曼滤波(unscented Kalman filter，UKF)算法中，通过约束迭代算法对阵列估计信号的导向角施加约束，实现约束UKF自适应波束形成算法的最优权值分配.仿真过程中，用所提算法与约束恒模迭代最小二乘算法和约束最小方差迭代最小二乘算法作对比，表明表明，该算法在收敛速度、信噪比、稳健性、跟踪性能方面具有较好的性能.

信号处理；自适应滤波；无迹卡尔曼滤波；波束形成；约束优化技术

在现代无线通讯系统中，通讯信号在干扰环境下的检测和提取是当今研究的热点问题.自适应波束形成技术，为解决这类问题提供一种有效的技术手段，它能在有用信号到达方向形成高增益，在干扰方向形成零陷，从而对干扰信号形成有效抑制，而被广泛应用于导航、雷达、声纳、医学图像等领域.盲自适应波束形成技术能在缺少很多先验信息的情况下形成有用信号的阵列响应.它在抑制干扰、提高追踪性能、降低计算复杂度方面具有重要的研究价值.附加约束条件的自适应波束形成算法按设计准则，主要分为两类：约束最小方差(constrained minimum variance, CMV)准则和约束恒模(constrained constant modulus, CCM)准则.CMV准则是对有用信号的波达角施加约束的情况下使得阵列输出方差最小.CCM准则是使得波束形成的模值输出偏差最小，且对有用信号的阵列输出施加条件约束.随机梯度法 (stochastic gradient, SG)是自适应波束形成的一种经典方法，它通过选取合适的步长，实现对最小均方目标函数的渐近逼近，不需要估计输入信号的协方差矩阵，能有效减少计算量.文献[1]对Bussgan类型的恒模算法给出了步长边界的计算公式；文献[2]利用CMV-SG算法，对CDMA通信系统进行多用户检测，并取得了较好的效果；文献[3]针对不同结构，提出了3种CCM-SG递推算法，并应用于CDMA系统中.该方法的不足是收敛速度慢，而且步长的选取相对困难，步长过小会导致较慢的收敛速度，步长过大又会造成偏差较大，甚至系统不稳定.

在众多自适应波束形成算法中，迭代最小二乘算法[4](recursive least squares, RLS)具有收敛快速性和算法实用性的特点.但该算法存在遗忘因子难以预先确定的问题.文献[5-6]针对这个问题提出可变遗忘因子的RLS算法；文献[7]将CCM-RLS算法应用于智能天线的波束形成中；文献[8]利用CMV-RLS和CCM-RLS算法消除CDMA系统的多径干扰.文献[9]利用UKF算法提出一种自适应盲波束形成算法，但未考虑附加约束条件.本文针对波达角先验信息已知的情况，利用信号的恒模特性，提出一种导向矢量约束的UKF自适应盲波束形成算法.

1　系统模型

1.1阵列模型

定义符号：(·)T为转置运算，(·)H为共轭转置运算，(·)*为共轭运算，In×n为n×n的单位阵，0m×n为m×n的零矩阵，R为实数域，C为复数域.

假设信号为点源、远场、传播介质均匀，且传播无衰减.L元传感器阵列如图1所示，信号的波达角为θi,i=0,1,…,N-1，N为信号数，阵列第k次快拍的接收向量xk∈CL×1为

信号的导向矢量矩阵为

波束形成的输出为

yk=wHxk.

式中w=[w(1),…,w(L)]T∈CL×1是复权值向量.

1.2状态空间模型

假设存在一个未知系统wk，输入信号xk的阵列响应具有恒模特性，模值为r，则状态空间模型为

(1)

(2)

式中p为一个常数，通常设置为2.

将式(1)、(2)变形可得

(3)

(4)

利用式(2)、(4)推得

由式(1)～(4)，得到状态空间模型的表达式为

2　约束自适应波束形成算法

约束自适应波束形成算法是通过最小化代价函数，对变量施加条件约束，以获取最优权矢量.写成表达式的形式为

式中：Jw为代价函数；B∈CL×q为约束矩阵；f∈Cq×1为约束向量.

2.1CMV算法

常用的CMV算法主要包括CMV-SG算法和CMV-RLS算法.CMV-SG算法的代价函数为

运用拉格朗日乘子法可求得递推公式

式中μ为算法步长.

CMV-RLS算法的代价函数为

式中α为RLS算法的遗忘因子.

运用拉格朗日乘子法可得

通过迭代运算，实现协方差矩阵Rk的求逆运算，由

运用矩阵求逆引理[10]得

若令

则CMV-RLS递推求解公式为

2.2CCM算法

恒模算法的代价函数为

式中：r为信号的恒模值，通常设置为1; p和q为非负常数，通常设置为2; yk表示阵列输出向量y的第k次快拍;E[·]表示期望.

与CMV算法类似，常用的CCM算法包括CCM-SG算法和CCM-RLS算法.CCM-SG算法的代价函数为

采用与CMV-SG相似的运算，可得

[xk-aH(θ0)xka(θ0)].

CCM-RLS的代价函数为

运用拉格朗日乘子法，可得

上式协方差矩阵Rk的求逆，可通过迭代实现，由于

由矩阵求逆引理[10]可得

式中

若令

则CCM-RLS递推求解公式为

更为详细的推算过程可参考文献[7-8].

3　恒模盲波束形成算法

3.1KF代价函数的最优解

因为直接利用1.2节建立的状态空间模型，构造约束UKF的代价函数求解比较困难，本文采用的思路是先对约束条件下的KF构造代价函数，求出最优估计的表达式，再将其推广到约束UKF波束形成算法中.

有关KF的详细递推公式可参考文献[11-13].利用拉格朗日乘子法，可求得约束状态估计及其协方差矩阵为

(5)

(6)

其中：

3.2UKF约束迭代算法

UKF是KF的非线性扩展，是通过UT变换，实现非线性系统线性化的一种近似.该迭代算法是约束KF向非线性约束状态估计的延伸，约束迭代过程如下.

1)计算sigma点：

(7)

2)传播sigma点：

(8)

式中ιj(·)表示(·)的第j列元素.

3)估计约束向量：

(9)

4)计算协方差矩阵：

(10)

(11)

3.3恒模盲波束形成算法的步骤

在自适应波束形成算法中，假设阵列波达角θ0为已知的先验信息，则约束条件为

即约束矩阵D=aH(θ0)，约束向量d=1.

附加导向矢量约束的UKF波束形成算法主要分为5个环节，第1个环节是迭代参数的初始化；第2个环节为系统的状态预测(步骤1～4)；第3个环节为系统的状态更新(步骤5～11)；第4个环节是约束迭代过程(步骤12)；第5个环节是最优权向量的输出(步骤13～14).详细的求解步骤如下.

步骤1选取sigma点：

其中Pk-1|k-1为先验协方差，χk-1|k-1为选取的sigma点.

步骤2传播sigma点：

步骤3估计系统状态：

步骤4估计误差协方差：

步骤5计算系统的输出sigma点：

式中函数g(·)表示对ιj(χk|k-1)中的每个元素的绝对值的p次方计算.

步骤6估计系统输出：

步骤7估计自协方差矩阵：

步骤8估计互协方差矩阵：

步骤9计算增益矩阵：

步骤10计算状态估计：

步骤11计算误差协方差矩阵：

步骤12利用式(5) ～(11)进行约束迭代.

步骤13计算系统的最优权值：

该算法的计算流程如图2所示.

图2　算法流程图

4　仿真验证

[7]的对比方式，以信干噪比(signaltointerferenceplusnoiseratio,SINR)为指标进行对比仿真，因为文献[2，7-8]已经验证CCM-RLS和CMV-RLS算法要好于CCM-SG和CMV-SG算法，所以本文仅用所提算法跟文献[7-8]的CCM-RLS，CMV-RLS作对比仿真，来验证算法的性能.

输出SINR的公式为

实验1给出了SINR随快拍数的变化，由图3可以看出，随着迭代次数的增加，所提算法最终收敛的SINR要高于CMV-RLS和CCM-RLS，而且当迭代次数到达200次时，SINR就可以完全收敛，其他两种算法需要迭代500次以上，由此可见，在收敛速度和收敛的SINR数值上都好于另外两种算法.该仿真也同时验证了CCM-RLS算法性能要好于CMV-RLS，该结论与文献[7-8]一致.

图3　实验1 SINR相对于迭代次数的变化

实验2为假定波达角DOA为11°，即存在1°的导向角失配情况，从图4中可以看出，3种算法存在不同程度的性能下降，但CMV-RLS和CCM-RLS算法性能下降更为明显，而所提算法仍能保持较好的性能，因而所提算法相比于其他两种算法具有较好的稳健性.

图4　实验2 SINR相对于迭代次数的变化

实验3为阵列维数较高时的信干噪比变化情况，此时设定L=60，从图5中可以看出，CCM-RLS算法收敛SINR为15 dB, CMV-RLS算法的收敛SINR为13 dB，相比于阵列维数较低时，存在着SINR下降的情况，而所提算法在收敛速度和SINR上仍保持较好的性能.

实验4为非平稳条件下的仿真，为了让仿真更加明显，此次仿真将迭代次数增加到了4 000次，并在迭代次数为2 000次时，增加两个INR为3 dB的干扰信号.从图6可以看出，3种算法在2 000次迭代处都发生了SINR的下降，但相比于其他两种算法，所提算法具有较快的恢复速度，体现出该算法具有较好的跟踪性能和稳态恢复能力.

图5　实验3 SINR相对于迭代次数的变化

图6　实验4 SINR相对于迭代次数的变化

5　结　论

1)提出一种新型的UKF波束形成算法.先对恒模算法的代价函数进行变形，然后在附加导向矢量约束的条件下，利用UKF的迭代更新等式，求解波束形成器的最优权矢量.

2)通过数值仿真实验可以得出，相比于经典的CCM-RLS和CMV-RLS算法，该算法具有如下优势：具有较高的输出SINR和较快的收敛速度；在导向角失配情况下表现出较好的稳健性；在阵列维数较高时，仍保持较好的性能指标；在非平稳情况下体现出较好的稳态恢复能力.

参考文献

[1] ABRAR S, NANDI A K. An adaptive constant modulus blind equalization algorithm and its stochastic stability analysis [J]. IEEE Signal Processing Letters, 2010, 17(1): 55-58. DOI: 10.1109/LSP.2009.2031765.

[2] XU Zhengyuan, TSATSANIS M K. Blind adaptive algorithms for minimum variance CDMA receivers [J]. IEEE Transactions on Communications, 2001, 49(1): 180-194. DOI: 10.1109/26.898261.

[3] MIGUEZ J, CASTEDO L. A linearly constrained constant modulus approach to blind adaptive multiuser interference suppression [J]. IEEE Communications Letters, 1998, 2(8): 217-219. DOI: 10.1109/4234.709436.

[4] CHEN Yuxin, LE NGOC T, CHAMPAGNE B, et al. Recursive least squares constant modulus algorithm for blind adaptive array [J]. IEEE Transactions on Signal Processing, 2004, 52(5):1452-1456. DOI: 10.1109/TSP.2004.826167.

[5] QIN Boya, CAI Yunlong, CHAMPAGNE B, et al. Low-complexity variable forgetting factor constant modulus RLS-based algorithm for blind adaptive beamforming [J]. Signal Processing, 2014, 105(1):277-282. DOI: 10.1016/ j.sigpro.2014.06.003.

[6] BHOTTO M Z A, ANTONIOU A. New improved recursive least-squares adaptive-filtering algorithms [J]. IEEE Transactions on Circuits and Systems I Regular Papers, 2013, 60(6): 1548-1558. DOI: 10.1109/TCSI.2012.2220452.

[7] WANG Lei, DE LAMARE R C. Blind adaptive beamforming based on constrained constant modulus RLS algorithm for smart antennas[EB/OL]. (2013-03-07)[2015-02-01]. http://arxiv.org/pdf/1303.1848v1/pdf.

[8] DE LAMARE R C, SAMPAIO NETO R. Blind adaptive MIMO receivers for space-time block-coded DS-CDMA systems in multipath channels using the constant modulus criterion [J]. IEEE Transactions on Communications, 2010, 58(1):21-27. DOI: 10.1109/TCOMM.2010.01.070549.

[9] BHOTTO M Z A, BAJIC I V. Constant modulus blind adaptive beamforming based on unscented Kalman filtering [J]. IEEE Signal Processing Letters, 2015, 22(4): 474-478. DOI: 10.1109/LSP.2014.2362932.

[10]西蒙.赫金. 自适应滤波器原理 [M].4版. 北京：电子工业出版社，2006：344-350.

HAYKIN S. Adaptive filter theory [M].4th ed. Beijing: Electronic Industry Press, 2006: 344-350.

[11]付梦印，邓志红，张继伟. Kalman滤波理论及其在导航系统中的应用[M]. 北京：科学出版社，2003:16-36.

FU Mengyin，DENG Zhihong，ZHANG Jiwei. Kalman filter theory with applications in navigation system [M]. Beijing: Science Press, 2003:16-36.

[12]JULIER S J, UHLMANN J K. Unscented filtering and nonlinear estimation [J]. Proceedings of the IEEE, 2004, 92(3): 401-422. DOI: 10.1109/JPROC.2003.823141.

[13]TEIXEIRA B O S, CHANDRASEKAR J, TORRES L A B, et al. State estimation for linear and non-linear equality-constrained systems [J]. International Journal of Control, 2009, 82(5): 918-936. DOI: 10.1080/002071708023 70033.

(编辑魏希柱)

A constant modulus blind beamforming algorithm based on constrained direction vectors

LIU Ke，QIAN Huaming，MA Junda

(College of Automation, Harbin Engineering University, Harbin 150001, China)

Aiming at the case of the known knowledge of the direction of arrival of the desired signal, by using the constant modulus feature, a new blind adaptive beamforming algorithm is proposed under the frame of the Kalman filter. According to the Lagrange multiplier method, the optimal estimation expression of the system states can be derived by minimizing the constrained cost function. Then, the optimal weight vector of the adaptive beamformer can be obtained by using the iteration and update equations of the unscented Kalman filter. In the simulation, the proposed algorithm is compared with the constrained constant modulus recursive least square (CCM-RLS) and constrained minimum variance recursive least square (CMV-RLS) to demonstrate its effectiveness in the terms of the convergence speed, signal to interference plus noise ratio, robustness to changeable environments and tracking capability in the non-stationary condition.

signal processing; adaptive filtering; unscented Kalman filter; beamforming; constrained optimization technique

10.11918/j.issn.0367-6234.2016.09.026

2015-03-31

国家自然科学基金(61573113)；哈尔滨市科技创新人才研究专项资金(优秀学科带头人)资助项目(2014RFXXJ074)作者简介: 刘可(1986—)，男，博士研究生；

钱华明，qianhuam@sina.com

TN911

A

0367-6234(2016)09-0151-06

钱华明(1965—)，男，教授，博士生导师