第4章快速傅里叶变换概要.pptxVIP

;2;3;4;;;;8;9;; ;12;13;14;;16;;;;20;;;;;; X1(k)、X2(k)都是 N/2 点的 DFT，它们各自又可分成 N/4 点的DFT，如此继续分下去，直至两点DFT。两点DFT 不需要乘法运算：;;28;29;;31;32;33; 从图可以看出，整个算法流图可以分为四段，（0）段为倒序重排，后面3段为3(log28=3)次迭代运算：首先计算2点DFT，然后将2点DFT的结果组合成4点DFT，最后将4点DFT组合为8点DFT。因此，对于N点FFT，只需要一列存储N个复数的存储器。 ;开始时，输入序列的N个数放于此N个存储器内，倒序重排后仍存于这N个存储器中，每一次迭代运算后的结果也仍然存于这N个存储器中，整个运算完成后，这N个存储器中即为所求的频谱序列X(k) （k = 0、1、…..、 N-1）。这就是所谓的同址计算，这样可以大大节约存储元件。 ;36; 倒位序：输入序列经过M级分解后，其排列顺序变为x(0) x(4) x(2) x(6) x(1) x(5) x(3) x(7) ，服从所谓的“码位倒置”的规律，称之为倒位序。 造成倒位序的原因： 将其按标号的偶奇的不断分组， 每次分解总是将偶序列放在上面， 把奇序列放在下面。 首先最低位按0、1分为偶、奇两组， 接着次低位也按0、1分组， 依此类推 右图为描述倒位序的树状图(N=8) ;可以看出，将输入序号的二进制表示（n2n1n0）位置颠倒，得到（n0n1n2），就是相应的倒序的二进制序号。因此，输入序列按倒序重排，实际上就是将序号为（n2n1n0）的元素与序号为（n0n1n2）的元素的位置相互交换。 ;39;40; 4) 旋转因子的变化规律 旋转因子和运算级数的关系 N=23=8时的各级旋转因子： 一般地，当N=2M （M为正整数）时，第L级蝶形运算共有 2L-1种不同的旋转因子，可表示为 这样，就可以确定第L级蝶形运算的旋转因子(程序中，L为循环变量)。 ;旋转因子的确定仅与指数P有关，当L一定时J可以确定，进而可以确定指数P。 对于同一P，参与本级蝶形运算的次数为2M-L，每级第一次调用 的蝶形结第一个输入节点的位置为J，第二次调用 的蝶形结第一个输入节点的位置为J+2L，即以2L为步进，搜索下一蝶形运算第一个输入节点的位置。以此类推，直至做完2M-L 个蝶形运算。;DIT-FFT原位计算步骤;关于FFT算法的计算机程序 在一般情况下，进行FFT运算的序列至少都有几百点的长度，因此需要编制FFT算法程序以便能够利用计算机来快速进行计算。可以参照下图来编制计算N (N必须等于2的正整数幂)点FFT的计算机程序，整个程序可以分为两部分：一部分是倒序重排，另一部分是用三层嵌套的循环来完成M=log2N次迭代。 ;三层循环的功能是：最里的一层循环完成蝶形运算，中间的一层循环完成因子WNk的变化，而最外的一层循环则是完成M次迭代过程。 倒序重排的程序是一段经典程序，它以巧妙的构思、简单的语句用高级编程语言来完成了倒序重排的功能。;46;47;48;各级中W的排列规律（自上而下）;50;;52;;; 图4-11 DIF―FFT一次分解运算流图（N=8）;56;57;58;;60;61;62;63;64;65;66;67;68;69;70;这些z平面上的抽样点所沿的周线是一条螺旋线，参数W0控制周线盘旋的倾斜率。若W0大于1，则随着k的增加，周线向内盘旋，趋向原点；若W0小于1，则随着k的增加，周线向外盘旋。若W0=1，螺旋线实际上是一段圆弧，而如果又有A0=1，则这段圆弧是单位圆的一部分。A0和θ0分别为第一个抽样点（k = 0）的模和幅角，其余抽样点沿螺旋周线按角度间隔φ0分布。;要计算： 式中N为序列x(n) 的长度。由恒等式： 并且令：