D,1)>>num=020.0000-29.0000-13.0000den=1-58-4②P361.5(3)已知系统传递函数阵:G(s)=求系统的状态空间模型num=[1422;3110];den=[1232];[A,B,C,D]=tf2ss(num,den)>>A=-2-3-2B=1C=2-10D=11000-5-8-630100③P562.3(3)已知系统方成如下:=x+uy=x求输入和初值为以下值时的状态响应和输出响应:u(t)=1(t),x(0)=a=[01;-6-5];b=[1;0];c=[1.50.5];d=0;G=ss(a,b,c,d);x0=[1;1];symsst;G0=inv(s*eye(size(a))-a);x1=ilaplace(G0)*x0G1=inv(s*eye(size(a))-a)*b;x2=ilaplace(G1/s)x=x1+x2y=c*xforI=1:61;tt=0.1*(I-1);xt(:,I)=subs(x(:),'t',tt);yt(I)=subs(y,'t',tt);end;plot(0:60,[xt;yt]);xlabel('t(s)');ylabel('ampitude');title('系统的状态响应和输出响应');运行结果如下:>>x1=[-3*exp(-3*t)+4*exp(-2*t)][-8*exp(-2*t)+9*exp(-3*t)]x2=[2/3*exp(-3*t)-3/2*exp(-2*t)+5/6][-1+3*exp(-2*t)-2*exp(-3*t)]x=[-7/3*exp(-3*t)+5/2*exp(-2*t)+5/6][-5*exp(-2*t)+7*exp(-3*t)-1]y=5/4*exp(-2*t)+3/4实验二系统的能空性、能观测性、稳定性分析及实现实验日期2013年5月21日成绩:______
全文预览
