Next Article in Journal
Heart Rhythm Analysis Using Nonlinear Oscillators with Duffing-Type Connections
Next Article in Special Issue
The Fractional Soliton Wave Propagation of Non-Linear Volatility and Option Pricing Systems with a Sensitive Demonstration
Previous Article in Journal
Fractional Transformation-Based Decentralized Robust Control of a Coupled-Tank System for Industrial Applications
Previous Article in Special Issue
Parameter Estimation in Rough Bessel Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Correlation Structure of Time-Changed Generalized Mixed Fractional Brownian Motion

1
Department of Mathematics, College of Science, Imam Abdulrahman Bin Faisal University, P.O. Box 1982, Dammam 31441, Saudi Arabia
2
Basic and Applied Scientific Research Center, Imam Abdulrahman Bin Faisal University, P.O. Box 1982, Dammam 31441, Saudi Arabia
Fractal Fract. 2023, 7(8), 591; https://doi.org/10.3390/fractalfract7080591
Submission received: 29 May 2023 / Revised: 5 July 2023 / Accepted: 19 July 2023 / Published: 30 July 2023
(This article belongs to the Special Issue Fractional Processes and Multidisciplinary Applications)

Abstract

:
The generalized mixed fractional Brownian motion (gmfBm) is a Gaussian process with stationary increments that exhibits long-range dependence controlled by its Hurst indices. It is defined by taking linear combinations of a finite number of independent fractional Brownian motions with different Hurst indices. In this paper, we investigate the long-time behavior of gmfBm when it is time-changed by a tempered stable subordinator or a gamma process. As a main result, we show that the time-changed process exhibits a long-range dependence property under some conditions on the Hurst indices. The time-changed gmfBm can be used to model natural phenomena that exhibit long-range dependence, even when the underlying process is not itself long-range dependent.

1. Introduction

The fractional Brownian motion (fBm) B H = { B t H , t 0 } with parameter H, is a centered Gaussian process with the covariance function
C o v ( B t H , B s H ) = 1 2 [ t 2 H + s 2 H | t s | 2 H ] , s , t 0 ,
where H is a real number in ( 0 , 1 ) , called the Hurst index. The case H = 1 2 corresponds to the Brownian motion (Bm).
An extension of the fBm was introduced by Cheridito [1], called the mixed fractional Brownian motion (mfBm), which is a linear combination of a Bm and an independent fBm of Hurst index H, with stationary increments, that exhibit a long-range dependence for H > 1 2 . A mfBm of parameters a 1 ,   a 2 and H is a process M H ( a 1 , a 2 ) = { M t H ( a 1 , a 2 ) , t 0 } , defined on some probability space ( Ω , F , P ) by
M t H ( a 1 , a 2 ) = a 1 B t + a 2 B t H , t 0 ,
where B = { B t , t 0 } is a Bm and B H = { B t H , t 0 } is an independent fBm of Hurst index H ( 0 , 1 ) . We refer also to [1,2,3,4] for further information on the mfBm process.
C. Elnouty [3] propose a generalization of the mfBm called fractional mixed fractional Brownian motion (fmfBm) of parameters a 1 , a 2 and H = ( H 1 , H 2 ) . A fmfBm is a process N H ( a 1 , a 2 ) = { N t H ( a 1 , a 2 ) , t 0 } , defined on some probability space ( Ω , F , P ) by
N t H ( a 1 , a 2 ) = a 1 B t H 1 + a 2 B t H 2 , t 0 ,
where B H i = { B t H i , t 0 } are independent fBms of Hurst indices H i ( 0 , 1 ) for i = 1 , 2 . In addition, the fmfBm was studied by Miao, Y et al. [5].
The fractional mixed fractional Brownian motion was further generalized by Thäle in 2009 [6] to the generalized mixed fractional Brownian motion. Let a 1 , , a n , n N * reel numbers and not all a i equals zero. A generalized mixed fractional Brownian motion (gmfBm) of parameters H = ( H 1 , H 2 , , H n ) and a = ( a 1 , a 2 , , a n ) ; H i ( 0 , 1 ) , a i R , n N * is a stochastic process Z H , a = { Z t H , a , t 0 } defined on some probability space ( Ω , F , P ) by
Z t H , a = i = 1 n a i B t H i ,
where { B t H i , t 0 } for i = 1 , 2 , , n are independent fBms of Hurst indices H i ( 0 , 1 ) . The gmfBm is a generalization of both fractional Brownian motion and subfractional Brownian motion. The gmfBm is a centered Gaussian process with stationary increments that exhibits long-range dependence if and only if there exists some k with H k > 1 2 , it can be used to model a wider range of natural phenomena than either fBm or sfBm. Internet traffic can be modeled using the gmfBm, as seen in [7]. The gmfBm market is a useful model for a variety of assets, including internet traffic. Internet traffic has been shown to exhibit long-range dependence, and the gmfBm model can be used to capture this dependence. The gmfBm market is a market where the underlying asset price satisfies the following stochastic differential equation:
d S t = a S t d t + b S t d Z t
where a and b stand for the standard deviation of the stock return and the volatility, see [8].
It should be noted that the gmfBm model is a generalization of all the fractional models studied in the literature. This generalized model degenerates into a single fBm model with n = 1 , a Bm model with n = 1 and H 1 = 1 2 , an mfBm model with n = 2 and H 1 = 1 2 and a fmfBm with n = 2 . For a detailed survey on the properties of the gmfBm, we refer to [6,9,10].
The time-changed generalized mixed fractional Brownian motion is defined as
T β H , a = { T β t H , a , t 0 } = { Z β t H , a , t 0 } ,
where the parent process T H , a is a gmfBm with parameters H = ( H 1 , H 2 , , H n ) , a = ( a 1 , a 2 , , a n ) , H i ( 0 , 1 ) , a i R , n N * and the internal process is a subordinator β = { β t , t 0 } assumed to be independent of B t H i , for i = 1 , 2 , , n . If H = ( 1 2 , 0 , , 0 ) and a = ( 1 , 0 , , 0 ) , the process T β H , a is called subordinated Brownian motion. Also, the process T β H , α , for H = ( H 1 , 0 , , 0 ) and a = ( 1 , 0 , , 0 ) is called subordinated fractional Brownian motion it was considered in [11,12].
A time-changed process is a stochastic process that is constructed by taking the superposition of two independent stochastic systems. The first system is called the external process, and the second system is called the subordinator. The evolution of time in the external process is replaced by the subordinator, which is a non-decreasing stochastic process. The resulting time-changed process often retains important properties of the base process; however, certain characteristics may change. The idea of subordination was introduced by Bochner in 1949 [13], and it has been explored in many papers since then (e.g., [14,15,16,17]). Subordination is a versatile tool that can be used to construct a wide variety of stochastic processes. It is a powerful tool for modeling real-world phenomena, and it has been used in many different fields, including finance, insurance, and physics.
In the case that H = ( 1 2 , H 2 , 0 , 0 ) and a = ( a 1 , a 2 , 0 , , 0 ) , the time-changed mixed fractional Brownian motion was discussed in [18] to present a stochastic model of the discounted stock price in some arbitrage-free and complete financial markets.
The time-changed processes have found many interesting applications, for example [18,19,20,21,22,23,24,25,26].
This study investigates the long-range dependency property of the time-changed gmfBm. We describe two processes that make up gmfBm’s “operational time”. In the first scenario, the internal process, which plays the role of time, is the tempered stable subordinator, whereas, in the second situation, it is the gamma process. As an application, we deduce the results concerned the long-range dependence property of the time-changed fBm by tempered stable subordinator and gamma process proved by Kumar et al. in [11,12], respectively.

