LMS自适应滤波器举例
N=500;
M=20;
n=1;
a1=-;
h=zeros(M,n+1,3);
e=zeros(M,n,3);
for d=1:3
if d==1 delta=;
else delta=*(d-1);
end;
for k=1:M
b=*randn(1,N);
y(1)=1;
for i=2:N
y(i)=-a1*y(i-1)+b(i);
end
for i=n+1:N
e(k,i,d)=y(i)-h(k,i,d)*y(i-1);
h(k,i+1,d)=h(k,i,d)+delta*y(i-1)*e(k,i,d);
end
end
end
for d=1:3
for i=1:N
em(i,d)=0;
hm(i,d)=0;
for j=1:M
em(i,d)=em(i,d)+e(j,i,d)^2;
hm(i,d)=hm(i,d)+h(j,i,d);
end
end
end
figure(1)
semilogy(1:150,em(1:150,1)),hold on
semilogy(1:150,em(1:150,2),'r'),hold on
semilogy(1:150,em(1:150,3),'g'),hold off
axis([0 150 1]),grid
title('Mean square error ')
xlabel('Samples')
gtext('\leftarrowd=');
gtext('\leftarrowd='
LMS自适应滤波器举例 来自淘豆网www.taodocs.com转载请标明出处.