ex3_8.m

来自「反馈控制系统分析与设计,虽然有点抽象」· M 代码 · 共 7 行

A=[-4.1053,6.6842,-2.3158,20; -0.5789,-0.7368,0.4632,4.2;
   -0.0526,-1.1579,0.6421,1.8;-1.2632,-1.1228,1.0772,7.2];
B=[1;0.2;0.1;0.4]; C=[0,0,8,-2]; D=0; G=ss(A,B,C,D);
W=diag([1,2,3,4]); V=lyap(G.a,-W); 
[key,sdet]=posdef(V)
eig(G.a)