2. Preliminaries

We define the tempered stable subordinator and gamma process in this section. Additionally, we quickly review the definitions of long-range dependence based on a process’s correlation function.
A subordinator is a process with stationary and independent non-negative increments starting at zero. Subordinators are a special class of Lévy processes taking values in [ 0 , ) and their sample paths are non-decreasing; this is a type of stochastic process that is used to model random phenomena that have jumps (see [27,28] for more details). Let β = { β t , t 0 } be a subordinator. The infinite divisibility of the law of β implies that its Laplace transform can be expressed in the form
E ( e λ β t ) = e t Φ ( λ ) , λ > 0 ,
where Φ : [ 0 , ) [ 0 , ) , called the Laplace exponent, is a Bernstein function. Such that the Laplace exponent Φ can be expressed as
Φ ( λ ) = a + b λ + ( 0 , ) ( 1 e λ τ ) σ ( d τ ) < ,
which is known as the Lévy-Khintchine formula for the subordinator β . Where a , b 0 and σ are a measure on the positive real half-line such that
( 0 , ) ( 1 τ ) σ ( d τ ) <

2.1. Tempered Stable Subordinator

Tempered stable subordinator, where index α ( 0 , 1 ) and tempering parameter λ > 0 (TSS) are the non-decreasing and non-negative Lévy process S λ , α = { S t λ , α , t 0 } with density function:
f λ , α ( x , t ) = e x p ( λ x + λ α t ) f α ( x , t ) , λ > 0 , α ( 0 , 1 ) ,
where
f α ( x , t ) = 1 π 0 e x y e t y α cos α π sin ( t y α sin α π ) d y .
More details about TSS can be found in [11].
Lemma 1.
(see [11] for the proof)
For q > 0 , the asymptotic behavior of q-th order moments of S t λ , α satisfies
E ( S t λ , α ) q ( α λ α 1 t ) q , a s t .

2.2. Gamma Subordinator

