Next Article in Journal
Results on Neutral Partial Integrodifferential Equations Using Monch-Krasnosel’Skii Fixed Point Theorem with Nonlocal Conditions
Next Article in Special Issue
Tychonoff Solutions of the Time-Fractional Heat Equation
Previous Article in Journal
A Note on Approximate Controllability of Fractional Semilinear Integrodifferential Control Systems via Resolvent Operators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Generalization of Multifractional Brownian Motion

1
Department of Mathematics, Indian Institute of Technology Ropar, Rupnagar 140001, India
2
Cardiff School of Mathematics, Cardiff University, Senghennydd Road, Cardiff CF24 4AG, UK
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(2), 74; https://doi.org/10.3390/fractalfract6020074
Submission received: 13 December 2021 / Revised: 20 January 2022 / Accepted: 27 January 2022 / Published: 30 January 2022

Abstract

:
In this article, some properties of multifractional Brownian motion (MFBM) are discussed. It is shown that it has persistence of signs long range dependence (LRD) and persistence of magnitudes LRD properties. A generalization called here nth order multifractional Brownian motion (n-MFBM) that allows to take the functional parameter H ( t ) values in the range ( n 1 , n ) is discussed. Two representations of the n-MFBM are given and their relationship with each other is obtained.

1. Introduction

The fractional Brownian motion (FBM) introduced in [1] is a generalization of Brownian motion which exhibit long range dependence (LRD) property with normally distributed marginal distributions. However, many real life time-series display both LRD and heavy tailed increments characteristics. An example of a stochastic process which exhibits both LRD and heavy tailed marginals is the linear fractional stable motion (LFSM). In LFSM, the integral is taken with respect to the α -stable random measure instead of Gaussian measure unlike in case of FBM (see e.g., [2]). The constant value of the Hurst exponent of FBM is a strong limitation for example in simulating artificial mountains [3], detection of cancer tumors in medical images [4]. There are examples where Hurst exponent H values above 1 have been observed as, for example, in Neil River data [5]. To allow the Hurst exponent in the interval ( n 1 , n ) , the nth order FBM (n-FBM) is discussed in [6]. A generalization of FBM is discussed in [7] where the Hurst parameter of FBM is replaced by a stochastic process. This generalization also have wide applications since the pointwise irregularity varies stochastically with time. Multifractional Brownian motion (MFBM) is another generalization of FBM in which the Hurst parameter is a deterministic function of time and was introduced in [8]. In MFBM, the local roughness varies which is measured by the local Hölder exponent and is defined by [8]
B H ( t ) ( t ) = 1 Γ ( H ( t ) + 1 / 2 ) 0 ( t s ) H ( t ) 1 / 2 ( s ) H ( t ) 1 / 2 d B ( s ) + 1 Γ ( H ( t ) + 1 / 2 ) 0 t ( t s ) H ( t ) 1 / 2 d B ( s ) = 1 Γ ( H ( t ) + 1 / 2 ) ( ( t s ) + ) H ( t ) 1 / 2 ( ( s ) + ) H ( t ) 1 / 2 d B ( s ) ,
where B ( t ) is the standard Brownian motion and the integration is taken in the mean square sense and H is a Hölder function of exponent β > 0 , such that,
H : [ 0 , ) [ a , b ] ( 0 , 1 )
and | H ( t ) H ( s ) | c | t s | β for some c > 0 . Further,
u + = u if u 0 , 0 if u < 0 .
Choosing appropriate values for the parameters (see e.g., [9,10]) the following functions are considered for Hölder exponent H ( t ) in literature:
linear H ( t ) = a t + b , t [ 0 , T ] ; logistic H ( t ) = c b 1 + exp ( d ( t t 0 T ) ) + b , t [ 0 , T ] ; periodic H ( t ) = a sin 4 π t T + b , t [ 0 , T ] ; arctan H ( t ) = 1 π arctan ( t ) + 1 2 , t > 0 ; H ( t ) = a t b + a t , t > 0 ; H ( t ) = a sin π 2 t α , 0 < a < 1 , 0 < α < 1 , t > 0 ; H ( t ) = 1 b + t α , 1 < b < , 0 < α < 1 , t > 0 .
The MFBM is also analogously defined with harmonizable representation of FBM as follows (see [11])
B ˜ H ( t ) ( t ) = e i t ξ 1 | ξ | H ( t ) + 1 / 2 d B ^ ( ξ ) ,
where B ^ denotes a complex-valued Gaussian measure. We refer to [12] for the precise definition and for the fact that B ˜ H ( t ) is real valued. Further, it is proved in [10] that two representations given in (1) and (4) neither have equality in distribution up to a multiplicative constant nor equality of the correlations which is claimed in Theorem 1 of [12]. Form Theorem 3.1 of [10], it follows
B H ( t ) ( t ) = cos ( ( H ( t ) + 1 / 2 ) π / 2 ) 2 π R e i t ξ 1 | ξ | H ( t ) + 1 / 2 d B ^ ( ξ ) + sin ( ( H ( t ) + 1 / 2 ) π / 2 ) 2 π R e i t ξ 1 i ξ | ξ | H ( t ) 1 / 2 d B ^ ( ξ ) .
Form (5), it is evident that processes B H ( t ) ( t ) and B ˜ H ( t ) ( t ) have the same laws up to a multiplicative deterministic function if H ( t ) = 1 / 2 . The MFBM is called a standard MFBM if its variance at time 1 is 1. The covariance function of standard MFBM is obtained in [13] and is given by
E ( B ˜ H ( t ) ( t ) B ˜ H ( s ) ( s ) ) = D ( H ( t ) , H ( s ) ) t H ( t ) + H ( s ) + s H ( t ) + H ( s ) | t s | H ( t ) + H ( s ) ,
where
D ( x , y ) = Γ ( 2 x + 1 ) Γ ( 2 y + 1 ) sin ( π x ) sin ( π y ) 2 Γ ( x + y + 1 ) sin ( π ( x + y ) / 2 ) .
Note that for x = y , D(x,y) = 1/2. This natural extension of FBM however results in some sense a loss of properties. Unlike FBM, the MFBM does not have stationary increments and the process is no more self-similar. However, the choice of H as a Hölder continuous function ensures continuity of its sample paths.
The LFSM is further generalized to linear multifractional stable motion which allows the scaling exponent (or functional exponent) to be a deterministic function of time as well as heavy tailedness of marginals, see e.g., [14]. In [15], the LFSM is generalized to tempered fractional stable motion. Moreover, the higher order fractional stable motion introduced in [16], allows the Hurst exponent to vary in the interval ( n 1 , n ) . The higher order fractional stable motion generalizes LFSM and n-FBM. In addition, MFBM is extended to so called generalized multifractional Brownian motion (GMFBM) [17]. Using the idea in [6,16], we generalize the MFBM to an nth order multifractional Brownian motion (n-MFBM) that allows the functional Hurst parameter H ( n 1 , n ) . In this article, we show that MFBM has persistence of signs LRD and persistence of magnitudes LRD properties. Further, two representations of the introduced process n-MFBM are given and their relationship with each other is obtained. The rest of the article is organized as follows. In Section 2, the LRD properties of MFBM are discussed. Section 3, deals with the n-MFBM and its properties like variance, covariance and harmonizable representation. We concludes the paper in Section 4. The last section discusses the possible generalizations of n-MFBM and other related processes.

