Economical Runge–Kutta methods for numerical solution of stochastic differential equations

For the numerical solution of stochastic differential equations an economical Runge–Kutta scheme of second order in the weak sense is proposed. Numerical stability is studied and some examples are presented to support the theoretical results. AMS subject classification (2000)  60H10     

正倒向随机微分方程组的数值解法

1990年,Pardoux和Peng(彭实戈)解决了非线性倒向随机微分方程(backward stochastic differential equation,BSDE)解的存在唯一性问题,从而建立了正倒向随机微分方程组(forward backward stochastic differential equations,FBSDEs)的理论基础;之后,正倒向随机微分方程组得到了广泛研究,并被应用于众多研究领域中,如随机最优控制、偏微分方程、金融数学、风险度量、非线性期望等.近年来,正倒向随机微分方程组的数值求解研究获得了越来越多的关注,本文旨在基于正倒向随机微分方程组的特性,介绍正倒向随机微分方程组的主要数值求解方法.我们将重点介绍讨论求解FBSDEs的积分离散法和微分近似法,包括一步法和多步法,以及相应的数值分析和理论分析结果.微分近似法能构造出求解全耦合FBSDEs的高效高精度并行数值方法,并且该方法采用最简单的Euler方法求解正向随机微分方程,极大地简化了问题求解的复杂度.文章最后,我们尝试提出关于FBSDEs数值求解研究面临的一些亟待解决和具有挑战性的问题.     

Existence and Hyers-Ulam stability of random impulsive stochastic functional differential equations with finite delays

ABSTRACT

In this paper, we investigate the existence and Hyers-Ulam stability for random impulsive stochastic functional differential equations with finite delays. Firstly, we prove the existence of mild solutions to the equations by using Krasnoselskii's fixed point. Then, we investigate the Hyers-Ulam stability results under the Lipschitz condition on a bounded and closed interval. Finally, an example is given to illustrate our results.     

Reflected backward doubly stochastic differential equations with discontinuous coefficients

In this paper, we study one-dimensional reflected backward doubly stochastic differential equations (RBDSDEs) with one continuous barrier and discontinuous (left or right continuous) generator. We obtain an existence theorem and a comparison theorem for solutions of the class of RBDSDEs.     

Global optimization and stochastic differential equations

Let ⁿ be then-dimensional real Euclidean space,x=(x ₁,x ₂, ...,x _n)^{T n} , and letf: ⁿ R be a real-valued function. We consider the problem of finding the global minimizers off. A new method to compute numerically the global minimizers by following the paths of a system of stochastic differential equations is proposed. This method is motivated by quantum mechanics. Some numerical experience on a set of test problems is presented. The method compares favorably with other existing methods for global optimization.This research has been supported by the European Research Office of the US Army under Contract No. DAJA-37-81-C-0740.The third author gratefully acknowledges Prof. A. Rinnooy Kan for bringing to his attention Ref. 4.     

A generalized existence theorem of backward doubly stochastic differential equations

In this paper, we deal with a class of one-dimensional backward doubly stochastic differential equations （BDSDEs）. We obtain a generalized comparison theorem and a generalized existence theorem of BDSDEs.     

Stationarity and stochastic differential equations

A linear differential equation of order N with stochastic process coefficients and excitation are studied. The objective of this paper is to demonstrate that, by using an expansion method, when the coefficients and excitation are strict sense stationary processes, the response is also a strict sense stationary process. Such problems occur frequently in the engineering sciences and are very important. Applications include parametric random vibrations, turbulent environment rotorcraft dynamics, and dynamics of axially loaded structural members, among others. An example application is provided.     

Global existence of solutions for stochastic impulsive differential equations

In this paper we obtain some results on the global existence of solution to Itô stochastic impulsive differential equations in M([0,∞),? ⁿ ) which denotes the family of ? ⁿ -valued stochastic processes x satisfying sup_t∈[0,∞) \(\mathbb{E}\)|x(t)|² < ∞ under non-Lipschitz coefficients. The Schaefer fixed point theorem is employed to achieve the desired result. An example is provided to illustrate the obtained results.     

