| 查看: 2909 | 回复: 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
» 猜你喜欢
遇见不省心的家人很难过
已经有16人回复
-
退学或坚持读
已经有25人回复
-
博士延得我，科研能力直往上蹿
已经有4人回复
-
免疫学博士有名额，速联系
已经有14人回复
-
面上基金申报没有其他的参与者成吗
已经有4人回复
-
多组分精馏求助
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程的求解
已经有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