2. Long Range Dependence (LRD)

A stationary finite variance process X ( t ) is said to exhibit short or long-range dependence (SRD or LRD) properties according as k = 0 γ k converges or diverges, where γ k = Cov ( X ( t ) , X ( t + k ) ) for k 1 , being the autocovariance at lag k [5]. A stationary process X ( t ) has persistence of signs LRD property if LRD holds for the process sign ( X ( t ) ) , where
sign ( x ) = + 1 , x > 0 , 0 , x = 0 , 1 , x < 0 ,
and persistence of magnitudes LRD property if LRD holds for the process | X ( t ) | [18]. The increments of FBM possess persistence of signs LRD, when 1 / 2 < H < 1 and also exhibit persistence of magnitudes LRD when 3 / 4 H < 1 , see e.g., [18]. Note that MFBM is neither a stationary process nor does it have stationary increments and hence the definition of LRD given above is not applicable. However, for a non-stationary process Y ( t ) an equivalent definition is given below.
Definition 1.
Let s > 0 be fixed and t > s . Then the process Y ( t ) is said to have LRD property if
Corr ( Y ( s ) , Y ( t ) ) c ( s ) t d ( s ) , as t ,
where c ( s ) is a constant depending on s and d ( s ) ( 0 , 1 ) . Further, f ( h ) g ( h ) , as h tends to infinity, when there exists 0 < a < b < such that for all sufficiently large h, a f ( h ) g ( h ) b .
In analogy to the definition of [18], we have following definitions for the persistence of signs LRD and persistence of magnitudes LRD.
Definition 2.
(i) A second order process X ( t ) is said to have persistence of signs LRD property if for fixed s and large t, there exist functions c 1 ( s ) and d 1 ( s ) ( 0 , 1 ) :
Corr ( sign X ( s ) , sign X ( t ) ) c 1 ( s ) t d 1 ( s ) as t .
(ii) It has persistence of magnitudes LRD property if for fixed s and large t, there exist functions c 2 ( s ) and d 2 ( s ) ( 0 , 1 ) :
Corr ( | X ( s ) | , | X ( t ) | ) c 2 ( s ) t d 2 ( s ) as t .
Proposition 1.
Let B ˜ H ( t ) ( t ) be a standard MFBM given in (4). Then,
(i)
For all admissible H ( t ) , MFBM has persistence of signs LRD in the sense of Definition 2.
(ii)
Also, it has persistence of magnitudes LRD taken in the sense of Definition 2, when
(a)
s > 0 , H ( t ) + H ( s ) < 1 and H ( t ) < 1 / 2 for large t.
(b)
H ( t ) + H ( s ) > 1 , H ( s ) > 1 / 2 for all s > 0 and large t.
Proof. 
The proof follows on the line of [18]. Let
X Y N 2 ( μ , Σ ) ,
be the bivariate normal vector with mean ( 0 , 0 ) T and covariance matrix
Σ = σ 1 2 ρ σ 1 σ 2 ρ σ 1 σ 2 σ 2 2 .
Then, we have (see e.g., [19,20])
P ( X > 0 , Y > 0 ) = P ( X < 0 , Y < 0 ) = 1 4 + 1 2 π sin 1 ρ ,
and
P ( X > 0 , Y < 0 ) = P ( X < 0 , Y > 0 ) = 1 4 1 2 π sin 1 ρ .
Further,
E ( | X | | Y | ) = 2 π σ 1 σ 2 ( 1 ρ 2 + ρ sin 1 ρ ) , E | X | = 2 π σ 1 , E | Y | = 2 π σ 2 ,
which implies
Cov ( | X | , | Y | ) = 2 π σ 1 σ 2 ( 1 ρ 2 1 + ρ sin 1 ρ ) .
Hence
Corr ( | X | , | Y | ) = 2 ( π 2 ) 1 ρ 2 1 + ρ sin 1 ρ .
In addition, E ( sign X ) = E ( sign Y ) = 0 and Var ( sign X ) = Var ( sign Y ) = 1 , see [18]. Moreover, using (14) and (15)
Cov ( sign X , sign Y ) = E ( sign X sign Y ) E ( sign X ) E ( sign Y ) = E ( sign X sign Y ) = 1 4 + 1 2 π sin 1 ρ + 1 4 + 1 2 π sin 1 ρ 1 4 + 1 2 π sin 1 ρ 1 4 + 1 2 π sin 1 ρ = 2 π sin 1 ρ .
It follows that
(i)
Corr ( sign B ˜ H ( s ) ( s ) , sign B ˜ H ( t ) ( t ) ) = Cov ( sign B ˜ H ( s ) ( s ) , sign B ˜ H ( t ) ( t ) ) ,
follows by the fact that B ˜ H ( t ) ( t ) N ( 0 , t 2 H ( t ) ) for all t 0 . We have by (14) and from the fact that Corr ( B ˜ H ( t ) ( t ) , B ˜ H ( s ) ( s ) ) 0 as t for all s 0 ,
Corr ( sign B ˜ H ( s ) ( s ) , sign B ˜ H ( t ) ( t ) ) = 2 π sin 1 ( ρ ( t , s ) ) 2 π ρ ( t , s ) a s t ,
where ρ ( t , s ) = Corr ( B ˜ H ( t ) ( t ) , B ˜ H ( s ) ( s ) ) . From Proposition 7 in [13], we have
Corr ( sign B ˜ H ( s ) ( s ) , sign B ˜ H ( t ) ( t ) ) t H ( t ) , if H ( t ) + H ( s ) < 1 t H ( s ) 1 , if H ( t ) + H ( s ) > 1 ,
for all fixed s > 0 and as t . Since H ( t ) and 1 H ( s ) both belong to ( 0 , 1 ) , we conclude that MFBM has persistence of signs LRD in the sense of Definition 2, for all values of H(t).
(ii)
In addition,
Corr ( | B ˜ H ( t ) ( t ) | , | B ˜ H ( s ) ( s ) | ) = 2 π 2 ( 1 ρ ( t , s ) 2 1 + ρ ( t , s ) sin 1 ρ ( t , s ) ) ρ ( t , s ) 2 π 2 ,
as t , where ρ ( t , s ) = Corr ( B ˜ H ( t ) ( t ) , B ˜ H ( s ) ( s ) ) . So we have,
C o r r ( | B ˜ H ( t ) ( t ) | , | B ˜ H ( s ) ( s ) | ) t 2 H ( t ) , if H ( t ) + H ( s ) < 1 t 2 H ( s ) 2 , if H ( t ) + H ( s ) > 1 .
If H ( t ) + H ( s ) < 1 and H ( t ) < 1 / 2 for large t and for all s > 0 , then we have persistence of magnitudes LRD property. In case H ( t ) + H ( s ) > 1 , persistence of magnitudes occurs when H ( s ) > 1 / 2 for all s > 0 .
Remark 1.
If for all s, H ( t ) + H ( s ) > 1 for all sufficiently large t, then by (19) the MFBM with functional parameter H has persistence of signs LRD in the sense of Definition 2, with functional LRD exponent α ( s ) = H ( s ) 1 . Similarly by (21), we have persistence of magnitude LRD in sense of Definition 2, with functional LRD exponent α ( s ) = 2 H ( s ) 2 , with H ( s ) < 1 / 2 .
In case H ( t ) + H ( s ) < 1 , MFBM may not have persistence of signs and persistence of magnitudes LRD properties when there exist a sequence { t n } tends to infinity s.t. H ( t n ) + H ( s ) < 1 , for every n and H ( t n ) does not have a limit. For example consider the functional parameter H ( t ) = 1 6 + 1 6 | sin ( π 2 t α ) | , α ( 0 , 1 ) . Here H ( t ) + H ( s ) 2 / 3 for all t and s , but lim n H ( t n ) does not exist for any sequence t n . In this case, there does not exist d ( s ) , s.t. C o r r ( B ˜ H ( t ) ( t ) , B ˜ H ( t ) ( s ) ) = O ( t d ( s ) ) . Similarly, when H ( t ) + H ( s ) < 1 MFBM will not have persistence of magnitude LRD in general in sense of Definition 2.
In definition of persistence of magnitudes, absolute values could be replaced by rth powers of absolute values for any r > 0 , similar to the result given in [18].
Proposition 2.
Let B ˜ H ( t ) ( t ) be a standard MFBM given in (4), then LRD holds in sense of Definition 2 for | X ( t ) | r , r > 0 . When,
(i)
s 0 , H ( t ) + H ( s ) < 1 and H ( t ) < 1 / 2 for sufficiently large t.
(ii)
H ( t ) + H ( s ) > 1 , H ( s ) > 1 / 2 for all s 0 and sufficiently large t.
Proof. 
The result again follows on the line of [18] for FBM. Let ( X , Y ) T be a bivariate Gaussian random vector as taken in Proposition 1. Then, for r > 0 (see e.g., [19]),
E | X Y | r = 2 r π σ 1 r σ 2 r Γ r + 1 2 2 F 1 1 2 r , 1 2 r ; 1 2 ; ρ 2 ,
where F 1 is the hypergeometric function [21]. In addition,
E | X | r = 2 r π σ 1 r Γ r + 1 2 , E | Y | r = 2 r π σ 2 r Γ r + 1 2 .
Moreover,
F 1 1 2 r , 1 2 r ; 1 2 ; ρ 2 = 1 + 1 2 r 2 ρ 2 + O ( ρ 4 ) a s ρ 0 ,
and Var ( | X | r ) = a r σ 1 2 r , Var ( | Y | r ) = a r σ 2 2 r , where,
a r = 2 r π Γ 2 r + 1 2 2 r π Γ r + 1 2 2 .
Thus,
Corr ( | B ˜ H ( t ) ( t ) | r , | B ˜ H ( s ) ( s ) | r ) = 2 r π a r Γ r + 1 2 2 F 1 1 2 r , 1 2 r ; 1 2 ; ρ ( t , s ) 2 1 ρ ( t , s ) 2 a s ρ ( t , s ) 0 ,
and the result follows with the help of Proposition 1(ii). □

