General nonlinear programming

Updated:


YALMIP really does not care if your general nonlinear program is convex or not. In the general case it will just be a nonlinear model and any nonlinear solver is used. If it happens to be convex, the solver might perform better, but that is not something YALMIP can influence.

However, general nonlinear models are not all the same in practice. There are nice convex models, and nasty nonconvex models, and everything in between from YALMIPs perspective.

Considerations in nonlinear models

A sure bet to see numerical issues in nonlinear solvers is model with singularities.

Analytic center of polytope - bad vs good forms

As an example, we will find the analytic center of a polytope \(Ax \leq b\). The analytic center is defined as the point which maximizes the expression \( \sum \log(b-Ax)\)

Define a polytope and the decision variable.

n = 5;m = 20;
A = randn(m,n);
b = rand(m,1)*m;
x = sdpvar(n,1);

Solve the problem using the overloaded concave log.

optimize(A*x <= b,-sum(log(b-A*x)))

If you have an exponential cone programming solver installed, this will solve nicely. However, if YALMIP has to revert to a standard nonlinear solver, it can easily fail. The reason is that these solvers often have problems with models where the objective function is undefined for infeasible points, or more generally have singularities.

If you have problems related to this on your model, we can use the exponential operator instead and solve an inverse formulation of the same problem.

y = sdpvar(m,1);
optimize(exp(y) <= b-A*x,-sum(y));

A completely equivalent problem with much better properties in a general nonlinear solver. Note that also this model will be solved with an exponential cone programming solver if available.

Hidden nonlinear equalities

Afraid of nnoconvex models?

By default, YALMIP allows nonconvex problems to be formulated. However, if you want to make sure that no nonconvex problems slips by, you can specify this. If we switch the sign on the objective a non-convex model is obtained, and YALMIP detects this.

optimize(A*x <= b,sum(log(b-A*x)),sdpsettings('allownonconvex',0))