Strong Stationarity Conditions for Optimal Control of Hybrid Systems

We present necessary and sufficient optimality conditions for finite time optimal control problems for a class of hybrid systems described by linear complementarity models. Although these optimal control problems are difficult in general due to the presence of complementarity constraints, we provide a set of structural assumptions ensuring that the tangent cone of the constraints possesses geometric regularity properties. These imply that the classical Karush-Kuhn-Tucker conditions of nonlinear programming theory are both necessary and sufficient for local optimality, which is not the case for general mathematical programs with complementarity constraints. We also present sufficient conditions for global optimality. We proceed to show that the dynamics of every continuous piecewise affine system can be written as the optimizer of a mathematical program which results in a linear complementarity model satisfying our structural assumptions. Hence, our stationarity results apply to a large class of hybrid systems with piecewise affine dynamics. We present simulation results showing the substantial benefits possible from using a nonlinear programming approach to the optimal control problem with complementarity constraints instead of a more traditional mixed-integer formulation.Comment: 30 pages, 4 figure

S-Lemma with Equality and Its Applications

Let $f(x)=x^TAx+2a^Tx+c$ and $h(x)=x^TBx+2b^Tx+d$ be two quadratic functions having symmetric matrices $A$ and $B$. The S-lemma with equality asks when the unsolvability of the system $f(x)<0, h(x)=0$ implies the existence of a real number $\mu$ such that $f(x) + \mu h(x)\ge0, ~\forall x\in \mathbb{R}^n$. The problem is much harder than the inequality version which asserts that, under Slater condition, $f(x)<0, h(x)\le0$ is unsolvable if and only if $f(x) + \mu h(x)\ge0, ~\forall x\in \mathbb{R}^n$ for some $\mu\ge0$. In this paper, we show that the S-lemma with equality does not hold only when the matrix $A$ has exactly one negative eigenvalue and $h(x)$ is a non-constant linear function ($B=0, b\not=0$). As an application, we can globally solve $\inf\{f(x)\vert h(x)=0\}$ as well as the two-sided generalized trust region subproblem $\inf\{f(x)\vert l\le h(x)\le u\}$ without any condition. Moreover, the convexity of the joint numerical range $\{(f(x), h_1(x),\ldots, h_p(x)):~x\in\Bbb R^n\}$ where $f$ is a (possibly non-convex) quadratic function and $h_1(x),\ldots,h_p(x)$ are affine functions can be characterized using the newly developed S-lemma with equality.Comment: 34 page

Convex inner approximations of nonconvex semialgebraic sets applied to fixed-order controller design

We describe an elementary algorithm to build convex inner approximations of nonconvex sets. Both input and output sets are basic semialgebraic sets given as lists of defining multivariate polynomials. Even though no optimality guarantees can be given (e.g. in terms of volume maximization for bounded sets), the algorithm is designed to preserve convex boundaries as much as possible, while removing regions with concave boundaries. In particular, the algorithm leaves invariant a given convex set. The algorithm is based on Gloptipoly 3, a public-domain Matlab package solving nonconvex polynomial optimization problems with the help of convex semidefinite programming (optimization over linear matrix inequalities, or LMIs). We illustrate how the algorithm can be used to design fixed-order controllers for linear systems, following a polynomial approach

An elementary proof of the Fritz-John and Karush-Kuhn-Tucker conditions in nonlinear programming

In this note we give an elementary proof of the Fritz-John and Karush-Kuhn-Tucker conditions for nonlinear finite dimensional programming problems with equality and/or inequality constraints.The proof avoids the implicit function theorem usually applied when dealing with equality constraints and uses a generalization of Farkas lemma and the Bolzano-Weierstrass property for compact sets

Necessary and Sufficient Optimality Conditions for Vector Equilibrium Problems on Hadamard Manifolds

The aim of this paper is to show the existence and attainability of Karush–Kuhn–Tucker optimality conditions for weakly efficient Pareto points for vector equilibrium problems with the addition of constraints in the novel context of Hadamard manifolds, as opposed to the classical examples of Banach, normed or Hausdorff spaces. More specifically, classical necessary and sufficient conditions for weakly efficient Pareto points to the constrained vector optimization problem are presented. The results described in this article generalize results obtained by Gong (2008) andWei and Gong (2010) and Feng and Qiu (2014) from Hausdorff topological vector spaces, real normed spaces, and real Banach spaces to Hadamard manifolds, respectively. This is done using a notion of Riemannian symmetric spaces of a noncompact type as special Hadarmard manifolds