Complex Boundary Value Problems of Nonlinear Differential Equations: Theory, Computational Methods, and Applications
View this Special IssueResearch Article  Open Access
Yanli Zhou, Yonghong Wu, Xiangyu Ge, B. Wiwatanapataphee, "A Robust Weak Taylor Approximation Scheme for Solutions of JumpDiffusion Stochastic Delay Differential Equations", Abstract and Applied Analysis, vol. 2013, Article ID 750147, 8 pages, 2013. https://doi.org/10.1155/2013/750147
A Robust Weak Taylor Approximation Scheme for Solutions of JumpDiffusion Stochastic Delay Differential Equations
Abstract
Stochastic delay differential equations with jumps have a wide range of applications, particularly, in mathematical finance. Solution of the underlying initial value problems is important for the understanding and control of many phenomena and systems in the real world. In this paper, we construct a robust Taylor approximation scheme and then examine the convergence of the method in a weak sense. A convergence theorem for the scheme is established and proved. Our analysis and numerical examples show that the proposed scheme of high order is effective and efficient for Monte Carlo simulations for jumpdiffusion stochastic delay differential equations.
1. Introduction
Stochastic delay differential equation models play a very important role in the study of many fields such as economics and finance, chemistry, biology, microelectronics, and control theories. Models of this type can more accurately describe many phenomena in the real world by taking into account the effect of time delay and random noise. By time delay, it means a certain period of time is required for the effect of an action to be observed after the moment when the action takes place. This phenomenon exists in most systems in almost any area of science; for example, a patient shows symptoms of an illness days or even weeks after he/she was infected. Similarly, random noises appear in almost all real world phenomena and systems, for example, the motion of molecules and the price of assets in financial markets. Hence, study of SDDEs is an important undertaking in order to understand real world phenomena and systems precisely.
Over the last couple of decades, a lot of work has been carried out to study differential equations with delay and/or random noises, for example, [1–5]. The bestknown and wellstudied theory and systems include the delay differential equations (DDEs) presented by Kolmanvskii and Myshkis [6] and their stochastic generalizations and the stochastic delay differential equations (SDDEs) established by Mohammed [7, 8], Mao [9, 10], and Mohammed and Scheutzow [11]. Other SDDEs theories of interest include, for instance, the socalled SDDEs with Markovian switching and Poisson jumps. These models have been investigated in the literature [12–14].
Analytical solutions of SDDEs can hardly be obtained. It is thus important to develop and study discretetime approximation methods for solving SDDEs. Discretetime approximations may be divided into two categories: weak approximations and strong approximations [15]. Some implicit and explicit numerical approximation methods for SDDEs in strong approximation sense were derived by Küchler and Platen [16]. Weak numerical methods for SDDEs have been studied by Küchler and Platen. The Monte Carlo simulation method has also been developed as a powerful simulation method for SDEs, while weak numerical approximations are required for Monte Carlo simulation [9, 17–20].
In this paper, we extend a fully weak approximation method for SDDEs to the type of jumpdiffusion SDDEs: subject to the initial condition where , is the time delay which is assumed to be constant at all time, is an adapted dimensional Wiener process, and denotes the Poisson random measure. Also here we denote by the almost sure lefthand limit of . The coefficient and are dimensional vectors of Borel measurable functions. Further, defined on is a matrix of Borel measurable function.
The main contributions of this work include construction of a numerical approximation method for weak solutions of SDDEs with jumpdiffusion, and establishment of a theoretical result for the convergence of the scheme.
The remaining part of this paper is organised as follows. In the following section, we give the conditions for the existence and uniqueness of the solution of the jumpdiffusion SDDEs (1) which is applicable to any cases where solution exists, and present various lemmas to be used later for the proof of the convergence theorem. We then introduce, in Section 3, a general weak approximation scheme, where the simplified stochastic Taylor approximation scheme with order is constructed; followed by a convergence theorem and its proof. In Section 4, we give a numerical example to demonstrate the application and the convergence of the numerical scheme.
2. Preliminaries
In this section, we present some basic concepts and definitions to be used in later sections, then establish the conditions for the existence and uniqueness of solution to the Jumpdiffusion SDDEs (1), and then give some lemmas to be used later for the proof of the convergence theorem. In this work, we assume that the dimensional vector valued function for the initial condition, , is right continuous and has lefthand limits.
From (1), we have the following integral form of the jumpdiffusion equation SDDE: where , for , are a sequence of pairs of jump times and corresponding values generated by the Poisson random measure .
Definition 1. Given a filtered probability space , a stochastic process given by is known as a solution of (1) subject to the initial condition (2) if is adapted, the integrals in the equation are well defined, and equalities (3) and (2) hold almost surely. Moreover, if any two solution processes are indistinguishable on with the same initial segment and the same path on , and where is the Euclidean norm, then if (1) has a solution, it is a unique solution for this initial value problem.
To guarantee the existence of a unique solution of the jumpdiffusion SDDE (1), we assume that the coefficients of (1) satisfy the following Lipschitz conditions: for and , as well as the growth conditions for , .
We denote by the Banach space of all dimensional continuous functions on equipped with the supremum norm . Furthermore, we suppose that the set of valued continuous process is measurable with
Following the work in [7, 10], the following theorem can be established for the existence and uniqueness of a solution to the problem defined by (1) and (2).
Theorem 2. Suppose that the Lipschitz conditions and the growth conditions are satisfied, and the initial condition is in . Then (1) subject to the initial condition (2) admits a unique solution.
Now we present some lemmas to be used later for the proof of the convergence theorem. Consider a right continuous process . is called a discretetime numerical approximation with maximum step size , if it is obtained by using a time discretization , and the random variable is measurable for . Further, can be expressed as a function of and the discretetime .
Because of dealing with the approximation of solutions of jumpdiffusion SDDEs, we introduce a concept of weak order convergence due to Kloeden and Platen [15].
Definition 3. A discretetime approximation converges weakly towards at time with order if for each there is a constant , independent of , such that where denotes the set of all polynomials .
We now give some auxiliary results to prepare for the proof of the weak convergence theorem to be presented.
Lemma 4. For and , we have for and .
The proof of the lemma for the case with no delay was established by Platen and BrutiLiberati [21], and a similar procedure can be used for the proof of this lemma.
Lemma 5. Given , there is a bounded constant satisfying for and .
The proof of (10) can be obtained by following that of a lemma for SDEs with jumps but with no delay in [15]. The following results are similar to what was given in Mikulevičius and Platen [22].
Lemma 6. Given , there is a finite constant satisfying for every .
Lemma 7. Given , there is and a bounded constant satisfying for each , , , and , where for all and .
The proof of estimate (12) can be established by following Itô’s formula for SDEs with jumps but with no delay as in [15].
Lemma 8. For , there exist and a finite constant satisfying for each , , and .
3. The JumpAdapted Weak Taylor Approximation Scheme
In Monte Carlo simulations for functionals of jumpdiffusion SDDEs, one uses numerical approximations evaluated only at discretization time. Here, we first give a jump adapted weak approximation Taylor scheme of order , and then study the basic properties of the discrete Taylor approximation in a weak order sense.
First we define the jumpadapted time discretization. Let . The jump adapted time discretization used throughout this paper is with the maximum step size . We choose the time discretization in such a way that all jump times are at the nodes of the time discretization. If the discretization node is not a jump time, then is measurable. Otherwise, is measurable. Also, throughout the paper, we denote the set of all multiindices by where the element is called a multiindex of length and has zero length. In the following, by a component of a multiindex we refer to the integration with respect to the th Wiener process in a multiple stochastic integral. A component with corresponds to integration with respect to time .
Now define the following operators for the coefficient functions:
A subset is the hierarchical set, and its corresponding remainder set is defined by . For each , we can then define the hierarchical set . The weak Taylor method of order is then constructed as follows: where
The multiple stochastic integral is then defined recursively as follows: where is obtained from by deleting its last component, while is obtained from by deleting its first component.
Now, we give the weak convergence theorem of the Taylor approximation with order .
Theorem 9. Given , let be the results obtained from the Taylor scheme (17) corresponding to with maximum step size . Suppose that for , and converges to weakly with order . Assume that the coefficients, , , , are in the space , for and , and the coefficient functions , with , satisfy the growth condition , with , for all , , and . Then for any there is a positive constant , which does not depend on , such that
Proof. For and , consider the Itô process below:
Then we can get .
Define also the process , by
for , , and for .
By the definition of the functional and the terminal condition of the stochastic process , we have
Since converges towards weakly with order , one has
By Lemma 4, we can write
From the properties of stochastic integrals, we obtain
Therefore, we have
where
In the following, we proceed to estimate and in Steps 1 and 2, respectively, and then complete the proof in Step 3.
Step 1. Let us assume that is so smooth that the deterministic Taylor expansion may be applied. Hence, by expanding in , we get
where the remainder term is
where is a diagonal matrix with for and and , respectively.
Therefore, according to the properties of expectation and absolute value, we get
By (31), the Hölder inequality, and Lemma 7, we get
Similarly, by Lemma 5 and the CauchySchwarz inequality, we have
Now, from the CauchySchwarz inequality, Lemmas 8 and 6, and inequalities (33) and (34), we obtain
Step 2. Now we estimate the term in inequality (27). By the jump coefficient and the smooth function , applying the Taylor expansion yields
Similarly, we can estimate the reminders as follows:
Then, by applying the Hölder inequality, Lemmas 8 and 6, inequalities (37) to estimate the inequality above, we get
Step 3. Finally, by inequalities (27) and (35) as well as (38), we have
4. A Numerical Example
Here we give an illustrative example to demonstrate the application and the convergence of the proposed numerical scheme. We consider the following linear SDDE with Poisson jumps: where , , and are, respectively, the drift coefficient, the diffusion coefficient and the jump coefficient; is the jump intensity; and and are the delay coefficients.
By the method for solving linear stochastic differential equations in [23], the analytical solution for is derived as follows: where
According to the weak Taylor approximation scheme (17)–(19) proposed in Section 3, we now expand it with weak order 1 (well known as Euler scheme): Here we have used the jump adapted time discretization, and is the maximum step size.
For higher accuracy and efficiency, one needs to construct higher order numerical schemes. We now give a Taylor scheme of weak order two below:
Next, we study the convergence of the two numerical schemes presented above by using the weak errors measured by and compare the results obtained from these two schemes to the explicit exact solution. We estimate the weak errors by running a very large number of simulations. The exact number depends on the implemented scheme. We use the following parameters: , , , , , and .
In Figure 1, we give the sample paths under the two approximation schemes and the numerical explicit solution of (40). We can see from the figure that the weak Taylor scheme path is closer to the analytical solution line than the Euler scheme.
Now we present the numerical errors generated by the two numerical schemes presented above. From Table 1, we notice that, for all the step sizes used in the numerical experiments, the weak Taylor method is more accurate. Moreover, the errors of the weak order two Taylor method decrease faster than the Euler scheme.

