Stopping simulation early on event

9 visualizaciones (últimos 30 días)
Abed Alnaif
Abed Alnaif el 25 de Nov. de 2019
Comentada: Abed Alnaif el 3 de Dic. de 2019
Hello,
Is there a good way to stop a simulation early if an event occurs; i.e., stop the simulation as soon as the event is triggered? Also, I'd like to do this in a way that causes the objective function value to be infinite when using 'sbiofit'.
Thank you,
Abed

Respuesta aceptada

Arthur Goldsipe
Arthur Goldsipe el 26 de Nov. de 2019
Hi Abed,
Interesting question. The only way an event can stop a simulation early is if the event errors. sbiofit would handle such a siutation by returning Inf for the objective function value. Since you mentioned fitting, I generally recommend setting MaximumWallClock during fitting to ensure that things don't get stuck if a fit encounters a bad parameter value that results in a really long running simulation.
Ok, back to a simulation that errors. The only tricky thing with that is creating an event that is compatible with SimBiology's acceleration (that is, with the subset of the MATLAB language that can be used for code generation). Here's an example of modifying the Lotka-Volterra model to end early due to an event:
First, save the following function to myfunction.m:
function time = myfunction(time)
if time >= 0
error("Ending simulation at time " + time);
end
Then, execute the following:
sbioloadproject('lotka');
threshold = m1.addparameter('threshold',1e6);
m1.addevent('y1 > threshold', 'z = myfunction(time)');
sbioaccelerate(m1); % Ensure this model can be accelerated
simdata = sbiosimulate(m1); % Runs to completion
threshold.Value = 1e3;
simdata = sbiosimulate(m1); % Simulation errors

Más respuestas (0)

Comunidades de usuarios

Más respuestas en  SimBiology Community

Categorías

Más información sobre Scan Parameter Ranges en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

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

Start Hunting!

Translated by