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 编辑 ] 用Maple解的结果,等于5054
但是把方程中的0.0072换成0.0064,发现结果出不来了forimg.php?url=http://www.simwe.com/forum/images/smilies/cry.gif 算的结果是疲劳寿命,但是结果这么不一致,很头疼,方程到底应该怎么解??
大家帮帮忙 应该是精度误差判断引起的,如下:
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 确实如此,本来函数值就小,默认精度0.001根本不够,改成1e-9后差别就很小了
谢谢楼上 但是Maple为什么只列出方程,而不能计算结果呢? fsolve(0.76341e-2/n^.10448+.575/n^.6 = 0.6415e-2, n, complex)
5419.531242
页:
[1]