一种改进的变步长LMS自适应滤波算法及性能分析

  • 格式:pdf
  • 大小:1.01 MB
  • 文档页数:3

收稿日期:2010 07 28;修回日期:2010 08 31作者简介:任自钊(1983 ),男,河南禹州人,硕士研究生,主要研究方向为基于SOPC IP 的设计方法与设计技术、射频功率放大器的研究(z i zhao218@126.co m );徐建城(1957 ),男,四川成都人,教授,硕导,主要研究方向为计算机网络、微电子技术;闫永鹏(1986 ),男,硕士研究生,主要研究方向为集成电路设计.一种改进的变步长L MS 自适应滤波算法及性能分析任自钊,徐建城,闫永鹏(西北工业大学电子信息学院,西安710129)摘 要:针对现有L M S (leastm ean squar e)算法不能同时提高收敛速度及降低稳态误差的矛盾,提出一种改进的变步长L M S 算法,建立了步长参数 (n)与误差信号e(n)之间的一种新的非线性函数关系:与现有的算法相比,同时引入记忆因子 和控制函数取值的参数 (n ),使当前步长与上一次迭代所得步长及前M 个误差的平方相关。

理论分析和计算机仿真结果表明,与现有几种常见的L M S 算法相比,改进的算法收敛速度和稳态误差的性能指标得到了提高。

关键词:自适应滤波;变步长;最小均方误差算法中图分类号:T N 911.72;TP301 6 文献标志码:A 文章编号:1001 3695(2011)03 0954 03do:i 10.3969/.j i ssn .1001 3695.2011.03.046I mproved variab le step size L M S adapti ve filteri ng a l gorith m andits perfor m ance anal ysi sREN Z i zhao ,XU Jian cheng,YAN Y ong peng(C olle ge of E lec t ronics&Infor ma tion,N orthw est ern P ol y t echnical University ,X i an 710129,Ch i na )Abstract :For existi ng LMS al gorith m could not sm i u ltaneousl y m i prove t he convergence speed and l ow er steady state error ofcontradicti ons .By bu il d i ng a nonli near functi onal relationsh i p bet w een (n)and e(n ),t h i s paper proposed an m i proved vari able step size LMS adapti ve filteri ng algori th m.Co m pared w it h existi ng al gorit h m s ,wh ile the introduction ofm e m ory factor and control f uncti ons val ues of the para m eters (n ),so that the current iteration step w ere related t o t he previous step and the for m er M the square of error .Theoretical anal ysis and compu ter sm i ulati ons show that several comm on w ith exi sti ng L M S al go rith m,t he m i proved convergence rate and steady state error perfor m ance is m i proved .Key words :adapti ve filteri ng;vari able step size ;L M S algorith m0 引言自适应滤波广泛应用于自适应控制、噪声消除、系统辨识、雷达和信号处理等领域。

由W i drow 和H off 提出的最小均方误差(LM S)算法,因其结构简单、鲁棒性好、易于实现,成为应用最广泛的一种自适应算法。

收敛速度、时变跟踪能力、稳态失调与算法复杂度是衡量自适应滤波算法性能的四个主要技术指标。

但是,传统的算法使用固定步长,导致收敛速度和稳态误差性不能同时得到满足。

为克服这一缺点,人们提出许多改进的LM S 算法,其中大部分是变步长L M S 算法。

文献[1]给出了S i g m o i d 函数变步长L M S 算法,该算法虽能同时获得较快的收敛速度、跟踪速度和较小的稳态误差,但当达到稳态时仍有较大步长变化,不具有缓慢变化的特性;文献[2~4]在si gm oid 函数基础上分别对e (n)引入了2次方、|e(n)e(n -1)|及3次方,克服了SVSLM S 算法在自适应稳态阶段步长调整过程中的不足;文献[5~7]提出了三种与误差信号成非线性关系的步长设计方法,这些算法均引入多个调整参数,因此步长因子不易设计和控制;文献[8]在文献[6]的基础上引入了测量噪声v 的方差,使得步长能够跟随测量噪声及时变化。

文献[9]引入记忆因子,使步长与前n 个误差值相关;文献[10]提出了一种新的基于双曲正切函数的变步长算法。

文献[8~10]都得到了较好的性能,但是这些算法在实现时,增加了计算复杂度。

为克服上述变步长L M S 自适应滤波器存在的不足,本文提出了一种新的变步长L M S 自适应滤波算法。

该算法使用的变步长因子兼顾了上一次迭代的步长和前M 个误差信号的具体变化过程,能够很好地匹配自适应滤波器的收敛特性,而且求变步长因子时计算量较小。

1 改进的变步长L M S 自适应滤波算法自适应滤波的原理如图1所示。

