SDELab: A package for solving stochastic differential equations in MATLAB

We introduce SDELab, a package for solving stochastic differential equations (SDEs) within MATLAB. SDELab features explicit and implicit integrators for a general class of Itô and Stratonovich SDEs, including Milstein's method, sophisticated algorithms for iterated stochastic integrals, and flexible plotting facilities.     

Pathwise accuracy and ergodicity of metropolized integrators for SDEs

Metropolized integrators for ergodic stochastic differential equations (SDEs) are proposed that (1) are ergodic with respect to the (known) equilibrium distribution of the SDEs and (2) approximate pathwise the solutions of the SDEs on finite‐time intervals. Both these properties are demonstrated in the paper, and precise strong error estimates are obtained. It is also shown that the Metropolized integrator retains these properties even in situations where the drift in the SDE is nonglobally Lipschitz, and vanilla explicit integrators for SDEs typically become unstable and fail to be ergodic. © 2009 Wiley Periodicals, Inc.     

Weak local linear discretizations for stochastic differential equations: Convergence and numerical schemes

Weak local linear (WLL) discretizations are playing an increasing role in the construction of effective numerical integrators and inference methods for stochastic differential equations (SDEs) with additive noise. However, due to limitations in the existing numerical implementations of WLL discretizations, the resulting integrators and inference methods have either been restricted to particular classes of autonomous SDEs or showed low computational efficiency. Another limitation is the absence of a systematic theoretical study of the rate of convergence of the WLL discretizations and numerical integratos. This task is the main purpose of the present paper. A second goal is introducing a new WLL scheme that overcomes the numerical limitations mentioned above. Additionally, a comparative analysis between the new WLL scheme and some conventional weak integrators is also presented.     

Different types of spdes in the eyes of girsanov's theorem

We prove Girsanov's theorem for continuous orthogonal martingale measures. We then define space-time SDEs, and use Girsanov's theorem to establish a oneto- one correspondence between solutions of two space-time SDEs differing only by a drift coefficient. For such stochastic equations, we give necessary conditions under which the laws of their solutions are absolutely continuous with respect to each other. Using Girsanov's theorem again, we prove additional existence and uniqueness results for space-time SDEs. The same one-to-one correspondence and absolute continuity theorems are also proved for the stochastic heat and wave equations     

Strong comparison result for a class of re-ected stochastic differential equations with non-Lipschitzian coeffcients

In this paper, we are concerned with a class of reflected stochastic differential equations (reflected SDEs) with non-Lipschitzian coeffcients. Under the same coeffcients assumptions as Fang and Zhang [Probab. Theory Relat. Fields, 2005, 132(3): 356 390] for a class of SDEs, we establish the pathwise uniqueness for the reflected SDEs. Furthermore, a strong comparison theorem is proved for the reflected SDEs in a onedimensional case.     

非Lipschitz条件下半鞅随机微分方程解的唯一性(英文)

本文研究了非Lipschitz条件下半鞅随机微分方程.利用It(o)分析和Gronwall不等式,探讨了随机微分方程无爆炸解,并证明了随机微分方程解的唯一性.     

Strong comparison result for a class of reflected stochastic differential equations with non-Lipschitzian coefficients

In this paper, we are concerned with a class of reflected stochastic differential equations (reflected SDEs) with non-Lipschitzian coefficients. Under the same coefficients assumptions as Fang and Zhang [Probab. Theory Relat. Fields, 2005, 132(3): 356–390] for a class of SDEs, we establish the pathwise uniqueness for the reflected SDEs. Furthermore, a strong comparison theorem is proved for the reflected SDEs in a one-dimensional case.      

Volume preserving RK methods for linear systems

