lms.m

matlab实现陷波器

function [w,y,e]=lms(x,d,u,N)
M=length(x);
y=zeros(1,M);
w=zeros(1,N);
e=zeros(1,M);
for n=N:M
    x1=x(n:-1:n-N+1);
    y(n)=w*x1';
    e(n)=d(n)-y(n);
    w=w+u*e(n)*x1;
end