Rigorous numerics in floquet theory: Computing stable and unstable bundles of periodic orbits

Roberto Castelli, Jean Philippe Lessard

Research output: Contribution to JournalArticleAcademicpeer-review

Abstract

In this paper, a rigorous method to compute Floquet normal forms of fundamental matrix solutions of nonautonomous linear differential equations with periodic coefficients is introduced. The Floquet normal form of a fundamental matrix solution F(t) is a canonical decomposition of the form F(t) = Q(t)eRt, where Q(t) is a real periodic matrix and R is a constant matrix. To rigorously compute the Floquet normal form, the idea is to use the regularity of Q(t) and to simultaneously solve for R and Q(t) with the contraction mapping theorem in a Banach space of rapidly decaying coefficients. The explicit knowledge of R and Q can then be used to construct, in a rigorous computer-assisted way, stable and unstable bundles of periodic orbits of vector fields. The new proposed method does not require rigorous numerical integration of the ODE.

Original languageEnglish
Pages (from-to)204-245
Number of pages42
JournalSIAM Journal on Applied Dynamical Systems
Volume12
Issue number1
DOIs
Publication statusPublished - 2013

Keywords

  • Contraction mapping theorem
  • Floquet theory
  • Fundamental matrix solutions
  • Periodic orbits
  • Rigorous numerics
  • Tangent bundles

Fingerprint

Dive into the research topics of 'Rigorous numerics in floquet theory: Computing stable and unstable bundles of periodic orbits'. Together they form a unique fingerprint.

Cite this