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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 求助！！对于调用龙格库塔法如何对结果判断反馈
71/1返回列表
查看: 1837  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

瀚海行舟

新虫 (初入文坛)

[求助] 求助！！对于调用龙格库塔法如何对结果判断反馈 已有2人参与

在matlab里调用龙格库塔算法，如何实现对每步的结果进行判断来选择不同的微分方程组？？具体情况如图：

求助！！对于调用龙格库塔法如何对结果判断反馈
360截图20140519204025928.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-05-20 15:10:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色

【答案】应助回帖

感谢参与，应助指数 +1
把你的两个方程组写成两个function，然后通过if判断来分别求解
赞一下 回复此楼
我同意叔本华的观点，人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活，逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
2楼2014-05-20 15:12:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

瀚海行舟

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by baobiao007 at 2014-05-20 15:12:24
把你的两个方程组写成两个function，然后通过if判断来分别求解

可是龙格库塔法是直接调用的，比如让算0到1之间的，在这个中间有判断，这时[T,X]=ode45('vdp1000',[0:0.001:1],[0,0,0])里面已经直接设定了起点终点，并不晓得从哪跳转啊？？（我表达能力有点差）
赞一下 回复此楼
3楼2014-05-20 15:19:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

reko34

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
瀚海行舟: 金币+5, ★★★很有帮助, 给出了思路 2014-06-10 15:48:59
龙哥库塔法我是只记得名字了，方程就看不懂了。程序上可以试试在调用命令的那行加个断点，然后step in进到调用的子程序，有些MATLAB的自带函数可以这样破解，复制个备份加你的判断。
赞一下(1人) 回复此楼
4楼2014-05-20 16:38:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

瀚海行舟

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by reko34 at 2014-05-20 16:38:51
龙哥库塔法我是只记得名字了，方程就看不懂了。程序上可以试试在调用命令的那行加个断点，然后step in进到调用的子程序，有些MATLAB的自带函数可以这样破解，复制个备份加你的判断。

非常感谢你的回答，作为一个初学者，对你的答案还不能完全明白，我先按你说的摸索一下，估计还会有问题，希望能够指教一二！！
赞一下 回复此楼
5楼2014-05-20 20:06:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

瀚海行舟

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by reko34 at 2014-05-20 16:38:51
龙哥库塔法我是只记得名字了，方程就看不懂了。程序上可以试试在调用命令的那行加个断点，然后step in进到调用的子程序，有些MATLAB的自带函数可以这样破解，复制个备份加你的判断。

你好，我按你说的还是没能搞定，我自己编了个简单函数，在程序里加判断，可是在判断的那个地方误差很大，能不能帮我分析下（是不是因为迭代的区间变了，该如何改呢？），奖励就给你了吧，看这样子是没人再回复了。具体如下：
function dy=vdp1000(t,y);
dy=zeros(1,1);
if y(1)<2
    dy(1)=1;
else
    dy(1)=2;
end
输入[T,X]=ode45('vdp1000',[0:0.001:4],[0]);
希望能帮帮忙
赞一下 回复此楼
6楼2014-05-21 20:33:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

reko34

木虫 (正式写手)
引用回帖:
6楼: Originally posted by 瀚海行舟 at 2014-05-21 20:33:55
你好，我按你说的还是没能搞定，我自己编了个简单函数，在程序里加判断，可是在判断的那个地方误差很大，能不能帮我分析下（是不是因为迭代的区间变了，该如何改呢？），奖励就给你了吧，看这样子是没人再回复了。 ...

龙哥库塔我没实际用过，学过就忘了，只记得是个迭代方法，在ode45里靠近结束的地方有个while循环嵌套，应该是主要的计算步骤，估计判断应该加在这里的循环内部。我理论不通所以看不懂哪步在干什么，还需要你自己分析。
赞一下(1人) 回复此楼
7楼2014-05-21 20:53:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 瀚海行舟 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学 +31 妮妮ninicgb 2026-04-04 35/1750 2026-04-09 21:06 by zhouxiaoyu
[考研] 085601初试330分找调剂 +8 流心奶黄包l 2026-04-09 8/400 2026-04-09 19:19 by 5268321
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 327求调剂 +12 Xxjc1107. 2026-04-06 12/600 2026-04-08 16:46 by luoyongfeng
[考博] 申博 +8 IQwQl 2026-04-04 8/400 2026-04-08 09:43 by 0608104024
[考研] 263分B区求调剂 +6 李nihao 2026-04-08 6/300 2026-04-08 09:38 by 南开小綦
[考研] 338求调剂 +5 小猪红色 678 2026-04-06 6/300 2026-04-07 21:18 by 乔哒哒哒
[考研] 材料调剂 +11 一样YWY 2026-04-07 11/550 2026-04-07 15:13 by shdgaomin
[考研] 材料调剂 +17 小刘同学吖吖 2026-04-06 18/900 2026-04-07 11:41 by 诗与自由
[考研] 一志愿太原理工大学计算机技术专硕348，求调剂指导 +3 nexious 2026-04-05 3/150 2026-04-07 08:19 by jp9609
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
[考研] 调剂 一志愿吉林大学357分 +5 .Starry. 2026-04-04 5/250 2026-04-06 09:28 by cql1109
[考研] 308求调剂 +4 maverick^_^ 2026-04-03 4/200 2026-04-05 19:08 by 蓝云思雨
[考研] 0860 求调剂 一志愿国科大 348 分 +3 WiiiP 2026-04-03 3/150 2026-04-05 17:43 by Ecowxq666！
[考研] 319求调剂 +4 星星不眨眼喽 2026-04-03 4/200 2026-04-04 16:25 by 中飞院空管学院�
[考研] 309求调剂 +4 快乐的小白鸽 2026-04-04 5/250 2026-04-04 15:55 by cql1109
[考研] 一志愿中国石油大学化学工程323分求调剂 +4 化工专硕323分 2026-04-03 6/300 2026-04-03 22:12 by dongzh2009
[考研] 315求调剂 +6 顺理成张 2026-04-03 8/400 2026-04-03 14:04 by 百灵童888
[考研] 调剂 +7 祉岷. 2026-04-02 7/350 2026-04-03 09:11 by 花呗还欠600
信息提示
关闭
请填处理意见
关闭