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