第四章快速傅理叶变换FFT01-11课题.ppt

第四章  快速傅里叶变换（FFT） 快速傅里叶变换 （FFT） §4.1 DFT运算的特点： §4.2 按时间抽取的FFT (3) 和 可以继续分解下去，可将N/2点的 子序列再按奇、偶项分解，一直到最后分解成  两两点的DFT为止。 时间抽取法FFT的运算特点： （1）蝶形运算 算法的运算速度 例：以N=8 按频率抽取的FFT 另外一种IFFT算法的实现 §4.5 任意基数的FFT算法 §4.6 Chirp-z变换（简称CZT) 将x(n)按n的顺 序前后对半分开： §4.3 按频率抽取的FFT 频率抽取法是N=2M情况下的另外一种FFT算法 按时间抽取FFT算法的基本思路 按频率抽取FFT算法的基本思路 x(n)按奇、偶一级级分开， X(k)按前一半、后一半一级级分开。 X(k)按奇、偶一级级分开， x(n)按前后对半一次次分开。 一.按频率抽取法的FFT算法原理 按k的奇、偶把X(k)进一步分解为两部分： a) k为偶数，令k=2r b) k为奇数，令k=2r+1 按频率抽取的蝶形运算流图： 把8点DFT分解成两个4点的DFT K为偶数 X(2r) K为奇数 X(2r＋1) 将两个4点DFT再分解成四个2点DFT 由于N/2仍是偶数可以继续分解，直到分解为两两点的DFT为止。 下图是N=8按频率抽取的FFT流图 频率抽取法的分解过程： (1)把输入序列前后对半分； (2)输出按k值的奇、偶将X(k)分解为奇数组和偶数组两部分；(3)输入如果按顺序输入，输出则是码位倒置形式。 按频率抽取法FFT的运算特点： （1）蝶形运算的运算量 其运算量与按时间抽取运算量相同。即： 复乘： 复加： （2）输入如果是自然顺序，输出则是码位倒置形式。如果要求输出为自然顺序，则需要进行变址运算。按频率抽取的FFT流图可以由按时域抽取的FFT流图转置 而得到。转置指将流图中所有支路的方向反 ,并将输入变输出，输出变输入。 频率抽取法的基本蝶形图与时间抽取法的基本蝶形图有所不同。时间抽取法是先乘后加、减，频率抽取法是先加、减，再乘系数 按频率抽取的 蝶形运算流图 按时间抽取的 蝶形运算流图 按时间抽取的FFT算法以及按频率抽取的FFT算法是DFT的两种快速算法，这些算法同样可以用於IDFT，我们称为快速傅立叶反变换(IFFT)。我们把IDFT公式与DFT公式作以比较： (1) DFT中的系数 (2) IDFT中要乘以常数 §4.4 IDFT的快速算法IFFT 以上讨论的时间抽取或频率抽取的FFT运算均可直接用于IFFT运算。但需要将蝶形系数 改为 ，每级乘以系数1/2 。 IFFT计算分为三步： ① 将X(k)取共轭(虚部乘以-1) ; ② 对 直接作FFT ； ③ 对FFT的结果取共轭并乘以1/N，就得到x(n)。 按上述方法求得的IFFT，完全不需要改动FFT程序，FFT和 IFFT可以共用一个程序。 总结FFT和IFFT蝶形图的表示方法： 按时间抽取的FFT 变到IFFT时，对应的流图应该是按频率抽取的IFFT。因为： 按时间抽取的FFT 输入变量是x(n) x(n)按奇、偶分解 按频率抽取的IFFT 输入变量是X(k) X(k)按奇、偶分解 同理，按频率抽取的FFT 变到IFFT时，对应的流图应该是按时间抽取的IFFT。 （1）输入按奇、偶分解： （2）输入按前、后分解： 前面讨论的FFT算法，序列的长度均满足N=2M ，这种情况实际使用得最多，因为它程序简单，效率 高，使用方便。如果不满足N=2M, 处理方法有两种: ① 补零: 将x(n)补零,使N增加到临近的一个值， 满足N=2M 。 例如：N=30, 补上x(30)=x(31)=0 两个零点, 使 得满足N=25=32, 这样可直接调用以2为基数的FFT 运算。有限长序列补零后并不影响其频谱 ， 补零只是频谱的采样点数增加了。 * ? 快速傅里叶变换（FFT）是计算DFT的一种快速高效的算法。 有限长序列在数字技术中占有很重要的地位，主要原因是由于其频谱可以离散化。有限长序列的DFT本身可以完全表达序列的频谱，所以DFT也可以直接对信号进行频谱分析。 谱分析在通信技术、图象传输、语音压缩、生物医学等领域都得到应用。 虽然频谱分析和DFT运算很重要，但在很长一段时间里，由于DFT运算量太大，因此它并没有得到真正的运用，频谱分