Can one run fminsearch for a function defined by cases?
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
alpedhuez
el 28 de Jul. de 2020
Comentada: Matt J
el 28 de Jul. de 2020
First, I have a function defined by cases like
if x<1 f(x)=x
if x>1 f(x)=sqrt(x)
Second, suppose I want to calculate the maximium of f for 0<x<2 (this problem itself is trivial, but as an example).
Third, can one run fminsearch for such problem?
0 comentarios
Respuesta aceptada
Más respuestas (1)
John D'Errico
el 28 de Jul. de 2020
Editada: John D'Errico
el 28 de Jul. de 2020
Can you? Yes.
Should you? This may sometimes be a poor idea, although it depends on the specific segment functions. Are the segments continuous across the break? Are they differentiable across the break?
Fminsearch does not require a several times continuously differentiable objective. It is even reasonably robust in that respect. But if your function is discontinuous across a break, then expect that it can easily fail. The nastier is your function, the more chance of failure.
What does failure mean here? It does not mean the algorithm will fail due to an error message. It just means fminsearch may not consistently find the optimal solution. It could get more easily stuck in a sub-optimal solution. Note there is NO assurance of convergence to the globally optimal solution for fminsearch (or virtually any optimizer) even on a more well-behaved problem, just that if you make things messy it may fail more easily.
If your problem is a 1-dimensional one, then it is advisable to best break down the problem into sub-intervals, optimizing over each interval, and to then use fminbnd for each interval. This will yield the best possile result. Using fminbnd is good here, because it employs a bracket to search over. You then choose the overall best choice of all solutions returned. In the case of a bracket that extends to -inf or +inf, just choose a reasonably large number, but not excessively large, as that can cause numerical problems.
Ver también
Categorías
Más información sobre Optimization en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!