limlms.m

自适应滤波器程序

%LMS算法，选用电力系统谐波信号作为输入信号。参考书目：现代数字信号处理，皇甫堪等，电子工业出版社。99页和107页
clear all;close all;
f1=50;ts=1/1000;n=1:256;N=256;M=10;
d=2*sin(2*pi*f1*n*ts)+1.4*sin(2*pi*f1*7*n*ts);
x=sin(2*pi*f1*n*ts+pi/2);w=zeros(1,M);w=w';u=0.005;
for n=1:N
    for m=1:M
        if n<=m-1
            x1(m)=0;
        else
            x1(m)=x(n-m+1);
        end
    end
        X(:,n)=x1';
        y(n)=w(:,n)'*X(:,n);
        e(n)=d(n)-y(n);
        w(:,n+1)=w(:,n)+2*u*e(n)*X(:,n);
end
    plot(y)