数字信号处理课件(第4章 快速傅里叶变换补充).pptVIP

* 3.7 基4 FFT和分裂基FFT算法 1。Radix 4算法 (1)算法原理：设 N=43=64，用DIT—FFT做第一次分解， 令N=16×4, 0≤m0≤15，0≤n2≤3 0≤ ≤15，0≤k2≤3 点DFT 点DFT 点DFT 点DFT … … … … x(4m0) x(4m0+1) x(4m0+2) x(4m0+3) ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ … … … … ○ ○ ○ ○ ○ ○ ○ ○ X(0) X(15) X(16) X(31) X(32) X(47) X(48) X(63) n2=1 n2=0 n2=2 n2=3 =0,1,2…15 X( ) X( + ) X( + ) X( + ) (2) 基4算法特点： ① 若点数N=4m, 则基4分解的级数是m级，且每级间都是 同位运算； ② 乘法次数：基4蝶形运算单元为： ? 每一基4蝶形运算 单元 需三次复乘 （用于级之间的旋转 因子）； ? 每级共有 个基4 蝶形； ? 共有m级， ○ ○ ○ ○ -j -1 j -1 -1 j -1 -j ∴ 基4运算所需的复乘次数是： ×(基2复乘次数) 结论：基4运算比基2运算节省 25% 的复乘次数，但复加次 数增加 50%； ③ 基4运算的基2排序： 对：N=4, 用基2倒位序画出的流图和基4比较 (p.145图3.7.1) 点DFT 点DFT 点DFT 点DFT … … … … ○ ○ ○ ○ ○ ○ ○ ○ … … … … X(0) X(15) X(16) X(31) X(32) X(47) X(48) X(63) =0,1,2…15 o o o o o o o o o o o o o o o o o o o o o o o o o o o －j －1 －1 －1 －1 x(4m0) x(4m0+2) x(4m0+1) x(4m0+3) n2=2 n2=0 n2=1 n2=3 X( ) X( + ) X( + ) X( + ) 2。分裂基FFT算法(Split—Radix FFT Algorithm) 1)分裂基FFT算法的提出（p.145）及分裂基FFT算法的 定义（p.146） 当N=2（2r+1）=25=32就难以用基4计算； 2）DIT—FFT分裂基算法原理： x(n)分成三个子序列：x1(r)=x(2r) ——0≤r≤ －1 x2( )=x(4 +1) —0≤ ≤ －1 x3( )=x(4 +3) —0≤ ≤ －1 X(k)=X1((k))N/2+ X2((k))N/4+ X3((k))N/4 ——0≤k≤N－1 ? 第一次分解的流图如p.147图3.7.2 ? 8点DIT分裂基FFT算法流图 o————o——o————o——o————o o————o——o————o——o————o o————o——o————o——o————o o————o——o————o——o————o o————o——o————o——o————o o————o——o————o——o————o o————o——o————o——o————o o————o——o————o——o————o ——(0≤k≤ －1) X(0) X(4) X(2) X(6)