Computation of Lyapunov Functions under State Constraints using Semidefinite Programming Hierarchies

We provide algorithms for computing a Lyapunov function for a class of systems where the state trajectories are constrained to evolve within a closed convex set. The dynamical systems that we consider comprise a differential equation which ensures continuous evolution within the domain, and a normal cone inclusion which ensures that the state trajectory remains within a prespecified set at all times. Finding a Lyapunov function for such a system boils down to finding a function which satisfies certain inequalities on the admissible set of state constraints. It is well-known that this problem, despite being convex, is computationally difficult. For conic constraints, we provide a discretization algorithm based on simplicial partitioning of a sim-plex, so that the search of desired function is addressed by constructing a hierarchy (associated with the diameter of the cells in the partition) of linear programs. Our second algorithm is tailored to semi-algebraic sets, where a hierarchy of semidefinite programs is constructed to compute Lyapunov functions as a sum-of-squares polynomial.

[1]  Didier Henrion,et al.  Optimization on linear matrix inequalities for polynomial systems control , 2013, 1309.3112.

[2]  Vincent Acary,et al.  Nonsmooth Modeling and Simulation for Switched Circuits , 2010, Lecture Notes in Electrical Engineering.

[3]  G. Chesi.  Homogeneous Polynomial Forms for Robustness Analysis of Uncertain Systems , 2009 .

[4]  Mirjam Dür,et al.  Algorithmic copositivity detection by simplicial partition , 2008 .

[5]  Antonis Papachristodoulou,et al.  Robust Stability Analysis of Nonlinear Hybrid Systems , 2009, IEEE Transactions on Automatic Control.

[6]  Amir Ali Ahmadi,et al.  SOS-Convex Lyapunov Functions and Stability of Difference Inclusions , 2018, ArXiv.

[7]  Peter J Seiler,et al.  SOSTOOLS: Sum of squares optimization toolbox for MATLAB , 2002 .

[8]  Amir Ali Ahmadi,et al.  Sum of Squares Certificates for Stability of Planar, Homogeneous, and Switched Systems , 2017, IEEE Transactions on Automatic Control.

[9]  M. Kanat Camlibel,et al.  Lyapunov Stability of Complementarity and Extended Systems , 2006, SIAM J. Optim..

[10]  S. Adly.  A Variational Approach to Nonsmooth Dynamics , 2017 .

[11]  Bernard Brogliato,et al.  Dynamical Systems Coupled with Monotone Set-Valued Operators: Formalisms, Applications, Well-Posedness, and Stability , 2020, SIAM Rev..

[12]  Bernard Brogliato,et al.  Well-Posedness and Output Regulation for Implicit Time-Varying Evolution Variational Inequalities , 2016, SIAM J. Control. Optim..

[13]  V. Powers,et al.  An algorithm for sums of squares of real polynomials , 1998 .

[14]  P. Parrilo.  Structured semidefinite programs and semialgebraic geometry methods in robustness and optimization , 2000 .

[15]  J. Lasserre.  An Introduction to Polynomial and Semi-Algebraic Optimization , 2015 .

[16]  Xinzhen Zhang,et al.  A Complete Semidefinite Algorithm for Detecting Copositive Matrices and Tensors , 2017, SIAM J. Optim..

[17]  Mirjam Dür,et al.  An Adaptive Linear Approximation Algorithm for Copositive Programs , 2009, SIAM J. Optim..

[18]  N. Wouw,et al.  Stability and Convergence of Mechanical Systems with Unilateral Constraints , 2008 .

[19]  A. Garulli,et al.  Positive Polynomials in Control , 2005 .

[20]  Daniel Goeleven,et al.  Stability and instability matrices for linear evolution variational inequalities , 2004, IEEE Transactions on Automatic Control.