请教Riccati方程的解为何总是出错

tiange

syms x12 x11 x22 x21 k11 k12 k13 k14 k21 k22 k23 k24;
L=1.8;
a1=1.875/1.8;
a2=4.964/1.875;
Ka=0.0002646;
A=[0 0 1 0;0 0 0 1;-12.7449 0 -0.07 0;0 -495.50 0 -0.445];
b112=sin(a1*x12)-sinh(a1*x12)+(sin(a1*L)+sinh(a1*L))*(cosh(a1*x12)-cos(a1*x12))/(cos(a1*L)+cosh(a1*L));
b111=sin(a1*x11)-sinh(a1*x11)+(sin(a1*L)+sinh(a1*L))*(cosh(a1*x11)-cos(a1*x11))/(cos(a1*L)+cosh(a1*L));
B11=b112-b111;
b122=sin(a1*x22)-sinh(a1*x22)+(sin(a1*L)+sinh(a1*L))*(cosh(a1*x22)-cos(a1*x22))/(cos(a1*L)+cosh(a1*L));
b121=sin(a1*x21)-sinh(a1*x21)+(sin(a1*L)+sinh(a1*L))*(cosh(a1*x21)-cos(a1*x21))/(cos(a1*L)+cosh(a1*L));
B12=b122-b121;
b212=sin(a2*x12)-sinh(a2*x12)+(sin(a2*L)+sinh(a2*L))*(cosh(a2*x12)-cos(a2*x12))/(cos(a2*L)+cosh(a2*L));
b211=sin(a2*x11)-sinh(a2*x11)+(sin(a2*L)+sinh(a2*L))*(cosh(a2*x11)-cos(a2*x11))/(cos(a2*L)+cosh(a2*L));
B21=b212-b211;
b222=sin(a2*x22)-sinh(a2*x22)+(sin(a2*L)+sinh(a2*L))*(cosh(a2*x22)-cos(a2*x22))/(cos(a2*L)+cosh(a2*L));
b221=sin(a2*x21)-sinh(a2*x21)+(sin(a2*L)+sinh(a2*L))*(cosh(a2*x21)-cos(a2*x21))/(cos(a2*L)+cosh(a2*L));
B22=b222-b221;
B=Ka*[0 0; 0 0;B11 B12;B21 B22];
K=[k11 k12 k13 k14;k21 k22 k23 k24];
AA=A-B*K;
BB=[0 0;0 0;0 0;0 0];
R=[1 0;0 1];
Q=[12.78 0 0 0;0 495.89 0 0;0 0 1 0;0 0 0 1];
[X]=care(AA,BB,Q,R)

[ 本帖最后由 eight 于 2007-11-22 16:31 编辑 ]

eight

请先耐心看看置顶帖《聚宝盆》，然后重新整理你的问题

[ 本帖最后由 ChaChing 于 2010-6-14 00:47 编辑 ]