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] & H+ U: g0 t- m* L  N+ NK=[k1 -k1;-k1 k1+k2] ' W' v# @4 v6 ^- BC=[c1 -c1;-c1 c1+c2] 7 s  D# B+ v  x$ b2 t. oKw=inv(M.^(1/2))*K*inv(M.^(1/2)) % N/ C# _) y& f' i# K# jCw=inv(M.^(1/2))*C*inv(M.^(1/2)) [P,Z]=eig(Kw) ' ]( M$ n- R2 u9 sP1=P.' ( Q1 b7 {! \1 e+ K/ FKh=P1*Kw*P Ch=P1*Cw*P S=inv(M.^(1/2))*P S1=S.' KS=S1*K*S # q6 P2 u2 d8 D  u# K% ICS=S1*C*S " Y2 k' ]6 L7 }7 {m1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; 0 `7 a0 T8 s$ J& N1 T- N9 M: EKS=simplify(subs(KS)) CS=simplify(subs(CS)) 2 E( {8 k6 z+ ^. P$ J % ?' \; D+ i5 X+ x) ] , [/ l- U  }6 L8 Wclc; ' o- ?$ j1 T- K6 d# D) T# pclear; 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 / r% ]( N3 l$ P1 n7 H1 P m1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; 9 p2 [2 e, j) GM=[m1 0;0 m2] K=[k1 -k1;-k1 k1+k2] : }/ \$ p0 i/ BC=[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)) ' T7 _; ]* O; o: |[P,Z]=eig(Kw) / k  v+ w8 ~& ~- e7 W$ |! HP1=P.' : Q# g$ e: s, N( M$ h* `Kh=P1*Kw*P ) h$ s' H7 F9 d0 g' `Ch=P1*Cw*P $ b; \/ O: T4 y4 S9 p7 SS=inv(M.^(1/2))*P S1=S.' KS=S1*K*S $ y4 `' R( v$ S! iCS=S1*C*S 疑问： 前一个程序得到的结果化简后： " v3 d$ s+ q! k! f4 i' XKS = 8 v3 _' N, o* o; j34.5311         0          0   12.1726 CS = * r. B2 I4 u3 k& i: U8 B; w% u   v2 u, B7 d# }* l$ V0 H, T    3.4531         0 : C3 s- K  |. U) G4 {+ Z         0    1.2173 后一个程序得到： * F* y4 z# h: t! j& \KS =   0.2482   -0.0000    -0.0000   12.0851 CS =    0.0248   -0.0000 0 m$ H. u! S6 K0 e   -0.0000    1.2085 两个程序大致相同，为什么会出现结果不同的情况

shapeofyou888

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

modengxian111

统一用数值算法！