#### 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.