| 查看: 2247 | 回复: 10 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求教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的计算流程是什么样的呢？ |
回复此楼» 猜你喜欢
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有60人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有0人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
11楼2018-01-15 20:38:39
2楼2014-11-27 20:11:20
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68126.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
3楼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的变化情况。  方程时间序列.jpg |
4楼2014-11-27 21:20:06
10