Gamma subordinator Γ = { Γ t , t } is a stationary and independent increments process with gamma distribution. More precisely, the increment Γ t + s Γ s has the density function
f ( x , t ) = 1 Γ ( t / ν ) x ( t / ν ) 1 e x , x > 0 , ν > 0 .
More details about gamma subordinator can be found in [12].
Lemma 2.
(see [12] for the proof)
For q > 0 , the asymptotic behavior of q-th order moments of Γ t satisfies
E ( Γ t ) q t ν q , a s t .
Lemma 3.
(see [12] for the proof)
The covariance of Γ t is
C o v ( Γ s , Γ t ) = t ν + t 2 s 2 ν 2 , w h e r e s < t .
Then for fixed s and t , it follows that
C o v ( Γ s , Γ t ) t 2 ν 2 .

2.3. Long-Range Dependence

Notation 1.
Let X and Y be two random variables defined on the same probability space ( Ω , F , P ) . We denote the correlation coefficient C o r r ( X , Y ) by
C o r r ( X , Y ) = C o v ( X , Y ) V a r ( X ) V a r ( Y ) .
Definition 1.
Please note that a finite variance stationary process { X t , t 0 } is said to have a long-range dependence property (Cont and Tankov [29]), if k = 0 γ k = , where
γ k = C o v ( X t , X t + k ) .
In the following definition, we give the equivalent definition for a non-stationary process { X t , t 0 } .
Definition 2.
Let s > 0 be fixed and t > s . The process { X t , t 0 } is said to have a long-range dependence property if
C o r r ( X t , X s ) c ( s ) t d , a s t ,
where c ( s ) is a constant depending on s and d ( 0 , 1 ) . An equivalent definition given in [30].
Let 0 < s < t and s be fixed. Assume a stochastic process { X t , t 0 } has the correlation function C o r r ( X ( s ) , X ( t ) ) that satisfies c 1 ( s ) t d C o r r ( X s , X t ) c 2 ( s ) t d for large t , d > 0 , c 1 ( s ) > 0 and c 2 ( s ) > 0 . , i.e.,
l i m t C o r r ( X s , X t ) t d = c ( s )
for some c ( s ) > 0 and d > 0 . We say { X t , t 0 } has the long-range dependence property if d ( 0 , 1 ) and has the short-range dependence property if d ( 1 , 2 ) .
Proposition 1.
The TSS with index α ( 0 , 1 ) and tempering parameter λ > 0 has LRD property.
Proof. 
First, we compute the covariance function using the subordinator’s independent increment characteristic. For 0 < s < t , we have
C o v ( S s λ , α , S t λ , α ) = C o v ( S s λ , α , ( S t λ , α S s λ , α ) + S s λ , α ) = C o v ( S s λ , α , ( S t λ , α S s λ , α ) ) + C o v ( S s λ , α , S s λ , α ) = V a r ( S s λ , α ) .
Thus, the correlation function is given by
C o r r ( S s λ , α , S t λ , α ) = C o v ( S s λ , α , S t λ , α ) V a r ( S s λ , α ) 1 2 V a r ( S t λ , α ) 1 2 = V a r ( S s λ , α ) 1 2 V a r ( S t λ , α ) 1 2 = s 1 2 t 1 2
Hence,
l i m t C o r r ( S s λ , α , S t λ , α ) t 1 2 = s 1 2
Therefore, the TSS { S t λ , α , t 0 } has an LRD property. □
Similar to the proof of Proposition 1, we obtain
Proposition 2.
The gamma process has a long-range dependence property.
Definition 3.
Let a 1 , , a n , reel numbers and not all a i equals zero. A generalized mixed fractional Brownian motion (gmfBm) of parameters H = ( H 1 , H 2 , , H n ) and a = ( a 1 , a 2 , , a n ) ; H i ( 0 , 1 ) , a i R , n N * is a stochastic process Z H , a = { Z t H , a , t 0 } defined on some probability space ( Ω , F , P ) by
Z t H , a = i = 1 n a i B t H i ,
where { B t H i , t 0 } for i = 1 , 2 , , n are independent fractional Brownian motions of Hurst indices H i ( 0 , 1 ) .
Lemma 4.
(see [6] for the proof) The gmfBm has stationary increments and exhibits a long-range dependence property if, and only if, there exist some k { 1 , , n } with H k > 1 2 .

3. gmfBm Time-Changed by Tempered Stable Subordinator

