convol.c

TI公司TMS320VC5402 DSP 上的卷积（Convolve）算法实验（前面传错了

/********************************************
*					   					   *
* 	 卷积（Convolve）算法实验（C语言版）   *
*					  					   *
********************************************/
#include <math.h>
#define  N   100
unsigned int N1,N2,n;
int i,k;
float sum_real,t_real;
float x_real[N];
float h_real[N];
float 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]=2000.0;
         h_real[i]=2000.0;         
      }
      else
      {
         x_real[i]=0;
         h_real[i]=0;
      }
   }   

   for(i=0;i<n;i++)
   {
      sum_real=0;      
      t_real=0;      

      for(k=0;k<=i;k++)
      {
         t_real=x_real[k]*h_real[i-k];         
         sum_real=sum_real+t_real;
       
      }

      y_real[i]=sum_real;                 
   }
   
   for(;;)
   {
     t++;
   }
}