MSP430 FFT算法

张衍波

/*****************main   programe********************/
#include  <msp430x14x.h>

#include       <math.h>     
#include       <stdio.h>     
#include       <stdlib.h>     
#include      "typedef.h "
      

float       result[8];   
struct     compx   s[8];      
int       Num=8;
const   float   pp=3.141592653589793;
extern void FFT(struct   compx *s,int num);

double mpr[8],mpi[8],mfr[8],mfi[8];
int n,k,l,il;
extern void kbfft();
extern void FFT_dingdian();
main()
{
    //   ????
    int   i;
    FILE   *fin,*fout;
    if((fin=fopen( "indata.txt ", "r "))==NULL)
    {  
         printf( "can 't   open   infile ");
        exit(0);
    };
    for(i=0;i <Num;i++)
    {   
         //fscanf(fin, "%lf ",&s.real);
        s.real=cos(2*3.1415926*(30)*(i)/80.0);//
         s.imag=(double)0;
        mpr=s.real;
         mpi=0;        
         
     };
    fclose(fin);
    //////////////
      
     //FFT(s,Num);
   
    //n=8;k=3;l=0;il=0;
     //kbfft(mpr,mpi,n,k,mfr,mfi,l,il);

   
    FFT_dingdian();
      
     //????
    if((fout=fopen( "fftresult.txt ", "w "))==NULL)
    {
         printf( "can 't   open   outfile ");
        exit(0);
    };
   
    for(i=0;i <Num;i++)
    {
        // result=sqrt(pow(s.real,2)+pow(s.imag,2));
        fprintf(fout, "%lf+(%lf)i   ",s.real,s.imag);
    };
   
    fclose(fout);
}