数字信号处理--实验五-用DFT(FFT)对信号进行频谱分析.docxVIP

学生实验报告 开课学院及实验室： 电子楼317 2013年4月29日 学院 机械与电气工 程学院 年级、专 业、班 姓名 学号 实验课程名称 数字信号处理实验 成绩 实验项目名称 实验五用DFT( FFT)对信号进行频谱分析 指导 老师 、实验目的 学习DFT的基本性质及对时域离散信号进行频谱分析的方法，进一步加深对频域概念和数字频 率的理解，掌握 MATLAB函数中FFT函数的应用。 二、实验原理 离散傅里叶变换(DFT)对有限长时域离散信号的频谱进行等间隔采样，频域函数被离散化了, 便于信号的计算机处理。 设x(n)是一个长度为 M的有限长序列，x(n)的N点傅立叶变换: X(k) N 1 j三 kn DFT[x(n)]N x(n)e N 0 k N 1 n 0 其中WN e .2 j N ，它的反变换定义为: 1 X(n) Nk N 1 nk X(k)WN 0 令z WNk， X(z z WNk 则有: N 1 x( n)Wjk n 0 可以得到，X(k) X(Z)Z WNk Z Wn *是Z平面单位圆上幅角为 2 k N 的点，就是将单 位圆进行N等分以后第 K个点。所以, X(K)是Z变换在单位圆上的等距采样，或者说是序列傅立叶 变换的等距采样。时域采样在满足 Nyquist定理时，就不会发生频谱混叠。 DFT是对序列傅立叶变 换的等距采样，因此可以用于序列的频谱分析。 如果用FFT对模拟信号进行谱分析，首先要把模拟信号转换成数字信号，转换时要求知道模拟 信号的最高截至频率，以便选择满足采样定理的采样频率。 般选择采样频率是模拟信号中最高频率 的3~4倍。另外要选择对模拟信号的观测时间，如果采样频率和观测时间确定，则采样点数也确定 了。这里观测时间和对模拟信号进行谱分析的分辨率有关, 最小的观测时间和分辨率成倒数关系。 最 小的采样点数用教材相关公式确定。要求选择的采样点数和观测时间大于它的最小值。 如果要进行谱分析的模拟信号是周期信号, 最好选择观测时间是信号周期的整数倍。 如果不知道 信号的周期，要尽量选择观测时间长一些，以减少阶段效应的影响。 用FFT对模拟信号作谱分析是一种近似的谱分析，首先一般模拟信号（除周期信号以外）的频 谱是连续谱，而用 FFT作谱分析得到的是数字谱，因此应该取 FFT的点数多一些，用它的包络作为 模拟信号的近似谱。另外，如果模拟信号不是严格的带限信号， 会因为频谱混叠现象引起谱分析的误 差，这种情况下可以预先将模拟信号进行滤波，或者尽量采样频率取高一些。 一般频率混叠发生在折叠频率附近， 分析时要注意因频率混叠引起的误差。 最后要注意一般模拟 信号是无限长的，分析时要截断，截断的长度和分辨率有关， 但也要尽量取长一些，取得太短会截断 引起的误差会很大。举一个极端的例子，一个周期性正弦波，如果所取观察时间太短，例如取小于一 个周期，它的波形和正弦波相差太大，肯定误差很大，但如果取得长一些，即使不是周期的倍数，这 种截断效应也会小一些。 如同理论课教材所讨论的，在运用 DFT进行频谱分析的时候可能有三种误差，即： （1 ）混叠现象 当采样率不满足 Nyquist定理，经过采样就会发生频谱混叠。这导致采样后的信号序列频谱不 能真实的反映原信号的频谱。所以，在利用 DFT分析连续信号频谱的时候，必须注意这一问题。避 免混叠现象的唯一方法是保证采样的速率足够高， 使频谱交叠的现象不出现。 这告诉我们，在确定信 号的采样频率之前，需要对频谱的性质有所了解。在一般的情况下，为了保证高于折叠频率的分量不 会出现，在采样之前，先用低通模拟滤波器对信号进行滤波。 （2 ）泄漏现象 实际中的信号序列往往很长，甚至是无限长。为了方便，我们往往用截短的序列来近似它们。 这样可以使用较短的 DFT来对信号进行频谱分析。这种截短等价于给原信号序列乘以一个矩形窗函 数。值得一提的是，泄漏是不能和混叠完全分离开的，因为泄漏导致频谱的扩展，从而造成混叠。为 了减少泄漏的影响，可以选择适当的窗函数使频谱的扩散减到最小。 （3 ）栅栏效应 因为DFT是对单位圆上Z变换的均匀采样，所以它不可能将频谱视为一个连续函数。 这样就产 生了栅栏效应，从某种角度看， 用DFT来观看频谱就好像通过一个栅栏来观看一幅景象， 只能在离 散点上看到真实的频谱。这样的话就会有一些频谱的峰点或谷点被 “栅栏”挡住，不能被我们观察到。 减小栅栏效应的一个方法是在源序列的末端补一些零值，从而变动 DFT的点数。这种方法的实质是 改变了真是频谱采样的点数和位置， 相当于搬动了 “栅栏”的位置，从而使得原来被挡住的一些频谱 的峰点或谷点显露出来。注意，这时候每根谱线所对应的频和原来的已经不相同了。 从上面的分析过程可以看出， DFT可以用于信号的频谱分析，