ex651.m

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

A=[-0.4 0 -0.01;1 0 0;-1.4 9 -0.02];B=[6; 0;8];G=[1;0;0];
C=[0 0 1];H=0;D=0;R=1;
Q=eye(3,3);Q0=1;R0=1;Xi=7e-4;theta=1e-8;
W=[Q,zeros(3,1);zeros(1,3),R];
V=[Xi*G*G',zeros(3,1);zeros(1,3),theta];
[Af,Bf,Cf,Df]=lqg(A,B,C,D,W,V);
disp('LQG控制器的状态方程为')
Gc=ss(Af,Bf,Cf,Df)