数学建模社区-数学中国
标题:
请教一个解方程组的问题
[打印本页]
作者:
prince12th
时间:
2012-8-27 17:00
标题:
请教一个解方程组的问题
用solve语句解这样一个方程组
8 [6 O8 W2 X( p) V' k
(x1 - a1)^2 == b^2*(3*x1^6 + 6*x1^2*x2^4)
, G: }' x2 Z; z- e6 i1 p
(x2 - a2)^2 == b^2*(3*x2^6 + 6*x2^2*x1^4)
2 {8 ^! ]' ] F8 T. n1 F; k9 N) H- w
a1,a2已知，b是0.00171053
- n/ |: _$ g! Q# n# D
, O1 s8 a, V. X' k, [- f* h
当设a1,a2都为1时，结果出来的是
. o9 @3 Q$ b# A2 k. n9 J
{{{x1 -> 1373.46, x2 -> 13.4299},
! T! e- o4 g" p. b% ~9 C! e
{x1 -> 10987.7, x2 -> 15.3663},
1 `2 Q I; O: x4 F; p
{x1 -> 1.05482*10^6, x2 -> 17.8499},
3 S/ P6 g3 z+ G+ g1 f% s$ a
{x1 -> 2.81284*10^6, x2 -> 17.8426}}}
6 F( u( K% U; ~8 s* E) k2 q
1 s- e% t( X% s) K1 K8 W# O, X
x1,x2不是应该对称的么，而且把结果代回原来的式子里也不正确，为什么呢？
作者:
sycxyc
时间:
2012-8-27 18:39
有小数点的计算精度是有限的，并且高次方程没有精确解，如果一开始b的值就是以小数的形式，则中间的计算过程误差会越来越大，你需要将0.00171053化成分数的形式
1 ?* P2 b0 @2 D" Q& M0 R4 X" @
Rationalize[0.00171053, 10^-16]
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5