Journal of Applied MathematicsVolume 2012 (2012), Article ID 305945, 14 pagesdoi:10.1155/2012/305945

Research Article

Analytic Approximation of the Solutions of Stochastic Differential Delay Equations with Poisson Jump and Markovian Switching

Hua Yang¹ and Feng Jiang²

¹School of Mathematics and Computer, Wuhan Polytechnic University, Wuhan 430023, China²School of Statistics and Mathematics, Zhongnan University of Economics and Law, Wuhan 430073, China

Received 2 May 2011; Revised 2 July 2011; Accepted 5 July 2011

Academic Editor: Ying U. Hu

Copyright © 2012 Hua Yang and Feng Jiang. 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.

Abstract.

We are concerned with the stochastic differential delay equations with Poisson jump and Markovian switching (SDDEsPJMSs). Most SDDEsPJMSs cannot be solved explicitly as stochastic differential equations. Therefore, numerical solutions have become an important issue in the study of SDDEsPJMSs. The key contribution of this paper is to investigate the strong convergence between the true solutions and the numerical solutions to SDDEsPJMSs when the drift and diffusion coefficients are Taylor approximations.

1. Introduction

Recently there has been an increasing interest in the study of stochastic differential delay equations with Poisson jump and Markovian switching (SDDEsPJMSs). Svīshchuk and Kazmerchuk [1] investigated stability of stochastic differential delay equations with linear Poisson jumps and Markovian switchings. While Luo [2] discussed the comparison principle and several kinds of stability of Itô stochastic differential delay equations with Poisson jump and Markovian switching. Besides, Li and Chang [3] discussed the convergence of the numerical solutions of stochastic differential delay equations with Poisson jump and Markovian switching. In the present paper we will further research this topic and our focus is on the convergence of numerical solution to stochastic differential delay equation with Poisson jump and Markovian switching when the coefficients are Taylor approximations.

Stochastic differential delay equation with Poisson jump and Markovian switching may be considered as extension of stochastic differential delay equation with Poisson jump. Of course, it may also be regarded as an generalization of stochastic differential delay equation with Markovian switching. Similar to stochastic differential delay equations with Poisson jump, explicit solutions can hardly be obtained for the stochastic differential delay equations with Poisson jump and Markovian switching. Thus appropriate numerical approximation schemes such as the Euler (or Euler-Maruyama) are needed if we apply them in practice or to study their properties. There is an extensive literature concerning the approximate schemes for either stochastic differential delay equations with Poisson jump or stochastic differential delay equations with Markovian switching [4–12].

However, the rate of convergence to the true solution by the numerical solution is different for different numerical schemes [13]. Recently, Janković and Ilić [14] have investigated the rate of convergence between the true solution and numerical solution of the stochastic differential equations in the sense of the -norm when the drift and diffusion coefficients are Taylor approximations, up to arbitrary fixed derivatives. Moreover, Jiang et al. [15] generalized the Taylor method to stochastic differential delay equations with Poisson jump. Since the rate of convergence for such a numerical method is faster than the result obtained in [13], in this paper, we intend to generalize the above method to the SDDEsPJMSs case and consider the strong convergence between the true solution and numerical; solution to SDDEsPJMSs if the drift and diffusion coefficients are Taylor approximations, up to arbitrary fixed derivatives. To the best of our knowledge, so far there seem to be no existing results. Therefore, the aim of this paper is to close this gap.

In Section 2, we introduce necessary notations and approximation scheme. Then comes our main result that the Taylor approximate solutions will converge to the true solutions of SDDEsPJMSs. The proof of this main result is rather technical so we present several lemmas in Section 3 and then complete the proof in Section 4.

2. Approximation Scheme and Hypotheses

Throughout this paper, we let be a complete probability space with a filtration satisfying the usual conditions (i.e., it is increasing and right-continuous while contains all -null sets). Let be the family of continuous function from to with the norm . We also denote by the family of all bounded, -measurable, -valued random variables and denote by the family of all -measurable, -valued random variables satisfying , where is a positive constant.

