应用快速傅里叶变换对信号进行频谱分析实验报告.pdfVIP

应用快速傅里叶变换对信号进行频谱分析 2.1 实验目的 1、通过本实验，进一步加深对DFT 算法原理和基本性质的理解，熟悉 FFT 算 法原理和 FFT 子程序应用 2、掌握应用FFT 对信号进行频谱分析的方法。 3、通过本次实验进一步掌握频域采样定理。 4、了解应用FFT 进行信号频谱分析过程中可能出现的问题，以便在实际中正。 确应用 FFT。 2.2 实验原理与方法 对于有限长序列我们可以使用离散傅里叶变换（DFT）。这一变换不但 可以好地反映序列的频域特性，而且易于用快速傅里叶变换在计算机上实 现 当序列 x(n)的长度为 N 时，它的离散傅里叶变换为： N 1 X (k )  DF T [x (n)]  x (n)Wkn N n0 其中W  e j (2 / N ) ，它的反变换定义为： N 1 N 1 kn x (n)  IDF T [X (k )]  X (k )W N N k 0 比较 Z 变换公式，令z  Wk 则 N N 1 X (z ) |  k  x (n)Wnk  DF T [x (n)] z W N N n0 因此有X (k )  X (z ) |  k 。 z W N 所以 X(k)是 x(n)的 Z 变换在单位圆上的等距采样，或者说是序列傅 里叶变换的等距采样。 DFT 是对序列傅里叶变换的等距采样，因此可以用于对序列的频谱分 析。在运用 DFT 进行频谱分析的过程中有可能产生三种误差： 1、混叠现象 序列的频谱是原模拟信号频谱的周期延拓，周期为2 / T 。因此，当采样 频率小于两倍信号的最大频率时，经过采样就会发生频谱混叠，使采样后 的信号序列频谱不能真实反映原信号的频谱。 2、泄漏现象 实际中信号序列往往很长，常用截短的序列来近似它们，这样可以用较短 的 DFT 对信号进行频谱分析，这种截短等价于给原信号序列乘以一个矩形 函数。这样得到的频谱会将原频谱扩展开。 3、栅栏效应 DFT 是对单位圆上 Z 变换的均匀采样，所以它不可能将频谱视为一个连续 函数。因此，只能看到离散点无法看出真实频谱。 2.3 实验内容及步骤 1、观察高斯序列的时域和幅频特性 (1)固定信号x (n)  e[(np )2 / q ] ,0  n  N  1 中参数 p=16，信号长度N=32，改变q a 的值，分别等于 2,10,30，观察他们的时域和幅频特性（作 FFT 时，点数 大于 256），了解当 q 取不同的值时，对信号序列的时域和幅频特性的影响。 MATLAB 程序： N=32;FFTN=256;p=16;q1=2;q2=