实数FFT算法的设计及其C语言实现

　　dataR[k+b]=TR-dataR[k+b]*cos_tab[p]-dataI[k+b]*sin_tab[p];

　　dataI[k+b]=TI+temp*sin_tab[p]-dataI[k+b]*cos_tab[p];

　　} /* END for (3) */

　　} /* END for (2) */

　　} /* END for (1) */

　　for(i=0;i32;i++){ /* 只需要32次以下的谐波进行分析 */

　　w[i]=sqrt(dataR[i]*dataR[i]+dataI[i]*dataI[i]);

　　w[i]=w[i]/64;}

　　w[0]=w[0]/2;

　　} /* END FFT */