Research Article | Open Access

Peng Jiang, Xiaofeng Ju, Dan Liu, Shaoqun Fan, "Exact Finite-Difference Schemes for -Dimensional Linear Stochastic Systems with Constant Coefficients", *Journal of Applied Mathematics*, vol. 2013, Article ID 830936, 6 pages, 2013. https://doi.org/10.1155/2013/830936

# Exact Finite-Difference Schemes for -Dimensional Linear Stochastic Systems with Constant Coefficients

**Academic Editor:**Jesus Vigo-Aguiar

#### Abstract

The authors attempt to construct the exact finite-difference schemes for linear stochastic differential equations with constant coefficients. The explicit solutions to Itô and Stratonovich linear stochastic differential equations with constant coefficients are adopted with the view of providing exact finite-difference schemes to solve them. In particular, the authors utilize the exact finite-difference schemes of Stratonovich type linear stochastic differential equations to solve the Kubo oscillator that is widely used in physics. Further, the authors prove that the exact finite-difference schemes can preserve the symplectic structure and first integral of the Kubo oscillator. The authors also use numerical examples to prove the validity of the numerical methods proposed in this paper.

#### 1. Introduction

Certainty barely exists as a natural phenomenon in isolation. On the contrary, the phenomenon of certainty is an integral part of a complex environment full of interrelations and interactions. Physics is the general analysis of nature, including elements of matter, motion, space, and time, all of which are relevant to the concept of certainty and uncertainty. In order to capture a clear understanding of physical movements in the natural world by using modeling techniques, it is inevitable to take the effects, as the result of the uncertainty upon any development course of any entities, into full consideration. The stochastic differential equations are one of the best known techniques to depict physical movement in precision. For instance, statistical physics is based on the ergodicity assumption of the development of a system and stochastic differential equations can give perfect solutions to such scenarios. Moreover, the stochastic differential equations and relevant applications have been applied in such a variety of scientific fields such as stochastic control, stochastic neural network, financial economic studies, species dynamics, and electronic engineering [1–5]. In this regard, it is of crucial importance to explore stochastic differential equations. Major research findings of the stochastic differential equations can be found in [4, 5].

Although Mao [4] and Kloeden and Platen [6] have managed to calculate the expression of the explicit solution to certain types of stochastic differential equations, it is not always possible to derive the explicit solution to equations of such kind. Strong nonlinearity and coupling characteristics of stochastic differential equations are the major reason and, in this sense, it is of significant importance to establish an effective numerical method [6–15] that is able to make rational judgment by using computer simulations. In saying so, a new branch of computational mathematics is emerging as the numerical solutions to stochastic differential equations. The fundamental theories in this area can be found in [6, 15].

In recent years, the exact finite-difference schemes of ordinary differential equations have captured more and more attention from the academia. Such exact finite-difference schemes are the same to the solutions to the original equations, whilst being able to retain the characteristics of the solutions. Further, such exact finite-difference schemes can be utilized to rationalize the construction of nonstandardized finite-difference schemes [16], for example, the exponentially fitted method [17–19]. As is well known, the linear ordinary differential equations with constant coefficients involve an exact finite-difference scheme (see [16, 20, 21]) and some scholars have discussed exact discretizations of other systems, for example, Vigo-Aguiar and Ferràndiz [22], Cieśliński [23, 24], Sakamoto et al. [25], Mickens and Washington [26], Roeger et al. [27–29], and so on. However, there are only a few published papers [30–32] discussing whether an exact finite-difference scheme exists for widely applied linear stochastic differential equations with constant coefficients. In order to bridge the gap, this paper aims to calculate the exact finite-difference schemes of general -dimensional linear stochastic differential equations with constant coefficients and hence apply the results to the Kubo oscillator so as to prove the structure-preserving property of the exact finite-difference schemes, which is a stochastic sample from [24].

The following part of the paper is organized as follows. In Section 2, the authors will discuss the exact finite-difference schemes for Itô and Stratonovich type stochastic differential equations. In Section 3, the authors demonstrate how to apply the findings generalized from Section 2 to the Kubo oscillator. Finally, the authors will use numerical examples to illustrate the validity of the findings.

#### 2. The Exact Finite-Difference Schemes for Linear Stochastic Differential Equations with Constant Coefficients