In this section, we will investigate the gmfBm time-changed by tempered stable subordinator.
Definition 4.
Let n N * . Let Z H , a = { Z t H , a , t 0 } be a gmfBm of parameters a = ( a 1 , a 2 , , a n ) and H = ( H 1 , , H n ) , H i ( 0 , 1 ) , a i R , n N * . Let S λ , α = { S t λ , α , t 0 } be a TSS with index α ( 0 , 1 ) and tempering parameter λ > 0 . The time-changed process of T H , a by means of S λ , α is the process T S λ , α H , a = { T S t λ , α H , a , t 0 } defined by:
T S t λ , α H , a = i = 1 n a i B S t λ , α H i ,
where the subordinator S t λ , α is assumed to be independent of all B H i for i = 1 , 2 , , n .
Proposition 3.
Let T H , a = { T t H , a , t 0 } be a gmfBm of parameters a = ( a 1 , , a n ) , and H = ( H 1 , , H n ) , H i ( 0 , 1 ) , a i R , n N * . Let T S λ , α H , a be the gmfBm time-changed by S λ , α . Then by Taylor’s expansion we obtain, for fixed s and large t,
C o v ( T S t λ , α H , a , T S s λ , α H , a ) i = 1 n a i 2 H i s ( α λ α 1 ) 2 H i t 2 H i 1 , a s t .
Proof. 
Let s > 0 be fixed and let s t . The covariance function of T S t λ , α H , a and T S s λ , α H , a is defined by
C o v ( T S t λ , α H , a , T S s λ , α H , a ) = E ( T S t λ , α H , a T S s λ , α H , a ) E ( T S t λ , α H , a ) E ( T S s λ , α H , a )
by observing that E [ T S t λ , α H , a ] = 0 , t 0 and using ([11], p. 195), the process T S t λ , α H , a follows
C o v ( T S t λ , α H , a , T S s λ , α H , a ) = E [ T S t λ , α H , a T S s λ , α H , a ] = 1 2 E [ ( T S t λ , α H , a ) 2 + ( T S s λ , α H , a ) 2 ( T S t λ , α H , a T S s λ , α H , a ) 2 ] = 1 2 E [ ( Z S t λ , α H , a ) 2 + ( Z S s λ , α H , a ) 2 ( Z S t λ , α H , a Z S s λ , α H , a ) 2 ] = 1 2 E [ ( i = 1 n a i B S t λ , α H i ) 2 + ( i = 1 n a i B S s λ , α H i ) 2 ] 1 2 E [ ( i = 1 n a i ( B S t λ , α H i B S s λ , α H i ) ) 2 ] .
Since the fractional Brownian motion has stationary increments, then
C o v ( T S t λ , α H , a , T S s λ , α H , a ) = 1 2 E [ ( i = 1 n a i B S t λ , α H i ) 2 + ( i = 1 n a i B S s λ , α H i ) 2 ( i = 1 n a i B S t s λ , α H i ) 2 ] = 1 2 E [ i = 1 n ( a i B S t λ , α H 1 ) 2 + 2 i j a i a j B S t λ , α H i B S t λ , α H j ] + 1 2 E [ i = 1 n ( a i B S s λ , α H i ) 2 + 2 i j n a i b j B S s λ , α H i B S s λ , α H j ] 1 2 E [ i = 1 n ( a i B S t s λ , α H i ) 2 + 2 i j a i b j B S t s λ , α H i B S t s λ , α H j ] .
By the independence of the fBms’ B t H i for i = 1 , , n and their independence of the S λ , α , we have
E [ B S t λ , α H k B S t λ , α H l ] = E [ E ( B r H k B r H l | S t λ , α ) ] = E [ B r H k B r H l ] f S t λ , α ( d r ) = 0
where f S t λ , α ( . ) is the distribution function of S t λ , α .
Thus, we obtain
C o v ( T S t λ , α H , a , T S s λ , α H , a ) = i = 1 n a i 2 2 [ E ( B S t λ , α H i ) 2 + E ( B S s λ , α H i ) 2 E ( B S t s λ , α H i ) 2 ] = i = 1 n a i 2 2 E ( B H i ( 1 ) ) 2 [ E ( S t λ , α ) 2 H i + E ( S s λ , α ) 2 H i E ( S t s λ , α ) 2 H i ] .
Hence for large t and using Lemma 1, we have
C o v ( T S t λ , α H , a , T S s λ , α H , a ) i = 1 n a i 2 2 [ ( α λ α 1 ) 2 H i t 2 H i + E ( S s λ , α ) 2 H i ( α λ α 1 ) 2 H i ( t s ) 2 H i ] = i = 1 n a i 2 2 ( α λ α 1 ) 2 H i t 2 H i ( 2 H i s t + E ( S s λ , α ) 2 H i t 2 H i + O ( t 2 ) ) i = 1 n a i 2 H i s ( α λ α 1 ) 2 H i t 2 H i 1 .
Proposition 4.
Let T H , a = { T t H , a , t 0 } be a gmfBm of parameters a = ( a 1 , , a n ) , and H = ( H 1 , , H n ) , H i ( 0 , 1 ) , a i R , n N * . Let S λ , α = { S t λ , α , t 0 } be the TSS with index α ( 0 , 1 ) and tempering parameter λ > 0 and let T S λ , α H , a be the gmfBm time-changed process by means of S λ , α . Then for fixed s > 0 and t , we obtain
E [ ( T S t λ , α H , a T S s λ , α H , a ) 2 ] i = 1 n a i 2 H i ( α λ α 1 ) 2 H i t 2 H i i = 1 n 2 a i 2 H i s ( α λ α 1 ) 2 H i t 2 H i 1 + i = 1 n a 2 H i ( α λ α 1 ) 2 H i s 2 H i .
Proof. 
Let s > 0 be fixed and t . Then, using Equation (4), we have
E [ ( T S t λ , α H , a T S s λ , α H , a ) 2 ] = E [ ( T S t λ , α H , a T S s λ , α H , a ) ( T S t λ , α H , a T S s λ , α H , a ) ] = E [ ( T S t λ , α H , a ) 2 T S t λ , α H , a T S s λ , α H , a T S t λ , α H , a T S s λ , α H , a + ( T S s λ , α H , a ) 2 ] = E [ ( T S t λ , α H ) 2 2 T S t λ , α H , a T S s λ , α H , a + ( T S s λ , α H , a ) 2 ] i = 1 n a i 2 2 [ ( α λ α 1 ) 2 H i t 2 H i + E ( S s λ , α ) 2 H i ( α λ α 1 ) 2 H i ( t s ) 2 H i ] = i = 1 n a i 2 2 ( α λ α 1 ) 2 H i t 2 H i ( 2 H i s t + E ( S s λ , α ) 2 H i t 2 H i + O ( t 2 ) ) i = 1 n a i 2 H i s ( α λ α 1 ) 2 H i t 2 H i 1 .
Now we discuss the long-range dependence behavior of T S λ , α H , a .
Theorem 1.
Let T H , a = { T t H , a , t 0 } be a gmfBm of parameters a = ( a 1 , , a n ) and H = ( H 1 , , H n ) , H i ( 0 , 1 ) , a i R , n N * with H i < H n for i = 1 , 2 , , n 1 . Let S λ , α = { S t λ , α , t 0 } be the TSS with index α ( 0 , 1 ) and tempering parameter λ > 0 . Then the time-changed gmfBm by means of S λ , α exhibits a long-range dependence property for all Hurst indices satisfying 0 < 2 H i H n < 1 .
Proof. 
Let n N * . Let T H , a = { T t H , a , t 0 } be a gmfBm of parameters a = ( a 1 , , a n ) and H = ( H 1 , , H n ) , H i ( 0 , 1 ) , a i R , n N * with H i < H n for i = 1 , 2 , , n 1 . Let S λ , α = { S t λ , α , t 0 } be the TSS with index α ( 0 , 1 ) and tempering parameter λ > 0 . The process T S λ , α H , a is not stationary, hence the Definition 2 will be used to establish the long-range dependence property.
Using Equations (2), (4) and by Taylor’s expansion we obtain, as t
C o r r ( T S t λ , α H , a , T S s λ , α H , a ) i = 1 n a i 2 H i s ( α λ α 1 ) 2 H i t 2 H i 1 i = 1 n a i 2 H i ( α λ α 1 ) 2 H i t 2 H i E ( T S s λ , α H , a ) 2 = i = 1 n a i 2 H i s ( α λ α 1 ) 2 H i t 2 H i 1 a n 2 H n ( α λ α 1 ) 2 H n σ 2 t 2 H n i = 1 n 1 ( H i a i 2 a n 2 H n ( α λ α 1 ) 2 H i 2 H n t 2 H i 2 H n ) + 1 i = 1 n a i 2 H i s ( α λ α 1 ) 2 H i t 2 H i 1 a n 2 H n σ 2 ( α λ α 1 ) 2 H n t 2 H n , s i n c e H i < H n , = i = 1 n 1 a i 2 H i H n 1 2 s ( α λ α 1 ) H n 2 H i t 2 H i H n 1 | a n | σ + | a n | H n 1 2 s ( α λ α 1 ) H n t H n 1 σ .
where σ 2 = E ( T S s λ , α H , a ) 2 . Then, for fixed s > 0 and t . For H i < H n < 1 for i = 1 , , n 1 , the correlation function is given by
C o r r ( T S t λ , α H , a , T S s λ , α H , a ) c 1 t 2 H i H n 1 + c 2 t H n 1 .
Therefore, for 0 < 2 H i H n < 1 the correlation function of T Γ t H , a decays like a t ( 1 H n ) for all 0 < H n < 1 . Then, in the sense of Definition 2, the time-changed process T S λ , α H , a exhibits a long-range dependence property for all 0 < 2 H i H n < 1 , i = 1 , 2 , , n 1 . □
Remark 1.
When n = 1 in Equation (4) and using Equation (2), we obtain
C o v ( T S t λ , α H , a , T S s λ , α H , a ) H 1 s ( α λ α 1 ) 2 H 1 t 2 H 1 1 , a s t , C o r r ( T S t λ , α H , a , T S s λ , α H , a ) H 1 s 1 H 1 t H 1 1 , a s t .
Hence we obtain the following results, proven in [11].
Corollary 1.
The fractional Brownian motion time-changed by TSS has long-range dependence for every H ( 0 , 1 ) .

