| 查看: 2156 | 回复: 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人回复
shikang999
新虫 (小有名气)
- 计算强帖: 1
- 应助: 27 (小学生)
- 金币: 1454.7
- 红花: 7
- 帖子: 182
- 在线: 206.4小时
- 虫号: 2201981
- 注册: 2012-12-23
- 性别: GG
- 专业: 固体力学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
brucexmc: 金币+30, ★★★很有帮助 2014-11-29 19:17:09
感谢参与，应助指数 +1
brucexmc: 金币+30, ★★★很有帮助 2014-11-29 19:17:09
|
1、ode45计算流程数值算法书上已经有，这里不详细叙述 2、你函数里记录的t和你外面输出的那个t不存在线性相关性（特别是一些算法做了技巧性处理，比如加入自适应或者复化分段自控误差等等），所以不影响你外面t按大小输出 3、看你这个数据记录情况，很可能t在1附近时，内部计算误差比较大，为了控制误差技巧性地自适应了一个区间（区间前后扩移了一下），进而得到相应点的值 4、当然，我并不了解ode45内部技巧性处理算法，仅是凭自己写的经验来判断！ 5、希望能帮到你！ |
» 本帖已获得的红花（最新10朵）
brucexmc
9楼2014-11-29 13:51:10