chap快速付里叶变换FFT2解读.ppt

对于任何流图，只要保持各节点所连的支路及传输系数不变，则不论节点位置怎么排列所得流图总是等效的，所得最后结果都是x(n)的DFT的正确结果，只是数据的提取和存放的次序不同而已。这样就可得到按时间抽取的FFT算法的若干其他形式流图。 六、按时间抽取的FFT算法的其他形式流图 ① 数据存放的方式不同 ②取用系数的顺序不同 图 时间抽取、 输入自然顺序、 输出倒位序的FFT流图 第四节 按频率抽取（DIF）的基 -2 FFT算法 一、 算法原理 仍设序列点数为N=2M，M为正整数。在把输出X(k)按k的奇偶分组之前，先把输入序列按前一半、后一半分开（不是按偶数、 奇数分开）， 把N点DFT写成两部分。 k=0, 1, …, N-1 k=0, 1, …, N-1 当k为偶数时，(-1)k=1；k为奇数时，(-1)k=-1。因此，按k的奇偶可将X(k)分为两部分: 所以为X(2r)为前一半输入与后一半输入之和的N/2点DFT， X(2r+1)为前一半输入与后一半输入之差再与WNn之积的N/2点DFT。 令: 图 频率抽取法蝶形运算单元 把一个N点DFT按k的奇偶分解为两个N/2点的DFT了。 N=8时，上述分解过程示于下图。 图 按频率抽取，将N点DFT分解为两个N/2点DFT的组合 与时间抽取法的推导过程一样，由于N=2M，N/2仍是一个偶数，因而可以将每个N/2点DFT的输出再分解为偶数组与奇数组，这就将N/2点DFT进一步分解为两个N/4 点DFT。 这两个N/4点DFT的输入也是先将N/2点DFT的输入上下对半分开后通过蝶形运算而形成的，下图出了这一步分解的过程。 图按频率抽取，将N点DFT分解为两个N/4点DFT的组合 这样的分解可以一直进行到第M次（N=2M），第M次实际上是做两点DFT，这N/2个两点DFT的N个输出就是x(n)的N点DFT的结果X(k)。下图表示一个N=8 完整的按频率抽取的基-2 FFT运算结构。 图 按频率抽取的FFT（N=8）信号流图 初看起来，DIF法与DIT法的区别是： DIF输入是自然顺序，输出是倒位序的，这与DIT法正好相反。但这不是实质性的区别，因为DIF法与DIT法一样，都可将输入或输出进行重排，使二者的输入或输出顺序变成自然顺序或倒位序顺序。DIF的基本蝶形与DIT的基本蝶形则有所不同，这才是实质的不同，DIF的复数乘法只出现在减法之后，DIT则是先作复乘后再作加减法。 频率抽取法(DIF)蝶形运算单元 时间抽取法(DIT)蝶形运算单元 但是，DIF与DIT就运算量来说则是相同的，即都有M级（列）运算，每级运算需N/2个蝶形运算来完成，DIF法与DIT法都可进行原位运算。频率抽取FFT算法的输入是自然顺序，输出是倒位序的。因此运算完毕后，要通过变址计算将倒位序转换成自然位序，然后再输出。转换方法与时间抽取法相同。 第五节 离散傅立叶反变换（IDFT）的快速计算方法 0≤k≤N-1 0≤n≤N-1 需要稍稍改动FFT的程序和参数才能实现IFFT算法。只要将DFT运算中的每个系数 换成 ，最后再乘以常数 ，则以上所有按时间抽选或按频率抽选的FFT都可以拿来计算IFFT。 即先将X(k)取共轭，就可直接利用FFT子程序，最后再将运算结果取一次共轭，并乘以1/N,即得到x(n)值。因此FFT运算和IFFT运算就可以共用一个子程序块。 4.10 线性卷积与线性相关的FFT算法 线性卷积的FFT算法——快速卷积 以FIR滤波器为例，因为它的输出等于有限长单位脉冲响应h(n)与有限长输入信号x(n)的离散线性卷积。 设x(n)为L点，h(n)为M点， 输出y(n)为 y(n)也是有限长序列，其点数为L+M-1 点。 下面首先讨论直接计算线性卷积的运算量。 由于每一个x(n)的输入值都必须和全部的h(n)值相乘一次，因而总共需要LM次乘法，这就是直接计算的乘法次数，以md表示为 md=LM 用FFT算法也就是用圆周卷积来代替这一线性卷积时，为了不产生混叠，其必要条件是使x(n)，h(n)都补零值点，补到至少N=M+L-1， 即: 0≤n≤L-1 L≤n≤N-1 0≤n≤M-1 M≤n≤N-1 然后计算圆周卷积 N 这时，y(n)就能代表线性卷积的结果。