第十三章 数字信号处理实验.docVIP

实验二 [实验目的步骤 x3(n) = x4(n) = cos(πn /4) x5(n) = sin(πn /8) x6(t) = cos8πt + cos16πt + cos20πt 应当注意，如果给出的是连续信号xa(t)，则首先要根据其最高频率确定抽样频率fs以及由频率分辨率选择抽样点数N，然后对其进行软件抽样(即计算 x(n)=xa(nT)，0≤n≤N-1)，产生对应序列 x(n)。对信x6(t)，频率分辨率的选择要以能分辨开其中的三个频率对应的谱线为准则。对周期序列，最好截取周期的整数倍进行谱分析，否则有可能产生较大的分析误差。请实验者根据DFT的隐含周期性思考这个问题。 (4) 编写主程序。图2.1 给出了主程序框图，供参考。 图 2.1 主程序框图 (5) 按实验内容要求，上机实验，并写出实验报告。 [实验内容] 对上面所给出的信号逐个进行谱分析。下面给出针对各信号的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 [思考题] (1) 在N=8时，x2(n)和x3(n)的幅频特性会相同吗？为什么？N=16呢？ (2) 如果周期信号的周期预先不知道，如何用FFT进行谱分析？ [实验报告要求] (1) 简述实验原理及目的； (2) 结合实验中所得给定典型序列幅频特性曲线，与理论结果比较，并分析 说明误差产生的原因以及用FFT作谱分析时有关参数的选择方法； (3) 总结实验所得主要结论； (4) 简要回答思考题。 开始 结束 读入长度N 调用信号产生子程序产生实验信号 调用绘图子程序(函数)绘制时间序列波形图 调用FFT子程序(函数)计算信号的DFT 调用绘图子程序(函数)绘制|X(k)|曲线