3. The nth Order Multifractional Brownian Motion (n-MFBM)

Similar to nth-order fractional Brownian motion (n-FBM) [6], we define nth-order multifractional Brownian motion (n-MFBM) by subtracting n terms in the definition of MFBM. Precisely, that will allow functional parameter H in the range ( n 1 , n ) . Using [6], for the harmonizable representation of n-FBM denoted by B ˜ H n ( t ) , we have
E ( B ˜ H n ( t ) B ˜ H n ( s ) ) = 1 2 π 1 | ξ | 2 H + 1 e i t ξ k = 0 n 1 ( i t ξ ) k k ! e i s ξ k = 0 n 1 ( i s ξ ) k k ! d ( ξ ) = ( 1 ) n C H n 2 | t s | 2 H j = 0 n 1 ( 1 ) j 2 H j t s j | s | 2 H + s t j | t | 2 H ,
where C H n = 1 Γ ( 2 H + 1 ) | sin π H | . Analogous to the definition of n-FBM, we define n-MFBM as follows.
Definition 3.
Let H : [ 0 , ) [ a , b ] ( n 1 , n ) be a Hölder function of exponent β > 0 . For t 0 the following random function, denoted by B H ( t ) n ( t ) , is called nth-order multifractional Brownian motion with functional parameter H:
B H ( t ) n ( t ) = 1 Γ ( H ( t ) + 1 2 ) { 0 [ ( t u ) H ( t ) 1 / 2 ( u ) H ( t ) 1 / 2 ( H ( t ) 1 2 ) ( u ) H ( t ) 3 / 2 t ( H ( t ) 1 2 ) ( H ( t ) 2 n 3 2 ) ( u ) H ( t ) n + 1 / 2 t n 1 ( n 1 ) ! ] d B ( u ) + 0 t ( t u ) H ( t ) 1 / 2 d B ( u ) } . = R f + ( t , u ) d B ( u ) ,
where
f + ( t , u ) = 1 Γ ( H ( t ) + 1 2 ) ( t u ) + H ( t ) 1 / 2 k = 0 n 1 ( u ) + H ( t ) k 1 / 2 t k k ! i = 1 k ( H ( t ) i + 1 / 2 ) .
Remark 2.
It is easy to show that
f + ( t , u ) H ( t ) 1 2 n t n ( u ) H ( t ) 1 2 n , as u , t H ( t ) 1 2 , as u 0 , ( t u ) + H ( t ) 1 2 , as u t
Thus, we have | f + ( t , u ) | 2 d u < .
Remark 3.
The following Hölder functions can be used for H ( t ) ( n 1 , n ) , t > 0 :
H ( t ) = ( n 1 ) + 1 π arctan ( t ) + 1 2 , t > 0 ; H ( t ) = ( n 1 ) + a t b + a t , t > 0 ; H ( t ) = ( n 1 ) + a sin π 2 t α , 0 < a < 1 , 0 < α < 1 , t > 0 .
In addition, the definition of n-MFBM can be given in following harmonizable representation form.
Definition 4.
Let H : [ 0 , ) [ a , b ] ( n 1 , n ) be a Hölder function of exponent β > 0 . For t 0 the following random function, denoted by B ˜ H ( t ) n ( t ) , we call nth-order multifractional Brownian motion with functional parameter H:
B ˜ H ( t ) n ( t ) = R 1 | ξ | H ( t ) + 1 / 2 e i t ξ k = 0 n 1 ( i t ξ ) k k ! d B ^ ( ξ ) ,
where B ^ denotes a complex-valued Gaussian measure.
For n = 1 , the definitions given in Definitions 3 and 4 correspond to the definition of MFBM given in Definitions 1 and 4.
Proposition 3.
Let f ^ + ( ξ ) = F ( f + ) ( ξ ) be the Fourier transform of the function f + ( u ) L 2 ( d u ) , which is defined in (26). For all t R , H ( t ) ( n 1 , n ) , we have
f + ^ ( t , ξ ) = e i s i g n ( ξ ) ( H ( t ) + 1 / 2 ) π / 2 ( 2 π ) | ξ | H ( t ) + 1 / 2 e i t ξ k = 0 n 1 ( i t ξ ) k k ! , ξ R \ { 0 } .
Proof. 
The Fourier transform of function g ( u ) = u + j , u R , j ( 1 / 2 , 1 / 2 ) is given by (see [22], p. 171)
F ( u + j ) ( ξ ) = Γ ( j + 1 ) ( 2 π ) e i s i g n ( ξ ) ( j + 1 ) π / 2 , ξ 0 .
This is shown for 1 < R ( j ) < 0 by regularizing the function u + j at infinity by considering u + j e τ u taking the Fourier transform and letting τ 0 . This result can be extend to R ( j ) R { 1 , 2 , 3 , } due to analytic continuation.
By Using (29) with H ( t ) ( k , k + 1 ) , we have
F ( u + ) H ( t ) k 1 / 2 ( ξ ) = Γ ( H ( t ) k + 1 / 2 ) ( 2 π ) e i s i g n ( ξ ) ( H ( t ) + 1 / 2 ) π / 2 | ξ | H ( t ) + 1 / 2 ( i ξ ) k .
Then we get (see Proposition 3.1 [10])
F ( f + ( t , u ) ) ( ξ ) = 1 Γ ( H ( t ) + 1 / 2 ) [ F ( ( t u ) + H ( t ) 1 / 2 ) ( ξ ) k = 0 n 1 F ( ( u ) + H ( t ) k 1 / 2 ) ( ξ ) t k k ! i = 1 k ( H ( t ) i + 1 / 2 ) ] = e i s i g n ( ξ ) ( H ( t ) + 1 / 2 ) π / 2 ( 2 π ) | ξ | H ( t ) + 1 / 2 e i t ξ k = 0 n 1 ( i ξ t ) k k ! , ξ 0 .
Proposition 4.
The moving average representation of the n-MFBM (25) can also be represented as
B H ( t ) n ( t ) = cos ( ( H ( t ) + 1 / 2 ) π / 2 ) 2 π R 1 | ξ | H ( t ) + 1 / 2 e i t ξ k = 0 n 1 ( i ξ t ) k k ! d B ^ ( ξ ) + sin ( ( H ( t ) + 1 / 2 ) π / 2 ) 2 π R 1 i ξ | ξ | H ( t ) 1 / 2 e i t ξ k = 0 n 1 ( i ξ t ) k k ! d B ^ ( ξ ) .
Proof. 
By using the Proposition 2.1 in [10], the result follows. □
Remark 4.
Alternatively, using the definition of FBM given in ([23], p. 328), the n-MFBM can be defined as
B ¯ H ( t ) n ( t ) = R 1 i ξ | ξ | H ( t ) 1 / 2 e i t ξ k = 0 n 1 ( i t ξ ) k k ! d B ^ ( ξ ) .
The representations for n-MFBM given in (28) and (32) are equivalent since both are centered Gaussian process and
E ( B ¯ H ( t ) n ( t ) B ¯ H ( s ) n ( s ) ) = E ( B ˜ H ( t ) n ( t ) B ˜ H ( s ) n ( s ) ) = R 1 | ξ | H ( t ) + H ( s ) + 1 e i t ξ k = 0 n 1 ( i t ξ ) k k ! e i s ξ k = 0 n 1 ( i s ξ ) k k ! d ξ .
Remark 5.
When H ( t ) = n 1 / 2 , the moving average and harmonizable representations of n-MFBM have the same laws up to a multiplicative deterministic function. For n = 2 m + 1 , sin ( 2 m + 1 ) π / 2 = ( 1 ) m and cos ( 2 m + 1 ) π / 2 = 0 , which gives
B 2 m + 1 / 2 n ( t ) = ( 1 ) m 2 π B ¯ 2 m + 1 / 2 ( t ) .
For n = 2 m , we have sin ( m π ) = 0 and cos ( m π ) = ( 1 ) m , which leads to
B 2 m 1 / 2 n ( t ) = ( 1 ) m 2 π B ˜ 2 m 1 / 2 ( t ) .

