ODE Events In SimBiology toolbox

Question

Alexander 2019-4-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/453687-ode-events-in-simbiology-toolbox

评论： Alexander 2019-4-1

采纳的回答： Arthur Goldsipe

I would like to define a custom ode stopping event for the solver of a SimBiology model. How to achieve that?

Say, I want to write some event function, like described in

https://de.mathworks.com/help/matlab/math/ode-event-location.html

and set isterminal conditioned at specific values of my concentrations, the model of which is defined by the SimBiology toolbox.

Where is the option located (or, how to export the code in the needed granularity) where I can provide my event function?

2 个评论
显示无隐藏无

zmi zmi 2019-4-1

You might be interested in "addevent" method of the model object and "trigger" property of the event object in SimBiology.

Alexander 2019-4-1

I'm aware of addevent. But it seems, that there is no opportunity to stop simulation by this event handler. At least, I didn't find any documentation about something analogous to the "isterminal" setting.

请先登录，再进行评论。

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

Answer 1

Arthur Goldsipe 2019-4-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/453687-ode-events-in-simbiology-toolbox#answer_368449

Currently (as of R2019a), SimBiology events cannot directly stop a simulation. I usually handle this sort of situation by post-processing simulation results that continue past the event of interest.

That said, there are 3 stop criteria: simulation stop time, maximum number of logs, and maximum (wall clock) simulation time. If you don't care about how long a simulation takes, you could use an event indirectly end the simulation by triggering the wall clock condition. For example, if you set the maximum wall clock time to 30 seconds and create an event that pauses for 30 seconds, then the simulation should stop after the event is triggered.

You also asked about export options. SimBiology supports export to SBML. You could also use the output from getequations as a starting point for converting your model to MATLAB code that can be used with ode15s or one of the other solvers in MATLAB.