lms.m

运用MATLAB设计一个LMS自适应滤波器

% w:滤波器加权系数， y:输出信号， e:误差， x:输入信号， d:期望输入， u:步长， N:权矢量个数

function [w,y]=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