数字信号处理第2版张小虹4数字信号处理1课件教学.pptVIP

例 N=210=1024 例 N=211=2048 ， Nlog2N N2 M N = 10 1024 = =102.4 N2 =4,194,304 Nlog2N=22,528 Nlog2N N2 M N = 11 2049 = =186.2 完成了。FFT算法的重大突破意义十分明显。 若直接计算要180分钟，则利用FFT只用不到一分钟就 4.2.2、基2时选 FFT运算规律 1、同址计算（原位） 可见N 越大，效率越高。 在整个计算中共有M级，每级有N/2 个蝶形 ，第L级 时选蝶形 L?1级 XL-1(p) XL-1(q) L级 XL(p) XL(q) 基本运算蝶形如图4.2-7所示。 WN p L级基本蝶形的运算关系为 从图4.2-7可见，为了计算第L列的p,q位置上的复数 节点值，只需要L -1列的p,q位置上的复数节点值。 如果计算得到的XL(p)、XL(q )又分别存在原来XL-1(p) XL-1(q )同一存储器内，则为完成整个计算实际只需要 一列(N个)复数存储单元。 XL (q) =XL?1 (p)? WN XL?1(q) p XL (q) =XL?1 (p)+ WN XL?1(q) p (4.2-17) 每个蝶形计算的输入、输出节点水平相邻时，该流 图对应的是同址计算。否则就需要两列复数存储器。 因为只要一组(N个)复数存储器，因此可以减少内 存，降低成本。从蝶形图看，只有当节点排列成使 因子的指数。从图4.2-6可见，运算级数不同，旋转因 子和循环方式就不同。我们要确定各级的旋转因子和 2、时选FFT系数的运算规律 用L表示从左到右的运算级数(L=1, 2 ,…, M) 。 从图4.2-6可见，第L级共有2L-1个不同的旋转因子。 N点FFT运算流图共有M级，每级有N/2个蝶形。 每个蝶形都要乘以因子WN ，称为旋转因子，p为旋转 p 循环方式，就要找出WN 旋转因子与运算级数的关系。 p 例N=8=23，M=3，各级旋转因子表示如下： L=1 L=2 ， L=3 ， WN = WN/4 ， p=0，一个旋转因子 r p WN = WN/2 ， p=0，1两个旋转因子 r p WN = WN ， p=0，1，2，3四个旋转因子 r p 其旋转因子为WN 0 其旋转因子为WN 、 WN 0 2 其旋转因子为WN 、 WN 、 WN、 WN 0 3 1 2 所以第L级的旋转因子为 推广到N=2M时的一般情况，因为2L =2M · 2L-M =N 2L-M， 可由(4.2-18)、(4.2-18)确定第L级的旋转因子。 在此基础上可以得到第L级蝶形一般运算规律。 (4.2-18) p =0,1,2,…, 2L?1 ?1 (4.2-19) 式中 r =p·2M?L WN = W = WN = W r 2L N·2L?M p p·2M?L p 第L级蝶形一般运算规律 p =0,1,2,…, 2L-1 ?1 式中 r =p·2M-L B =2L-1 L=0,1,2,…, M 设序列x(n)经时间抽取后，存入数组X中。如果蝶形 运算的两个输入数据相距B点，应用同址计算，则第 L级蝶形运算输出表示 XL (q) =XL (p+B) =XL?1 (p)? WN XL?1(p+B) r XL (q) =XL?1 (p)+ WN XL?1(p+B) r 从图4.2-6及上面的分析，可以归纳出时选FFT的运算 形。根据上述规律，从输入端（第一级）开始，逐级 进行，共作M级运算。在作第L级运算时，依次求出 FFT运算。程序流图如图4.2-8所示。 规律：第L级每个蝶形的两个输入数据相距 B =2L-1 个点；同一旋转因子对应着间隔为 2L点的 2M-L 个蝶 它对应的所有2M-L个蝶形。这样可以用三重循环实现 个不同的旋转因子，每求出一个旋转因子，就计算完 结束 输 出 ② 倒 序 开始 ① 输入x(n),M N=2M B?2L-1 L=1,M ② ① X(k) ?X (k)+WN X(k+B) p X(k+B) ?X (k)?WN X(k+B) p J=0,B?1 k=J,N?1,2L P=2M?LJ 的，实际上的输入是“倒序位”存储。仍以8点的DFT流图 序中将长点的DFT（按奇、偶）分解为短点DFT引起 而是x(0), x(4) ,x(2),… x(N-1) 。这是由时选法不断在时 按顺序输出。但输入的x(n)不是按这种自然顺序排列的 正好顺序放着X(0), X(1) ,X(2),… X(N-1) 。因此可以直接 由时间抽取FFT可见，当全部FFT完成后，存储单元中 4、变址 N=8=23 要用三位二进制码标明数列的次