fft_test.c

一个基于dsp的fft实验程序

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

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

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


#define NX 64

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;
double f[32];   //取放频谱

/*****************************************************************************/
/* MAIN																		 */
/*****************************************************************************/

void main()
{   int i,m;
    
  while(1)
  { 
   m=0;
  //正变换
    cbrev(x,y,NX/2);
    rfft(y,64,0);
   //求频谱
      
     
   } 
       
     return;
}