Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 678976, 6 pages

http://dx.doi.org/10.1155/2014/678976

Research Article

Impact of Correlated Noises on Additive Dynamical Systems

Chujin Li¹ and Jinqiao Duan1,2

¹School of Mathematics and Statistics, Huazhong University of Science and Technology, Wuhan 430074, China

²Department of Applied Mathematics, Illinois Institute of Technology, Chicago, IL 60616, USA

Received 14 April 2014; Accepted 4 August 2014; Published 3 September 2014

Academic Editor: Wuquan Li

Copyright © 2014 Chujin Li and Jinqiao Duan. 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.

Impact of correlated noises on dynamical systems is investigated by considering Fokker-Planck type equations under the fractional white noise measure, which correspond to stochastic differential equations driven by fractional Brownian motions with the Hurst parameter . Firstly, by constructing the fractional white noise framework, one small noise limit theorem is proved, which provides an estimate for the deviation of random solution orbits from the corresponding deterministic orbits. Secondly, numerical experiments are conducted to examine the probability density evolutions of two special dynamical systems, as the Hurst parameter varies. Certain behaviors of the probability density functions are observed.

1. Introduction

Dynamical systems arising from financial, biological, physical, or geophysical sciences are often subject to random influences. These random influences may be modeled by various stochastic processes, such as Brownian motions, Lévy motions, or fractional Brownian motions. A fractional Brownian motion , , in a probability space , with Hurst parameter , is a continuous-time Gaussian process with mean zero, starting at zero and having the following correlation function:

In particular, when it is just the standard Brownian motion. The time derivative of a fractional Brownian motion, , as a generalized stochastic process, has nonvanishing correlation [1, 2] and it is thus called a correlated noise or colored noise. In the special case of , this noise is uncorrelated and thus is called white noise [3]. Correlated noises appear in the modeling of some geophysical systems [4–6].

For systematic discussions about fractional Brownian motions and their stochastic calculus, we refer to [7–12] and the references therein. Fractional Brownian motions have stationary increments and are Hölder continuous with exponent less than , but they are no longer semimartingales, even no longer Markovian. They possess some other significant properties such as long range dependence and self-similarity which result in wide applications in fields such as hydrology, telecommunications, and mathematical finance. During the last decade or so, several reasonable stochastic integrations with respect to fractional Brownian motions were developed. See, for example, Lin [13], Duncan et al. [14], Decreusefond and Üstunel [15], and the references mentioned therein. Stochastic differential equations (SDEs) driven by fractional Brownian motions also have been attracting more attention recently [1, 10, 16–18].

In this paper, we consider the following scalar stochastic differential equation (SDE):

where the drift is a Lipschitz continuous function on , is the noise intensity, is a fractional Brownian motion with , and the initial state value is assumed to be independent of the natural filtration of . Since this system has a unique solution [17, 19], here we intend to understand some impact of correlated noises on this additive dynamical system as the Hurst parameter varies.

This paper is organized as follows. In Section 2, we set up a fractional white noise analysis framework which makes correlated noises as functionals of standard white noises and prove a small noise limit theorem which implies the stochastic continuity of the system with respect to noise intensity. In Section 3, we show that the probability density function of satisfies a Fokker-Planck type partial differential equation with respect to the fractional white noise measure. Then, we implement numerical experiments to examine the probability density evolutions as the Hurst parameter varies. As to one linear system and one double-well system, certain behaviors of the probability density functions are observed.

2. Analysis Framework and Small Noise Limit

2.1. Analysis Framework

White noise framework is one natural and flexible stochastic analysis thoughtway, and fractional white noise analysis takes correlated noise as functionals of standard white noise. This approach has shown to be very effective in investigating distributions and path properties of stochastic processes. In the following, we describe the fractional white noise analysis framework.

Let be the Schwartz space of rapidly decreasing smooth functions on and the space of tempered distributions. And denote by the dual pairing on . For , define

where is beta function; , , .

Lemma 1. For , let

Then, for ,

that is, is the dual map of .

Now we can only prove the linear map is continuous from to . Since is not continuous from to (even not a proper operator in ), we could not obtain a dual map from to by duality. By using Itô’s regularization theorem, we construct a unique -valued random variable such that

which extends the map in view of (5).

Theorem 2. Let be the image measure of induced by the map T. Then, for any , the distribution of under is the same as under . In particular,

