- Research
- Open Access
- Published:

# Modeling of applied problems by stochastic systems and their analysis using the moment equations

*Advances in Difference Equations*
**volume 2013**, Article number: 152 (2013)

## Abstract

The paper deals with systems of linear differential equations with coefficients depending on the Markov process. Equations for particular density and the moment equations for given systems are derived and used in the investigation of solvability of initial problems and stability. Results are illustrated by examples.

**MSC:**34K50, 60H10, 60H30, 65C30.

## 1 Introduction

Most of the notable achievements in theoretical economics in the last fifty years were related to finances. The first Nobel Memorial Prize in Economic Sciences was awarded in 1969 jointly to Frisch and Tinbergen ‘for having developed and applied dynamic models for the analysis of economic processes’. Many of their works were devoted to the development of mathematical methods to the analysis of economic processes, including mathematical modeling of financial processes [1]. The second Nobel Memorial Prize in Economic Sciences was awarded in 1970 to Samuelson ‘for the scientific work through which he has developed static and dynamic economic theory and actively contributed to raising the level of analysis in economic science’. In his works he studied the role of expectations in the theory of finance.

The growing importance of finance theory in economics is linked to two trends: still wider use of mathematics in the modeling of economic processes, and using the results of theoretical economics in practice. The both trends have a close relationship to finance. Mathematical modeling assumes the exact determination of the parameters - they usually are expressed in the finances; application of the theory in practice assumes description of the cash flows and the risk of using models.

Significant development of the theory of finance, which includes the theory of corporate finance and the theory of investment, occurred in the twentieth century. Until then, the theory of finance was developed as a theory of state finance, in the twentieth century became the theory of capital markets. The amount of significant works on the theory of finance were written in the years 1950 to 2000. Bachelier, the founder of the modern theory of finance, has merit that the theory of finance received a mathematical basics. He anticipated many of the ideas of the twentieth century in his works: the relationship between random and diffusion processes, Markov processes, the theory of Brownian motion and much more than today lies not only in the investment theory. One of the first models of the offers loan funds was built in the early twentieth century by Fisher. Equations to balance between savings and investments (known as IS-LM model, and Mandella-Fleming model) are the basis of the modern macroeconomics. Its authors Hicks and Mundell are Nobel Prize winners. Mundell, in addition, created the theory of optimum currency areas, which allows to call him the father of the euro. In the theory of financial investment, there is no concept that would be such widely verified and so little credible as ‘efficient markets’. The so-called efficient market hypothesis performs a primary function - to justify the use of probabilistic calculation in the analysis of capital markets. But if markets are ‘nonlinear stochastic dynamical systems’, the use of standard statistical analysis can lead to erroneous results, especially if they are based on the model of random walks.

One of the methods that permit to examine the stability of stochastic systems is a traditional method of Lyapunov functions, which was developed, for example, in the works by Barbashin [2], Hasminski [3], Valeev [4], Zubov [5] and others.

Investigating the mean stability or mean square stability of solutions of differential equations with random coefficients depending on Markov process is a current problem. The theory of Markov processes was studied in the works by Chung [6], Davis [7], Dynkin [8, 9], Kolmogorov [10], Lèvy [11], Skorohkod [12] and others. The use of the theory of Markov processes to the study of various economic processes can be found in the works by Elliot, Kopp [13], Malliaris, Brock [14] and Williams [15].

Dynamic systems considered in the present paper belong to the class of the so-called systems with random states. The works by Artem’ev [16], Katz, Krasovskii [17] and others are dedicated to such systems.

We offer a new approach to simulation by creating algorithms for the construction of moment equations and their quantification. The origin of the theory of moment equations and their use in the examination of the stability can be found in the works by Valeev [4] and his scientific school (*e.g.*, [18]).

In the present paper we derive the functional equations for particular density functions and the moment equations for the system which are used in the investigation of solvability and mean square stability. There is shown the application of the results to solve various problems of practice.

## 2 Statement of the problem

Let $(\mathrm{\Omega},\mathcal{F},\mathbb{P})$ be a probability space (see, for example, [19]). On the probability space, we consider the initial value problem formulated for the stochastic system

where *A* is an $m\times m$ matrix with random elements, *B* is an *m*-dimensional column vector function whose elements are random variables, $\phi :\mathrm{\Omega}\to {\mathbb{R}}^{m}$, $\phi \in C(\mathrm{\Omega})$, $\xi (t)$ is a random Markov process with a finite number of states ${\theta}_{k}$, $k=1,2,\dots ,q$, the probabilities of which are

and satisfy the system of linear differential equations

with the transition matrix ${({\pi}_{ks}(t))}_{k,s=1}^{q}$.

