信号与系统实验二分析.doc

南京工程学院 通信工程学院 实 验 报 告 课 程 名 称 信号与线性系统 实验项目名称 线性系统的时域分析 实验学生姓名 实验学生班级 实验学生学号 实 验 时 间 实 验 地 点 信号与系统实验室 实验成绩评定 指导教师签字 年 月 日 一、实验目的： 1、掌握连续时间系统的零输入响应，单位冲激响应，零状态响应和全响应的求解方法。 2、掌握离散时间系统的零输入响应，单位冲激响应，零状态响应和全响应的求解方法。 3、掌握连续时间信号与离散时间信号的卷积运算方式，分析建立信号波形间的联系。 二、实验原理： （1）连续时间系统调用函数 连续系统的零输入响应：initial(A,B,C,D,zi,t) 冲激响应： impulse(b,a) impulse(b,a,t) impulse(b,a,t1:p:t2) y=impulse(b,a,t1:p:t2)(数值解) 阶跃响应 step(b,a) step(b,a,t) step(b,a,t1:p:t2) y=step(b,a,t1:p:t2) 连续系统的零状态响应： lsim(A,B,C,D,x,t,0) lsim(b,a,x,t) 连续系统的全响应： lsim(b,a,x,t,zi) （2）离散时间系统调用函数 离散系统的零输入响应：dinitial(A,B,C,D,zi) 离散系统的冲激响应： impz(b,a); impz(b,a,n); impz(b,a,n1:n2); y=impz(b,a,n1:n2); dimpulse(A,B,C,D) dimpulse(b,a) y=filter(b,a,Impseq) 离散系统的阶跃响应： stepz(b,a); stepz(b,a,n); stpez(b,a,n1:n2); y=stepz(b,a,n1:n2); dstep(A,B,C,D); dstep(b,a); filter(b,a,stepfun); 离散系统的零状态响应： dlsim(A,B,C,D,x) y=filter(b,a,X) 离散系统的全响应： dlsim(A,B,C,D,x,zi) filter(b,a,x,zi) 三、实验内容： 1、已知微分方程: ，，， 激励信号为。 人工解出系统的零输入响应、单位冲激响应、单位阶跃响应、零状态响应以及全响应，利用Matlab软件画出所有响应的波形图。 t=0:0.01:10; yzi=3*exp(-1*t)+(-1)*exp(-2*t); subplot(2,3,1); plot(t,yzi); xlabel(时间(t));ylabel(yzi(t));title(零输入响应); t=0:0.01:10; yzs=(-2)*exp(-1*t)+exp(-2*t)+1; subplot(2,3,2); plot(t,yzs); xlabel(时间(t));ylabel(yzs(t));title(零状态响应); y=yzi+yzs; subplot(2,3,3) plot(t,y); xlabel(时间(t));ylabel(y(t));title(全响应); t=0:0.01:10; h=2*exp(-1*t)+(-2)*exp(-2*t); subplot(2,3,4); plot(t,h); xlabel(时间(t));ylabel(h(t));title(单位冲激响应); t=0:0.01:10; s=(-2)*exp(-1*t)+exp(-2*t)+1; subplot(2,3,5); plot(t,s); xlabel(时间(t));ylabel(s(t));title(单位阶跃响应); （2）利用函数initial(A，B，C，D，zi，t) 或者lsim(sys，x，t，zi)（当x=0时） 求解零输入响应画波形，与（1）的结果比较，看是否相同。 t=0:0.01:10; a=[1 3 2]; b=[2]; zi=[1 2]; [A,B,C,D]=tf2ss(b,a);