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

justlikethisis

clc; - M6 U2 u  B8 n  d( j$ yclear; 3 s  d5 Z& ]' ~, tsyms 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 1 r0 C" d/ ^- A5 c. q0 S9 fM=[m1 0;0 m2] 7 e3 A; a  s2 v& n6 D/ _K=[k1 -k1;-k1 k1+k2] 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)) 8 a& Y. e5 y6 o" p+ u1 o0 y9 r[P,Z]=eig(Kw) P1=P.' Kh=P1*Kw*P 3 a0 I; O* W+ [, d! kCh=P1*Cw*P ! w" ^# J2 K: j* q7 K* ES=inv(M.^(1/2))*P ; h+ G: N( L" K- `9 h( t3 RS1=S.' KS=S1*K*S CS=S1*C*S ) N. U0 Y4 T( p! t m1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; ! ]. H: R  P* e0 OKS=simplify(subs(KS)) CS=simplify(subs(CS)) * F+ |% E& D" U8 h* l7 e 9 F, @% ?. h& R( W1 ^6 S) }) B2 R, |clc; ) F- m+ l0 P( u7 y6 @0 aclear; 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 ) c/ I! ]6 ]& a$ } . H- z+ R* N! g; S4 z9 b9 Im1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; M=[m1 0;0 m2] K=[k1 -k1;-k1 k1+k2] % q& ]5 \- `& p1 \. q* K; vC=[c1 -c1;-c1 c1+c2] ) G5 z" q" {0 w# a3 ?) J! ^  |) `Kw=inv(M.^(1/2))*K*inv(M.^(1/2)) Cw=inv(M.^(1/2))*C*inv(M.^(1/2)) ' T% q5 P5 q4 Y" z- r/ L/ D[P,Z]=eig(Kw) P1=P.' Kh=P1*Kw*P / k9 N! L+ Q* e' |6 iCh=P1*Cw*P S=inv(M.^(1/2))*P S1=S.' - d- K2 F" Q' T1 o" F) t( Z! i8 o% EKS=S1*K*S % t0 V6 U# p, z$ NCS=S1*C*S + k/ \' o, m/ f9 s3 E+ j8 X7 K 疑问： 前一个程序得到的结果化简后： KS =   e, f+ K  V$ E34.5311         0 7 j, k. M. K% ?* j         0   12.1726 4 W- V- E5 g" X7 c* g/ E. j: O: `& fCS = 4 x5 l8 d' R0 a7 t' c  I % n* v$ \. x0 d7 z    3.4531         0          0    1.2173 ! Q9 g) p% {0 M后一个程序得到： KS = & F9 h+ \, E' ~1 b4 |0 M& N  0.2482   -0.0000 1 L( ]* |; A4 S9 P. F   -0.0000   12.0851 CS = . u5 y* R1 W% c0 _3 K   0.0248   -0.0000 2 x6 D& V4 g1 D; U: Y! ?7 h! i! Y   -0.0000    1.2085 两个程序大致相同，为什么会出现结果不同的情况

shapeofyou888

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

modengxian111

统一用数值算法！