**Definition 1** The *m*-dimensional random vector function $x(t)$, the components of which are random variables is called a solution of the initial value problem (1), (2) if $x(t)$ satisfies (1) and initial condition (2) in the meaning of strong solution (defined in [20]) of the initial Cauchy problem.

Our task is to obtain a reliable and simple method for investigating the stability of solutions of this class of systems. To solve this task, we present below the method of moment equations. On a series of examples, we demonstrate that the method is effective and useful.

**Definition 2** Let $x\in {\mathbb{R}}^{m}$ be a continuous random variable depending on a random Markov process $\xi (t)$ with *q* possible states ${\theta}_{k}$, $k=1,2,\dots ,q$. The matrices

where

are called moments of the first or second order of the random variable *x* respectively. The values ${E}^{(k)}(t)$ and ${D}^{(k)}(t)$, $k=1,2,\dots ,q$, are called particular moments of the first or second order respectively.

The ${\mathbb{E}}_{m}$ in Definition 2 denotes an *m*-dimensional Euclidian space, functions ${f}_{k}(t,x)$, $k=1,2,\dots ,q$ are the particular density functions of the random variable *x*.

**Remark 1** The moments of the random variable *x* in a scalar case, $x\in \mathbb{R}$, are defined for any $s=1,2,\dots $ , and are called moments of the *s* th order. The particular moments are defined by the formula

Several different stability statements are possible. We here recall mean square stability definition, which is based on that given in [3].

**Definition 3** The trivial solution of the associated homogenous system to system (1) is said to be mean square stable on the interval $[0,\mathrm{\infty})$ if for each $\epsilon >0$ there exists $\delta >0$ such that any solution $x(t)$ of the associated system, corresponding to the initial data ${x}_{0}$, exists for all $t\ge 0$ and the mathematical expectation

## 3 Moment equations for the linear differential equations

Before the initial value problem (1), (2) formulated in the previous section will be investigated, a simpler problem will be studied. First we derive the moment equations in the scalar case of system (1), that is, if instead of the system there is an equation. In the first part of this section, the linear homogenous differential equation, the coefficient of which depends on a random Markov process, with two states only is considered. In the second part, the moment equations are derived for nonhomogenous linear differential equations with *q* possible states of a random process, on which the coefficients depend.

### 3.1 Homogenous linear differential equations

On the probability space $(\mathrm{\Omega},\mathcal{F},\mathbb{P})$, we consider initial value problem (1), (2) where instead of system (1) there is a stochastic linear homogenous differential equation of the first order of the form

where *a* is a scalar function of a random variable. We suppose that the function *a* depends on the random Markov process $\xi (t)$, which has only two states ${\theta}_{1}$, ${\theta}_{2}$ with probabilities

that satisfy the system of linear differential equations

In the following, we use the denotations

**Theorem 1** *Moment equations of any order* $s=0,1,2,\dots $ *for equation* (5) *are of the form*

*Proof* We divide the time line $[0,\mathrm{\infty})$ into intervals of length *h*. Next we replace the considered system of differential equations (5) by an approximated system of difference equations. If we denote ${t}_{n}=nh$, $h>0$, $n=1,2,\dots $ , and approximate $dx({t}_{n+1})/dt$ with $(x({t}_{n+1})-x({t}_{n}))/h$, then the approximated system to system (5) can be written in the form

or the approximated system to system (6) is of the form

In accordance with the formula for total probability, we obtain relationships for the particular density functions ${f}_{k}({t}_{n},x)$, $k=1,2$, which satisfy the following system of functional equations:

Rename ‘${t}_{n}$’ to ‘*t*’ and suppose that the particular density functions can be expressed in powers of parameter *h* by the Taylor formula. Let functions in (9) be represented as

where *O* is Landau order symbol. Now, using the obtained expressions and comparing the left-hand side to the right-hand side of (9) and assuming $h\to 0$, we get

Similarly, decomposition of the particular density functions in (10) gives the second equation

Finally, multiplying equations (11), (12) by ${x}^{s}$, $s=0,1,2,\dots $ and integrating them by parts from −∞ to ∞, in accordance with Definition 2, a system of linear differential equations with constant coefficients (7) can be obtained. □

Let us note that moment equations (7) can be derived in a different way. If system (8) of difference equations for probabilities is known, then the particular moments of the *s* th order satisfy the following relations:

Particular moments contained in the first equation of (13) can be expressed in powers of parameter *h* by the Taylor formula:

If we put the obtained expressions into the first equation of (13), then under assumption $h\to 0$, we get the first equation of system (7). In the same way, using the second equation of (13), the second equation of system (7) can be constructed.

