第3章节续 快速傅里叶变换(FFT)1.ppt

图4.4.2 分裂基FFT算法L形排列示意图与结构示意图 (a)分裂基FFT算法L形排列示意图； (b)分裂基FFT算法运算流图结构示意图 图4.4.3 16点分裂基第一次分解L形流图(图中省去箭头) 4.4.2 分裂基FFT算法的运算量 设第j级有lj个L形，j=1,2,…,M-1,M=log2N，则有l1=N/4。由图4.4.2(b)可见，第j-1列中的L形包含了第j列中的一部分结点的计算，即空白部分，所占结点数刚好等于第j-1列中所有L形对应结点的一半，所以第j列L形个数就减少l j-1/2个，即 由于每个L形有两次复(数)乘运算，所以全部复 乘次数为 (4.4.5) DIF与DIT比较2 不同点： （1）DIF与DIT两种算法结构倒过来。 DIF为输入顺序，输出倒序。运算完毕再运行“二进制倒读”程序。 DIT为输入倒序，输出顺序。先运行“二进制倒读”程序，再进行求DFT。 （2）DIF与DIT根本区别：在于蝶形结不同。 DIT：W因子在蝶形输入端 DIF：W因子在蝶形输出端 五、变形流图 对于任何流图，只要保持各节点所连接的支路及传输系数不变，则不论节点位置如何排列，所得流图总是等效的，由此可得DIT基2-FFT算法的若干变形流图。 变形流图主要有以下三种： （1）输入顺序，输出倒序 （2）输入输出都顺序 （3）各级的蝶形和蝶距都是相同的 图4.2.14 DIT―FFT的一种变形运算流图 图4.2.15 DIT―FFT的一种变形运算流图 4.2.6 IDFT的高效算法 上述FFT算法流图也可以用于离散傅里叶逆变换(Inverse Discrete Fourier Transform，简称IDFT)。比较DFT和IDFT的运算公式： 图4.2.16 DIT―IFFT运算流图 图4.2.17 DIT―IFFT运算流图(防止溢出) 如果希望直接调用FFT子程序计算IFFT，则可用下面的方法： 由于 对上式两边同时取共轭，得 4.3 进一步减少运算量的措施 回顾 减少运算量的途径 利用DFT运算系数 的周期性和对称性 将长序列分为短序列: 基2时域抽取法(DIT) 基2频域抽取法(DIF) 特点:算法简单,编程效率高 以程序的复杂度来换取进一步减少运算量: 蝶形的旋转因子: 值为 分裂基(混合基) 一、 多类蝶形单元运算 一类蝶形单元运算：包含所有旋转因子 二类蝶形单元运算：去掉 的旋转因子 三类蝶形单元运算：去掉 的旋转因子 四类蝶形单元运算：去掉 的旋转因子 多类蝶形单元  二、减少的计算量: 由DIT―FFT运算流图已得出结论， N=2M点FFT共需要 MN/2 （N/2Log2N）次复数乘法。 ?当L=1时，只有一种旋转因子 ，则第一级不需要乘法运算。 当L=2时，共有两个旋转因子 则第二级不需要乘法运算。 综上所述，先除去第一、二两级后，所需复数乘法次数应是 当L=3时，每级都包含有 同一旋转因子对应 蝶形 从L=3至L=M共减少复数乘法次数为 因此，　DIT―FFT的复乘次数降至 当L=3时，都包含有 一般实现一次复数乘法运算等于 四次实数乘法, 二次实数加法 但对于 任一复数与其相乘 只需两次实数乘法,二次实数加法 从实数运算考虑，计算N=2M点DIT―FFT所需实数乘法次数为 (4.3.4) 从L=3至L=M, 节省的实数法乘同 一样为 三、 旋转因子的生成 在FFT运算中，旋转因子WmN的产生： 每级运算直接产生