20.05.2014 Views

link to my thesis

link to my thesis

link to my thesis

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Chapter 4<br />

Global optimization of multivariate<br />

polynomials<br />

Finding the global minimum of a real-valued multivariate polynomial is a problem<br />

which has several useful applications in systems and control theory as well as in many<br />

other quantitative sciences including statistics, mathematical finance, economics, systems<br />

biology, etc. Multivariate global polynomial optimization is a hard problem because<br />

of the non-convexity of the problem and the existence of local optima. In this<br />

chapter we focus on computing the global optimum of a special class of polynomials.<br />

The class of polynomials discussed in this chapter is the class of so-called Minkowski<br />

dominated polynomials in several variables. When working with a polynomial p λ (x 1 ,<br />

...,x n ) from this class, the problem of finding its global minimum can be reformulated<br />

as an eigenvalue problem by applying the Stetter-Möller matrix method introduced in<br />

the previous chapter. This is possible because the system of first-order conditions of<br />

such a polynomial is always in Gröbner basis form with a special structure, showing<br />

that the number of solutions is finite. Applying the Stetter-Möller matrix method<br />

yields a set of real non-symmetric large and sparse commuting matrices A T x 1<br />

,...,A T x n<br />

.<br />

Using the same approach, it is also possible <strong>to</strong> construct a matrix A T p λ (x which<br />

1,...,x n)<br />

has some useful properties: it turns out that the smallest real eigenvalue of the matrix<br />

A T p λ x 1,...,x n)<br />

is of interest for computing the global optimum of the polynomial<br />

p λ (x 1 ,...,x n ).<br />

Various extensions and generalizations of these techniques exist including the optimization<br />

of rational or logarithmic functions [63] and the use of systems that are<br />

not in Gröbner basis form [12]. However, these generalizations do not satisfy all the<br />

properties described below and therefore require adjustments <strong>to</strong> the techniques used<br />

here.<br />

Numerical experiments using the approach described in this chapter are discussed<br />

in Chapter 7.<br />

43

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!