**Example 1** Let us establish conditions for *s*-mean stability of linear differential equation (5). The characteristic equation for the system of moment equations (7) is written as follows:

Therefore, the conditions of asymptotic stability of solutions of moment equations (7), in accordance with the Hurwitz criterion, are of the following form (assume $s\ne 0$, the case $s=0$ is considered below):

Let us use the denotations

where ${p}_{k}^{0}={lim}_{t\to +\mathrm{\infty}}{p}_{k}(t)$, $k=1,2$ and *a* is mean value of coefficients ${a}_{1}$, ${a}_{2}$. It allows us to derive a simpler form of the above conditions:

The domains of stability for moments of various order are determined by their boundaries as it is shown in Figure 1. Any domain of stability includes the third quadrant where the values of coefficients ${a}_{1}$, ${a}_{2}$ are negative, *i.e.*, ${a}_{1}<0$, ${a}_{2}<0$.

Using moment equations, it is also possible to determine the domain of stability for the deterministic equation

where *a* is independent of a random variable $\xi (t)$. This case corresponds to the moment equations of the zeroth order, *i.e.*, if $s=0$.

### 3.2 Nonhomogenous linear differential equation

We have derived the system of moment equations for a linear homogenous equation with random coefficient under assumptions that the random variable can only be in two states. It was a simple enough case that allowed us to understand the process of deriving the system of moment equations. Now we establish a system of moment equations in the same way for linear the nonhomogeneous differential equation

where $\xi (t)$ is the Markov process which has *q* possible states ${\theta}_{1},{\theta}_{2},\dots ,{\theta}_{q}$, with probabilities ${p}_{k}(t)=P\{\xi (t)={\theta}_{k}\}$, $k=1,2,\dots ,q$. We suppose that the probabilities satisfy the system of linear differential equations

where the transition matrix ${({\pi}_{ks}(t))}_{k,s=1}^{q}$ satisfies the following relationships:

Since the coefficients of studied system (14) depend on *t*, we can denote

**Theorem 2** *Moment equations of any order* $s=1,2,\dots $ *for equation* (14) *are of the form*

*Proof* By dividing the time line into intervals of length *h*, we obtain the approximated system

to the considered system (14) and

to system (15).

Particular probability density functions ${f}_{k}({t}_{n},x)$ satisfy, in this case, the system of difference equations

Similarly as in the proof of Theorem 1, we assume that the particular density functions can be represented in powers of parameter *h* by the Taylor formula, and by the same way as in the proof of Theorem 1, we get

The system of moment equations (16) can be derived from the last system for particular probability density functions by using the same modifications as in the proof of Theorem 1. □

## 4 Moment equations for the linear differential system

Now we come back to the initial problem (1), (2) that we have formulated in Section 2. We also suppose that the matrix *A* and vector *B* depend on a random Markov process $\xi (t)$ with *q* possible states, the probabilities of which (3) satisfy the system of linear differential equations (4).

Moreover, we use the denotations

**Theorem 3** *Moment equations of the first and second order respectively for system* (1) *are of the form*

*Proof* The philosophy of the proof is the same as in the proof of Theorem 1, only the calculations are more complicated, because now we work with the matrix case. In a similar way, by dividing the time line into intervals of length *h*, for the particular density functions ${f}_{k}(t,x)$, $k=1,2,\dots ,q$, we get the system of equations

where

Assume that the particular density functions can be expressed in powers of parameter *h* by the Taylor formula. If we put ${t}_{n}=t$, then decompositions of the functions on the left-hand side and on the right-hand side in (20) are equal to

where

Using obtained expressions, next comparing the left-hand side to the right-hand side of equation (20) and assuming $h\to 0$, we get the system of differential equations for the particular density functions

Finally, multiplying equation (21) by *x* and integrating it by parts on the Euclidean space ${\mathbb{E}}_{m}$, in accordance with Definition 2, we obtain a system of linear equations for the particular moments of the first order in the form (18). The particular moments of the second order satisfy the matrix system of differential equations (19) which we get in the same way. The difference is that (21) is multiplied by the matrix $x{x}^{\ast}$, next it is integrated over the Euclidean space ${\mathbb{E}}_{m}$. □

**Remark 2** The moment equations (18), (19) are deterministic and can be solved by using usual methods, *e.g.*, [21].

The following examples illustrate the use of moment equations for the investigation of stability.

**Example 2** Let us investigate the mean square stability of solutions of the homogenous linear differential equation

where the Markov process $\xi (t)$ can be in two states ${\theta}_{1}$, ${\theta}_{2}$ with probabilities

which satisfies the system of differential equations

where $\lambda \ge 0$. We establish the system of moment equations for the considered system. Let the values of the matrix $A(\xi (t))$ corresponding with the states of the Markov process be

