数字信号处理第4章剖析.ppt

　　形成倒序J后，将原数组A中存放的输入序列重新按倒序排列。 (1) 原输入序列x(I)先按自然顺序存入数组A中。 对N=8，A(0)，A(1)，…，A(7)中依次存放着x(0)，x(1)，x(2)，…，x(7) (2) 对x(n)重新排序(倒序)：规律如图4.2.8所示： ① 第一个序列值x(0)和最后一个序列值x(N－1)不需要重排； ② 每计算出一个倒序值J，便与循环语句自动生成的顺序I 比较：当I=J时，不需要交换数据； 当I≠J时, A(I)与A(J)交换数据。 ③为了避免再次调换前面已调换过的一对数据，只对IJ的 情况调换A(I)和A(J)的内容。 图4.2.8 倒序规律 图4.2.9 倒序程序框图 4.2.5 频域抽取法FFT(DIF-FFT) 　　设序列x(n)长度为N=2M，首先将x(n)前后对半分开，得到两个子序列，其DFT如下： 式中 将X(k)分解成偶数组与奇数组： 当k取偶数(k=2m, m=0, 1, …, N/2－1)时 (4.2.14) 当k取奇数(k=2m+1, m=0, 1, …, N/2－1)时， 令 将x1(n)和x2(n)分别代入(4.2.14)和(4.2.15)式，可得 　　　　　　　　　　　　　　　 (4.2.16) ——X(k)按k值的奇偶分为两组：偶数组是x1(n)的N/2点DFT， 奇数组是x2(n)的N/2点DFT。 x1(n)、x2(n)和x(n)之间的关系也可用蝶形运算流图符号表示 图4.2.11 DIF-FFT第一次分解运算流图（N=8） 　　由于N=2M，N/2仍然是偶数，继续将N/2点DFT分成偶数组合奇数组，这样每个N/2点DFT又可由两个N/4点DFT形成，其输入序列分别是x1(n)和x2(n)按上下对半分开形成的四个子序列。图4.2.12示出了N=8时第二次分解运算流图。 这样继续分解下去，经过M－1次分解，最后分解为2M－1个两点DFT，两点DFT就是一个基本蝶形运算流图。当N=8，经两次分解，便分解为四个两点DFT。N = 8的完整DIF-FFT运算流图如图4.2.13所示。 图4.2.12 DIF-FFT第二次分解运算流图（N = 8） 图4.2.13 DIF-FFT运算流图（N =8） 　　这种算法是对X(k)进行奇偶抽取分解的结果，所以称之为频域抽取法FFT。观察图4.2.13可知： DIF-FFT算法与DIT-TTF算法类似，可以原位计算； DIF-FFT算法共有M级运算，每级共有N/2个蝶形运算，所以两种算法的运算次数亦相同； 不同的是DIF-FFT算法输入序列为自然顺序，而输出为倒序排列； M级运算完后，要对输出数据进行倒序才能得到自然顺序的X(k)； 蝶形运算略有不同，DIT-FFT蝶形先乘后加(减)，而DIF-FFT蝶形先加(减)后相乘。 4.2.6 IDFT的高效算法 　　上述FFT算法流图也可以用于计算IDFT。比较DFT和IDFT的运算公式： 只要将DFT运算式中的系数　　　改变为　　　，最后乘以1/N，就是IDFT运算公式； 现在流图的输入是X(k)，输出就是x(n)； 原来的DIT-FFT→DIF-IFFT； DIF-FFT→DIT-IFFT。 　　如果希望直接调用FFT子程序计算IFFT，则可用下面的方法： 　　由于 1、将X(k)取复共轭； 2、直接调用FFT子程序，或者送入FFT专用硬件设备进行DFT运算； 3、最后取复共轭并乘以1/N得到序列x(n)。 这种方法虽然用了两次取共轭运算，但可以与FFT共用同一子程序，因而用起来很方便。 4.3 进一步减少运算量的措施 4.3.1 多类蝶形单元运算 DIT-FFT运算：N=2M点FFT共需要MN/2次复数乘法 由(4.2.12)式， 当L=1时，只有一种旋转因子 ，第一级不需要乘法运算； 当L=2 时，共有两个旋转因子： 和 ，第二级也不需要乘法运算； 依此类推 在DFT中，又称其值为±1和±j的旋转因子为无关紧要的旋转因子，如 ， ， 等。 综上所述： (1) 先除去第一、二两级后，所需复数乘法次数应是 　　　　　　　　　　　　　　　　　　(4.3.1) (2) 当L=3时，有两