フィルターのクリア

Can I use fmincon with a whole model as function?

11 ビュー (過去 30 日間)
Timo Beune
Timo Beune 2024 年 6 月 21 日
コメント済み: Timo Beune 2024 年 6 月 24 日 11:13
I am currently working on an optimization for a flow and electric current problem. The parameters to be optimized are geometric parameters like lengths, diameters, thicknesses of channels etc. I want to use fmincon, but to my understanding I would need to define a function that is dependent on my variables. However the problem itself needs to iterate before a solution is found. Therefore, fmincon would need to put the values of x into the flow and current model, calculate values for current efficiency and pressure drop, and then use these values to optimize. Is that in someway possible? and if yes, how would I define fmincon and the function itself? The model itself is too big to put in, but i will try to give an example:
The value of Dp depends on the input parameters defined in x and the same goes for Ic,k but cannot immediatly be linked to them since they first need to iterate. Fmincon then needs to find an optimum for the system efficiency, which is defined as . So basically my function would be something like , where x is a vector containing values for channel length, diameter, depth etc. Hopefully this is somewhat clear.
  2 件のコメント
Sam Chak
Sam Chak 2024 年 6 月 21 日
Hi @Timo Beune, the MATLAB code for writing the objection function has been provided by @Torsten.
Because fmincon() only minimizes the objection function, if you wish to maximize the system efficiency, then the objection function should be written as
,
unless I misunderstand the concept of optimality in system efficiency.
Torsten
Torsten 2024 年 6 月 21 日
Thanks for the minus - I added it below.

サインインしてコメントする。

採用された回答

Torsten
Torsten 2024 年 6 月 21 日
編集済み: Torsten 2024 年 6 月 21 日
I don't see your problem.
You call "fmincon" as
sol = fmincon(@fun,x0)
and write a function "fun" as
function obj = fun(x)
% Compute P_p from x (maybe by iteration or even software different from MATLAB)
% Compute eta_i from x (maybe by iteration or even software different from MATLAB)
obj = -eta_i/P_p
end
  1 件のコメント
Timo Beune
Timo Beune 2024 年 6 月 24 日 11:13
Thanks. It was namely a matter of code construction, not so much a problem. Thanks for your help

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeSurrogate Optimization についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by