test16.c

来自「这里面有很多dsp应用程序源码」· C语言 代码 · 共 49 行

#include <math.h>

#define pi 3.1415926  

int N,order;

float x[500],w[50],d[500],e[500],y[500];

main()
{ 
   int i,j;
   float sum; 
   float step=0.0000001;
   
   N=500;
   order=3;                                                                        
   

   for(i=0;i<N;i++)
   { 
      x[i]=0;
      d[i]=0;
      e[i]=0;
      y[i]=0;
   }                                                                   

   for(i=0;i<N;i++)
   { 
      x[i]=(float)100*sin(pi*i/20);
      d[i]=x[i-2];
   }

   for(i=0;i<order;i++)
      w[i]=0;

   for(i=order;i<N;i++)
   {        
      sum=0;
      for(j=0;j<order;j++)
         sum=sum+x[i-j]*w[j];
      y[i]=sum;
      e[i]=d[i]-y[i];
      for(j=0;j<order;j++)
         w[j]=w[j]+2*step*e[i]*x[i-j];
   }   
}