7-2.m

matlab数字信号与应用源代码7

%例程7－2 LMS算法

Function[h,y]=lms(x,d,u,N)
%LMS算法的实现
%[h,y]=lms(x,d,u,N)
%    h=估计的FIR滤波器
%    y=输出数组y(n)
%    x=输入数组x(n)
%    d=预期数组d(n),其长度应与x相同
%    u=步长
%    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+u*e*x1;
end