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

justlikethisis

clc; clear; 7 D3 ]* j# x+ h5 v4 }3 _( N! z2 \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] + f4 X* z1 v. E( n) O7 FC=[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) P1=P.' Kh=P1*Kw*P ' ?. I4 c0 U# C; d: L; BCh=P1*Cw*P S=inv(M.^(1/2))*P + Q- t% x5 O0 G# hS1=S.' 6 ?) T' r( k( Q9 ]) e# IKS=S1*K*S CS=S1*C*S m1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; $ L- Z# I& ~; e2 m. H" u4 dKS=simplify(subs(KS)) 6 x! R' v6 o* x2 ICS=simplify(subs(CS)) 8 c7 K8 k8 t9 {7 {  C/ @; L/ w8 P 8 y& |5 R* f& ~: w8 O/ x" qclc; clear; ! ~2 L7 F8 [* z3 z, ssyms 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 ) N* H8 Z; J- U- m1 G4 v m1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; " X% D5 I4 G& f# W6 B" I0 Q. BM=[m1 0;0 m2] K=[k1 -k1;-k1 k1+k2] C=[c1 -c1;-c1 c1+c2] Kw=inv(M.^(1/2))*K*inv(M.^(1/2)) 0 T# w$ ?& q* O$ d: y+ a3 H- p$ w' M, FCw=inv(M.^(1/2))*C*inv(M.^(1/2)) $ \7 H+ l  U) M2 P5 P) F, q# n[P,Z]=eig(Kw) ' X# v+ X; V! m- MP1=P.' ' w1 h; V% }- A9 M6 n  pKh=P1*Kw*P Ch=P1*Cw*P S=inv(M.^(1/2))*P S1=S.'   M& |. Z3 i3 _) k5 c2 b/ ]KS=S1*K*S CS=S1*C*S , e) x2 u* T4 |, H- g# M 疑问： 前一个程序得到的结果化简后： KS = 34.5311         0          0   12.1726 CS = . [% S8 ~. O; B# a: h     3.4531         0 8 B; C7 F- Z) h1 K         0    1.2173 2 K& }# O; d0 v- T7 S后一个程序得到： KS = - p4 {0 [: J! V( ~  0.2482   -0.0000 3 x- I0 n& u3 z* W   -0.0000   12.0851 8 e( L9 j1 }: P. u  v- ^CS = 5 j! A2 {, S$ {" X. ^$ s) e  t   0.0248   -0.0000    -0.0000    1.2085 ; ~, [  |0 r3 O3 |两个程序大致相同，为什么会出现结果不同的情况

shapeofyou888

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

modengxian111

统一用数值算法！