FFT谱分析实验.docVIP

实验报告 课程名称： 数字信号处理 实验题目： 用FFT作谱分析 院 系： 学院 班 级： 姓 名： 学 号： 指导教师： 实验时间： 年 月 XXXX大学 实验一： 用FFT作谱分析 1.实验目的 (1) 进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法， 所以FFT的运算结果必然满足DFT的基本性质。 熟悉FFT算法原理和FFT子程序的应用。 (3) 学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。 实验步骤 复习DFT的定义、 性质和用DFT作谱分析的有关内容。 复习FFT算法原理与编程思想， 并对照DIT-FFT运算流图和程序框图。 编制信号产生子程序， 产生以下典型信号供谱分析用： (4) 编写主程序。 按实验内容要求， 上机实验， 并写出实验报告。 上机实验内容 (1) 画出2 中所给出的信号，并逐个进行谱分析。 下面给出针对各信号的FFT变换区间N以及对连续信号x6(t)的采样频率fs， 供实验时参考。 x1(n), x2(n), x3(n), x4(n), x5(n): N=8,16 x6(t): fs=64(Hz), N=16, 32, 64（n=0:1:69） 作图及谱分析如下： X1(n): 8点 时域 幅频 相频 X2(n): 8点 时域 幅频 相频 X3(n): 8点 时域 幅频 相频 X4(n): 8点 时域 幅频 相频 X5(n): 8点 时域 幅频 相频 X6(n): 8点 时域 幅频 相频 由于16点的与上述方法相同，故画在一张图中： 时域： 幅频： 相频： (2) 令x(n)=x4(n)+x5(n)， 用FFT计算 8 点和 16 点离散傅里叶变换， X(k)=DFT［x(n)］ 结果如下： 8点 时域： 幅频： 相频： 16点 时域： 幅频： 相频： 令x(n)=x4(n)+jx5(n)， 重复(2)。 结果如下： 8点 时域（取模）： 幅频： 相频： 16点 时域（取模）： 幅频： 相频： 由于x6(n)比较特殊，故在此另作一部分画出其图像以及频率特性： 时域波形： 16点 幅频： 相频： 32点 幅频： 相频： 64点 幅频： 相频： 4. 思考题 (1) 在N=8时， x2(n)和x3(n)的幅频特性会相同吗? 为什么? N=16呢? 在N=8时，x2(n)和x3(n)的幅频特性相近，因为此时变换的点数不多，两函数变换后相似。 在N=8时，x2(n)和x3(n)的幅频特性区别开来，因为此时变换点数增加，突出了信号幅频的细节。 (2) 如果周期信号的周期预先不知道， 如何用FFT进行谱分析? 估计信号上限频率，先从较大的采样频率开始，若特性图只有一小部分有，则说明采样频率过大，应该适当减小采样频率。若最高频率有值则说明采样值不够，需增大采样频率。 5. 体会与问题 1.在实验前应当巩固本部分知识，熟悉MATLAB的使用； 2.实验中当点数较少时，变换后特性粗糙；当点数增加则图像变得平滑； 3.FFT是DFT的快速算法，所以其结果也和DFT的变换结果一致，但计算量大大减小。 6. 代码附录 8点变换的所有代码: function [x1,x2,x3,x4,x5,x6]= Signal(N) i=8; x1=zeros(1,i); for n=1:i if (n=1)(n=4) x1(n)=1; else x1(n)=0; end end x2=zeros(1,i); for n=1:i if (n=1)(n=4) x2(n)=n+1-1; else if (n=5n=8) x2(n)=8-(n-1); else x2(n)=0; end end end x3=zeros(1,i); for n