Computing Lyapunov constants for random recurrences with smooth coefficients

Abstract

In recent years, there has been much interest in the growth and decay rates (Lyapunov constants) of solutions to random recurrences such as the random Fibonacci sequence . Many of these problems involve non-smooth dynamics (nondifferentiable invariant measures), making computations hard. Here, however, we consider recurrences with smooth random coefficients and smooth invariant measures. By computing discretised invariant measures and applying Richardson extrapolation, we can compute Lyapunov constants to ten digits of accuracy. In particular, solutions to the recurrence , where the are independent standard normal variables, increase exponentially (almost surely) at the asymptotic rate . Solutions to the related recurrences , and (where the are also independent standard normal variables) increase (decrease) at the rates and respectively.