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

justlikethisis · 发表于 2022-7-19 10:33

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

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 x0 dx0 r0 dr0
M=[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))
Cw=inv(M.^(1/2))*C*inv(M.^(1/2))
[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

m1=9;m2=1;c1=0.3;c2=0.9;k1=3;k2=9;x0=[1;0];dx0=[0;0];

KS=simplify(subs(KS))
CS=simplify(subs(CS))

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

m1=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]
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))
[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

疑问：
前一个程序得到的结果化简后：
KS =
34.5311       0
      0 12.1726
CS =

3.4531       0
      0 1.2173
后一个程序得到：
KS =
  0.2482 -0.0000
-0.0000 12.0851
CS =
0.0248 -0.0000
-0.0000 1.2085
两个程序大致相同，为什么会出现结果不同的情况

shapeofyou888 · 发表于 2022-7-19 13:18

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

modengxian111 · 发表于 2022-7-19 15:22

统一用数值算法！

帐号		自动登录	找回密码
密码			注册

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

EDA365欢迎您登录！

浏览过的版块

推荐内容 /1