ODE Event Finding Accuracy

7 次查看（过去 30 天）
Grant Hecht
Grant Hecht 2021-2-19
评论： Grant Hecht 2021-2-19
I'm working on a project that requires that I determine when an ODE event occurs with a very high degree of accuracy. Looking at MATLAB's documentation, I've been unable to find any information on how to set the error tollerance for the root finding method or exactly how accurate the root finding method used for determining event locations is.
Is it somehow possible to set abs/rel tollerance on the root finding method that MATLAB uses to determine where ODE events occur? I've recently been using the Julia programming language with it's DifferentialEquations.jl package which allows me to explicitly define the absolute and relative error tolorance for the Regula-Falsi root finding method used to determine where events occur. Ideally, I would like to do this exact thing with MATLABs ODE event handling.
  1 个评论
Bjorn Gustavsson
Bjorn Gustavsson 2021-2-19
Have you checked that it isn't good enough? For ballode I get the event-time to be accurate with a relative error of 2.75e-14...

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Grant Hecht
Grant Hecht 2021-2-19
编辑：Grant Hecht 2021-2-19
I ended up answering my own question... digging into the source code a little I found that MATLAB uses Regula Falsi with a tollerance computed as 128*max(eps(t),eps(tnew)) to determine where an event occured. I wish MATLAB had a way of specifying the tollerance used in the root finding routine via odeset() so I could decrease it to around 10*eps(t) (which is what I used in Julia) but for the time being, I looks like I'm stuck with the default.
If anyone knows of a way to specify a different tollerence please let me know. I attempted to edit the souce code in odezero.m but realized the source code is unfortunately read only.
  2 个评论
Bjorn Gustavsson
Bjorn Gustavsson 2021-2-19
If you write your own odezero function with identical behaviour you can at least hard-code the tolerance you need, and if this function is found before the odezero of Mathworks then yours will be used. This is obviously system-programming-wise a very dubious thing to do since you'll have to manually control which function to use/have been used. The legalities I don't know about, but as long as you dont sell or redustribute and don't copy the interior of the original function bu only the input-output?
Grant Hecht
Grant Hecht 2021-2-19
I'll have to try this out! I just assumed that MATLAB's source code would take precidence over any user defined functions on MATLAB's path but thinking about it further, it wouldn't make much sence for MATLAB to behave this way. Thank you for the suggestion!

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

产品


版本

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by