DFT是重要的变换.pptVIP

§5.7 FFT 一.DFT的计算工作量 两者的差别仅在指数的符号和因子1/N. 通常x(n)和 都是复数,所以计算一个 X(k)的值需要N次复数乘法运算,和 次 复数加法运算.那么,所有的X(k)就要N2次复 数乘法运算,N(N-1)次复数加法运算.当N很 大时,运算量将是惊人的,如N=1024,则要完 成1048576 次(一百多万次)运算.这样,难以做到实时处理. 二.改进的途径 1. 的对称性和周期性 利用上述特性,可以将有些项合并,并 将DFT分解为短序列,从而降低运算次数,提 高运算速度.1965年,库利(cooley)和图基 (Tukey)首先提出FFT算法.对于N点DFT,仅需 (N/2)log2N 次复数乘法运算.例如N=1024=210 时， 需要(1024/2)log2 210 =512*10=5120次。 5120/1048576=4.88% ,速度提高20倍 一.算法原理(基2FFT) (一)N/2点DFT 1.先将 按n的奇偶分为两组作DFT,设N=2L ,不足时,可补些零。这样有: n为偶数时: n为奇数时: 由于: 所以,上式可表示为: 其中, 2.两点结论: (1) X (k),X (k)均为N/2点的DFT。 (2) X(k)=X (k)+W X (k)只能确定出 X(k)的k= 个； 即前一半的结果。 (奇中偶) (奇中奇) 同样对n为奇数时，N/2点分为两个N/4点 的序列进行DFT，则有: 例如,N=8时的DFT可分解为四个N/4的DFT, 具体步骤如下: (1) 将原序列x(n)的“偶中偶”部分: 构成N/4点DFT，从而得到X3(0), X3(1)。 7.共轭对称特性之五、六 8.X(k)圆周共轭对称分量与圆周共轭反对称分量的对称性 9.实、虚序列的对称特性 当x(n)为实序列时，根据特性之三，则 X(k)=Xep(k) 又据Xep(k)的对称性： 当x(n)为纯虚序列时，根据特性之四，则 X(k)=Xop(k) 又据Xop(k)的对称性： 四.圆周卷积和 1.时域卷积定理 设 和 均为长度为N的有限长 序列，且 ， 如果 ，则 N N 证明： 相当于将 作周期卷积和后， 再取主值序列。 将 周期延拓： 则有： 在主值区间 ，所以： N 同样可证： N 2.时域圆周卷积过程 N-1 0 n N-1 0 n 0 m 0 m 0 m 0 m 0 2 3 3 2 1 1 N-1 n N 最后结果： 五.有限长序列的线性卷积与圆周卷积 1.线性卷积 的长度为 的长度为 它们线性卷积为 的非零区间为 的非零区间为 两不等式相加得 也就是 不为零的区间. 例如: 1 0 1 2 n 1 0 1 2 n 3 m -1 -2 -3 m m 1 0 1 2 m m n 2 1 0 3 1 4 5 2 3 3 2 1 1 0 1 2 m 2.用圆周卷积计算线性卷积 圆周卷积是线性卷积的周期延拓序列的主值序列. 的长度为 , 的长度为 ,先构造长度均为L长的序列, 即将 补零点;然后再对它们进行周期延拓 ，即 所以得到周期卷积： 可见,周期卷积为线性卷积的周期延拓,其周期为L.由于 有 个非零值,所以周期L必须满足: 又由于圆周卷积是周期卷积的主值序列,所以圆周卷积是线性卷积的周期延拓序列的主值序列,即 一个X(k)的值的工作量,如X(1) 得: 对称性: 周期性: 因此， (n为偶数) (n为奇数) 1 2 1 2 k N 同理, 这就是说,X1(k),X2(k)的