数字信号处理_9.pptVIP

离散傅里叶变换快速算法(FFT);学 习 要 求;重 点 和 难 点 ;FFT: Fast Fourier Transform;直接计算DFT的问题及改进途径; 通常x[k]和WNkm都是复数，所以计算一个 X[m]的值需要N次复数乘法运算和N-1次 复数加法运算。;所有的X[m]就要N2次复数乘法运算，N(N-1)次复数加法运算。;小结运算量;二.改进的途径;2. WNkm的可约性和特殊值; 利用上述特性，可以将有些项合并，并将DFT长序列分解为短序列，从而降低运算次数，提高运算速度。;FFT算法分类:;3.1 按时间抽取(DIT)的基-2FFT 算法 —库利-图基算法;;椅糙鉴蚂吉轴惫丢旗桶肇考胁垄约怂嚷洪凋椰坡咯胚竹刨上糠液酬鞋檬尧数字信号处理_9数字信号处理_9;请注意，这里m取值范围！;利用系数周期性;请注意，这里m取值范围！;综合上述，已知：;;例如 N=8 时的DFT，可以分解为两个 N/2=4点的DFT。具体步骤如下:;(2). k为奇数时，即x[1]， x[3]， x[5]， x[7] 分别记作: ;(3).对X1[m]和X2[m]进行蝶形运算，前半部为X[0]~ X[3],后半部分为X[3]~X[7]。; 由于N=2L，所以 N/2仍为偶数，可以进一步把每个N/2点的序列再按其奇偶部分分解为两个N/4的子序列。 例如，k为偶数时的N/2点分解为:;从而可得到：;;例如，N=8时的DFT可分解为四个N/4的DFT, 具体步骤如下:;(2) 将原序列x[k]的“偶中奇”部分:;(3) 将原序列x[k]的“奇中偶”部分:;（5）由 X3[0], X3[1], X4[0], X4[1]进行碟形运算, 得到X1[0], X1[1], X1[2], X1[3]。;（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] 。;对于N=8时DFT，N/4点即为两点DFT,因此;;;基2时间抽取FFT算法流图; 这种FFT算法,是在时间上对输入序列的次序是属于偶数还是属于奇数来进行分解的，所以称作按时间抽取的算法(DIT)。 ;基2时间抽取FFT算法流图;4点基2时间抽取FFT算法流图;4点基2时间抽取FFT算法流图;8点基2时间抽取FFT算法流图;;第一级;二、运算量;x[0] x[4];复数乘法：;;三、算法特点;X2[0]; 所以，当n=1时，则有(前两个蝶形) ; 当n=2时,则有（前两个蝶形）;在某列进行蝶形运算的任意两个节点(行)m和j的节点变量Xn-1[m]，Xn-1[j]就完全可以确定蝶形运算的结果Xn[m]，Xn[j]。 ??其它行(节点)无关。 蝶形运算的两个输出值仍可放回蝶形运算的两个输入所在的存储器中，即实现所谓原位运算。 每一组(列)有N/2个蝶形运算，所以只需N个存储单元，可以节 省存储单元。;2.倒位序规律;0; 3.倒位序实现;0 0 0 0 0 0 0 0;x[0] x[1] x[2] x[3] x[4] x[5] x[6] x[7];FFT算法流图旋转因子 规律;例：试利用N=4基2时间抽取的FFT流图计算8点序列x[k]={1, -1, 1, -1, 2, 1, 1, 2}的DFT。(书上作业3-5);例：试利用N=4基2时间抽取的FFT流图计算8点序列x[k]={1, -1, 1, -1, 2, 1, 1, 2}的DFT。;作业