4. gmfBm Time-Changed by the Gamma Subordinator

This section looks into generalized mixed fractional Brownian motion time-changed by the gamma process.
Definition 5.
Let Z H , a = { Z t H , a , t 0 } be a gmfBm of parameters a = ( a 1 , , a n ) and H = ( H 1 , , H n ) , H i ( 0 , 1 ) , a i R , n N * . Let Γ = { Γ t , t 0 } be a gamma process. The time-changed process of Z H , a by means of Γ is the process T Γ H , a = { T Γ t H , a , t 0 } defined by:
T Γ t H , a = i = 1 n a i B Γ t H i ,
where the process Γ t is assumed to be independent of B t H i for i = 1 , , n .
Proposition 5.
Let Z H , a = { Z t H , a , t 0 } be a gmfBm of parameters a = ( a 1 , , a n ) and H = ( H 1 , , H n ) , H i ( 0 , 1 ) , a i R , n N * . Let T Γ H , a be the gmfBm time-changed by Γ. Then we have
  • For t > s , the covariance function for the process T Γ t H , a follows
    C o v ( T Γ t H , a , T Γ s H , a ) = i = 1 n a i 2 2 Γ ( 2 H i + t / ν ) Γ ( t / ν ) + Γ ( 2 H i + s / ν ) Γ ( s / ν ) Γ ( 2 H i + ( t s ) / ν ) Γ ( ( t s ) / ν )
  • For fixed s and large t, the process Y Γ t H follows
    C o v ( T Γ t H , a , T Γ s H , a ) i = 1 n 2 a i 2 H i s ν 2 H i t 2 H i 1 .
