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 ' C* o4 [% Y" {8 S5 KM=[m1 0;0 m2] 3 R) `, t+ C* h5 sK=[k1 -k1;-k1 k1+k2] C=[c1 -c1;-c1 c1+c2] Kw=inv(M.^(1/2))*K*inv(M.^(1/2)) 8 ?3 T1 A! M" MCw=inv(M.^(1/2))*C*inv(M.^(1/2)) [P,Z]=eig(Kw) : w; q8 A' ?$ u8 B8 X- yP1=P.' Kh=P1*Kw*P 6 A: `7 \% j8 J2 @. c; w! F7 CCh=P1*Cw*P 0 Z/ q* a9 `: LS=inv(M.^(1/2))*P 4 @; N6 B8 p# ES1=S.' KS=S1*K*S 8 W; E5 n% P3 A& yCS=S1*C*S 1 y5 S2 d& T9 Z6 g/ n9 Fm1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; KS=simplify(subs(KS)) ! q* i$ D: j) ?5 Z6 JCS=simplify(subs(CS)) & c8 Z& J6 u* X- `& R0 M * i& B5 b% B+ p  T9 r+ q* s" j # c& G( z  r; |7 Dclc; 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 2 _8 a* Z. }3 r6 M2 s3 R 2 F2 }4 ~% N% k! T6 c) Q- tm1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0]; ! X2 R9 \  a+ t' n! n M=[m1 0;0 m2] # c+ E, T, ^1 A' W* _K=[k1 -k1;-k1 k1+k2] ( r+ K2 `8 O9 t8 t5 d. UC=[c1 -c1;-c1 c1+c2] Kw=inv(M.^(1/2))*K*inv(M.^(1/2)) ' h, ~6 G  D& h$ }1 KCw=inv(M.^(1/2))*C*inv(M.^(1/2)) ) l- N" x5 K9 u# [5 H[P,Z]=eig(Kw) & Z; G- F2 C6 z* Z0 N) Q2 C" JP1=P.' 0 P, Q3 ^9 {/ DKh=P1*Kw*P ( g$ j5 K# T: ~$ q7 bCh=P1*Cw*P S=inv(M.^(1/2))*P + K9 K. ]" C1 x: n8 o( QS1=S.' + @( W* ?  [$ G& lKS=S1*K*S 4 \, E3 ~. v  t1 _- L8 Y, Y/ QCS=S1*C*S 8 P' [( r4 S9 _疑问： 7 Y3 C3 T9 ~. H8 r6 H" W6 f0 W前一个程序得到的结果化简后： KS = 34.5311         0 3 A) a3 ?/ O& B- n         0   12.1726 , m; a! z8 g! O6 e0 P2 s8 T% z1 k* xCS =     3.4531         0 $ Y  \& c. q; q# \* l! Z1 x         0    1.2173 后一个程序得到： KS =   0.2482   -0.0000    -0.0000   12.0851 " T$ s3 |4 f+ @CS = ( M7 }" t# `6 Q$ w   0.0248   -0.0000    -0.0000    1.2085 5 E( B* s2 d: {; n$ j' T. T  c两个程序大致相同，为什么会出现结果不同的情况

( G& u- s+ r1 c. h

shapeofyou888

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

modengxian111

统一用数值算法！