How to handle stiffness parameter using Try and Catch?

Question

Jinsu Kim 2020-9-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/588277-how-to-handle-stiffness-parameter-using-try-and-catch

评论： Mohammad Sami 2020-9-3

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.

Regards,

Kim.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Mohammad Sami 2020-9-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/588277-how-to-handle-stiffness-parameter-using-try-and-catch#answer_488944

编辑：Mohammad Sami 2020-9-3

在 MATLAB Online 中打开

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)
    try
        somefunction(parameters{i,1},parameters{i,2},parameters{i,3},parameters{i,4},parameters{i,5});
        % or short form
        %somefunction(parameters{i,:});
    catch ME
        errors{i} = ME;
    end
end

2 个评论
显示无隐藏无

Jinsu Kim 2020-9-3

Thanks for your opinion. Actually it's good way to try.

I have one more question. I want to handle exception case more automatically. How can I apply your idea when using GA (Genetic Algorithm) ?

This is the full code of GA.m

fitness = @Main ;

parameter = 5 ;

LB = [ 1*10^(-7), 40, 1*10^(-7), 40, 1*10^(-7) ] ;

UB = [ 1.0*10^(-5), 100, 1.0*10^(-5), 100, 1.0*10^(-5) ] ;

[x, fval] = ga (fitness, parameter, [], [], [], [], LB, UB ) ;

Mohammad Sami 2020-9-3

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.

请先登录，再进行评论。