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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 关于使用maltab中的龙格库塔方法解常微分方程组：
汕头大学海洋科学接受调剂
81/1返回列表
查看: 1828  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wangdexina

金虫 (正式写手)

[求助] 关于使用maltab中的龙格库塔方法解常微分方程组：

关于使用maltab中的龙格库塔方法解常微分方程组：
例如下面这个问题，这是个方程组，是分开求解还是用一个function呢？
目前我用的是一个function，可是在求解x的二阶导数的时候会调用y的二阶导数，这个真的没遇见过。希望大神们可以给予帮助~！谢谢

关于使用maltab中的龙格库塔方法解常微分方程组：
wt.jpg
回复此楼

» 猜你喜欢
1楼 2017-04-22 17:12:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cmaik
顶一下

发自小木虫Android客户端

» 本帖已获得的红花（最新10朵）

wangdexina
2楼2017-04-22 17:22:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangdexina

金虫 (正式写手)
送红花一朵
引用回帖:
2楼: Originally posted by cmaik at 2017-04-22 17:22:41
顶一下

谢谢
回复此楼
3楼2017-04-22 17:23:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

somomo91

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
wangdexina(月只蓝代发): 金币+15, 感谢应助 2017-04-23 10:58:16
wangdexina: 金币+20, ★★★★★最佳答案 2017-04-24 09:00:04
这个方程组，可以简化为 四个一介方程组：解高阶方程组都是用这个 降介 以及 分解的方法
然后直接用 ode45 就可以
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

wangdexina
4楼2017-04-22 17:45:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangdexina

金虫 (正式写手)
送红花一朵
引用回帖:
4楼: Originally posted by somomo91 at 2017-04-22 17:45:50
这个方程组，可以简化为 四个一介方程组：解高阶方程组都是用这个 降介 以及 分解的方法
然后直接用 ode45 就可以

首先非常感谢你的回答，其次，我就是降阶之后做的，具体如下图
我的做法是将dy(2)在方程外用符号表达，例如YY=c*y(1)+d*y(2)
程序中的dy(2)=YY;
在dy(4)中dy(4)=a*y(3)+(b+YY)*y(4)
不知道这么做对不对，但是算出来的结果和文献上不一样
关于使用maltab中的龙格库塔方法解常微分方程组：-1
wt2.jpg
赞一下 回复此楼
5楼2017-04-22 18:39:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

somomo91

专家顾问 (职业作家)
引用回帖:
5楼: Originally posted by wangdexina at 2017-04-22 18:39:49
首先非常感谢你的回答，其次，我就是降阶之后做的，具体如下图
我的做法是将dy(2)在方程外用符号表达，例如YY=c*y(1)+d*y(2)
程序中的dy(2)=YY;
在dy(4)中dy(4)=a*y(3)+(b+YY)*y(4)
不知道这么做对不对，但是算 ...

你图片里不是很清楚，
正常
dy(1)=y(2) ;
dy(2)=c*y(1)+d*y(2) ;
dy(3)=y(4) ;
dy(4)=a*y(3)+(b+ dy(2))*y(4) ;
结果之所以对不上，是因为这里有好几个因素：
1. 时间的步长
2. 二阶导 中 的 常数 dy = y+c，这里有好几个地方都可以加 c 常数
具体得看文献中给出的数据够不够具体
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

wangdexina
6楼2017-04-22 19:14:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangdexina

金虫 (正式写手)
送红花一朵
引用回帖:
6楼: Originally posted by somomo91 at 2017-04-22 19:14:37
你图片里不是很清楚，
正常
dy(1)=y(2) ;
dy(2)=c*y(1)+d*y(2) ;
dy(3)=y(4) ;
dy(4)=a*y(3)+(b+ dy(2))*y(4) ;
结果之所以对不上，是因为这里有好几个因素：
1. 时间的步长
2. 二阶导 中 的 常数 dy = y+ ...

我的思路是没有错的吗？
我总感觉是我的做法有错误呢？
如果没有错，那就是文献中给的数据不充分？
我在看看，仔细检查一下，谢谢
赞一下 回复此楼
7楼2017-04-22 19:38:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangdexina

金虫 (正式写手)
引用回帖:
6楼: Originally posted by somomo91 at 2017-04-22 19:14:37
你图片里不是很清楚，
正常
dy(1)=y(2) ;
dy(2)=c*y(1)+d*y(2) ;
dy(3)=y(4) ;
dy(4)=a*y(3)+(b+ dy(2))*y(4) ;
结果之所以对不上，是因为这里有好几个因素：
1. 时间的步长
2. 二阶导 中 的 常数 dy = y+ ...

我能再问一个问题吗？
如果是这样的方程呢：
dy(1)=y(2);
dy(2)=a*y(1)+(b+dy(4))*y(2);
dy(3)=y(4);
dy(4)=c*y(3)+(d+dy(2))*y(4);
这样的方程，两个方程中互有对方的二阶导数，
不需要相减，化简约分吗？
赞一下 回复此楼
8楼2017-04-24 19:39:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangdexina 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +12 璃茉一定上岸 2026-04-10 13/650 2026-04-14 00:08 by Equinoxhua
[考研] 材料相关专业344求调剂双非工科学校或课题组 +17 hualkop 2026-04-12 18/900 2026-04-13 22:42 by pies112
[考研] 271求调剂 +24 2261744733 2026-04-11 24/1200 2026-04-13 22:24 by pies112
[考研] 调剂 +27 不逢春 2026-04-07 28/1400 2026-04-13 20:19 by xxy88888
[考研] 一志愿双非085400电子信息344 求调剂，对材料和化学方向也感兴趣 +12 无情的小羊 2026-04-09 13/650 2026-04-13 14:17 by 张zhihao
[考研] 化工学硕 285求调剂 +29 Wisjxn 2026-04-07 29/1450 2026-04-13 13:47 by floriea
[考研] 本科南方医科大学 一志愿985 药学学硕284分 求调剂 +5 弱水听文 2026-04-09 5/250 2026-04-12 13:16 by shengxi123
[考研] 调剂求助 +6 果然有我 2026-04-11 7/350 2026-04-11 16:22 by 明月此时有
[考研] 298求调剂 +9 钉叮咚冬瓜 2026-04-07 11/550 2026-04-11 09:35 by zhq0425
[考研] 284求调剂 +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[考研] 化学工程与技术324调剂 +23 孙常华 2026-04-09 25/1250 2026-04-11 00:07 by 骑牛渡寒江
[考研] 一志愿北理工298英一数二已上岸，感谢各位老师 +14 Reframe 2026-04-10 16/800 2026-04-10 23:07 by caotw2020
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 计算机类求调剂，22408-274分 +7 上岸de小虫 2026-04-09 8/400 2026-04-10 19:56 by fxue1114
[考研] 293调剂 +25 yj1221 2026-04-08 26/1300 2026-04-10 15:02 by 柴小白
[考研] 一志愿沪9，326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[考研] 070300化学 求调剂 +13 73372112 2026-04-08 13/650 2026-04-09 20:22 by maddjdld
[考研] 材料工程322 +18 哈哈哈吼吼吼哈 2026-04-07 19/950 2026-04-09 10:44 by cymywx
[考研] 生物学学硕，初试351分，求调剂 +4 …～、王…～ 2026-04-08 5/250 2026-04-08 21:49 by limeifeng
[考研] 259求调剂 +5 就爱吃土豆呀呀 2026-04-07 5/250 2026-04-07 22:40 by JourneyLucky
信息提示
关闭
请填处理意见
关闭