A regression-based numerical scheme for backward stochastic differential equations

Based on Fourier cosine expansion, two approximations of conditional expectations are studied, and the local errors for these approximations are analyzed. Using these approximations and the theta-time discretization, a new and efficient numerical scheme, which is based on least-squares regression, for forward–backward stochastic differential equations is proposed. Numerical experiments are done to test the availability and stability of this new scheme for Black–Scholes call and calls combination under an empirical expression about volatility. Some conclusions are given.     

Energy methods for stochastic differential equations

This article is devoted to the existence of strong solutions to stochastic differential equations (SDEs). Compared with Ito's theory, we relax the assumptions on the volatility term and replace the global Lipschitz continuity condition with a local Lipschitz continuity condition and a Hoelder continuity condition. In particular, our general SDE covers the Cox–Ingersoll–Ross SDE as a special case. We note that the general weak existence theory presumably extends to our general SDE (although the explicit time dependence of the drift term and the volatility term might require some extra considerations). However, avoiding weak existence theory we prove the existence of a strong solution directly using a priori estimates (the so-called energy estimates) derived from the SDE. The benefit of this approach is that the argument only requires some basic knowledge about stochastic and functional analysis. Moreover, the underlying principle has developed to become one of the cornerstones of the modern theory of partial differential equations (PDEs). In this sense, the general goal of this article is not just to establish the existence of a strong solution to the SDE under consideration but rather to introduce a new principle in the context of SDEs that has already proven to be successful in the context of PDEs.     

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     

A self-dual variational approach to stochastic partial differential equations

Unlike many of their deterministic counterparts, stochastic partial differential equations are not amenable to the methods of calculus of variations à la Euler–Lagrange. In this paper, we show how self-dual variational calculus leads to variational solutions of various stochastic partial differential equations driven by monotone vector fields. We construct solutions as minima of suitable non-negative and self-dual energy functionals on Itô spaces of stochastic processes. We show how a stochastic version of Bolza's duality leads to solutions for equations with additive noise. We then use a Hamiltonian formulation to construct solutions for non-linear equations with non-additive noise such as the stochastic Navier–Stokes equations in dimension two.     

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.     

Exponential mean square stability of numerical methods for systems of stochastic differential equations

This paper is concerned with exponential mean square stability of the classical stochastic theta method and the so called split-step theta method for stochastic systems. First, we consider linear autonomous systems. Under a sufficient and necessary condition for exponential mean square stability of the exact solution, it is proved that the two classes of theta methods with θ≥0.5$\theta \ge 0.5$ are exponentially mean square stable for all positive step sizes and the methods with θ<0.5$\theta <0.5$ are stable for some small step sizes. Then, we study the stability of the methods for nonlinear non-autonomous systems. Under a coupled condition on the drift and diffusion coefficients, it is proved that the split-step theta method with θ>0.5$\theta >0.5$ still unconditionally preserves the exponential mean square stability of the underlying systems, but the stochastic theta method does not have this property. Finally, we consider stochastic differential equations with jumps. Some similar results are derived.     

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.     

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.     

Two-step numerical methods for parabolic differential equations

A family of two-stepA-stable methods of maximal order for the numerical solution of ordinary differential systems is developed. If these methods are applied to the stiff, large systems which originate from linear parabolic differential equations they yield a large, sparse set of linear algebraic equations of special form. This set is considerably easier to solve than the algebraic equations which are obtained when using diagonal Obrechkoff methods, which are one-step,A-stable and of maximal order     

Three-stage stochastic Runge–Kutta methods for stochastic differential equations

In this paper we discuss three-stage stochastic Runge–Kutta (SRK) methods with strong order 1.0 for a strong solution of Stratonovich stochastic differential equations (SDEs). Higher deterministic order is considered. Two methods, a three-stage explicit (E3) method and a three-stage semi-implicit (SI3) method, are constructed in this paper. The stability properties and numerical results show the effectiveness of these methods in the pathwise approximation of several standard test problems.     

Analysis of multiscale methods for stochastic differential equations

We analyze a class of numerical schemes proposed [26] for stochastic differential equations with multiple time scales. Both advective and diffusive time scales are considered. Weak as well as strong convergence theorems are proven. Most of our results are optimal. They in turn allow us to provide a thorough discussion on the efficiency as well as optimal strategy for the method. © 2005 Wiley Periodicals, Inc.     

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.     

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

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

Box-spline–based approach to the formulation of numerical methods for partial differential equations

A multivariate box spline framework for the formulation of numerical methods for partial differential equations has been constructed. In particular, a fourth-order Galerkin method and a second-order collocation method were derived and applied to a test problem (classical Poisson equation on a square). The examples indicate that accuracy compares favorably with standard methods and the success of iterative schemes suggests an underlying stabilizing effect.     

Mean-square exponential dichotomy of numerical solutions to stochastic differential equations

We present the ability of numerical simulations to reproduce the mean-square exponential dichotomy of stochastic differential equations. Under some conditions, we show that the mean-square exponential dichotomy of stochastic differential equations is equivalent to that of the numerical method for sufficient small step sizes     

A weak second-order split-step method for numerical simulations of stochastic differential equations

In an analogy from symmetric ordinary differential equation numerical integrators, we derive a three-stage, weak 2nd-order procedure for Monte-Carlo simulations of Itô stochastic differential equations. Our composite procedure splits each time step into three parts: an \(h/2\) -stage of trapezoidal rule, an \(h\) -stage martingale, followed by another \(h/2\) -stage of trapezoidal rule. In \(n\) time steps, an \(h/2\) -stage deterministic step follows another \(n-1\) times. Each of these adjacent pairs may be combined into a single \(h\) -stage, effectively producing a two-stage method with partial overlap between successive time steps.     

Trees and numerical methods for ordinary differential equations

This paper presents a review of the role played by trees in the theory of Runge–Kutta methods. The use of trees is in contrast to early publications on numerical methods, in which a deceptively simpler approach was used. This earlier approach is not only non-rigorous, but also incorrect. It is now known, for example, that methods can have different orders when applied to a single equation and when applied to a system of equations; the earlier approach cannot show this. Trees have a central role in the theory of Runge–Kutta methods and they also have applications to more general methods, involving multiple values and multiple stages.