%%(2)RLS 滤波
W=[0;0]; %设置权系数初始值
X=[x(2);x(1)]; %设置信号初始值
R=X*X'; %计算 x(n)相关矩阵
T=inv(R); %计算 T 初始值
P=[0;0];
for i=L+1:m
X=[x(i-1),x(i-2)]';
R=lamuda*R+ X*X';
e=x(i)-W'*X;
W=W+inv(R)*X*e;
RLS_a1(i)=W(1); %RLS 算法中权系数 a1 的提取
RLS_a2(i)=W(2); %RLS 算法中权系数 a2 的提取
end
figure(1)
plot(x,'r-');
title('x(n)曲线绘图');xlabel('迭代次数 n');ylabel('x(n)');grid on;
figure(2)
plot(LMS_a1,'r'); hold on;
plot(RLS_a1,'g'); hold on
plot(LMS_a2,'b'); hold on;
plot(RLS_a2,'m');
legend('LMS(a1)','RLS(a1)','LMS(a2)','RLS(a2)');
title('自适应权系数的过渡过程(RLS 和 LMS 算法比较)'); xlabel('迭代次数 m'); ylabel('权系
数'); grid on;