串行fft遞歸算法蝶式遞歸計算原理求傅里葉變換-文庫吧資料
2025-07-02 03:16本頁面
【正文】 mp。 /*size/2/l是該級W的相鄰上標差，l是該級該組取的W總個數(shù)*/add(x[j+k],product,amp。k++) //算出第i級內(nèi)j組蝶形單元的結(jié)果{ //算出j組中第k個蝶形單元mul(x[j+k+l],W[(size_x/2/l)*k],amp。j+=2*l){for(k=0。for(j=0。ilog(size_x)/log(2)。change()。*/voidfft(){inti=0,j=0,k=0,l=0。} /*進行基2 FFT運算,蝶形算法。 //利用fft快速算法進行DFT變化output()。initW()。x[i].real,amp。isize_x。printf(Please input the data in x[N]:(such as:5 6)\n)。 /*輸入序列的長度*/scanf(%d,amp。PI=atan(1)*4。 //pi的值int main(){inti。 /*輸出序列的值*/intsize_x=0。 /*復數(shù)除法*/void output()。 /*復數(shù)乘法*/void sub(plex,plex,plex *)。 //碼元位置倒置函數(shù)函數(shù)聲明void add(plex,plex,plex *)。 //快速傅里葉變換函數(shù)聲明void initW()。}plex。 輸出: B=(b0,b1,…,bn1) Procedure RFFT(a,b) begin if n=1 then b0=a0 else (1)RFFT(a0,a2,…,an2, u0,u1,…,un/21) (2)RFFT(a1,a3,…,an1, v0,v1,…,vn/21) (3)z=1 (4)for j=0 to n1 do ()bj=uj mod n/2+zvj mod n/2 ()z=zω endfor endifend 注: (1)算法時間復雜度t(n)=2t(n/2)+O(n) t(n)=O(nlogn) n=8的FFT蝶式計算圖: FFT蝶式計算圖n=6的FFT遞歸計算流程圖: FFT遞歸計算流程圖 開始計算出前size_x/2個exp(j*2π*k/size_x)個值，即W的值輸入序列對應值（例如5+j3，輸入5 3）輸入序列長度size_x飛級數(shù)i=？級數(shù)i加1 是 輸出fft結(jié)果序列結(jié)束 否 該級該組起
外語相關(guān)推薦
鄂ICP備17016276號-1