function [W,Y]=lms(X,d,n,u)
% X: input data
% d: reference value
% n: adaptive filter order
% u: factor
% Wi: weight
% Y: actual output

format long;

m=length(X);
W=zeros(n,m);
Y=zeros(1,m);

for j=1:m
    for i=1:n
        if(j-i>0) 
            Y(j)=Y(j)+W(i,j)*X(j-i);
        end 
    end
    e(j)=d(j)-Y(j);
    for i=1:n
        if(j-i>0)
            W(i,j+1)=W(i,j)+2*u*e(j)*X(j-i);
        end
    end
end
% by yyu 2006.3.19