FFT算法介绍[修复的]概要.pptxVIP

;本章内容：;本讲学习目标;本章作业练习 ;第四章 快速傅里叶变换;4.2.1直接计算DFT的特点及减少运算量的基本途径;运算量;FFT算法分类:;4.2.2、时间抽取法基-2FFT算法基本思想 （基-2 Decimation-In-Time FFT）;则x(n)的DFT:; ; 同理, 这就是说,X1(k),X2(k)的后一半,分别 等于其前一半的值。; 可见,X(k)的后一半，也完全由X1(k), X2 (k)的前一半所确定。 * N点的DFT可由两个N/2点的DFT来计算。;实现上式运算的流图称作蝶形运算;图4.2.1 蝶形运算符号 ; x1(0)=x(0) x1(1)=x(2) N/2点 x1(2)=x(4) DFT x1(3)=x(6) x2(0)=x(1) x2(1)=x(3) N/2点 x2(2)=x(5) DFT x2(3)=x(7) ;6. N/2分解后的运算量：; 由于N=2 M ,所以 N/2仍为偶数,可以进 一步把每个N/2点的序列再按其奇偶部分 分解为两个N/4的子序列。例如，n为偶数时的 N/2点分解为:;同理可将奇数序号组成的N/2点序列进行分解得：;这样逐级分解，直到2点DFT 当N = 8时，即分解到X3(k)，X4(k)，X5(k)，X6(k)，k = 0, 1 ; (0) (4) (2) (6) (1) (5) (3) (7);4.2.3 基2DIT-FFT算法与直接DFT运算量的比较;;4.2.4 DIF-FFT的运算规律 及编程思想; (0)=A0(0) A1(0) A2(0) A3(0)=X(0) (4)=A0(1) A1(1) A2(1) A3(1)=X(1) (2)=A0(2) A3(2)=X(2) (6)=A0(3) A3(3)=X(3) (1)=A0(4) A1(4) A2(4) A3(4)=X(4) (5)=A0(5) A3(5)=X(5) (3)=A0(6) A3(6)=X(6) (7)=A0(7) A1(7) A2(7) A3(7)=X(7);开始时，输入序列的N个数放于此N个存储器内，倒序重排后仍存于这N个存储器中，每一次迭代运算后的结果也仍然存于这N个存储器中，整个运算完成后，这N个存储器中即为所求的频谱序列X(k) （k = 0、1、…..、 N-1）。这就是所谓的同址计算，这样可以大大节约存储元件。 ; N点DIT―FFT运算流图中，每级都有N/2个蝶形。每个蝶形都要乘以因子WpN，称其为旋转因子，p称为旋转因子的指数。 ;3）蝶形运算规律;旋转因子的确定仅与指数P有关，当L一定时J可以确定，进而可以确定指数P。对于同一P，参与本级蝶形运算的次数为 ，每级第一次调用 的蝶形结第一个输入节点的位置为J，第二次调用 的蝶形结第一个输入节点的位置为J+2L，即以2L为步进，搜索下一蝶形运算第一个输入节点的位置。以此类推，直至做完 个蝶形运算。;DIT-FFT原位计算步骤;; 4） 编程思想及程序框图 在一般情况下，进行FFT运算的序列至少都有几百点的长度，因此需要编制FFT算法程序以便能够利用计算机来快速进行计算。输入倒位序，输出顺序的DIT-FFT的编程思想， N必须等于2的正整数幂，FFT的计算程序可以分为两部分：一部分是倒序重排，另一部分是用三层嵌套的循环来完成M=log2N次迭代。 ; 三层循环的功能是：最里的一层循环完成相同WNP 的蝶形运算，中间的一层循环完成因子WNP的变化，而最外的