filter_nlms.c

归一化最小均方（LMS）自适应数字滤波

#include "math.h"
void nlms(x,d,n,w,m,mu,sigma2,a,px)
int m,n;
double a,mu,sigma2,d[],x[],w[],px[];
{int i,k;
double e,tmp;
for(k=0;k<n;k++)
{px[0]=x[k];
x[k]=0.0;
for(i=0;i<m;i++)
{x[k]+=px[i]*w[i];}
e=d[k]-x[k];
sigma2=a*px[0]*px[0]+(1.0-a)*sigma2;
tmp=2*mu/(m*sigma2);
for(i=0;i<m;i++)
{w[i]+=tmp*e*px[i];}
for(i=(m-1);i>=1;i--)
{px[i]=px[i-1];}
}
}