数字信号处理实验指导书(更新).docVIP

《数字信号处理》 实 验 指 导 书 吴启宏 编 华北电力大学 2007 年 3 月 前 言 实验总体目标 通过实验，学生应能深入了解“数字信号处理”这门课程的知识点，加深对基本概念的理解，同时基本掌握MATLAB语言，为今后进一步学习和研究打下一个良好的基础。 ⒉ 适用专业 全校工科各专业 ⒊ 先修课程 信号与系统 ⒋ 实验课时分配 实验项目 学时 实验一 熟悉MATLAB环境 2 实验二 快速傅里叶变换（FFT）及其应用 2 实验三 IIR数字滤波器的设计 2 实验四 FIR数字滤波器的设计 2 ⒌ 实验环境 普通配置计算机，安装MATLAB6.5或以上版本。 ⒍ 实验总体要求 自行编写MATLAB程序，完成上机实验内容，并提交实验报告。 ⒎ 本实验的重点、难点及教学方法建议 本实验的重点为“数字信号处理”的基本理论及应用；难点在MATLAB语言的熟练使用。 教学方法建议：学生在实验前先自学MATLAB语言，指导教师将实验中用到的主要的MATLAB函数列出，提前给学生，供参考。 目 录 实验一、熟悉MATLAB环境 …………………………………… 3 实验二、快速傅里叶变换（FFT）及其应用 …………………………………… 5 实验三、IIR数字滤波器的设计 …………………………………… 10 实验四、FIR数字滤波器的设计 …………………………………… 13 实验一 熟悉MATLAB环境 一、实验目的 熟悉MATLAB环境的主要操作命令。 学会简单的矩阵输入和数据读写。 掌握简单的绘图命令 用MATLAB编程并学会创建函数。 观察离散系统的频率响应。 二、实验类型 验证型 三、实验仪器 装有MATLAB语言的计算机 四、实验原理 MATLAB的基本命令和函数 五、实验内容和要求 ?数组的加、减、乘、除和乘方运算。输入A=[1 2 3 4]，B=[3 4 5 6]，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.^B，并用stem语句画出A、B、C、D、E、F、G。 用MATLAB实现下列序列： 将（3）中的扩展为以16为周期的函数，绘出四个周期。 将（3）中的扩展为以10为周期的函数，绘出四个周期。 3．，产生并绘出下列序列的样本： 4．绘出下列时间函数的图形，对轴、轴以及图形上方均须加上适当的标注： 5．编写函数stepshift（n0，n1，n2）实现u（n-n0），n1n0n2，绘出该函数的图形，起点为n1，终点为n2。 6．给定一因果系统，求出并绘制的幅频响应和相频响应。 7．计算序列{8 -2 -1 2 3}和序列{2 3 -1 3}的离散卷积，并作图表示卷积结果。 8．求以下差分方程所描述系统的单位脉冲响应 六、思考题 比较实验内容第2题中（4）和（5）两小题的结果，试说明对于周期性信号，应当如何采样，才能保证周期扩展后与原信号保持一致？ 对于有限长序列，如何用MATLAB计算其DTFT？ 对于由两个子系统级联或并联的系统，如何用MATLAB计算它们的幅频响应与相频响应？ 七、实验报告要求 简述实验目的及原理。 按实验步骤附上实验程序。 按实验步骤附上有关离散系统的频率特性曲线。 简要回答思考题。 备注：本实验用到的主要函数 h=freqz(num,den,w) y=filter(num,den,x) c=conv(a,b) 实验(FFT)及其应用 在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉MATLAB中的有关函数。 熟悉应用FFT对典型信号进行频谱分析的方法。 了解应用FFT进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用FFT。 熟悉应用FFT实现两个序列的线性卷积和相关的方法。 二、实验类型 演示型 三、实验仪器 装有MATLAB语言的计算机 四、实验原理 ?在各种信号序列中，有限长序列信号处理占有很重要地位，对有限长序列，我们可以使用离散Fouier变换(DFT)。这一变换不但可以很好的反映序列的频谱特性，而且易于用快速算法在计算机上实现，当序列x(n)的长度为N时，它的DFT定义为： 反变换为： ??? 有限长序列的DFT是其Z变换在单位圆上的等距采样，或者说是序列Fourier变换的等距采样，因此可以用于序列的谱分析。 ??? FFT并不是与DFT不同的另一种变换，而是为了减少DFT运算次数的一种快速算法。它是对变换式进行一次次分解，使其成为若干小点数的组合，从而减少运算量。常用的FFT是以2为基数的，其长度 。它的效率高，程序简单，使用非常方便，当要变换的序列长度不等于2的整数次方时，为了使用以2为基数的