rls自适应滤波器代码（matlab）

clear all; 
clc; 
m=500;%将输入量赋值 
u=0.002;%将输入量赋值 
a1=-1.6;a2=0.8; 
w1(1)=0;w2(1)=0; 
w1(2)=0;w2(2)=0; 
n(1)=1;n(2)=2; 
%%%%%%%% 
% true sequence x(n) 
rd=randn(1,m); 
x(1)=rd(1);x(2)=rd(2); 
for k=3:m 
    x(k)=rd(k)-a1*x(k-1)-a2*x(k-2); 
end 
R=[x(1)*x(1) 0;0 0]; 
T=0; 
e(1)=0; 
W=[0;1]; 
X=[x(2);x(1)]; 
e(2)=x(2)-W'*X; 
%%%%%%%%%%   
  for  i=3:1:m   %权系数迭代m次    
       R=[x(i-1)*x(i-1) x(i-1)*x(i-2);x(i-1)*x(i-2) x(i-2)*x(i-2)];%列出自相关矩阵 
        T=1/(T+trace(R)); %求出迹的值，为后续u的判断做准备 
        if(u>T) 
            error('u is larger than 1/t[R]');%判断u的值是否小于迹的倒数 
        end 
       W=W+2*u*X*e(i-1);   %LMS算法的权系数迭代公式 
       X=[x(i-1)  x(i-2)]';  %LMS算法中输入信号矢量的递推