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

justlikethisis

clc; , r0 T5 U( Q# B; Kclear; 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] K=[k1 -k1;-k1 k1+k2] 6 d3 R- x1 x) M5 w& v+ }: C* ^5 ^C=[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)) - X( B5 S$ S" d$ i( V4 c: `[P,Z]=eig(Kw) $ A% [8 i6 O3 d# f& N- M$ Y+ a5 f+ FP1=P.' Kh=P1*Kw*P ' A* j1 t" [4 U7 H2 F* t6 @9 F2 ECh=P1*Cw*P S=inv(M.^(1/2))*P S1=S.' KS=S1*K*S CS=S1*C*S ! C; P# U. F, j7 \ 5 l+ c# E( r9 `& |' A4 [m1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; * d0 i- Y& S6 z6 G9 W KS=simplify(subs(KS)) 9 M  ]+ T7 J$ P5 B$ S. ^/ }CS=simplify(subs(CS))   c% q7 }. J! ~6 |0 |clc; ! G( V6 P5 j4 U" I/ Wclear; - j7 T) v. u8 Fsyms 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 0 N# D$ e) C+ s( P5 G& Y 9 Z5 s, R* N  N" Q0 fm1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; M=[m1 0;0 m2] , N& V3 b/ }4 u! E+ `K=[k1 -k1;-k1 k1+k2] 1 i# g4 k4 K. y2 Q$ o7 t3 K$ s. nC=[c1 -c1;-c1 c1+c2] / D- r, N4 _# S1 Z# _Kw=inv(M.^(1/2))*K*inv(M.^(1/2)) Cw=inv(M.^(1/2))*C*inv(M.^(1/2)) [P,Z]=eig(Kw) ; p6 f9 ~; Q( sP1=P.' Kh=P1*Kw*P Ch=P1*Cw*P 0 q3 q* q2 F0 V  NS=inv(M.^(1/2))*P / w) L# _. {" kS1=S.' KS=S1*K*S CS=S1*C*S # r' k/ c+ K9 l( A" [! T+ o疑问： 前一个程序得到的结果化简后： KS = " ?9 O7 ~+ o+ C2 I& e+ [1 ]1 R34.5311         0 6 g& h* E# w# v7 i3 G/ W" T5 R+ r         0   12.1726 CS =     3.4531         0 3 T1 ~5 [3 U/ Y' A6 w         0    1.2173 后一个程序得到： KS = 4 I' b- F2 K0 j0 K% L( p5 z2 w6 z5 L  0.2482   -0.0000 6 p4 D& L" ]! r0 _+ F- A   -0.0000   12.0851 1 f% F  Z. P1 j6 s8 LCS =    0.0248   -0.0000 9 [; ]6 M% w7 c2 \9 Y' {: H   -0.0000    1.2085 2 g/ `9 b0 \" V3 x' p: L7 |两个程序大致相同，为什么会出现结果不同的情况

shapeofyou888

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

modengxian111

统一用数值算法！