24小时热门版块排行榜    

查看: 1679  |  回复: 2

康雅思

新虫 (初入文坛)

[求助] 关于龙格库塔法的求解时间步长设定问题 已有1人参与

新手编fortran程序,解一个微分方程,用龙格库塔法,关于时间步长,到底取多少合适,取多少才可以稳定呢,不知,请问怎么确定,有哪些教材里面有关于时间步长的取法解释呢?求高手相助,不甚感激。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongcailiu

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
详见李荣华等著的《微分方程数值解法》(第四版)中 P43 表5.2给出了绝对收敛域的取值范围,它和龙格库塔法的级数有关
2楼2015-02-25 09:41:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

康雅思

新虫 (初入文坛)

您好,谢谢您的答复,我是这么个方程,如下,打算用四级显示龙格库塔法去解,我看了您说的那个教材,但是因为数学基础差,看不懂,那个教材的表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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 康雅思 的主题更新
信息提示
请填处理意见