| 查看: 1759 | 回复: 2 | ||
[求助]
关于龙格库塔法的求解时间步长设定问题 已有1人参与
|
新手编fortran程序，解一个微分方程，用龙格库塔法，关于时间步长，到底取多少合适，取多少才可以稳定呢，不知，请问怎么确定，有哪些教材里面有关于时间步长的取法解释呢？求高手相助，不甚感激。 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
找导师
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有11人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 如何求解多元二阶微分方程组
已经有7人回复
- 二阶龙格库塔λ是否都为非负？
已经有3人回复
- 帮忙：求解一组微分代数方程，奖励100金币
已经有7人回复
- 动力学方程组参数求解问题
已经有9人回复
- 求助！！对于调用龙格库塔法如何对结果判断反馈
已经有6人回复
- matlab解方程组求助攻
已经有4人回复
- Fluent求解瞬态问题时，中途改变时间步长，会不会对后来求解形成大的误差？
已经有7人回复
- matlab 非线性微分方程求解
已经有3人回复
- 简单四阶导数求解！！
已经有7人回复
- 拼了。。700金币求高手。最优控制，动态规划 matlab 仿真。。
已经有17人回复
- 6元方程求解----绘制等高线
已经有5人回复
- 求解一组四元一次方程
已经有4人回复
- Matlab用龙格库塔法求二阶非线性方程
已经有13人回复
- 【求助】非线性方程组的求解问题
已经有6人回复
- 分数阶微分方程的求解（用差分法，matlab编程求解）
已经有9人回复
- 龙格库塔法 matlab
已经有6人回复
- 如何求解描述振动的二阶微分方程
已经有21人回复
- matlab解微分方程组
已经有15人回复
- LABVIEW四阶龙格库塔法求解一阶微分方程组
已经有12人回复
- 【求助】急请微分方程高手看下这个方程组！
已经有14人回复
yongcailiu
金虫 (小有名气)
- 应助: 63 (初中生)
- 金币: 2283.9
- 红花: 10
- 帖子: 181
- 在线: 250.4小时
- 虫号: 1406563
- 注册: 2011-09-18
- 性别: GG
- 专业: 计算数学与科学工程计算
2楼2015-02-25 09:41:34
|
您好，谢谢您的答复，我是这么个方程，如下，打算用四级显示龙格库塔法去解，我看了您说的那个教材，但是因为数学基础差，看不懂，那个教材的表5.2说四级的绝对稳定域是-2.78<uh<0,这里h为时间步长，u为常数。请您帮我看看我这个方程到底时间步长取多少，哪里是我的u呢。 ∂T/∂t=m{0.94q+5.4999×【10】^(-8)∙[1-0.261∙exp(-0.74×【10】^(-4)× y^2 ) ]∙(1+0.017C^2 )∙【(273+y)】^4-5.4999×【10】^(-8)∙【(273+T)】^4-(4.324+0.2162【W】^2 )∙(T-y)-(9.2+0.46W^2 )∙{exp[2.3026∙(7.5T/(T+237.3)+0.6609)]- exp[2.3026∙(7.5d/(d+237.3)+0.6609)]}} 上述t为时间，T为我要求的变量， 其它参数都是常数。 |
3楼2015-02-25 17:07:10
20