Let , be a one-dimensional Brownian motion defined on the probability space and let , be a scalar Poisson process with intensity which is independent of . Let , be a right-continuous Markov chain on the probability space taking values in a finite state space with the generator given by where . Here is the transition rate from to if while . We assume that the Markov chain is independent of the Brownian motion and Poisson process . It is well known that almost every sample path of is a right-continuous step function with finite number of simple jumps in any finite subinterval of .

Consider stochastic differential delay equations with Poisson jump and Markovian switching of the form on the time interval   with initial data being independent of and , and , where

In the paper, , and denote, respectively, th-order partial derivatives of and with respect to . Further, to ensure the existence and uniqueness of the solution to (2.2), we impose the following hypotheses.

        (H₁) , and satisfy Lipschitz and linear growth condition; that is, there exists a positive constant such that

for and .

        (H₂) There exist constants and such that, for all and ,

        (H₃) , and have Taylor approximations in the second argument, up to th, th, and th derivatives, respectively.

        (H₄) Partial derivatives of the order , , and of the functions , , and , , and , are uniformly bounded; that is, there exist positive constants , and obeying

For some sufficiently large integer , we define the time step by , where . Then, the approximation solution to (2.2) is computed by on , and for any , where, for with integer ,

In order to show the strong convergence of the numerical solutions and the exact solutions to (2.2), we will need the following assumption.

        (H₅) There exists a positive constant such that for any , where .

Remark 2.1. shows that the exact solutions to (2.2) admit finite moments; see [15]. If , then shows also that the numerical solutions and the exact solutions admit finite moments (see, [1, 3]).We can now state our main result of this paper.

Theorem 2.2. Under assumptions–, for any  , then

The proof of this theorem is rather technical. We will present a number of useful lemmas in Section 3 and then complete the proof in Section 4.

3. Lemmas

Throughout our analysis, denote generic constants, independent of . In order to prove the main theorem, the following lemmas are useful.

Lemma 3.1. If assumptions, and hold, then, for, where is a positive constant independent of .

Proof. For notation simplicity reason, let us denote that Obviously, for any , there exists an integer such that . Then, by (2.7) and (2.8) we obtain Hence, we have By the Hölder inequality, we have By the Burkholder-Davis-Gundy inequality and the Hölder inequality, we yield, for some positive constant , For the jump integral, we convert to the compensated Poisson process , which is a martingale with see, for example, [1, 4, 12, 16]. By the Burkholder-Davis-Gundy inequality and the Hölder inequality, for some positive constant , we then obtain Now, by the mean value theorem there is a such that This, together with –yields that Similarly, we can also show that there are two positive constant and for which Next, since the boundedness of , , and , we then have some positive constant , independent of , such that The desired assertion is complete.

Lemma 3.2. If assumptions hold, then, for , where and is a positive constant independent of .

Proof. Clearly, for any , there exists an integer such that . In what follows, we split the following three cases to complete the proof. Case 1. If , we then have, from , Case 2. If , then, by Lemma 3.1, we have Case 3. If , note that Then, by the above two cases, it follows easily that Now, combining the three cases altogether, for and , where is a positive constant independent of . The proof is complete.

Lemma 3.3. If assumptions and hold, then, for , where , and are positive constants dependent on , but independent of .

Proof. Let be the integer part of . Then with being . By and , we derive Now, by the Markov property [6], we have So, (3.19) is complete. Similarly, we can show (3.20) and (3.21).

4. Proof of Theorem 2.2

Let us now begin to prove our main result Theorem 2.2. Clearly, So, for any , by the Hölder inequality, By the Burkholder-Davis-Gundy inequality, for some positive constant , we have By Doob’s martingale inequality and the compensated Poisson integral, for some positive constant , we yield Now, by virtue of the assumption , we obtain that

Moreover, by and Lemma 3.2, we have and by Lemma 3.1 and , there exists a such that Hence, by (4.6) and (4.7), together with Lemma 3.3, we yield Similarly, we have Now we use to denote a generic positive constant that may change between occurrences. Hence, by (4.2)–(4.4) and (4.8)–(4.10), we yield that

