快速傅里叶算法数字信号处理重点分析.ppt

第一级 第二级 第三级 000 001 010 011 100 101 110 111 000 100 010 110 001 101 011 111 0 0 0 0 0 0 0 0 自然顺序 二进制k2k1k0 倒位序二进制k0k1k2 倒位顺序 1 0 0 1 1 0 0 4 2 0 1 0 0 1 0 2 3 0 1 1 1 1 0 6 4 1 0 0 0 0 1 1 5 1 0 1 1 0 1 5 6 1 1 0 0 1 1 3 7 1 1 1 1 1 1 7 倒位序规律 由按时间抽选法FFT运算流图(书140页图3-6)可知，输出X[m]按正常顺序排列在存储单元，而输入x[k]是按顺序: 0 1 0 1 0 1 0 1 以N=8为例， 说明如下: k0=1 (奇) k0=0 (偶) k1=0 k1=1 k1=0 k1=1 k0 k1 k2 x[000] 0 x[100] 4 x[010] 2 x[110] 6 x[001] 1 x[101] 5 x[011] 3 x[111] 7 序号k x[0] x[1] x[2] x[3] x[4] x[5] x[6] x[7] x[0] x[4] x[2] x[6] x[1] x[5] x[3] x[7] 变址处理方法 存储单元 自然顺序 变址 倒位序 A(1) A(2) A(3) A(4) A(5) A(6) A(7) A(8) 不必调换数值 交换x[k]与x[ ]之间数值 1.原位计算 n表示第n级迭代，m，j表示数据所在的行数 输入数据、中间运算结果和最后输出均用同一存储器 Xn-1[m] Xn-1[j] 第二级的蝶形系数为 ，蝶形节点的距离为2。 第一级的蝶形系数均为 ，蝶形节点的距离为1。 第三级的蝶形系数为 ，蝶形节点的距离为4。 第M级的蝶形系数为 ，蝶形节点的距离为N /2。 级数 碟形运算的数量 序列倒序 序列原位 旋转因子分布 例：试利用N=4基2时间抽取的FFT流图计算8点序列x[k]={1, -1, 1, -1, 2, 1, 1, 2}的DFT。(书上作业3-5) 解： 根据基2时间抽取FFT算法原理，8点序列的DFT X[m]可由两个4点序列的DFT X1[m]和X2[m]表达。如果按照序列x[k]序号的奇偶分解为x1[k]和 x2[k]，则存在 其中x1[k]={1, 1, 2, 1}，x2[k]={-1, -1, 1, 2}，X1[m]和X2[m]可通过4点的FFT来计算。 例：试利用N=4基2时间抽取的FFT流图计算8点序列x[k]={1, -1, 1, -1, 2, 1, 1, 2}的DFT。 解： X1[m]={5, -1, 1, -1}, X2[m]={1, -2+3j, 1, -2-3j} 利用上述公式，可得序列x[k]的DFT X[m]为 X[m]={6, -0.293+3.535j, 1+j, -1.707 + 3.535j, 4, -1.707-3.535j, 1-j, -0.293-3.535j} 快速傅里叶变换(FFT) 基2时间抽取FFT算法 与DFT的运算量对比 算法的特点 课本: P94 ~108 参考书: P94~108 练习: P114 3-1, 3-2 * * * * * * * 离散傅里叶变换快速算法(FFT) 快速傅里叶变换(FFT) 问题的提出 基2时间抽取FFT算法 基2频率抽取FFT算法 IFFT算法的实际应用 掌握基2时间抽取FFT算法的基本思想和方法 掌握基2频率抽取FFT算法的基本思想和方法 掌握实序列F