| 查看: 2148 | 回复: 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论文润色/翻译怎么收费?
已经有206人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有19人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
关于含量计算的公式问题
已经有6人回复- 在线求教 计算表达式化简
已经有21人回复
- 求教关于CCSD计算单点能问题
已经有19人回复
- 求教沼气集中供气工程工艺设计各部分的具体计算过程
已经有9人回复
- 求教LCALCPOL计算极化的结果怎么处理的，多出来一行
已经有4人回复
- 求教个ASPEN PLUS的问题！
已经有5人回复
- 高金币求教高手用FCclasses计算FC因子
已经有11人回复
- 求教一道关于ΔH和ΔS计算的物理化学问题
已经有8人回复
- 求教matlab线性优化求最小值的问题
已经有7人回复
- 求教吉布斯函数关于平衡的计算问题
已经有10人回复
- QE-PWscf计算声子谱流程。
已经有10人回复
- 求教关于锂电池rate capacity的计算问题
已经有3人回复
- 求教分子轨道系数的计算
已经有3人回复
- 求45号钢加热过程的对流换热系数
已经有4人回复
- 求教一个小白问题
已经有10人回复
- 求教Gaussian用TDDFT计算单重态多重态的问题
已经有6人回复
- 在做fluent例子的时候遇到了问题，发现无法进行迭代计算，求教大家
已经有9人回复
- 【求助】求教DPD高手怎样用amorphous cell计算Flory huggins参数
已经有8人回复
- 【求助】MS计算CPU占用率问题，求教
已经有17人回复
8楼2014-11-28 09:15:41
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