数字信号处理--快速傅里叶变换(FFT).ppt

　　形成倒序J后，将原数组A中存放的输入序列重新按倒序排列。 (1) 原输入序列x(I)先按自然顺序存入数组A中。 对N=8，A(0)，A(1)，…，A(7)中依次存放着x(0)，x(1)，x(2)，…，x(7) (2) 对x(n)重新排序(倒序)：规律如图4.2.8所示： ① 第一个序列值x(0)和最后一个序列值x(N－1)不需要重排； ② 每计算出一个倒序值J，便与循环语句自动生成的顺序I 比较：当I=J时，不需要交换数据； 当I≠J时, A(I)与A(J)交换数据。 ③为了避免再次调换前面已调换过的一对数据，只对IJ的 情况调换A(I)和A(J)的内容。 图4.2.8 倒序规律 图4.2.9 倒序程序框图 4.2.5 频域抽取法FFT(DIF-FFT) 　　设序列x(n)长度为N=2M，首先将x(n)前后对半分开，得到两个子序列，其DFT如下： 式中 将X(k)分解成偶数组与奇数组： 当k取偶数(k=2m, m=0, 1, …, N/2－1)时 (4.2.14) 当k取奇数(k=2m+1, m=0, 1, …, N/2－1)时， 令 将x1(n)和x2(n)分别代入(4.2.14)和(4.2.15)式，可得 　　　　　　　　　　　　　　　 (4.2.16) ——X(k)按k值的奇偶分为两组：偶数组是x1(n)的N/2点DFT， 奇数组是x2(n)的N/2点DFT。 x1(n)、x2(n)和x(n)之间的关系也可用蝶形运算流图符号表示 图4.2.11 DIF-FFT第一次分解运算流图（N=8） 第4章 快速傅里叶变换(FFT) 第4章 快速傅里叶变换(FFT) 4.1 引言 4.2 基2FFT算法 4.3 进一步减少运算量的措施 4.4 其他快速算法简介 习题与上机题 4.1 引 言 DFT是数字信号分析与处理中的一种重要变换 直接计算N点DFT：与变换区间长度N的平方成正比 N较大时，计算量太大——难以进行实时处理 快速傅里叶变换FFT(Fast Fourier Transform) 高效计算DFT的算法 使DFT运算效率提高1~2个数量级 4.2 基2FFT算法 4.2.1 直接计算DFT的特点及减少运算量基本途径 一、直接计算DFT　　 有限长序列x(n)的N点DFT为 　　　　　　　　　　　　　　　　　 x(n)为复数序列的一般情况： 计算X(k)的1个值： 需要 次复数乘法和 次复数加法； 计算X(k)的所有N个值： 共需N2次复数乘法和N(N－1)次复数加法运算 （N1时， N(N－1)= N2 ）。 (4.2.1) N (N－1) N较大时，运算量很大： 例如N=1024时，N2=1 048 576 实时处理信号难以实现，对运算速度要求高 二、减少运算量的基本途径： 1、N较大，运算量大 减小N，把N点DFT分解为几个较短的DFT 2、旋转因子具有周期性和对称性 周期性： 　　　　　　　　　　　　　　　 　　 对称性： （4.2.2） 或者 (4.2.3a) (4.2.3b) 2：利用　　的周期性和对称性来减少DFT的运算次数 减少运算量的基本途径： 1：不断地把长序列的DFT分解成几个短序列的DFT 4.2.2 时域抽取法基2FFT基本原理 基2FFT算法： 频域抽取法FFT (DIF－FFT) 时域抽取法： 　　设序列x(n)的长度为N，且满足N=2M，M为自然数。 （若不满足可补零） 按n的奇偶分解x(n) 为两个N/2点的子序列如下： 时域抽取法FFT(DIT－FFT ) 则x(n)的N点DFT为 而 (4.2.4) 其中Ｘ1(k)和X2(k)分别为x1(r)和x2(r)的N/2点DFT： (4.2.5) (4.2.6) 所以 由于X1(k)和X2(k)均以N/2为周期，且 (4.2.7) (4.2.8) 1、X1(k)和X2(k)两个N/2点DFT 2、上式的运算 图4.2.1 蝶形运算符号 因此X(k)又可表示为