fft_test.c

DSP课程课堂作业程序

/*****************************************************************************/
/* fft-test.c																	 */ 
/* 																			 */	
/* FFT example													 */
/*																			 */
/*****************************************************************************/

#include <type.h>
#include <math.h>
#include <tms320.h>
#include <dsplib.h>

/*****************************************************************************/
/* Function Prototypes														 */
/*****************************************************************************/


#define NX 64
#pragma DATA_SECTION(x,"data_buf")

DATA x[NX]={ 
       13,    -32,    -31,     22,     52,    -16,    -84,     -9,    117,
      55,   -142,   -128,    147,    225,   -115,   -339,     30,    455,
     123,   -551,   -359,    597,    691,   -549,  -1143,    341,   1775,
      176,  -2825,  -1578,   5900,  13543,  13543,   5900,  -1578,  -2825,
     176,   1775,    341,  -1143,   -549,    691,    597,   -359,   -551,
     123,    455,     30,   -339,   -115,    225,    147,   -128,   -142,
       55,    117,     -9,    -84,    -16,     52,     22,    -31,    -32,
      13

       };
DATA y[NX];
double p,q,temp;



u16 f[NX/2+1];   //取放频谱
/*****************************************************************************/
/* MAIN																		 */
/*****************************************************************************/

void main()
{   int i,m;
    
 // while(1)
  //{ 
     
  //正变换
    cbrev(x,y,NX/2);
    rfft(y,NX,0);
   //求频谱
   f[0]=abs(y[0]);
   f[NX/2]=abs(y[1]);
   m=1;
      for(i=2;i<NX;i=i+2)
    {
     p=y[i];
     q=y[i+1];
     temp=sqrt(p*p+q*q);
     f[m]=temp;
     m++;
     }
     
   
     //反变换 ,如果进行反变换，打开以下3行   
     // unpacki(y,NX);
     //cbrev(y,y,NX/2);
     //rifft(y,NX,1);
     
   } 
       
    // return;*/
//}