24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1748  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

瀚海行舟

新虫 (初入文坛)

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

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

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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的回帖
查看全部 7 个回答

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的自带函数可以这样破解,复制个备份加你的判断。
4楼2014-05-20 16:38:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见