.DFT的計算工作量
?兩者的差別僅在指數的符號和因子1/N. 通常x(n)和?都是復數,所以計算一個
? X(k)的值需要N次復數乘法運算,和??? 次
? 復數加法運算.那么,所有的X(k)就要N2次復
? 數乘法運算,N(N-1)次復數加法運算.當N很
? 大時,運算量將是驚人的,如N=1024,則要完
? 成1048576 次(一百多萬次)運算.這樣,難以做到實時處理.
二.改進的途徑
? 1.?? 的對稱性和周期性

???? 利用上述特性,可以將有些項合并,并
將DFT分解為短序列,從而降低運算次數,提
高運算速度.1965年,庫利(cooley)和圖基
(Tukey)首先提出FFT算法.對于N點DFT,僅需
(N/2)log2N 次復數乘法運算.例如N=1024=210 時，
需要(1024/2)log2 210 =512*10=5120次。
5120/1048576=4.88% ,速度提高20倍
?按時間抽取(DIT)的FFT算法?????? —庫利-圖基算法
一.算法原理(基2FFT)
(一)N/2點DFT
1.先將??? 按n的奇偶分為兩組作DFT,設N=2L ,不足時,可補些零。這樣有:
? n為偶數時:
? n為奇數時:由于:?
?所以,上式可表示為:
?其中,
2.兩點結論:
?(1) X (k),X (k)均為N/2點的DFT。
?(2) X(k)=X (k)+W? X (k)只能確定出
?? X(k)的k=?????????? 個；
即前一半的結果。
3.X(k)的后一半的確定
?同理,
?這就是說,X1(k),X2(k)的后一半,分別
?等于其前一半的值。

????????? 可見,X(k)的后一半，也完全由X1(k), X2 (k)的前一半所確定。
???? *N點的DFT可由兩個N/2點的DFT來計算。

實現上式運算的流圖稱作蝶形運算

(1)N/2點的DFT運算量:復乘次數:
?????復加次數:
(2)兩個N/2點的DFT運算量:復乘次數:
?????復加次數:
(3)N/2個蝶形運算的運算量:復乘次數:
?????復加次數:

??? 例如 N=8 時的DFT,可以分解為兩個
? N/2=4點的DFT.具體方法如下:
?? (1)n為偶數時,即
? 分別記作: