数字信号处理Digital Signal ProcessingLecture 13 快速傅立叶变换.ppt

4）需要的存储单元 输入序列 x(n)：N 个存储单元 系数 ： N / 2 个存储单元 4、DIT算法的其他形式流图 输入倒位序输出自然序 输入自然序输出倒位序 输入输出均自然序 相同几何形状 四 、IFFT算法 比较： IDFT: DFT: 利用FFT程序计算IFFT，需做一些小的改动 可以不做任何改动，直接调用FFT程序就实现IFFT吗？ 共轭 FFT 共轭 乘1/ N 直接调用FFT子程序计算IFFT的方法： 常系数线性差分方程 数字信号处理DIGITAL SIGNAL PROCESSING 第四章 快速傅立叶变换（FFT） Lecture 4-1 快速傅里叶变换的由来 Lecture 4-2 DFT算法的快速实现方法之一：时间抽取（DIT）基-2 FFT算法 离散傅立叶变换 离散傅立叶变换在数字信号处理中的作用： 物理学、电子类学科、数论、组合数学、信号处理、概率论、统计学、密码学、声学、光学、海洋学、结构动力学等 离散傅立叶变换的应用情况 效果好 计算时间长 直接计算DFT的问题及改进途径 运算量 复数乘法 复数加法 单个X(k) N N – 1 N 个X(k) (k = 0,1,…, N-1) N 2 N (N – 1) 实数乘法 实数加法 一次复乘 4 2 一次复加 2 单个 X(k) 4N 2N + 2 (N – 1) =2 (2N – 1) N 个 X(k) (k=0,1,…, N-1) 4N 2 2N (2N – 1) 傅立叶变换的直接计算方法 结论——直接计算计算量巨大 FFT： Fast Fourier Transform 1965年，Cooley, Tukey在期刊《Mathematics of Computation》上发表著名论文： 《机器计算傅里叶级数的一种算法》 使得傅立叶变换获得实际的应用 问题： 什么样的方法可以大大降低傅立叶变换的计算量？入口在哪里？？？ 计算的重复性！ FFT算法分类: 1、时间抽选法 DIT: Decimation-In-Time 2、频率抽选法 DIF: Decimation-In-Frequency 按时间抽选的基-2 FFT算法 算法原理 设序列点数 N = 2L，L 为整数。若不满足，则补零 N 为 2 的整数幂的 FFT 算法称 基-2 FFT算法。 将序列 x(n) 按 n 的奇偶顺序分成两组： 则 x(n) 的 DFT: 利用周期性，X(k) 可 分两部分进行求解 分解后的运算量： 复数乘法 复数加法 一个N / 2点DFT (N / 2)2 N / 2 (N / 2 –1) 两个N / 2点DFT N 2 / 2 N (N / 2 –1) 一个蝶形 1 2 N / 2个蝶形 N / 2 N 总计 运算量减少近一半！ 若N / 2仍为偶数，进一步分解：N / 2 N / 4 同理： 其中： 这样逐级分解，直到 2 点DFT。 当 N = 8 时，即分解到 X3(k)，X4(k)，X5(k)，X6(k)，k = 0, 1 运算量 当 N = 2L 时，共有 L 级蝶形，每级 N / 2 个蝶形， 每个蝶形有 1 次复数乘法 2 次复数加法。 复数乘法： 复数加法： DFT与FFT的计算量之比 算法特点 1）原位计算 m 表示第 m 级迭代，k，j 表示数据所在的行数 2）倒位序 倒位序 自然序 000 0 0 000 100 4 1 001 010 2 2 010 110 6 3 011 001 1 4 100 101 5 5 101 011 3 6 110 111 7 7 111 n0 n1 n2 0 0 0 1 1 0 1 1 0 0 1 1 0 1 3）蝶形运算 对 N = 2L点FFT，输入倒位序，输出自然序， 第 m 级运算每个蝶形的两节点距离为 2m–1 第 m 级运算： 对 m 级运算，蝶形运算的两节点之间的距离为2m-1,设第一个节点的标号值为 k ，将其乘以 2L-m 并取整，取该整数的低 L 位二进制数，该数即为 r 。 清华大学《计算机文化基础》电子教案 2003年3月 * 页 常系数线性差分方程