信号实验(二).docVIP

实验二：线性系统的时域分析 一、实验目的： 1、掌握连续时间系统的零输入响应，单位冲激响应，零状态响应和全响应的求解方法。 2、掌握离散时间系统的零输入响应，单位冲激响应，零状态响应和全响应的求解方法。 3、掌握连续时间信号与离散时间信号的卷积运算方式，分析建立信号波形间的联系。 二、课前验证实验内容： 1、练习并验证实验指导书上实验五----八的内容。 2、连续时间系统响应 （1）冲激响应： ·impulse(b,a) ·Impulse(b,a,t) ·Impulse(b,a,t1:p:t2) ·y=impulse(b,a,t1:p:t2) （2）阶跃响应 ·Step(b,a) ·Step(b,a,t) ·Step(b,a,t1:p:t2) （3）LTI连续系统的响应 ·连续系统的零输入响应：initial(sys,[y’(0-)，y(0-)…],t) ·连续系统的零状态响应：lsim(b,a,X,t) lsim(sys,X,t,0) ·连续系统的全响应： lsim(sys,X,t,zi) 例题： 已知微分方程: ，， （1）写出人工求解的零输入响应表达式并画波形；利用零输入响应函数initial(A,B,C,D,y(0-),t) 求解画波形；利用lsim函数在激励为零的情况下求解画波形，比较三种结果是否一致。 （2）写出人工求解的冲激响应表达式并画波形，利用冲激响应函数impulse(sys,t) 求解画波形，比较结果。 （3）若激励信号为，写出人工求解的零状态响应表达式并画波形；利用阶跃响应函数step(sys,t) 求解画波形；利用零状态响应函数lsim求解画波形；比较结果。 （4）利用零输入响应与零状态响应叠加求解全响应并画波形；利用初始条件以及lsim函数求解全响应并画波形；比较结果。 解：由人工计算得：  （1） t=0:0.1:5; yzi=-3*exp(-2*t)+5*exp(-t); subplot(331) plot(t,yzi); axis([0,5,0,3]); grid on a=[1 3 2]; b=[1]; zi=[1 2]; [A,B,C,D]=tf2ss(b,a); sys=ss(A,B,C,D); yzi2=initial(sys,zi,t); subplot(334) initial(sys,zi,t); axis([0,5,0,3]); grid on f=zeros(1,length(t)); yzi3=lsim(sys,f,t,zi); subplot(337) lsim(sys,f,t,zi); axis([0,5,0,3]); grid on （2） h1=-exp(-2*t)+exp(-t) subplot(332) plot(t,h1); axis([0,5,0,0.4]); grid on h2=impulse(sys,t) subplot(335) plot(t,h2) axis([0,5,0,0.4]); grid on (3) yzs=1/2*exp(-2*t)-exp(-t)+1/2*stepfun(t,0); subplot(333) plot(t,yzs); axis([0,5,0,0.5]); grid on yzs2=step(b,a,t); subplot(336); plot(t,yzs2) axis([0,5,0,0.5]); grid on x2=stepfun(t,0); yzs3=lsim(A,B,C,D,x2,t); subplot(339); plot(t,yzs3) axis([0,5,0,0.5]) （）y1=yzi+yzs; subplot(211) plot(t,y); axis([0,5,0,3]); grid on f2=stepfun(t,0); y2=lsim(A,B,C,D,f2,t,zi); subplot(212) plot(t,y2) axis([0,5,0.5,3]); grid on 3、离散时间系统的响应 ·离散差分方程： ·冲激响应：impz(b,a); impz(b,a,n); impz(b,a,n1:n2); dimpulse(A,B,C,D) dimpulse(b,a) ·阶跃响应：stepz(b,a); stepz(b,a,n); stpez(b,a,n1:n2); dstep(A,B,C,D); dstep(b,a);