CONVERGENCE ANALYSIS OF PARAREAL ALGORITHM BASED ON MILSTEIN SCHEME FOR STOCHASTIC DIFFERENTIAL EQUATIONS

In this paper, we propose a parareal algorithm for stochastic differential equations (SDEs), which proceeds as a two-level temporal parallelizable integrator with the Milstein scheme as the coarse propagator and the exact solution as the fine propagator. The convergence order of the proposed algorithm is analyzed under some regular assumptions. Finally, numerical experiments are dedicated to illustrating the convergence and the convergence order with respect to the iteration number $k$, which show the efficiency of the proposed method.     

Implicit Milstein method for stochastic differential equations via the Wong-Zakai approximation

The aim of this paper is to derive a numerical scheme for solving stochastic differential equations (SDEs) via Wong-Zakai approximation. One of the most important methods for solving SDEs is Milstein method, but this method is not so popular because the cost of simulating the double stochastic integrals is high. For overcoming this complexity, we present an implicit Milstein scheme based on Wong-Zakai approximation by approximating the Brownian motion with its truncated Haar expansion. The main advantages of this method lie in the fact that it preserves the convergence order and also stability region of the Milstein method while its simulation is much easier than Milstein scheme. We show the convergence rate of the method by some numerical examples.     

The tamed Milstein method for commutative stochastic differential equations with non-globally Lipschitz continuous coefficients

For stochastic differential equations (SDEs) with a superlinearly growing and globally one-sided Lipschitz continuous drift coefficient, the classical explicit Euler scheme fails to converge strongly to the exact solution. Recently, an explicit strongly convergent numerical scheme, called the tamed Euler method, has been proposed in [8] for such SDEs. Motivated by their work, we here introduce a tamed version of the Milstein scheme for SDEs with commutative noise. The proposed method is also explicit and easily implementable, but achieves higher strong convergence order than the tamed Euler method does. In recovering the strong convergence order one of the new method, new difficulties arise and kind of a bootstrap argument is developed to overcome them. Finally, an illustrative example confirms the computational efficiency of the tamed Milstein method compared to the tamed Euler method.     

Mean-square stability of Milstein method for linear hybrid stochastic delay integro-differential equations

In this paper we study the mean-square (MS) stability of the Milstein method for linear stochastic delay integro-differential equations (SDIDE) with Markovian switching by extending the techniques of [Z. Wang, C. Zhang, An analysis of stability of Milstein method for stochastic differential equations with delay, Computers and Mathematics with Applications 51 (2006) 1445–1452; L. Ronghua, H. Yingmin, Convergence and stability of numerical solutions to SDDEs with Markovian switching, Applied Mathematics and Computation 175 (2006) 1080–1091]. It is established that the Milstein method is MS-stable for linear stochastic delay differential equations (Wang and Zhang (2006); in the above reference). Here we prove that it is MS-stable for linear SDIDE with Markovian switching also under suitable conditions on the integral term. A numerical example is provided to illustrate the theoretical results.     

The composite Milstein methods for the numerical solution of Stratonovich stochastic differential equations

In this paper, we present two composite Milstein methods for the strong solution of Stratonovich stochastic differential equations driven by d-dimensional Wiener processes. The composite Milstein methods are a combination of semi-implicit and implicit Milstein methods. The criterion for choosing either the implicit or the semi-implicit method at each step of the numerical solution is given. The stability and convergence properties of the proposed methods are analyzed for the linear test equation. It is shown that the proposed methods converge to the exact solution in Stratonovich sense. In addition, the stability properties of our methods are found to be superior to those of the Milstein and the composite Euler methods. The convergence properties for the nonlinear case are shown numerically to be the same as the linear case. Hence, the proposed methods are a good candidate for the solution of stiff SDEs.     

Newton–Milstein scheme for stochastic differential equations and its fast uniform convergence

