自适应滤波LMS算法及RLS算法及其仿真

1、自适应滤波第1章绪论11.1自适应滤波理论发展过程11 .2自适应滤波发展前景21. 2.1小波变换与自适应滤波21 .2.2模糊神经网络与自适应滤波3第2章线性自适应滤波理论42. 1最小均方自适应滤波器42.1.1最速下降算法42.1.2最小均方算法62.2递归最小二乘自适应滤波器7第3章仿真123.1 基于LMS算法的MATLAB真123.2 基于RLS算法的MATLAB真15组别：第二小组组员：黄亚明李存龙杨振第1章绪论从连续的(或离散的)输入数据中滤除噪声和干扰以提取有用信息的过程称为滤波。相应的装置称为滤波器。实际上，一个滤波器可以看成是一个系统，这个系统的目的是为了从含有噪声的数

2、据中提取人们感兴趣的、或者希望得到的有用信号，即期望信号。滤波器可分为线性滤波器和非线性滤波器两种。当滤波器的输出为输入的线性函数时，该滤波器称为线性滤波器，当滤波器的输出为输入的非线性函数时，该滤波器就称为非线性滤波器。自适应滤波器是在不知道输入过程的统计特性时，或是输入过程的统计特性发生变化时，能够自动调整自己的参数，以满足某种最佳准则要求的滤波器。1.1自适应滤波理论发展过程自适应技术与最优化理论有着密切的系。自适应算法中的最速下降算法以及最小二乘算法最初都是用来解决有/无约束条件的极值优化问题的。1942年维纳(Wiener)研究了基于最小均方误差(MMSE涯则的在可加性噪声中信号的最

3、佳滤波问题。并利用Wiener.Hopf方程给出了对连续信号情况的最佳解。基于这准则的最佳滤波器称为维纳滤波器。20世纪60年代初，卡尔曼(Kalman)突破和发展了经典滤波理论，在时间域上提出了状态空间方法，提出了一套便于在计算机上实现的递推滤波算法，并且适用于非平稳过程的滤波和多变量系统的滤波，克服了维纳(Wiener)滤波理论的局限性，并获得了广泛的应用。这种基于MMSE则的对于动态系统的离散形式递推算法即卡尔曼滤波算法。这两种算法都为自适应算法奠定了基础。从频域上的谱分析方法到时域上的状态空间分析方法的变革，也标志着现代控制理论的诞生。最优滤波理论是现代控制论的重要组成部分。在控制论的

4、文献中，最优滤波理论也叫做Kalman波理论或者状态估计理论。从应用观点来看，Kalman波的缺点和局限性是应用Kalman波时要求知道系统的数学模型和噪声统计这两种先验知识。然而在绝大多数实际应用问题中，它们是不知道的，或者是近似知道的，也或者是部分知道的。应用不精确或者错误的模型和噪声统计设计Kalma砒波器将使滤波器性能变坏，导致大的状态估计误差，甚至使滤波发散。为了解决这个矛盾，产生了自适应滤波。最早的自适应滤波算法是最小JY(LMS)算法。它成为横向滤波器的一种简单而有效的算法。实际上，LMSf法是一种随机梯度算法，它在相对于抽头权值的误差信号平方幅度的梯度方向上迭代调整每个抽头权值

5、。1996年Hassibi等人证明了LMSf法在H。准则下为最佳，从而在理论上证明了LMSf法具有孥实性。自Widrow等人197阿提出LM咱适应滤波算法以来，经过30多年的迅速发展，已经使这一理论成果成功的应用到通信、系统辨识、信号处理和自适应控制等领域，为自适应滤波开辟了新的发展方向。在各种自适应滤波算法中，LMST法因为其简单、计算量小、稳定性好和易于实现而得到了广泛应用。这种算法中，固定步长因子以对算法的性能有决定性的影响。若以较小时，算法收敛速度慢，并且为得到满意的结果需要很多的采样数据，但稳态失调误差较小：当“较大时，该算法收敛速度快，但稳态失调误差变大，并有可能使算法发散。收敛速

6、度与稳态失调误差是不可兼得的两个指标。以往的文献对LMSJ法的性能和改进算法已经做出了相当多的研究，并且至今仍然是一个重要的研究课题。另一类重要的自适应算法是最小二乘（LS）算法。LSB法早在179成就由高斯提出来了，但LSJT法存在运算量大等缺点，因而在自适应滤波中一般采用其递推形式一一递才t最小二乘（RLS）算法，这是一种通过递推方式寻求最佳解的算法，复杂度比直接LS#法小，因而获彳#了广泛应用。199的SayedffiKailath建立了Kalman波和RLSS法之间的对应关系，证明了RLST法事实上是Kalman波器的一种特例，从而使人们对RLST法有了进一步的理解，而且Kalman波

7、的大量研究成果可应用于自适应滤波处理。这对自适应滤波技术起到了重要的推动作用。基本上，现有的参考文献都是基于这两种算法进行改进的。而这两种算法又可以简单的用以下语句来描述：LMS（抽头权向量更新值尸（老的抽头权向量值）+（学习速率）（抽头输出向量）（误差信号）RLS（状态递推值尸（旧的状态值）+（卡尔曼增益）（新息向量）以往的研究多集中在线性滤波方面，非线性滤波理论还有待于进一步的研究开发。1. 2自适应滤波发展前景现代信号处理理论为自适应滤波技术的发展提供了广阔的空间。尤其是小波技术和人工智能理论的发展，更是推动和加快了自适应滤波技术的前进。1.2. 1小波变换与自适应滤波小波变换是由法国地

8、球物理学家Morlet于80年代初在分析地球物理信号时作为种信号分析的数学工具提出来的。通俗地讲，小波是一种短期波。在积分变换中，小波作为核函数的用法大体与傅立叶分析中的正弦和余弦函数或与沃耳什（Walsh）分析中的沃耳什函数的用法相同。目前，小波分析主要用于信号处理、图像压缩、次能带编码、医学显像、数据压缩、地震分析、消除噪声数据、计算机图像、声音合成等领域。小波变换的基本特点是多分辨率或多标度的观点，目的是“既要看到森林（信号的概貌）,又要看到树木（信号的细节）”。借助于小波的精辟理论，自适应滤波技术有了新的发展方向，这也引起了信号处理领域许多学者专家的浓厚兴趣和热切关注。基于小波变换的自

9、适应滤波技术是未来自适应滤波发展的方向之一，有着广阔的应用前景。目前还有许多问题亟待解决，例如不同形式的小波滤波器的滤波效果研究；在时变信号滤波方面的应用研究以及对于失调噪声的滤波等等。1. 2.2模糊神经网络与自适应滤波神经网络作为一种新的计算方法，已经引起了人们广泛的研究兴趣。神经网络可以认为是一种由许多称为“神经元”（neuron）的基本计算单元通过广泛的连接所组成的网络。它是在现代神经科学研究成果的基础上提出来的，反映了人脑功能的基本特征。网络的信息处理由神经元之间的相互作用来实现，网络的学习与识别决定于各神经元之间联接权系数的动态变化过程。人工神经网络是模仿和延伸人脑认知功能的新型智

10、能信息处理系统，由于神经元本身具有高度自适应性，因而由大量神经元组成的神经网络具有自学习性、自组织性、巨量并行性、存储分布性、结构可变性等特点，能解决常规信息处理方法难以解决或无法解决的问题。模糊技术也是现代智能理论的一个重要方面。利用模糊技术我们可以很容易的将人们熟悉的语言描述应用到自动控制或信号处理中来，然而，模糊理论与神经网络都存在着各自的优缺点。模糊逻辑和神经网络在许多方面具有关联性和互补性。它们的交叉研究正是基于二者互补性和关联性的结合。首先，两者具有互补性。一方面，模糊技术的特长在于模糊推理能力，容易进行高阶的信息处理。将模糊技术引入神经网络，可以大大拓宽神经网络处理信息的范围和能

11、力，使其不仅能够处理精确信息，也能够处理模糊信息或其它不精确信息，不仅能够实现精确性联想及映射，也能够实现不精确性联想及映射，特别是模糊联想和模糊映射。另一方面，神经网络在自学习和自动模式识别方面有极强的优势，采取神经网络技术进行模糊信息处理，则使得模糊规则的自动提取及隶属函数的全自动生成有可能得以解决。其次,两者具有关联性，有许多共同点。它们都着眼于模拟人的思维，都是为了处理实际中不确定性、不精确性等引起的系统难以控制等问题。两者在形式上有不少相似之处，其信息都是分布式存储于其结构之中，从而都具有好的容错能力。不管是神经网络还是模糊逻辑，都不需要建立数学模型，只需根据输入的采样数据去获取所需

12、的结论，也就是模型无关估计器。另外，神经网络的映射功能早已得到证明。近年来，Kosko、L。X.Wan第证明了模糊系统能以任意精度逼近紧密集上的实连续函数，这也说明了二者之间有着密切的关系。因此，将两者融合在一起的模糊神经网络可以有效的克服两者的缺点，提高整个网络的性能。譬如；神经网络可以降低透明程度，使它们更接近于模糊系统；而模糊系统可以提高自适应性，更接近于神经网络。神经网络与模糊神经网络有力的推动了自适应滤波技术，特别是非线性自适应滤波技术的发展。事实上，一个神经网络或模糊神经网络本身就可以看作是一个自适应滤波器。第2章线性自适应滤波理论线性自适应滤波理论比非线性自适应滤波理论发展的比较

13、早、比较成熟。这里只重点介绍基本的最小均方(LMS)年tl最小二乘(LS)自适应滤波器。2. 1最小均方自适应滤波器本节所讨论的lmst法是应用最广泛的一类算法，从理论体系上来看，最优化方法中的最速下降法是随机梯度信息处理的一种递归算法，在误差性能未知的情况下，它可以寻求误差曲面的最小点，可为平稳随机条件下的lmst法提供若干启发性思路，并且是分析算法性能的重要基础。2. 1.1最速下降算法最速下降法是一种不用求矩阵逆来解正规方程组的方法。它通过递推方式寻求加权矢量的最佳值。虽然在自适应滤波中很少直接使用最速下泽算法，但它构成了其它自适应算法、特别是LMST法的基础。图2.1表示一个横向滤波器

14、。图2.1横向滤波器输入矢量为x(n)=lx(n,x(n-1),.,x(n-M+1)(21)加权矢量(即滤波器参数矢量)为w(n)=l.w0(n),w(n)1,.,Wm(n)T(22)滤波器输出为my(n)=£Wj(nk(n-i+1)=wT(nk(n)(23)i1利用图2-1中输出信号y(n)与期望信号的d(n)的关系，误差序列e(n)可以写成:e(n)=d(n)-y(n)(24)自适应滤波器就是根据误差序列8(行)按照某种准则和算法对其系数w(n)进行调整，最终使自适应滤波器的目标(代价)函数达到最小，即最佳滤波状态。按照均方误差(MSE潍则所定义的目标函数是F(e(n)=己(n)

15、=Ee2(n)=Ed2(n)-2d(n)y(n)+y2(n)(2-5)将式(23)代入上式，得到E(n)=Ed2(n)-2Ed(n)wT(n)x(n)+EwT(n)x(n)xT(n)w(n)(2-6)固定滤波器系数，则目标函数(2.6)可写为(n)=Ed2(n)-2wTp+wTw(2.7)式中，R=Ex(n)xT(n)是输入信号的自相关矩阵；P=Ed(n)x(n)是期望信号与输入信号的互相关矢量。假设输入自相关矩阵R为非奇异的，当R与P均已知时，将式(2.7)对W求导数，并令其等于零，可得到使目标函数最小的最佳滤波参数W为W=R1P(2-8)这个解称为维纳解，即最佳滤波参数值。从式(2.7)可

16、以看出，自适应滤波器的目标函数是滤波参数W勺二次函数，因此形成了一个多维的超抛物曲面，二维时好像是一个碗状的曲面且具有唯一的碗底最小点。这个多维的超抛物曲面通常称之为自适应滤波器的误差性能曲面。位于该曲面上的某一点，经过自适应调节过程，可以朝碗底最小方向移动，最终到达最小点。最速下降法就是实现这种从初始值到最佳值搜索的一种优化技术，它利用梯度信息分析白适应滤波性能和最佳滤波状态，避免了对输入信号自相关矩阵R直接求逆。可以想象，沿着己减小的方向调整权值W应该可以找到最佳值W,因为梯度的方向是增长最快的方向，所以负的梯度方向就是己减少最快的方向，这样，就可以采用如下的递推公式来调整WZ寻找W:(2

17、-9)w(h+1)=w(n)+式中，()代表珂时刻孝的梯度，是一个MX1维的矢量，这里浊滤波器滤波权系数的数目；u是一个正实常数，通常被称为步长或步长因子。根据梯度定义，(n)可以写成：»爵e喘(MO)(Ml)=-£2e(n)x(H)J将式(2.7)对W(偏导，可以得至IJ:(神)m-2P+2Rw(n)将式(211)代入式(2-9)可得到最速下降法的递推公式：w(n-hl)=w(n)+2P-RwQO(2-12)有关最速下降法的收敛性，这里只给出其结果。当最速下降法满足下面条件时是收敛的：0<<Xmax式中，入ma是自相关矩阵R的最大特征根O2.1.2最小均方算法

18、如上节所述，要使用最速下降法，就要知道均方误差性能函数的梯度的精确值，见式(211),这就要求输入信号x(n)和期望信号d(n)平稳，并且要求它们的二阶统计特性已知。而这是相当复杂的，很多情况下它们是未知的或不完全知道的，因此一般采用梯度的估计值(n)来代替梯度(n),即(M4)+1)=w(w)-4十(叫最小均方(LMS)算法就是使用瞬时输出误差功率的梯度Ee2(n)作为均方误差梯度Ee2(n)的估计值，也就是令(2-15)(2-16)行(切二干讥/伽月=上“9将式(215)代入式(2-14),有w(w+1)=w(n)一=w(n)-口“(")-wr(«)x(rt)2=vv(

19、n)+2e(w)x(n)式(2.16)就是LMST法的递推公式。自适应lmst法简单，既不需要计算输入信号的相关函数，又不需要求矩阵的逆，因而得到了广泛应用。但是，由于LMSJ法采用梯度矢量的瞬时估计，它有很大的方差，以至不能获得最优滤波性能。下面给出完整的算法步骤：有关参量：M-滤波器抽头数；a步长；0<n<(MRO-1,Pn=Ex1(n)2初始条件：w(O)=0或者由先验知识确定计算步骤：对于H=1,2,?取x(n),d(n)(2)滤波y(n尸WT(n)x(n)估计误差e(n)=d(n)-y(n)更新权向量w(H+1)=w(n)+e(n)x(n)当滤波器的输入信号为有色随机过程

20、时，特别是输入信号为高度相关的情况，大多数自适应滤波算法的收敛速度都要下降，对于上述典型LMSJ法。此问题更加突出。为了提高收敛速度，相继提出了多种改进的算法。下面给出几种改进的权值递推公式：归一化LMSf法：w(n+1)=w(n)+中而rx(»)(2-17)式中，0<n<1为控制失调的固定收敛因子；丫是为了避免XT(n)x(n)过小而导致步长值太大而设置的正常数。简化的lmst法:w(w+l)=w()+2附即(乂5)卜町(2'18>或+1)=w(h)+2jtix(M)5J?4(?7)(2'19)或wS+1)=w(rt)+2?igwx(n)i(2-2

21、0)式中，signx)="2。-10产或gn(x)=*忖，工*00/=0(3)MLMSU法:w(w+i)=w(rt)+S+l)£w(n+l)x(«+1)(2-2JJ式中，乂诃=川'1+卬6+1)、(界+1):(n+1)-d(H+l)-w,S)x(fi+D以上这几种改进LM算法中，前两种都是交步长方法，后一种采用的是改变梯度估值的方法。除了这几种方法以外，还有很多改进的算法，比如变换域LM算法，该算法适用于输入信号具有高度相关性的情况下；频域LM算法，该算法与经典梯度下降法相比较有着更好的收敛性；分块LMSJ法，与其它方法相比该算法计算量大为减少；最小高阶均

22、方(LMK)算法，该算法是LMSJ法的扩展,或者说LMS算法是该算法的特例，当系统噪声为非高斯分布时，LMKT法要比LMSJ法收敛精度高；Q的解LMSJ法，该算法适用于横向延迟线抽头数目比较小的情况。其它方法这里不再列举。2.2递归最小二乘自适应滤波器最小二乘(LS)法是一种典型的有效的数据处理方法，既可用于静态系统，又可用于动态系统：既可用于线性系统，又可用于非线性系统；既可用于离线估计，又可用于在线估计。递归最小二乘(RecursiveLeastSquare简写RLS)!最小乘法的一种快速算法，它包含时间递归最小二乘(TRLS)算法和阶数递归最小二乘(ORES算法两方面内容，一般前者适用于

23、动态系统辨识和在线估计，后者适用于静态系统辨识和离线估计。与LMST法相比，RLST法有着非常快的收敛速度。图2-2自适应RUM法采用RUM法的自适应滤波器如图2-2所示。输入信号x(n)，含有Nt已知样本x(1),x(2),x(3)-x(N),期望输出为d(n)=d(1),d(2),，d(N)。滤波器输出为MMRxZ吗工5-A+l)(月=12，N(2-22)M伪滤波器长度，且M<No误差信号为e(n)=d(n)-y(n)(2-23)按照最小平方准则，设计滤波器目标函数为00=£e2(i)=XK0-M02(2-24)J-Ijnl将式(2-22)代入式(224)得到45)=Zd&

24、#39;-2£呜S)£-左+1)f=ii=iMMN+X£/匕x(i-k+1)x("掰+1)(2-25)定义如下参量：一左)式,-w)：k.m=01-1(2-26)为确定性自相关函数，表示输入信号在抽头k与抽头於问两信号的相关性：何此k)=玳;k=OL、材-1(2-27)为确定性互相关函数，表示期望响应与在抽头嗣入信号之间的互相关性(2-28)45)=%(?)()为期望响应序列的能量。将式(2-26)(2-28)代入式(2.25),目标函数可写为&5)=Edn-2工w*5*(N;*-1)4«1(2-29)十£X"%(汽

25、源(N法一1,第一1)*-T用1因此有舞=询97噂水加)(2-30)令式(2-30)等于零，得到031)£为6心(”;4一Lm-i)=#一】卜A=,M明1二1把它简化成矩阵形式，有O(n)w(n)Me(/j)(2-32)式中，w(n)为MX1维最小平方估计的滤波器系数矢量，即wS)=吗5),w25),(w)f(2-33)(n)为延迟线抽头输入信号的确定性自相关函数，是一个MXMS矩阵即<前明0.0)。(胃:0.1)0(n;O,M-1)(2-34)奴风0)/(矶1)椁*只/一1)0(/1)二.*二百(w;A/-1,0)M-1J)故-1)9(n)为脉冲响应序列与输入信号之间的确定性

26、互相关函数MX1维矢量，即0(«)=阳(厘;0)-1),(2-35)假定矩阵(n)是非奇异的，其逆矩阵存在，则由(2-32)可求得w(w)=O)1(rt)0(77)(2-36)式f2,36)就是最+-乘法自适应滤波算法的正规方程。从式中可以看出，它要求(n)是可逆的矩阵，虽然如此，但是对大多数应用来说这都是成立的。若对于某应用(n)为降秩，即为不可逆矩阵，则式(2-36可理解为采用了伪逆矩阵。下面推导RUM法由式(2,36)有w(n-1)=n7)0伽-1)(2-37)根据式(2.26)可得到写成矩阵形式<!>(»)=。(万-i)+xS)x,()利用矩阵求逆定理(

27、A十BCD)-1=%"一47刁(。7+DA-lByDA'再由式(2-39)可求得力T(界-DxS)x丁也)T51)i+界)“伽-i)x(m(2-38)(2-39)(2-40)(2*41)引入MXM隹矩阵。42)P5)=T(/t)K5)=1+xr(«)P(«-l)x(n)(2-43)和N隹矢量K(n)被称为增益系数。将式(2.42)和式(2-43)代入(2-41)得到P(w)=PS1)-K5X，(h)P5-1)(2-44)上式两边右乘x(n)P(rt)x(n)=P(nl)x(w)-K(w)xr-l)x(n)(2-45)用式(2.43)中分母多项式右乘上式两边

28、，整理得到K(«)x;(n)P(rt-l)x(w)-P(/J-l)x(rt)-K(«)(2-46)(2,47)把上式带入式(2.45)中，可化简成下式K5)=pg)xs)另外，根据式(2.27)有如下递推关系式8®k)=e(k-1;*)+d(n)x(n-k)(2-48)写成矩阵形式Q(n)=e(K-D+45)x5)(2-49)式中MX1维矢量d(n)x(n)代表递归计算的更新校正项,把式(242)代入式(2-36)得到：w(n)=P(n)e(n)(2-50)再将式(2-49)代入上式，并由式(2-47)可得到w(n)=P(n)0(MT)+P(?i)d(/i)x(n

29、)P(n)6(n-1)+KSWpi)(2-51)将式(2-44)代入式(2-51)中，且根据式(2-50)有w(w)=P(n-l)0(w-1)一K5)x/5)PS-l)g(理一】)+K(n)d(m=w(-1)+K5)a(，)一xr(n)w(n-1)(2-52)上式就是rum法的递推公式。与LMsr法相比较，两者的主要差别在于增益系数，LMSJ法简单的利用输入矢量乘上常数，而RLSJ法则用较复杂的K(n)。下面给出RLSJ法的初始条件以及运算步骤：初始条件：w(O)=0;P(O)=-1(0)=6-1I为小的正实数运算步骤：对于n=1,2,?(1)取d(n),x(n)更新增益矢量K(n)=詈T)x

