Next Article in Journal
On Sharp Bounds of Local Fractional Metric Dimension for Certain Symmetrical Algebraic Structure Graphs
Next Article in Special Issue
A New Fractional Curvature of Curves and Surfaces in Euclidean Space Using the Caputo’s Fractional Derivative
Previous Article in Journal
An Adaptive Mesh Refinement–Rotated Lattice Boltzmann Flux Solver for Numerical Simulation of Two and Three-Dimensional Compressible Flows with Complex Shock Structures
Previous Article in Special Issue
Fractional Integrals Associated with the One-Dimensional Dunkl Operator in Generalized Lizorkin Space
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Existence and Averaging Principle for a Class of Fractional Hadamard Itô–Doob Stochastic Integral Equations

by
Mohamed Rhaima
1,
Lassaad Mchiri
2 and
Abdellatif Ben Makhlouf
3,*
1
Department of Statistics and Operations Research, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
2
Panthéon-Assas University Paris II, 92 Rue d’Assas, 75006 Paris, France
3
Department of Mathematics, Faculty of Sciences of Sfax, University of Sfax, Road of Soukra, Sfax 1171, Tunisia
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(10), 1910; https://doi.org/10.3390/sym15101910
Submission received: 25 August 2023 / Revised: 10 October 2023 / Accepted: 11 October 2023 / Published: 12 October 2023
(This article belongs to the Special Issue Applied Mathematics and Fractional Calculus III)

Abstract

:
In this paper, we investigate the existence and uniqueness properties pertaining to a class of fractional Hadamard Itô–Doob stochastic integral equations (FHIDSIE). Our study centers around the utilization of the Picard iteration technique (PIT), which not only establishes these fundamental properties but also unveils the remarkable averaging principle within FHIDSIE. To accomplish this, we harness powerful mathematical tools, including the Hölder and Gronwall inequalities.

1. Introduction

The theory of fractional calculus is an ancient subject that appeared in the same century as classical calculus. The question of fractional derivatives was raised as early as 1695 by Leibniz in a letter to l’Hospital, but when Leibniz asked what the derivative of a function means when the order is 1 2 , Leibniz replied that this leads to a paradox. This letter to l’Hospital was considered the first incident of fractional calculus theory. Later, Euler’s functions made the definition of fractional operators possible. However, a systematic development of the subject did not occur until the nineteenth century with Riemann, Liouville, Grünwald, and Letnikov (see [1,2,3]).
Fractional calculus has practical applications across a range of disciplines, including biology, machine learning, and physics. A notable example is found in the study of anomalous diffusion phenomena, where fractional variants of the diffusion equation serve as indispensable tools for investigation (as highlighted in references such as [4,5,6]). These applications underscore the versatility and potency of fractional calculus in unraveling intricate phenomena that conventional calculus may struggle to address effectively.
Caputo–Hadamard and Hadamard derivatives are crucial classes of fractional derivatives. The scientific community has extensively studied Hadamard fractional integral equations (HFIE); for further details, readers can consult [7,8,9]. The FHIDSIE represents a significant category within HFIE, drawing the attention of numerous researchers (see [10,11]). In [12], the authors analyzed the existence and uniqueness of solutions using PIT and the averaging principle of FHIDSIE.
Inspired by the analysis above, this paper investigates the averaging principle for FHIDSIE. In the literature, most published articles that study the averaging principle focus on ordinary stochastic differential equations. In this sense, this paper extends the results in [12] to the neutral delay FHIDSIE. Therefore, in this paper, our main aims are:
  • Analyze the existence and uniqueness of the solution of FHIDSIE by employing the PIT;
  • Investigate the averaging principle of FHIDSIE using the famous stochastic calculus techniques.
The structure of this article is summarized as follows. Section 2 is devoted to several definitions and classical results. In Section 3, we prove the existence and uniqueness result using the PIT. In Section 4, we present two theorems related to demonstrating the averaging principle of FHIDSIE. Section 5 is dedicated to showcasing an application related to the studied problem. Finally, we provide a conclusion in Section 6.

2. Basic Notions

Let D = { N , Z , ( Z ϰ ) ϰ 1 , P } be a complete probability space and W ( ϰ ) be a standard Brownian motion with respect the space D . Denote by C ( [ 1 ϕ , 1 ] ; R n ) the family of functions φ from [ 1 ϕ , 1 ] to R n that are right-continuous and have limits on the left. C ( [ 1 ϕ , 1 ] ; R n ) is equipped with the norm φ = sup 1 ϕ r 1 | φ ( r ) | and | y | = y T y for any y R n . Denote by L F ϰ 2 ( [ 1 ϕ , 1 ] ; R n ) , ϰ 1 the set of all F ϰ -measurable, C ( [ 1 ϕ , 1 ] ; R n ) -valued random variables φ = { φ ( τ ) : 1 ϕ τ 1 } satisfying E sup 1 ϕ τ 1 | φ ( τ ) | 2 < .
Definition 1. 
The Hadamard fractional integral of order ς for a function h is defined as
I ς h ( ϑ ) = 1 Γ ( ς ) 1 ϑ log ϑ ν ς 1 h ( ν ) ν d ν , ς > 0 ,
provided the integral exists.
For a more comprehensive and in-depth understanding of fractional calculus, we highly recommend consulting the authoritative work in [1,13,14], which offers a wealth of detailed insights and thorough explanations on this subject matter.
Let
H : C ( [ 1 ϕ , 1 ] , R n ) R n ,
b : [ 1 , T ] × C ( [ 1 ϕ , 1 ] , R n ) R n ,
σ 1 : [ 1 , T ] × C ( [ 1 ϕ , 1 ] , R n ) R n ,
σ 2 : [ 1 , T ] × C ( [ 1 ϕ , 1 ] , R n ) R n
all be Borel measurable. Consider the following FHIDSIE:
ζ ( ϰ ) H ( ζ ϰ ) = ζ ( 1 ) H ( ζ 1 ) + 1 ϰ b s , ζ s d s + 1 ϰ σ 1 s , ζ s d W ( s ) + ι 1 ϰ log ϰ s ι 1 σ 2 s , ζ s s d s
with initial condition
ζ 1 = α = { α ( λ ) ; 1 ϕ λ 1 } L Z 1 2 [ 1 ϕ , 1 ] , R n ,
1 2 < ι < 1 and ϕ > 0 .
Take the following hypotheses:
H 1
: There exists K ̲ > 0 that satisfies
| b ( ρ , θ ) b ( ρ , θ ˜ ) | 2 | σ 1 ( ρ , θ ) σ 1 ( ρ , θ ˜ ) | 2 | σ 2 ( ρ , α 1 ) σ 2 ( ρ , θ ˜ ) | 2 K ̲ | | θ θ ˜ | | 2 ,
for all ( ρ , θ , θ ˜ ) [ 1 , T ] × C ( [ 1 ϕ , 1 ] , R n ) × C ( [ 1 ϕ , 1 ] , R n ) .
H 2
: There exists K > 0 that satisfies
| b ( ρ , θ ) | 2 | σ 1 ( ρ , θ ) | 2 | σ 2 ( ρ , θ ) | 2 K ( 1 + | | θ | | 2 ) ,
for all ( ρ , θ ) [ 1 , T ] × C ( [ 1 ϕ , 1 ] , R n ) .
H 3
: Assume that H ( 0 ) = 0 and there is ϱ ̲ ( 0 , 1 ) such that
| H ( θ ) H ( θ ˜ ) | ϱ ̲ | | θ θ ˜ | | ,
for all ( θ , θ ˜ ) C ( [ 1 ϕ , 1 ] , R n ) × C ( [ 1 ϕ , 1 ] , R n ) .

