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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 关于使用maltab中的龙格库塔方法解常微分方程组：
81/1返回列表
查看: 1820  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 307求调剂 +10 冷笙123 2026-03-17 10/500 2026-03-21 01:54 by JourneyLucky
[考研] 297求调剂 +9 戏精丹丹丹 2026-03-17 9/450 2026-03-21 01:49 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 299求调剂 +4 shxchem 2026-03-20 5/250 2026-03-21 00:24 by 星空星月
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 265求调剂 +12 梁梁校校 2026-03-19 13/650 2026-03-20 21:01 by 无际的草原
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 广西大学家禽遗传育种课题组2026年硕士招生（接收计算机专业调剂） +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 生物学071000 329分求调剂 +3 我爱生物生物爱� 2026-03-17 3/150 2026-03-18 10:12 by macy2011
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 考研化学学硕调剂，一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 11408 一志愿西电，277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
信息提示
关闭
请填处理意见
关闭