Numerical Approximation of Fractal Dimension of Gaussian Stochastic Processes

In this paper we propose a numerical method to estimate the fractal dimension of stationary Gaussian stochastic processes using the random Euler numerical scheme and based on an analytical formulation of the fractal dimension for filtered stochastic signals. The discretization of continuous time processes through this random scheme allows us to find, numerically, the expected value, variance and correlation functions at any point of time. This alternative method for estimating the fractal dimension is easy to implement and requires no sophisticated routines. We use simulated data sets for stationary processes of the type Random Ornstein Uhlenbeck to graphically illustrate the results and compare them with those obtained whit the box counting theorem.


Introduction
A common practice when trying to measure irregular shapes, such as the perimeter of an island or the length of a coastline, is to use Euclidean geometry, ignoring the fact that these shapes do not correspond to the ones of ideal objects, such as polygons and circles, whose dimension is an integer.In contrast, many forms observed in nature have a special feature that makes their resulting measurement depend on the measuring scale, such that the lower the scale, the higher the value of the measurement, to the point where it becomes infinite or indeterminate.The dimension of irregular shapes is a non-integer dimension known as the Hausdorff-Besicovitch dimension.In the early twentieth century, Lewis Fry Richardson found a relationship that allows determining the value of a constant that indicates the degree of roughness of a coast or a geographical border, and used it to calculate the border length between several countries [1].Based on Richardson's findings, Benoit Mandelbrot [2] found that such irregular shapes have the property of statistical self-similarity, which states that any portion of an irregular shape can be considered a fully scaled image of the entire shape.Later, Mandelbrot coined this property as the fractal dimension to describe sets in d  .There are many studies in the literature in which maps, border areas and natural and physical phenomena are modeled and simulated with fractal features.In [3] for example, a proposal to simulate real maps using the fractal dimension of Gaussian stochastic processes is presented and some reports modeling terrestrial surfaces, irregular oceanic coasts and rivers, can be found in [4]- [9].
There are several approaches to estimating of the fractal dimension.In fact in [10] different types of fractals that appear in scientific research are examined, application in nonlinear dynamical systems is discussed, and an exahustive review of numerical methods for their estimation is done.[11] proposed another method for estimating the fractal dimension of subset of d  , as well as an additional correction to address the problem of resolution, allowing obtaining not only the fractal dimension, but also a confidence interval for this estimation, thus improving the box counting method.They also show some experimental results based on simulated data and real data from the coastline of Norway.The techniques for the measurement of fractal dimension can be divided into two groups: those based on vectors, which include structured algorithms, and those based on matrices, for the analysis of image systems that include the distance transformation method reported in [12].Moreover, [13] introduces the roughness method in order to determine the fractal dimension of the Sierpinski carpet, triangle and tetrahedron, the standard Cantor set and the Menger sponge.Meanwhile, in the study of stochastic signals filtering by [14] is obtained the analytical formulation of the fractal dimension to investigate the effects of a linear filter and the regularity in a particular mono-fractal stochastic process, based on an explicit, continuous formulation, linking the spectral properties of a given filter and the fractal properties of a given process, knowing the function of probability density and output, departing from the correlation coefficient of the filter.With regard to the arguments established by [14], the aim of our study is to estimate, numerically, the fractal dimension of ergodic stationary Gaussian stochastic processes using the random Euler scheme.
This paper is organized as follows.Section 2 presents the Euler random scheme and builds the numerical solution, the expected value and the crossed expected value for a more general stochastic processe that allows calculating the variance and correlation of the solution of problems of random initial value; the definition of fractal dimension and the box counting theorem are described in Section 3. In Section 4 we present a numericalanalytical formulation for calculating length and fractal dimension for stationary Gaussian stochastic processes.A method to calculate the length and fractal dimension of stationary Gaussian stochastic processes that include processes such as random Ornstein Uhlenbeck is presented in Section 5. Section 6 illustrates graphically, some experimental numerical results associated with the calculation of the fractal dimension of ergodic stationary Gaussian stochastic processes, compared with the results obtained from using the box count theorem.