Proof. 
1.
Let s fixed. Let s < t . We use a similar procedure as in the proof of Equation (5). By the independence of the fBms’ B t H i for i = 1 , , n and their independence of the Γ t , we obtain
C o v ( T Γ t H , a , T Γ s H , a ) = E ( T Γ t H , a T Γ s H , a ) = i = 1 n a i 2 2 E ( B Γ t H i ) 2 + E ( B Γ s H i ) 2 E ( B Γ t s H i ) 2 = i = 1 n a i 2 2 Γ ( 2 H i + t / ν ) Γ ( t / ν ) + Γ ( 2 H i + s / ν ) Γ ( s / ν ) Γ ( 2 H i + ( t s ) / ν ) Γ ( ( t s ) / ν ) .
2.
Let g ( x ) = Γ ( x + 2 H i ) / Γ ( x ) and f ( x ) = Γ ( x + 1 ) / Γ ( x ) . By Taylor expansion and [12] we have
g ( x + h ) g ( x ) = 1 + 2 H i ( h / x ) + H ( 2 H i 1 ) ( h / x ) 2 + O ( x 3 ) ,
and
f ( x + h ) f ( x ) = 1 + ( h / x ) + O ( x 2 ) .
For fixed s and large t, using Equations (7)–(9), T Γ t H follows
C o v ( T Γ t H , a , T Γ s H , a ) = i = 1 n a i 2 2 g ( t / ν ) 1 + g ( s / ν ) g ( t / ν ) g ( ( t s ) / ν ) g ( t / ν ) = i = 1 n a i 2 2 ( t / ν ) 2 H i [ 1 + g ( s / ν ) g ( t / ν ) ( 1 2 H i s t + H i ( 2 H i 1 ) s 2 t 2 + O ( t 3 ) ) ] i = 1 n 2 a i 2 H i s ν 2 H i t 2 H i 1 .
Theorem 2.
Let Z H , a = { Z t H , a , t 0 } be a gmfBm of parameters a = ( a 1 , , a n ) and H = ( H 1 , , H n ) , H i ( 0 , 1 ) , a i R , n N * with H i < H n for i = 1 , 2 , , n 1 . Let Γ = { Γ t , t 0 } be a gamma process with parameter ν > 0 . Let T Γ H , a be the gmfBm time-changed by Γ. Then, the time-changed gmfBm by means of Γ has la ong-range dependence property for all Hurst indices satisfying 0 < 2 H i H n < 1 .
Proof. 
Let n N * . Let Z H , a = { Z t H , a , t 0 } be a gmfBm of parameters a = ( a 1 , , a n ) and H = ( H 1 , , H n ) , H i ( 0 , 1 ) , a i R with H i < H n for i = 1 , 2 , , n 1 . Let Γ = { Γ t , t 0 } be a gamma process with parameter ν > 0 . The process T Γ H , a is not stationary, hence the Definition 2 will be used to establish the long-range dependence property.
Using Equations (2), (7) and by Taylor’s expansion we obtain, as t
C o r r ( T Γ t H , a , T Γ s H , a ) i = 1 n 2 a i 2 H i s ν 2 H i t 2 H i 1 ( i = 1 n 2 a i 2 H i ν 2 H i t 2 H i ) E ( T Γ s H , a ) 2 = i = 1 n 2 a i 2 H i s ν 2 H i t 2 H i 1 2 a n 2 H n t 2 H n ν 2 H n [ i = 1 n 1 2 a i 2 H i ν 2 H i t 2 H i 2 H n + 1 ] E ( T Γ s H , a ) 2 i = 1 n 2 a i 2 H i s ν 2 H i t 2 H i 1 2 a n 2 H n ν 2 H n t 2 H n E ( T Γ s H , a ) 2 , s i n c e H i < H n = i = 1 n 1 2 a i 2 H i s ν 2 H i 2 a n 2 H n ν 2 H n E ( T Γ s H , a ) 2 t 2 H i H n 1 + | a n | ( 2 H n ) 1 2 ν H n s E ( T Γ s H , a ) 2 t H n 1 .
Then, for fixed s > 0 and t . For H i < H n < 1 for i = 1 , , n 1 , the correlation function is given by
C o r r ( T Γ t H , a , T Γ s H , a ) c 1 t 2 H i H n 1 + c 2 t H n 1 .
Therefore, for 0 < 2 H i H n < 1 the correlation function of T Γ t H , a decays like a t ( 1 H n ) for all 0 < H n < 1 . Then, in the sense of Definition 2 the time-changed process T S λ , α H , a exhibits the long-range dependence property for all 0 < 2 H i H n < 1 for i = 1 , 2 , , n 1 . □
Hence we obtain the following results, proved in [12].
Corollary 2.
The fractional Brownian motion time-changed by the gamma process has long-range dependence for all H ( 0 , 1 ) .

