alp.c

LMS自适应滤波算法,采用LMS算法来完成自适应滤波.

#define twomu 96.0/32768.0
void alp(float *in,float *y,float *e,float *x,float *w,unsigned int Ns,unsigned int N)
{
 unsigned int n;
 int i;
 float temp;
 float uen;
 for(n=0;n<Ns;n++)
 {
  temp=0.0;
  for (i=N-1;i>=0;i--)
   temp+=(w[i]*x[i]);
  y[n]=temp;
  e[n]=in[n]-y[n];
  uen=twomu*e[n];
  for(i=N-1;i>=0;i--)
   w[i]+=uen*x[i];
  for(i=N-1;i>0;i--)
   x[i]=x[i-1];
  x[0]=in[n];
  }
 }