离散傅里叶变换的计算精要.pptVIP

可见,上面两式均为N/2的DFT。 k为奇数时： 3.蝶形运算 -1 4.N=8时,按k的奇偶分解过程 先蝶形运算，后DFT: -1 -1 -1 -1 W W W W N N N N 0 1 2 3 5.仿照DIT的方法 再将N/2点DFT按k的奇偶分解为两个 N/4点的DFT,如此进行下去,直至分解为 2点DFT。 (0) X(0) (1) X(4) (2) X(2) (3) X(6) (4) X(1) (5) X(5) (6) X(3) (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 例如 N=8时DIF的FFT流图如下： 二.原位运算 每级(列)都是由N/2个蝶形运算构成,即 -1 W N r 三.蝶形运算两节点的距离 一般公式为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。 四. 的计算 由于DIF蝶形运算的两节点的距离为 N/2m ,所以蝶形运算可表为： r的求法： k=(n2 n1 n0 ) ,左移m-1位,右边空出补零, 得(r)2 ,亦即(r)2 =(k)2 2m-1 . 例如,N=8: (1)m=1,k=2, k=(010)2 左移0位,(r)2=(010)2=2; (2)m=2,k=1, k=(001)2 左移1位,(r)2=(010)2=2; (3)m=2,k=5, k=(101)2 左移1位,(r)2=(010)2=2 . 1.相同点 (1)进行原位运算； (2)运算量相同,均为（N/2) Log2N 次复乘,N Log2N次复加。 2.不同点 (1)DIT输入为倒位序,输出为自然顺序；DIF正好与此相反。但DIT也有输入为自然顺序,输出为倒位序的情况。 五.DIF法与DIT法的异同 (2)蝶形运算不同 a.(DIT) 用矩阵表示 = 1 1 b.(DFT) 用矩阵表示 = 1 1 (3)两种蝶形运算的关系 互为转置（矩阵）； 将流图的所有支路方向都反向,交换 输入和输出，即可得到另一种蝶形。 a.(DIT) b.(DFT) 1 1 1 1 IFFT算法 一.稍微变动FFT程序和参数可实现IFFT 另外,可以将常数1/N分配到每级运算中, ∵ ,也就是每级蝶形运算均乘 以1/2。 比较两式可知,只要DFT的每个系数 换 成 ,最后再乘以常数1/N就可以得到IDFT 的快速算法--IFFT。 二.不改(FFT)的程序直接实现IFFT 这就是说,先将X(k)取共轭,即将X(k)的 虚部乘-1，直接利用FFT程序计算DFT；然后 再取一次共轭；最后再乘1/N,即得x(n)。所 以FFT,IFFT可用一个子程序。 FFT FFT IFFT x x(n) h(n) y(n) X(k) H(k) Y(k) 流程图 二.用FFT算y(n)的步骤： 三.几点说明 1. 当 M=L 时,用圆周卷积计算线性 卷积的速度快,点数越多,速度越快, N≥64时,速度增加明显. 2. LM 时,速度增加不太明显,为了 增加速度,采用 (1)重叠相加法 (2) 重叠保留法(从略). 9.6 用卷积实现DFT 有些情况下卷积比FFT的计算量小 9.6.1 Winograd 傅里叶变换 其计算量：乘法的计算量正比于N。但加法的计算量多，因此适合乘法