Closed physical systems eventually come to rest, the reason beingthat due to friction of some kind they continuously lose energy.The mathematical extension of this principle is the concept of aLyapunov function. A Lyapunov function for a dynamical system, ofwhich the dynamics are modelled by an ordinary differentialequation (ODE), is a function that is decreasing along anytrajectory of the system and with exactly one local minimum. Thisimplies that the system must eventually come to rest at thisminimum. Although it has been known for over 50 years that theasymptotic stability of an ODE's equilibrium isequivalent to the existence of a Lyapunov function for the ODE,there has been no constructive method for non-local Lyapunovfunctions, except in special cases. Recently, a novel method toconstruct Lyapunov functions for ODEs via linear programming waspresented [5], [6], which includes an algorithmicdescription of how to derive a linear program for a continuousautonomous ODE, such that a Lyapunov function can be constructedfrom any feasible solution of this linear program. We will showhow to choose the free parameters of this linear program,dependent on the ODE in question, so that it will have a feasiblesolution if the equilibrium at the origin is exponentially stable.This leads to the first constructive converse Lyapunov theorem inthe theory of dynamical systems/ODEs.
Read full abstract