Then the system of moment equations of the second order (19) for equation (22) is of the form

Denote

and rewrite the last system of moment equations into the scalar form

The obtained system of moment equations is a system of ordinary linear differential equations. Its stability is determined by the eigenvalues of the matrix system. The characteristic equation

can be transformed into the following equation:

the roots of which are

It is easy to see that the real parts of all eigenvalues are negative or equal to zero. Therefore, the solutions of the system of equations (22) are stable in the mean square.

## References

- 1.
Tinbergen J:

*Income Distribution: Analysis and Policies*. American Elsevier, New York; 1975. - 2.
Barbashyn EA:

*Lyapunov Functions*. Nauka, Moscow; 1970. (in Russian) - 3.
Hasminski RZ Monographs and Textbooks on Mechanics of Solids and Fluids: Mechanics and Analysis 7. In

*Stochastic Stability of Differential Equations*. Sijthoff & Noordhoff, Alphen aan den Rijn; 1980. - 4.
Valeev KG, Dzhalladova IA:

*Optimization of Random Process*. KNEU, Kiev; 2006. (in Russian) - 5.
Zubov VI:

*Methods of A. M. Lyapunov and Their Application*. P. Noordhoff, Groningen; 1964. - 6.
Chung KL:

*Lectures from Markov Processes to Brownian Motion*. Springer, Berlin; 1982. (in Russian) - 7.
Davis MNA:

*Markov Models and Optimization*. Chapman & Hall, London; 1993. - 8.
Dynkin EB: Markov Processes, vol. I. Translated with the authorization and assistance of the author by J. Fabius, V. Greenberg, A. Maitra, G. Majone. Die Grundlehren der Mathematischen Wissenschaften Vol. 121. Academic Press, New York; 1965.

- 9.
Dynkin, EB: Markov Processes, vol. II. Translated with the authorization and assistance of the author by J. Fabius, V. Greenberg, A. Maitra, G. Majone. Die Grundlehren der Mathematischen Wissenschaften, vol. 122. Academic Press, New York (1965)

- 10.
Kolmogorov AN:

*Selected Works of A.N. Kolmogorov (Probability Theory and Mathematical Statistics)*. Springer, New York; 1992. - 11.
Lèvy P:

*Stochastic Processes and Brownian Motion*. Nauka, Moscow; 1972. (in Russian) - 12.
Skorokhod AV Translations of Mathematical Monographs 78. In

*Asymptotic Methods in the Theory of Stochastic Differential Equations*. Am. Math. Soc., Providence; 1989. Translated from the Russian by H. H. McFaden - 13.
Elliot RJ, Kopp PE:

*Mathematics of Financial Markets*. Springer, Berlin; 1999. - 14.
Malliaris AG, Brock WA:

*Stochastic Methods in Economics and Finance*. North-Holland, Amsterdam; 1982. - 15.
Williams D:

*Diffusions, Markov Processes and Martingales*. Wiley, New York; 1979. - 16.
Artem’ev VM, Kazakov IE:

*Handbook on the Theory of Automatic Control*. Nauka, Moscow; 1987. (in Russian) - 17.
Katz IY, Krasovskii NN: On stability of systems with random parameters.

*Prikl. Mat. Meh.*1960, 24(5):809–823. (in Russian) - 18.
Dzhalladova IA:

*Optimization of Stochastic System*. KNEU, Kiev; 2005. (in Russian) - 19.
Jacod J, Shiryaev AN:

*Limit Theorems for Stochastic Processes*. Springer, New York; 1987. - 20.
Gihman II, Skorohod AV:

*Controlled Stochastic Processes*. Springer, Berlin; 1979. - 21.
Samoilenko AM, Perestyuk NA:

*Impulsive Differential Equations*. Vyshcha Shkola, Kiev; 1987. (in Russian)

## Acknowledgements

The first author was supported by the Grant No. P201/11/0768 of the Grant Agency of the Czech Republic (GA CR). The third and fourth authors were supported by the Grant No. 081ŽU-4/2011 of the Grant Agency of the Slovak Republic (KEGA).

## Author information

### Affiliations

### Corresponding author

## Additional information

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

The authors have made the same contribution. All authors read and approved the final manuscript.

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

**Open Access**
This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (
https://creativecommons.org/licenses/by/2.0
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

### Cite this article

Diblík, J., Dzhalladova, I., Michalková, M. *et al.* Modeling of applied problems by stochastic systems and their analysis using the moment equations.
*Adv Differ Equ* **2013, **152 (2013). https://doi.org/10.1186/1687-1847-2013-152

Received:

Accepted:

Published:

### Keywords

- stochastic systems
- Markov process
- moment equations
- solvability
- stability