ex641.m

来自「详细介绍了Mathlab的使用方法与技巧」· M 代码 · 共 7 行

A=[0 1 0;0 0 1;-5 -10 -18];B=[0;0;1];
C=[1 0 0];D=0;
Q=eye(3,3);R=1;N=[0;0;1];
[K,S,e]=lqr(A,B,Q,R,N);
disp('要求的性能指标求出的LQ控制器的状态反馈矩阵K');K
disp('Riccati的正定解为');S