4-快速傅里叶变换分析.ppt

§4-1引言 一.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倍 §4-2 按时间抽取(DIT)的FFT算法 —库利-图基算法 一.算法原理(基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= 个； 即前一半的结果。 §4-6 CZT（线性调频z变换算法） CZT的特点 CZT可计算单位圆上任一段曲线上的Z 变换，可任意给定起止频率； 作变换时输入的点数N和输出点数M可以不相等； 可达到频域“细化”的目的。 CZT的计算： 0 1 1 2 3 2 3 。 。 0 1 2 3 。 0 1 1 2 3 2 3 。 0 1 2 3 4 0 1 2 3 4 5 0 1 1 2 3 2 3 。 。 0 1 2 3 4 5 1 3 3 5 6 6 。 二.用FFT算 的步骤： FFT FFT IFFT x x(n) h(n) y(n) X(k) H(k) Y(k) 流程图 三.几点说明 1. 当 M=L 时,用圆周卷积计算线性 卷积的速度快,点数越多,速度越快, N≥64时,速度增加明显. 2. LM 时,速度增加不太明显,为了 增加速度,采用 (1)重叠相加法 (2) 重叠保留法(从略). 其中： Z在其 ROC 内取值，现为Z指定一离散的路径： Z变换： 做DFT时，Z变换在单位圆上的等分的 N个点上取值。 CZT时，离散路径可在单位圆内、外，或圆上。 CZT在Z平面上的变换 路径是一条螺旋线 决定CZT的起点； 决定变换路径如何倾斜 决定变换的步长。 信号的点数 N 和变换路径的点数 M 可以不相等。 CZT变成了DFT 时，起点在单位圆外， 反之，在圆内； 时，内旋，反之外旋； 时, CZT变换路径 为单位园上一段弧， 所以,当m=1时,则有（前两个蝶形） 当m=2时,则有（前两个蝶形） 当m=3时,则有（前两个蝶形） 可见,在某列进行蝶形运算的任意两个节点(行)k和j的节点变量 就完全可以确定蝶形运算的结果 ,与其它行(节点)无关。 这样,蝶形运算的两个输出值仍可放回蝶形运算的两个输入所在的存储器中,即实现所谓原位运算。每一组(列)有N/2个蝶形运算,所以只需N个存储单元，可以节 省存储单元。 2 倒位序规律 由图可知,输出X(k)按正常顺序排 列在存储单元,而输入是按顺序: 这种顺序称作倒位序,即二进制数 倒位。 n =0 0 n =1 0 n =0 1 n =1 1 n =0 1 n =1 1 0 1 0 1 0 1 0 1 (n2) x(000) 0 乾 x(100) 4 兑 x(010) 2 离 x(110) 6 震 x(001) 1巽 x(101) 5 坎 x(011) 3 艮 x(111) 7 坤 （偶） （奇） 这是由奇偶分组造成的,以N=8为例 说明如下: 3.倒位序实现 输入序列先按自然顺序存入存储