The continuous Gronwall inequality then gives

This proof is therefore complete.

Remark 4.1. If , then the approximate scheme (2.7) is reduced to Euler-Maruyama method for stochastic differential delay equations with Poisson Jump and Markovian switching, which has been discussed in [3].

Remark 4.2. As is well known, Taylor approximation is effectively applicable in engineering if the equations can be solved explicitly. If not, since polynomials are very useful analytic functions, the approximation in the paper can be useful in other applications of stochastic Taylor expansion, especially in the construction of various time discrete approximations of Itô processes by using Itô-Taylor expansion such as Euler-Maruyama approximation and Milstein approximation, which has order 1/2 and 1, respectively [3, 8, 13]. All these show that numerical methods based on Taylor expansions of higher degrees could be improved by combining them with analytic approximations presented in this paper.

Acknowledgments

The project reported here was supported by the Foundation of Wuhan Polytechnic University and Zhongnan University of Economics and Law for Zhenxing.

References

A. V. Svīshchuk and Y. I. Kazmerchuk, “Stability of stochastic Itô equations with delays, with jumps and Markov switchings, and their application in finance,” Theory of Probability and Mathematical Statistics, no. 64, pp. 141–151, 2002.

J. Luo, “Comparison principle and stability of Ito stochastic differential delay equations with Poisson jump and Markovian switching,” Nonlinear Analysis: Theory, Methods & Applications, vol. 64, no. 2, pp. 253–262, 2006.

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.

D. J. Higham and P. E. Kloeden, “Strong convergence rates for backward Euler on a class of nonlinear jump-diffusion problems,” Journal of Computational and Applied Mathematics, vol. 205, no. 2, pp. 949–956, 2007.

J. Bao and Z. Hou, “An analytic approximation of solutions of stochastic differential delay equations with Markovian switching,” Mathematical and Computer Modelling, vol. 50, no. 9-10, pp. 1379–1384, 2009.

X. Mao and C. Yuan, Stochastic Differential Equations with Markovian Switching, Imperial College Press, London, UK, 2006.

X. Mao, C. Yuan, and G. Yin, “Approximations of Euler-Maruyama type for stochastic differential equations with Markovian switching, under non-Lipschitz conditions,” Journal of Computational and Applied Mathematics, vol. 205, no. 2, pp. 936–948, 2007.

Y. Saito and T. Mitsui, “Stability analysis of numerical schemes for stochastic differential equations,” SIAM Journal on Numerical Analysis, vol. 33, no. 6, pp. 2254–2267, 1996.

B. Li and Y. Hou, “Convergence and stability of numerical solutions to SDDEs with Markovian switching,” Applied Mathematics and Computation, vol. 175, no. 2, pp. 1080–1091, 2006.

L.-S. Wang, C. Mei, and H. Xue, “The semi-implicit Euler method for stochastic differential delay equations with jumps,” Applied Mathematics and Computation, vol. 192, no. 2, pp. 567–578, 2007.

C. Yuan and W. Glover, “Approximate solutions of stochastic differential delay equations with Markovian switching,” Journal of Computational and Applied Mathematics, vol. 194, no. 2, pp. 207–226, 2006.

F. Jiang, Y. Shen, and J. Hu, “Stability of the split-step backward Euler scheme for stochastic delay integro-differential equations with Markovian switching,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 2, pp. 814–821, 2011.

P. E. Kloeden and E. Platen, Numerical Solution of Stochastic Differential Equations, vol. 23 of Applications of Mathematics, Springer, Berlin, Germany, 1992.

S. Janković and D. Ilić, “An analytic approximation of solutions of stochastic differential equations,” Computers & Mathematics with Applications, vol. 47, no. 6-7, pp. 903–912, 2004.

F. Jiang, Y. Shen, and L. Liu, “Taylor approximation of the solutions of stochastic differential delay equations with Poisson jump,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 2, pp. 798–804, 2011.

D. J. Higham and P. E. Kloeden, “Numerical methods for nonlinear stochastic differential equations with jumps,” Numerische Mathematik, vol. 101, no. 1, pp. 101–119, 2005.