matlab中计算矩阵时直接赋值和算出符号表达式后再代入后结果不同

justlikethisis

clc; clear; syms m1 m2 k1 k2 c1 c2 x0 dx0 K M C Kw Cw Kh Ch P P1 Z S S1 KS CS SF x0 dx0 r0 dr0 M=[m1 0;0 m2] ! d5 s' A& V- \" pK=[k1 -k1;-k1 k1+k2] 2 @* v) K  D* b6 S: ]% A2 cC=[c1 -c1;-c1 c1+c2] 1 L0 O. z3 O1 h% CKw=inv(M.^(1/2))*K*inv(M.^(1/2)) # q2 m5 T9 S2 ]5 |8 T$ P% UCw=inv(M.^(1/2))*C*inv(M.^(1/2)) [P,Z]=eig(Kw) # l3 x. t  C4 L" c1 UP1=P.' ; ]0 y! C& K7 z# h( @' l  HKh=P1*Kw*P Ch=P1*Cw*P ' b; Q* W- T% M% L% _8 `) I' uS=inv(M.^(1/2))*P S1=S.' 7 _1 J( M. w- aKS=S1*K*S CS=S1*C*S 4 m2 u4 L9 }& c! S$ \' K& L4 b m1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; KS=simplify(subs(KS)) " `1 K/ [8 q  t' H5 D2 gCS=simplify(subs(CS)) ( X. l  C" l9 c( U% hclc; clear; syms m1 m2 k1 k2 c1 c2 x0 dx0 K M C Kw Cw Kh Ch P P1 Z S S1 KS CS SF r0 dr0 x0 dx0 3 ~7 z* ]5 t  w m1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; M=[m1 0;0 m2] K=[k1 -k1;-k1 k1+k2] 2 @+ _" S6 {0 O1 k3 MC=[c1 -c1;-c1 c1+c2] Kw=inv(M.^(1/2))*K*inv(M.^(1/2)) Cw=inv(M.^(1/2))*C*inv(M.^(1/2)) [P,Z]=eig(Kw) 4 a6 m& [9 E2 G" U' _/ x! t- e: j  VP1=P.' 7 y2 I7 K$ z) PKh=P1*Kw*P . K9 p5 S1 O: g6 c# Y1 iCh=P1*Cw*P S=inv(M.^(1/2))*P 9 f% l! ]7 ~$ d6 L" {' o5 W1 bS1=S.' / Q8 N0 c, B  v. X5 r6 }+ Q+ t% b0 ZKS=S1*K*S CS=S1*C*S 9 d1 G" \4 g, K( D% ]% D; _ * t, `0 K# g2 W# D, @ 6 m( Y" I2 O' E3 c疑问： 前一个程序得到的结果化简后： KS = 34.5311         0 / Y, P  F& i. B% E9 Y; g! Y7 T         0   12.1726 1 h6 h# c6 i* l- f; |" nCS =     3.4531         0          0    1.2173 9 k: C& H9 Z1 t$ g$ E8 O后一个程序得到： " V6 D! p* P. h9 w$ JKS = % |0 Q5 C2 U' {. p+ l- q- ]  0.2482   -0.0000    -0.0000   12.0851 / G. {2 l; l" b+ }, t4 ]" K+ PCS = 7 `2 ]3 i% ?6 L  A: z6 ^   0.0248   -0.0000 " x) l0 d* J+ [( t9 D& u$ z   -0.0000    1.2085 两个程序大致相同，为什么会出现结果不同的情况

. K; E$ K) V4 }6 f. _1 f0 a

shapeofyou888

一个是符号算法、一个是数值算法。 结果可能不一样。

modengxian111

统一用数值算法！