2012年第31卷第3期 传感器与微系统(Transducer and Microsystem Technologies) 15 t. t声 t 研究与探讨(i 、 一种新的LMS自适应滤波算法分析仿真研究 贺洪江,王春霞 (河北工程大学信息与电气工程学院,河北邯郸056038) 摘要:传统变步长最小均方(LMS)算法存在收敛速度慢、易受噪声干扰等缺点,为了提高算法的性能, 通过对变步长LMS算法进行分析研究,在步长因子 (n)与误差信号e(n)的相关统计量之间建立一种新 的非线性函数关系,提出了一种新的变步长LMS自适应滤波算法。该算法采用误差信号的自相关时间均 值来调节步长,并用绝对估计误差的扰动量以加快自适应滤波器抽头权向量的收敛。理论分析与计算机 仿真结果表明:与SVSLMS和G—SVSLMS算法比较,该算法具有较快的收敛速度、较小的稳态误差以及较 强的抗干扰能力。 关键词:变步长;最小均方算法;收敛速度;稳态误差;噪声抑制能力 中图分类号:TP301.6 文献标识码:A 文章编号:1000-9787(2012)03-0015-03 Analysis and Anal3"and simulation study on a new LMSstudy ue se’l。f 、 adaptive filtering algorithm HE Hong-jiang.WANG Chun—xi (School of Information and Electrical Engineering,Hebei University of Engineering,Handan 056038,China) Abstract:The common variable step size LMS algorithm has many weakness,such as poor convergence speed and sensitivity to noise.A uew variable step size adaptive filter algorithm is presented.A new non—linear function between step factor and error signal is established.In this algorithm,the step size factor is adjusted by the absolute value of the product of the current and former errors.The algorithm also introduces the disturbance of the absolute estimation error to update the tapping vector of the self-adaptive filter.The theoretical analysis and simulation results shows that compared with SVSLMS and G—SVSLMS algorithm,the new algorithm has faster convergence speed,lower steady state elror and better performance of noise suppression. Key words:variable step size;LMS algorithm;convergence speed;steady—state error;performance of noise SH— ppression 0引 言 度和稳态误差2个性能指标之间权衡。为此,人们提出了 最小均方(LMS)算法是由Widrow和Hoff两人在1960 多种变步长LMS自适应滤波算法 。 年提出的,由于其计算复杂度低、易于实现等优点 ,广泛 文献[4]中提出的SVSLMS算法的收敛速度较快,但 应用于自适应控制、系统辨识、信号处理和噪声抵消等领 在自适应稳态阶段的步长变化较大,稳态误差也较大。文 域,同时自适应滤波器在通信、雷达、工业控制、地震预报及 献[5]中提出的G—SVSLMS算法是在SVSLMS的基础上进 生物医学电子学等领域也已经有了越来越广泛的应用。初 行了改进,稳态步长变化比较缓慢,收敛速度也较快。文 始收敛速度、对时变系统跟踪能力、稳态误差以及抗噪声干 献[6]提出一种新的变步长函数,具有在初始阶段和未知 扰能力是衡量LMS算法优劣的重要性能指标 。传统 系统时变时步长自动增大而在稳态时步长很小的特点,但 LMS算法中固定步长的取值不同会影响算法的性能。减 收敛速度需要进一步提高。SVSLMS和G—SVSLMS算法的 小步长取值可以降低稳态误差,但会降低算法的收敛速度 步长受误差调节,抗噪声能力比较低。 和对时变系统的跟踪能力;增大步长取值可以提高收敛速 本文算法采用输入信号与误差信号不相关的特点,用 度,但会增大稳态误差。传统的固定步长的LMS算法的大 误差信号的自相关时间均值来调节步长,并用绝对估计误 收敛速度和小稳态误差不能同时满足,这就要求在收敛速 差的扰动量以更新自适应滤波器的抽头向量,算法性能较 收稿日期:2011-07-24 16 传感器与微系统 第3l卷 上述算法有较大的优越性。 1 自适应滤波算法原理 LMS算法的基本原理是基于最速下降法,即沿着权值 的梯度估值的负方向搜索,达到权值最优,实现滤波器的输 出信号与期望输出信号之间的LMS误差。自适应滤波的 原理框图 如图l所示。 图1 自适应滤波器原理框图 Fig 1 Principle diagram of adaptive filter 图1中, (n)为n时刻的输入信号矢量; (n)为输出 信号;v(n)为噪声信号;d(n)为期望输出信号;P(n)为d(n) 和.),(n)之间的误差信号估计;通过误差信号P(n)调节自适 应滤波器抽头权向量,使自适应滤波器收敛至稳定状态。 基于最速下降法的LMS算法公式 如下 e(n):d( )一 。(n)w(n), (1) II,(n+1)=’‘r(n)+ F(n) (n), (2) 式中 (n)为 时刻Ⅳ阶自适应滤波器的权系数; 为控 制稳定性和收敛性能的参量即步长因子。LMS算法收敛 时,步长因子 的取值范围为:0< <1/3.…,Amox为输入信 号自相关矩阵的最大特征值。 2新的变步长LMS算法与性能分析 变步长LMS算法的基本思想 。川是:在初始收敛阶段 或者系统参数发生时变的时候,最优权值与自适应滤波器 的权值相距较远,为了保证能有较快的收敛速度和对时变 系统的跟踪速度,选取较大的步长/Z;在算法接近收敛时, 滤波器的权值接近最优权值,为了减少算法的稳态误差,选 取较小的步长肛。基于这种思想,本文在前人研究的基础 上,通过在步长参数 (n)与误差信号e(n)之问建立一种 新的非线性函数关系来调节步长,提出一种新的变步长 LMS算法,本文算法步长因子函数和权系数更新公式如下 (n) {l—exp[一 I e(n)e(n--1)I]}, (3) w(n+1)=w(n)+ (n)P(n) (n)+ (n)(I (n)1—1 P(n-1)I), (4) (n)=py(n-1)一・=p A(n), (5) 式中 为控制函数形状的常数,决定曲线上升的快慢; 为控制函数取值范围的常数。Ot, 分别变化时步长因子 (n)和误差e(n)的关系曲线如图2和图3所示。 由图2可知,当Ot固定时, 取值越大,步长初始值也 越大,收敛速度也越快,同时稳态误差也越大;卢取值越小, 步长初始值也越小,收敛速度也越慢,稳态误差也越小。由 图3可知,当 固定时, 取值变化时也具有同样的特点, e(n) 图2 =1,卢变化时 ( 与e(,1)的关系曲线 Fig 2 Relation curYes with肛( )and el,I)when 01=1 and diferent P( ) 图3 p=1, 变化时 (,1)与e(n)的关系曲线 Fig 3 Relation curves with (,1)and P(n)when卢=1 and diferent 同时,在I e(n)I>1.3时步长变化不明显,但 越大时步长 也越大,收敛速度也越快。e(n)趋近于0时,步长变化也较 缓慢。由以上分析可知,新算法稳态时步长变化平滑,克服 了SVSLMS算法e(n)接近0时步长变化太大的缺点。算法 性能由 和口共同决定,当要求较快的收敛速度时,2个参 数的取值都应该较大,当要求较小的稳态误差时,2个参数 的取值都应该较小,在实际应用中根据环境的不同来确定 最佳取值。 式(4)中的ky(n)(f e(n)I—I e(n-1)I)为本文算法 绝对估计误差的扰动量。传统算法采用随机梯度调整抽头 权向量,每次迭代时的估计误差e(n)的正负并不确定,抽 头权向量的调整只能围绕 震荡。本文用绝对估计误差 能更好地表示估计信号偏离期望信号的程度,扰动量可以 对w(n)进行正反向调节来加快算法的初始收敛速度。通 过 (n)扰动量幅度因子,把扰动量对 (n)的调节控制在 一个最佳水平。 (n)以式(5)的指数形式衰减,其中,P< l。在迭代次数比较少时,扰动量对 (n)的影响较大,随着 迭代次数的增加, (n)趋近于0,e( )的波动几乎不会对 ’.,(n)造成影响,从而使稳态误差抑制在较低的水平。为了 加强算法的抗噪性能,本文算法中用式(3)的误差向量自相 关值e(n)e(n--1)来调节步长。设自适应滤波器抽头权向 量的维纳解为W。,令△w='.,(n)一w。为n时刻权系数矢量 与最佳值之差,期望信号的介入噪声为 (n),与 (n)不相 关,则有 d(n)= (n)w0+ (n), (6) d(n)=P(n)+ (n)’I,(n). (7) 第3期 贺洪江,等:一种新的LMS自适应滤波算法分析仿真研究 17 由式(6)、式(7)得 e(n)= (n)一 (n)Aw, (8) e(n)e(n--1)=孝( )孝(n--1)一 (n)x (n-1)Aw— Aw x(n) (n--1)+Aw (,Ox (n—1)Aw. (9) 设噪声功率为 ,由系统特性可知 E Ex(n)g(m)] 0, (10) E Eg(n)g(m)]= 。艿(n--m), (11) E[e(n)e(n-1)]=E[Aw X(n)x ( 一1)Aw],(12) E[e ( )]= [Aw X(n)x (n一1)Aw]+ . (13) 由式(12)和式(13)可知,误差自相关函数只与输入信 号有关,而与噪声无关。SVSLMS算法用误差调节步长,G— SVSLMS算法用误差功率调节步长,导致在低信噪比条件下 性能恶化。本文算法的噪声抑制能力优于上述算法。 3算法仿真结果分析比较 3.1仿真条件 1)自适应滤波器阶数L=2; 2)未知系统的FIR系数为W =[0.8,0.5r; 3)参考输入信号 (n)是均值为0,方差为1的高斯白 噪声; 4)v(n)为与x(n)不相关的高斯白噪声,其均值为0, 方差为l; 5)在本文实验条件下对三种算法进行了大量的仿真, 测定了三种算法性能最优时的参数取值范围,参数取值如 下:SVSLMS算法中OL=150, =0.07;G—SVSLMS算法中 =450,卢=0.06;本文算法中 =1500,卢=0.05,',(0)= 0.4,P=0.4: 6)采样点数为1000,分别做500次独立仿真,然后通过 求其统计平均,得出学习曲线。 3.2仿真结果分析比较 本文算法和SVSLMS算法及G—SVSLMS算法收敛曲线 的比较如图4所示。 迭代次数 图4三种不同算法收敛曲线 Fig 4 Convergence curve of three diferent algorithm 图4中的三条收敛曲线从上到下依次为SVSLMS算 法、G—SVSLMS算法和本文算法的收敛曲线。由图3可知, 本文算法的收敛速度快于SVSLMS算法和G—SVSLMS算 法,稳态误差小于SVSLMS算法和G—SVSLMS算法的稳态 均方误差较接近,本文算法有效抑制了随机噪声对信号的 干扰,在综合性能上优于SVSLMS算法和G—SVSLMS算法, 验证了前文对算法性能的分析。 4系统发生时变时算法分析 系统发生时变时的总采样点数为1000,系统在第 500个采样点时刻未知系统发生时变,未知系统的FIR系数 由原来的W =[o.8,0.5 r变为W :[0.7,0.5,0.2r,做 500次独立仿真,求其统计平均,得出系统发生时变时本文 算法收敛曲线如图5所示。 迭代次数 图5本文算法发生时变时收敛曲线 Fig 5 Convergence curve when this algorithm have time—varying 由图可以看出:本文算法在系统发生时变时同样具有 很快的收敛速度,因此,本文算法在综合性能上优于SVS— LMS算法和G—SVSLMS算法的同时仍具有较好地对时变系 统的跟踪能力。 5结论 本文通过在步长因子与误差信号之间建立一种新的非 线性函数数,同时引入绝对估计的扰动量,提出了一种新的 变步长LMS自适应滤波算法,较好地解决了收敛速度和稳 态误差之间的矛盾。理论分析与仿真验证显示,本文算法 不仅初始收敛速度、稳态误差和抗噪性能均优于传统的几 种算法,而且具有较好的对时变系统的跟踪能力。 参考文献: [1]Diniz P S R.自适应滤波算法与实现[M].刘郁林,译,北京: 电子工业出版社,2004:48. [2]Diniz P S R.Adaptive filtering:Algorithms and practical imple— mentation[M].2nd ed.USA:Spring,2002. [3]叶水生,余荣贵,吴霄,等.一种新的自适应最小均方算法 及其应用研究[J].电测与仪表,2008,45(511):9. [4]覃景繁,欧阳景正.一种新的变步长自适应滤波算法[J].数 据采集与处理,1997,12(3):171--174. [5]高 鹰,谢胜利.一种变步长LMS自适应滤波算法及分 析【J].电子学报,2001,29(8):1094--1097. [6] 李方伟,张浩.一种新的变步长LMS自适应滤波算法及其 仿真[J].重庆邮电大学学报:自然科学版,2009,21(5): 591-594. [7]Irusta U,Ruiz de G S,Ruiz J,et a1.A variable step size LMS algo— rithm for the suppression of the CPR artifact from a VF signal[J]. Computers in Cardiology,2005,48(4):179--182. (下转第21页) 第3期 唐骥锋,等:基于多传感器融合的机器人蒙特卡洛定位决策 21 g g。\删曾}咖丌J匠 辑 仿真实验证明:SR—UKF—MCL算法在非线性复杂环境中能 霸 ∞ ∞z 极大提高机器人定位的抗干扰能力,减小定位误差。 删 5 4 3 2 O 参考文献: -叵 椒 暴 [1]宋骊平,姬红兵.多被动传感器UKF与EKF算法的应用与比 较[J].系统工程与电子技术,2009,31(5):1083--1086. 时间,s [2] Julier S,Uhhnann J,Hush F.A new method for the nonlinear 图3沿 轴方向的位置误差 transformation of means and covariance in filters and estimator. Fig 3 Position error in X-axis direction s[J].IEEE Transaction on Automatic Control,2000,45(3): 222--226. [3] 吴盘龙,孔建寿.基于平方根UKF的水下纯方位目标跟 踪[J].南京理工大学学报:自然科学版,2009,33(6):751— 755. [4] 陈小宁,黄玉清,杨佳.多传感器信息融合在移动机器人定 时间/s 位中的应用[J].传感器与微系统,2008,27(6):110--113. 图4沿l,轴方向的位置误差 [5]赵璇,何波,吉得志,等.基于粒子滤波器的机器人定位 及动态目标追踪[J].系统仿真学报,2008,2O(23):3715— Fig 4 Position error in Y-axis direction 3723. e 。 [6]周 华.传感器融合技术在移动机器人定位中的应用研 究[D].武汉:武汉理工大学,2009. .叵 5 [7] 赵璇.基于粒子滤波的机器人定位及动态目标追踪[D]. 青岛:中国海洋大学,2008. 0 [8]Delleart F,Fox D,Burgard W,et a1.Monte Carlo localization for 时间/s mobile robots『J].IEEE Transactions on Robotics and Automa— 图5定位的方向角误差 tion,1999,2(5):10-15. Fig 5 Direction angle error of localization [9] 潘勃,冯金富,李骞,等.毫米波/红外多传感器融合跟踪 算法研究[J].红外与毫米波学报,2010,29(3):230-235. [1O]Van der,Merwe R,Wan E A.The square-root unscented Kalman filter for state and parameter—estimation[c]∥Processing of the IEEE International Conference on Acoustics,Speech,and Signal 5O 1OO Processing,NJ,USA:IEEE,2001:3461--3464. 迭代次数 [11]Sadhu S,Modndal S,Srinivasan M,et a1.Sigma point Kalman fil— 图6机器人NEES平均值变化曲线 ter for beating only tracking[J].Signal Processing,2006, Fig 6 Average value change curve of the robot NEES 86(12):3769-3777. 由图6可知,SR.UKF—MCL算法的NEES平均值基本 [12]Mortemerlo M,Thrnn S,Koller S T D,et a1.FastSLAM:A fac- 上都在要求之下,符合一致性条件,并且其一致性要优于 tored solution to the simuhaneous focalization and mapping prob— EKF—MCL和UKF—MCL。 lem[C]∥Proc of the AAAI National Conf on Artiifcial Intelli- 5结束语 gence,Edmonton,Canada:AAAI,2002:593--598. 通过将多传感器信息进行平方根无迹卡尔曼滤波融合 作者简介: 之后的机器人蒙特一卡洛定位(SR—UKF—MCL)决策,较之经 唐骥锋(1986一),男,河南济源人,硕士,主要从事基于多传 典卡尔曼滤波决策,在机器人定位精度上有很大提高。 感器融合移动机器人决策与规划的研究。 、t; 、0; 0 、 、0 、 、 j ) ) ) > ) ) (上接第17页) 算法[J].电子学报,2006,34(6):1123--1126. [8]Mohammad z U R,Ahamed S R,Koti R D V R.Cancellation of [1O]Givens M.Enhanced convergence normalized LMS algorithm[J]. artifacts in ECG signals using sign based normalized adaptive fil— IEEE Signal Processing Magazine,2009,26(3):81-95. teifng technique[C]//2009 IEEE Symposium on Industrila Elec— tronics and Applications,ISIEA 2009,Kuala Lumpur,Malaysia, 作者简介: 2009:442—445. 贺洪江(1964一),男,河北邯郸人,教授,硕士生导师,主要研 [9] 罗小东,贾振红,王强.一种新的变步长LMS自适应滤波 究方向为计算机控制技术、传感器技术和信号处理技术。