t=0:0.01:3; sys1=ss(sys); y=[1 1 ];
u=zeros(1,length(t)); rzi=lsim(sys1,u,t,y);%
subplot(311);plot(t,rzi);title('零输入响应yzi(t)');ylabel('rzi(t)') %% 零状态响应
syms s f
f=ilaplace(3/((s+2)*(s+3)));t=0:0.01:3; rzs=3*exp(-2*t)-3*exp(-3*t); subplot(312);plot(t,rzs)
title('零状态响应');ylabel('rzs(t)') %% 全响应=零状态响应+零输入响应
r=rzi+rzs';
subplot(313);plot(t,r);
title('全响应');ylabel('r(t)');xlabel('时间（秒）');ylabel('r(t)')
21