（解ＲＩＣＡＴＴＩ方程）D=0.1;H=3.2;Ch=0.3;Kw=1;Kp=0.05;Tt=0.3;W1=1;W2=1;g=1/0.2;A=[-D/H,1/H;-Ch*Kw/Tt,-1/Tt];B1=[-1/H;-Ch*Kp/Tt];B2=[0;-Ch/Tt];B=[B1,B2];C=[W1,0;0,W2];m1=size(B1,2);m2=size(B2,2);R=[-g^2*eye(m1)zeros(m1,m2);zeros(m2,m1

Question

（解ＲＩＣＡＴＴＩ方程）
D=0.1;H=3.2;Ch=0.3;Kw=1;Kp=0.05;Tt=0.3;W1=1;W2=1;g=1/0.2;
A=[-D/H,1/H;-Ch*Kw/Tt,-1/Tt];
B1=[-1/H;-Ch*Kp/Tt];
B2=[0;-Ch/Tt];
B =[B1,B2];
C=[W1,0;0,W2];
m1=size(B1,2);
m2=size(B2,2);
R =[-g^2*eye(m1) zeros(m1,m2);zeros(m2,m1) eye(m2)];
X=care(A,B,C'*C,R)
运行错误提示：
Attempt to execute SCRIPT care as a function.
Error in ==> care at 16
X=care(A,B,C'*C,R)

Answer 1

用我的代码试试,应该是好的：
clear all;clc;
D=0.1;
H=3.2;
Ch=0.3;
Kw=1;
Kp=0.05;
Tt=0.3;
W1=1;
W2=1;
g=1/0.2;
A=[-D/H,1/H;-Ch*Kw/Tt,-1/Tt];
B1=[-1/H;-Ch*Kp/Tt];
B2=[0;-Ch/Tt];
B =[B1,B2];
C=[W1,0;0,W2];
m1=size(B1,2);
m2=size(B2,2);
R=[-g^2*eye(m1),zeros(m1,m2);zeros(m2,m1),eye(m2)];
X=care(A,B,C'*C,R)
X =
4.4269 0.3423
0.3423 0.1774