f8.m

利用lms算法实现自适应滤波器的matlab仿真

clear
t=0:0.1:999.9;
xs=5*sin(0.05*t);
actual_isi=[0.05 -0.063 0.088 -0.126 -0.25 0.9047 0.25 0.0126 0.038 0.088];
N=length(t);
K=13;
sigma=1;
delta=0.0006
xs1=filter(actual_isi,1,xs);
noise=wgn(1,N,sigma);
y1=xs1+noise;
estimated_c=zeros(1,K);
y_k=zeros(1,K);
y=zeros(1,N);
e_k=zeros(1,N);
for n=K+1:N;
y_k=y1(n-K:n-1);
y(n)=estimated_c*y_k';
e_k(n)=xs(n)-y(n);
estimated_c=estimated_c+2*delta*e_k(n)*y_k;
end
subplot(3,1,1);
plot(t,xs,'r');grid;
ylabel('幅度');
xlabel('时间');
title('输入信号');
subplot(3,1,2);
plot(t,y1,'black');
grid;
ylabel('幅度');
xlabel('时间');
title('加入噪声后的信号');
subplot(3,1,3);
plot(t,y);grid;
ylabel('幅度');
xlabel('时间');
title('维纳滤波后的输出信号');
axis([0 1000 -6 6]);