dsplibconvol.c

这里面有很多dsp应用程序源码

#include <math.h>

#define  N   100
unsigned int N1,N2,n;
int i,k;
short sum_real,t_real;
short x_real[N];
short h_real[N];
short y_real[2*N];

extern c54init();

main()
{	
   unsigned int t;
     
   N1=N;                           /* x 长度*/
   N2=N;                           /* h 长度*/
   n=N1+N2-1;                      /* 输出 y 长度*/  
   
   c54init();

   for(i=0;i<N;i++)              /* 初始化数组 */
   {
     x_real[i]=0;     
     h_real[i]=0;
   }
   for(i=0;i<2*N;i++)
     y_real[i]=0;

   for(i=0;i<N;i++)              
   {
      if(i<N/2)      
         x_real[i]=i;
      else     
         x_real[i]=N-i;
      if(i<N/2)
         h_real[i]=i;       
      else     
         h_real[i]=N-i;
     
   }    
  
  /*
  **************************************************************************
  DSPLIB库函数short convol (DATA *x, DATA *y, DATA *r, ushort ny, ushort nr);
  的调用。
  x[nx]       Pointer to real input vector a of  nx real elements
  y[ny]       Pointer to real input vector y of  ny real elements
  r[nx+ny-1]  Pointer to real output vector h of nx+ny-1 real elements
  ny          Number of elements in vector y
  nr          Number of real elements in vector r
  
  这里nx=ny;
  **************************************************************************
  */  
    convol((short *)x_real,(short *)h_real,(short *)y_real,N,2*N-1);
   
   for(;;)
   {
     t++;
   }
}