//********************************************************** // 函數名: 快速傅立葉變換（來源《C常用算法集》） // 本函數測試OK,可以在TC2.0,VC++6.0,Keil C51測試通過。 // 如果你的MCS51系統有足夠的RAM時,可以驗證一下用單片機處理FFT有多么的慢。 // // 入口參數： // l: l = 0, 傅立葉變換; l = 1, 逆傅立葉變換 // il: il = 0,不計算傅立葉變換或逆變換模和幅角；il = 1,計算模和幅角 // n: 輸入的點數，為偶數，一般為32，64，128，...,1024等 // k: 滿足n=2^k(k>0),實質上k是n個采樣數據可以分解為偶次冪和奇次冪的次數 // pr[]: l=0時，存放N點采樣數據的實部 // l=1時, 存放傅立葉變換的N個實部 // pi[]: l=0時，存放N點采樣數據的虛部 // l=1時, 存放傅立葉變換的N個虛部 // // 出口參數： // fr[]: l=0, 返回傅立葉變換的實部 // l=1, 返回逆傅立葉變換的實部 // fi[]: l=0, 返回傅立葉變換的虛部 // l=1, 返回逆傅立葉變換的虛部 // pr[]: il = 1,i = 0 時，返回傅立葉變換的模 // il = 1,i = 1 時，返回逆傅立葉變換的模 // pi[]: il = 1,i = 0 時，返回傅立葉變換的輻角 // il = 1,i = 1 時，返回逆傅立葉變換的輻角 // data: 2005.8.15,Mend Xin Dong #include <math.h>void kkfft(double pr[], double pi[], int n, int k, double fr[], double fi[], int l, int il) {int it,m,is,i,j,nv,l0;double p,q,s,vr,vi,poddr,poddi;for (it=0; it<=n-1; it++){ m = it; is = 0;for(i=0; i<=k-1; i++){ j = m/2; is = 2*is+(m-2*j); m = j;}fr[it] = pr[is]; fi[it] = pi[is];}//----------------------------pr[0] = 1.0; pi[0] = 0.0;p = 6.283185306/(1.0*n);pr[1] = cos(p); pi[1] = -sin(p);if (l!=0) pi[1]=-pi[1];for (i=2; i<=n-1; i++){ p = pr[i-1]*pr[1]; q = pi[i-1]*pi[1];s = (pr[i-1]+pi[i-1])*(pr[1]+pi[1]);pr[i] = p-q; pi[i] = s-p-q;}for (it=0; it<=n-2; it=it+2){ vr = fr[it]; vi = fi[it];fr[it] = vr+fr[it+1]; fi[it] = vi+fi[it+1];fr[it+1] = vr-fr[it+1]; fi[it+1] = vi-fi[it+1];}m = n/2; nv = 2;for (l0=k-2; l0>=0; l0--){ m = m/2; nv = 2*nv;for(it=0; it<=(m-1)*nv; it=it+nv)for (j=0; j<=(nv/2)-1; j++){ p = pr[m*j]*fr[it+j+nv/2];q = pi[m*j]*fi[it+j+nv/2];s = pr[m*j]+pi[m*j];s = s*(fr[it+j+nv/2]+fi[it+j+nv/2]);poddr = p-q; poddi = s-p-q;fr[it+j+nv/2] = fr[it+j]-poddr;fi[it+j+nv/2] = fi[it+j]-poddi;fr[it+j] = fr[it+j]+poddr;fi[it+j] = fi[it+j]+poddi;}}if(l!=0)for(i=0; i<=n-1; i++){ fr[i] = fr[i]/(1.0*n);fi[i] = fi[i]/(1.0*n);}if(il!=0)for(i=0; i<=n-1; i++){ pr[i] = sqrt(fr[i]*fr[i]+fi[i]*fi[i]);if(fabs(fr[i])<0.000001*fabs(fi[i])){ if ((fi[i]*fr[i])>0) pi[i] = 90.0;else pi[i] = -90.0;}elsepi[i] = atan(fi[i]/fr[i])*360.0/6.283185306;}return; }
與50位技術專家面對面20年技術見證，附贈技術全景圖

總結

以上是生活随笔為你收集整理的fft快速傅利叶变的C实现的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。