Error analysis of a hybrid method for computing Lyapunov exponents

In a previous paper (Beyn and Lust in Numer Math 113:357–375, 2009) we suggested a numerical method for computing all Lyapunov exponents of a dynamical system by spatial integration with respect to an ergodic measure. The method extended an earlier approach of Aston and Dellnitz (Comput Methods Appl Mech Eng 170:223–237, 1999) for the largest Lyapunov exponent by integrating the diagonal entries from the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$QR$$\end{document}-decomposition of the Jacobian for an iterated map. In this paper we provide an asymptotic error analysis of the method for the case in which all Lyapunov exponents are simple. We employ Oseledec multiplicative ergodic theorem and impose certain hyperbolicity conditions on the invariant subspaces that belong to neighboring exponents. The resulting error expansion shows that one step of extrapolation is enough to obtain exponential decay of errors.