Cone-Copositive Lyapunov Functions for Complementarity Systems: Converse Result and Polynomial Approximation

This article establishes the existence of Lyapunov functions for analyzing the stability of a class of state-constrained systems, and it describes algorithms for their numerical computation. The system model consists of a differential equation coupled with a set-valued relation that introduces discontinuities in the vector field at the boundaries of the constraint set. In particular, the set-valued relation is described by the subdifferential of the indicator function of a closed convex cone, which results in a cone-complementarity system. The question of analyzing the stability of such systems is addressed by constructing cone-copositive Lyapunov functions. As a first analytical result, we show that exponentially stable complementarity systems always admit a continuously differentiable cone-copositive Lyapunov function. Putting some more structure on the system vector field, such as homogeneity, we can show that the aforementioned functions can be approximated by a rational function of cone-copositive homogeneous polynomials. This latter class of functions is seen to be particularly amenable for numerical computation as we provide two types of algorithms for precisely that purpose. These algorithms consist of a hierarchy of either linear or semidefinite optimization problems for computing the desired cone-copositive Lyapunov function. Some examples are given to illustrate our approach.

[1]  Katta G. Murty,et al.  Some NP-complete problems in quadratic and nonlinear programming , 1987, Math. Program..

[2]  Jong-Shi Pang,et al.  A Posteriori Error Bounds for the Linearly-Constrained Variational Inequality Problem , 1987, Math. Oper. Res..

[3]  O. Mangasarian,et al.  Lipschitz continuity of solutions of linear inequalities, programs and complementarity problems , 1987 .

[4]  K. Schmüdgen.  TheK-moment problem for compact semi-algebraic sets , 1991 .

[5]  L. Rosier.  Homogeneous Lyapunov function for homogeneous continuous vector field , 1992 .

[6]  F. Clarke,et al.  Subgradient Criteria for Monotonicity, The Lipschitz Condition, and Convexity , 1993, Canadian Journal of Mathematics.

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

[8]  Yu. S. Ledyaev,et al.  Asymptotic Stability and Smooth Lyapunov Functions , 1998 .

[9]  A. Teel,et al.  A Smooth Lyapunov Function from a Class-kl Estimate Involving Two Positive Semideenite Functions , 1999 .

[10]  Clyde F. Martin,et al.  A Converse Lyapunov Theorem for a Class of Dynamical Systems which Undergo Switching , 1999, IEEE Transactions on Automatic Control.

[11]  Etienne de Klerk,et al.  On Copositive Programming and Standard Quadratic Optimization Problems , 2000, J. Glob. Optim..

[12]  P. Parrilo,et al.  Semidefinite programming based tests for matrix copositivity , 2000, Proceedings of the 39th IEEE Conference on Decision and Control (Cat. No.00CH37187).

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

[14]  B. Reznick,et al.  A new bound for Pólya's Theorem with applications to polynomials positive on polyhedra , 2001 .

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

[16]  Bernard Brogliato,et al.  Some perspectives on the analysis and control of complementarity systems , 2003, IEEE Trans. Autom. Control..

[17]  Daniel Liberzon,et al.  Switching in Systems and Control , 2003, Systems & Control: Foundations & Applications.

[18]  F. Facchinei,et al.  Finite-Dimensional Variational Inequalities and Complementarity Problems , 2003 .

[19]  Johannes Schumacher,et al.  Copositive Lyapunov functions , 2004 .

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

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

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

[23]  Yacine Chitour,et al.  Common Polynomial Lyapunov Functions for Linear Switched Systems , 2006, SIAM J. Control. Optim..

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

[25]  Matthew M. Peet,et al.  Exponentially Stable Nonlinear Systems Have Polynomial Lyapunov Functions on Bounded Regions , 2007, IEEE Transactions on Automatic Control.

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

[27]  Richard W. Cottle,et al.  Linear Complementarity Problem , 2009, Encyclopedia of Optimization.

[28]  M. Laurent.  Sums of Squares, Moment Matrices and Optimization Over Polynomials , 2009 .

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

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

[31]  Alberto Seeger,et al.  A Variational Approach to Copositive Matrices , 2010, SIAM Rev..

[32]  Mirjam Dür,et al.  Copositive Programming – a Survey , 2010 .

[33]  H. Brezis.  Functional Analysis, Sobolev Spaces and Partial Differential Equations , 2010 .

[34]  Immanuel M. Bomze,et al.  Copositive optimization - Recent developments and applications , 2012, Eur. J. Oper. Res..

[35]  Ricardo G. Sanfelice,et al.  Hybrid Dynamical Systems: Modeling, Stability, and Robustness , 2012 .

[36]  P. Olver.  Nonlinear Systems , 2013 .

[37]  Bernard Brogliato,et al.  Stability and Observer Design for Lur'e Systems with Multivalued, Nonmonotone, Time-Varying Nonlinearities and State Jumps , 2014, SIAM J. Control. Optim..

[38]  C. Kellett.  Classical Converse Theorems in Lyapunov's Second Method , 2015, 1502.04809.

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

[40]  Luigi Iannelli,et al.  Cone-Copositive Piecewise Quadratic Lyapunov Functions for Conewise Linear Systems , 2015, IEEE Transactions on Automatic Control.

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

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

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

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

[45]  F. Vasca,et al.  Linear Complementarity Systems and Cone-Copositive Lyapunov Stability , 2019, IEEE Control Systems Letters.

[46]  L. Zaccarian,et al.  Smooth Approximation of Patchy Lyapunov Functions for Switched Systems , 2019, IFAC-PapersOnLine.

[47]  Amir Ali Ahmadi,et al.  On Algebraic Proofs of Stability for Homogeneous Vector Fields , 2018, IEEE Transactions on Automatic Control.

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

[49]  L. Iannelli,et al.  Convergence of proximal solutions for evolution inclusions with time-dependent maximal monotone operators , 2019, Mathematical Programming.