Give gradient to fmincon
2 次查看(过去 30 天)
显示 更早的评论
Hey,
so I am using fmincn to solve a problem. Im somehow not to satisfioed with the results because the solver seems to bug itsself into the lower boundary limit and ste sthe solution for 1 of the 3 parameters to zero therefore. Lokking at the gradient output it seems like there is a potential to optimize, because the gradient for 2 of the variables is like 10^-2, but for the one peramater where the solution is on the boundary, the gradient is like 10^4.
So I know that there is the possibility to set DerivativeCheck on and to give it a user suplied gradient. But I can only activate SpecifyObjectiveGradient to on or off. I am not sure wethere the user suplied gradient should represent a treshold that the output gradient should satisfy or if its an analytical gradient funtion that is required. Unfortunately I can not provide an analytical solution because the function is just to complex, filling multiple pages..
So is there another possibility to stop fmincon from bugging one of the parameters into the boundary since fmincon clearly knows that the gradient at the solution point is pretty bad?
Best regards
6 个评论
Torsten
2022-10-11
And what about the value of the objective ? Which one is "better" - yours or the one from the commercial software ?
But in any case: You should check whether the problem formulations in the two softwares are identical.
采纳的回答
Matt J
2022-10-11
编辑:Matt J
2022-10-11
If the problem is constrained, the gradient in a local extremum usually is not 0.
In fact, if the gradient is non-negative on the lower boundary, it is a very strong sign that a valid local minimum has been found. If you think this is a sub-optimal local minimum, you should try other initial guesses x0.
If the function is too complex to know if it is even differentiable, you may even want to try a non-derivative-based solver like ga, if you have the Global Optimization Toolbox, or fminsearch if you have a small number (<6) of unknnown parameters.
2 个评论
Matt J
2022-10-12
its A=[0,1,1], b=1; but x(2)+x(3) in the solution are always larger 1
That in itself doesn't imply a violation. It depends on the ConstraintTolerance parameters that you've set.
You also need to check exitflags. Possibly, ga could not find a feasible solution because there is none.
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Introduction to Installation and Licensing 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!