RLS和LMS自适应算法分析报告.doc

wordRLS和LMS自适应算法分析摘要：本文主要介绍了自适应滤波的两种算法：最小均方(LMS, Least Mean Squares)和递推最小二乘(RLS, Recursive Least Squares)两种根本自适应算法我们对这两种根本的算法进展了原理介绍，并进展了Matlab仿真通过仿真结果，我们对两种自适应算法进展了性能分析，并对其进展了比拟用Matlab求出了LMS自适应算法的权系数，与其学习过程曲线，和RLS自适应权系数算法的学习过程关键词：自适应滤波、LMS、RLS、Matlab仿真Abstract: this article mainly introduces two kinds of adaptive filtering algorithms: Least Mean square (LMS), further Mean Squares) and Recursive Least Squares (RLS, Recursive further Squares) two basic adaptive algorithm. Our algorithms of these two basic principle is introduced, and Matlab simulation. Through the simulation results, we have two kinds of adaptive algorithm performance analysis, and carries on the parison. Matlab calculate the weight coefficient of the LMS adaptive algorithm, and its learning curve, and the RLS adaptive weight coefficient algorithm of the learning process.Keywords:, LMS and RLS adaptive filter, the Matlab simulation课题简介：零均值、单位方差的白噪声通过一个二阶自回归模型产生的AR过程。

AR模型的系统函数为： H(Z)=假设=-1.6,=0.8将系统函数转化为差分方程为：其中w(n)为白噪声，参数=-1.6,=0.8激励源是白噪声w(n)本文用Matlab仿真做出了模型系数的收敛过程与平均的学习曲线分别用LMS算法和RLS算法，分别做出了模型系数的收敛过程与学习曲线，还对两种算法的特性进展了比照引言：由于随机信号的未知性和随时间变化的统计特性，需要设计参数随时间变化的滤波器算法，即所谓的自适应滤波它是利用前一时刻以获得的滤波器参数的结果，自动的调节现时刻的滤波器参数，以适应信号和噪声未知的或随时间变化的统计特性，从而实现最优滤波自适应滤波器的特性变化是由自适应算法通过调整滤波器系数来实现的不同的自适应滤波器算法，具有不同的收敛速度、稳态失调和算法复杂度自适应滤波算法中利用了输出反应，属于闭环算法其优点是能在滤波器输入变化时保持最优的输出，而且还能在某种程度上补偿滤波器元件参数的变化和误差以与运算误差但其缺点是存在稳定性问题以与收敛速度不高所以探讨如何提高收敛速度、增强稳定性以满足信号处理的高效性、实时性，一直是人们研究的重点和热点本文基比照研究了两类根本的自适应算法LMS和RLS，并对它们权系数的收敛过程与学习过程进展了分析。

LMS原理分析：LMS算法是自适应滤波器中常用的一种算法与维纳算法不同的是其系统的系数随输入序列而改变维纳算法中截取输入序列自相关函数的一段构造系统的最优系数而LMS算法如此是对初始化的滤波器系数依据最小均方误差准如此进展不断修正来实现的因此理论上讲LMS算法的性能在同等条件下要优于维纳算法但是LMS算法是在一个初始化值得根底上进展逐步调整得到的因此在系统进入稳定之前有一个调整的时间这个时间受到算法步长因子的控制在一定值X围内增大会减小调整时间但超过这个值X围时系统不再收敛的最大取值为R的迹LMS采用平方误差最小的原如此代替均方误差最小的原如此，信号根本关系如下：写成矩阵形式为：式中W(n)为n时刻自适应滤波器的权矢量，N为自适应滤波器的阶数X(n)为n时刻自适应滤波器的参考输入矢量，由最近的N个信号的采样值构成，d(n)是期望的输出值；e(n)为自适应滤波器的输出误差调节信号；μ是控制自适应速度与稳定性的增益常数LMS的算法流程图：读取x(n)和d(n)初始化w(n)计算误差e(n)=d(n)-y(n) 计算因子更新权RLS算法原理分析：为遗忘因子，它是小于1的正数：参考信号或期望信号第n次迭代的权值均方误差按照如下准如此：越旧的数据对的影响越小。

