2007dsp第十三讲.ppt.pptVIP

第十三讲 4.2.5 频域抽取法FFT(DIF-FFT) 4.2.5 频域抽取法FFT(DIF-FFT)(桑德—图基算法) 1、算法原理 N /2仍为偶数，进一步分解：N /2 N /4 逐级分解，直到2点DFT 2、算法特点 1）原位计算 2）蝶形运算 对N=2M点FFT，输入自然序，输出倒位序时 两节点距离：2M-L=N / 2L 3、DIT与DIF的异同 基本蝶形不同 离散傅里叶变换对为： 式中 。 下面要用矩阵来表示DFT关系。 如果用i (0≤i≤N-1)表示这两个N阶方阵的行号，用j(0≤j≤N-1)表示这两个N阶方阵的列号，那末很容易看出，WN 方阵中i行j列的元素为 ，而方阵中i行j列的元素为 。 于是DFT可以写成： 而IDFT式可以写成：  由DIT―FFT运算流图已得出结论，N=2M点FFT共需要MN/2次复数乘法。 当L=1时，只有一种旋转因子 ，所以，第一级不需要乘法运算。 综上所述，先除去第一、二两级后，所需复数乘法次数应是 从L=3至L=M共减少复数乘法次数为 从实数运算考虑，计算N=2M点DIT―FFT所需实数乘法次数为 4.3.2 旋转因子的生成 （略) 4.3.3 实序列的FFT算法 设x(n)为N点实序列，取x(n)的偶数点和奇数点分别作为新构造序列y(n)的实部和虚部，即 由于x(n)为实序列，所以X(k)具有共轭对称性，X(k)的另外N/2点的值为 2.特殊的复数运算 任何复数(x+jy)与其相乘有如下规律 只需两次实加和两次实乘，节约两次实乘。 在DIT-FFT运算流图中，从L=3到L=M级，每一级都包含旋转因子 ，第L级中， 对应 个碟形运算。因此从第三级至最后一级，旋转因子 节省的实数乘次数相同 。 对y(n)进行N/2点FFT，输出Y(k)，则 根据DIT―FFT的思想及式(4.2.7)和(4.2.8)，可得到 * 设序列点数N=2M，M为整数。 将X(k)按k的奇偶分组前，先将输入x(n)按n的顺序分成前后两半： 该算法的思想是把X(k)分成N/2、N/4……2的子序列，以这种方式为基础的FFT算法称为频率抽取。 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。 x3(0) x3(1) -1 -1 x4(0) x4(1) N/4点 DFT N/4点 DFT x1(0) x1(1) x1(2) x1(3) X3(0)=X1(0)=X(0) X4(0)=X1(1)=X(2) X3(1)=X1(2)=X(4) X4(1)=X1(3)=X(6) 同理： 其中： 当N=8时，即分解到x3(n)，x4(n)，x5(n)，x6(n)，n=0,1 (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 -