30、5)一(»IX，"更新滤波器参数=wS-D+K5)d(用-xJ(n)w(w-1)更新逆矩阵P5)=P(7-K5)x75)PS-1)第3章仿真3.1基于LMSJ法的MATLA的真在matlab中输入以下程序，可得出仿真结果g=100;%统计仿真次数为gN=1024；%输入信号抽样点数k=128;%滤波器阶数pp=zeros(g,N-k);%将每次独立循环的误差结果存于矩阵pp中，以便后面对其平均u=0.00026;%滤波器收敛因子forq=1:gt=1:N;a=1;s=a*sin(0.05*pi*t);%输入单信号sfigure(1);subplot(411)plot(s);

31、%信号s时域波形title('信号s时域波形');xlabel('n');axis(0,N,-a-1,a+1);xn=awgn(s,5);%加入均值为零的高斯白噪声%设置初值y=zeros(1,N);%输出信号yy(1:k)=xn(1:k);%将输入信号xn的前k个值作为输出y的前k个值w=zeros(1,k);%设置抽头加权初值e=zeros(1,N);%误差信号%用LMS算法迭代滤波fori=(k+1):NXN=xn(i-k+1):(i);y(i)=w*XN'e(i)=s(i)-y(i);w=w+u*e(i)*XN;endpp(q,:)=(e(k+1

32、):N)A2;endsubplot(412)plot(xn);title('信号s加噪声后的时域波形');subplot(413)plot(y);title('自适应滤波后的输出时域波形');subplot(414)plot(e);title('估计误差e(t)的时域波形');forb=1:N-kbi(b)=sum(pp(:,b)/g;endfigure(2);t=1:N-k;plot(t,bi,'r');title('收敛曲线)；holdon%信号s时域波形%信号s时域波形%求误差统计平均%算法收敛曲线%将每次循环的图形显示结果保存下来正弦信号5加噪%后的时域波形oViywvVWWVAK/yvV/WvNVwV1-5J11JL020040060080010001200自适应滤波后的输出时域波形nI!0ZljlVvA/WWWWWXAAA/VW收敛曲线01002003004005006007008009003.2基于RLSJ法的MATLA时真程序代码：clearN=1000;Fs=500;n=0:N-1;t=n/Fs;xs=(sin(2*pi*5*t)'figure；subplot(211);plot(t,xs);grid;ylabel('幅