On the construction of Lyapunov functions with computer assistance

Computer assisted procedures of Lyapunov functions defined in given neighborhoods of fixed points for flows and maps are discussed. We provide a systematic methodology for constructing explicit ranges where quadratic Lyapunov functions exist in two stages; negative definiteness of associating matrices and direct approach. We note that the former is equivalent to the procedure of cones describing enclosures of the stable and the unstable manifolds of invariant sets, which gives us flexible discussions of asymptotic behavior not only around equilibria for flows but also fixed points for maps. Additionally, our procedure admits a re-parameterization of trajectories in terms of values of Lyapunov functions. Several verification examples are shown for discussions of applicability.

[1]  J. Smoller.  Shock Waves and Reaction-Diffusion Equations , 1983 .

[2]  Ja. G. Sinai,et al.  Discovery of closed orbits of dynamical systems with the use of computers , 1980 .

[3]  Kaname Matsue.  Rigorous numerics for stationary solutions of dissipative PDEs - Existence and local dynamics - , 2013 .

[4]  P. Zgliczynski,et al.  Covering relations, cone conditions and the stable manifold theorem , 2009 .

[5]  Daniel Wilczak.  The Existence of Shilnikov Homoclinic Orbits in the Michelson System: A Computer Assisted Proof , 2006, Found. Comput. Math..

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

[7]  Kazuaki Tanaka,et al.  Numerical validation of blow-up solutions of ordinary differential equations , 2016, J. Comput. Appl. Math..

[8]  Maciej J. Capinski,et al.  Computer Assisted Existence Proofs of Lyapunov Orbits at L2 and Transversal Intersections of Invariant Manifolds in the Jupiter-Sun PCR3BP , 2012, SIAM J. Appl. Dyn. Syst..

[9]  Peter Giesl,et al.  Computation and Verification of Lyapunov Functions , 2015, SIAM J. Appl. Dyn. Syst..

[10]  Nobito Yamamoto,et al.  Some remarks on numerical verification of closed orbits in dynamical systems , 2015 .

[11]  Pablo A. Parrilo,et al.  SOSTOOLS Version 3.00 Sum of Squares Optimization Toolbox for MATLAB , 2013, ArXiv.

[12]  K. Matsue.  Rigorous numerics for fast-slow systems with one-dimensional slow variable: topological shadowing approach , 2015, 1507.01462.

[13]  P. Giesl,et al.  Review on computational methods for Lyapunov functions , 2015 .

[14]  D. Chillingworth.  DYNAMICAL SYSTEMS: STABILITY, SYMBOLIC DYNAMICS AND CHAOS , 1998 .

[15]  Piotr Zgliczynski,et al.  C1 Lohner Algorithm , 2002, Found. Comput. Math..

[16]  H. Kokubu,et al.  Rigorous verification of cocoon bifurcations in the Michelson system , 2007 .

[17]  D. Wilczak,et al.  Topological method for symmetric periodic orbits for maps with a reversing symmetry , 2004, math/0401145.

[18]  Ja. G. Sinaĭ,et al.  Hyperbolicity conditions for the Lorenz model , 1981 .