Optimization with affine homogeneous quadratic integral inequality constraints

G. Fantuzzi, A. Wynn, P. J. Goulart and A. Papachristodoulou

BibTeX  Preprint  Code 

  author = {G. Fantuzzi and A. Wynn and P. J. Goulart and A. Papachristodoulou},
  title = {Optimization with affine homogeneous quadratic integral inequality constraints},
  year = {2016}

We introduce a new technique to optimize a linear cost function subject to an affine homogeneous quadratic integral inequality, i.e., the requirement that a given homogeneous quadratic integral functional, affine in the optimization variables, is non-negative over a space of functions defined by homogeneous boundary conditions. Problems of this type often arise when studying dynamical systems governed by partial differential equations. First, we derive a hierarchy of outer approximations for the feasible set of a homogeneous quadratic integral inequality in terms of linear matrix inequalities (LMIs), and show that a convergent, non-decreasing sequence of lower bounds for the optimal cost value can be computed by solving a sequence of semidefinite programs (SDPs). Second, we obtain inner approximations in terms of LMIs and sum-of-squares constraints, which enable us to formulate SDPs to compute upper bounds for the optimal cost, as well as to compute a strictly feasible point for the integral inequality. To aid the formulation and solution of our SDP relaxations, we implement our techniques in QUINOPT, an open-source add-on to the MATLAB optimization toolbox YALMIP. We demonstrate our techniques by solving typical problems that arise in the context of stability analysis for dynamical systems governed by PDEs.