数字信号处理实验第一次报告实验三快速傅立叶变换及其应用.pdf

《数字信号处理》 Matlab 实验报告 实验三 快速傅立叶变换及其应用 姓名： 学号： 一． 实验平台 二． 实验目的： （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 可以很快 的计算互相关函数。 四． 实验内容： 实验中用到的函数序列： (a) Gaussian序列 exp(-(n-p).^2)/q), 0=n=15 X (n)= a 0, 其他 （b）衰减正弦序列 1 《数字信号处理》 Matlab 实验报告 exp(-an)*sin(2pi*fn),0=n=15 X(b)= 0, 其他 (c)三角波序列 n, 0=n=3 X b(n)= 8-n, 4=n=7 0, 其他 (d)反三角波序列 4-n, 0=n=3 X c(n)= n-4, 4=n=7 0, 其他 五． 上机实验内容： 1．观察高斯序列的时域