ch3_4_1.m

来自「这是书籍:Matlab控制系统与应用的源码,」· M 代码 · 共 19 行

A=[0 1 0;0 0 1;-14 -11 -6];
B=[0;0;1];
C=[1 0 0];
D=0;
sys=ss(A,B,C,D);
Q=[200 0 0;0 1 0;0 0 1];
R=1;
N=[0;0;1];
[K,P,E]=lqr(A,B,Q,R,N)
AC=A-B*K;BC=B*K(1);CC=C;DC=0;
sysc=ss(AC,BC,CC,DC);
[y,t,x]=step(sysc,8);
figure(1)
plot(t,x)
grid
figure(2);
plot(t,y)
grid