Covariance Structure of n-MFBM

MFBM is a zero mean non-stationary Gaussian process with continuous sample paths. From definition it is obvious that n-MFBM is also a zero mean Gaussian process. Since n-MFBM is Gaussian, we will only be interested in statistical properties up to the second order, namely, the covariance function and variance of the process.
Proposition 5.
Let B ˜ H ( t ) n ( t ) be the n-MFBM given in (28) with functional parameter H ( t ) . Then,
Cov ( B ˜ H ( t ) n ( t ) , B ˜ H ( s ) n ( s ) ) = ( 1 ) n C ( t , s ) 2 [ | t s | H ( t ) + H ( s ) j = 0 n 1 ( 1 ) j H ( t ) + H ( s ) j × t s j | s | H ( t ) + H ( s ) + s t j | t | H ( t ) + H ( s ) ] ,
where C ( t , s ) is given by
C ( t , s ) = 2 π Γ ( H ( t ) + H ( s ) + 1 ) | sin π 2 ( H ( t ) + H ( s ) ) | .
Proof. 
By definition,
E ( B ˜ H ( t ) n ( t ) B ˜ H ( s ) n ( s ) ) = E [ R 1 | ξ | H ( t ) + 1 / 2 e i t ξ k = 0 n 1 ( i t ξ ) k k ! d B ^ ( ξ ) × R 1 | ξ | H ( s ) + 1 / 2 e i s ξ k = 0 n 1 ( i s ξ ) k k ! d B ^ ( ξ ) ] = R 1 | ξ | H ( t ) + H ( s ) + 1 e i t ξ k = 0 n 1 ( i t ξ ) k k ! e i s ξ k = 0 n 1 ( i s ξ ) k k ! d ξ .
Covariance function of standard n-FBM B H n ( t ) with Hurst parameter H, is given in (24). For fixed t and s, let B H n ( t ) be a standard n-FBM with Hurst parameter H = H ( t ) + H ( s ) 2 . Using (24), it follows
R 1 | ξ | H ( t ) + H ( s ) + 1 e i t ξ k = 0 n 1 i t ξ k ! e i s ξ k = 0 n 1 ( i s ξ ) k ! d ξ = ( 1 ) n C ( t , s ) 2 | t s | H ( t ) + H ( s ) j = 0 n 1 ( 1 ) j H ( t ) + H ( s ) j t s j | s | H ( t ) + H ( s ) + s t j | t | H ( t ) + H ( s ) ,
where
C ( t , s ) = 2 π Γ ( H ( t ) + H ( s ) + 1 ) | sin π 2 ( H ( t ) + H ( s ) ) | .
Using (35) and (36), the desired result follows. □
Remark 6.
For n = 1 and H ( t ) = H for all t > 0 , we have
E ( B ˜ H ( t ) n ( t ) B ˜ H ( s ) n ( s ) ) = C ( t , s ) 2 t 2 H + s 2 H | t s | 2 H
where C ( t , s ) = 2 π H Γ ( 2 H ) sin ( H π ) , which is the covariance function of FBM [6].
Remark 7.
Using (33), it follows
Var ( B ˜ H ( t ) n ( t ) ) = ( 1 ) n 1 C ( t , t ) | t | 2 H ( t ) j = 0 n 1 ( 1 ) j 2 H ( t ) j = C ( t , t ) | t | 2 H ( t ) 2 H ( t ) 1 n 1 ,
where C ( t , t ) = π H ( t ) Γ ( H ( t ) ) | sin ( π H ( t ) ) | .
Proposition 6.
Let B ˜ H ( t ) n ( t ) be the n-MFBM defined in (28). Then, for all admissible H ( t ) , the n-MFBM has LRD property in the sense of Definition 1.
Proof. 
We have
Corr ( B ˜ H ( t ) n ( t ) ) , B ˜ H ( s ) n ( s ) ) ) = 1 C ( t , t ) C ( s , s ) 2 H ( t ) 1 n 1 2 H ( s ) 1 n 1 t H ( t ) s H ( s ) [ t H ( t ) + H ( s ) ( 1 s t ) H ( t ) + H ( s ) j = 0 n 1 ( 1 ) j H ( t ) + H ( s ) j t j s H ( t ) + H ( s ) j + s j t H ( t ) + H ( s ) j ] = D ( t , s ) t s H ( s ) 1 s t H ( t ) + H ( s ) j = 0 n 1 ( 1 ) j H ( t ) + H ( s ) j s t H ( t ) + H ( s ) j + s t j = D ( t , s ) t s H ( s ) j = n ( 1 ) j H ( t ) + H ( s ) j s t j j = 0 n 1 ( 1 ) j H ( t ) + H ( s ) j s t H ( t ) + H ( s ) j s t H ( t ) n + 1 if H ( t ) + H ( s ) 2 n 1 / 2 s t n H ( s ) if H ( t ) + H ( s ) 2 n 1 / 2 .

