根據DFT的基二分解方法，可以發現在第L（L表示從左到右的運算級數，L＝1，2，3…M）級中，每個蝶形的兩個輸入數據相距B=2^(L-1)個點，同一旋轉因子對應著間隔為2^L點的2^(M-L)個蝶形。從輸入端開始，逐級進行，共進行M級運算。在進行L級運算時，依次求出個2^(L-1)不同的旋轉因子，每求出一個旋轉因子，就計算完它對應的所有的2^(M-L)個蝶形。因此我們可以用三重循環程序實現FFT變換。同一級中，每個蝶形的兩個輸入數據只對本蝶形有用，而且每個蝶形的輸入、輸出數據節點又同在一條水平線上，所以輸出數據可以立即存入原輸入數據所占用的存儲單元。這種方法可稱為原址計算，可節省大量的存儲單元。附件包含算法流程圖和源程序。

點擊標簽查看更多相關資源：