24小时热门版块排行榜    

查看: 1355  |  回复: 5

ayi1021

新虫 (初入文坛)

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

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

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

郑州苗鹏

金虫 (小有名气)

百度上有,可以用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的回帖

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的回帖

taylor66

新虫 (初入文坛)

号高深
6楼2015-12-25 15:17:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ayi1021 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 17:01 by BCdvBAEpt3zm
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 16:58 by BCdvBAEpt3zm
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-01 16:48 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:45 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 14/700 2026-11-01 16:32 by BCdvBAEpt3zm
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-01 15:04 by BCdvBAEpt3zm
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 14:55 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 14:24 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 12:38 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:53 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 X5bHxOULxFqb 2026-10-30 26/1300 2026-11-01 10:22 by BCdvBAEpt3zm
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
请填处理意见