从deSolve R-package中使用ode()函数时,在每个集成步骤中提取本地截断错误(LTE)的值

时间:2019-04-09 14:31:19

标签: r ode numerical-integration lte

deSolve R包中的ode()函数使用的默认数字方案是lsode方法,该方法实现BDF和隐式Adams线性多步方案来求解ODE系统。集成以可变步长完成,该步长通过估计每个步骤的本地截断误差(LTE)来控制。实际上,以使得LTE的估计小于某个预设值rtol和atol的方式来调整步长。 ode方法中的默认值为:

defaultValue

并且可以适应。

但是,实际评估的LTE总是高于这些值。有没有办法从求解器中提取此值?

下面是一个示例系统:

ode(rtol = 1e-6, atol = 1e-6)

0 个答案:

没有答案