Letbe a complete probability space with a filtration. The filtrationis increasing and right continuous, andcontains all-null sets. Letbe a standard-dimensional Brownian motion defined on, whose increment() is a Gaussian random variable.

Then we discuss numerical methods for strong solutions to stochastic differential equations: where the deterministic termis the drift coefficient and the stochastic terms() are diffusion coefficients. The solution of (1) can be written as Given an equidistant discretization of the intervalwith grid points() and letting, the stochastic integrals in (2) can be calculated by the limit of the approximating sums in the mean square sense, as, The values of stochastic integrals depend on the choice ofand the Itô integrals when. The corresponding Itô type stochastic differential equations are the equations using the usual notation (1) and the Stratonovich integrals when. The corresponding Stratonovich type stochastic differential equations are denoted by The relationships between these two types of stochastic differential equations are given by

The exact finite-difference schemes of-dimensional linear stochastic differential equations with constant coefficients of Itô type can be discussed as whereand() arematrices and() are independent one-dimensional Brownian motion. Suppose that the matricesare commutative; that is, Letbe the time step and.() denotes the increment of Brownian motion. Then the following theorem can be derived.

Theorem 1. *Equation (6) admits the exact finite-difference schemes given by
**
whereis the exact value of solutionto (6) at.*

*Proof. *Using the theories in [4, Chapter 3], the general solution of (6) can be expressed in terms ofas follows:
Insertingandinto (9), it can be deduced that
which ends the proof.

Consider a-dimensional linear stochastic differential equation with constant coefficient of Stratonovich type driven byindependent one-dimensional Brownian motion: whereand() arematrices. It is assumed that the matricesare commutative; that is, The exact finite-difference schemes can be constructed by using the general solution of (11).

Theorem 2. *Equation (11) admits the exact finite-difference schemes given by
**
whereis the exact value of solutionto (11) at.*

*Proof. *Firstly, we prove that the general solution to (11) is given by
Equation (14) follows (9) by a known relation (5) between Itô and Stratonovich type stochastic differential equations. Taking into account thatand, in particular,, it can be obtained that
This completes the proof of the theorem.

*Remark 3. *Although the existence of the exact solution of (6) and (11) is known, it is impossible to compute the exact value of the solution at a fixed time, because of the randomness of Brownian motion. The author of [33] has used a Matlab program to simulate the exact solution of one-dimensional linear stochastic differential equations. The codes in the program produce a discretized Brownian path() with the help of the computer. The value of Brownian motionatis obtained by. Then it is possible to derive the value of exact solutionat. Such sum must be calculated in every step of the iteration, increasing the amount of computation needed. However, the scheme proposed in this paper can overcome such a problem.

#### 3. An Application to the Kubo Oscillator

Consider the Kubo oscillator [34]: whereandare constants andis a standard one-dimensional Brownian motion. The small circle “” beforedenotes stochastic differential equations of Stratonovich type. Using (14), the exact solution to (16) is illustrated as follows.

Theorem 4. *The exact solution of (16) is given by
**
Further, the exact finite-difference scheme for (16) can be constructed by applying (17). Assumeis the time increment andis independent-distributed Gaussian random variables.and() are the exact discrete values toandwhen. The exact finite-difference scheme of (16) is given by
**
whereand.**Due to
**
thenis conservative along the phase flow of (16). That is,is a first integral of (16), indicating that a phase trajectory of (16) is a circle with the center atand with the radius. It can be proved that scheme (18) can preserve the first integralexactly.*

Theorem 5. *Scheme (18) for solving (16) has the propertyfor any.*

*Proof. *Substitution of (18) intoyields
This completes the proof.

It is obvious that (16) is a stochastic Hamiltonian system withand, and thus the phase flow of (16) preserves the symplectic structurefor all. A good analytical and numerical study of stochastic Hamiltonian systems can be found in [34, 35]. The following theorem elicits whether scheme (18) is symplectic.

Theorem 6. *Scheme (18) for solving (16) preserves the symplectic structure; that is,.*

*Proof. *Differentiating (18), it is known that
From the exterior product, it can be derived as
This completes the proof.

#### 4. Numerical Experiments

A linear stochastic differential equation of Itô type can be written as The exact solution to (23) is given by By (8), the exact finite-difference scheme of (23) is Figure 1 exhibits the exact solution of a sample phase trajectory of (23) simulated by (25).

