数字信号处理实验教案探究.doc

数字信号处理实验教案 （1112341用 目 录 数字信号处理实验要求 1 预备实验 2 实验一 系统响应及系统稳定性 4 实验二 时域采样与频域采样 9 实验三 用FFT对信号作频谱分析 13 实验四 IIR数字滤波器设计及软件实现 18 实验五 FIR数字滤波器设计与软件实现 23 附 录 31 附录1 ：MATLAB编程及绘图方法 31 附录2 ：MATLAB矩阵及矩阵的运算 38 信号处理实验要求 一、实验课程任务与要求 实验目的 信号处理实验教学是为了将学生的计算机操作能力、分析能力、软件设计能力与应用实践结合起来，引导学生由浅入深地掌握信号处理理论与开发工具，具备实际应用的信号处理软件开发2. 实验基本要求（以软件设计为主要表现形式） （1）上机前应准备好实验的程序设计算法描述与关键分析内容； （2）准备好程序测试数据和设备操作步骤，上机调试、运行； （3）完成每个实验后进行数据与程序对比分析； （4）递交实验结果的可执行程序、源程序并演示实验结果预备实验（3学时） （3学时） （3学时） （3学时） （3学时） （3学时） （1）实验预习报告； （2）实验签到； （3）上机实际操作； （4）实验设计报告； 五、教材及参考书 1．教材： [1]高西全，丁玉美.数字信号处理（第三版）西电出版社，2008.3 2. 参考书： [1]张志勇等. MATLAB教程R2010a. 北京航空航天大学出版社. 2010.8 [2]陈怀琛等. MATLAB及在电子信息课程中的应用.电子工业出版社. 2002 [3]程佩清.数字处理出版社2003 预备实验 用于信号处理的MATLAB基本操作 1．实验目的： 学习和掌握MATLAB最基本的矩阵运算与绘图工具。 学习用MATLAB实现信号处理的基本操作。 2．实验步骤： 熟悉MATLAB的使用环境和方法。 练习使用基本的向量生成、矩阵运算、绘图等语句。 利用冒号（：）生成向量 X1=[1 2 3 4 5] X2=[1.000 1.500 2.000 2.500] X3=[5 4 3 2 1] 分别生成3*3，3*4的全0矩阵，全1矩阵和随机矩阵 分别输入矩阵 1 2 3 1.0 1.1 1.2 A= 4 5 6 B= 2.0 2.1 2.2 7 8 9 3.0 3.1 3.2 分别计算A+B，A-B，A+3，A-4，A*B，A.*B，C=inv（A）,A/B,A./B 分别计算sin(x1),cos(x1),exp(x1),lg(x2),sqrt(x2) 3.实验说明 向量的生成 利用冒号（：）生成向量 X=j:k (jk) 生成X=[j,j+1,j+2,…,k-1,k] X=j:I:k 如果I0 且jk或I0且jk,则生成向量 X=[j,j+I,j+2I,…,k-I,k] 如果I0 且jk,则生成向量 X=[j,j+I,j+2I,…,k] 例 X1=1：5 X1=[1 2 3 4 5] X2=1：0.5：3 X2=[1.000 1.500 2.000 2.500 3.000] X3=5：-1：1 X3=[5 4 3 2 1] 矩阵的生成 Zeros 生成全0阵 B= Zeros(n) 生成n*n的全0矩阵 B= Zeros(m,n) 生成m*n的全0矩阵 B= Zeros(Size(A)) 生成与与矩阵A大小相同的全0矩阵 注：m,n必须是一个非负数 ones生成全1阵 rand生成均匀分布的随机阵 二 矩阵的算术运算 加法和减法 对于同维矩阵指令的 A+B B 对于矩阵和标量（一个数）的加减运算，指令为： A+3 9 乘法和除法运算 A*B 是数学中的矩阵乘法，遵循矩阵乘法规则 A．*B 是同维矩阵对应位置元素做乘法 B=inv（A）是求矩阵的逆 A/B 是数学中的矩阵除法，遵循矩阵除法规则 A．/B 是同维矩阵对应位置元素相除 另 表示矩阵的转置运算 三．数组函数 下面列举一些基本函数，他们的用法和格式都相同。 Sin(A),cos(A),exp(A),log(A)(相当于ln) Sgrt(A)开平方 abs(A)求模 real(A)求实部 imag(A)求虚部 式中A可以是标量也可以是矩阵 例： 利用等差向量产生一个正弦值向量 t=0:0.1:10 A=sin(t) PLot(A) 这时候即可看到一个绘有正弦曲线的窗口弹