北京交通大学dsp研究性学习报告频谱计算.pdf

《数字信号处理》课程研究性学习报告 DFT 近似计算信号频谱专题研讨 姓 名 李 帆 学 号 同组成员 张 静 林 恒 王亚君 李亚伟 指导教师 薛 健 时 间 2013 年 5 月 8 日 利用 DFT 近似计算信号频谱专题研讨 【目的】 (1) 掌握利用 DFT 近似计算不同类型信号频谱的原理和方法； (2) 理解误差产生的原因及减小误差的方法； (3) 研究用 DFT 近似计算连续周期信号的方法； (4) 培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。 【研讨内容 】 基本题 基本题是课程的基本要求，所有的人都需完成。 问题一 已知某离散序列为 x[k ] sin( 0.2 πk), k 0,1, ,31 (1) 用 L =32 点 DFT 计算该序列的频谱，求出频谱中谱峰的频率； (2)对序列进行补零，然后分别用 L=64 、128、256 、512 点 DFT 计算该序列的频谱，求出频谱中谱峰 的频率； (3)讨论所获得的结果，从中你能得到了什么结论？该结论对序列的频谱计算有何指导意义？ 【题目分析 】 本题讨论补零对离散序列频谱计算的影响。 补零可以使 DFT计算得出的频谱更加细致，但是不能改变序列的 DTFT 【温磬提示】 在计算离散非周期序列频谱时常用 / 作为横坐标，称 / 为归一化频率 normalized frequency) 。 在画频谱时需给出横坐标。每幅图下都需给出简要的文字说明。 由于离散非周期序列频谱是周期的，所以在计算时不必用 fftshift 函数对 fft 计算的结果进行重 新排列。 【序列频谱计算的基本方法 】 在 MATLAB 中，用函数 fft （x ， N）可以计算 X[k] 序列的 N 点 DFT 【仿真结果 】 L=32 20 10 0 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 L=64 20 10 0 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 L=128 20 10 0 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 L=256 20 10 0 0 0.2 0.4 0.6 0.8 1 1.2