数字信号处理第四章 fft.pptVIP

第四章 快速傅立叶变换（FFT) 本章主要内容 掌握FFT算法基本思想和方法 掌握基2DIT-FFT算法、规律及流图 掌握基2DIF-FFT算法和流图 掌握利用DFT进行计算 利用DFT对信号进行谱分析 第四章 快速傅立叶变换（FFT) 概述 DFT是数字信号中的一种重要变换，但从DFT定义可以容易得到直接计算一个N点的DFT需要：N2次复数乘法；N(N-1)次复数加法。即其运算量随着N按平方增加，当N较大时，其计算量非常大，直接用DFT进行实时计算或谱分析是不切实际的。 1965年库利(J.W.Cooley)和图基(J.W.Tukey)发现DFT的快速算法后，DFT才得到实际的应用。 自1965年后，DFT的快速计算算法的研究得到空前的发展，除了Cooley-Tukey算法；Sande-Tukey算法外，还有许多其它算法，如：Winograd算法；余弦变换快速算法；Walsh变换；数论变换等 第四章 快速傅立叶变换（FFT) 基2FFT算法 FFT的基本思想 长为N的序列x(n)的DFT定义： 式中： 旋转因子 旋转因子的周期性和对称性 周期性： 对称性： 基2FFT算法 FFT的基本思想： 利用 的周期性和对称性，可使DFT运算中的某些项合并； 因为DFT的运算量与N2成正比，若将长序列DFT运算尽可能地分解成几个短序列的DFT，这样可以降低运算量 基2时域抽取FFT(Cooley-Tukey算法，DIT-FFT) 基2FFT ：通过补零使N满足： ，M为自然数 时域抽取原理 按n的奇偶将x(n)分解为两个N/2点的子序列： 基2时域抽取FFT(Cooley-Tukey算法，DIT-FFT) 则x(n)的DFT可写作： 再由的 周期性和对称性可求的DFT的后一半： 由周期性： 得： 和 基2时域抽取FFT(Cooley-Tukey算法，DIT-FFT) 再由对称性： 从而有： 这样，一个N点的DFT被分解成了两个N/2点的DFT线性组合： 将DFT分解M次，最后为2点DFT，完成FFT分解。 蝶形运算表示 上式定义的运算称为蝶形运算(Bullerfly Computation)，它可由图4.2.1形象表示，利用蝶形运算符号可将FFT运算用流图描述。 基2时域抽取FFT(Cooley-Tukey算法，DIT-FFT) 一个蝶形运算由一次复乘法；两次复加法实现。向上加；向下减。 N=8的Cooley-Tukey法示例 一个N点基2FFT算法可以通过分解M次，每次用N/2个蝶形运算表示。图4.2.2和图4.2.3分别给出了8点时域抽取FFT的一、二次分解过程，图4.2.4为分解完成后的8点时域FFT流图。 图4.2.3 N点DFT的第二次时域抽取分解图（N=8） 图4.2.4 N点DIT-FFT运算流图（N=8） 基2时域抽取FFT(Cooley-Tukey算法，DIT-FFT) Cooley-Tukey算法的规律及特点 运算次数 基2FFT流图共有M级蝶形；每级有N/2个蝶形；每个蝶形最多需要一次复乘法、二次复加法。这样一个N点基2FFT的运算量最多为： 复乘法： 复加法： 与直接运算比较： 复乘法： 复加法： N越大，FFT效率越高，由图4.2.5显见。 Cooley-Tukey算法的规律及特点 原位运算 在运算中无需中间寄存器。仅需（N+N/2）个存储单元。 蝶距 定义：蝶形输入两信号点间的节点数称为蝶距。 式中：N为点数；M为级数；l为级号。 旋转因子 各级蝶形有 组 ；每组有 个 ，而且组中 的幂m按差 由0递增。 输入序列的倒序 按M位二进制“码位倒置”规律扰乱输入序列的角标顺序。 例： （N＝8见表4.2.1 ） 基2FFT算法 IDFT的运算 将 改为 ，计算完后再乘1/N。此法需要修改FFT子程序； 由IDFT表达式有： 先将X(k)取共轭，然后直接调用FFT子程序（或用与正变换相同的专用硬件），再将结果取共轭并乘以1/N。 基2频域抽取FFT(Sande-Tukey算法 ，DIF-FFT) 抽取原理 将x(n)分成前N/2和后N/2两半，即： 基2频域抽取FFT(Sande-Tukey算法 ，DIF-FFT) 这样，DFT可