新数字信号处理 教学课件 尹为民 14 4.1改进途经+4.2基2时抽FFT.pptVIP

第四章 快速傅里叶变换(Fast? Fourier Transform) 4.1 直接计算DFT的问题及改进途径 4.2 按时间抽取的基-2FFT算法 4.3 按频率抽取的基-2FFT算法 4.4 N为复合数的FFT算法 4.5 线性调频z变换算法 §4.1 直接计算DFT的问题及改进途径 一、直接计算DFT的运算量 二、改善途径 二、改善途径 § 4.2 按时间抽取的基-2FFT算法 一、算法原理 一、算法原理 一、算法原理 一、算法原理 一、算法原理 二、运算量 二、运算量 三、按时间抽选的FFT算法的特点 三、按时间抽选的FFT算法的特点 四、其他形式流图 由于计算机上乘法运算所需的时间比加法运算所需的时间多得多，故以乘法为例，比较计算量大小。 当点数N越大时，FFT的优点更为明显。 上述按时间抽选的FFT算法的一种有益的附带结果是：该流图提出存储数据以及存储中间各列计算结果的有效方式。 讨论算法特点的目的是寻求FFT流图的一般规律，以便给出流图结构，为FFT软件或硬件实现作准备。 将Xm( k )和Xm( j )分别存放在原存放Xm-1(k)和Xm-1( j )的同一存储寄存器中，则实现全部运算只需要一列存储N个复数的寄存器。 由表4-1可见，自然顺序数I增加1，是在顺序数的二进制数最低位加1，向左进位。而倒序数J则是在二进制数最高位加1， 逢2向右进位。 用这种算法，可以从当前任一倒序值求得下一个倒序值。 流图中特别重要的是连接节点的支路与每个支路的传输比。只要节点间的连接和连接的传输比维持不变，则无论在流图中诸节点如何重新排列，它始终表示相同的计算。只是数据的提取和存放的次序不同而已。这样就可得到按时间抽取的FFT算法的若干其他形式流图。 将图3-5中和x(4)水平相连的所有节点和x(1)水平相连的所有节点位置对调，再将和x(6)水平相连的所有节点与和x(3)水平相连的所有节点对调，其余诸节点保持不变，可得图3-8的流图。 图3-8与图3-5的蝶形相同，运算量也一样，不同点是：① 数据存放的方式不同； ②取用系数的顺序不同。 该流图就是Cooley和Tukey(1965)提出的按时间抽选算法的形式。 只有当重新排列节点使得每个蝶形运算的输入节点和输出节点呈水平相邻时才可以得到对应的同址运算的流图，否则需要两列复数寄存器。 有必要讲一道例题,演示FFT过程. 海军工程大学电气与信息工程学院 ■ ▲ ■ * ■ ▲ ■ ▲ 在线教务辅导网： 教材其余课件及动画素材请查阅在线教务辅导网 QQ:349134187 或者直接输入下面地址： 直接计算DFT的运算量 改善途径 复数加法 N(N-1) 复数乘法 N 2 如何提高DFT的运算效率? 利用旋转因子 的周期性、对称性、可约性。 (1) 周期性 (2) 对称性 (3) 可约性 合并项 长序列分解为短序列 将时域序列逐次分解为一组子序列，利用旋转因子的特性，由子序列的DFT来实现整个序列的DFT。 基2时间抽取(Decimation in time)FFT算法 基2频率抽取(Decimation in frequency)FFT算法 FFT基本思想 算法原理 运算量 按时间抽取的FFT算法的特点 其他形式流图 重点：理解算法原理，掌握8点FFT蝶形运算流图 难点：由短序列的DFT表达相应长序列的DFT的基本原理及方法 分解过程遵循规则：对时间进行奇偶分解 设序列x(n)长度为N，且满足N=2L, L为正整数。按n的奇偶把x(n)分解为两个N/2点的子序列： 对于X(k)后N/2点的值 时域奇偶分解结果：对频域进行前后分解 求一个N点DFT简化为求两个N/2点DFT 时间抽取法蝶形运算单元 按时间( n 的奇偶)抽取的第一次分解 (N=8) X (0) X (1) X (2) X (3) X (4) X (5) X (6) X (7) N点DFT共有2个N/2点DFT， N/2个蝶形运算 复乘数 复加数 运算工作量减少一半！ 分解后运算量： 进一步分解： 按时间( n 的奇偶)抽取的第二次分解 (N=8) -1 -1 -1 -1 按时间( n 的奇偶)抽取的第三次分解 (N=8) 第一级 第二级 第三级 8点基2时间抽取FFT算法流图 当N=2L时，共有L级蝶形， 每级都由N/2个蝶形运算组成，共有蝶形L? N/2个蝶形。 复乘数 复加数 直接计算DFT与FFT算法的复乘计算量之比为： 复乘次数 N N 2 1965，Cooley和Tukey在《计算数学》(Mathematics of Computation)上发表“用机器计算复序列傅里叶级数的一种