How can I use Genetic Algorithm toolbox to calibrate a model?

1 view (last 30 days)
Hi, this is my first time I use Optimize Live Editor Task.
I have a function, which calculates some outputs throughot the time,
= solvemodel_v12(m_start,m_stop,initfile,'lake',inputfile,'timeseries', parafile,'lake');
Each output is calculated throughout the time step (if it is a year, we would have 365 values for every output).
I have the observation of an output Tzt, and now I want to calibrate my model to minimize RMSQ, and during the calibration process, I would like parameters in the parafile excelsheet to be changed, which its maximum and minimum values are also indicated in that excelsheet. I don't know how I can use FA in live editor task to calibrate this model.

Answers (1)

Alan Weiss
Alan Weiss on 28 Feb 2022
Edited: Alan Weiss on 28 Feb 2022
It is entirely possible that I have misunderstood your question in some way. Feel freee to ask again with more explanation if my answer is not relevant.
First, it is not at all clear to me that ga is the best tool for this task. Why not use lsqnonlin, which is generally faster and more accurate?
But if you have a good reason for using ga then you can write the objective as something like
function Tzt = ...
myfun(m_start,m_stop,initfile,'lake',inputfile,'timeseries', parafile,'lake');
[~,~,~,~,~,~,Tzt,~,~,~,~,~,~,~,~,~,~,~,~,~,~,~,~] = ...
solvemodel_v12(m_start,m_stop,initfile,'lake',inputfile,'timeseries', parafile,'lake');
Call ga like this:
fun = @(parafile)sum((myfun(m_start,m_stop,initfile,'lake',inputfile,'timeseries', parafile,'lake') ...
- (theExpectedValue)).^2);
[sol,fval,efag,output] = ga(fun,nvar,A,b,Aeq,beq,lb,ub,nlcon)
Put in any necessary constraints in the ga call as well.
The solution sol has the new values for the parafile variables. You can write this out to the Excel file after ga runs.
Good luck,
Alan Weiss
MATLAB mathematical toolbox documentation
Alan Weiss
Alan Weiss on 3 Mar 2022
Yes, you have to get the data into MATLAB.
Alan Weiss
MATLAB mathematical toolbox documentation

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!

Translated by