对滤波器的系数w求偏导，并令结果等于0知整理得到标准方程为：定义:标准方程可以简化为：经求解可以得到迭代形式：定义：，如此可知T的迭代方程为：系数的迭代方程为其中增益和误差的定义分别为：参数递推估计，每取得一次新的观测数据后，就在前次估计结果的根底上，利用新引入的观测数据对前次估计的结果，根据递推算法进展修正，减少估计误差，从而递推地得出新的参数估计值这样，随着新观测数据的逐次引入，一次接一次地进展参数估计，直到参数估计值达到满意的准确程度为止RLS算法流程图：初始化；计算T(n),w(n),k(n),e(n|n-1)计算误差e(n)=d(n)-y(n)更新权LMS算法程序：clearclose allclca1=-1.6; a2=0.8; n=1000; P=50;e=zeros(1,n);ep=zeros(1,n);ee=zeros(1,n);x=zeros(1,n)'; w=randn(1,n)'; %算法for p=1:P x(1)=w(1); x(2)=-a1*x(1)+w(2); for i=3:n x(i)=-a1*x(i-1)-a2*x(i-2)+w(i);end L=2; u=0.0005; wL=zeros(L,n);for i=(L+1):n X=x(i-L:1:(i-1)); y(i)=X'*wL(:,i); %i时刻输出信号 e(i)=x(i)-y(i); %i时刻误差信号 wL(:,(i+1))=wL(:,i)+2*u*e(i)*X; %i时刻滤波器的权值 ee(i)=e(i)^2;end ep=ep+ee; endeq=ep/P;a1L=-wL(2,1:n); % a1在LMS算法下值的变化，wL矩阵中第一行的1到n个数 a2L=-wL(1,1:n); % a2在LMS算法下值的变化 ，wL矩阵中第二行的1到n个数 %画图subplot(2,2,1);plot(1:n,x); title('高斯白噪声w');subplot(2,2,2); plot(1:n,a1L,'r-',1:n,a1,'k-'); title('a1的学习过程');subplot(2,2,3); plot(1:n,a2L,'r-',1:n,a2,'k-'); title('a2的学习过程');subplot(2,2,4);plot(1:n,eq);title('50次平均后的学习过程'); 图1：步长因子μ 图2：步长因子μ图3：步长因子μ结果分析：1. 在仿真过程中可以看到，图形的收敛时间随着步长因子μ的增大而减小。

说明步长因子μ与收敛时间成反比，其决定了LMS算法学习过程的快慢2. 由上图比照可知，当步长因子μ增大时，收敛时间减少，但会使失调增大，当μ等于0.0005与0.001时图形没有失调，但当μ等于0.005时，就会发现图形失调严重3. 综上所述可得出结论：控制失调与加快收敛速度矛盾 LMS与RLS比照程序：程序：clear;close all;clc;a1=-1.6; a2=0.8; n=1000; x=zeros(1,n)'; w=randn(1,n)'; x(1)=w(1); x(2)=-a1*x(1)+w(2); for i=3:n x(i)=-a1*x(i-1)-a2*x(i-2)+w(i); end; %LMS滤波 L=2; u=0.001; wL=zeros(L,n);for i=(L+1):n X=x(i-1:-1:(i-L)); y(i)=X'*wL(:,i); e(i)=x(i)-y(i); wL(:,(i+1))=wL(:,i)+2*u*e(i)*X; end; a1L=-wL(1,1:n); a2L=-wL(2,1:n); %RLS滤波 L=2; namuta=0.98; wR=zeros(L,n);T=eye(L,L)*10;% %RLS算法下T参数的初始化,T初始值为10 for i=(L+1):n X=x(i-1:-1:(i-L)); K=(T*X)/(namuta+X'*T*X);%i时刻增益值 e1=x(i)-wR(:,i-1)'*X; wR(:,i)=wR(:,i-1)+K*e1; %i时刻权值 y(i)=wR(:,i)'*X; e(i)=x(i)-y(i); T=(T-K*X'*T)/namuta; %i时刻的维纳解 end; a1R=-wR(1,1:n); a2R=-wR(2,1:n); %画图subplot(2,1,1); plot(1:n,a1L,'r-',1:n,a1R,'g-',1:n,a1,'k-'); title('LMS与RLS算法a1权系数收敛过程比照'); subplot(2,1,2); plot(1:n,a2L,'r-',1:n,a2R,'g-',1:n,a2,'k-'); title('LMS与RLS算法a2权系数收敛过程比照');图4：LMS与RLS仿真图形比照结果分析：1. RLS算法在算法的稳态阶段即算法的后期收敛阶段其性能和LMS算法相差不明显但在算法的前期收敛段RLS算法的收敛速度要明显高于LMS算法。

但是RLS算法复杂度高计算量比拟大RLS算法与LMS比照：。