按时间抽选的基2 FFT算法课堂.pptVIP

第三节 按时间抽选的基 2-FFT 算法 1 、算法原理 设输入序列长度为 N=2 M (M 为正整数，将该序列 按时间顺序的奇偶分解 为越来越短的子序列，称为 基 2 按时间抽取的 FFT 算法。也称为 Coolkey-Tukey 算法。 其中 基 2 表示： N=2 M ， M 为整数 . 若不满足这个 条件，可以人为地加上若干零值（ 加零补长 ）使其 达到 N=2 M 。 1 先将 x(n) 按 n 的奇偶分为两组，作变量置换 : 当 n= 偶数时，令 n=2r; 当 n= 奇数时，令 n=2r+1; ? 分组，变量置换 2 、算法步骤 1 0 ) ( )] ( [ ) ( 1 0 ? ? ? ? ? ? ? ? N k W n x n x DFT k X N n nk N 得到： 1 ,..., 0 ) ( ) 1 2 ( ) ( ) 2 ( 2 2 1 ? ? ? ? ? ? ? ? N r r x r x r x r x 2 ? 带入 DFT 中 ? ? ? ? ? 1 0 ) ( )] ( [ ) ( N n nk N W n x n x DFT k X ? ? ? ? ? ? ? ? ? ? 1 2 / 0 ) 1 2 ( 1 2 / 0 2 ) 1 2 ( ) 2 ( N r k r N N r rk N W r x W r x ? ? ? ? ? ? ? ? 1 0 1 0 ) ( ) ( N n nk N N n nk N n n W n x W n x 为 奇 数 为 偶 数 ? ? ? ? ? ? ? ? 1 2 / 0 2 2 1 2 / 0 2 1 ) ( ) ( N r rk N k N N r rk N W r x W W r x 3 所以 ? ? ? ? ? ? ? ? 1 2 / 0 2 2 1 2 / 0 2 1 ) ( ) ( ) ( N r rk N k N N r rk N W r x W W r x k X 由于 n N n N j n N j n N W e e W 2 / 2 / 2 2 2 2 ? ? ? ? ? － － ? ? ? ? ? ? ? ? 1 2 / 0 2 / 2 1 2 / 0 2 / 1 ) ( ) ( N r rk N k N N r rk N W r x W W r x ) ( ) ( 2 1 k X W k X k N ? ? 1 , , 1 , 0 2 ? ? N k ? 1 ? N k ? , , 1 , 0 2 ? N ? 4 X 1 (k) 、 X 2 (k) 只有 N/2 个点，以 N/2 为周期；而 X (k) 却有 N 个点，以 N 为周期。要用 X 1 (k) 、 X 2 (k) 表达全部的 X (k) 值，还必须利用 W N 系数的 周期特性 。 rk N k N r N W W 2 / ) 2 / ( 2 / ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 2 / 0 2 / 1 1 2 / 0 ) 2 / ( 2 / 1 1 ) ( ) ( ) 2 / ( N r rk N N r k N r N W r x W r x k N X ? ? ? ? ? ? ? ? ) ( ) 2 / ( ) ( ) 2 / ( 2 2 1 1 k X k N X k X k N X 5 后半部分 前半部分 又考虑到 的对称性： k N W k N k N N N k N N W W W W ? ? ? ? ? 2 / ) 2 / ( ) 2 / ( ) 2 / ( ) 2 / ( 2 ) 2 / ( 1 k N X W k N X k N X k N N ? ? ? ? ? ? 有： 1 , , 1 , 0 ) ( ) ( ) ( 2 2 1 ? ? ? ? N k N k k X W k X k X ? 1 , , 1 , 0 ) ( ) ( 2 2 1 ? ? ? ? N k N k k X W k X ? 6 后半部分 前半部分 ) ( ) ( ) 2 / ( 2 1 k X W k X k N X k N ? ? ? ) ( ) ( ) ( 2 1 k X W k X k X k N ? ? 1 , , 1 , 0 2 ? ? N k ? ) ( 1 k X ) ( 2 k X k N W ) ( ) ( 2 1 k X W k X k N ? ) ( ) ( 2 1 k X W k X k N ? 蝶形运算流图符号 说明： (1) 左边两路为输入 (2) 右边两路为输出 (3) 中间以一个小圆表示加、 减运算（右上路为相加 输出、右下路为相减输 出 ) 1 个蝶形运算需要 1 次复乘， 2 次复加 7 复数乘法 复