数学建模社区-数学中国

标题: 请教一个解方程组的问题 [打印本页]

作者: prince12th    时间: 2012-8-27 17:00
标题: 请教一个解方程组的问题
用solve语句解这样一个方程组, O# h7 {0 R  ~
(x1 - a1)^2 == b^2*(3*x1^6 + 6*x1^2*x2^4)% z  f' N+ r/ C6 R$ p2 E# J$ p/ \
(x2 - a2)^2 == b^2*(3*x2^6 + 6*x2^2*x1^4)
; R. t$ v% F* x  @a1,a2已知,b是0.001710531 S8 e) a4 V+ z/ S
7 j; I4 b, q$ E/ N
当设a1,a2都为1时,结果出来的是. W$ G0 F+ |. a4 E1 g
{{{x1 -> 1373.46, x2 -> 13.4299},
  b  j5 r7 \0 B  H( m( ^% O% E{x1 -> 10987.7,    x2 -> 15.3663}, # J" t$ C% U" X, K4 X- Z1 [/ _* B
{x1 -> 1.05482*10^6,   x2 -> 17.8499}, : R, _4 _8 s# X' A
{x1 -> 2.81284*10^6, x2 -> 17.8426}}}
$ d2 c* |3 O" v9 K$ O, Q
' S( h- R- j- R: f+ O: Zx1,x2不是应该对称的么,而且把结果代回原来的式子里也不正确,为什么呢?
作者: sycxyc    时间: 2012-8-27 18:39
有小数点的计算精度是有限的,并且高次方程没有精确解,如果一开始b的值就是以小数的形式,则中间的计算过程误差会越来越大,你需要将0.00171053化成分数的形式
# t. N+ p  V$ }. {Rationalize[0.00171053, 10^-16]




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5