function [h,y] = lms(x,d,delta,N)
%系數調整的LMS算法
% ----------------------------------------
%  [h,y] = lms(x,d,delta,N)
%     h = 估計的 FIR 濾波器
%     y = 輸出數組 y(n)
%     x = 輸入數組 x(n)
%     d = 預期數組 d(n), 其長度應與 x 相同
%  delta = 步長
%     N =  FIR濾波器的長度
%
M = length(x); y = zeros(1,M);
h = zeros(1,N);
for n = N:M
    x1 = x(n:-1:n-N+1);
     y = h * x1';
     e = d(n) - y;
     h = h + delta*e*x1;
end