lms.m

来自「It is a GUI designed for the learning of」· M 代码 · 共 28 行

function [e,W]=lms(mu,N,u,d);


% Input arguments:
% mu = step size, dim 1x1
% M = filter length, dim 1x1
% u = input signal, dim Nx1
% d = desired signal, dim Nx1
%
% Output arguments:
% e = estimation error, dim Nx1
% w = final filter coefficients, dim Mx1

w=zeros(N,1);
M=length(u);
u=u(:);
d=d(:);
W=[];
%LMS starts here
for n=N:M
    uvec=u(n:-1:n-N+1);
    e(n)=d(n)-w'*uvec;
    w=w+2*mu*uvec*conj(e(n));
    W=[W w];
end
e=cumsum((e(:)).^2)./cumsum(ones(length(e),1));
% e=(e(:)).^2;