致各位高手--求RLS算法（matlab）
各位前辈们好，本人现在正在研究自适应滤波器的除噪效果由于比较菜,在做RLS自适应滤波器时遇到一些困难，现想求一个RLS自适应滤波器的matlab实现程序
以下是基于LMS算法的程序代码
clear all;
t=0:0.0001:1.9999;
N=20000;
%生成有用信号s
s=3*sin(2*pi*100*t);
%figure(1);
%subplot(3,1,1);
%plot(t,real(s));
%title('有用信号s');
%xlabel('n');
%ylabel('s');
%生成两个相关的噪声n1,n2
n1=wgn(1,N,3);
n2=zeros(1,N);
n2(1:18000)=n1(1:18000);
n3=wgn(1,N,3);
n2(18001:N)=n3(18001:N);
%figure(1);
%subplot(3,1,2);
%plot(t,real(n2));
%title('噪声参考输入');
%xlabel('n');
%ylabel('n2');
%输入端期望信号d=s+n1
d=s+n1;
subplot(3,1,3);
plot(t,real(d));
title('含有噪声的信号d');
xlabel('n');
ylabel('d');
%初始化各参量
M=8;%滤波器阶数
w=zeros(M,1); %滤波器权系数
mu=0.01;%更新步长
iter=2000; %迭代次数
x=zeros(M,1); %输入矢量
%LMS循环算法
for i=1:iter
i
x(1)=n2(i);
y(i)=dot(w,x);
e=d(i)-y(i);
w=w+mu*e*x;
for j=M:-1:2
x(j)=x(j-1);
end
error(i)=e; %保存输出误差
cost(i)=e*e;
end
figure(2);
subplot(1,1,1);
plot(error);
title('除噪后的信号');
xlabel('n');
ylabel('error');
求助各位高手把核心的lms循环算法改成RLS的算法，小弟将不胜感激
页:
[1]