[物理]第4章 离散傅里叶变换.pptVIP

* * 令 ，则 上式表明，要计算x(n)在半径为r的圆上的等间隔频谱分量，可以先对x(n)乘以r-n，再计算N点DFT(FFT)即可得到。若要求x(n)分布在该圆的有限角度内N点等间隔频谱分量，可以通过尾部补零的方法，仍按式(3-75)用DFT分析整个圆上的等间隔频谱，最后只取所需角度内的频谱分量即可。 显然这种方法的计算量大，下面要介绍的Chirp-Z变换可使这种谱分析的运算量大大减少。 3-75 * * 3. ChirpZ变换 设序列x(n)长度为N，要分析z平面上M点频谱采样值，分析点为zk， k=0, 1, 2, …, M-1。 则 zk=AW-k, 0≤k≤M-1 式中A和W为复数， 用极坐标形式表示为 式中A0和W0为实数。 当k=0时有 3-76 * 西安建筑科技大学信于与控制学院 * 将zk代入Z变换公式得到 利用下面的关系式： 得到: 令 * * 图 3-17 Chrip-Z变换分析频率点分布图 3-77 * * 图 3-18 Chirp z变换计算框图 * * 图 3-19 Chirp-Z变换中hL(n)序列的形成 * * 由(3-66)式知，y(n)×h(n)是V(n)的周期延拓序列的主值序列， 延拓周期为L， 即 综上所述， 可归纳出具体计算步骤如下： (1) 形成hL(n)序列 (2) (3) * * (4) (5) 计算 (6) (7) * * 与标准DFT(FFT)算法相比较， Chirp-Z变换有以下特点： (1) 输入序列长度N和输出序列长度不需要相等， 且二者均可以素数。 (2) 分析频率点zk的起始点z0及相邻两点的夹角φ0是任意的(即频率分辨率是任意的)，因此可从任意频率上开始，对输入数据进行窄带高分辨率的谱分析。 (3) 谱分析路径可以是螺旋形的。 (4) 当A=1，M=N， 时， zk均匀分布在单位圆上，此时Chirp-Z变换就是序列的DFT。 * * 4. 用DFT进行谱分析的误差问题 DFT(实际中用FFT计算)可用来对连续信号和数字信号进行谱分析。 (1) 混叠现象。 (2) 栅栏效应。 (3) 截断效应。 根据傅里叶变换的频域卷积定理有 * * 幅度谱|RN(ω)|~ω曲线如图3-20所示(|RN(ω)|以2π为周期，只画低频部分)。图中|ω|2π/N的部分称为主瓣，其余部分称为旁瓣。 例如， x(n)=cos(ω0n), ω0=π/4其频谱为 其中 * * 图 3-20 矩形窗函数的幅度谱 * * 图 3-21 加矩形窗前后的频谱 * * 令 则 3-60 3-61 式(3-61)称为用X(k)表示X(z)的内插公式，φk(z)称为内插函数。当z=ejω时，(3-60)式和(3-61)式就成为x(n)的傅里叶变换变换X(ejω)的内插函数和内插公式，即 * * 进一步化简可得 在数字滤波器的结构与设计中，我们将会看到，频域采样理论及有关公式可提供一种有用的滤波器结构和滤波器设计途径。 3-62 3-63 * * 3.4 DFT的应用举例 DFT的快速算法FFT的出现，使DFT在数字通信、语音信号处理、图象处理、功率谱估计、仿真、系统分析、雷达理论、光学、医学、地震以及数值分析等各个领域都得到广泛应用。 然而，各种应用一般都以卷积和相关运算的具体处理为依据，或者以DFT作为连续傅里叶变换的近似为基础。 本节主要介绍用DFT计算卷积和相关系数的基础原理以及用DFT对连续信号和序列进行谱分析等最基本的应用。只要掌握了这两种基本应用的原理，就为用DFT解决数字滤波和系统分析等问题打下了基础。 * * 3.4.1用DFT计算线性卷积 如果 且 则由时域循环卷积定理有 0≤k≤L-1 由此可见，循环卷积既可在时域直接计算，也可以按照图3-10所示的计算框图，在频域计算。由于DFT有快速算法FFT，当N很大时，在频域计算的速度快得多，因而常用DFT(FFT)计算循环卷积。 * * 图3-10 用DFT计算循环卷积 在实际应用中，为了分析时域离散线性非时变系统或者对序列进行滤波处理等，需要计算两个序列的线性卷积。和计算循环卷积一样，为了提高运算速度，也希望用DFT(FFT)计算线性卷积。而DFT只能用来计算循环卷积，为此导出线性卷积和循环卷积之间的关系以及循环卷积和线性卷积相等的条件。 * * 假设h(n)和x(n)都是有限长序列，长度分别是N和M。它们的线性卷