基于DSP的数字信号实验二.docVIP

实验二 用FFT对信号进行频谱分析 一、实验目的 1、加深对离散信号DFT的理解； 2、掌握FFT算法的流程及其MATLAB实现； 3、利用FFT对典型信号进行频谱分析； 4、结合理论知识，对频谱分析中出现的有关现象进行理论分析。 二、实验原理 采样序列的频谱是被采样模拟信号频谱的周期延拓，当采样频率不满足奈奎斯特采样定理时，就会发生频谱的混叠。解决混叠问题的唯一办法是保证采样频率足够高，以防止频谱混叠。 对一个时间无限的信号，虽然频带有限，但在实际DFT运算中，时间长度总是取有限值，在将信号截短的过程中，出现了分散的扩展谱线的现象，称为频谱泄漏或功率泄漏。泄漏也会引起混叠。 从某种意义上来讲，用DFT来观察频谱就如同通过一个栅栏来观看景象一样，只能在离散点上看到真实的频谱，这样一些频谱的峰点或谷点就可能被“尖桩的栅栏”挡住，也就是正好落在两个离散采样点之间，不能被观察到。减小栅栏效应的一个方法是在原序列的末端填补一些零值，从而变动DFT的点数。 三、实验内容 1、对高斯序列xa(n)，令p=8， q分别等于2、4和8，观察q值的改变对高斯序列时域特性和幅频特性的影响，并给出理论解释： 2、对衰减正弦序列xb(n)，a=0.1，f=0.0625，观察其时域和幅频特性，检查谱峰位置是否正确；改变f，使f分别等于0.4375和0.5625，比较这两种情况下频谱形状和谱峰位置，并给出理论解释： 3、观察三角波序列xc(n)和反三角波序列xd(n)的时域特性和幅频特性，用8点FFT分析幅频特性，观察两者的序列形状和频谱曲线有什么异同？为什么？在xc(n)和xd(n)的末尾补零，用32点FFT分析幅频特性，观察幅频特性发生了什么变化？这些变化说明了什么？ 4、对连续单频周期信号，按采样频率进行采样，截取长度分别选N=16和N＝20，比较其幅频特性，并给出理论解释；当截取长度选N＝164时，与N＝20比较，频谱发生了哪些变化？为什么？ 5、一个连续信号含两个频率分量，经采样得： 令N=16，Δf分别为1/16和1/64，观察其频谱，两个谱峰是否可以区分开？为什么？当N=128时，Δf不变，结果有何不同？为什么？ 四、实验结果与分析 1、因为 ; p=8 程序见附录程序一，由matlab作图可得下图，其中左边是时域特性，右边是频域特性 分析：由高斯序列表达式n=p为其对称轴，当p取固定值时，时域图都关于n=8对称截取长度为周期的整数倍，没有发生明显的泄露现象，但存在混叠，当q由2增加至8的过程中，时域图形变化越来越平缓，中间包络越来越大，可能函数周期开始增加，频率降低，渐渐小于fs/2，混叠减弱。 2、程序见附录程序二，下图中，左边为时域特性，右边为频域特性 分析： 当f= 0.0625时，谱峰位置出现正确，存在混叠现象，时域采样为一周期，不满足采样定理。 当f=0.4375和0.5625时，时域图像关于Y轴对称，频域完全相同。这是因为频域图是取绝对值的结果，所以完全相同。另外由于时域采样为6个半周期，满足采样定理，无混叠；但由于截取长度不是周期的整数倍，出现泄漏。 3、程序见附录程序三 N=8时的时域和幅度频谱图： 分析： 由图知，三角波序列和反三角波序列的时域图像成镜像关系，但频域图像完全一样，只是因为幅频图是对X(k） 分析： 由实验所得的图形知，N=32点时，Xc(n)和Xd(n)的幅频特性都更加密集，更多离散点的幅值显示，“栅栏效应”减小，分辨率提高，而对于Xd(n)来说变化更加明显。在原序列的末端填补零值，变动了DFT的点数，人为的改变了对真实频谱采样的点数和位置，相当于搬动了“尖桩栅栏”的位置，从而使得频谱的峰点和谷点暴露出来。N=32时，Xc(n)和Xd(n)的频谱差别较大，但总体趋势仍然都是中间的最小，两侧呈对称。 4、先求离散序列： 其对应程序见附录程序四，所得图像为： 分析：因为所选取的采样频率fs=8fa，82，满足奈奎斯特采样频率，由此可以得到单一谱线的DFT结果，N=16与N=20之间的频谱差距较大，是因为N=20所选的时间长度不是周期的整数倍。而当N=164时，其图形在形状上与N=20时差不多，但是在X（k） 分析： 由图可以看出N=16时，当△f由1/16减小为1/64时，频谱图出现失真，可能是△f的改变引起的周期变化导致混叠。 当N增加至128时，频谱更加密集，分辨率明显提高，混叠现象消失。 五、实验小结 此次实验主要研究利用DFT做连续频谱的分析，并验证了其中的混叠、泄漏、栅栏效应等现象的产生、减小措施等，这次的实验明显比第一次要难，做起来也没有那么顺手，但是经过一遍遍的看、改，也总算是勉强完成了，准确率不知道怎么样，但收获却是颇丰，实验将课本知识与实践中的现象联合在了一起，显得更加形