How do I help quadprog converge?

5 次查看（过去 30 天）

Leigh Sneddon 2019-8-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/476207-how-do-i-help-quadprog-converge

编辑： Matt J 2019-8-16

Quadprog quits with exitFlag = -2, and the message below. How do I help it converge to a feasible point?

"quadprog optimization failed: Converged to an infeasible point.

quadprog stopped because the size of the current step is less than the default value of the step size tolerance but constraints are not satisfied to within the selected value of the constraint tolerance.

Stopping criteria details:

Optimization stopped because the relative changes in all elements of x are less than options.StepTolerance = 1.000000e-12, but the relative maximum constraint violation, 1.902013e-14, exceeds options.ConstraintTolerance = 1.000000e-06.

Optimization Metric Options

max(abs(delta_x./x)) = 1.82e-13 StepTolerance = 1e-12 (default)

relative max(constraint violation) = 1.90e-14 ConstraintTolerance = 1e-06 (selected)"

9 个评论
显示 7更早的评论隐藏 7更早的评论

Walter Roberson 2019-8-15

https://www.mathworks.com/help/optim/ug/when-the-solver-fails.html

I cannot find the information on the problem that came to mind; unfortunately the bug reports are now difficult to search :(

Leigh Sneddon 2019-8-15

OK. Thanks.

请先登录，再进行评论。

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

回答（1 个）

Matt J 2019-8-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/476207-how-do-i-help-quadprog-converge#answer_387757

编辑：Matt J 2019-8-15

在 MATLAB Online 中打开

Is checking that the constraints are satisfied and the optimality measure is low a good rule of thumb for making this decision?

Checking the first order KKT conditions would be the best test, assuming your quadratic is convex. The final output argument of quadprog gives the solver's idea of the optimal Lagrange multipliers,

[x,fval,exitflag,output,lambda] = quadprog(___)

But I would first recommend upgrading to a Matlab version that doesn't have this bug.

2 个评论
显示无隐藏无

Leigh Sneddon 2019-8-16

Optimizer precision limitations will mean that none of the conditions is satisfied exactly. Is there a way of knowing, given the tolerances used, how much mismatch is consistent with a correct solution?

Matt J 2019-8-16

编辑：Matt J 2019-8-16

What does "consistent with a correct solution" mean to you? Even if quadprog's exit message had been a proper one, what is the deviation distance from the true optimum that your application can tolerate, and how would you have known that the result is within that distance if quadprog had behaved normally?

请先登录，再进行评论。

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