Impact of Correlated Noises on Additive Dynamical Systems

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 parameterH > 1/2. 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 H varies. Certain behaviors of the probability density functions are observed.


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    ,  ≥ 0, in a probability space (Ω, F, ), with Hurst parameter  ∈ (0, 1), is a continuous-time Gaussian process with mean zero, starting at zero and having the following correlation function: In particular, when  = 1/2 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  = 1/2, this noise is uncorrelated and thus is called white noise [3].Correlated noises appear in the modeling of some geophysical systems [4][5][6].
For systematic discussions about fractional Brownian motions and their stochastic calculus, we refer to [7][8][9][10][11][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][17][18].
In this paper, we consider the following scalar stochastic differential equation (SDE): where the drift (⋅) is a Lipschitz continuous function on ,  > 0 is the noise intensity,    is a fractional Brownian motion with  > 1/2, 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.

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 S() be the Schwartz space of rapidly decreasing smooth functions on  and S  () the space of tempered distributions.And denote by ⟨⋅, ⋅⟩ the dual pairing on S  () × S().For 1/2 <  < 1, define where Then, for ,  ∈ S(), Now we can only prove the linear map Γ  is continuous from S() to  2 ().Since Γ  is not continuous from S() to S() (even not a proper operator in S()), we could not obtain a dual map from S  () to S  () by duality.By using Itô's regularization theorem, we construct a unique S  ()valued random variable  : S  () → S  () such that which extends the map Γ *  in view of (5).
Theorem 2. Let   =  ∘  −1 be the image measure of  induced by the map T.Then, for any  ∈ S(), the distribution of ⟨⋅, ⟩ under   is the same as ⟨⋅, Γ  ⟩ under .In particular, is a fractional Brownian motion with Hurst constant .Moreover, where   () ≡ ⟨, 1 [0,] ⟩ is the standard Brownian motion.
(See proof in [20].) Let {F  ,  ∈  + } and {F   ,  ∈  + } be the filtrations generated by {  } and {   }, respectively.Then, in view of (8), we have , where ( * )() := ().So, the filtrated probability space (S  (), F  , ) is the extension of (S  (), F   ,   ).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.

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 [0, ].We have the following result.Hence, for any small enough  > 0, we have which completes the proof when  → 0. 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 E in the above theorem corresponds to the fractional white noise measure.And, henceforth, we take all expectations E with respect to the fractional white noise measure (i.e., for simplicity, we omit the subscript  mentioned above).

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 Mathematical Problems in Engineering 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 E[ℎ(  )] = ∫ R ℎ()(, ); by integration by parts and  = 0 at  = ±∞, we obtain that is, In the following, we numerically simulate this partial differential equation for two special cases: () =  −  3 and () = , with finite noise intensity (for simplicity we take  = 1).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, (, 0) = (1/ √ 2) − 2 /2 .
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 () =  −  2 , this method also applies smoothly.
3.1.Numerical Simulation: () =  −  3 .We first simulate the dynamical evolutions of the probability density function (, ) for the corresponding stochastic differential equation ( 2) with the double-well drift () =  −  3 , for various values of  > 1/2.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  = 0.5.

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.