数字信号处理时间抽取FFT分析报告.ppt

数字信号处理(Digital Signal Processing) 离散傅里叶变换快速算法(FFT) 问题的提出 解决问题的思路与方法 基2时间抽取FFT算法 基2频率抽取FFT算法 FFT算法的实际应用—— 实序列的DFT计算，IDFT的快速计算方法 问题的提出 解决问题的思路 旋转因子 的性质 解决问题的方法 基2时间抽取FFT算法 基2时间抽取FFT算法推导 基2时间抽取FFT算法推导 基2时间抽取FFT算法流图 4点基2时间抽取FFT算法流图 4点基2时间抽取FFT算法流图 8点基2时间抽取FFT算法流图 FFT算法流图旋转因子 规律 8点基2时间抽取FFT算法流图 算法的计算复杂度 原位运算 倒序 倒序的实现——变址 补零,插零序列的DFT 时间抽取FFT * * 4点序列{2，3，3，2} DFT的计算复杂度 复数加法 N(N-1) 复数乘法 N 2 如何提高DFT的运算效率? 1. 将长序列DFT分解为短序列的DFT 2. 利用旋转因子 的周期性、对称性、可约性。 (1) 周期性(以N为周期) (2) 对称性 (3) 可约性 将时域序列逐次分解为一组子序列，利用旋转因子的特性，由子序列的DFT来实现整个序列的DFT。 基2时间抽取(Decimation in time)FFT算法 基2频率抽取(Decimation in frequency)FFT算法 基2时间抽取FFT算法推导 基2时间抽取FFT算法流图 基2时间抽取FFT算法的计算复杂度 基2时间抽取FFT算法流图规律 因此有: 由于X1[m] 和X2[m]隐含有周期性，可得 N=2 x[k]={x[0], x[1]} x[0] x[2] x[1] x[3] X1[0] X1[1] X2[0] X2[1] 2点DFT 2点DFT -1 -1 -1 -1 X [0] X [1] X [2] X [3] 4点DFT 4点DFT x[0] x[2] x[4] x[6] x[1] x[3] x[5] x[7] X1[0] X1[1] X1[2] X1[3] X2[0] X2[1] X2[2] X2[3] X [0] X [1] X [2] X [3] X [4] X [5] X [6] X [7] -1 -1 -1 -1 4点DFT 4点DFT x[0] x[2] x[4] x[6] x[1] x[3] x[5] x[7] X1[0] X1[1] X1[2] X1[3] X2[0] X2[1] X2[2] X2[3] X [0] X [1] X [2] X [3] X [4] X [5] X [6] X [7] -1 -1 -1 -1 8点基2时间抽取FFT算法流图 第一级 第二级 第三级 8点基2时间抽取FFT算法流图 第二级的蝶形系数为 ，蝶形节点的距离为2。 第一级的蝶形系数均为 ，蝶形节点的距离为1。 第三级的蝶形系数为 ，蝶形节点的距离为4。 第M级的蝶形系数为 ，蝶形节点的距离为N /2。 第一级 第二级 第三级 复乘次数 复乘次数 N N 2 x[0] x[4] x[2] x[6] x[1] x[5] x[3] x[7] A(1) A(2) A(3) A(4) A(5) A(6) A(7) A(8) 输入序列 存储单元 第一级输出 第二级输入 第二级输出 第三级输入 X1[0] X1[1] X2[0] X2[1] X3[0] X3[1] X4[0] X4[1] A(1) A(2) A(3) A(4) A(5) A(6) A(7) A(8) X5[0] X5[1] X5[2] X5[3] X6[0] X6[1] X6[2] X6[3] A(1) A(2) A(3) A(4) A(5) A(6) A(7) A(8) X [0] X [1] X [2] X [3] X [4] X [5] X [6] X [7] A(1) A(2) A(3) A(4) A(5) A(6) A(7) A(8) 第三级输出 k0 k1 k2 x [ k 2 k 1 k 0 ] x [00 0] x [10 0] x [01 0] 0 1 0 1 1] 1 2 x [ k k 0] x [ k 2 k 1 0 1 x [11 0] x [001 ] x [101 ] x [01 1] x [111 ] 0 1 0 1 0 1 0 1 A(1) A(2)