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);

}