第三章 离散傅里叶变换和其快速算法.ppt

计算 的信号流图如下：类似的，可得：注： 均是 点的DFT 上面时间抽选FFT流图具有三个特点： ①基本计算单元为一碟形； ②输入为“混序”排列；输出为正序排列； ③具有“同址计算”特性。 1.碟形计算 对于任意 ，总可以通过M级分解成2点DFT计算，每次由 个碟形计算而成。如下图所示，计算方程为： 一个 的DFT流程图如下图所示: 例：当 时， 为： 如上图红线所示。 谢 谢！ 3.6.2 （Decimation-In-Time FFT－DIT FFT） 基本出发点：利用 的周期性和对称性，将DFT的计算分解成一些逐次减小的DFT计算； 分解规则：（1）对时间进行偶奇分， （2）对频率进行前后分。 M 设 N ? 2 M ? log N 。 2 为讨论方便，以时间信号序列为例： x ( n ) ? { x ( 0 ), x ( 1 ), x ( 2 ), x ( 3 ), x ( 4 ), x ( 5 ), x ( 6 ), x ( 7 )} …(3.6.2a) 由规则（1），将 x(n)按序号偶奇分，得： { x ( 0 ), x ( 2 ), x ( 4 ), x ( 6 ) | x ( 1 ), x ( 3 ), x ( 5 ), x ( 7 )} 即： g ( r ) ? x ( 2 r ), even N , r ? 0 , 1 , 2 ,..., - 1 h ( r ) ? x ( 2 r ? 1 ), odd 2 …(3. 6.2b) 时间抽取FFT算法 M：正整数,则： 由DFT定义得： N - 1 kn kn kn X ( k ) ? x ( n ) W ? x ( n ) W ? x ( n ) W N N N n ? 0 n ? even n ? odd N N - 1 - 1 2 2 2 rk ( 2 r ? 1 ) k ? x ( 2 r ) W ? x ( 2 r ? 1 ) W N N r ? 0 r ? 0 N N - 1 - 1 2 2 2 kr k 2 kr ? g ( r ) ( W ) ? W h ( r ) ( W ) N N N r ? 0 r ? 0 N N - 1 - 1 2 2 rk k rk ? g ( r ) W ? W h ( r ) W N N N 2 2 r ? 0 r ? 0 k ? G ( k ) ? W H ( k ), k ? 0 , 1 , 2 ..., N - 1 …（3.6.2c） N 由规则（2），将 X (k )分为前后两组： X ( k ) ? { X ( 0 ), X ( 1 ), X ( 2 ), X ( 3 ) | X ( 4 ), X ( 5 ), X ( 6 ), X ( 7 )} 由式（3.6.2c），前4个k值的 X (k ) 为： N k X ( k ) ? G ( k ) ? W H ( k ), k ? 0 , 1 , 2 ,... - 1 …（3.6.2d） N 2 后4个k值的 X (k ) 为： N k ? N N N 2 X ( k ? ) ? G ( k ? ) ? W H ( k ? ) …（3.6.2e） N 2 2 2 N k ? G ( k ) - W H ( k ), k ? 0 , 1 , 2 ...