lms.m

正弦干扰陷波LMS算法

function [w,y,EE,e,r]=LMS(x,d,u,N)
%LMS算法

M=length(x);
y=zeros(1,M);
w=zeros(1,N);

%计算权矢量w
r=0;
e=x-y;
for n=N:M
    x1=x(n:-1:n-N+1);
    y(n)=w*x1';%x'是x的转置
    e(n)=d(n)-y(n);
    w=w+u*e(n)*x1;
    EE(r+1)=e*e'/(r+1);  %EE是均方误差
    r=r+1;
end