天津大学测试技术课件-Cha5.pptVIP

在时域可以看到，如果一个周期信号的周期趋于无穷大，则周期信号将演变成一个非周期信号；反过来，如果将任何非周期信号进行周期性延拓，就一定能形成一个周期信号。 我们把非周期信号看成是周期信号在周期趋于无穷大时的极限，从而考查连续时间傅立叶级数在 T趋于无穷大时的变化，就应该能够得到对非周期信号的频域表示方法。 §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= 个； 即前一半的结果。 2.用圆周卷积计算线性卷积 圆周卷积是线性卷积的周期延拓序列的主值序列. 的长度为 , 的长度为 ,先构造长度均为L长的序列, 即将 补零点;然后再对它们进行周期延拓 ，即 所以得到周期卷积： 可见,周期卷积为线性卷积的周期延拓,其周期为L.由于 有 个非零值,所以周期L必须满足: 又由于圆周卷积是周期卷积的主值序列,所以圆周卷积是线性卷积的周期延拓序列的主值序列,即 一个X(k)的值的工作量,如X(1) 得: 对称性: 周期性: 因此， (n为偶数) (n为奇数) 1 2 1 2 k N 同理, 这就是说,X1(k),X2(k)的后一半,分别 等于其前一半的值。 3.X(k)的后一半的确定 由于 （周期性）,所以: 可见,X(k)的后一半，也完全由X1(k), X2 (k)的前一半所确定。 *N点的DFT可由两个N/2点的DFT来计算。 又由于 ,所以 实现上式运算的流图称作蝶形运算 前一半 4.蝶形运算 后一半 （N/2个蝶形) (前一半) (后一半) 1 1 1 1 -1 由X1(k)、X 2(k)表示X(k)的运算是一种特殊的运算-碟形运算 (1)N/2点的DFT运算量:复乘次数: 复加次数: (2)两个N/2点的DFT运算量:复乘次数: 复加次数: (3)N/2个蝶形运算的运算量:复乘次数: 复加次数: 5.计算工作量分析 复乘: 复加: 总共运算量: 按奇、偶分组后的计算量： *但是,N点DFT的复乘为N2 ;复加N(N-1);与 分解后相比可知,计算工作点差不多减少 一半。 例如 N=8 时的DFT,可以分解为两个 N/2=4点的DFT.具体方法如下: (1)n为偶数时,即 分别记作: 四.离散时间、离散频率的傅氏变换--DFT x(nT)=x(n) t 0 T 2T 1 2 N n 0 0 1 2 3 k NT 由上述分析可知，要想在时域和频域都是离散的，那么两域必须是周期的。 时域信号 频域信号 离散的 周期的 周期