《数字信号处理教程——MATLAB释义与实现》第四章分析.ppt

第四章 信号频谱的高效计算 信号频谱的高效计算 4.1 各种傅立叶变换及其关系 4.2 快速付立叶变换(FFT) 4.3 用FFT计算序列的频谱 4.4 时域采样中的频谱变换 4.5 连续信号的频谱计算 4.6 用反变换从频谱计算信号 4.7 用FFT计算能量 4.8 小结 4.1 各种傅立叶变换及其关系 （1）时域的周期性对应于频域的离散化。 （2）?时域的离散化对应于频域的周期性。其主周期就是乃奎斯特频率范围[-π/T,π/T]。 （3） 周期性的离散序列将对应于离散并周期性的频谱。即离散傅立叶级数(DFS)。离散有利于数值计算，但信号无穷延伸又不利于计算。 （4）把时域和频域数据长度都限于主周期，并且使之相等，形成离散傅立叶变换(DFT)。它既离散，又长度有限，适合于计算机数值计算。 各种傅立叶变换的特点 各种傅立叶变换及其相互关系 对离散傅立叶变换(DFT)，人们开发了可以高效地进行计算的方法，称为快速傅立叶变换(FFT)。 人们想尽量利用FFT来解决其他类型信号的频谱计算问题。所以要充分弄清各种傅立叶变换之间的关系。本章就讨论这个主题。先把主要结果列出，其中有些结论已经讨论过，与采样定理有关的结论将在后面讨论，读者可先接受下来。目的是走通下图的路线，完成时频域傅立叶变换的数值计算。 各种傅立叶变换及其相互关系 模拟信号 时域采样 周期延拓 主值区间 时域xa(t) —— x(n) —— ——x(n)·RN(n) | | (FFT) 频域Xa(jΩ)——X(jΩ)—— ——X(k)·RN(k) CTFT DTFT DFS DFT ICTFT IDTFT 频域采样 IDFT 各种傅立叶变换及其相互关系 (1). DFT与DFS的关系： 或 (4.1.7) (2). DFT与DTFT的关系 : 采样 插补 各种傅立叶变换及其相互关系 (3). DTFT与CTFT的关系:用采样定理或乃奎斯特定理建立关系。由CTFT求DTFT 由DTFT求CTFT:频率泄漏可以忽略不计时近似解。 各种傅立叶变换及其相互关系 DFT与CTFT的相互关系 : 由CTFT求DFT只是取出其一部分，没有误差问题。而由DFT求CTFT要先经过由DFT求DTFT，再经过由DTFT求CTFT的两步。由离散求连续不一定有解，只有在信号的频谱足够窄，频率泄漏可以忽略不计的条件才可能近似由下式求得。 4.2 快速付立叶变换(FFT) 计算DFT的运算次数按N2快速增长。设N可以被2整除，把x(n)分成两个子序列x1(n)和x2(n)， 则原序列的DFT可写成(设N1=N/2)： 快速付立叶变换(FFT) 设它们的傅立叶变换分别为X1(m)和X2(m)，其周期为N1=N/2： 则x(n)的傅立叶变换X(m)可表为： 如果N=8,则N1=4, 故X(m)周期是8 ,而X1(m)和X2(m)周期是4, 即X1(4)=X1(0), X1(5)=X1(1),… 。依次取m=0,1,…7，上式对应于右方的运算图。 快速付立叶变换(FFT) X1(m)和X2(m)是N1点的DFT，它们的计算又可以用类似方法化为两个更短的N2=N1/2点的DFT,一直分解下去，直到2为止，这就构成了上述FFT的全部运算图。 粗算其中每一根线条代表一次乘法，线条的合成点代表一次加法。则每一级要N次乘法和加法。N=8时，需log28=3级，故共要24次乘法和加法。原来要N2=64次。若N=1024，需要10242次乘法，而用FFT, 需分解为log21024=9级，只需1024×9次乘法，加快了100多倍。 快速付立叶变换(FFT) 把上述运算次数的估计精确化： (1). 每次分解的乘法次数为N，共log2N次， 乘法次数＝N·log2N (2). 考虑到蝶形运算，又把乘法次数减半。 在FFT运算图中，基本单元为如右图所示的蝶形结构，实际上不需要 四次乘法，而只需 要两次乘法即可完 成。 蝶形运算节省一半乘法 考虑到 在算 时，把m和N/2+m两点成组来进行，即构成上图的蝶形运算，就节省一半乘法 快速付立叶变换(FFT) 用这些措施后总的乘法次数约为（当N很大时）： 当N=1024时，结果为5120，与1