贵州大学实验报告 学院:计信 专业: 计科 班级:课程名称:数字信号处理 姓名 学号 实验序号 2 实验时间 指导教师 成绩 实验项目名称 离散系统的差分方程、冲激响应和卷积分析 实验目的 加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。 实验内容 编制程序求解下列两个系统的单位冲激响应和阶跃响应,并绘出其图形。 1]}4[]3[]2[]1[{25.0][nxnxnxnxny 2]1[][]2[125.0]1[75.0][nxnxnynyny 实验环境 MATLAB 7.4.0(R2007a) 实验 原理 离散系统 ][nx][nyDiscrete-timesystme 其输入、输出关系可用以下差分方程描述: MkkNkkknxpknyd00][][ 输入信号分解为冲激信号,mmnmxnx][][][。记系统单位冲激响应 ][][nhn ,则系统响应为如下的卷积计算式: mmnhmxnhnxny][][][][][ 当Nkdk,...2,1,0时,h[n]是有限长度的(n:[0,M]),称系统为FIR 系统;反之,称系统为IIR 系统。 在 MATLAB 中,可以用函数 y=Filter(p,d,x) 求解差分方程,也可以用函数 y=Conv(x,h)计算卷积。 实验编程及运行结果 1.输入为单位冲击函数)(t引起的零状态响应称为单位冲击响应。 2.输入为单位阶跃函数)(t引起的零状态响应称为单位阶跃响应。 3.在MATLAB中,可用dirac函数表示单位冲击函数,用heaviside函数表示阶跃函数,用函数y=Filter(b,a,x) 求解差分方程。其中b表示该差分方程的右边系数,a表示该差分方程的左边系数, x表示输入的序列函数。 4.对于系统1:y[n]=0.25{x[n-1]+x[n-2]+x[n-3]+x[n-4]} 可知在输入为单位冲击函数)(t时,有 b=[0 0.25 0.25 0.25 0.25 ]; a=[1]; x=dirac(n); 在输入为单位阶跃函数)(t时,有 b=[0 0.25 0.25 0.25 0.25 ]; a=[1]; x=heaviside (n); 所以系统1 的单位冲激响应和阶跃响应的编写代码及输出图形如下图所示: %y[n]=0.25{x[n-1]+x[n-2]+x[n-3]+x[n-4]}单位冲击响应 n=1:30; b=[0 0.25 0.25 0.25 0.25 ]; a=[1]; x=dirac(n); y=filter(b,a,x); stem(n,y) title('y[n]=0.25{x[n-1]+x[n-2]+x[n-3]+x[n-4]}单位冲击响应') xlabel('n') 图1-1 %y[n]=0.25{x[n-1]+x[n-2]+x[n-3]+x[n-4]}单位阶跃响应 n=1:30; b=[0.25 0.25 0.25 0.25 ]; a=[1]; x=heaviside(n); y=filter(b,a,x); stem(n,y) title('y[n]=0.25{x[n-1]+x[n-2]+x[n-3...