本帖最后由 Happy.Q 于 2013-5-12 09:28 编辑

sbit led0 = P2^0;
sbit led1 = P2^1;
sbit led2 = P2^2;
sbit led3 = P2^3;
sbit led4 = P2^4;
sbit led5 = P2^5;
sbit led6 = P2^6;
sbit led7 = P2^7;
uint8 COUNT=15,COUNT1=0,ADC_Count=0,LINE=15,G=0,T;
uint8 i,j,k,b,p,anum,led=0;                 
int Temp_Real,Temp_Imag,temp;                // 中间临时变量  
uint16 TEMP1,max=0;
int xdata Fft_Real[128];
int xdata Fft_Image[128];               // fft的虚部
uint8 xdata LED_TAB[64];    //记录FFT计算出的各点幅值大小
void FFT()//基2fft
{                 
    for( i=1; i<=7; i++)                            /* for(1) */
    {
        b=1;
        b <<=(i-1);                                 //碟式运算，用于计算隔多少行计算例如 第一极 1和2行计算
        for( j=0; j<=b-1; j++)                      /* for (2) */
        {
            p=1;
            p <<= (7-i);            
            p = p*j;
            for( k=j; k<128; k=k+2*b)               /* for (3)  */
            {
                Temp_Real=Fft_Real[k];
    Temp_Imag=Fft_Image[k];
    temp=Fft_Real[k+b];
                Fft_Real[k]=Fft_Real[k]+((Fft_Real[k+b]*COS_TAB[p])>>7)+((Fft_Image[k+b]*SIN_TAB[p])>>7);
                Fft_Image[k]=Fft_Image[k]-((Fft_Real[k+b]*SIN_TAB[p])>>7)+((Fft_Image[k+b]*COS_TAB[p])>>7);
                Fft_Real[k+b]=Temp_Real-((Fft_Real[k+b]*COS_TAB[p])>>7)-((Fft_Image[k+b]*SIN_TAB[p])>>7);
                Fft_Image[k+b]=Temp_Imag+((temp*SIN_TAB[p])>>7)-((Fft_Image[k+b]*COS_TAB[p])>>7);          // 移位.防止溢出. 结果已经是本值的 1/64               
                Fft_Real[k] >>= 1;            
                Fft_Image[k] >>= 1;
                Fft_Real[k+b] >>= 1;                 
                Fft_Image[k+b] >>= 1;                                                               
            }     
        }
    }