实验4-fft频谱分析.pdfVIP

实验 4 FFT 频谱分析 一、实验目的 应用离散傅里叶变换的快速算法 FFT 分析信号的频谱。深刻理解应用 FFT 分析离散和 连续信号的原理，掌握分析过程中出现的现象及解决方法。 二、实验原理 1. 离散周期信号频谱的分析 ~ ~ 周期为 N 的离散信号(序列) x (k ) 的频谱函数X (m) 定义为： 式中：N 是序列的周期；k 为离散的时间变量；m 为离散的频率变量；2 m / N 是次谐波的 数字频率。 ~ 离散周期信号的频谱 X (m) 也是周期为 N 的离散谱，谱线间隔为 2 / N 。利用 MATLAB 提供的 fft 函数可以计算离散周期信号的频谱。对于离散周期序列，只需对周期序 列上一个周期内的数值x (k ) 进行N 点的 FFT 运算，就可准确地得到其频谱在一个周期上的 ~ N 个数值X (m) ，且有X (m)  X (m) ，m [0, N  1] 。其分析步骤如下。 ~ (1) 确定离散周期序列x (k ) 的基本周期为N 。 ~ (2) 利用 fft 函数对序列x (k ) 一个周期内进行N 点FFT 计算，得到X (m) 。 ~ (3) X (m)  X (m) 。 2. 离散非周期信号频谱的分析  离散非周期信号x (k ) 的频谱函数X (ej  ) 为：X (ej  )   x (k )ej k 。 k  利用 MATLAB 提供的 fft 函数可以计算离散非周期信号的频谱。当序列长度有限时， 可以求得准确的序列频谱X (ej  ) 的样点值。若序列很长或无限长时，则由于截短产生泄漏 误差，计算的结果只能是序列频谱X (ej  ) 样点值的近似。分析步骤如下。 (1) 确定序列x (k ) 的长度 M 及窗函数的类型。当序列为无限长时，需要根据能量分布， 利用窗函数进行截短。 (2) 确定做 FFT 的点数N 。根据频域采样定理，为使时域波形不产生混叠，必须取N ≥ M 。 (3) 使用 fft 函数做N 点 FFT 以计算X (m) 。 3. 连续周期信号频谱的分析 ~ 周期为T 的连续时间信号x (t) 的频谱函数X (n ) 定义为： 0 0 1 T ~ 0 j n t X (n )  x (t)e 0 dt 0 T 0 0 2 式中：T 是信号的周期；   2 f 称为信号的基频(基波) ；n 称为信号的谐频。 0 0 0 0 T 0 连续周期信号的