ex6_5.m

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

A=[-0.3,0.1,-0.05; 1,0.1,0; -1.5,-8.9,-0.05]; 
B=[2; 0; 4]; x0=zeros(3,1);
Q=eye(3,3); R=1; C=[1,2,3]; D=0; Kc=lqr(A,B,Q,R)

[y,t,x]=step(ss(A-B*Kc,B,C,D)); 
plot(t,x,':'), figure;  plot(t,y)