DSP_13离散傅立叶变换-FFT要点.pptVIP

数字信号处理 陈鹏 3299340 快速傅立叶变换 Fast Fourier Transform (FFT) Lesson 13 复习提问 如何利用循环卷积计算线性卷积？ 什么是频率取样，其特点是？ DFT的计算量是多少？ DFT的计算量 DFT的计算量 DFT的计算量 DFT的计算量 DFT的计算量 DFT的计算量 DFT的计算量 时间抽取基2 FFT算法(库里－图基算法) 时间抽取基2 FFT算法(库里－图基算法) 时间抽取基2 FFT算法(库里－图基算法) 时间抽取基2 FFT算法(库里－图基算法) 时间抽取基2 FFT算法(库里－图基算法) 时间抽取基2 FFT算法(库里－图基算法) 时间抽取基2 FFT算法(库里－图基算法) 时间抽取基2 FFT算法(库里－图基算法) 时间抽取基2 FFT算法(库里－图基算法) 时间抽取基2 FFT算法(库里－图基算法) 时间抽取基2 FFT算法(库里－图基算法) 时间抽取基2 FFT算法(库里－图基算法) 时间抽取基2 FFT算法(库里－图基算法) 时间抽取基2 FFT算法(库里－图基算法) 离散傅立叶变换(DFT)在实际应用中是非常重要的，利用它可以计算信号的频谱、功率谱和线性卷积等。但当 N 很大时，利用DFT直接计算运算量太大，因此，如何提高计算DFT的速度成为很重要的课题。1965年库利(Cooley)和图基(Tukey)在前人研究成果的基础上提出了快速计算DFT的算法，之后，又出现了各种各样快速计算DFT的方法，这些方法统称为快速傅立叶变换(Fast Fourier Transform 或 FFT)。 介绍FFT之前，先估计一下直接计算DFT所需的运算量。先展开计算DFT为方程组形式如下 从上式可以看出，假设 是复数，则计算一个 需要 N 次复数乘法和 (N-1) 次复数加法。 所以计算 N 个 值，共需 次复数乘法和 次复数加法。而且我们知道： 即每次复数乘法包括 4 次实数乘法和 2 次实数加法，每次复数加法包括 2 次实数加法，因此计算 N 点 DFT 共需要 次实数乘法和 次实数加法。当 N 很大时，这是一个非常大的计算量。 N (N – 1) N 2 N个X(k)(N点DFT) N – 1 N 一个X(k) 复数加法 复数乘法 2N+2 (N – 1)=2 (2N – 1) 4N 一个X (k) 2N (2N – 1) 4N 2 N个X (k)(N点DFT) 2 一次复加 2 4 一次复乘 实数加法 实数乘法 FFT算法主要利用了 的两个性质： 对称性 周期性 r 为整数。左图是 N=8 时这两个性质的情况。 FFT算法有很多种，我们只介绍其中两种，即时间抽选基2FFT算法和频率抽选基2FFT算法，这两个名字的含义在介绍它们时解释。这两种算法特别适用于 N 等于 2 的幂次的情况。 这种算法简称为时间抽取FFT算法，其基本出发点是，利用旋转因子 的对称性和周期性，将一个大的DFT分解成一些逐次变小的DFT来计算。分解过程遵循两条规则：(1) 对时间进行奇偶分解；(2) 对频率进行前后分解。 设 ，M 为正整数。为了推导方便，取 ，即离散时间信号为 按照规则(1)，将序列 分为两组，一组序号为偶数，另一组序号为奇数，即 并分别表示为 并利用性质 重写DFT如下： * *