4. Conclusions

Brownian motion is a central object in stochastic calculus and stochastic differential equations. However due to independent increments, Brownian motion has limitations in modeling random phenomena exhibiting dependence. Fractional Brownian motion (FBM) is a generalization of Brownian motion possessing LRD, persistence of signs LRD and persistence of magnitude LRD properties. The pointwise irregularities in FBM do not vary with time and hence it cannot be used to simulate artificial mountains for example. To overcome these limitations MFBM was introduced in literature where pointwise irregularities represented by a Hölder continuous function. In this article, we generalize the MFBM to n-MFBM which allows the Hölder continuous function to take values in the interval ( n 1 , n ) . The basic properties like covariance and variance of the introduced process are obtained. It is shown that n-MFBM has LRD property. Further, the harmonizable representation of n-MFBM is given. We believe that n-MFBM will be useful in modeling of natural phenomena where the pointwise irregularities vary with time and the Hurst parameter does not necessarily belong to the interval ( 0 , 1 ) .

5. Discussion

Note that the tempered fractional stable motion introduced in [15] adds an exponential tempering to the power-law kernel in the LFSM. In the line of higher order fractional stable motion [16], it is very natural to define a higher order tempered fractional stable motion. Furthermore, the higher order fractional stable motion and the tempered fractional stable motion can be extended further by choosing the Hurst parameter as a stochastic process instead of deterministic function of time as done in [7]. In addition, the GMFBM can be further theorized to higher order GMFBM and higher order generalized MFBM.

