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

3 次查看(过去 30 天)
Hi, this is my first time I use Optimize Live Editor Task.
I have a function, which calculates some outputs throughot the time,
[zz,Az,Vz,tt,Qst,Kzt,Tzt,Czt,Szt,Pzt,Chlzt,PPzt,DOPzt,DOCzt,Qzt_sed,lambdazt,...
P3zt_sed,P3zt_sed_sc,His,DoF,DoM,MixStat,Wt]...
= 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.

回答(1 个)

Alan Weiss
Alan Weiss 2022-2-28
编辑:Alan Weiss 2022-2-28
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');
end
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
  6 个评论
Sahar khalili
Sahar khalili 2022-3-1
Thank you so much. In the function, paraval I think is a MATLAB variable, while in order to run the function solvemodel_v12 I need to put this as an excel file, in the 'lake' sheet. So, I have to use xlsread and xlswrite?
Alan Weiss
Alan Weiss 2022-3-3
Yes, you have to get the data into MATLAB.
Alan Weiss
MATLAB mathematical toolbox documentation

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Genetic Algorithm 的更多信息

产品


版本

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by