| 查看: 2611 | 回复: 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的计算流程是什么样的呢？ |
回复此楼» 猜你喜欢
青基已中
已经有3人回复
-
散金币祈福
已经有186人回复
-
物理学I论文润色/翻译怎么收费?
已经有113人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有110人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有66人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于含量计算的公式问题
已经有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、希望能帮到你！ |
brucexmc
9楼2014-11-29 13:51:10
2楼2014-11-27 20:11:20
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68175.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.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
