ex4_19.m

反馈控制系统设计与分析（matlab程序）

G1=tf([5,100],[1,8,32,80,100]); G=ss(G1);
u0=[2,-1,-24]; Y0=[0.3; 0.5; 0; 0];
hatY0=Y0-[0; G.c*G.b*u0(1); 
     G.c*G.b*u0(2)+G.c*G.a*G.b*u0(1); 
     G.c*G.b*u0(3)+G.c*G.a*G.b*u0(2)+G.c*G.a^2*G.b*u0(1)];
Wo=obsv(G); X0=inv(Wo)*hatY0;  
cc=1; dd=[-1,1,0,5];
[y,t,x]=ss_simu(G,X0,cc,dd,0.05,100);
plot(t,y); axis([0,5,0,1.2])