is a fractional Brownian motion with Hurst constant . Moreover,

where is the standard Brownian motion. (See proof in [20].)

Let and be the filtrations generated by and , respectively. Then, in view of (8), we have

(1), for all ;(2)for any , a.s. , where . So, the filtrated probability space is the extension of . Thus the stochastic analysis with respect to measure could be reduced to the standard white noise framework naturally. Therefore, we choose the standard white noise measure as the reference measure rather than , and this treatment is more useful and more convenient for applications. For more details, we refer to [20] and the reference therein.

2.2. Small Noise Limit

Now, we consider the SDE (2) in fractional white noise framework

And to investigate the impact of noise on deterministic dynamical system

which is solvable on any finite time interval . We have the following result.

Theorem 3. The solution of (2) converges in probability to the solution of (10) uniformly on any finite time interval .

Proof. Firstly, we rewrite the equation as

Then, by assuming the Lipschitz condition on with Lipschitz constant , it follows from the Gronwall inequality that

Hence, for any small enough , we have

which completes the proof when . In the final step, we have used the self-similarity of the fractional Brownian motion

This theorem provides an estimate for the deviation of random solution orbits from the corresponding deterministic orbits. Note that the expectation in the above theorem corresponds to the fractional white noise measure. And, henceforth, we take all expectations with respect to the fractional white noise measure (i.e., for simplicity, we omit the subscript mentioned above).

3. Probability Density Evolution

For SDE, such as (2), the probability density function of the solution carries significant dynamical information. This is considered here by examining a fractional Fokker-Planck type equation. The key step in the derivation of this Fokker-Planck type equation is the application of Ito’s formula for SDEs driven by fractional Brownian motion, under fractional white noise analysis framework [1, 10, 16, 20, 21]. We sketch the derivation here.

By Ito’s formula [10], Theorem 6.3.6, for a second order differentiable function with compact support, we have

Taking expectations on both sides yields

Let be the probability density function of the solution of the system (2). Recall that ; by integration by parts and at , we obtain

that is,

In the following, we numerically simulate this partial differential equation for two special cases: and , with finite noise intensity (for simplicity we take ). Through these two special cases, we expect to illustrate the impact of correlated noises on additive dynamical systems as the Hurst parameter varies.

Here, we perform the popular Crank-Nicolson scheme in Matlab for (17) with zero boundary values;, the grid size is 0.05, total grid points are 801, and the time step size is 0.01. And the initial probability density function is taken to be standard normal; that is, .

Since the system is tridiagonal, we could solve it using Thomas Algorithm efficiently. Moreover, for other initial conditions and other drift coefficients, for instance, the initial uniform distribution or , this method also applies smoothly.

3.1. Numerical Simulation:

We first simulate the dynamical evolutions of the probability density function for the corresponding stochastic differential equation (2) with the double-well drift , for various values of . The double-well dynamics is a rich and typical model for understanding numerous physical or geophysical systems [22, 23], focusing on the maxima (minima), symmetry, kurtosis, and so forth.

As observed in Figure 1, the probability density function evolves from the unimodal (one peak) to the flat top and then to the bimodal (two peaks) shape for various Hurst parameter values , as time increases. Simultaneously, the effect of Hurst parameter on the dynamics is significant. As value increases, the plateau for becomes lower when time exceeds .

(a)

(b)

(c)

(d)

Figure 1: Plot of with , at .

3.2. Numerical Simulation:

Now, for comparison we investigate the dynamical evolutions of the probability density function of the corresponding stochastic differential equation (2) with the linear drift , which is a rich toy example for understanding dynamical systems.

Also as observed in Figure 2, at given time instants, ’s peak becomes higher as increases. This illustrates the significant and distinguishing influence of Hurst parameter on the dynamics when time evolves. The bigger makes the solution of (2) has more centralized value, but the long time effect shows that the values of the solution distribute more scatteredly.

(a)

(b)

(c)

(d)

Figure 2: Plot of with : , , , and .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

Chujin Li acknowledges the support of Self-Renovation Project HUST2013QN171.

References

C. Bender, “An Itô formula for generalized functionals of a fractional Brownian motion with arbitrary Hurst parameter,” Stochastic Processes and their Applications, vol. 104, no. 1, pp. 81–106, 2003.