Next, the exact finite-difference scheme (18) can be applied to solve the Kubo oscillator (16). The coefficients of (16) are chosen as,,,,, and. Figure 2 exhibits the numerical solutions of a sample phase trajectory of (16) simulated by (18). Figure 3 shows that the numerical solutions created by (18) could preserve the first integralof (16).

#### 5. Conclusions

In this paper, the authors extend the exact finite-difference schemes to linear stochastic differential equations with constant coefficients. The exact finite-difference schemes have been calculated for general-dimensional Itô and Stratonovich type stochastic differential equations. By using the exact finite-difference schemes to solve the Kubo oscillator, the authors have proven that the findings illustrated in this paper can preserve the symplectic structure and first integral. Numerical examples demonstrate the validity of the exact finite-difference schemes in this paper.

#### Acknowledgment

This work is supported by the National Science Council of the Republic of China, under Project No. NSC 71031003.

#### References

- L. P. Blenman, R. S. Cantrell, R. E. Fennell et al., “An alternative approach to stochastic calculus for economic and financial models,”
*Journal of Economic Dynamics & Control*, vol. 19, no. 3, pp. 553–568, 1995. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - F. G. Ball and O. D. Lyne, “Optimal vaccination policies for stochastic epidemics among a population of households,”
*Mathematical Biosciences*, vol. 177-178, pp. 333–354, 2002. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. Liu and K. Wang, “Global stability of a nonlinear stochastic predator-prey system with Beddington-DeAngelis functional response,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 16, no. 3, pp. 1114–1121, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - X. Mao,
*Stochastic Differential Equations and Their Applications*, Horwood Publishing Series in Mathematics & Applications, Horwood Publishing, Chichester, UK, 1997. View at: MathSciNet - B. Øksendal,
*Stochastic Differential Equations*, Springer, Berlin, Germany, 2005. - P. E. Kloeden and E. Platen,
*Numerical Solution of Stochastic Differential Equations*, vol. 23 of*Applications of Mathematics*, Springer, Berlin, Germany, 1992. View at: MathSciNet - A. Rathinasamy, B. Yin, and B. Yasodha, “Numerical analysis for stochastic age-dependent population equations with Poisson jump and phase semi-Markovian switching,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 16, no. 1, pp. 350–362, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - P. Wang, “Three-stage stochastic Runge-Kutta methods for stochastic differential equations,”
*Journal of Computational and Applied Mathematics*, vol. 222, no. 2, pp. 324–332, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - D. J. Higham, X. Mao, and A. M. Stuart, “Strong convergence of Euler-type methods for nonlinear stochastic differential equations,”
*SIAM Journal on Numerical Analysis*, vol. 40, no. 3, pp. 1041–1063, 2002. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - I. Shoji, “A note on convergence rate of a linearization method for the discretization of stochastic differential equations,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 16, no. 7, pp. 2667–2671, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - L. Hu and S. Gan, “Convergence and stability of the balanced methods for stochastic differential equations with jumps,”
*International Journal of Computer Mathematics*, vol. 88, no. 10, pp. 2089–2108, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - X. Wang and S. Gan, “Compensated stochastic theta methods for stochastic differential equations with jumps,”
*Applied Numerical Mathematics*, vol. 60, no. 9, pp. 877–887, 2010. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - F. Jiang, Y. Shen, and F. Wu, “A note on order of convergence of numerical method for neutral stochastic functional differential equations,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 17, no. 3, pp. 1194–1200, 2012. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - G. N. Milstein,
*Numerical Integration of Stochastic Differential Equations*, vol. 313 of*Mathematics and Its Applications*, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1995. View at: MathSciNet - R. E. Mickens,
*Nonstandard Finite Difference Models of Differential Equations*, World Scientific Publishing, River Edge, NJ, USA, 1994. View at: MathSciNet - T. E. Simos and J. Vigo-Aguiar, “A dissipative exponentially-fitted method for the numerical solution of the Schrödinger equation and related problems,”
*Computer Physics Communications*, vol. 152, no. 3, pp. 274–294, 2003. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - T. E. Simos and J. Vigo-Aguiar, “An exponentially-fitted high order method for long-term integration of periodic initial-value problems,”
*Computer Physics Communications*, vol. 140, no. 3, pp. 358–365, 2001. View at: Publisher Site | Google Scholar - A. Tocino and J. Vigo-Aguiar, “Symplectic conditions for exponential fitting Runge-Kutta-Nyström methods,”
*Mathematical and Computer Modelling*, vol. 42, no. 7-8, pp. 873–876, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - R. B. Potts, “Differential and difference equations,”
*The American Mathematical Monthly*, vol. 89, no. 6, pp. 402–407, 1982. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - R. P. Agarwal,
*Difference Equations and Inequalities: Theory, Methods, and Applications*, vol. 228 of*Monographs and Textbooks in Pure and Applied Mathematics*, Marcel Dekker, New York, NY, USA, 2nd edition, 2000. View at: Zentralblatt MATH | MathSciNet - J. Vigo-Aguiar and J. M. Ferrándiz, “A general procedure for the adaptation of multistep algorithms to the integration of oscillatory problems,”
*SIAM Journal on Numerical Analysis*, vol. 35, no. 4, pp. 1684–1708, 1998. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J. L. Cieśliński, “Comment on ‘conservative discretizations of the Kepler motion’,”
*Journal of Physics A*, vol. 43, no. 22, Article ID 228001, 4 pages, 2010. View at: Publisher Site | Google Scholar | MathSciNet - J. L. Cieśliński, “On the exact discretization of the classical harmonic oscillator equation,”
*Journal of Difference Equations and Applications*, vol. 17, no. 11, pp. 1673–1694, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - T. Sakamoto, N. Hori, and Y. Ochi, “Exact linearization and discretization of nonlinear systems satisfying a lagrange pde condition,”
*Transactions of the Canadian Society for Mechanical Engineering*, vol. 35, no. 2, pp. 215–228, 2011. View at: Google Scholar - R. E. Mickens and T. M. Washington, “A note on exact finite difference schemes for the differential equations satisfied by the Jacobi cosine and sine functions,”
*Journal of Difference Equations and Applications*, vol. 19, no. 6, pp. 1042–1047, 2013. View at: Publisher Site | Google Scholar - L.-I. W. Roeger and R. E. Mickens, “Exact finite difference and non-standard finite difference schemes for $dx/dt=-\lambda {x}^{\alpha}$,”
*Journal of Difference Equations and Applications*, vol. 18, no. 9, pp. 1511–1517, 2012. View at: Publisher Site | Google Scholar - L.-I. W. Roeger and R. E. Mickens, “Exact finite-difference schemes for first order differential equations having three distinct fixed-points,”
*Journal of Difference Equations and Applications*, vol. 13, no. 12, pp. 1179–1185, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - L.-I. W. Roeger, “Exact finite-difference schemes for two-dimensional linear systems with constant coefficients,”
*Journal of Computational and Applied Mathematics*, vol. 219, no. 1, pp. 102–109, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - T. Simos, “The exact discrete model of a third-order system of linear stochastic differential equations with observable stochastic trends,”
*Macroeconomic Dynamics*, vol. 13, no. 5, pp. 656–672, 2009. View at: Publisher Site | Google Scholar - H. Schurz, “Numerical analysis of stochastic differential equations without tears,” in
*Handbook of Stochastic Analysis and Applications*, V. Lakshmikantham and D. Kannan, Eds., pp. 237–359, Dekker, New York, NY, USA, 2002. View at: Google Scholar | MathSciNet - R. García-Rubio, “Accurate numerical solution of Black-Scholes option pricing equations,”
*International Journal of Financial Markets and Derivatives*, vol. 2, no. 3, pp. 236–243, 2011. View at: Publisher Site | Google Scholar - D. J. Higham, “An algorithmic introduction to numerical simulation of stochastic differential equations,”
*SIAM Review*, vol. 43, no. 3, pp. 525–546, 2001. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - G. N. Milstein, Yu. M. Repin, and M. V. Tretyakov, “Numerical methods for stochastic systems preserving symplectic structure,”
*SIAM Journal on Numerical Analysis*, vol. 40, no. 4, pp. 1583–1604, 2002. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - G. N. Milstein, Yu. M. Repin, and M. V. Tretyakov, “Symplectic integration of Hamiltonian systems with additive noise,”
*SIAM Journal on Numerical Analysis*, vol. 39, no. 6, pp. 2066–2088, 2002. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet

#### Copyright

Copyright © 2013 Peng Jiang 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.