Table of Contents
Fetching ...

Dynamical systems with fast switching and slow diffusion: Hyperbolic equilibria and stable limit cycles

Nguyen H. Du, Alexandru Hening, Dang H. Nguyen, George Yin

Abstract

We study the long-term qualitative behavior of randomly perturbed dynamical systems. More specifically, we look at limit cycles of stochastic differential equations (SDE) with Markovian switching, in which the process switches at random times among different systems of SDEs, when the switching is fast and the diffusion (white noise) term is small. The system is modeled by $$ dX^{ε,δ}(t)=f(X^{ε,δ}(t), α^ε(t))dt+\sqrtδσ(X^{ε,δ}(t), α^ε(t))dW(t) , \ X^ε(0)=x, $$ where $α^ε(t)$ is a finite state space Markov chain with irreducible generator $Q=(q_{ij})$. The relative changing rates of the switching and the diffusion are highlighted by the two small parameters $ε$ and $δ$. We associate to the system the averaged ODE \[ d\bar X(t)=\bar f(\bar X(t))dt, \ X(0)=x, \] where $\bar f(\cdot)=\sum_{i=1}^{m_0}f(\cdot, i)ν_i$ and $(ν_1,\dots,ν_{m_0})$ is the unique invariant probability measure of the Markov chain with generator $Q$. Suppose that for each pair $(ε,δ)$ of parameters, the process has an invariant probability measure $μ^{ε,δ}$, and that the averaged ODE has a limit cycle in which there is an averaged occupation measure $μ^0$ for the averaged equation. We are able to prove that if $\bar f$ has finitely many unstable or hyperbolic fixed points, then $μ^{ε,δ}$ converges weakly to $μ^0$ as $ε\to 0$ and $δ\to 0$. Our results generalize to the setting of state-dependent switching \[ \mathbb{P}\{α^ε(t+Δ)=j~|~α^ε=i, X^{ε,δ}(s),α^ε(s), s\leq t\}=q_{ij}(X^{ε,δ}(t))Δ+o(Δ),~~ i\neq j \] as long as the generator $Q(\cdot)=(q_{ij}(\cdot))$ is bounded, Lipschitz, and irreducible for all $x\in\mathbb{R}^d$. We conclude our analysis by studying a predator-prey model.

Dynamical systems with fast switching and slow diffusion: Hyperbolic equilibria and stable limit cycles

Abstract

We study the long-term qualitative behavior of randomly perturbed dynamical systems. More specifically, we look at limit cycles of stochastic differential equations (SDE) with Markovian switching, in which the process switches at random times among different systems of SDEs, when the switching is fast and the diffusion (white noise) term is small. The system is modeled by where is a finite state space Markov chain with irreducible generator . The relative changing rates of the switching and the diffusion are highlighted by the two small parameters and . We associate to the system the averaged ODE where and is the unique invariant probability measure of the Markov chain with generator . Suppose that for each pair of parameters, the process has an invariant probability measure , and that the averaged ODE has a limit cycle in which there is an averaged occupation measure for the averaged equation. We are able to prove that if has finitely many unstable or hyperbolic fixed points, then converges weakly to as and . Our results generalize to the setting of state-dependent switching as long as the generator is bounded, Lipschitz, and irreducible for all . We conclude our analysis by studying a predator-prey model.

Paper Structure

This paper contains 10 sections, 19 theorems, 210 equations, 5 figures.

Key Result

Theorem 1.1

Suppose Assumptions asp1 and asp2 hold. The family of invariant probability measures $(\mu^{\varepsilon,\delta})_{\varepsilon>0}$ converges weakly to the measure $\mu^0$ given by e:mu0 in the sense that for every bounded and continuous function $g:\mathbb{R}^d\times\mathcal{M}\to \mathbb{R}$, where $T_\Gamma$ is the period of the limit cycle, $y\in\Gamma$ and $\overline g(x)=\sum_{i\in\mathcal{M}

Figures (5)

  • Figure 1: From left to right: Graphs of the $x^{\varepsilon,\delta}(t)$ component of \ref{['nex1']} with $(\varepsilon,\delta)=(0.001, 0.001)$, $(\varepsilon,\delta)=(0.00005, 0.00005)$ and $x(t)$ of the averaged system \ref{['nex2']} respectively.
  • Figure 2: From left to right: Graphs of the $y^{\varepsilon,\delta}(t)$ component of \ref{['nex1']} with $(\varepsilon,\delta)=(0.001, 0.001)$, $(\varepsilon,\delta)=(0.00005, 0.00005)$ and $y(t)$ of the averaged system \ref{['nex2']} respectively.
  • Figure 3: Phase portraits of \ref{['nex1']} for different values of $\varepsilon$ and $\delta$.
  • Figure 4: The two figures illustrate an approximate density of the occupation measure on the time interval $[0,1000]$ with $\varepsilon=\delta=0.001$ on the right and $\varepsilon=\delta=0.0001$ on the left. The step size is $h=0.0001$. We divide the domain in $50\times 50$ cells, approximate the density by the frequency of the process staying in each cell, and then interpolate. The simulations support the theoretical results that the occupation measures converge to the invariant probability measure.
  • Figure 5: The left figure is a sample path starting at the equilibrium $(0,0,0)$ of the limit system. We see that this path has a tendency of approaching $(0,0,1)$ before it is finally attracted to the proximity of the limit cycle. The right figure is an approximated density of the occupation measure of two components $x,y$ on the time interval $[0,1000]$.

Theorems & Definitions (42)

  • Remark 1.1
  • Remark 1.2
  • Remark 1.3
  • Theorem 1.1
  • Remark 1.4
  • Remark 1.5
  • Theorem 1.2
  • Remark 1.6
  • Lemma 2.1
  • Lemma 2.2
  • ...and 32 more