Author Contributions

Conceptualization, N.G., A.K. and N.L.; Methodology, A.K. and N.L.; Writing—Original Draft Preparation, N.G., A.K. and N.L.; Writing—Review & Editing, N.G., A.K. and N.L. All authors have read and agreed to the published version of the manuscript.

Funding

N.G. would like to thank Council of Scientific and Industrial Research (CSIR), India for supporting her research under the fellowship award number 09/1005(0021)2018-EMR-I. Further, A.K. would like to express his gratitude to Science and Engineering Research Board (SERB), India for the financial support under the MATRICS research grant MTR/2019/000286. N.L. was supported by the Grant of London Mathematical Society, Ref. 42007.

Acknowledgments

N.G. would like to thank Council of Scientific and Industrial Research (CSIR), India, for the award of a research fellowship. N.L. was partially supported by LMS grant 42997 and Isak Newton Institute (Cambridge) under the Program “Fractional Differential Equations”.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
MFBMMultifractional Brownian motion
FBMFractional Brownian motion
LRDLong range dependence
SRDShort range dependence
n-MFBMnth order multifractional Brownian motion.

References

  1. Mandelbrot, B.B.; Van Ness, J.W. Fractional Brownian motions, fractional Gaussian noises and applications. SIAM Rev. 1968, 10, 422–437. [Google Scholar] [CrossRef]
  2. Grahovac, D.; Leonenko, N.N.; Taqqu, M.S. Scaling properties of the empirical structure function of linear fractional stable motion and estimation of its parameters. J. Stat. Phys. 2015, 158, 105–119. [Google Scholar] [CrossRef] [Green Version]
  3. Echelard, A.; Véhel, J.L.; Barriére, O. Terrain Modeling with Multifractional Brownian Motion and Self-regulating Processes. In Proceedings of the International Conference on Computer Vision and Graphics (IICCVG 2010), Warsaw, Poland, 20–22 September 2010; pp. 342–351. [Google Scholar]
  4. Islam, A.; Reza, S.M.; Iftekharuddin, K.M. Multifractal texture estimation for detection and segmentation of brain tumors. IEEE Trans. Biomed. Eng. 2013, 60, 3204–3215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Beran, J. Statistics for Long-Memory Processes; Chapman & Hall: New York, NY, USA, 1994. [Google Scholar]
  6. Perrin, E.; Harba, R.; Berzin, C.; Iribarren, I.; Bonami, A. nth-Order Fractional Brownian Motion and Fractional Gaussian Noises. IEEE Trans. Signal Process. 2001, 49, 1049–1059. [Google Scholar] [CrossRef]
  7. Ayache, A.; Taqqu, M.S. Multifractional processes with random exponent. Publ. Mat. 2005, 49, 459–486. [Google Scholar] [CrossRef] [Green Version]
  8. Peltier, R.; Lévy-Vehel, J. Multifractional Brownian Motion: Definition and Preliminary Results; Res. Rept. 2645; INRIA: Le Chesnay-Rocquencourt, France, 1995; Available online: https://hal.inria.fr/inria-00074045/document (accessed on 26 January 2021).
  9. Chan, G.; Wood, A.T. Simulation of multifractional Brownian motion. In COMPSTAT; Springer: Berlin/Heidelberg, Germany, 1998; pp. 233–238. [Google Scholar]
  10. Stoev, S.A.; Taqqu, M.S. How rich is the class of multifractional Brownian motions? Stoch. Process. Appl. 2006, 116, 200–221. [Google Scholar] [CrossRef] [Green Version]
  11. Benassi, A.; Jaffard, S.; Roux, D. Elliptic Gaussian random processes. Rev. Mat. Iberoam. 1997, 13, 19–90. [Google Scholar] [CrossRef] [Green Version]
  12. Cohen, S. From self-similarity to local self-similarity: The estimation problem. In Fractals: Theory and Applications in Engineering; Dekking, M., Lévy-Véhel, J., Lutton, E., Tricot, C., Eds.; Springer: London, UK, 1999. [Google Scholar]
  13. Ayache, A.; Cohen, S.; Vehel, J.L. The covariance structure of multifractional Brownian motion with application to long-range dependence. In Proceedings of the 2000 IEEE International Conference on Acoustics, Speech, and Signal Processing, Istanbul, Turkey, 5–9 June 2000. [Google Scholar]
  14. Stoev, S.; Taqqu, M.S. Stochastic Properties of the Linear Multifractional Stable Motion. Adv. Appl. Probab. 2004, 36, 1085–1115. [Google Scholar] [CrossRef]
  15. Meerschaert, M.M.; Sabzikar, F. Tempered Fractional Stable Motion. J. Theor. Probab. 2016, 29, 681–706. [Google Scholar] [CrossRef]
  16. Kawai, R. Higher Order Fractional Stable Motion: Hyperdiffusion with Heavy Tails. J. Stat. Phys. 2016, 165, 126–152. [Google Scholar] [CrossRef]
  17. Ayache, A.; Lévy-Véhel, J. Generalized multifractional Brownian motion: Definition and preliminary results. In Fractals: Theory and Applications in Engineering; Dekking, M., Lévy-Véhel, J., Lutton, E., Tricot, C., Eds.; Springer: London, UK, 1999; pp. 17–32. [Google Scholar]
  18. Heyde, C.C. On Modes of Long-Range Dependence. J. Appl. Prob. 2002, 39, 882–888. [Google Scholar] [CrossRef]
  19. Johnson, N.L.; Kotz, S. Distributions in Statistics: Continuous Multivariate Distributions; John Wiley: New York, NY, USA, 1972. [Google Scholar]
  20. Chandramouli, R.; Ranganathan, N. Computing the bivariate Gaussian probability integral. IEEE Signal Process. Lett. 1999, 6, 129–131. [Google Scholar] [CrossRef] [Green Version]
  21. Nabeya, S. Absolute moments in 2-dimensional normal distribution. Ann. Inst. Stat. Math. 1951, 3, 2–6. [Google Scholar] [CrossRef]
  22. Gelfand, I.M.; Shilov, G.E. Generalized Functions: Properties and Operations; Academic Press: New York, NY, USA, 1964. [Google Scholar]
  23. Samoradnitsky, G.; Taqqu, M.S. Stable Non-Gaussian Random Processes: Stochastic Models with Infinite Variance; Chapman and Hall/CRC: Boca Raton, FL, USA, 1994. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gupta, N.; Kumar, A.; Leonenko, N. A Generalization of Multifractional Brownian Motion. Fractal Fract. 2022, 6, 74. https://doi.org/10.3390/fractalfract6020074

AMA Style

Gupta N, Kumar A, Leonenko N. A Generalization of Multifractional Brownian Motion. Fractal and Fractional. 2022; 6(2):74. https://doi.org/10.3390/fractalfract6020074

Chicago/Turabian Style

Gupta, Neha, Arun Kumar, and Nikolai Leonenko. 2022. "A Generalization of Multifractional Brownian Motion" Fractal and Fractional 6, no. 2: 74. https://doi.org/10.3390/fractalfract6020074

APA Style

Gupta, N., Kumar, A., & Leonenko, N. (2022). A Generalization of Multifractional Brownian Motion. Fractal and Fractional, 6(2), 74. https://doi.org/10.3390/fractalfract6020074

Article Metrics

Back to TopTop