Construction of Lyapunov functions for nonlinear planar systems by linear programming

Recently the authors proved the existence of piecewise affine Lyapunov functions for dynamical systems with an exponentially stable equilibrium in two dimensions (Giesl and Hafstein, 2010 [7], ). Here, we extend these results by designing an algorithm to explicitly construct such a Lyapunov function. We do this by modifying and extending an algorithm to construct Lyapunov functions first presented in Marinosson (2002) [17], and further improved in Hafstein (2007) [10]. The algorithm constructs a linear programming problem for the system at hand, and any feasible solution to this problem parameterizes a Lyapunov function for the system. We prove that the algorithm always succeeds in constructing a Lyapunov function if the system possesses an exponentially stable equilibrium. The size of the region of the Lyapunov function is only limited by the region of attraction of the equilibrium and it includes the equilibrium.