can I use a compiled (.exe) file to generate my objective functions while running optimisation in MATLAB?

Question

Mohsen2015 2015-6-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/224110-can-i-use-a-compiled-exe-file-to-generate-my-objective-functions-while-running-optimisation-in-ma

评论： Kurt Stewart 2019-9-18

I have developed a Fortran code which simulates a thermodynamic cycle. I want to optimise the cycle for few objectives (like cost, efficiency etc.). My Fortran program reads design variables [X] from a file and writes the objectives values F[X] in a separate file. I can call the .exe file of my program in MATLAB and it communicates with the input and output files without problem. I want to see whether it is possible to call this .exe file to generate objective functions from [X] in each iteration of optimisation procedure? If yes, which types of optimisation solvers I may choose?

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

请先登录，再进行评论。

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

Answer 1

Andrew Schenk 2015-6-17

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/224110-can-i-use-a-compiled-exe-file-to-generate-my-objective-functions-while-running-optimisation-in-ma#answer_183041

在 MATLAB Online 中打开

Optimization functions like fminunc simply need a MATLAB function of the form y = f(x) to operate on. To interface with your external code, you should create a MATLAB function similar to the one below:

 function y = myfun(x)
%write the x data
xh = fopen('x.txt', 'w');
fprintf(xh, '%f', x);
fclose(xh);
 %run your program
system('yourProgram.exe');
 %read the f data
fh = fopen('f.txt');
y = fscanf(fh, '%f');
fclose(fh);

Then use fminunc:

minunc(@myfun,x0);

6 个评论
显示 4更早的评论隐藏 4更早的评论

Mohsen2015 2015-6-18

在 MATLAB Online 中打开

Thanks Andrew, but when I use the function as you suggested, it returns an error:

At line 9 of file ../test2.f (unit = 11, file = 'x.txt')
Fortran runtime error: Bad real number in item 1 of list input
Error using fminunc (line 341)
User supplied objective function must return a scalar value.

I tested the function below as a sample to see if I can replace it with a compiled fortran program.

 function f = myfun(x)
f = 3*x(1)^2 + 2*x(1)*x(2) + x(2)^2;    % Cost function

This is the fortran program (test2.exe) to read x1 and x2 from x.txt and write the function result to f.txt:

      program test2
      implicit none
      real y
      real x
      dimension x(2)
      open (unit = 11, file = "x.txt", status= "old")
      open (unit = 21, file = "f.txt", status= "replace")
      read(11,*) x(1)
      read(11,*) x(2)
      y = 3*x(1)**2 + 2*x(1)*x(2) + x(2)**2
      write(21,*) y
      end program test2

and then used it as you suggested:

function y = myfun(x)
%write the x data
xh = fopen('x.txt', 'w');
fprintf(xh, '%f', x);
fclose(xh);
%run your program
system('test2.exe');
%read the f data
fh = fopen('f.txt');
y = fscanf(fh, '%f');
fclose(fh);

Any ideas what should I do?

Mohsen2015 2015-6-22

Thanks

I used fminsearch instead and it worked with no problem either defining the function in MATLAB or calling a compiled Fortran .exe file. I also ran a two objective function optimisation and used "gamultiobj" (the example function given in MATLAB help). No problems and it works properly. Just two questions are still remaining:

1. why I can get correct answer with fminunc regardless of x0 value when define the function directly in MATLAB, but it has problems when calling .exe file? It only gives correct minimum when I choose x0 close enough.

2. Using .exe file increased the calculation time, especially in "gamultiobj". How can I reduce the time?

I'm going to ask this in a separate question.

Kurt Stewart 2019-9-18

fminunc is a local solver - thus it gets stuck in local minima, that is why you would need to have the initial conditions very close to the solution to actually find it in a non-convex problem. gamultiobj uses a genetic algorithm. They just inherently take quite a long time. Try ParetoSearch and see how that works, probably faster than gamultiobj

请先登录，再进行评论。

can I use a compiled (.exe) file to generate my objective functions while running optimisation in MATLAB?

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

采纳的回答

6 个评论
显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

can I use a compiled (.exe) file to generate my objective functions while running optimisation in MATLAB?

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

采纳的回答

6 个评论 显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

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

6 个评论
显示 4更早的评论隐藏 4更早的评论