DSP研究性学习分析方案频谱计算96565.docVIP

《近代数字信号处理》课程研究性学习报告 姓名 学号 同组成员 指导教师 时间 DFT近似计算信号频谱专题研讨 【目的】 (1) 掌握利用DFT近似计算不同类型信号频谱的原理和方法。 (2) 理解误差产生的原因及减小误差的方法。 (3) 培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。 【研讨题目】 基本题 M2-0 已知一离散序列为 (1)用L=32点DFT计算该序列的频谱，求出频谱中谱峰的频率； (2)对序列进行补零，然后分别用L=64、128、256、512点DFT计算该序列的频谱，求出频谱中谱峰的频率； (3)讨论所获得的结果，给出你的结论。该结论对序列的频谱计算有何指导意义？ 【题目分析】 本题讨论补零对离散序列频谱计算的影响。 【温磬提示】 在计算离散非周期序列频谱时常用?/?作为横坐标，称?/?为归一化频率?M2-3 已知一离散序列为 x [k]=Acos?0k+Bcos?????0+??)k)。用长度N=64的哈明窗对信号截短后近似计算其频谱。试用不同的A和B的取值，确定用哈明窗能分辩的最小的谱峰间隔中c的值。 【题目分析】 本题讨论用哈明窗计算序列频谱时的频率分辨率。 【仿真结果】 【结果分析】 将实验结果与教材中定义的窗函数的有效宽度相比较，发表你的看法。 【自主学习内容】 【阅读文献】 【发现问题】 (专题研讨或相关知识点学习中发现的问题)： 【问题探究】 在离散序列频谱计算中为何要用窗函数？用不同的窗函数对计算结果有何影响？与矩形窗相比哈明窗有何特点?如何选择窗函数？ 【仿真程序】 M2-4 已知一离散序列为 x[k]=cos(?0k)+0.75cos(?1k), 0( k ( 63 其中?0=0.4?, ?1=?0+?/64 (1) 对x[k]做64点FFT, 画出此时信号的频谱。 (2) 如果(1)中显示的谱不能分辨两个谱峰，是否可对(1)中的64点信号补零而分辨出两个谱峰。通过编程进行证实，M2-6 试用DFT近似计算高斯信号的频谱抽样值。高斯信号频谱的理论值为 通过与理论值比较，讨论信号的时域截取长度和抽样频率对计算误差的影响。 【题目分析】 连续非周期信号频谱计算的基本方法。计算中出现误差的主要原因及减小误差的方法。 【仿真结果】 【结果分析】 由于信号及频谱都有理论表达式，在进行误差分析时希望给出一些定量的结果。 【自主学习内容】 【阅读文献】 【发现问题】 (专题研讨或相关知识点学习中发现的问题)： 【问题探究】 【仿真程序】 扩展题 M2-7 本题研究连续周期信号频谱的近似计算问题。 周期为T0的连续时间周期信号x(t)可用Fourier级数表示为 其中 X(n?0)称为连续时间周期信号x(t)的频谱函数。称为信号的基频（基波），称为信号的谐波。 如果信号x(t)函数表达式已知，则可由积分得出信号的频谱。如果信号x(t)函数表达式未知，或者x(t)函数表达式非常复杂，则很难由积分得信号的频谱。本题的目的就是研究如何利用DFT近似计算连续时间周期信号的频谱。 (1)若在信号x(t)的一个周期T0内抽样N个点，即， T为抽样周期(间隔)，可获得序列x[k] 试分析序列x[k]的DFT与连续时间周期信号x(t)的频谱X(n?0)的关系； (2)由(1)的结论，给出由DFT近似计算周期信号频谱X(n?0)的方案； (3)周期信号x(t)的周期T0=1，x(t)在区间[0，1]的表达式为 x(t)=20t2(1?t)4cos(12?t) (a)试画出信号x(t)在区间[0，1]的波形； (b)若要用10次以内的谐波近似表示x(t)，试给出计算方案，并计算出近似表示的误差。讨论出现误差的原因及减小误差的方法。 【题目分析】 【理论推导】 DFT计算所得结果X[m]与连续周期信号频谱X(n?0)的关系。 【计算方案】 根据理论推导结果设计近似计算方案。分析产生误差的主要原因。 【扩展分析】 如果周期信号x(t)是带限信号