根据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变换。同一级中，每个蝶形的两个输入数据只对本蝶形有用，而且每个蝶形的输入、输出数据节点又同在一条水平线上，所以输出数据可以立即存入原输入数据所占用的存储单元。这种方法可称为原址计算，可节省大量的存储单元。附件包含算法流程图和源程序。