3. Existence and Uniqueness of Solution

Theorem 1. 
Under assumptions H 1 H 3 , there exists a unique solution ζ ( ϰ ) of system (1) that satisfies ζ ( ϰ ) M 2 [ 1 ϕ , T ] ; R n .
Before demonstrating Theorem 1, we require the following lemma.
Lemma 1. 
Assume that assumptions H 2 and H 3 hold. If ζ ( ϰ ) is a solution to Equation (1) with initial condition (2), then
E sup 1 ϕ ϰ T | ζ ( ϰ ) | 2 L 1 ,
where
L 1 = 1 + 4 + ϱ ̲ ϱ ̲ ( 1 ϱ ̲ ) ( 1 ϱ ̲ ) E | | α | | 2 exp 3 K ( T 1 ) T + 3 + 3 ι 2 K ( log T ) 2 ι 1 2 ι 1 ( 1 ϱ ̲ ) ( 1 ϱ ̲ ) .
Moreover, ζ ( ϰ ) M 2 [ 1 ϕ , T ] ; R n .
Proof. 
For any n N * , denote by τ n = T inf { ϰ [ 1 , T ] ; | | ζ ϰ | | n } a stopping time. It is obvious that τ n T a.s. Let ζ n ( ϰ ) = ζ ( ϰ τ n ) , for ϰ [ 1 ϕ , T ] . Thus, for 1 ϰ T ,
ζ n ( ϰ ) = H ( ζ ϰ n ) H ( α ) + O n ( ϰ ) ,
where
O n ( ϰ ) = α ( 1 ) + 1 ϰ b s , ζ s n 1 [ 1 , τ n ] ( s ) d s + 1 ϰ σ 1 s , ζ s n 1 [ 1 , τ n ] ( s ) d W ( s )
+ ι 1 ϰ log ϰ τ n s ι 1 σ 2 s , ζ s n s 1 [ 1 , τ n ] ( s ) d s .
By Lemma 2.3 in [15] and assumption H 3 , we can derive that
| ζ n ( ϰ ) | 2 1 ϱ ̲ | H ( ζ ϰ n ) H ( α ) | 2 + 1 1 ϱ ̲ | O n ( ϰ ) | 2 ϱ ̲ | | ζ ϰ n α | | 2 + 1 1 ϱ ̲ | O n ( ϰ ) | 2 ϱ ̲ | | ζ ϰ n | | 2 + ϱ ̲ 1 ϱ ̲ | | α | | 2 + 1 1 ϱ ̲ | O n ( ϰ ) | 2 .
Consequently,
E sup 1 s ϰ | ζ n ( s ) | 2 ϱ ̲ E sup 1 ϕ s ϰ | ζ n ( s ) | 2 + ϱ ̲ 1 ϱ ̲ E | | α | | 2 + 1 1 ϱ ̲ E sup 1 s ϰ | O n ( s ) | 2 .
It is not hard to see that sup 1 ϕ s ϰ | ζ n ( s ) | 2 | | α | | 2 + sup 1 s ϰ | ζ n ( s ) | 2 . Then,
E sup 1 ϕ s ϰ | ζ n ( s ) | 2 E | | α | | 2 + E sup 1 s ϰ | ζ n ( s ) | 2 .
Therefore,
E sup 1 ϕ s ϰ | ζ n ( s ) | 2 ϱ ̲ E sup 1 ϕ s ϰ | ζ n ( s ) | 2 + 1 + ϱ ̲ ϱ ̲ 1 ϱ ̲ E | | α | | 2 + 1 1 ϱ ̲ E sup 1 s ϰ | O n ( s ) | 2 .
Hence,
E sup 1 ϕ s ϰ | ζ n ( s ) | 2 1 + ϱ ̲ ϱ ̲ ( 1 ϱ ̲ ) 2 E | | α | | 2 + 1 ( 1 ϱ ̲ ) ( 1 ϱ ̲ ) E sup 1 s ϰ | O n ( s ) | 2 .
Using the Cauchy–Schwartz inequality and assumption H 2 , we obtain
| O n ( ϰ ) | 2 3 | α ( 1 ) | 2 + 3 ( ϰ 1 ) 1 ϰ | b s , ζ s n | 2 d s + 3 | 1 ϰ σ 1 s , ζ s n 1 [ 1 , τ n ] ( s ) d W ( s ) | 2 + 3 ι 2 1 ϰ 1 s 2 log ϰ s 2 ι 2 d s 1 ϰ | σ 2 s , ζ s n | 2 d s 3 | α ( 1 ) | 2 + 3 K ( ϰ 1 ) 1 ϰ ( 1 + | | ζ s n | | ) 2 d s + 3 | 1 ϰ σ 1 s , ζ s n 1 [ 1 , τ n ] ( s ) d W ( s ) | 2 + 3 ι 2 K 2 ι 1 ( log ϰ ) 2 ι 1 1 ϰ ( 1 + | | ζ s n | | ) 2 d s .
Therefore, applying Theorem 7.2 (p. 40 in [15]) and hypothesis H 2 , one can obtain
E sup 1 s ϰ | O n ( s ) | 2 3 E | | α | | 2 + 3 K ( T 1 ) 1 ϰ ( 1 + E | | ζ s n | | 2 ) d s + 12 E 1 ϰ | σ 1 s , ζ s n | 2 1 [ 1 , τ n ] ( s ) d s + 3 ι 2 K 2 ι 1 ( log T ) 2 ι 1 1 ϰ ( 1 + E | | ζ s n | | 2 ) d s 3 E | | α | | 2 + 3 K T + 3 + 3 ι 2 K ( log T ) 2 ι 1 2 ι 1 1 ϰ ( 1 + E | | ζ s n | | 2 ) d s .
Plugging (12) into (10), we can obtain
E sup 1 ϕ s ϰ | ζ n ( s ) | 2 4 + ϱ ̲ ϱ ̲ ( 1 ϱ ̲ ) ( 1 ϱ ̲ ) E | | α | | 2 + 3 K T + 3 + 3 ι 2 K ( log T ) 2 ι 1 2 ι 1 ( 1 ϱ ̲ ) ( 1 ϱ ̲ ) 1 ϰ ( 1 + E | | ζ s n | | 2 ) d s .
Hence,
1 + E sup 1 ϕ s ϰ | ζ n ( s ) | 2 1 + 4 + ϱ ̲ ϱ ̲ ( 1 ϱ ̲ ) ( 1 ϱ ̲ ) E | | α | | 2 + 3 K T + 3 + 3 ι 2 K ( log T ) 2 ι 1 2 ι 1 ( 1 ϱ ̲ ) ( 1 ϱ ̲ ) 1 ϰ 1 + E sup 1 ϕ l s | ζ n ( l ) | 2 d s .
By Gronwall’s inequality, one has
1 + E sup 1 ϕ s ϰ | ζ n ( s ) | 2 1 + 4 + ϱ ̲ ϱ ̲ ( 1 ϱ ̲ ) ( 1 ϱ ̲ ) E | | α | | 2 exp 3 K ( T 1 ) T + 3 + 3 ι 2 K ( log T ) 2 ι 1 2 ι 1 ( 1 ϱ ̲ ) ( 1 ϱ ̲ ) .
Consequently,
E sup 1 ϕ s ϰ | ζ n ( s ) | 2 1 + 4 + ϱ ̲ ϱ ̲ ( 1 ϱ ̲ ) ( 1 ϱ ̲ ) E | | α | | 2 exp 3 K ( T 1 ) T + 3 + 3 ι 2 K ( log T ) 2 ι 1 2 ι 1 ( 1 ϱ ̲ ) ( 1 ϱ ̲ ) .
Therefore, we can deduce the desired result when n . □
Proof of Theorem 1:
Proof. 
Uniqueness:
Consider ζ ( ϰ ) and ζ ¯ ( ϰ ) to be two solutions to Equation (1). Using Lemma 1, ζ ( ϰ ) and ζ ¯ ( ϰ ) belong to M 2 [ 1 ϕ , T ] ; R n . Thus,
ζ ( ϰ ) ζ ¯ ( ϰ ) = H ( ζ ϰ ) H ( ζ ¯ ϰ ) + O ( ϰ ) ,
where
O ( ϰ ) = 1 ϰ b s , ζ s b s , ζ ¯ s ) d s + 1 ϰ σ 1 s , ζ s σ 1 s , ζ ¯ s ) d W ( s )
+ ι 1 ϰ log ϰ s ι 1 σ 2 s , ζ s σ 2 s , ζ ¯ s ) s d s .
By Lemma 2.3 in [15] and assumption H 3 , we can see that
| ζ ( ϰ ) ζ ¯ ( ϰ ) | 2 ϱ ̲ | | ζ ϰ ζ ¯ ϰ | | 2 + 1 1 ϱ ̲ | O ( ϰ ) | 2 .
Hence,
E sup 1 s ϰ | ζ ( s ) ζ ¯ ( s ) | 2 ϱ ̲ E sup 1 s ϰ | ζ ( s ) ζ ¯ ( s ) | 2 + 1 1 ϱ ̲ E sup 1 s ϰ | O ( s ) | 2 .
Using (15), it yields that
E sup 1 s ϰ | ζ ( s ) ζ ¯ ( s ) | 2 1 ( 1 ϱ ̲ ) 2 E sup 1 s ϰ | O ( s ) | 2 .
Moreover, by Hölder’s inequality, Theorem 7.2 (p. 40 in [15]), and H 1 , using a similar method as in the proof of Lemma 1, we obtain
E sup 1 s ϰ | O ( s ) | 2 3 K ̲ T + 4 + ι 2 ( log T ) 2 ι 1 2 ι 1 1 ϰ E | | ζ s ζ ¯ s | | 2 d s 3 K ̲ T + 4 + ι 2 ( log T ) 2 ι 1 2 ι 1 1 ϰ E sup 1 w s | ζ ( w ) ζ ¯ ( w ) | 2 d s .
Therefore,
E sup 1 s ϰ | ζ ( s ) ζ ¯ ( s ) | 2 3 K ̲ T + 4 + ι 2 ( log T ) 2 ι 1 2 ι 1 ( 1 ϱ ̲ ) 2 1 ϰ E sup 1 w s | ζ ( w ) ζ ¯ ( w ) | 2 d s .
Using the Gronwall inequality, we can derive that: E sup 1 ϰ T | ζ ( ϰ ) ζ ¯ ( ϰ ) | 2 = 0 . Consequently, for all 1 ϕ ϰ T , ζ ( ϰ ) = ζ ¯ ( ϰ ) , a.s.
Existence: Let λ be sufficiently small such that
ϑ = ϱ ̲ + 3 K ̲ λ + 5 + ι 2 ( log ( λ + 1 ) ) 2 ι 1 2 ι 1 λ 1 ϱ ̲ < 1 ,
and write [ 1 , T ] = k = 1 r [ 1 + ( k 1 ) λ , 1 + k λ ] [ 1 + r λ , T ] . We split the proof into two steps:
Step 1: Let ζ 1 0 = α and ζ 0 ( ϰ ) = α ( 1 ) for 1 ϰ 1 + λ . For each n = 1 , 2 , 3 , , let ζ 1 n = α , and we define, using the Picard iterations:
ζ n ( ϰ ) H ( ζ ϰ n 1 ) = α ( 1 ) H ( α ) + 1 ϰ b s , ζ s n 1 d s + 1 ϰ σ 1 s , ζ s n 1 d W ( s ) + ι 1 ϰ log ϰ s ι 1 σ 2 s , ζ s n 1 s d s .
It is not hard to see that ζ n ( · ) M 2 [ 1 ϕ , 1 + λ ] ; R n . Noting that for 1 ϰ 1 + λ ,
ζ 1 ( ϰ ) ζ 0 ( ϰ ) = ζ 1 ( ϰ ) α ( 1 ) = H ( ζ ϰ 0 ) H ( α ) + 1 ϰ b s , ζ s 0 d s + 1 ϰ σ 1 s , ζ s 0 d W ( s ) + ι 1 ϰ log ϰ s ι 1 σ 2 s , ζ s 0 s d s .
Proceeding as in the proof of the uniqueness, we have
E sup 1 s ϰ | ζ 1 ( ϰ ) ζ 0 ( ϰ ) | 2 ϱ ̲ E sup 1 s ϰ | | ζ ϰ 0 α | | 2 + 3 K T + 4 + ι 2 ( log T ) 2 ι 1 2 ι 1 1 ϱ ̲ E 1 ϰ 1 + | | ζ ϰ 0 | | 2 d ϰ C 1
where
C 1 = 2 ϱ ̲ E | | α | | 2 + 3 K λ + 5 + ι 2 ( log ( λ + 1 ) ) 2 ι 1 2 ι 1 1 ϱ ̲ 1 + E | | α | | 2 λ .
Noting that for n 1 and 1 ϰ 1 + λ , we have
ζ n + 1 ( ϰ ) ζ n ( ϰ ) = H ( ζ ϰ n ) H ( ζ ϰ n 1 ) + 1 ϰ b s , ζ s n b s , ζ s n 1 d s + 1 ϰ σ 1 s , ζ s n σ 1 s , ζ s n 1 d W ( s ) + ι 1 ϰ log ϰ s ι 1 σ 2 s , ζ s n σ 2 s , ζ s n 1 s d s .
Proceeding as in the proof of the uniqueness and using (18), we can obtain
E sup 1 ϰ 1 + λ | ζ n + 1 ( ϰ ) ζ n ( ϰ ) | 2 ϱ ̲ E sup 1 ϰ 1 + λ | ζ n ( ϰ ) ζ n 1 ( ϰ ) | 2 + 3 K ̲ λ + 5 + ι 2 ( log ( λ + 1 ) ) 2 ι 1 2 ι 1 1 ϱ ̲ 1 ϰ E sup 1 s ϰ | ζ n ( s ) ζ n 1 ( s ) | 2 d ϰ ϑ E sup 1 ϰ 1 + λ | ζ n ( ϰ ) ζ n 1 ( ϰ ) | 2 ϑ n E sup 1 ϰ 1 + λ | ζ 1 ( ϰ ) ζ 0 ( ϰ ) | 2 C 1 ϑ n .
By (16), we can prove from (19) that there exists a solution 1 ζ to Equation (1) on [ 1 ϕ , 1 + λ ] in the same way as in the proof of Theorem 3.1 in [12].
Step 2: For 2 k r and ϰ [ 1 + ( k 1 ) λ , 1 + k λ ] : let ζ 0 ( ϰ ) = k 1 ζ ( ϰ ) for 1 ϕ ϰ 1 + ( k 1 ) λ and ζ 0 ( ϰ ) = k 1 ζ ( 1 + ( k 1 ) λ ) for 1 + ( k 1 ) λ ϰ 1 + k λ . For each n = 1 , 2 , 3 , , let ζ n ( ϰ ) = k 1 ζ ( ϰ ) for ϰ [ 1 ϕ , 1 + ( k 1 ) λ ] . We define the following Picard iterations:
ζ n ( ϰ ) H ( ζ ϰ n 1 ) = α ( 1 ) H ( α ) + 1 1 + ( k 1 ) λ b s , k 1 ζ s d s + 1 1 + ( k 1 ) λ σ 1 s , k 1 ζ s d W ( s ) + ι 1 1 + ( k 1 ) λ log ϰ s ι 1 σ 2 s , k 1 ζ s s d s + 1 + ( k 1 ) λ ϰ b s , ζ s n 1 d s + 1 + ( k 1 ) λ ϰ σ 1 s , ζ s n 1 d W ( s ) + ι 1 + ( k 1 ) λ ϰ log ϰ s ι 1 σ 2 s , ζ s n 1 s d s .
Proceeding as in step 1, we obtain the existence of the solution k ζ on [ 1 ϕ , 1 + k λ ] .
Repeating this on [ 1 + r λ , T ] , we can observe that there exists a solution to Equation (1) on the entire interval [ 1 ϕ , T ] , which complete the proof. □

