FFT算法C语言程序代码.doc

《FFT算法C语言程序代码.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、DIT-基2FFT的浮点C语言程序：1、生成旋转因子，复数结构,旋转因子Wn=exp(-j*2*pi/N)//twiFactor——指向旋转因子矩阵的指针//wLen——FFT的长度StructcomplexData{//定义一个复数结构floatre;floatim;};Voidgen_w_r2(structcomplexData*twiFactor,intwLen){intiFactor;floatstepFactor;stepFactor=2.0*pi/wLen;for(iFactor=0;iFactor<(

2、wLen>>1);iFactor++){twiFactor[iFactor].re=cos(stepFactor*iFactor);twiFactor[iFactor].im=sin(stepFactor*iFactor);//W[n]=exp(j*2*pi*n/N),n=0,1,…,(N/2-1)}}2、在运行FFT之前，对输入序列进行倒序变换，代码如下：//bitRevData——指向位变换序列的指针//revLen——FFT长度Voidbit_rev(structcomplexData*bitRevData,

3、intrevLen){structcomplexDatatempRev;intiRev,jRev,kRev,halfLen;halfLen=revLen>>1;jRev=0;for(iRev=0;iRev<(revLen-1);iRev++){If(iRev

4、v-kRev;kRev=kRev>>1;}}}3、FFT计算。有3个循环体，分别为a内循环,b中间循环,c外循环,内循环实现蝶形结计算，循环a和b完成所有的蝶形结运算，而循环c则表示完成FFT算法所需要的级数。//x——输入数据的指针//w——旋转因子指针//n——FFT的长度voidsp_cfftr2_dit(float*x,float*w,shortn){shortn2,ie,ia,i,j,k,m;floatrtemp,itemp,c,s;n2=n;ie=1;for(k=n;k>1;k>>=1)//loopc{

5、n2>>=1;ia=0;for(j=0;j