lms.m

自适应算法 主要用于自适应天线的算法防真程序

function [e,h]=lms(y,fl,mu);
% [e,h]=lms(y,fl,mu)
%	y = input signal
%	fl = filter length
%	mu = convergence factor (0<mu<1)
%	e = error signal
%	h = filter coefficients
%
% ref _Adaptive_Filters_ by Cowan and Grant, p44

% cjm 15mar95
% Chris McCormack / Radiation Laboratory / University of Michigan
%                                            [e-mail]

if mu==0; mu=.001; end
if fl<1;
   e='lms : Illegal filter length'
  else
   sl=size(y);
   if sl(1)<sl(2); y=y'; sl=sl(2); end
   sl=sl(1);			% set up length of signal
   e=zeros(sl,1);		% initialize
   h=zeros(fl,1);
   phi=zeros(fl,1);
   for n=1:sl;
      e(n)=y(n)-phi'*h;
      h=h+2*mu*e(n)*phi;
      phi=[y(n); phi(1:fl-1)];
   end
end