J. Duan, C. Li, and X. Wang, “Modeling colored noise by fractional Brownian motion,” Interdisciplinary Mathematical Sciences, vol. 8, pp. 119–130, 2009.

L. Arnold, Stochastic Differential Equations, John Wiley & Sons, New York, NY, USA, 1974.

A. Du and J. Duan, “A stochastic approach for parameterizing unresolved scales in a system with memory,” Journal of Algorithms & Computational Technology, vol. 3, no. 3, pp. 393–405, 2009.

J. Duan, “Stochastic modeling of unresolved scales in complex systems,” Frontiers of Mathematics in China, vol. 4, no. 3, pp. 425–436, 2009.

B. Chen and J. Duan, “Stochastic quantification of missing mechanisms in dynamical systems,” Interdisciplinary Mathematical Sciences, vol. 8, pp. 67–76, 2009.

D. Nualart, “Stochastic calculus with respect to the fractional Brownian motion and applications,” Contemporary Mathematics, vol. 336, pp. 3–39, 2003.

E. Alòs and D. Nualart, “Stochastic integration with respect to the fractional Brownian motion,” Stochastics and Stochastics Reports, vol. 75, no. 3, pp. 129–152, 2003.

A. Chronopoulou and F. Viens, “Hurst index estimation for self-similar processes with long-memory,” Interdisciplinary Mathematical Sciences, vol. 8, pp. 91–118, 2009.

F. Biagini, Y. Hu, B. Øksendal, and T. Zhang, Stochastic Calculus for Fractional Brownian Motion and Applications, Springer, New York, NY, USA, 2008.

Y. S. Mishura, Stochastic Calculus for Fractional Brownian Motion and Related Processes, Springer, New York, NY, USA, 2008.

C. A. Tudor and F. G. Viens, “Variations and estimators for self-similarity parameters via Malliavin calculus,” The Annals of Probability, vol. 37, no. 6, pp. 2093–2134, 2009.

S. J. Lin, “Stochastic analysis of fractional Brownian motions,” Stochastics and Stochastics Reports, vol. 55, no. 1-2, pp. 121–140, 1995.

T. E. Duncan, Y. Hu, and B. Pasik-Duncan, “Stochastic calculus for fractional Brownian motion. I. Theory,” SIAM Journal on Control and Optimization, vol. 38, no. 2, pp. 582–612, 2000.

L. Decreusefond and A. S. Üstunel, “Stochastic analysis of the fractional Brownian motion,” Potential Analysis, vol. 10, no. 2, pp. 177–214, 1999.

F. Baudoin and L. Coutin, “Operators associated with a stochastic differential equation driven by fractional Brownian motions,” Stochastic Processes and Their Applications, vol. 117, no. 5, pp. 550–574, 2007.

M. A. Diop and Y. Ouknine, “A linear stochastic differential equation driven by a fractional Brownian motion with Hurst parameter ,” Statistics & Probability Letters, vol. 81, no. 8, pp. 1013–1020, 2011.

B. Saussereau, “A stability result for stochastic differential equations driven by fractional Brownian motions,” International Journal of Stochastic Analysis, vol. 2012, Article ID 281474, 13 pages, 2012.

L. Denis, M. Erraoui, and Y. Ouknine, “Existence and uniqueness for solutions of one dimensional SDE’s driven by an additive fractional noise,” Stochastics and Stochastics Reports, vol. 76, no. 5, pp. 409–427, 2004.

Z. Huang and C. Li, “On fractional stable processes and sheets: white noise approach,” Journal of Mathematical Analysis and Applications, vol. 325, no. 1, pp. 624–635, 2007.

G. Ünal, “Fokker-Planck-Kolmogorov equation for fBM: derivation and analytical solution,” in Proceedings of the 12th Regional Conference, pp. 53–60, Islamabad, Pakistan, 2006.

D. Farrelly and J. E. Howard, “Double-well dynamics of two ions in the Paul and Penning traps,” Physical Review A, vol. 49, no. 2, pp. 1494–1497, 1994.

E. Kierig, U. Schnorrberger, A. Schietinger, J. Tomkovic, and M. K. Oberthaler, “Single-particle tunneling in strongly driven double-well potentials,” Physical Review Letters, vol. 100, no. 19, Article ID 190405, 2008.