lms.m

Matlab 6.0图形图像处理所附带的光盘

function [h,y] = lms(x,d,delta,N)

% LMS Algorithm for Coefficient Adjustment

% ----------------------------------------

% [h,y] = lms(x,d,delta,N)

%     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