4. Averaging Principle

Take the following standard form of Equation (1):
ζ ϵ ( ϰ ) = ζ ( 1 ) H ( ζ 1 ) + H ( ζ ϵ , ϰ ) + ϵ 1 ϰ b s , ζ ϵ , s d s + ϵ 1 ϰ σ 1 s , ζ ϵ , s d W ( s ) + ϵ ι ι 1 ϰ log ϰ s ι 1 σ 2 s , ζ ϵ , s s d s ,
where b , σ 1 , and σ 2 satisfy H 1 and H 2 , and ϵ ( 0 , ϵ 0 ] , with ϵ 0 ( 0 , 1 2 ) , is a fixed number. Using Theorem 1, Equation (21) has a unique global solution ζ ϵ ( ϰ ) , ϰ [ 1 ϕ , T ] , for any ϵ ( 0 , ϵ 0 ] .
Suppose that the following hypothesis holds true:
H 4
: Let the measurable functions b ¯ ( ζ ) : C ( [ 1 ϕ , 1 ] , R n ) R n , σ 1 ¯ ( ζ ) : C ( [ 1 ϕ , 1 ] , R n ) R n and σ 2 ¯ ( ζ ) : C ( [ 1 ϕ , 1 ] , R n ) R n exist, satisfying H 1 and H 2 . For any ( ζ , T 1 ) C ( [ 1 ϕ , 1 ] , R n ) × [ 1 , T ] , we have
  • 1 T 1 1 T 1 | b ( s , ζ ) b ¯ ( ζ ) | 2 d s ψ 1 ( T 1 ) ( 1 + | | ζ | | 2 ) ;
  • 1 T 1 1 T 1 | σ 1 ( s , ζ ) σ 1 ¯ ( ζ ) | 2 d s ψ 2 ( T 1 ) ( 1 + | | ζ | | 2 ) ;
  • 1 T 1 1 T 1 | σ 2 ( s , ζ ) σ 2 ¯ ( ζ ) | 2 d s ψ 3 ( T 1 ) ( 1 + | | ζ | | 2 ) ;
