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

justlikethisis

clc; clear; 5 i! s5 G7 l6 @3 ]" S, Dsyms 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] C=[c1 -c1;-c1 c1+c2] , r! H/ A4 z# N# a. }  YKw=inv(M.^(1/2))*K*inv(M.^(1/2)) # V3 `' S. S: P5 ]Cw=inv(M.^(1/2))*C*inv(M.^(1/2)) 9 S- D  T0 L' e7 N[P,Z]=eig(Kw) 4 Z! [( b- a+ ?7 f% P3 r  G- [P1=P.' Kh=P1*Kw*P Ch=P1*Cw*P ! N. _3 ]7 v2 @S=inv(M.^(1/2))*P . C. n. Y# a% W' E7 CS1=S.' KS=S1*K*S ! x% \1 z# H' _9 h3 |* _5 sCS=S1*C*S m1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; 3 n, p0 S; x- d- \2 p % k# z- |  h4 l& XKS=simplify(subs(KS)) CS=simplify(subs(CS)) / i( \9 ?- t, a$ _& X 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 r0 dr0 x0 dx0 ! l1 _4 b- w& q& C m1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; / ]1 a  f. `) Q. s  c) n6 d6 e . a: T3 a9 ^0 D4 pM=[m1 0;0 m2] & f7 R+ P! ^/ t5 N2 sK=[k1 -k1;-k1 k1+k2] 6 G3 J5 h/ s7 b' P9 vC=[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)) / j3 s5 g/ N5 L! j7 E[P,Z]=eig(Kw) P1=P.' Kh=P1*Kw*P Ch=P1*Cw*P S=inv(M.^(1/2))*P S1=S.' KS=S1*K*S CS=S1*C*S 2 r  B( _* N* H  r* U" w' [. @ & q" d/ k! P8 B: U; J / N9 k' ?) l# V; @) S6 U; k疑问： 前一个程序得到的结果化简后： 1 J% ]4 S% {  `/ o( DKS = 34.5311         0          0   12.1726 CS = " D, B; q5 M" R8 ^ % T: k! K) P" w+ y; T2 M$ L$ f) {    3.4531         0 8 ~* A, @: g# j         0    1.2173 后一个程序得到： - V: O/ A! ?. s) D: y* a4 [: CKS = ; `4 e, O# z. y4 G& Z9 I! N  0.2482   -0.0000    -0.0000   12.0851 CS =    0.0248   -0.0000 : K- z& t/ x" l   -0.0000    1.2085 4 V4 A$ ]( S! y% ~2 W4 }6 t两个程序大致相同，为什么会出现结果不同的情况

# K+ T7 F0 a* c1 a7 [( {' U+ v

shapeofyou888

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

modengxian111

统一用数值算法！