Revised CPA method to compute Lyapunov functions for nonlinear systems

The CPA method uses linear programming to compute Continuous and Piecewise Affine Lyapunov function for nonlinear systems with asymptotically stable equilibria. In it was shown that the method always succeeds in computing a CPA Lyapunov function for such a system. The size of the domain of the computed CPA Lyapunov function is only limited by the equilibrium’s basin of attraction. However, for some systems, an arbitrary small neighborhood of the equilibrium had to be excluded from the domain a priori. This is necessary, if the equilibrium is not exponentially stable, because the existence of a CPA Lyapunov function in a neighborhood of the equilibrium is equivalent to its exponential stability as shown in. However, if the equilibrium is exponentially stable, then this was an artifact of the method. In this paper we overcome this artifact by developing a revised CPA method. We show that this revised method is always able to compute a CPA Lyapunov function for a system with an exponentially stable equilibrium. The only conditions on the system are that it is C² and autonomous. The domain of the CPA Lyapunov function can be any a priori given compact neighborhood of the equilibrium which is contained in its basin of attraction. Whereas in a previous paper we have shown these results for planar systems, in this paper we cover general n-dimensional systems.

[1]  G. M.,et al.  Partial Differential Equations I , 2023, Applied Mathematical Sciences.

[2]  Antonis Papachristodoulou,et al.  A converse sum-of-squares Lyapunov result: An existence proof based on the Picard iteration , 2010, 49th IEEE Conference on Decision and Control (CDC).

[3]  Peter Giesl,et al.  Construction of Lyapunov functions for nonlinear planar systems by linear programming , 2012 .

[4]  C. Conley.  Isolated Invariant Sets and the Morse Index , 1978 .

[5]  Peter Giesl,et al.  Existence of piecewise affine Lyapunov functions in two dimensions , 2010 .

[6]  Wolfgang Hahn,et al.  Stability of Motion , 1967 .

[7]  Sigurður F. Marinósson,et al.  Lyapunov function construction for ordinary differential equations with linear programming , 2002 .

[8]  Conrad Sanderson,et al.  Armadillo: An Open Source C++ Linear Algebra Library for Fast Prototyping and Computationally Intensive Experiments , 2010 .

[9]  Sigurdur Hafstein,et al.  A constructive converse Lyapunov theorem on asymptotic stability for nonlinear autonomous ordinary differential equations , 2005 .

[10]  Douglas E. Norton.  The fundamental theorem of dynamical systems , 1995 .

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

[12]  John N. Tsitsiklis,et al.  Introduction to linear optimization , 1997, Athena scientific optimization and computation series.

[13]  L. Grüne.  Asymptotic Behavior of Dynamical and Control Systems under Perturbation and Discretization , 2002 .

[14]  Yasuaki Oishi,et al.  Probabilistic design of a robust state-feedback controller based on parameter-dependent Lyapunov functions , 2003, 42nd IEEE International Conference on Decision and Control (IEEE Cat. No.03CH37475).

[15]  C. Hang,et al.  An algorithm for constructing Lyapunov functions based on the variable gradient method , 1970 .

[16]  Tor Arne Johansen,et al.  Computation of Lyapunov functions for smooth nonlinear systems using convex optimization , 2000, Autom..

[17]  Sigurður F. Marinósson,et al.  Stability Analysis of Nonlinear Systems with Linear Programming , 2002 .

[18]  Konstantin Mischaikow,et al.  An Algorithmic Approach to Chain Recurrence , 2005, Found. Comput. Math..

[19]  S. Sastry.  Nonlinear Systems: Analysis, Stability, and Control , 1999 .

[20]  A. M. Lyapunov.  The general problem of the stability of motion , 1992 .

[21]  Mohammad Rezaiee-Pajand,et al.  Estimating the Region of Attraction via collocation for autonomous nonlinear systems , 2012 .

[22]  Sigurdur Hafstein,et al.  A CONSTRUCTIVE CONVERSE LYAPUNOV THEOREM ON EXPONENTIAL STABILITY , 2004 .

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

[24]  Richard H. Bartels,et al.  Algorithm 432 [C2]: Solution of the matrix equation AX + XB = C [F4] , 1972, Commun. ACM.

[25]  Peter Giesl,et al.  Existence of piecewise linear Lyapunov functions in arbitrary dimensions , 2012 .

[26]  V. Zubov.  Methods of A.M. Lyapunov and their application , 1965 .

[27]  Ali Karimpour,et al.  Discontinuous piecewise quadratic Lyapunov functions for planar piecewise affine systems , 2013 .

[28]  A. Papachristodoulou,et al.  On the construction of Lyapunov functions using the sum of squares decomposition , 2002, Proceedings of the 41st IEEE Conference on Decision and Control, 2002..

[29]  P. Julián,et al.  A parametrization of piecewise linear Lyapunov functions via linear programming , 1999 .

[30]  A. Rantzer,et al.  On the computation of piecewise quadratic Lyapunov functions , 1997, Proceedings of the 36th IEEE Conference on Decision and Control.

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

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

[33]  Lars Grüne,et al.  Linear programming based Lyapunov function computation for differential inclusions , 2011 .

[34]  P. Giesl.  Construction of Global Lyapunov Functions Using Radial Basis Functions , 2007 .

[35]  Stefan Ratschan,et al.  Providing a Basin of Attraction to a Target Region of Polynomial Systems by Computation of Lyapunov-Like Functions , 2010, SIAM J. Control. Optim..