lms.m

自适应中关于格型算法和LMS的仿真程序

clc
clear all
u=0.0004;
N=10000;
x(1)=0;
x(2)=0;
v=randn(1,N); 
for n=3:N
    x(n)=1.2728*x(n-1)-0.81*x(n-2)+v(n);
end
w1(2)=0; 
w2(2)=0;
w1(3)=0;
w2(3)=0;
for n=3:N-1
    y(n)=w1(n)*x(n-1)+w2(n)*x(n-2);
    e(n)=x(n)-y(n);
    w1(n+1)=w1(n)+2*u*e(n)*x(n-1);
    w2(n+1)=w2(n)+2*u*e(n)*x(n-2);
end
subplot(2,1,1);
plot(w1);
title('w1(n)-n');
subplot(2,1,2);
plot(w2);
title('w2(n)-n');