The Numerical Invariant Measure of Stochastic Differential Equations With Markovian Switching

The existence and uniqueness of the numerical invariant measure of the backward Euler-Maruyama method for stochastic differential equations with Markovian switching is yielded, and it is revealed that the numerical invariant measure converges to the underlying invariant measure in the Wasserstein metric. Under the polynomial growth condition of drift term the convergence rate is estimated. The global Lipschitz condition on the drift coefficients required by Bao et al., 2016 and Yuan et al., 2005 is released. Several examples and numerical experiments are given to verify our theory.Comment: 25 pages, 4 figure

Approximate solutions of stochastic differential delay equations with Markovian switching

Our main aim is to develop the existence theory for the solutions to stochastic differential delay equations with Markovian switching (SDDEwMSs) and to establish the convergence theory for the Euler-Maruyama approximate solutions under the local Lipschitz condition. As an application, our results are used to discuss a stochastic delay population system with Markovian switching

Almost sure and moment exponential stability of Euler-Maruyama discretizations for hybrid stochastic differential equations

Positive results are derived concerning the long time dynamics of numerical simulations of stochastic differential equation systems with Markovian switching. Euler-Maruyama discretizations are shown to capture almost sure and momente xponential stability for all sufficiently small timesteps under appropriate conditions

Convergence Rate of Numerical Solutions for Nonlinear Stochastic Pantograph Equations with Markovian Switching and Jumps

The sufficient conditions of existence and uniqueness of the solutions for nonlinear stochastic pantograph equations with Markovian switching and jumps are given. It is proved that Euler-Maruyama scheme for nonlinear stochastic pantograph equations with Markovian switching and Brownian motion is of convergence with strong order 1/2. For nonlinear stochastic pantograph equations with Markovian switching and pure jumps, it is best to use the mean-square convergence, and the order of mean-square convergence is close to 1/2

Switching and diffusion models for gene regulation networks

We analyze a hierarchy of three regimes for modeling gene regulation. The most complete model is a continuous time, discrete state space, Markov jump process. An intermediate 'switch plus diffusion' model takes the form of a stochastic differential equation driven by an independent continuous time Markov switch. In the third 'switch plus ODE' model the switch remains but the diffusion is removed. The latter two models allow for multi-scale simulation where, for the sake of computational efficiency, system components are treated differently according to their abundance. The 'switch plus ODE' regime was proposed by Paszek (Modeling stochasticity in gene regulation: characterization in the terms of the underlying distribution function, Bulletin of Mathematical Biology, 2007), who analyzed the steady state behavior, showing that the mean was preserved but the variance only approximated that of the full model. Here, we show that the tools of stochastic calculus can be used to analyze first and second moments for all time. A technical issue to be addressed is that the state space for the discrete-valued switch is infinite. We show that the new 'switch plus diffusion' regime preserves the biologically relevant measures of mean and variance, whereas the 'switch plus ODE' model uniformly underestimates the variance in the protein level. We also show that, for biologically relevant parameters, the transient behaviour can differ significantly from the steady state, justifying our time-dependent analysis. Extra computational results are also given for a protein dimerization model that is beyond the scope of the current analysis