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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求教ode45计算流程问题
51/1返回列表
查看: 2147  |  回复: 10
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

brucexmc

新虫 (正式写手)

[求助] 求教ode45计算流程问题已有2人参与

大家好，本人最近利用Matlab的ODE45函数求解微分方程[t,y]=ode45(@fun,tspan,t0)，在fun函数中用写文件命令fsprint()记录了一些需要使用的参数，最后发现ODE45计算过程中时间项t不是单调增加，而是有时候会往回走一段儿，比如[0, 0.1, 0.2, 0.3, 0.25, 0.28, 0.3, 0.35, 0.4]，但是最后输出的[t,y]中的t确实单调增大的。
所以想问一问各位有没有谁比较懂这个，ode45的计算流程是什么样的呢？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-11-27 19:32:20
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

brucexmc

新虫 (正式写手)
送红花一朵
引用回帖:
9楼: Originally posted by shikang999 at 2014-11-29 13:51:10
1、ode45计算流程数值算法书上已经有，这里不详细叙述
2、你函数里记录的t和你外面输出的那个t不存在线性相关性（特别是一些算法做了技巧性处理，比如加入自适应或者复化分段自控误差等等），所以不影响你外面t按大 ...

非常感谢您细心的回复，我也在研究ode45的时间步长调整方法，过程比较复杂，目前还没有完全理解。但你说的很有道理，是其算法自身的一些技巧导致了这样的现象，因为在t=1时，二阶导数不连续了，可能会导致计算误差增大的现象。
赞一下 回复此楼
10楼2014-11-29 19:15:54
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

brucexmc

新虫 (正式写手)
求大神指教啊，给自己顶一下下吧！
赞一下 回复此楼
2楼2014-11-27 20:11:20
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
你说的这种情况比较少见。
麻烦给出完整代码。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2014-11-27 20:29:38
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

brucexmc

新虫 (正式写手)
引用回帖:
3楼: Originally posted by 月只蓝 at 2014-11-27 20:29:38
你说的这种情况比较少见。
麻烦给出完整代码。

非常感谢您的回复，我的调用函数代码为：
clear all;
clc;
options=odeset('reltol',1e-5);
[t,y]=ode45(@fun,[0 2],[1 -1],options);
plot(t,y(:,1),t,y(:,2));

函数fun的定义为
function dx=fun(t,x)

if t<1
    a=0.1;
else
    a=0.2;
end

%数据写入文件
fdata=fopen('data.txt','a');
fprintf(fdata,'%f\n',t);
fclose(fdata);

dx=[x(2);1+a*x(2)];

经过计算之后，在t=1左右就会出现t值先变小，之后又变大的现象，下图所示为记录得到的t的变化情况。
求教ode45计算流程问题
方程时间序列.jpg
赞一下 回复此楼
4楼2014-11-27 21:20:06
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答
信息提示
关闭
请填处理意见
关闭