with ψ j ( T 1 ) as positive bounded functions satisfying lim x ψ j ( x ) = 0 for j = 1 , 2 , 3 .
Our objective is to establish that the solution ζ ϵ ( ϰ ) can be approximated by the solution y ϵ ( ϰ ) of the next equation
y ϵ ( ϰ ) = ζ ( 1 ) H ( ζ 1 ) + H ( y ϵ , ϰ ) + ϵ 1 ϰ b ¯ y ϵ , s d s + ϵ 1 ϰ σ ¯ 1 y ϵ , s d W ( s ) + ϵ ι ι 1 ϰ log ϰ s ι 1 σ ¯ 2 y ϵ , s s d s ,
for ϰ [ 1 , T ] .
Now, we will show the first main theorem in this section.
Theorem 2. 
Assume that H 1 H 4 hold. For a given arbitrary small constant π 1 > 0 , there are two constants L > 0 and ϵ 1 ( 0 , ϵ 0 ] that satisfy ϵ ( 0 , ϵ 1 ] ,
E sup ϰ [ 1 , L ϵ ] | ζ ϵ ( ϰ ) y ϵ ( ϰ ) | 2 π 1 .
Proof. 
Using Lemma 2.3 in [15] and assumption H 3 , one can derive
E sup ϰ [ 1 , u ] | ζ ϵ ( ϰ ) y ϵ ( ϰ ) | 2 E sup ϰ [ 1 , u ] | ζ ϵ ( ϰ ) y ϵ ( ϰ ) H ( ζ ϵ , ϰ ) H ( y ϵ , ϰ ) | 2 ( 1 ϱ ̲ ) 2 .
Note that
ζ ϵ ( ϰ ) y ϵ ( ϰ ) H ( ζ ϵ , ϰ ) H ( y ϵ , ϰ ) = ϵ 1 ϰ b s , ζ ϵ , s b ¯ y ϵ , s d s + ϵ 1 ϰ σ 1 s , ζ ϵ , s σ ¯ 1 y ϵ , s d W ( s ) + ι ϵ ι 1 ϰ log ϰ s ι 1 σ 2 s , ζ ϵ , s σ ¯ 2 y ϵ , s s d s .
According to the elementary inequality, one has
E sup ϰ [ 1 , u ] | ζ ϵ ( ϰ ) y ϵ ( ϰ ) H ( ζ ϵ , ϰ ) H ( y ϵ , ϰ ) | 2 ϵ 2 E sup ϰ [ 1 , u ] | 1 ϰ b s , ζ ϵ , s b ¯ y ϵ , s d s | 2 + 3 ϵ E sup ϰ [ 1 , u ] | 1 ϰ σ 1 s , ζ ϵ , s σ ¯ 1 y ϵ , s d W ( s ) | 2 + 3 ι 2 ϵ 2 ι E sup ϰ [ 1 , u ] | 1 ϰ log ϰ s ι 1 σ 2 s , ζ ϵ , s σ ¯ 2 y ϵ , s s d s | 2
= J 1 + J 2 + J 3 ,
where
J 1 = 3 ϵ 2 E sup ϰ [ 1 , u ] | 1 ϰ b s , ζ ϵ , s b ¯ y ϵ , s d s | 2 6 ϵ 2 E sup ϰ [ 1 , u ] | 1 ϰ b s , ζ ϵ , s b s , y ϵ , s d s | 2 + 6 ϵ 2 E sup ϰ [ 1 , u ] | 1 ϰ b s , y ϵ , s b ¯ y ϵ , s d s | 2 6 ϵ 2 J 11 + 6 ϵ 2 J 12 .
By the Hölder inequality and H 1 , one has
J 11 = E sup ϰ [ 1 , u ] | 1 ϰ b s , ζ ϵ , s b s , y ϵ , s d s | 2 ( u 1 ) E 1 u | b s , ζ ϵ , s b s , y ϵ , s | 2 d s K ̲ u E 1 u | | ζ ϵ , s y ϵ , s | | 2 d s K ̲ u 1 u E sup v [ 1 , s ] | ζ ϵ ( v ) y ϵ ( v ) | 2 d s .
By H 4 , one can obtain
J 12 = E sup ϰ [ 1 , u ] | 1 ϰ b s , y ϵ , s b ¯ y ϵ , s d s | 2 u E sup ϰ [ 1 , u ] 1 ϰ | b s , y ϵ , s b ¯ y ϵ , s | 2 d s u 2 E sup ϰ [ 1 , u ] 1 ϰ 1 ϰ | b s , y ϵ , s b ¯ y ϵ , s | 2 d s .
As in the proof of Theorem 1 in [16], one can derive
J 12 u 2 sup ϰ [ 1 , u ] ψ 1 ( ϰ ) 1 + E sup ϰ [ 1 ϕ , u ] | y ϵ ( ϰ ) | 2 .
Therefore,
J 1 6 ϵ 2 K ̲ u 1 u E sup v [ 1 , s ] | ζ ϵ ( v ) y ϵ ( v ) | 2 d s + 6 ϵ 2 u 2 sup ϰ [ 1 , u ] ψ 1 ( ϰ ) 1 + E sup ϰ [ 1 ϕ , u ] | y ϵ ( ϰ ) | 2 .
J 2 = 3 ϵ E sup ϰ [ 1 , u ] | 1 ϰ σ 1 s , ζ ϵ , s σ ¯ 1 y ϵ , s d W ( s ) | 2 = 3 ϵ E sup ϰ [ 1 , u ] | 1 ϰ σ 1 s , ζ ϵ , s σ 1 s , y ϵ , s + σ 1 s , y ϵ , s σ ¯ 1 y ϵ , s d W ( s ) | 2 6 ϵ E sup ϰ [ 1 , u ] | 1 ϰ σ 1 s , ζ ϵ , s σ 1 s , y ϵ , s d W ( s ) | 2 + 6 ϵ E sup ϰ [ 1 , u ] | 1 ϰ σ 1 s , y ϵ , s σ ¯ 1 y ϵ , s d W ( s ) | 2 6 ϵ J 21 + 6 ϵ J 22 .
By employing Theorem 7.2 (p. 40 in [15]) and H 1 , one can obtain
J 21 = E sup ϰ [ 1 , u ] | 1 ϰ σ 1 s , ζ ϵ , s σ 1 s , y ϵ , s d W ( s ) | 2 4 E 1 u | σ 1 s , ζ ϵ , s σ 1 s , y ϵ , s | 2 d s 4 K ̲ E 1 u | | ζ ϵ , s y ϵ , s | | 2 d s 4 K ̲ 1 u E sup v [ 1 , s ] | ζ ϵ ( v ) y ϵ ( v ) | 2 d s .
Using H 4 , one can derive
J 22 = E sup ϰ [ 1 , u ] | 1 ϰ σ 1 s , y ϵ , s σ ¯ 1 y ϵ , s d W ( s ) | 2 4 E 1 u | σ 1 s , y ϵ , s σ ¯ 1 y ϵ , s | 2 d s 4 u E 1 u 1 u | σ 1 s , y ϵ , s σ ¯ 1 y ϵ , s | 2 d s 4 u sup ϰ [ 1 , u ] ψ 2 ( ϰ ) 1 + E sup ϰ [ 1 ϕ , u ] | y ϵ ( ϰ ) | 2 .
Consequently,
J 2 24 ϵ K ̲ 1 u E sup v [ 1 , s ] | ζ ϵ ( v ) y ϵ ( v ) | 2 d s + 24 ϵ u sup ϰ [ 1 , u ] ψ 2 ( ϰ ) 1 + E sup ϰ [ 1 ϕ , u ] | y ϵ ( ϰ ) | 2 .
Moreover,
J 3 = 3 ι 2 ϵ 2 ι E sup ϰ [ 1 , u ] | 1 ϰ log ϰ s ι 1 σ 2 s , ζ ϵ , s σ ¯ 2 y ϵ , s s d s | 2 6 ι 2 ϵ 2 ι E sup ϰ [ 1 , u ] | 1 ϰ log ϰ s ι 1 σ 2 s , ζ ϵ , s σ 2 s , y ϵ , s s d s | 2 + 6 ι 2 ϵ 2 ι E sup ϰ [ 1 , u ] | 1 ϰ log ϰ s ι 1 σ 2 s , y ϵ , s σ ¯ 2 y ϵ , s s d s | 2 6 ι 2 ϵ 2 ι J 31 + 6 ι 2 ϵ 2 ι J 32 .
By Hölder’s inequality and H 1 , one can derive
J 31 = E sup ϰ [ 1 , u ] | 1 ϰ log ϰ s ι 1 σ 2 s , ζ ϵ , s σ 2 s , y ϵ , s s d s | 2 sup ϰ [ 1 , u ] 1 ϰ log ϰ s 2 ι 2 1 s d s E 1 u | σ 2 s , ζ ϵ , s σ 2 s , y ϵ , s | 2 d s 1 2 ι 1 ( log u ) 2 ι 1 K ̲ E 1 u | | ζ ϵ , s y ϵ , s | | 2 d s 1 2 ι 1 ( log u ) 2 ι 1 K ̲ 1 u E sup v [ 1 , s ] | ζ ϵ ( v ) y ϵ ( v ) | 2 d s .
Using Hölder’s inequality and H 4 , one has
J 32 = E sup ϰ [ 1 , u ] | 1 ϰ log ϰ s ι 1 σ 2 s , y ϵ , s σ ¯ 2 y ϵ , s s d s | 2 u 2 ι 1 ( log u ) 2 ι 1 E 1 u 1 u | σ 2 s , y ϵ , s σ ¯ 2 y ϵ , s | 2 d s u 2 ι 1 ( log u ) 2 ι 1 ψ 3 ( u ) 1 + E sup ϰ [ 1 ϕ , u ] | y ϵ ( ϰ ) | 2 .
Therefore,
J 3 6 ι 2 ϵ 2 ι 2 ι 1 ( log u ) 2 ι 1 K ̲ 1 u E sup v [ 1 , s ] | ζ ϵ ( v ) y ϵ ( v ) | 2 d s + 6 ι 2 ϵ 2 ι u 2 ι 1 ( log u ) 2 ι 1 ψ 3 ( u ) 1 + E sup ϰ [ 1 ϕ , u ] | y ϵ ( ϰ ) | 2 .
Plugging (31)–(39) into (25), using Lemma 1 and (23), one can derive that
E sup ϰ [ 1 , u ] | ζ ϵ ( ϰ ) y ϵ ( ϰ ) | 2
1 ( 1 ϱ ̲ ) 2 6 ϵ K ̲ ϵ u + 4 + ι 2 ϵ 2 ι 1 u 2 ι 1 ( log u ) 2 ι 1 1 u E sup v [ 1 , s ] | ζ ϵ ( v ) y ϵ ( v ) | 2 d s + 6 ϵ u ϵ u sup l [ 1 , u ] ψ 1 ( l ) + 4 sup l [ 1 , u ] ψ 2 ( l ) + ι 2 ϵ 2 ι 1 u 2 ι 1 ( log u ) 2 ι 1 ψ 3 ( u ) × 1 + E sup ϰ [ 1 ϕ , u ] | y ϵ ( ϰ ) | 2 × 1 ( 1 ϱ ̲ ) 2 ϵ L 3 + ϵ L 2 1 u E sup v [ 1 , s ] | ζ ϵ ( v ) y ϵ ( v ) | 2 d s ,
with
L 2 = 6 K ̲ ϵ T + 4 + ι 2 ϵ 2 ι 1 u 2 ι 1 ( log T ) 2 ι 1 ( 1 ϱ ̲ ) 2 ,
and
L 3 = 6 T ϵ T sup l [ 1 , T ] ψ 1 ( l ) + 4 sup l [ 1 , T ] ψ 2 ( l ) + ι 2 ϵ 2 ι 1 2 ι 1 ( log T ) 2 ι 1 sup l [ 1 , T ] ψ 3 ( l ) 1 + L 1 ( 1 ϱ ̲ ) 2 .
By Gronwall’s inequality, one can obtain
E sup ϰ [ 1 , u ] | ζ ϵ ( ϰ ) y ϵ ( ϰ ) | 2 ϵ L 3 e ϵ L 2 u .
Therefore, given any number π 1 > 0 , there are two constants L > 0 and ϵ 1 ( 0 , ϵ 0 ] that satisfy, for every ϵ ( 0 , ϵ 1 ] ,
E sup ϰ [ 1 , L ϵ ] | ζ ϵ ( ϰ ) y ϵ ( ϰ ) | 2 π 1 .
Next, we show the convergence in probability between ζ ϵ ( ϰ ) and y ϵ ( ϰ ) .
Theorem 3. 
Suppose that the FHIDSIE (21) and (22) both satisfy H 1 H 4 . For a given arbitrary small constant π 1 > 0 , there are two constants L > 0 and ϵ 1 ( 0 , ϵ 0 ] that satisfy ϵ ( 0 , ϵ 1 ] ,
lim ϵ 0 P sup ϰ [ 1 , L ϵ ] | ζ ϵ ( ϰ ) y ϵ ( ϰ ) | > π 1 = 0 .
Proof. 
Given any constant π 1 > 0 , using Theorem 4.1 and the inequality of Chebyshev, we can obtain
P sup ϰ [ 1 , L ϵ ] | ζ ϵ ( ϰ ) y ϵ ( ϰ ) | > π 1 1 π 1 2 E sup ϰ [ 1 , L ϵ ] | ζ ϵ ( ϰ ) y ϵ ( ϰ ) | 2 1 π 1 2 ϵ L 3 e ϵ L 2 T .
Letting ϵ 0 , (42) holds. □

