non linear least squares problem using BFGS

1 次查看(过去 30 天)
Hello,
I need to find solution to a non linear least squares problem using Gauss-Newton method, however I am only able to compute Jacobian matrix for a simple model, so I need to use updating methods like BFGS instead full Jacobian computation. Does anybody know how to solve this problem using Matlab?
Thanks
  2 个评论
Matt J
Matt J 2013-1-21
编辑:Matt J 2013-1-21
That sounds like it could have doubtful convergence properties. Gauss-Newton is already not guaranteed to converge, even in its pure form with the exact Jacobian. Seems like using an approximate Jacobian might only make things worse. If Benji's answer doesn't work for you (because it's not Gauss-Newton), maybe you could elaborate on why Gauss-Newton is your only option.
David Kusnirak
David Kusnirak 2013-1-22
Thanks for hint. As I'm trying to follow 'recipe' from an arcticle, I'd prefer the Gauss-Newton method. I've tried to write my own script, but it yields quite poor results, so I would like to compare it somehow with Matlab function to improve the results.

请先登录,再进行评论。

采纳的回答

Shashank Prasanna
Shashank Prasanna 2013-1-21
Do you want to code this yourself or is it OK to use something that comes with MATLAB?
For the FMINUNC function in the Optimization Toolbox, if set the Algorithm to Medium Scale it will use BFGS Quasi-Newton method with a cubic line search procedure.
Take a look at the following documentation page and scroll down to the section call Algorithm > Medium-Scale Optimization
  2 个评论
David Kusnirak
David Kusnirak 2013-1-22
Thank you. I'm going to try this way, but till now I didn't find out how to set the algorithm for the medium scale..
Shashank Prasanna
Shashank Prasanna 2013-1-22
Do the following:
>> options = optimset('LargeScale','off');
>> x = fminunc(fun,x0,options)
I urge you to read the FMINUNC documentation:

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Solver Outputs and Iterative Display 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by