数字信号处理4～分析报告.ppt

第4章 快速傅立叶变换(FFT) 直接计算DFT的问题和改进的途径 按时间抽取(DIT)的基2 FFT算法 按频率抽取(DIF)的基2 FFT算法(选学) 线性调频z变换（Chirp-z变换）算法(选学) 线性卷积与线性相关的FFT算法 数字信号处理的实现（DSP课,选学） 1.知识回顾 有限长序列的离散傅立叶变换(DFT) 2.直接计算DFT的问题和改进的途径 (1)、DFT的计算工作量 以正变换为例： 2.直接计算DFT的问题和改进的途径 在DSP芯片中，没有复数运算，因此，复数运算需转换成实数运算实现。 2.直接计算DFT的问题和改进的途径 (2)、改进的途径 a、 的性质 2.直接计算DFT的问题和改进的途径 b、有关系数关系 3.按时间抽取(DIT)的基2 FFT算法 一、算法原理 (一) N/2点DFT 1、先将x(n)按n的奇偶分为两组作DFT，设N=2L ,不足时可补零。这样有: n为偶数时: n为奇数时: 3.按时间抽取(DIT)的基2 FFT算法 2、两点结论: (1) 、 均为N/2点的DFT。 (2) 只能确定出 的N/2个值。 3.按时间抽取(DIT)的基2 FFT算法 3、X(k)的后一半的确定 由于 （周期性）,所以: 3.按时间抽取(DIT)的基2 FFT算法 3.按时间抽取(DIT)的基2 FFT算法 4、蝶形运算 3.按时间抽取(DIT)的基2 FFT算法 5、运算量分析 按奇、偶分组后的计算量： 3.按时间抽取(DIT)的基2 FFT算法 总共运算量: 3.按时间抽取(DIT)的基2 FFT算法 例如 N=8时的DFT,可以分解为两个N/2=4点的DFT。 具体方法如下: (1)n为偶数时,即 分别记作: 3.按时间抽取(DIT)的基2FFT算法 (2) n为奇数时,分别记作: 3.按时间抽取(DIT)的基2 FFT算法 整个过程如下图所示: 3.按时间抽取(DIT)的基2 FFT算法 (二) N/4点DFT 3.按时间抽取(DIT)的基2 FFT算法 3.按时间抽取(DIT)的基2 FFT算法 同样对n为奇数时，N/2点分为两个N/4点的序列: 3.按时间抽取(DIT)的基2 FFT算法 从而有： 3.按时间抽取(DIT)的基2 FFT算法 3.按时间抽取(DIT)的基2 FFT算法 例如，N=8时的DFT可分解为四个N/4的DFT， 具体步骤如下: （1）序列分解。 3.按时间抽取(DIT)的基2 FFT算法 同样： 3.按时间抽取(DIT)的基2 FFT算法 （2）蝶形运算 由 X3(0), X3(1), X4(0), X4(1)进行碟形运算, 得到X1(0), X1(1), X1(2), X1(3)。 由 X5(0), X5(1), X6(0), X6(1)进行碟形运算,得到X2(0), X2(1), X2(2), X2(3)。 由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) 3.按时间抽取(DIT)的基2 FFT算法 整个过程如下图所示: 3.按时间抽取(DIT)的基2 FFT算法 这样,又一次分解,得到四个N/4点DFT,两级蝶形运算,其运算量有大约减少一半。 对于N=8时DFT,N/4点即为两点DFT,也可以用蝶形运算实现，如下所示： 3.按时间抽取(DIT)的基2 FFT算法 也即： 3.按时间抽取(DIT)的基2 FFT算法 因此，8点DFT的FFT的运算流图如下 3.按时间抽取(DIT)的基2 FFT算法 注意系数关系： （1）第一级，加权系数只有1个 ，没有乘法。 （2）第二级，系数为2个， ，也可以不用乘法计算。 （3）在最后一级(L级），系数为N/2个， 3.按时间抽取(DIT)的基2 FFT算法 这种FFT算法,是在时间上对输入序列的次序是属于偶数还是属于奇数来进行分 解的，所以称作按时间抽取的算