5. Application

In this section, we will present an application that illustrates the key findings outlined in the preceding section.
Consider the following perturbed Malthusian FHIDSIE model of population growth (see [15,17]):
ζ ϵ ( ϰ ) = ζ ( 1 ) H ( ζ 1 ) + H ( ζ ϵ , ϰ ) + ϵ 1 ϰ b s , ζ ϵ , s d s + ϵ 1 ϰ σ 1 s , ζ ϵ , s d W ( s ) + ϵ ι ι 1 ϰ log ϰ s ι 1 σ 2 s , ζ ϵ , s s d s ,
with initial condition
ζ 1 = α = { α ( λ ) ; 1 ϕ λ 1 } L Z 1 2 [ 1 ϕ , 1 ] , R n ,
where
H ( ζ ϵ , ϰ ) = 0.1 ζ ϵ ( ϕ ) , b s , ζ ϵ , s = c 1 ζ ϵ ( ϕ ) , σ 1 s , ζ ϵ , s = c 2 ζ ϵ ( ϕ ) , σ 2 s , ζ ϵ , s = c 3 cos 2 s ζ ϵ ( ϕ ) .
The assumptions H 1 , H 2 , and H 3 are satisfied for K ̲ = K = c 1 2 c 2 2 c 3 2 and ϱ ̲ = 0.1 .
We obtain the corresponding averaged equation
y ϵ ( ϰ ) = ζ ( 1 ) H ( ζ 1 ) + H ( y ϵ , ϰ ) + ϵ 1 ϰ b ¯ y ϵ , s d s + ϵ 1 ϰ σ ¯ 1 y ϵ , s d W ( s ) + ϵ ι ι 1 ϰ log ϰ s ι 1 σ ¯ 2 y ϵ , s s d s ,
where
b ¯ y ϵ , s = c 1 y ϵ ( ϕ ) , σ ¯ 1 y ϵ , s = c 2 y ϵ ( ϕ ) , σ ¯ 2 y ϵ , s = 0.5 c 3 y ϵ ( ϕ ) .
Now, we present a numerical simulation of Equations (44) and (46), initialized with the condition α ( ϰ ) = exp ( ϰ 2 / 2 ) , while employing the selected parameters ϕ = 0.5 and ϵ = 0.04 . Subsequently, in Figure 1 and Figure 2, we offer a comparative analysis between the precise solution ζ ϵ for Equation (44) and the averaged solution y ϵ for Equation (46), each evaluated for two distinct values of ι .

