We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
This work concerns stochastic differential equations with jumps. We prove convergence for solutions to a sequence of (possibly degenerate) stochastic differential equations with jumps when the coefficients converge in some appropriate sense. Then some special cases are analyzed and some concrete and verifiable conditions are given.
This paper focuses on the averaging principle concerning the fast–slow McKean–Vlasov stochastic differential equations driven by mixed fractional Brownian motion with Hurst parameter $\tfrac{1}{2} < H < 1$. The integral associated with Brownian motion is the standard Itô integral, while the integral with respect to fractional Brownian motion is a generalized Riemann–Stieltjes integral. Under the non-Lipschitz condition and certain appropriate assumptions regarding the coefficients, we initially establish the existence and uniqueness theorem for the fast–slow McKean–Vlasov stochastic differential equation driven by mixed fractional Brownian motion. Subsequently, we demonstrate the averaging principle of the fast–slow McKean–Vlasov stochastic differential equations, signifying that the slow stochastic differential equation converges to the associated averaged equation in terms of mean-square convergence.
where $b\,:\, \mathbb{R}^d \rightarrow \mathbb{R}^d$ is a Lipschitz-continuous function, $A \in \mathbb{R}^{d \times d}$ is a positive-definite matrix, $(Z_t)_{t\geqslant 0}$ is a d-dimensional rotationally symmetric $\alpha$-stable Lévy process with $\alpha \in (1,2)$ and $x\in\mathbb{R}^{d}$. We use two Euler–Maruyama schemes with decreasing step sizes $\Gamma = (\gamma_n)_{n\in \mathbb{N}}$ to approximate the invariant measure of $(X_t)_{t \geqslant 0}$: one uses independent and identically distributed $\alpha$-stable random variables as innovations, and the other employs independent and identically distributed Pareto random variables. We study the convergence rates of these two approximation schemes in the Wasserstein-1 distance. For the first scheme, under the assumption that the function b is Lipschitz and satisfies a certain dissipation condition, we demonstrate a convergence rate of $\gamma^{\frac{1}{\alpha}}_n$. This convergence rate can be improved to $\gamma^{1+\frac {1}{\alpha}-\frac{1}{\kappa}}_n$ for any $\kappa \in [1,\alpha)$, provided b has the additional regularity of bounded second-order directional derivatives. For the second scheme, where the function b is assumed to be twice continuously differentiable, we establish a convergence rate of $\gamma^{\frac{2-\alpha}{\alpha}}_n$; moreover, we show that this rate is optimal for the one-dimensional stable Ornstein–Uhlenbeck process. Our theorems indicate that the recent significant result of [34] concerning the unadjusted Langevin algorithm with additive innovations can be extended to stochastic differential equations driven by an $\alpha$-stable Lévy process and that the corresponding convergence rate exhibits similar behaviour. Compared with the result in [6], our assumptions have relaxed the second-order differentiability condition, requiring only a Lipschitz condition for the first scheme, which broadens the applicability of our approach.
We consider time-inhomogeneous ordinary differential equations (ODEs) whose parameters are governed by an underlying ergodic Markov process. When this underlying process is accelerated by a factor $\varepsilon^{-1}$, an averaging phenomenon occurs and the solution of the ODE converges to a deterministic ODE as $\varepsilon$ vanishes. We are interested in cases where this averaged flow is globally attracted to a point. In that case, the equilibrium distribution of the solution of the ODE converges to a Dirac mass at this point. We prove an asymptotic expansion in terms of $\varepsilon$ for this convergence, with a somewhat explicit formula for the first-order term. The results are applied in three contexts: linear Markov-modulated ODEs, randomized splitting schemes, and Lotka–Volterra models in a random environment. In particular, as a corollary, we prove the existence of two matrices whose convex combinations are all stable but are such that, for a suitable jump rate, the top Lyapunov exponent of a Markov-modulated linear ODE switching between these two matrices is positive.
The aim of this article is to study the asymptotic behaviour of non-autonomous stochastic lattice systems. We first show the existence and uniqueness of a pullback measure attractor. Moreover, when deterministic external forcing terms are periodic in time, we show the pullback measure attractors are periodic. We then study the upper semicontinuity of pullback measure attractors as the noise intensity goes to zero. Pullback asymptotic compact for a family of probability measures with respect to probability distributions of the solutions is demonstrated by using uniform a priori estimates for far-field values of solutions.
Continuous-time Markov chains are frequently used to model the stochastic dynamics of (bio)chemical reaction networks. However, except in very special cases, they cannot be analyzed exactly. Additionally, simulation can be computationally intensive. An approach to address these challenges is to consider a more tractable diffusion approximation. Leite and Williams (Ann. Appl. Prob.29, 2019) proposed a reflected diffusion as an approximation for (bio)chemical reaction networks, which they called the constrained Langevin approximation (CLA) as it extends the usual Langevin approximation beyond the first time some chemical species becomes zero in number. Further explanation and examples of the CLA can be found in Anderson et al. (SIAM Multiscale Modeling Simul.17, 2019).
In this paper, we extend the approximation of Leite and Williams to (nearly) density-dependent Markov chains, as a first step to obtaining error estimates for the CLA when the diffusion state space is one-dimensional, and we provide a bound for the error in a strong approximation. We discuss some applications for chemical reaction networks and epidemic models, and illustrate these with examples. Our method of proof is designed to generalize to higher dimensions, provided there is a Lipschitz Skorokhod map defining the reflected diffusion process. The existence of such a Lipschitz map is an open problem in dimensions more than one.
We consider the stochastic volatility model obtained by adding a compound Hawkes process to the volatility of the well-known Heston model. A Hawkes process is a self-exciting counting process with many applications in mathematical finance, insurance, epidemiology, seismology, and other fields. We prove a general result on the existence of a family of equivalent (local) martingale measures. We apply this result to a particular example where the sizes of the jumps are exponentially distributed. Finally, a practical application to efficient computation of exposures is discussed.
We give an extension of Lê’s stochastic sewing lemma. The stochastic sewing lemma proves convergence in $L_m$ of Riemann type sums $\sum _{[s,t] \in \pi } A_{s,t}$ for an adapted two-parameter stochastic process A, under certain conditions on the moments of $A_{s,t}$ and of conditional expectations of $A_{s,t}$ given $\mathcal F_s$. Our extension replaces the conditional expectation given $\mathcal F_s$ by that given $\mathcal F_v$ for $v<s$, and it allows to make use of asymptotic decorrelation properties between $A_{s,t}$ and $\mathcal F_v$ by including a singularity in $(s-v)$. We provide three applications for which Lê’s stochastic sewing lemma seems to be insufficient. The first is to prove the convergence of Itô or Stratonovich approximations of stochastic integrals along fractional Brownian motions under low regularity assumptions. The second is to obtain new representations of local times of fractional Brownian motions via discretization. The third is to improve a regularity assumption on the diffusion coefficient of a stochastic differential equation driven by a fractional Brownian motion for pathwise uniqueness and strong existence.
We study a skew Ornstein–Uhlenbeck process with zero being a sticky reflecting boundary, which is defined as the weak solution to a stochastic differential equation (SDE) system involving local time. The main results obtained include: (i) the existence and uniqueness of solutions to the SDE system, (ii) the scale function and speed measure, and (iii) the distributional properties regarding the transition density and the first hitting times. On the application side, we apply the process to interest rate modeling and obtain the explicit pricing formula for zero-coupon bonds. Numerical examples illustrate the impacts on bond yields of skewness and stickiness parameters.
This paper is concerned with stochastic Schrödinger delay lattice systems with both locally Lipschitz drift and diffusion terms. Based on the uniform estimates and the equicontinuity of the segment of the solution in probability, we show the tightness of a family of probability distributions of the solution and its segment process, and hence the existence of invariant measures on $l^2\times L^2((-\rho,\,0);l^2)$ with $\rho >0$. We also establish a large deviation principle for the solutions with small noise by the weak convergence method.
This research studies the robustness of permanence and the continuous dependence of the stationary distribution on the parameters for a stochastic predator–prey model with Beddington–DeAngelis functional response. We show that if the model is extinct (resp. permanent) for a parameter, it is still extinct (resp. permanent) in a neighbourhood of this parameter. In the case of extinction, the Lyapunov exponent of predator quantity is negative and the prey quantity converges almost to the saturated situation, where the predator is absent at an exponential rate. Under the condition of permanence, the unique stationary distribution converges weakly to the degenerate measure concentrated at the unique limit cycle or at the globally asymptotic equilibrium when the diffusion term tends to 0.
We establish the exponential nonuniform Berry–Esseen bound for the maximum likelihood estimator of unknown drift parameter in an ultraspherical Jacobi process using the change of measure method and precise asymptotic analysis techniques. As applications, the optimal uniform Berry–Esseen bound and optimal Cramér-type moderate deviation for the corresponding maximum likelihood estimator are obtained.
Based on biochemical kinetics, a stochastic model to characterize wastewater treatment plants and dynamics of river water quality under the influence of random fluctuations is proposed in this paper. This model describes the interaction between dissolved oxygen (DO) and biochemical oxygen demand (BOD), and is in the form of stochastic differential equations driven by multiplicative Gaussian noises. The stochastic persistence problem for the model of the system is analysed. Further, a numerical simulation of the stationary probability distributions of BOD and OD by approximations of the stochastic process solution is presented. These results have implications for the prediction and control of pollutants.
We continue with the systematic study of the speed of extinction of continuous-state branching processes in Lévy environments under more general branching mechanisms. Here, we deal with the weakly subcritical regime under the assumption that the branching mechanism is regularly varying. We extend recent results of Li and Xu (2018) and Palau et al. (2016), where it is assumed that the branching mechanism is stable, and complement the recent articles of Bansaye et al. (2021) and Cardona-Tobón and Pardo (2021), where the critical and the strongly and intermediate subcritical cases were treated, respectively. Our methodology combines a path analysis of the branching process together with its Lévy environment, fluctuation theory for Lévy processes, and the asymptotic behaviour of exponential functionals of Lévy processes. Our approach is inspired by the last two previously cited papers, and by Afanasyev et al. (2012), where the analogue was obtained.
This paper analyzes the training process of generative adversarial networks (GANs) via stochastic differential equations (SDEs). It first establishes SDE approximations for the training of GANs under stochastic gradient algorithms, with precise error bound analysis. It then describes the long-run behavior of GAN training via the invariant measures of its SDE approximations under proper conditions. This work builds a theoretical foundation for GAN training and provides analytical tools to study its evolution and stability.
In this paper, we prove a support theorem of Stroock–Varadhan type for pinned diffusion processes. To this end, we use two powerful results from stochastic analysis. One is quasi-sure analysis for Brownian rough path. The other is Aida–Kusuoka–Stroock’s positivity theorem for the densities of weighted laws of non-degenerate Wiener functionals.
We consider parallel single-server queues in heavy traffic with randomly split Hawkes arrival processes. The service times are assumed to be independent and identically distributed (i.i.d.) in each queue and are independent in different queues. In the critically loaded regime at each queue, it is shown that the diffusion-scaled queueing and workload processes converge to a multidimensional reflected Brownian motion in the non-negative orthant with orthonormal reflections. For the model with abandonment, we also show that the corresponding limit is a multidimensional reflected Ornstein–Uhlenbeck diffusion in the non-negative orthant.
Let X be a d-dimensional diffusion and M the running supremum of its first component. In this paper, we show that for any $t>0,$ the density (with respect to the $(d+1)$-dimensional Lebesgue measure) of the pair $\big(M_t,X_t\big)$ is a weak solution of a Fokker–Planck partial differential equation on the closed set $\big\{(m,x)\in \mathbb{R}^{d+1},\,{m\geq x^1}\big\},$ using an integral expansion of this density.
For linear stochastic differential equations with bounded coefficients, we establish the robustness of nonuniform mean-square exponential dichotomy (NMS-ED) on $[t_{0},\,+\infty )$, $(-\infty,\,t_{0}]$ and the whole ${\Bbb R}$ separately, in the sense that such an NMS-ED persists under a sufficiently small linear perturbation. The result for the nonuniform mean-square exponential contraction is also discussed. Moreover, in the process of proving the existence of NMS-ED, we use the observation that the projections of the ‘exponential growing solutions’ and the ‘exponential decaying solutions’ on $[t_{0},\,+\infty )$, $(-\infty,\,t_{0}]$ and ${\Bbb R}$ are different but related. Thus, the relations of three types of projections on $[t_{0},\,+\infty )$, $(-\infty,\,t_{0}]$ and ${\Bbb R}$ are discussed.