24小时热门版块排行榜    

查看: 1291  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ayi1021

新虫 (初入文坛)

[求助] 龙格-库塔法解三阶微分方程的步骤 已有2人参与

近日遇到一个三阶微分方程需要求解,所以想要了解怎么用龙格库塔方法求解,可以随便举个例子说明,我主要是想知道四阶龙格库塔公式里面的f(x,y)怎么确定。

龙格-库塔法解三阶微分方程的步骤
25](7ALX9@FT1CDCXQLB_F7.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

例如:
y''(x)=f(x,y),y(x0)=a,y'(x0)=b.
改为:
y'(x)=z(x)
z'=f(x,y)
下面用显式Euler方法计算:
y0=a;
z0=b;
y1=y0+h*z0
z1=z0+h*f(x0,y0)
y2=z1+h*z1
z2=z1+h*f(x1,y1)
.......
h为步长
y1,y2....就是一y(x1),y(x2)....的数值解!
善恶到头终有报,人间正道是沧桑.
5楼2015-12-24 17:43:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

郑州苗鹏

金虫 (小有名气)

百度上有,可以用MATLAB中ode45求解!

发自小木虫Android客户端
天行健君子以自强不息
2楼2015-10-14 12:27:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ayi1021

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 郑州苗鹏 at 2015-10-14 12:27:44
百度上有,可以用MATLAB中ode45求解!

据说是把高阶变低阶,你能举个例子说明一下吗?

发自小木虫Android客户端
3楼2015-10-16 10:08:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伫立的灯塔

铁虫 (正式写手)

【答案】应助回帖

z=[y y' y''],z(1)=y,z'(1)=z(2)=y',z'(2)=z(3)=y'',z'(3)=y'''=表达式,把未知量化为一阶变量
4楼2015-12-24 16:48:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 14:24 by jCd0dEvKHShX
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 14:15 by jCd0dEvKHShX
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 14:03 by jCd0dEvKHShX
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:02 by jCd0dEvKHShX
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:51 by jCd0dEvKHShX
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 11:49 by jCd0dEvKHShX
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +14 晴天加油 2026-08-26 15/750 2026-08-29 11:06 by biily
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
信息提示
请填处理意见