lms_demo.m

卡尔曼滤波和自适应滤波 y

M=100;
for i=1:M
    d(i)=sqrt(2.0)*sin(2*pi*i/20.0);
end
d=d+0.5*randn(1,M);
delay=3;
X(1+delay:M)=d(1:M-delay);
n=20;
u=0.0005;

[W,Y]=lms(X,d,n,u);

%   输入信号 
plot(X,':')    
hold on 
% 　参考信号
plot(d,'r')  
%　 实际输出
plot(Y,'m')
clear
% by yyu 2006.3.19