When does fsolve use Levenberg-Marquardt algorithm?
8 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Manuela Gräfe
el 16 de Mayo de 2017
Comentada: John D'Errico
el 16 de Mayo de 2017
Hi,
when does Matlabs fsolve method uses the Levenberg-Marquardt algorithm? The warning tells me, that the trust-region-dogleg algorithm cannot handle non-square systems.
What are the features/characteristics of a non-square system? Do I have more unknowns than equations? What else can be the 'problem'?
Thanks!
0 comentarios
Respuesta aceptada
John D'Errico
el 16 de Mayo de 2017
With fewer unknowns than equations, nothing you do can "solve" the problem, because there will be normally be infinitely many solutions (if any exist). At best, you can find a completely arbitrary solution among many. The solution you obtain will be completely dependent on your start point. Change the start point, and you will find some completely different solution.
With more unknowns than equations, you will have no EXACT solution in general. At best, you will need to find a set of parameters that minimizes the sum of squares of residuals. This is what algorithms like Levenberg-Marquardt are designed to solve, and why that algorithm must be invoked. At that point, you might as well have decided to used lsqnonlin instead of fsolve, since lsqnonlin is designed to solve that class of problem anyway.
As for what else can be the problem, it tells you exactly what the problem was! You have a non-square system.
2 comentarios
John D'Errico
el 16 de Mayo de 2017
If you change the number of fixed parameters, fsolve must change the way it solves the problem. Not all such problems will have an exact solution.
Más respuestas (0)
Ver también
Categorías
Más información sobre Nonlinear Dynamics 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!