| 查看: 2927 | 回复: 2 | ||
[求助]
mathematica 使用NDSolve时出现的问题 已有1人参与
|
| 在mathematica 9.0中，使用NDSolve 时，只是将边界条件中的0改为0.03后，为什么会出现图中的问题呢？ 在 x == 0.00116` 处，步长实际上为零；可能存在奇点或者刚性系统。怎么解决？？ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : mathematica.png
2014-04-12 17:22:45, 34.12 K
» 猜你喜欢
窗边初夏的小雨
已经有10人回复
-
护理论文 晋升
已经有4人回复
-
西南大学考核制博士
已经有5人回复
-
求碳排放博导；方向是LCA、生命周期可持续发展以及碳排放
已经有8人回复
-
评审感受-评审感受-评审感受
已经有11人回复
-
26年申博自荐-计算机视觉
已经有5人回复
-
导师各种操作恶心咋办
已经有12人回复
-
现在不知道怎么办，感觉很痛苦
已经有5人回复
-
生物学博士
已经有6人回复
-
2026年申博-电池方向
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程的求解
已经有9人回复
- 请教mathematic问题？关于解微分方程组的和画图的一些问题。
已经有19人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kexiaolong: 金币+30 2014-12-11 09:22:01
kexiaolong: 金币+30 2014-12-11 09:22:01
|
0.03的问题我没仔细研究，但是从你的提问来看，你是因为0处失败了才转向了0.03吧？那么，对0处的失败，可以添加选项MaxSteps->Infinity来解决。不过，你确定你这方程没问题吗？： eqn = With[{i = 0}, Rationalize[{y''[x] - y[x]/0.00008585917334 + 1/30.02 Exp[38 + 32.633/(2.586 + 138 x + 25.655 x^0.732)] == 0, y' == 100/30.02 y, y[11.6 10^-4] == 0}, 0]]; sol = NDSolve[eqn, y, {x, 0, 400}, MaxSteps -> Infinity, WorkingPrecision -> 16]; y[100] /. sol (* {-9.6701192450*10^4701} *) 调了下精度发现冒似没什么改善……直觉上你最好再分析下你的方程。 |
2楼2014-07-03 18:51:18
3楼2014-07-03 18:52:49