levenberg-marquardt with linear constraints

7 次查看（过去 30 天）
Eze Far
Eze Far 2023-8-8
评论： Eze Far 2023-8-8
Hi, I just see the the box constraints are now allowed for the levenberg-marquardt algorithm in lsqnonlin. However it seems that linear constraints are still not implemented. Is that correct? Are there any plans to implement that in future realeases?
Cheers
  6 个评论
显示 4更早的评论
Bruno Luong
Bruno Luong 2023-8-8
编辑：Bruno Luong 2023-8-8
AFAIK (As Far As I Know) it is not trivial to introduce constraints in LM algorithm. There are few papers on the topic but probably not mature enough to be implemented at the presence time.

请先登录，再进行评论。

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

回答（2 个）

Torsten
Torsten 2023-8-8
编辑：Torsten 2023-8-8
x = lsqnonlin(fun,x0,lb,ub,A,b,Aeq,beq,nonlcon)
A,b,Aeq,beq define the linear constraints.
  1 个评论
Eze Far
Eze Far 2023-8-8
Yes, that's what I thought, but in https://nl.mathworks.com/help/optim/ug/lsqnonlin.html one can read:
"The 'interior-point' algorithm is the only algorithm that can solve problems with linear or nonlinear constraints. If you include these constraints in your problem and do not specify an algorithm, the solver automatically switches to the 'interior-point' algorithm."
I also tried myself and it switches automatically to interior-point when the linear constraints are enabled. That doesn't happen when only box constraints are set.

请先登录，再进行评论。

John D'Errico
John D'Errico 2023-8-8
编辑：John D'Errico 2023-8-8
I'm sorry, but nobody can tell you IF there are any plans to introduce some behavior in a future release, until it is released.
If someone does know, then they will be under a non-disclosure agreement. If they don't know, well what do you expect them to be able to tell you?
At best, you can submit a feature request to MathWorks. DIRECTLY TO THEM, NOT THROUGH ANSWERS. Answers is not technical support. Ask for this feature to be included in a future release. Remember that TMW will generally not write code that nobody ever wants or hopes to use. Instead, they will direct their efforts at things for which people have shown a need.
https://servicerequest.mathworks.com/mysr/cp_case_new1?cc=us
The nice thing is, I have always found them to be responsive to reasonable requests. In fact, if they do add the idea as a feature request, they will also attach your e-mail to it, so when it does come alive, you will be notified directly.
  3 个评论
显示 1更早的评论
John D'Errico
John D'Errico 2023-8-8
编辑：John D'Errico 2023-8-8
@Eze Far
Geez, RELAX!
  1. TMW employees will be under an NDA. How can you not expect that?
  2. I put it in caps, since far too many people seem to think Answers is indeed official tech support, so I wanted to make that part clear, drawing explicit attention to it. I could have made it bold face or italicized it instead. Would that seriously have made you happier? Really? Did my use of caps truly upset you? relax. relax. Do those make you happier?
  3. I did not say that nobody will need that option, only that they may not have decided to expend the effort to add that option to the LM codes. Mybe they did not think it was necessary. Remember that every hour spent on one option is an hour not invested in development of other tools. And since they only have so many employees, this is a tradeoff they must always make. If they already offered a capability under one sub-optimizer in lsqnonlin, I can easily understand they might have chosen not to put it into all sub-optimzers in lsqnonlin. Honestly, it has been too many years since I worked seriously with Levenberg-Marquardt and writing optimization codes, that I don't know if there may be some explicit reason why they did not put it into the LM code. That they did not is a suggestion there was some reason not to do so, if it were easy.
Please develop a thicker skin, and don't get all upset when someone gives you some information that you clearly needed. As for your advice, I don't give a tinker's damn. In fact, I felt I was being quite reasonably polite. My feelings towards you are not so much anymore, but I doubt you care anymore than I care about yours. The fact is, you obviously learned several things from my response.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Solver Outputs and Iterative Display 的更多信息

标签

产品


版本

R2023a

Community Treasure Hunt

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

Start Hunting!

Translated by