6. Conclusions

In this work, we delve into the intricacies surrounding the existence and uniqueness of FHIDSIE with an order denoted by ι , where ι belongs to the interval ( 1 2 , 1 ) . Our investigation in this paper harnesses the power of the PIT method to unravel these fundamental aspects of FHIDSIE. Moreover, we establish the profound averaging principle associated with FHIDSIE, leveraging moment inequalities to provide a rigorous and comprehensive proof.

Author Contributions

Conceptualization, L.M.; methodology, M.R.; writing—original draft preparation, A.B.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by “Researchers Supporting Project number (RSPD2023R683), King Saud University, Riyadh, Saudi Arabia”.

Data Availability Statement

No data were used to support this study.

Acknowledgments

The authors extend their appreciation to King Saud University in Riyadh, Saudi Arabia for funding this research work through Researchers Supporting Project number (RSPD2023R683).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  2. Podlubny, I. Fractional Differential Equations; Academic Press: Cambridge, MA, USA, 1999. [Google Scholar]
  3. Ceng, L.C.; Cho, S.Y. On approximate controllability for systems of fractional evolution hemivariational inequalities with Riemann-Liouville fractional derivatives. J. Nonlinear Var. Anal. 2022, 6, 421–438. [Google Scholar]
  4. Atanackovic, T.M.; Pilipovic, S.; Stankovic, B.; Zorica, D. Fractional Calculus with Applications in Mechanics; Wiley-ISTE: Indianapolis, IN, USA, 2014. [Google Scholar]
  5. Baleanu, D.; Machado, J.A.; Luo, A.C. Fractional Dynamics and Control; Springer Science and Business Media: New York, NY, USA, 2011. [Google Scholar]
  6. Pedjeu, J.C.; Ladde, G.S. Stochastic fractional differential equations: Modeling, method and analysis. Chaos Solitons Fractals 2012, 45, 279–293. [Google Scholar] [CrossRef]
  7. Abbas, S.; Benchohra, M.; Lazreg, J.E.; Zhou, Y. A survey on Hadamard and Hilfer fractional differential equations: Analysis and stability. Chaos Solitons Fractals 2017, 102, 47–71. [Google Scholar] [CrossRef]
  8. Almeida, R. Caputo-Hadamard fractional derivatives of variable order. Numer. Funct. Anal. Optim. 2017, 38, 1–19. [Google Scholar] [CrossRef]
  9. Makhlouf, A.B.; Mchiri, L. Some results on the study of Caputo-Hadamard fractional stochastic differential equations. Chaos Solitons Fractals 2022, 155, 111757. [Google Scholar] [CrossRef]
  10. Abouagwa, M.; Liu, J.; Li, J. Caratheodory approximations and stability of solutions to non-Lipschitz stochastic fractional differential equations of Itô-Doob type. Appl. Math. Comput. 2018, 329, 143–153. [Google Scholar] [CrossRef]
  11. Wang, W.; Guo, Z. Optimal index and averaging principle for Itô–Doob stochastic fractional differential equations. Stoch. Dyn. 2022, 22, 2250018. [Google Scholar] [CrossRef]
  12. Ben Makhlouf, A.; Mchiri, L.; Arfaoui, H.; Dhahri, S.; El-Hady, E.S.; Cherif, B. Hadamard Itô-Doob Stochastic Fractional Order Systems. Discret. Contin. Dyn. Syst.-S 2023, 16, 2060–2074. [Google Scholar] [CrossRef]
  13. Hattaf, K. A New Class of Generalized Fractal and Fractal-Fractional Derivatives with Non-Singular Kernels. Fractal Fract. 2023, 7, 395. [Google Scholar] [CrossRef]
  14. Tatom, F.B. The Relationship between fractional calculus and fractals. Fractals 1995, 3, 217–229. [Google Scholar] [CrossRef]
  15. Mao, X. Stochastic Differential Equations and Applications; Ellis Horwood: Chichester, UK, 1997. [Google Scholar]
  16. Stoyanov, I.M.; Bainov, D.D. The averaging method for a class of stochastic differential equations. Ukr. Math. J. 1974, 26, 186–194. [Google Scholar] [CrossRef]
  17. Lewis, E.R. Delay-line models of polpulation growth. Ecology 1972, 53, 797–807. [Google Scholar] [CrossRef]
