应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 4阶龙格库塔方法疑问
71/1返回列表
查看: 1301  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

人民海军

木虫 (职业作家)

[求助] 4阶龙格库塔方法疑问

本人在做应用问题的时候遇到一个常微分方程组，打算用4阶R-K方法求数值解。非常奇怪的是，当我把迭代步长设为某个值的时候，得到的数值解都是正的（由于是实际问题，必须为正），然而当我把步长增加或者减小的时候，得到的数值解结果却会出现负值，求高人解释！！！谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
Letbygonesbebygones.
1楼 2012-07-16 21:23:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aborna

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
人民海军: 金币+10, ★有帮助 2012-07-17 11:42:01
这个可能与你的方程有关，比方说你的方程有刚性。这样就会造成很大的误差
赞一下 回复此楼
2楼2012-07-16 21:30:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

人民海军

木虫 (职业作家)
引用回帖:
2楼: Originally posted by aborna at 2012-07-16 21:30:51
这个可能与你的方程有关，比方说你的方程有刚性。这样就会造成很大的误差

你的意思是我所得到的全部为正的数值解实际上是误差太大的缘故么？
赞一下 回复此楼
Letbygonesbebygones.
3楼2012-07-16 23:21:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aborna

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
人民海军: 金币+5, ★有帮助 2012-07-18 11:36:34
不太清楚你的步长小的时候有多小，过小的话，可能累计误差太大。不敢说你的正数值解是否正确。建议你用多种算法试验下
赞一下 回复此楼
4楼2012-07-17 08:16:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lwingkei

禁虫 (小有名气)
本帖内容被屏蔽
5楼2012-07-19 09:45:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

人民海军

木虫 (职业作家)
引用回帖:
5楼: Originally posted by lwingkei at 2012-07-19 09:45:15
你改变时间步的时候同时记得改下网格，时间步跟网格粗细有关系。

网格？是常微分方程组，不是偏微分方程组！！！你的意思是？？
赞一下 回复此楼
Letbygonesbebygones.
6楼2012-07-19 11:20:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lwingkei

禁虫 (小有名气)
本帖内容被屏蔽
7楼2012-07-19 16:17:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 人民海军 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 6/300 2026-08-31 01:48 by hZiFeudZyoGR
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-31 01:27 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-31 00:45 by hZiFeudZyoGR
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-31 00:32 by hZiFeudZyoGR
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 22:24 by xiangy672
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:30 by hZiFeudZyoGR
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 zICmwzsBXjbN 2026-08-29 5/250 2026-08-30 20:08 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:57 by hZiFeudZyoGR
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:35 by hZiFeudZyoGR
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 19:01 by hZiFeudZyoGR
[博后之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 18:28 by hZiFeudZyoGR
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 17:36 by hZiFeudZyoGR
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +8 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 17:11 by hZiFeudZyoGR
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
关闭
请填处理意见
关闭