minimization fmincon with ode

3 次查看（过去 30 天）
AM
AM 2019-1-8
评论： Torsten 2019-1-8
Hello, I am using fmincon to fit my model to data.
I have a set of differential equations that give me the model results:
dy/dt=f(x,y,t)
where y are the points I want to compare to experimental points, x are my decision variables and t is time.
My objective function objfun is the sum of the squares of the residuals (y-data)^2
I profiled my code and I saw that it spends the most time solving the ode. I wanted to know if there is a way, numerically, to use my set of ODE dy/dt to determine the gradient dobjfun/dx so I can give it to fmincon beforehand instead of it using finite differences to determine it.
Thank you,

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

采纳的回答

Torsten
Torsten 2019-1-8
编辑：Torsten 2019-1-8
To get dobjectfun/dx numerically, you had to solve even more ODEs:
https://de.mathworks.com/help/optim/ug/optimizing-a-simulation-or-ordinary-differential-equation.html
Section: Use a Gradient Evaluation Function.
I wouldn't advice you to do so.
Best wishes
Torsten.
  2 个评论
AM
AM 2019-1-8
I see, thank you!
Does that mean that if all the other aspects of my code were optimized (which isn't the case yet), there is nothing else that can be done regarding this bottleneck to improve overall efficiency? I can't provide an anlytical expression of the gradient so I wanted to see if it could be done numerically.
Torsten
Torsten 2019-1-8
You can use derivative-free optimization methods, e.g. fminsearch.
For methods that need the Jacobian, nothing can be done.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

产品


版本

R2017b

Community Treasure Hunt

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

Start Hunting!

Translated by