/*****************************************************************************/
/* 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);
   //求頻譜
      for(i=0;i<NX;i=i+2)
    {
     p=y[i];
     q=y[i+1];
     temp=p*p+q*q;
     f[m]=sqrt(temp);
     m++;
     }
    
   //反變換    
    unpacki(y,NX);
     cbrev(y,y,NX/2);
     rifft(y,NX,1);
     
   } 
       
     return;
}