Ulam's method for computing stationary densities of invariant measures for piecewise convex maps with countably infinite number of branches

Md Shafiqul Islam,Paweł Góra,A H M Mahbubur Rahman
2024-05-05
Abstract:Let $\tau: I=[0, 1]\to [0, 1]$ be a piecewise convex map with countably infinite number of branches. In \cite{GIR}, the existence of absolutely continuous invariant measure (ACIM) $\mu$ for $\tau$ and the exactness of the system $(\tau, \mu)$ has been proven. In this paper, we develop an Ulam method for approximation of $f^*$, the density of ACIM $\mu$. We construct a sequence $\{\tau_n\}_{n=1}^\infty$ of maps $\tau_n: I\to I$ s. t. $\tau_n$ has a finite number of branches and the sequence $\tau_n$ converges to $\tau$ almost uniformly. Using supremum norms and Lasota-Yorke type inequalities, we prove the existence of ACIMs $\mu_n$ for $\tau_n$ with the densities $f_n$. For a fixed $n$, we apply Ulam's method with $k$ subintervals to $\tau_n$ and compute approximations $f_{n,k}$ of $f_n$. We prove that $f_{n,k}\to f^*$ as $n\to \infty, k\to \infty,$ both a.e. and in $L^1$. We provide examples of piecewise convex maps $\tau$ with countably infinite number of branches, their approximations $\tau_n$ with finite number of branches and for increasing values of parameter $k$ show the errors $\|f^*-f_{n,k}\|_1$.
Dynamical Systems
What problem does this paper attempt to address?
This paper aims to solve the problem of how to calculate the stationary density of the invariant measure of piecewise convex mappings with countably infinite branches. Specifically, the paper focuses on the piecewise convex mapping $\tau$ defined on the interval $I = [0,1]$, which has countably infinitely many branches. Previous work has proven that such mappings have absolutely continuous invariant measures (ACIM), and the system $(\tau, \mu)$ is exact. However, there is still a lack of effective methods for effectively approximating the density function $f^*$ of these invariant measures, especially when the mapping has infinitely many branches. ### Main contributions 1. **Method development**: The paper proposes a numerical algorithm based on the Ulam method for approximately calculating the density function $f^*$ of the invariant measure of the piecewise convex mapping $\tau$. 2. **Theoretical proof**: The author constructs a sequence of mappings $\{\tau_n\}_{n = 1}^{\infty}$, where each $\tau_n$ has a finite number of branches, and this sequence converges almost uniformly to $\tau$. By using the supremum norm and the Lasota - Yorke - type inequality, it is proven that $\tau_n$ has an invariant measure $\mu_n$ with a density function $f_n$. 3. **Convergence analysis**: Further application of the Ulam method, by dividing $\tau_n$ into $k$ sub - intervals to calculate the approximate value $f_{n,k}$ of $f_n$, and it is proven that $f_{n,k}$ converges to $f^*$ almost everywhere and in the $L^1$ norm as $n$ and $k$ tend to infinity. ### Key results - **Construction of approximate mappings**: By constructing a sequence of mappings $\{\tau_n\}$ with a finite number of branches, such that $\tau_n$ converges almost uniformly to the original mapping $\tau$. - **Application of the Ulam method**: Use the Ulam method to approximate each $\tau_n$ to obtain the approximate density function $f_{n,k}$. - **Convergence proof**: Prove that $f_{n,k}$ converges to $f^*$ almost everywhere and in the $L^1$ norm as $n$ and $k$ tend to infinity. ### Numerical experiments The paper also provides specific numerical experiments to demonstrate the effectiveness of this method. By calculating the approximate density functions $f_{n,k}$ under different parameters $n$ and $k$, the correctness of the theoretical results is verified, and the trend of the approximate error gradually decreasing as $n$ and $k$ increase is shown. ### Conclusion This paper successfully solves the problem of how to calculate the density function of the invariant measure of piecewise convex mappings with countably infinite branches, provides an effective numerical method, and verifies the effectiveness of the method through strict mathematical proofs and numerical experiments. This provides an important tool for studying the long - term behavior and chaotic properties of such mappings.