5. Conclusions
In this work, we have extended previous research on weak convergence to a more general class of stochastic differential equations involving both jumps and time delay. We proved that under the Poisson random measure and a fixed time delay, a simplified Taylor method gives weak convergence rate arbitrarily close to order . There is much scope for further work in the context of weak solution of jumpdiffusion SDDEs. For example, it is clearly of great importance to extend the weak convergence theory to the case where coefficients in the equations are not globally Lipschitz, and to develop and analyze new methods that maintain good properties of convergence and stability.
Acknowledgments
The authors thank the anonymous referees and the editors for their constructive comments and suggestions. Yanli Zhou acknowledges financial support from Curtin International Postgraduate Research Scholarship (CIPRS) and Chinese Scholarship Council (CSC). Acknowledgement is also made to the China National Social Science Foundation (SSF) for the research support (10BJY104).
References
 P. Hu and C. Huang, “Stability of stochastic $\theta $methods for stochastic delay integrodifferential equations,” International Journal of Computer Mathematics, vol. 88, no. 7, pp. 1417–1429, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 J. Tan and H. Wang, “Meansquare stability of the EulerMaruyama method for stochastic differential delay equations with jumps,” International Journal of Computer Mathematics, vol. 88, no. 2, pp. 421–429, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 N. Jacob, Y. Wang, and C. Yuan, “Numerical solutions of stochastic differential delay equations with jumps,” Stochastic Analysis and Applications, vol. 27, no. 4, pp. 825–853, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Zhang, S. Gan, and L. Hu, “The splitstep backward Euler method for linear stochastic delay differential equations,” Journal of Computational and Applied Mathematics, vol. 225, no. 2, pp. 558–568, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. Wang and Y. Chen, “Meansquare stability of semiimplicit Euler method for nonlinear neutral stochastic delay differential equations,” Applied Numerical Mathematics, vol. 61, no. 5, pp. 696–701, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. Kolmanovskii and A. Myshkis, Applied Theory of Fundamental Differential Equations, Kluwer, New York, NY, USA, 1992.
 S. E. A. Mohammed, Stochastic Functional Differential Equations, vol. 99 of Research Notes in Mathematics, Pitman, Boston, Mass, USA, 1984. View at: MathSciNet
 S. E. A. Mohammed, “Lyapunov exponents and stochastic flows of linear and affine hereditary systems,” in Diffusion Processesand Related Problems in Analysis, Volume II, vol. 27 of Progress in Probability, pp. 141–169, Birkhäuser, Boston, Mass, USA, 1992. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 X. Mao, Exponential Stability of Stochastic Differential Equations, vol. 182, Marcel Dekker, New York, NY, USA, 1994. View at: MathSciNet
 X. Mao, Stochastic Differential Equations and Applications, Horwood, New York, NY, USA, 1997.
 S. E. A. Mohammed and M. K. R. Scheutzow, “Lyapunov exponents and stationary solutions for affine stochastic delay equations,” Stochastics and Stochastics Reports, vol. 29, no. 2, pp. 259–283, 1990. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 A. V. Svishchuk and Y. I. Kazmerchuk, “Stability of stochastic Ito equations delays, with jumps and Markov switchings, and their application in finance,” Theory of Probability and Mathematical Statistics, no. 64, pp. 141–151, 2002. View at: Google Scholar
 J. Luo, “Comparison principle and stability of Ito stochastic differential delay equations with Poisson jump and Markovian switching,” Nonlinear Analysis: Theory, Methods and Applications A, vol. 64, no. 2, pp. 253–262, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. Li and Z. Chang, “Convergence of numerical solution to stochastic delay differential equation with Poisson jump and Markovian switching,” Applied Mathematics and Computation, vol. 184, no. 2, pp. 451–463, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. E. Kloeden and E. Platen, Numerical Solution of Stochastic Differential Equations, vol. 23 of Applied Mathematics, 3rd corrected printing, Springer, Berlin, Germany, 1999.
 U. Küchler and E. Platen, “Strong discrete time approximation of stochastic differential equations with time delay,” Mathematics and Computers in Simulation, vol. 54, no. 1–3, pp. 189–205, 2000. View at: Publisher Site  Google Scholar  MathSciNet
 C. T. H. Baker and E. Buckwar, “Introduction to the numerical analysis of stochastic delay differential equations,” Numberical Analysis Report 345, University of Manchester, Department of Mathematics, Manchester, UK, 1999. View at: Google Scholar
 H. Gilsing, “On the stability of the Euler scheme for an affine stochastic differential equations with time delay,” Diskussion Paper 20, SFB 373, Humboldt University, Berlin, Germany, 2001. View at: Google Scholar
 C. Tudor and M. Tudor, “On approximation of solutions for stochastic delay equations,” Studii şi Cercetări Matematice, vol. 39, no. 3, pp. 265–274, 1987. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 C. Tudor, “Approximation of delay stochastic equations with constant retardation by usual Itô equations,” Revue Roumaine de Mathématiques Pures et Appliquées, vol. 34, no. 1, pp. 55–64, 1989. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 E. Platen and N. BrutiLiberati, Numerical Solution of Stochastic Differential Equations with Jumps in Finance, vol. 64 of Stochastic Modelling and Applied Probability, Springer, Berlin, Germany, 2010. View at: Publisher Site  MathSciNet
 R. Mikulevičius and E. Platen, “Time discrete Taylor approximations for Itô processes with jump component,” Mathematische Nachrichten, vol. 138, pp. 93–104, 1988. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. L. Evans, An Introduction to Stochastic Differential Equations, Lecture Notes, University of California, Department of Mathematics, Berkeley, Calif, USA, 2006.
Copyright
Copyright © 2013 Yanli Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.