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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 关于使用maltab中的龙格库塔方法解常微分方程组：
81/1返回列表
查看: 1832  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿A区211，22408 321求调剂 +7 随心所欲☆ 2026-04-15 8/400 2026-04-21 08:22 by Equinoxhua
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 一志愿中科大材料与化工，353分还有调剂学校吗 +11 否极泰来2026 2026-04-15 13/650 2026-04-20 22:31 by Equinoxhua
[考研] 320求调剂 +6 深郊akm 2026-04-17 6/300 2026-04-20 18:57 by fs26jie
[论文投稿] 期刊推荐 +3 材料研究生 2026-04-15 5/250 2026-04-20 16:02 by 豆豆7758
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 5/250 2026-04-20 10:47 by YuY66
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 291求调剂 +12 关忆北. 2026-04-14 13/650 2026-04-19 16:50 by 中豫男
[考研] 085404 22408 309分求调剂 +10 lzmk 2026-04-14 11/550 2026-04-19 16:42 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 297，工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
信息提示
关闭
请填处理意见
关闭