stargazer 发表于 2007-1-22 23:47

MathCAD和Maple解方程的问题!

用MathCAD解方程,调用root函数,发现设不同初始值,得到结果差异很大,到底哪一个才对阿,抓狂!!!

令X0=5000,解得 x=5054;
令X0=6000,解得 x=6420;
令X0=8000,解得 x=3833;
令X0=10000,解得 x=6420;
令X0=12000,解得 x=5054;

用Maple解的结果,等于5054
但是把方程中的0.0072换成0.0064,发现结果出不来了forimg.php?url=http://forimg.php?url=http://forimg.php?url=http://forimg.php?url=http://www.simwe.com/forum/images/smilies/cry.gif



[ 本帖最后由 stargazer 于 2007-1-22 23:49 编辑 ]

stargazer 发表于 2007-1-22 23:50

用Maple解的结果,等于5054
但是把方程中的0.0072换成0.0064,发现结果出不来了forimg.php?url=http://www.simwe.com/forum/images/smilies/cry.gif

stargazer 发表于 2007-1-22 23:53

算的结果是疲劳寿命,但是结果这么不一致,很头疼,方程到底应该怎么解??
大家帮帮忙

dingd 发表于 2007-1-23 10:19

应该是精度误差判断引起的,如下:
x               f(x)
2527        5.5374201E-6
3285.3        1.8366879E-6
4043.6        4.3810763E-7
4801.9        2.0968216E-8
5560.2        6.7445897E-8
6318.5        3.4482538E-7
7076.8        7.3853063E-7
7835.1        1.1884995E-6
8593.4        1.6618859E-6
9351.7        2.1402621E-6
10110        2.6131922E-6

f(x)的变化很小,如终止判断误差取为1E-6左右,就可的任意多的解。实际上只有一个最优解。下面是用1stOpt算的:
1:Function 0.007634063*x^(-0.10448)+0.685*x^(-0.6)-0.007238722;
结果:x: 5053.95687520354

2:Function 0.007634063*x^(-0.10448)+0.685*x^(-0.6)-0.006415;
结果:x: 6969.7310450444

stargazer 发表于 2007-1-23 17:42

确实如此,本来函数值就小,默认精度0.001根本不够,改成1e-9后差别就很小了
谢谢楼上

stargazer 发表于 2007-1-23 17:43

但是Maple为什么只列出方程,而不能计算结果呢?

suffer 发表于 2007-2-4 18:36

fsolve(0.76341e-2/n^.10448+.575/n^.6 = 0.6415e-2, n, complex)

5419.531242
页: [1]
查看完整版本: MathCAD和Maple解方程的问题!