MathCAD和Maple解方程的问题！

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

用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，发现结果出不来了

登录/注册后可看大图

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

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

用Maple解的结果，等于5054
但是把方程中的0.0072换成0.0064，发现结果出不来了

登录/注册后可看大图

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

账号		自动登录	找回密码
密码			我要加入

[Mathcad] MathCAD和Maple解方程的问题！

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。