Geometric programming
This
example requires MOSEK,
GPPOSY or
fmincon Geometric programming solvers are capable of solving a sub-class of geometric problems where c≥0 with the additional constraint t≥0, so called posynomial geometric programming. The following example is taken from the MOSEK manual. (note, the positivity constraint on t will be added automatically)
If the geometric program violates the posynomial assumption, an error will be issued.
YALMIP will automatically convert some simple violations of the posynomial assumptions, such as lower bounds on monomial terms and maximization of negative monomials. The following small program maximizes the volume of an open box, under constraints on the floor and wall area, and constraints on the relation between the height, width and depth (example from [S. Boyd, S. Kim, L. Vandenberghe, A. Hassibi] ).
The posynomial geometric programming problem is not convex in its
standard formulation. Hence, if a general nonlinear solver is applied to
the problem, it will typically fail. However, by performing a suitable
logarithmic variable transformation, the problem is rendered convex.
YALMIP has built-in support for performing this variable change, and solve
the problem using the nonlinear solver
fmincon. To invoke this module in
YALMIP, use the solver
tag
The current
version of YALMIP has a bug that may cause problems if you have convex
quadratic constraints. To avoid this problem, use Never use the commands sqrt and cpower when working with geometric programs, i.e. always use the ^ operator. The reason is implementation issues in YALMIP. The commands sqrt and cpower are meant to be used in optimization problems where a conic model is derived using convexity propagation, see nonlinear operators. Generalized geometric programmingSome geometric programs, although not given in standard form, can still be solved using a standard geometric programming solver after some some additional variables and constraints have been introduced. YALMIP has built-in support for some of these conversion. To begin with, nonlinear operators can be used also in geometric programs, as in any other optimization problems (as long as YALMIP is capable of proving convexity, see the nonlinear operator examples)
Powers of posynomials are allowed in generalized geometric programs. YALMIP will automatically take care of this and convert the problems to a standard geometric programs. Note that the power has to be positive if used on the left-hand side of a <, and negative otherwise.
To understand how a generalized geometric program can be converted to a standard geometric program. the reader is referred to [S. Boyd, S. Kim, L. Vandenberghe, A. Hassibi] Mixed integer geometric programmingThe branch and bound solver in YALMIP is built in a modular fashion that makes it possible to solve almost arbitrary convex mixed integer programs. The following example is taken from [S. Boyd, S. Kim, L. Vandenberghe, A. Hassibi]. To begin with, define the data for the example.
Introduce symbolic expressions used in the model.
The objective function and constraints (notice the use of the nonlinear operator max in the objective).
Solve!
An alternative model is discussed in the paper, and is just as easy to define.
|