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
, 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:
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 be a complete probability space and be a standard Brownian motion with respect the space . Denote by the family of functions from to that are right-continuous and have limits on the left. is equipped with the norm and for any . Denote by , the set of all -measurable, -valued random variables satisfying .
Definition 1. The Hadamard fractional integral of order ς for a function h is defined as 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
all be Borel measurable. Consider the following FHIDSIE:
with initial condition
and .
Take the following hypotheses:
: There exists
that satisfies
for all
.
: There exists
that satisfies
for all
.
: Assume that
and there is
such that
for all
.
3. Existence and Uniqueness of Solution
Theorem 1. Under assumptions –, there exists a unique solution of system (1) that satisfies . Before demonstrating Theorem 1, we require the following lemma.
Lemma 1. Assume that assumptions and hold. If is a solution to Equation (1) with initial condition (2), then Moreover, .
Proof. For any
, denote by
a stopping time. It is obvious that
a.s. Let
, for
. Thus, for
,
where
By Lemma 2.3 in [
15] and assumption
, we can derive that
It is not hard to see that
. Then,
Using the Cauchy–Schwartz inequality and assumption
, we obtain
Therefore, applying Theorem
(p. 40 in [
15]) and hypothesis
, one can obtain
Plugging (
12) into (
10), we can obtain
By Gronwall’s inequality, one has
Therefore, we can deduce the desired result when . □
Proof of Theorem 1:
Proof. Uniqueness:
Consider
and
to be two solutions to Equation (
1). Using Lemma 1,
and
belong to
. Thus,
where
By Lemma 2.3 in [
15] and assumption
, we can see that
Using (
15), it yields that
Moreover, by Hölder’s inequality, Theorem
(p. 40 in [
15]), and
, using a similar method as in the proof of Lemma 1, we obtain
Using the Gronwall inequality, we can derive that: . Consequently, for all , , a.s.
Existence: Let
be sufficiently small such that
and write
. We split the proof into two steps:
Step 1: Let
and
for
. For each
, let
, and we define, using the Picard iterations:
It is not hard to see that
. Noting that for
,
Proceeding as in the proof of the uniqueness, we have
where
Noting that for
and
, we have
Proceeding as in the proof of the uniqueness and using (
18), we can obtain
By (
16), we can prove from (
19) that there exists a solution
to Equation (
1) on
in the same way as in the proof of Theorem 3.1 in [
12].
Step 2: For
and
: let
for
and
for
. For each
, let
for
. We define the following Picard iterations:
Proceeding as in step 1, we obtain the existence of the solution on .
Repeating this on
, we can observe that there exists a solution to Equation (
1) on the entire interval
, which complete the proof. □
4. Averaging Principle
Take the following standard form of Equation (
1):
where
, and
satisfy
and
, and
, with
, is a fixed number. Using Theorem 1, Equation (
21) has a unique global solution
, for any
.
Suppose that the following hypothesis holds true:
: Let the measurable functions , and exist, satisfying and . For any , we have
with as positive bounded functions satisfying for .
Our objective is to establish that the solution
can be approximated by the solution
of the next equation
for
.
Now, we will show the first main theorem in this section.
Theorem 2. Assume that – hold. For a given arbitrary small constant , there are two constants and that satisfy , Proof. Using Lemma 2.3 in [
15] and assumption
, one can derive
According to the elementary inequality, one has
where
By the Hölder inequality and
, one has
As in the proof of Theorem 1 in [
16], one can derive
By employing Theorem
(p. 40 in [
15]) and
, one can obtain
Using
, one can derive
By Hölder’s inequality and
, one can derive
Using Hölder’s inequality and
, one has
Plugging (
31)–(
39) into (
25), using Lemma 1 and (
23), one can derive that
with
and
By Gronwall’s inequality, one can obtain
Therefore, given any number
, there are two constants
and
that satisfy, for every
,
□
Next, we show the convergence in probability between and .
Theorem 3. Suppose that the FHIDSIE (21) and (22) both satisfy –. For a given arbitrary small constant , there are two constants and that satisfy , Proof. Given any constant
, using Theorem
and the inequality of Chebyshev, we can obtain
Letting
, (
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]):
with initial condition
where
The assumptions , , and are satisfied for and .
We obtain the corresponding averaged equation
where
Now, we present a numerical simulation of Equations (
44) and (
46), initialized with the condition
, while employing the selected parameters
and
. Subsequently, in
Figure 1 and
Figure 2, we offer a comparative analysis between the precise solution
for Equation (
44) and the averaged solution
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 . 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
- Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
- Podlubny, I. Fractional Differential Equations; Academic Press: Cambridge, MA, USA, 1999. [Google Scholar]
- 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]
- Atanackovic, T.M.; Pilipovic, S.; Stankovic, B.; Zorica, D. Fractional Calculus with Applications in Mechanics; Wiley-ISTE: Indianapolis, IN, USA, 2014. [Google Scholar]
- Baleanu, D.; Machado, J.A.; Luo, A.C. Fractional Dynamics and Control; Springer Science and Business Media: New York, NY, USA, 2011. [Google Scholar]
- Pedjeu, J.C.; Ladde, G.S. Stochastic fractional differential equations: Modeling, method and analysis. Chaos Solitons Fractals 2012, 45, 279–293. [Google Scholar] [CrossRef]
- 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]
- Almeida, R. Caputo-Hadamard fractional derivatives of variable order. Numer. Funct. Anal. Optim. 2017, 38, 1–19. [Google Scholar] [CrossRef]
- 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]
- 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]
- Wang, W.; Guo, Z. Optimal index and averaging principle for Itô–Doob stochastic fractional differential equations. Stoch. Dyn. 2022, 22, 2250018. [Google Scholar] [CrossRef]
- 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]
- Hattaf, K. A New Class of Generalized Fractal and Fractal-Fractional Derivatives with Non-Singular Kernels. Fractal Fract. 2023, 7, 395. [Google Scholar] [CrossRef]
- Tatom, F.B. The Relationship between fractional calculus and fractals. Fractals 1995, 3, 217–229. [Google Scholar] [CrossRef]
- Mao, X. Stochastic Differential Equations and Applications; Ellis Horwood: Chichester, UK, 1997. [Google Scholar]
- 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]
- Lewis, E.R. Delay-line models of polpulation growth. Ecology 1972, 53, 797–807. [Google Scholar] [CrossRef]
| 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. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).