图中,x (n)是时刻n 的输入信号矢量,通过未知系统后受到测量噪声v (n )的干扰,生成d (n )。

W (n )为时刻n 的自适应滤波器权矢量,为了辨识未知或时变系统的抽头值,将x (n )通过自适应滤波器,根据LM S第28卷第3期2011年3月计算机应用研究Application R esearc h of C o m putersV o.l 28N o .3M ar .2011算法,自适应地调整该滤波器抽头权值,使得输出误差e(n)的均方值最小化。

经迭代收敛后得到的W即为未知系统的辨识模型。

根据步长调整的原则:初始收敛阶段,步长较大,以得到较快的收敛速度;在收敛阶段,步长应较小,以达到较低的稳态失调。

本文提出的变步长LM S自适应滤波算法如下:(i)=exp(-2i)(i=0,1,2, ,M-1)(1)其中, 为指数形式的记忆因子。

e~2(n)= M-1i=0(i)e2(n-i)(2)式(2)利用记忆因子 对过去M个误差功率加权,越旧的信息对现在的影响就越小。

(n)= (n-1)+(1- )e~2(n)(3)式(3)是利用过去m个误差功率的加权求出 (n)。

(n)= (n)(1-exp(- (n)|e(n-1)e(n)|))(4)利用式(3)中求出的 (n),代入式(4),从而求出步长 (n)的值。

e(n)=d(n)-X T(n)W(n)(5)式(5)是求滤波值和期望值的误差W(n+1)=W(n)+ (n)e(n)X(n)(6)利用上面求出 (n)、e(n),代入式(6),从而求出滤波系数。

其中:W(n)为自适应滤波器在时刻n的权值向量;X(n)为输入信号向量;d(n)为期望信号;e(n)为误差信号; (n)是控制函数形状的函数; (n)是控制函数取值范围的函数; (0< < 1)是相关因子,用于描述 (n)及 (n-1)的相关度; (n)为步长因子,可控制稳态性能和收敛速度,由 (n)和 (n)共同控制。

综上所述,本算法用记忆因子和误差的乘累加来控制步长,这样首先保证步长不会在收敛前就降到最小,同时也具有抗噪性能。

2 算法分析2 1 变步长因子分析根据步长调整原则,初始阶段步长最大,设 (0)=opt,opt为定步长LM S算法达到最快收敛速度时的步长。

由式(4)可知, (n)由 (n-1)及误差功率共同求得。

在初始收敛阶段,由于W(n)与未知系统权值向量差别较大,使得e(n)较大,令 接近于1,如0.99,则由式(4)及 (0)可得到下式成立:(n)o pt(7)在自适应滤波迭代收敛阶段,由于W(n)与未知系统权值向量近似相等,使得e(n)很小,由式(2)~(4),可知下式成立:(n)<o pt (8)式(7)(8)说明,在自适应滤波器迭代的收敛阶段,本算法能够根据输出误差的大小,匹配地、自适应地减小变步长因子的值,从而达到较低的稳态失调。

2 2 收敛条件分析一般的,算法收敛的步长因子范围是0< <2/m ax,此处引用文献[11]中给出的更精确的L M S算法收敛条件:0< <(1/3)tr R(9)其中,tr R为输入信号向量X(n)自相关函数的迹。

由式(7)(8)可知,要使本算法收敛,只需要满足0<opt< (1/3)trR,即0< (0)<(1/3)trR(10)2 3 抗干扰性分析将误差式(5)改写为d(n)=e(n)+X T(n)W(n)(11)因为误差e(n)与输入信号X(n)相关,此处将参考信号d(n)表示为d(n)= (n)+X T(n)W*(n)(12)其中: (n)是零均值噪声,与输入噪声无关;W*(n)为自适应滤波器最优权矢量。

令:K(n)=W(n)-W*(n)(13) K(n)表示系数偏差矢量。

则由式(7)~(9)可得e(n)e(n-1)=K T(n)X(n) (n-1)+K T(n)X(n)X T(n-1)K(n)+(n)X T(n-1)K(n)+ (n) (n-1)(14)由于 (n)是均值为0的噪声,并且 (n)与X(n)无关,故有E[e(n)e(n-1)]=E[K(n)T X(n)X T(n-1)K(n)](15)由式(15)可看出,本文采用e(n)e(n-1)调整步长后,降低了噪声 (n)对步长的影响,噪声可以忽略不计。

由于改进算法步长只与输入信号X(n)有关,与噪声 (n)无关,具有收敛速度快、稳态误差小的优点,且在低信噪比的环境中仍保持较好的性能,故有更广泛的用途。

2 4 参数 和 的选择及其对算法性能的影响参数 和 的不同选择,将对算法的收敛速度和稳态失调产生影响。