Simulation of stochastic differential equations

A lot of discrete approximation schemes for stochastic differential equations with regard to mean-square sense were proposed. Numerical experiments for these schemes can be seen in some papers, but the efficiency of scheme with respect to its order has not been revealed. We will propose another type of error analysis. Also we will show results of simulation studies carried out for these schemes under our notion.     

Stochastic partitioned averaged vector field methods for stochastic differential equations with a conserved quantity

In this paper, stochastic differential equations in the Stratonovich sense with a conserved quantity are considered. A stochastic partitioned averaged vector field method is proposed and analyzed. We prove this numerical method is able to preserve the conserved quantity of the original system. Then the convergence analysis is carried out in detail and we derive the method is convergent with order $1$ in the mean-square sense. Finally some numerical examples are reported to verify the effectiveness and flexibility of the proposed method.     

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.     

On Runge-Kutta-type methods for two-dimensional stochastic differential equations

In the multidimensional case, second-order weak Runge-Kutta methods for stochastic differential equation (SDE) need simulation of correlated random variables, unless the diffusion matrix of SDE satisfies the commutativity condition. In this paper, we show that this can be avoided for some types of diffusion matrices and test functions important for applications. Published in Lietuvos Matematikos Rinkinys, Vol. 46, No. 3, pp. 403–412, July–September, 2006.     

Implicit elliptic differential equations

Let be a bounded domain in ⁿ (n 3) having a smooth boundary, letY be a closed, connected and locally connected subset of ^h , letf be a real-valued function defined on × ^h × ^nh ×Y, and letL be a linear, second-order elliptic operator. In this paper, the existence of strong solutionsu W ^2,p (, ^h ) W ₀ ^1,p (, ^h ) (n<p<+) to the implicit elliptic equationf(x, u, Du, Lu)=0, whereu=(u ₁,u ₂, ...,u _h ),Du=(Du ₁,Du ₂, ...,Du _h ) andLu=(Lu ₁,Lu ₂, ...,Lu _h ), is established. The abstract framework where the equation is studied is that of set-valued analysis.Dedicated to Professor G. Pulvirenti on the occasion of his sixtieth birthday     

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.      

Nonlocal stochastic differential equations with time-varying delay driven by G-Brownian motion

This article studies a class of nonlocal stochastic differential equations driven by G-Brownian motion (G-NSDEs for short). We show the existence and uniqueness results of solutions by means of fixed point theorem. In addition, exponential estimation of (1) has been discussed. Furthermore, we present global solution to Equation (1) with the help of G-Lyapunov functional and ψ-type function.     

Multistep methods for differential algebraic equations

Multistep methods for the differential/algebraic equations (DAEs) in the form of

Existence theory for first order discontinuous functional differential equations

We prove the existence of extremal solutions for a first order functional differential equation subject to nonlinear boundary conditions of functional type. Moreover, the functions that define our problem are allowed to be discontinuous. The proof of our main result is based on a generalized iterative technique.

     

Variational solutions for partial differential equations driven by a fractional noise

In this article we develop an existence and uniqueness theory of variational solutions for a class of nonautonomous stochastic partial differential equations of parabolic type defined on a bounded open subset D⊂R^d and driven by an infinite-dimensional multiplicative fractional noise. We introduce two notions of such solutions for them and prove their existence and their indistinguishability by assuming that the noise is derived from an L²(D)-valued fractional Wiener process W^H with Hurst parameter , whose covariance operator satisfies appropriate integrability conditions, and where γ∈(0,1] denotes the Hölder exponent of the derivative of the nonlinearity in the stochastic term of the equations. We also prove the uniqueness of solutions when the stochastic term is an affine function of the unknown random field. Our existence and uniqueness proofs rest upon the construction and the convergence of a suitable sequence of Faedo-Galerkin approximations, while our proof of indistinguishability is based on certain density arguments as well as on new continuity properties of the stochastic integral we define with respect to W^H.