快速傅里叶和小波变换.ppt

由于 故 因此 X(k)可表示为 2.N点DFT按k的奇偶分组可分为两个N/2的DFT 当k为偶数,即 k=2r时,(-1)k =1； 当k为奇数,即 k=2r+1 时 (-1)k =-1 。 这时 X(k) 可分为两部分： k为偶数时： 可见,上面两式均为N/2的DFT。 k为奇数时： -1 3.蝶形运算 -1 -1 -1 -1 W W W W N N N N 0 1 2 3 4.N=8时,按k的奇偶分解过程 先蝶形运算，后DFT: 5.仿照DIT的方法 再将N/2点DFT按k的奇偶分解为两个N/4点的DFT,如此进行下去,直至分解为2点DFT。 以下是8点的DIF-DFT流程： x(0) X(0) x(1) X(4) x(2) X(2) x(3) X(6) x(4) X(1) x(5) X(5) x(6) X(3) x(7) X(7) -1 -1 -1 -1 W W W W N N N N 0 1 2 3 -1 -1 -1 -1 W W W W N N N N 0 2 0 2 -1 -1 -1 -1 W W W W N N N N 0 0 0 0 B.原位运算 每级(列)都是由N/2个蝶形运算构成,即 -1 W N r C.蝶形运算两节点的距离 一般公式为2L-m =N/2m 例如 N=23 =8 ： (1)m=1 时的距离为 8/2=4； (2)m=2 时的距离为 8/4=2; (3)m=3 时的距离为 8/8=1。 1.相同点 (1)进行原位运算； (2)运算量相同,均为（N/2) Log2N次复乘、N Log2N次复加。 2.不同点 (1)DIT输入为倒位序,输出为自然顺序； DIF正好与此相反。但DIT也有输入为自然顺序,输出为倒位序的情况。 D.DIF法与DIT法的异同 (2)蝶形运算不同 A.DIT 用矩阵表示 = 1 1 B.DIF 用矩阵表示 = 1 1 (3)两种蝶形运算的关系－－互为转置（矩阵）； 将流图的所有支路方向都反向,交换输入和输出，即可得到另一种蝶形。 A. DIT B.DIF 1 1 1 1 FFT算法，同样可以适用于离散博里叶反变换(IDFT)运算，并简称为IFFT，即快速博里叶反变换，从IDFT公式看： 而DFT公式 比较以上两式可知，只要把DFT运算中的每一个系数 换成 ，并且最后再乘以常数1/N, 则 六、IDFT的快速算法--IFFT 以上所有时间抽取或频率抽取的FFT算法都可以拿来运算IDFT。 此外，有一种直接采用FFT程序计算IFFT的方法。 对 取共轭，则： 两边同时再取共轭，有： 上述情况表明：在进行IDFT运算的时候，可以先将X(k)取共轭，直接调用FFT子程序，然后再取共轭并乘以1/N就可得到序列x(