myc1main.c

基于dsp的快速傅立叶变化

 #include "math.h"
 #include "myc1.h"
 #include "stdio.h"
 #pragma DATA_SECTION(x,"FFTx")
 double x[32],y[32];
 double amp[32];
 signed int buffer[32];
 double angle[32];
/*extern void fft(x,y,n);
int n;
double x[],y[];*/


main()
 { int i,j,n;
   /*double a1,a2;*/
   
   
   n=32;
   /*a1=0.9;
   a2=0.3;
   x[0]=1.0;
   y[0]=0.0;
   
   for(i=1;i<n;i++)
   { x[i]=a1*x[i-1]-a2*y[i-1];
     y[i]=a2*x[i-1]+a1*y[i-1];
    }
    
   printf("\n COMPLEX INPUT SEQUENCE\n");
   for(i=0;i<n/2;i++)
    { for(j=0;j<2;j++)
       {printf(" %10.7f+j%10.7f",x[2*i+j],y[2*i+j]);
        printf("\n");
       }
     }
   n=8;
   x[0]=1020.0;
   x[1]=1021.0;
   x[2]=1022.0;
   x[3]=1023.0;
   x[4]=0.0;
   x[5]=0.0;
   x[6]=0.0;
   x[7]=0.0;
   y[0]=0.0;
   y[1]=0.0;
   y[2]=0.0;
   y[3]=0.0;
   y[4]=0.0;
   y[5]=0.0;
   y[6]=0.0;
   y[7]=0.0;*/
   for(i=0;i<n;i++)
   {x[i]=buffer[i];}
   
   
   
   for(i=0;i<n;i++)
   {y[i]=0;}
   
   
   
   
   fft(x,y,n);
   
   /*for(i=0;i<n;i++)
   { amp[i]=sqrt(x[i]*x[i]+y[i]*y[i]);} 
    
    
    for(i=0;i<n;i++)
   { angle[i]=atan(y[i]/x[i]);}*/ 
   
   
    printf("\n DISCRETE FOURISER TRANSFORM\n");
    for(i=0;i<n/2;i++)
    { for(j=0;j<2;j++)
       {printf(" %10.7f+j%10.7f",x[2*i+j],y[2*i+j]);
        printf("\n");
       }
     }
    
    /*printf("\n AMPLITUDE\n"); 
    for(i=0;i<n;i++)
    {printf("%f",amp[i]);
     printf("\n");
    }
    
    printf("\n ANGLE\n"); 
    for(i=0;i<n;i++)
    {printf("%f",angle[i]);
     printf("\n");
    }*/
    
    
    
   }