On the Random Euler Numerical Scheme
Consider the initial random value problem of the form1 , ( [ ] ( ) where 0 X is a random variable and therefore ( ) , f X t t are both unknown stochastic processes defined on the same probability space ( ) , , Ω   .The aim of this section is to use the Euler scheme to construct a numerical solution, the expected value and the crossed expected value to random initial value problem (1).
Note that the discretized solution of Equation ( 1) can be expressed in integral form as , d .
By using the approximation , the numerical scheme ( ) is obtained, where n t X , ( ) are second order random variables (2-r.v.) and In particular, consider a (2-r.v.) 0 X and a stochastic process ( ) b t second order (2-s.p.) m.s.integrable defined on the same probability space ( ) , , Let ( ) a t a continuous deterministic function defined on T .The idea is to solve numerically initial random value problems given by the expression Note that the Euler numerical scheme for problem (3)2 takes the form It is easy to show using the principle of mathematical induction, (see [19] and [20] pp.58-59), that n t X and n m t t X X can be expressed as where ( ) . From Equations ( 5), ( 6) and the results reported in [19], it follows that the expected

Fractal Dimension and the Box Counting Theorem
The fractal dimension is an objective way to compare fractals geometrically.This section introduces the concept of fractal dimension and presents the box counting theorem.
Let ( ) and center at a point x X ∈ .We wish to define an integer,

( )
, N A  , to be the least number of closed balls of radius  needed to cover the set .
A That is, ( ) In other words, surround every point where ( ) then A has fractal dimension D .The proof of this theorem and other important definitions can be found at [21].

Analytical Formulation of the Fractal Dimension for Stochastic Processes
An analytical expression for the fractal dimension of filtered stochastic signals with the purpose of studying the effects of a linear filter and the regularity of a given stochastic mono-fractal process is introduced in [14].They present some arguments using the Euclidean length of a deterministic curve ( ) x t with limited duration T , as starting point to estimate the length and fractal dimension of such stochastic signals.
However, stochastic processes can not be described by deterministic continuous curves and irregularities are present even when examined with higher resolutions.This leads to the paradoxical conclusion that the lengths of such signals are not finite.To overcome the fact that stochastic processes lack of a finite longitude, [2] assumes that the longitude x L of a fractal curve depends on the resolution and the length of the measurement ins- trument, λ . [1]proves that measuring ( ) x L λ as a function of λ allows to determine the fractal dimension i.e.: ( ) where x D is the fractal dimension of the curve and x K is a constant.Thus, the length x L can be written as a function of λ given by ( ) ( ) where ( ) ( ) and where the elementary length is given by ( ) ( ) ( ) This definition of length for continuous curves can be extended to ergodic stationary stochastic processes by replacing the arithmetic mean by the mathematical expectation {} .

( ) f  and ( )
G  be real valued functions of the positive real variable  .Note that only is solved numerically.In particular, for a high frequency sampling it can be assumed that h ∆  to obtain an analytical expression given by [22] ( )

Length and Fractal Dimension for Gaussian Stochastic Processes
Following the results reported in [14], we will make a brief extension that also requires a prior information such as the probability density function of the process and the variance and correlation functions.
Consider a stochastic process { } 0 t t X ≥ with values in  defined on a probability space ( ) , , Ω   .The average length of the process can be estimated by the expression with , ; , 0, . Thus, the expected value of the absolute value of the difference between ( ) x t and ( ) x t takes the form , d d , , x x p x t x t denotes the joint probability density function.
Now suppose that the stochastic process { } 0 t t X ≥ is a stationary Gaussian process.Consequently, ( ) x t and ( ) x t are also Gaussian stationary and therefore the joint probability density can be written as ( The probability density of the increment ( ) ( ) Analogously to [14], the probability density of the increment follows a normal distribution x t and ( ) x t belong to the same process { } 0 and the standard deviation is ( ) .
and thus Using Equation ( 9), the fractal dimension x D of the process { } 0 t t X ≥ can be approximated as Finally, Equation ( 18) becomes

Numerical Approach to Fractal Dimension for Stationary Gaussian Processes
In this section we present the numerical construction of the fractal dimension for stationary Gaussian stochastic processes including random Ornstein Uhlenbeck process, using the results of Section 4, the random Euler numerical scheme and the approximation of correlation function using the calculations for { } t t E X X described in Section 2.

Calculation of Fractal Dimension for the Random Ornstein Uhlenbeck Process
Consider the random initial value problem given by ( ) W t is a Gaussian white noise independent of 0 X for each t lying in the inter- val [ ] 0,T (defined as in [16] p. 121) and Note that this is stationary Gaussian and ergodic, and can be considered as a modification of the random walk in continuous time and analog to the order 1 autoregressive process, AR (1) with discrete time, that broadly describes the velocity of a Brownian mass particle under the influence of friction [23].Under random Euler numerical scheme (5), the Equation ( 21) can be expressed as .

Calculation of Fractal Dimension for a Random Mean Reverting Process
Consider the random initial value problem given by W t is a Gaussian white noise independent of 0 X for each t lying in the interval [ ] 0,T (defined as in [16] p. 121) and σ αµ stationary Gaussian and ergodic.When the initial condition is constant, this process is also known as mean reversion process with additive noise.This kind of process has been widely used in modeling the dynamic behavior of commodity prices and short term interest rates (see [24] and references therein), and the simulating of realistic irregular lattices in [3].Equation ( 24) can be written as Under random Euler numerical scheme, Equation (25) can be expressed as .
( ) It is easy to see that

Numerical Results
This section graphically illustrates some experimental numerical results related to the calculation of the fractal dimension for stationary Gaussian stochastic processes associated with random initial value problems ( 21) and (24).

Conclusions
This paper made a slight extension of the approach proposed in [14] for the numerical estimation of the fractal dimension of ergodic stationary Gaussian stochastic processes associated with the random Ornstein Uhlenbeck process.Using random Euler scheme performed numerical construction of the expectation value, variance and correlation functions at any point of time to initial random value problems general.The computational implementation of this alternative method to estimate the fractal dimension is easy and requires no sophisticated routines.It illustrates graphically the results obtained compared with box counting theorem using simulated data sets for stationary processes type random Ornstein Uhlenbeck.
A clear convergence of the proposed methods is shown in x A ∈ by an open ball of radius 0 >  to provide a cover of A by open sets.Because A is compact this cover possesses a finite subcover, consisting of an integer number, say M ˆ, of open balls.By taking the closure of each ball, thus obtained a cover consisting of M closed balls.Let C denote the set of covers of A by at most M closed balls of radius  .Then C contains at least one element.Let { } : 1, 2,3, , f C M →  defined by ( ) f c = number of balls in the cover c C ∈ .Then ( ) { } : f c c C ∈ is a finite set of positive integers.It follows that it contains a least integer, ( ) , N A  .The intuitive idea behind fractal dimension is that a set A has fractal dimension D if:

Figure 1 (
a) shows a comparison between the approximate average fractal dimension of 1000 simulated trajectories for Ornstein Uhlenbeck process using the Box Counting theorem and the numerical fractal dimension proposed in Equation (23) for different values of h .In Figure1(b) the dynamic behavior of the fractal dimension is shown for different values of the parameter α .

Figure 2 (Figure 1 .
Figure 1.(a) Numerical comparison of fractal dimension for Ornstein Uhlenbeck processes with = 1 T , = 10 α and = 20 σ for different values of h ; (b) Dynamic behavior of the fractal dimension for Ornstein Uhlenbeck processes with = 1 T ,

Figure 1 (Figure 2 (
a) and Figure2(a).In addition, the results show that the fractal dimension of Ornstein Uhlenbeck processes is independent of the parameter µ .This independence is maintained in the mean reverting process for the case where b) reflect the fact that the fractal dimension of this type of process is more sensitive to the parameter α .
then D is called the fractal dimension A .Also used the notation m A H ∈  , where the Euclidean metric is used.Cover m  by closed square boxes of side length ( ) 2π