快速傅里叶变换及其应用.pdfVIP

实验三 快速傅里叶变换及其应用 陈郁蕾 一． 实验目的 1. 在理论学习的基础上，通过本实验，加深对FFT 的理解，熟悉MATLAB 中的有关函数。 2. 应用FFT 对典型信号进行频谱分析。 3. 了解应用FFT 进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用FFT。 4. 应用FFT 实现序列的线性卷积和相关。 二． 实验原理 1. 混叠 采样序列的频谱是被采样信号频谱的周期延拓，当采样频率不满足奈奎斯特采样定 理的时候，就会发生混叠，使得刺痒后的序列信号的频谱不能真实的反映原采样信号的 频谱。 2. 泄漏 根据理论分析，一个时间的信号其频带宽度为无限，一个时间无限的信号其频带宽 度则为有限。因此对一个时间有限的信号，应用DFT 进行分析，频谱混叠难以避免。对 一个时间无限的信号虽然频带有限，但在实际运算中，时间总是取有限值，在将信号截 断的过程中，出现了展谱线的现象，称之为频谱泄漏或功率泄漏。 3. 栅栏效应 DFT 是对单位圆上Z 变换的均匀采样，所以它不可能将频谱视为一个连续函数，就 在一定意义上看，用DFT 来观察频谱就好象通过一个栅栏来观看一个景象一样，只能在 离散点上看到真实的频谱，这样就有可能发生一些频谱的峰点和谷点被“尖桩的栅栏” 所挡住，不能被我们观察到。 4. 圆周卷积 把序列X （N）分布在N 等份的圆周上，而序列Y （N）经反摺后也分布在另一个具 有 N 等份的同心圆的圆周上。两圆上对应的数两量两相乘求和，就得到全部卷积序列。 这个卷积过程称做圆周卷积。 5. 互相关函数 反映了两个序列X （N）和Y （N）的相似程度，用FFT 可以很快的计算互相关函数。 三． 实验内容及结果 1. 观察高斯序列的时域和幅频特性，固定信号xa(n) 中参数p=8，改变q 的 值，使 q 分别等于2，4 ，8，观察他们的时域和幅频特性，了解当q 取不同值时，对信号序列 的时域和幅频特性影响；改变p，使p 分别等于8，13，14，观察参数p 变化对信号序 列的时域和幅频特性影响，注意p 等于多少时，会发生明显的泄漏现象，混叠是否也随 之出现？记录实验中观察到的现象，绘出相应的时域序列和幅频特性曲线。 1.1 幅频特性曲线 1.1.1 p=8,q=2 1.1.2 p=8,q=4 1.1.3 p=8,q=8 1.1.4 p=13,q=8 1.1.5 p=14,q=8 1.2 结果分析 X （n）中的参数p 为高斯序列的峰值位置，q 则表示高斯序列峰的尖锐度，（即峰值 边沿的陡峭度）q 值越大，时域图中图象越平缓，序列变化越慢；其幅频特性图中高频 分量越少，频谱越窄，越不容易产生混叠。p 值越大，序列右移，在规定的窗口内有效 值被截断的越多。因为窗口截断会造成窗口泄漏，所以我们可以在幅频特性图中看到， 随着p 值的变大，高频分量会增加。易出现泄露，当p=13 时，特别是p=14 时，产生了 明显的泄漏与混叠。 2. 观察衰减正弦序列xb(n) 的时域和幅频特性，a=0.1,f=0.0625,检查谱峰出现位置是否正确， 注意频谱的形状，绘出幅频特性曲线，改变f ，使f 分别等于0.4375 和0.5625，观察这 两种情况下，频谱的形状和谱峰出现的位置，有无混叠和泄漏现象？说明产生现象的原 因。 2.1 幅频特性曲线 2.1.1 f=0.0625 2.1.2 f=0.4375 2.1.3 f=0. 5625 2.2 结果分析 图中的幅频特性图，当f=0.0625 时，没有产生明显的混叠和泄露，；当f=0.4375 和 f=0.5625 时，产生了混叠，是因为不满足奈奎斯特