lms.m

数字信号处理（机械工业出版社）的源码

function [h,y] = lms(x,d,delta,N)
% 	Least-Mean-Square(LMS) Algorithm for Coefficient Adjustment
%     h = estimated FIR filter
%     y = output array y(n)
%     x = input array x(n)
%     d = desired array d(n), length must be same as x
% 	  delta = step size
%     N = length of the FIR filter
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