Computational Stability Analysis of Lotka-Volterra Systems

Abstract

This paper concerns the computational stability analysis of locally stable Lotka-Volterra (LV) systems by searching for appropriate Lyapunov functions in a general quadratic form composed of higher order monomial terms. The Lyapunov conditions are ensured through the solution of linear matrix inequalities. The stability region is estimated by determining the level set of the Lyapunov function within a suitable convex domain. The paper includes interesting computational results and discussion on the stability regions of higher (3,4) dimensional LV models as well as on the monomial selection for constructing the Lyapunov functions. Finally, the stability region is estimated of an uncertain 2D LV system with an uncertain interior locally stable equilibrium point.