ex6_7.m

来自「反馈控制系统设计与分析（matlab程序）」· M 代码 · 共 9 行

A=[-0.2,0.5,0,0,0;0,-0.5,1.6,0,0;0,0,-14.3,85.8,0;
    0,0,0,-33.3,100;0,0,0,0,-10];
B=[0; 0; 0; 0; 30]; Q=diag([10,20,6,2,5]); R=1;
C=[1,0,0,0,0]; D=0; rho=[1,10,100,100,10000];
for i=1:length(rho)
   [K,P]=lqr(A,B,rho(i)*Q,R); Ac=A-B*K;
   step(ss(Ac,B,C,D)); hold on
end