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

justlikethisis

clc; ; U- S# G9 C' a6 f" Q3 Y4 Uclear; 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] . E5 T% q8 }5 R, @: CK=[k1 -k1;-k1 k1+k2] & F2 e1 y: Y& BC=[c1 -c1;-c1 c1+c2] ) q3 ~9 a6 N9 K9 k! `+ o! \) LKw=inv(M.^(1/2))*K*inv(M.^(1/2)) ) W4 G9 G, _1 ZCw=inv(M.^(1/2))*C*inv(M.^(1/2)) [P,Z]=eig(Kw) P1=P.'   P3 V) a3 ]( d+ Q8 _/ LKh=P1*Kw*P 9 Y. k, t+ K! KCh=P1*Cw*P ( n$ B7 W, \/ e, I! R! a  X6 RS=inv(M.^(1/2))*P S1=S.' KS=S1*K*S CS=S1*C*S 8 p! c& q# a( Y7 l5 i6 Z% x/ m2 { m1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; KS=simplify(subs(KS)) 5 e3 P  I6 G# tCS=simplify(subs(CS)) % E* J: v  g) K0 y3 o' A' m+ f 2 [  c& ~, ?/ T1 J% Z4 f9 z$ lclc; ) j" Z7 B* L. |6 v4 Y5 v9 {4 `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 1 f6 ~: `- ?* e& p0 E6 N2 Om1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; " B; _; q0 {7 R6 }7 [' a7 kM=[m1 0;0 m2] K=[k1 -k1;-k1 k1+k2] C=[c1 -c1;-c1 c1+c2] : s; ?3 V/ b4 k: _( `Kw=inv(M.^(1/2))*K*inv(M.^(1/2)) Cw=inv(M.^(1/2))*C*inv(M.^(1/2)) [P,Z]=eig(Kw) 8 s! q$ t+ e9 w5 HP1=P.' Kh=P1*Kw*P 9 K' D( s( G3 V8 w% R4 `Ch=P1*Cw*P 2 e, c$ Z' C! C; S! ^3 n+ C/ ZS=inv(M.^(1/2))*P S1=S.' ; J1 a! x; y+ f% K' RKS=S1*K*S CS=S1*C*S 1 @" }+ I& z% G! j& T+ d# B ; y' \: j, y* y7 W) J6 T" Z疑问： 前一个程序得到的结果化简后： KS =   N, A# s! J! i8 [  I2 o3 y$ N9 ]5 Q34.5311         0 5 h% @) ?! Y: y! X2 u% ^         0   12.1726 6 i. v4 m" p1 W" }8 I/ XCS =     3.4531         0 7 B  B4 I# W  Z9 n5 u. G# L8 M9 ^         0    1.2173 后一个程序得到： / K6 p% q, s8 R  V: ]KS = * T% Z- |; ?% E( M% ?/ y9 i  0.2482   -0.0000 3 b1 _! m5 t- h9 |) ]   -0.0000   12.0851 CS = ( Y  o. B/ N9 |" ^) H' y) J  m   0.0248   -0.0000 9 ]; i: X- k; x2 f' B6 @   -0.0000    1.2085 5 ~. K- }3 I/ D7 M; p% A, F两个程序大致相同，为什么会出现结果不同的情况

shapeofyou888

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

modengxian111

统一用数值算法！