Figure 1. Comparison of the exact solution (blue) and the approximate solution (red), with ι = 0.75 (left), ι = 0.95 (right), c 1 = 2 , c 2 = 1 , c 3 = 2 , and α ( ϰ ) = exp ( ϰ 2 / 2 ) , on the interval [ 0.5 , 3 ] .
Figure 1. Comparison of the exact solution (blue) and the approximate solution (red), with ι = 0.75 (left), ι = 0.95 (right), c 1 = 2 , c 2 = 1 , c 3 = 2 , and α ( ϰ ) = exp ( ϰ 2 / 2 ) , on the interval [ 0.5 , 3 ] .
Symmetry 15 01910 g001
Figure 2. Comparison of the exact solution (blue) and the approximate solution (red), with ι = 0.5 (left), ι = 0.65 (right), c 1 = 3 , c 2 = 1 , c 3 = 1 , and α ( ϰ ) = exp ( ϰ 2 / 2 ) , on the interval [ 0.5 , 3 ] .
Figure 2. Comparison of the exact solution (blue) and the approximate solution (red), with ι = 0.5 (left), ι = 0.65 (right), c 1 = 3 , c 2 = 1 , c 3 = 1 , and α ( ϰ ) = exp ( ϰ 2 / 2 ) , on the interval [ 0.5 , 3 ] .
Symmetry 15 01910 g002
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

Rhaima, M.; Mchiri, L.; Ben Makhlouf, A. The Existence and Averaging Principle for a Class of Fractional Hadamard Itô–Doob Stochastic Integral Equations. Symmetry 2023, 15, 1910. https://doi.org/10.3390/sym15101910

AMA Style

Rhaima M, Mchiri L, Ben Makhlouf A. The Existence and Averaging Principle for a Class of Fractional Hadamard Itô–Doob Stochastic Integral Equations. Symmetry. 2023; 15(10):1910. https://doi.org/10.3390/sym15101910

Chicago/Turabian Style

Rhaima, Mohamed, Lassaad Mchiri, and Abdellatif Ben Makhlouf. 2023. "The Existence and Averaging Principle for a Class of Fractional Hadamard Itô–Doob Stochastic Integral Equations" Symmetry 15, no. 10: 1910. https://doi.org/10.3390/sym15101910

APA Style

Rhaima, M., Mchiri, L., & Ben Makhlouf, A. (2023). The Existence and Averaging Principle for a Class of Fractional Hadamard Itô–Doob Stochastic Integral Equations. Symmetry, 15(10), 1910. https://doi.org/10.3390/sym15101910

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop