aanlms.m

来自「Adaptive Filtering Primer with MATLAB」· M 代码 · 共 22 行

function[w,y,e,J,w1]=aanlms(x,dn,mubar,M,c)
%all quantities are real-valued;
%x=input data to the filter;dn=desired signal;
%M=filter order;c=constant;mu-bar=step-size equivalent;
%x and dn must be of the same length;J=learning curve;
N=length(x);
y=zeros(1,N);
w=zeros(1,M);%initialized filter coefficient vector;
for n=M:N
   x1=x(n:-1:n-M+1);%for each n vector x1 is of length
   %M with elements from x in reverse order;
   y(n)=w*x1';
   e(n)=dn(n)-y(n);
   w=w+(mubar/(c+x1*x1'))*e(n)*x1;
   w1(n-M+1,:)=w(1,:);
end;
J=e.^2;
%the columns of the matrix w1 depict the history of the
%filter coefficients;