We show that a modified Milstein scheme combined with explicit Newton’s method enables us to construct fast converging sequences of approximate solutions of stochastic differential equations. The fast uniform convergence of our Newton–Milstein scheme follows from Amano’s probabilistic second-order error estimate, which had been an open problem since 1991. The Newton–Milstein scheme, which is based on a modified Milstein scheme and the symbolic Newton’s method, will be classified as a numerical and computer algebraic hybrid method and it may give a new possibility to the study of computer algebraic method in stochastic analysis.     

The composite Milstein methods for the numerical solution of Ito stochastic differential equations

In this paper, we present the composite Milstein methods for the strong solution of Ito stochastic differential equations. These methods are a combination of semi-implicit and implicit Milstein methods. We give a criterion for choosing either the implicit or the semi-implicit scheme at each step of our numerical solution. The stability and convergence properties are investigated and discussed for the linear test equation. The convergence properties for the nonlinear case are shown numerically to be the same as the linear case. The stability properties of the composite Milstein methods are found to be more superior compared to those of the Milstein, the Euler and even better than the composite Euler method. This superiority in stability makes the methods a better candidate for the solution of stiff SDEs.     

随机延迟微分方程的Milstein方法的非线性均方稳定性

本文针对一般的非线性随机延迟微分方程,证明了当系统理论解满足均方稳定性条件时,则当方程的漂移和扩散项满足一定的条件时,Milstein方法也是均方稳定的.数学实验进一步验证了我们的结论.     

中立型随机延迟微分方程分裂步θ方法的强收敛性

中立型随机延迟微分方程常出现在一些科学技术和工程领域中.本文在漂移系数和扩散系数关于非延迟项满足全局Lipschitz条件,关于延迟项满足多项式增长条件以及中立项满足多项式增长条件下,证明了分裂步θ方法对于中立型随机延迟微分方程的强收敛阶为1/2.数值实验也验证了这一理论结果.     

求解随机微分方程的二级Milstein方法(英文)

In this paper we discuss two-stage Miistein methods for solving Ito stochastic differential equations （SDEs）. Six fully explicit methods （TSM 1 -- TSM 6） are given in this paper. Their order of strong convergence is proved. The stability properties and numerical results show the effectiveness of these methods in the pathwise approximation of Ito SDEs.     

Continuous weak approximation for stochastic differential equations

A convergence theorem for the continuous weak approximation of the solution of stochastic differential equations (SDEs) by general one-step methods is proved, which is an extension of a theorem due to Milstein. As an application, uniform second order conditions for a class of continuous stochastic Runge–Kutta methods containing the continuous extension of the second order stochastic Runge–Kutta scheme due to Platen are derived. Further, some coefficients for optimal continuous schemes applicable to Itô SDEs with respect to a multi–dimensional Wiener process are presented.     

改进欧拉格式求解随机微分方程的收敛性

采用改进的欧拉格式求解随机微分方程,当方程的偏移系数和扩散系数均满足全局Lipschitz条件和线性增长条件时,证明改进格式的强收敛的阶是1/2.     

Strong convergence of a stochastic Rosenbrock-type scheme for the finite element discretization of semilinear SPDEs driven by multiplicative and additive noise

This paper aims to investigate the numerical approximation of a general second order parabolic stochastic partial differential equation(SPDE) driven by multiplicative and additive noise. Our main interest is on such SPDEs where the nonlinear part is stronger than the linear part, usually called stochastic dominated transport equations. Most standard numerical schemes lose their good stability properties on such equations, including the current linear implicit Euler method. We discretize the SPDE in space by the finite element method and propose a novel scheme called stochastic Rosenbrock-type scheme for temporal discretization. Our scheme is based on the local linearization of the semi-discrete problem obtained after space discretization and is more appropriate for such equations. We provide a strong convergence of the new fully discrete scheme toward the exact solution for multiplicative and additive noise and obtain optimal rates of convergence. Numerical experiments to sustain our theoretical results are provided.     

随机延迟微分方程平衡方法的均方收敛性与稳定性

