24小时热门版块排行榜    

查看: 226  |  回复: 2
当前主题已经存档。

PPCAI

至尊木虫 (著名写手)

小辣椒

[交流] 【求助】关于matlab计算终止条件设定

利用matlab的ODE工具计算二介微分方程,形式如下:
y1'=y2
y2'=f(t, y1,y2)
现在用的命令是:
Tspan = [0 0.02]; %时间的积分范围
Y0=[R0 0]; %初始值,即t=0, y1=R0, y2=0
options=odeset('RelTol',1e-4,'AbsTol',1e-4);
[T,Y]=ode45(@fun,Tspan,Y0,options);

现在想设定一个计算终止条件, 想在y1/R0<0.1时结束计算, 不知道该如何处理, 是通过odeset中的事件属性设定吗?如果是,该如何设定, 另外是否有其它实现方法, 请高手指点!不甚感激!
如果问题得到解决, 将给予5个金币的酬谢

[ Last edited by PPCAI on 2009-3-25 at 16:09 ]
回复此楼
爱工作!爱生活!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

windhair

铜虫 (小有名气)


sunxiao(金币+1,VIP+0):谢谢参与交流,欢迎常来仿真编程版 3-27 00:10
这个好象不行吧,因为方程的求解是由 RelTol 控制的,没法自己控制
2楼2009-03-26 17:16:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

PPCAI

至尊木虫 (著名写手)

小辣椒

应该是可以强制结束吧, 我看见好多这样的类似的处理方式
但是我照着做,却总是出错.
引用回帖:
Originally posted by windhair at 2009-3-26 17:16:
这个好象不行吧,因为方程的求解是由 RelTol 控制的,没法自己控制

爱工作!爱生活!
3楼2009-03-27 08:46:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 PPCAI 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见