5. Conclusions

The time-changed gmfBm is a versatile and powerful tool for modeling natural phenomena that exhibit long-range dependence. The ability to control the long-range dependence property through the Hurst indices is a key feature of this model, and it allows us to tailor the model to the specific characteristics of the phenomenon we are trying to model. In this paper, it is shown that the time-changed gmfBm exhibits a long-range dependence property under some conditions on the Hurst indices when it is time-changed by a tempered stable subordinator or a gamma process. This is a significant result, as it shows that the time-changed gmfBm can be used to model a wide variety of natural phenomena that exhibit long-range dependence, even when the underlying process is not itself long-range dependent. We deduce that the fractional Brownian motion time-changed by tempered stable subordinator or gamma process has long-range dependence for all H ( 0 , 1 ) .

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The author thank the reviewers for the careful reading of the manuscript and their helpful comments.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Cheridito, P. Mixed fractional Brownian motion. Bernoulli 7 2001, 2, 913–934. [Google Scholar] [CrossRef]
  2. Alajmi, S.; Mliki, E. Mixed generalized fractional Brownian motion. J. Stoch. Anal. 2021, 2, 2–14. [Google Scholar]
  3. El-Nouty, C. The fractional mixed fractional Brownian motion. Stat. Prob. Lett. 2003, 65, 111–120. [Google Scholar]
  4. Majdoub, M.; Mliki, E. Well-posedness for Hardy-Hénon parabolic equations with fractional Brownian noise. Anal. Math. Phys. 2021, 11, 1–12. [Google Scholar]
  5. Miao, Y.; Ren, W.; Ren, Z. On the fractional mixed fractional Brownian motion. App. Math. Sci 2008, 35, 1729–1938. [Google Scholar]
  6. Thäle, C. Further remarks on mixed fractional Brownian motion. Appl. Math. Sci. 2009, 3, 1885–1901. [Google Scholar]
  7. Drakakis, K.; Radulovic, D. A discretized version of the self-similar model for internet traffic. Appl. Math. Sci. 2008, 2, 2743–2756. [Google Scholar]
  8. Chena, W.; Bowen, Y.; Guanghua, C.; Zhangd, Y. Numerically pricing American options under the generalized mixed fractional Brownian motion model. Phys. A Stat. Mech. Its Appl. 2016, 451, 180–189. [Google Scholar]
  9. Herry, P.S. Generalized mixed fractional Brownian motion as a generalized white noise functional. J. Mat. Stat. 2011, 11, 10–19. [Google Scholar]
  10. Herry, P.; Gunarso, B. Self-intersection local times of generalized mixed fractional Brownian motion as white noise distributions. J. Phys. Conf. Ser. 2017, 855, 012050. [Google Scholar]
  11. Kumar, A.; Gajda, J.; Wyłomańska, A.; Połoczański, R. Fractional Brownian motion delayed by tempered and inverse tempered stable subordinators. Methodol. Comput. Appl. Probab. 2019, 21, 185–202. [Google Scholar] [CrossRef] [Green Version]
  12. Kumar, A.; Wylomańska, A.; Polozański, R.; Sundar, S. Fractional Brownian motion time-changed by gamma and inverse gamma process. Pys. A Stat. Mech. Its Appl. 2017, 468, 648–667. [Google Scholar]
  13. Bochner, S. Diffusion equation and stochastic processes. Proc. Nat. Acad. Sci. USA 1949, 35, 368–370. [Google Scholar] [PubMed] [Green Version]
  14. Mliki, E. Correlation structure of time-changed fractional mixed fractional Brownian motion. 2023. submitted for publication. [Google Scholar]
  15. Önalan, Ö. Time-changed generalized mixed fractional Brownian motion and application to arithmetic average Asian option pricing. Int. J. Appl. Math. Res. 2017, 6, 85–92. [Google Scholar]
  16. Teuerle, M.; Wyłomańska, A.; Sikora, G. Modeling anomalous diffusion by a subordinated fractional Lévy-stable process. J. Stat. Mech. Theory Exp. 2013, 2013, 5–16. [Google Scholar]
  17. Alajmi, S.; Mliki, E. On the mixed fractional Brownian motion time changed by inverse α-stable subordinator. Appl. Math. Sci. 2020, 14, 755–763. [Google Scholar]
  18. Guo, Z.; Yuan, H. Pricing European option under the time-changed mixed Brownian-fractional Brownian model. Phys. A Stat. Mech. Its Appl. 2014, 406, 73–79. [Google Scholar] [CrossRef]
  19. Shokrollahi, F. The evaluation of geometric Asian power options under time changed mixed fractional Brownian motion. J. Comput. Appl. Math. 2018, 344, 716–724. [Google Scholar]
  20. Gu, H.; Liang, J.R.; Zhang, Y.X. The time changed geometric fractional Brownian motion and option pricing with transaction costs. Phys. A Stat. Mech. Its Appl. 2012, 391, 3971–3977. [Google Scholar]
  21. Kim, K.; Kim, S.; Jo, H. Option pricing under mixed hedging strategy in time-changed mixed fractional Brownian model. J. Comput. Appl. Math. 2022, 416, 73–79. [Google Scholar]
  22. Melnikov, A.; Mishura, Y. On pricing in financial markets with long-range dependence. Math. Financ. Econ. 2011, 5, 29–46. [Google Scholar]
  23. Miao, J. Option Pricing with Transaction Costs under the Subdiffusive Mixed Fractional Brownian Motion. J. Phys. Conf. Ser. 2020, 1670, 012045. [Google Scholar]
  24. Magdziarz, M. Stochastic path properties of sub-diffusion, a martingale approach. Stoch. Models 2010, 26, 256–271. [Google Scholar] [CrossRef]
  25. Mijena, J.B.; Nane, E. Correlation structure of time-changed pearson diffusions. Stat. Probab. Lett. 2014, 90, 68–77. [Google Scholar]
  26. Mliki, E. On the fractional mixed fractional Brownian motion time changed by inverse α-Stable subordinator. Glob. Stoch. Anal. 2023, 10, No 1. [Google Scholar]
  27. Ber, J. Subordinators: Theory and Applications; Cambridge University Press: Cambridge, UK, 1996. [Google Scholar]
  28. Bertoin, J. Lévy Processes of Cambridge Tracts in Mathematics; Cambridge University Press: Cambridge, UK, 1996; Volume 121. [Google Scholar]
  29. Cont, R.; Tankov, P. Financial Modeling with Jump Processes; CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  30. Maheshwari, A.; Vellaisamy, P. On the long-range dependence of fractional poisson and negative binomial processes. J. Appl. Probab. 2016, 53, 989–1000. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mliki, E. Correlation Structure of Time-Changed Generalized Mixed Fractional Brownian Motion. Fractal Fract. 2023, 7, 591. https://doi.org/10.3390/fractalfract7080591

AMA Style

Mliki E. Correlation Structure of Time-Changed Generalized Mixed Fractional Brownian Motion. Fractal and Fractional. 2023; 7(8):591. https://doi.org/10.3390/fractalfract7080591

Chicago/Turabian Style

Mliki, Ezzedine. 2023. "Correlation Structure of Time-Changed Generalized Mixed Fractional Brownian Motion" Fractal and Fractional 7, no. 8: 591. https://doi.org/10.3390/fractalfract7080591

APA Style

Mliki, E. (2023). Correlation Structure of Time-Changed Generalized Mixed Fractional Brownian Motion. Fractal and Fractional, 7(8), 591. https://doi.org/10.3390/fractalfract7080591

Article Metrics

Back to TopTop