Related papers: Estimating generalised Lyapunov exponents for prod…
We give lower and upper bounds on both the Lyapunov exponent and generalised Lyapunov exponents for the random product of positive and negative shear matrices. These types of random products arise in applications such as fluid stirring…
The exact value of the Lyapunov exponents for the random matrix product $P_N = A_N A_{N-1}...A_1$ with each $A_i = \Sigma^{1/2} G_i^{\rm c}$, where $\Sigma$ is a fixed $d \times d$ positive definite matrix and $G_i^{\rm c}$ a $d \times d$…
For products $P_N$ of $N$ random matrices of size $d \times d$, there is a natural notion of finite $N$ Lyapunov exponents $\{\mu_i\}_{i=1}^d$. In the case of standard Gaussian random matrices with real, complex or real quaternion elements,…
Inspired by the latest developments in multilevel Monte Carlo (MLMC) methods and randomised sketching for linear algebra problems we propose a MLMC estimator for real-time processing of matrix structured random data. Our algorithm is…
Monte Carlo method is a broad class of computational algorithms that rely on repeated random sampling to obtain numerical results. They are often used in physical and mathematical problems and are most useful when it is difficult or…
Lyapunov exponents describe the asymptotic behavior of the singular values of large products of random matrices. A direct computation of these exponents is however often infeasible. By establishing a link between Lyapunov exponents and an…
We study the top Lyapunov exponents of random products of positive $2 \times 2$ matrices and obtain an efficient algorithm for its computation. As in the earlier work of Pollicott, the algorithm is based on the Fredholm theory of…
A Monte Carlo method for computing the action of a matrix exponential for a certain class of matrices on a vector is proposed. The method is based on generating random paths, which evolve through the indices of the matrix, governed by a…
We introduce and implement an importance-sampling Monte Carlo algorithm to study systems of globally-coupled oscillators. Our computational method efficiently obtains estimates of the tails of the distribution of various measures of…
This article concerns the non-asymptotic analysis of the singular values (and Lyapunov exponents) of Gaussian matrix products in the regime where $N,$ the number of term in the product, is large and $n,$ the size of the matrices, may be…
We present a sequential Monte Carlo sampler algorithm for the Bayesian analysis of generalised linear mixed models (GLMMs). These models support a variety of interesting regression-type analyses, but performing inference is often extremely…
We prove that the Lyapunov exponents of random products in a (real or complex) matrix group depends continuously on the matrix coefficients and probability weights. More generally, the Lyapunov exponents of the random product defined by any…
The cumulant expansion is used to estimate generalized Lyapunov exponents of the random-frequency harmonic oscillator. Three stochastic processes are considered: Gaussian white noise, Ornstein-Uhlenbeck, and Poisson shot noise. In some…
We discuss several numerical methods for calculating Lyapunov exponents (a quantitative measure of chaos) in systems of ordinary differential equations. We pay particular attention to constrained systems, and we introduce a variety of…
A random matrix with rows distributed as a function of their length is said to be isotropic. When these distributions are Gaussian, beta type I, or beta type II, previous work has, from the viewpoint of integral geometry, obtained the…
The goal of this paper is twofold. In the first part we discuss a general approach to determine Lyapunov exponents from ensemble- rather than time-averages. The approach passes through the identification of locally stable and unstable…
The problems that we consider in this paper are as follows. Let $A_1, \ldots, A_k$ be square matrices (over reals). Let $W=w(A_1, \ldots, A_k)$ be a random product of $n$ matrices. What is the expected absolute value of the largest (in the…
We discuss certain recent metric space methods and some of the possibilities these methods provide, with special focus on various generalizations of Lyapunov exponents originally appearing in the theory of dynamical systems and differential…
A generalized matrix-pencil approach is proposed for the estimation of complex exponential components with segmented signal samples, which is very efficient and provides super-resolution estimations. It is applicable to the signals sampled…
We consider orthogonally invariant probability measures on $\mathrm{GL}_n(\mathbb{R})$ and compare the mean of the logs of the moduli of eigenvalues of the matrices to the Lyapunov exponents of random matrix products independently drawn…