| 查看: 2146 | 回复: 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的计算流程是什么样的呢？ |
回复此楼» 猜你喜欢
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有285人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有19人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.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
- 金币: 68120.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