本文讨论求解刚性随机延迟微分方程的平衡方法.证明了随机延迟微分方程平衡方法的均方收敛阶为1/2.给出了线性随机延迟微分方程平衡方法均方稳定的条件.     

On the strong Feller property for stochastic delay differential equations with singular drift

In this paper, we prove the strong Feller property for stochastic delay (or functional) differential equations with singular drift. We extend an approach of Maslowski and Seidler to derive the strong Feller property of those equations, see Maslowski and Seidler (2000). The argumentation is based on the well-posedness and the strong Feller property of the equations’ drift-free version. To this aim, we investigate a certain convergence of random variables in topological spaces in order to deal with discontinuous drift coefficients.     

Split-step forward methods for stochastic differential equations

In this paper we discuss split-step forward methods for solving Itô stochastic differential equations (SDEs). Eight fully explicit methods, the drifting split-step Euler (DRSSE) method, the diffused split-step Euler (DISSE) method and the three-stage Milstein (TSM 1a-TSM 1f) methods, are constructed based on Euler-Maruyama method and Milstein method, respectively, in this paper. Their order of strong convergence is proved. The analysis of stability shows that the mean-square stability properties of the methods derived in this paper are improved on the original methods. The numerical results show the effectiveness of these methods in the pathwise approximation of Itô SDEs.     

Sufficient conditions for the eventual strong Feller property for degenerate stochastic evolutions

The strong Feller property is an important quality of Markov semigroups which helps for example in establishing uniqueness of invariant measure. Unfortunately degenerate stochastic evolutions, such as stochastic delay equations, do not possess this property. However the eventual strong Feller property is sufficient in establishing uniqueness of invariant probability measure. In this paper we provide operator theoretic conditions under which a stochastic evolution equation with additive noise possesses the eventual strong Feller property. The results are used to establish uniqueness of invariant probability measure for stochastic delay equations and stochastic partial differential equations with delay, with an application in neural networks.     

带有弱奇性核的多项分数阶非线性随机微分方程的改进Euler-Maruyama格式

针对一类带有弱奇性核的多项分数阶非线性随机微分方程构造了改进Euler-Maruyama (EM)格式,并证明了该格式的强收敛性.具体地,利用随机积分解的充分条件,将此多项分数阶随机微分方程等价地转化为随机Volterra 积分方程的形式,详细推导出对应的改进EM格式,并对该格式进行了强收敛性分析,其强收敛阶为α_m-α_m-1,其中α_i为分数阶导数的指标,且满足0<α₁<…<α_m-1<α_m<1.最后,通过数值实验验证了理论分析结果的正确性.     

A Runge–Kutta type scheme for nonlinear stochastic partial differential equations with multiplicative trace class noise

In this paper a new Runge–Kutta type scheme is introduced for nonlinear stochastic partial differential equations (SPDEs) with multiplicative trace class noise. The proposed scheme converges with respect to the computational effort with a higher order than the well-known linear implicit Euler scheme. In comparison to the infinite dimensional analog of Milstein type scheme recently proposed in Jentzen and Röckner (2012), our scheme is easier to implement and needs less computational effort due to avoiding the derivative of the diffusion function. The new scheme can be regarded as an infinite dimensional analog of Runge–Kutta method for finite dimensional stochastic ordinary differential equations (SODEs). Numerical examples are reported to support the theoretical results.     

A family of fully implicit Milstein methods for stiff stochastic differential equations with multiplicative noise

In this paper a family of fully implicit Milstein methods are introduced for solving stiff stochastic differential equations (SDEs). It is proved that the methods are convergent with strong order 1.0 for a class of SDEs. For a linear scalar test equation with multiplicative noise terms, mean-square and almost sure asymptotic stability of the methods are also investigated. We combine analytical and numerical techniques to get insights into the stability properties. The fully implicit methods are shown to be superior to those of the corresponding semi-implicit methods in term of stability property. Finally, numerical results are reported to illustrate the convergence and stability results.