1. IntroductionIn recent yearss there has been a great interest in constructing numerical integrationschemes for ODEs in such a way that some qualitative geometrical properties of the solutionof the ODEs are exactly preserved. R.th[ll and Feng Kang[2'31 has proposed symplectic algorithms for Hamiltollian systems, and since then st ruct ure s- preserving me t ho ds fordynamical systems have been systematically developed[4--7]. The symplectic algorithms forHamiltonian systems, the volume-pre…     

RELATIONS BETWEEN SOLUTIONS TO STOCHASTIC DIFFERENTIAL EQUATIONS DRIVEN BY SEMIMARTINGALE WITH NON-LIPSCHITZ COEFFICIENTS

In this paper, a class of stochastic differential equations (SDEs) driven by semi-martingale with non-Lipschitz coefficients is studied. We investigate the dependence of solutions to SDEs on the initial value. To obtain a continuous version, we impose the conditions on the local characteristic of semimartingale. In this case, it gives rise to a flow of homeomorphisms if the local characteristic is compactly supported.     

High strong order stochastic Runge-Kutta methods for Stratonovich stochastic differential equations with scalar noise

This paper concerns the stochastic Runge-Kutta (SRK) methods with high strong order for solving the Stratonovich stochastic differential equations (SDEs) with scalar noise. Firstly, the new SRK methods with strong order 1.5 or 2.0 for the Stratonovich SDEs with scalar noise are constructed by applying colored rooted tree analysis and the theorem of order conditions for SRK methods proposed by Rößler (SIAM J. Numer. Anal. 48(3), 922–952, 2010). Secondly, a specific SRK method with strong order 2.0 for the Stratonovich SDEs whose drift term vanishes is proposed. And another specific SRK method with strong order 1.5 for the Stratonovich SDEs whose drift and diffusion terms satisfy the commutativity condition is proposed. The two specific SRK methods need only to use one random variable and do not need to simulate the multiple Stratonovich stochastic integrals. Finally, the numerical results show that performance of our methods is better than those of well-known SRK methods with strong order 1.0 or 1.5.     

DISCRETIZATION OF JUMP STOCHASTIC DIFFERENTIAL EQUATIONS IN TERMS OF MULTIPLE STOCHASTIC INTEGRALS

1.IntroductionFOrthestrongdiscretizationofSDEs,anynumericalmethodwhichonlydependsonthevaluesofBrownianpathsorPoissonpathsatthepartitionnodescannotachieveanorderhigherthan0.5ingeneral[')'1'].Thereforetheevaluationofmultiplestochasticintegralsontheintervalsbetweennodesisamajorobstaclethatmustbeovercome.Someattemptshavebeenmadepreviouslyindifferentapproachestoapproximatemul-tiplestochasticintegrals.[2]suggestsanapproximationintermsofFourierGaussiancoefficientsoftheBrownianbridgeprocess.Asthel…     

BSDEs and SDEs with time-advanced and -delayed coefficients

ABSTRACT

This paper introduces a class of backward stochastic differential equations (BSDEs), whose coefficients not only depend on the value of its solutions of the present but also the past and the future. For a sufficiently small time delay or a sufficiently small Lipschitz constant, the existence and uniqueness of such BSDEs is obtained. As an adjoint process, a class of stochastic differential equations (SDEs) is introduced, whose coefficients also depend on the present, the past and the future of its solutions. The existence and uniqueness of such SDEs is proved for a sufficiently small time advance or a sufficiently small Lipschitz constant. A duality between such BSDEs and SDEs is established.     

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.     

SOLVXBILITY OF FORWARD－BACKWARD SDES AND THE NODAL SET OF HAMILTON－JACOBI－BELLMAN EQUATIONS

ＳＯＬＶＸＢＩＬＩＴＹＯＦＦＯＲＷＡＲＤ－ＢＡＣＫＷＡＲＤＳＤＥＳＡＮＤＴＨＥＮＯＤＡＬＳＥＴＯＦＨＡＭＩＬＴＯＮ－ＪＡＣＯＢＩ－ＢＥＬＬＭＡＮＥＱＵＡＴＩＯＮＳ￥ＭＡＪＩＮ；ＹＯＮＧＪＩＯＮＧＭＩＮＡｂｓｔｒａｃｔ：Ｔｈｅｓｏｌｖａｂｉｌｉｔｙｏｆ...     

SOLVABILITY OF FORWARD-BACKWARD SDES AND THE NODAL SET OF HAMILTON-JACOBI-BELLMAN EQUATIONS

The solvability of a class of forward-backward stochastic differential equations (SDEs for short) over an arbitrarily prescribed time duration is studied. The authors design a stochastic relaxed control problem, with both drift and difftusion all being controlled, so that the solvability problem is converted to a problem of finding the nodal set of the viscosity solution to a certain Hamilton-Jacobi-Bellman equation. This method overcomes the fatal difficulty encountered in the traditional contraction mapping approach to the existence theorem of such SDEs.     

Mean square polynomial stability of numerical solutions to a class of stochastic differential equations

The exponential stability of numerical methods to stochastic differential equations (SDEs) has been widely studied. In contrast, there are relatively few works on polynomial stability of numerical methods. In this letter, we address the question of reproducing the polynomial decay of a class of SDEs using the Euler–Maruyama method and the backward Euler–Maruyama method. The key technical contribution is based on various estimates involving the gamma function.     

High order local linearization methods: An approach for constructing A-stable explicit schemes for stochastic differential equations with additive noise

An approach for the construction of A-stable high order explicit strong schemes for stochastic differential equations (SDEs) with additive noise is proposed. We prove that such schemes also have the dynamical property that we call Random A-stability (RA-stability), which ensures that, for linear equations with stationary solutions, the numerical scheme has a random attractor that converges to the exact one as the step size decreases. Basically, the proposed integrators are obtained by splitting, at each time step, the solution of the original equation into two parts: the solution of a linear ordinary differential equation plus the solution of an auxiliary SDE. The first one is solved by the Local Linearization scheme in such a way that A-stability is guaranteed, while the second one is approximated by any extant scheme, preferably an explicit one that yields high order of convergence with low computational cost. Numerical integrators constructed in this way are called High Order Local Linearization (HOLL) methods. Various efficient HOLL schemes are elaborated in detail, and their performance is illustrated through computer simulations. Furthermore, mean-square convergence of the introduced methods is studied.     

半鞅非Lipschitz系数随机微分方程解的大偏差

建立了半鞅非Lipschitz系数随机微分方程, 研究了Freidlin-Wentzell型大偏差原理.     

A class of new Magnus-type methods for semi-linear non-commutative Itô stochastic differential equations

Numerical Algorithms - In this paper, a class of new Magnus-type methods is proposed for non-commutative Itô stochastic differential equations (SDEs) with semi-linear drift term and...     

SDEs Driven by a Time-Changed Lévy Process and Their Associated Time-Fractional Order Pseudo-Differential Equations

It is known that the transition probabilities of a solution to a classical It? stochastic differential equation (SDE) satisfy in the weak sense the associated Kolmogorov equation. The Kolmogorov equation is a partial differential equation with coefficients determined by the corresponding SDE. Time-fractional Kolmogorov-type equations are used to model complex processes in many fields. However, the class of SDEs that is associated with these equations is unknown except in a few special cases. The present paper shows that in the cases of either time-fractional order or more general time-distributed order differential equations, the associated class of SDEs can be described within the framework of SDEs driven by semimartingales. These semimartingales are time-changed Lévy processes where the independent time-change is given respectively by the inverse of a single or mixture of independent stable subordinators. Examples are provided, including a fractional analogue of the Feynman–Kac formula.