How to handle stiffness parameter using Try and Catch?

1 view (last 30 days)
Hi all,
I'm trying to do parameter estimation (using GA) in my stiffness 1D ODE system.
When GA automatically try various parameter values to fit data, some parameter cases are so sensitive that main function were shut down due to stiffness problem.
In this case, how to handling exception case using "try and catch" ?
In my system, I solve 95 ODEs simultaneously. Among them a parameter affects several ODEs and it cause stiffness problem.
I want to change setting parameter value as zero in case of stiffness problem.

Accepted Answer

Mohammad Sami
Mohammad Sami on 3 Sep 2020
Edited: Mohammad Sami on 3 Sep 2020
You can perhaps try a for loop over the various parameters.
% assuming you have 5 parameters needed for your function
% and you have 20 different combinations you want to try
parameters = cell(20,5);
% assign parameters
errors = cell(20,1);
for i = 1:size(parameters,1)
% or short form
catch ME
errors{i} = ME;
Mohammad Sami
Mohammad Sami on 3 Sep 2020
If the error is thrown by the function you are optimizing, you can put try catch inside the function. In the catch you can assign an output value such as 0 or Inf depending on your use case.

Sign in to comment.

More Answers (0)

Community Treasure Hunt

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

Start Hunting!

Translated by