基于MSP430的FFT算法源码

之前由于工作需要，在MSP430F169和IAR EW430v3.42环境下做了整型数的FFT算法，主要是为了追求快速计算，但实际实验中发现定点数算法在小信号的时候精度损失比较大，所以最近又重新作了单精度浮点数的FFT算法，同样只针对实数序列。
    经实验发现，MSP430的浮点数计算并没有想象中的那么慢，完成一次128点的FFT序列变换，大概需要130ms左右（8MHz主频），下面是我的源代码：

__no_init float fft_r[FFT_N]; //FFT输入实数序列，保存变换后的频域实部
__no_init float fft_i[FFT_N]; //保存变换后的频域虚部
__no_init float harm[HARM_M+1]; //各次谐波幅值
const Uint8 FFT_BIT[8]={ //供完成倒位序排列使用的位定义
0x01,0x02,0x04,0x08,0x10,0x20,0x40,0x80};

void fft(void)
{
    Uint16 i,j,k,L;
    Uint16 b,c,p,n;
    Uint8 x,xc;
    float Gr,Gi,Wr,Wi;

   
    for(i=0;i<FFT_N;i++)
    {
        x =i;
        xc=0;
        for(j=0;j<FFT_M;j++)
        {
            if((x&0x01)!=0)
            {
                xc|=FFT_BIT[(FFT_M-1)-j];
            }
            x>>=1;
        }
        fft_i[xc]=fft_r[i];