Next Article in Journal
Partial Derivatives Estimation of Multivariate Variance Function in Heteroscedastic Model via Wavelet Method
Previous Article in Journal
Nonlinear Tolerancing: Variation Simulation and Assembly Analysis with Regard to Contact Interaction of Parts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Existence and Averaging Principle for Caputo Fractional Stochastic Delay Differential Systems with Poisson Jumps

1
XJTLU Wisdom Lake Academy of Pharmacy, Xi’an Jiaotong-Liverpool University, Suzhou 215123, China
2
Department of Mathematics, Huaiyin Normal University, Huaian 223300, China
*
Author to whom correspondence should be addressed.
Submission received: 23 November 2023 / Revised: 12 January 2024 / Accepted: 16 January 2024 / Published: 20 January 2024
(This article belongs to the Special Issue Fractional and Stochastic Differential Equations in Mathematics)

Abstract

:
In this paper, we obtain the existence and uniqueness theorem for solutions of Caputo-type fractional stochastic delay differential systems(FSDDSs) with Poisson jumps by utilizing the delayed perturbation of the Mittag–Leffler function. Moreover, by using the Burkholder–Davis–Gundy inequality, Doob’s martingale inequality, and Hölder inequality, we prove that the solution of the averaged FSDDSs converges to that of the standard FSDDSs in the sense of L p . Some known results in the literature are extended.

1. Introduction

Fractional stochastic delay differential systems (FSDDs) are mathematical models that involve fractional derivatives, stochastic noise, and time delays. The fractional derivatives represent the memory effects and long-range dependence in the system, while the stochastic noise and delays account for the random fluctuations and time delays, respectively. FSDDs find applications in many fields, including physics, biology, finance, and engineering. They can be used to model systems with memory and randomness, such as anomalous diffusion processes, fractional-order control systems with stochastic disturbances, and biological systems with fractional-order kinetics and stochastic effects. They provide a powerful framework for understanding and predicting the behavior of complex systems with memory, randomness, and time delays. See, for example, [1,2,3,4,5,6] and the references cited therein.
The averaging principle is a mathematical tool used to simplify the analysis of dynamical systems with fast and slow time scales. It provides an approximate description of the system’s behavior. In 1968, Khasminskii [7] first used the average principle to prove that the solution of the average equation converges to the solution of the corresponding equation. In [8], the authors presented an averaging method for stochastic differential equations with non-Gaussian Lévy noise. Due to the importance of fractional calculus in theory and application, many works have emerged that apply the averaging principle to fractional stochastic differential equations (FSDEs). In [9], Xu et al. presented an averaging principle for Caputo FSDEs driven by Brown motion. In [10], Luo et al. established an averaging principle for the solution of a class of FSDEs with time delays. In the sense of a mean square, Ahmed and Zhu [11] studied the averaging principle for the Hilfer FSDEs with Poisson jumps. In [12], Ahmed investigated the periodic averaging method for impulsive and conformable FSDEs with Poisson jumps. In [13], Wang and Lin consider the following FSDEs
C D 0 α [ x ( t ) h ( t , x ( t ) ] = f ( t , x ( t ) ) + g ( t , x ( t ) ) d B t d t , t J = [ 0 , T ] , x ( 0 ) = x 0 ,
the main results obtained extend some of the works on the average principle of FSDES [9,10,14] from L 2 convergence to L p convergence ( p 2 ). In [15], Yang, et al. studied the averaging principle for a class of ψ -Caputo FSDDEs with Poisson jumps.
Recently, Li and Wang in [16] investigated the existence, uniqueness, and averaging principle for the following Caputo-type FSDDEs:
( C D 0 α Y ) ( t ) = A Y ( t ) + B Y ( t h ) + f ( t , Y ( t ) ) + σ ( t , Y ( t ) ) d W ( t ) d t , t J , Y ( t ) = Φ ( t ) , h t 0 , h > 0 .
Motivated by [11,13,16], we will study the following Caputo FSDDSs with Poisson jumps
( C D 0 α x ) ( t ) = A x ( t ) + B x ( t σ ) + g ( t , x ( t ) , x ( t σ ) ) + κ ( t , x ( t ) , x ( t σ ) ) d W ( t ) d t , + V f ( t , x ( t ) , x ( t σ ) , v ) N ¯ ( d t , d v ) , t J , x ( t ) = ϕ ( t ) , σ t 0 ,
where C D 0 α is the left Caputo fractional derivative with 1 2 < α < 1 , J = [ 0 , T ] , A , B R n × n are two constant matrices, the state vector x R n is a stochastic process, g : J × R n × R n R n , κ : J × R n × R n R n × m and f : J × R n × R n × V R n are measurable continuous functions, W ( t ) is an m-dimensional Brownian motion on the probability space ( Ω , F , P ) . Let ( V , Φ , λ ( d v ) ) be a σ -finite measurable space. Define N ¯ ( t , d v ) : = N ( t , d v ) t λ ( d v ) , where N ( t , d v ) is the counting measure of the stationary Poisson point process p t .
In this paper, we first prove the existence and uniqueness of solutions of Caputo-type FSDDSs (1) using the delayed perturbation of the Mittag–Leffler function and Banach fixed-point theorem; secondly, we prove the averaging principle for Caputo FSDDSs (1) in the sense of L p (pth moment) with inequality techniques. The main contributions and advantages of this paper are as follows:
(1)
The solution of the averaged FSDDSs converges to that of the standard FSDDSs in the sense of L p , which is a generalization of the existing result ( p = 2 ) of the averaging principle for FSDDSs;
(2)
Stochastic inequality, fractional calculus, and Hölder inequality are utilized to establish our results very effectively.
(3)
Our work in this article is innovative. Our result extends the main results of [16].
The remainder of this paper is arranged as follows. In Section 2, we give some definitions and preliminaries. In Section 3, we prove the existence and uniqueness of solutions for Caputo-type FSDDSs (1) with Poisson jumps. In Section 4, we prove that the solution of the FSDDSs (1) converges to that of the standard one in the L p sense. In Section 5, two examples are presented to illustrate our theoretical results. Finally, the paper is concluded in Section 6.

2. Preliminaries

Let Υ = L p ( Ω , F , P ) denote the space of all F ( t ) -measurable, p-square integrable functions x : Ω R n with x ( t ) p s : = i = 1 n E ( | x i ( t ) | p ) 1 / p , and x = i = 1 n | x i | and A = max 1 j n i = 1 n | a i j | be the vector norm and matrix norm, respectively. A process x : [ σ , T ] L p ( Ω , F , P ) is said to be F ( t ) -adapted if x ( t ) Υ .
Definition 1 ([17]).  
Let α > 0 and f be a real function defined on [ a , b ] . The left Riemann–Liouville fractional integral operator of order α is defined by
a I t α f ( t ) = 1 Γ ( α ) a t ( t s ) α 1 f ( s ) d s , t > a .
Definition 2 ([17]).  
Let n 1 < α < n and f C n ( [ a , b ] ) . The left Caputo fractional derivative of order α is defined by
a C D t α f ( t ) = ( a I t n α f ( n ) ) ( t ) = 1 Γ ( n α ) a t ( t s ) n α 1 f ( n ) ( s ) d s , t > a ,
where n = [ α ] + 1 .
Definition 3 ([18]).  
The coefficient matrices Q k ( s ) , k = 0 , 1 , 2 , . . . , satisfy the following multivariate determining matrix equation
Q 0 ( s ) = Q k ( τ ) = Θ , Q 1 ( 0 ) = I , k = 0 , 1 , 2 , , s = 0 , τ , 2 τ , ,
Q k + 1 ( s ) = A Q k ( s ) + B Q k ( s τ ) , k = 0 , 1 , 2 , , s = 0 , τ , 2 τ , ,
where I is an identity matrix and Θ is a zero matrix.
Definition 4 ([18]).  
Delayed perturbation of two parameter Mittag–Leffler-type matrix function X σ , α , β A , B generated by A , a n d   B is defined by
X σ , α , β A , B ( t ) : = Θ , t [ σ , 0 ) , I , t = 0 , i = 0 Q i + 1 ( 0 ) t i α + β 1 Γ ( i α + β ) + i = 1 Q i + 1 ( σ ) ( t σ ) i α + β 1 Γ ( i α + β ) + + i = 0 Q i + 1 ( p σ ) ( t p σ ) i α + β 1 Γ ( i α + β ) , p σ < t ( p + 1 ) σ .
From [18], we can easily obtain the following definition.
Definition 5. 
A R n -value stochastic process { x ( t ) : σ t T } is called a solution of (1) if x ( t ) satisfies the following form:
x ( t ) = X σ , α , 1 A , B ( t + σ ) ϕ ( σ ) + σ 0 X σ , α , α A , B ( t s ) [ C D σ + α ϕ ) ( s ) A ϕ ( s ) ] d s + 0 t X σ , α , α A , B ( t s ) g ( s , x ( s ) , x ( s σ ) ) d s + 0 t X σ , α , α A , B ( t s ) κ ( s , x ( s ) , x ( s σ ) ) d W ( s ) + 0 t X σ , α , α A , B ( t s ) V f ( s , x ( s ) , x ( s σ ) , v ) N ¯ ( d s , d v ) , t J , ϕ ( t ) , t [ σ , 0 ] ,
where x ( t ) is F ( t ) -adapted and E ( σ T x ( t ) p d t ) < .
Lemma 1 ([19]).  
For each t 0 , 0 < α < 1 , 0 < β 1 and α + β 1 , one has
X σ , α , β A , B ( t ) t β 1 E α , β ( ( A + B ) t α ) ,
where E α , β ( z ) = k = 0 z k Γ ( k α + β ) , z R is the Mittag–Leffler function.
Lemma 2. 
For any p 2 , α 1 1 p , 1 and μ > 0 , we have
0 t ( t s ) p α p E p α p + 1 , 1 ( μ s p α p + 1 ) d s Γ ( p α p + 1 ) μ E p α p + 1 , 1 ( μ t p α p + 1 ) ,
where Γ ( α ) : = 0 + s α 1 e s d s is the Gamma function.
Proof. 
Let μ > 0 be arbitrary. Consider the corresponding linear Caputo fractional differential equation of the following form
C D 0 + p α p + 1 x ( t ) = μ x ( t ) .
From [20], it is easy to know that the Mittag–Leffler function E p α p + 1 , 1 ( μ t p α p + 1 ) is a solution of (8). So, the following equality holds:
E p α p + 1 , 1 ( μ t p α p + 1 ) = 1 + μ Γ ( p α p + 1 ) 0 t ( t s ) p α p E p α p + 1 , 1 ( μ s p α p + 1 ) d s ,
which completes the proof.  □
Lemma 3 ([21,22]).  
Let ϕ : R + × V R n and assume that
0 t V | ϕ ( s , v ) | p λ ( d v ) d s < , p 2 .
Then, there exists D p > 0 such that
E sup 0 t u 0 t V ϕ ( s , v ) N ¯ ( d s , d v ) p D p E 0 u V | ϕ ( s , v ) | 2 λ ( d v ) d s p 2 + E 0 u V | ϕ ( s , v ) | p λ ( d v ) d s .
Lemma 4 ([23]).  
Let u , v be two integrable functions and g be continuously defined on the domain [ a , b ] . Suppose that
(1)
u and v are non-negative, and v is non-decreasing;
(2)
g is non-negative and non-decreasing.
If
u ( t ) v ( t ) + g ( t ) a t ( t τ ) α 1 u ( τ ) d τ ,
then
u ( t ) v ( t ) E α ( g ( t ) Γ ( α ) ( t a ) α ) , t [ a , b ] ,
where E α ( · ) is the Mittag–Leffler function.
To study the problem (1), we impose the following conditions:
(H1) For each x 1 , x 2 , y 1 , y 2 R n and t J , there exist two constants C 1 , C 2 > 0 such that g ( t , x 1 , y 1 ) g ( t , x 2 , y 2 ) p   κ ( t , x 1 , y 1 ) κ ( t , x 2 , y 2 ) p
V f ( t , x 1 , y 1 , v ) f ( t , x 2 , y 2 , v ) p λ ( d v ) C 1 p ( x 1 x 2 p +   y 1 y 2 p ) , where · is the norm of R n , x y = m a x { x , y } .
(H2) Let κ ( · , 0 , 0 ) and f ( · , 0 , 0 , 0 ) be essentially bounded, i.e.,
κ ( · , 0 , 0 ) : = e s s sup t [ 0 , ) κ ( t , 0 , 0 )   < + , f ( · , 0 , 0 , 0 ) : = e s s sup t [ 0 , ) f ( t , 0 , 0 , 0 ) < + ,
and g ( · , 0 , 0 ) is L p integrable, i.e.,
g L p = 0 T g ( t , 0 , 0 ) p d t < + .

3. Existence and Uniqueness Result

Let H p ( [ 0 , T ] ) be the space of all the processes x which are measurable, F ( t ) -adapted, and satisfied that x H p : = sup 0 t T x ( t ) p s < . Obviously, ( H p ( [ 0 , T ] ) , · H p ) is a Banach space. Set γ =   A + B . For each t [ σ , T ] and ϕ C ( [ σ , 0 ] , R n ) , we define an operator T : H p ( [ 0 , T ] ) H p ( [ 0 , T ] ) as follows:
( T x ) ( t ) = X σ , α , 1 A , B ( t + σ ) ϕ ( σ ) + σ 0 X σ , α , α A , B ( t s ) [ C D σ + α ϕ ) ( s ) A ϕ ( s ) ] d s + 0 t X σ , α , α A , B ( t s ) g ( s , x ( s ) , x ( s σ ) ) d s + 0 t X σ , α , α A , B ( t s ) κ ( s , x ( s ) , x ( s σ ) ) d W ( s ) + 0 t X σ , α , α A , B ( t s ) V f ( s , x ( s ) , x ( s σ ) , v ) N ¯ ( d s , d v ) .
Lemma 5. 
Let 1 1 p < α < 1 . Assume that (H1) and (H2) hold. Then, the operator T is well defined.
Proof. 
For any x H p ( [ 0 , T ] ) , by (10) and the following elementary inequality,
i = 1 m a i p m p 1 i = 1 m a i p , a i R n , i = 1 , 2 , . . . , m .
we have
( T x ) ( t ) p s p 5 p 1 E ( X σ , α , 1 A , B ( t + σ ) ϕ ( σ ) p ) + 5 p 1 E σ 0 X σ , α , α A , B ( t s ) [ ( C D σ + α ϕ ) ( s ) A ϕ ( s ) ] d s p + 5 p 1 E 0 t X σ , α , α A , B ( t s ) g ( s , x ( s ) , x ( s σ ) ) d s p + 5 p 1 E 0 t X σ , α , α A , B ( t s ) κ ( s , x ( s ) , x ( s σ ) ) d W ( s ) p + 5 p 1 E 0 t X σ , α , α A , B ( t s ) V f ( s , x ( s ) , x ( s σ ) , v ) N ¯ ( d s , d v ) p : = I 1 + I 2 + I 3 + I 4 + I 5 .
For I 1 , from Lemma 1, one has
I 1 = 5 p 1 E ( X σ , α , 1 A , B ( t + σ ) ϕ ( σ ) p ) 5 p 1 E ( X σ , α , 1 A , B ( t + σ ) p ϕ ( σ ) p ) 5 p 1 ϕ ( σ ) p ( E α , 1 ( γ ( T + σ ) α ) ) p .
For I 2 , by Lemma 1, Hölder inequality, and α > 1 1 p , we obtain
I 2 = 5 p 1 E σ 0 X σ , α , α A , B ( t s ) [ ( C D σ + α ϕ ) ( s ) A ϕ ( s ) ] d s p 5 p 1 σ 0 X σ , α , α A , B ( t s ) p d s · E σ 0 ( C D σ + α ϕ ) ( s ) A ϕ ( s ) q d s p 1 5 p 1 Ξ ( T + σ ) p α p + 1 p α p + 1 ( E α , α ( γ ( T + σ ) α ) ) p ,
where 1 p + 1 q = 1 and Ξ = σ 0 ( C D σ + α ϕ ) ( s ) A ϕ ( s ) q d s p 1 < .
For I 3 , applying (H1), (H2), Hölder inequality, Lemma 1 and Jensen inequality, one has
I 3 = 5 p 1 E 0 t X σ , α , α A , B ( t s ) g ( s , x ( s ) , x ( s σ ) ) d s p 5 p 1 0 t X σ , α , α A , B ( t s ) q d s p q · E 0 t g ( s , x ( s ) , x ( s σ ) ) g ( s , 0 , 0 ) + g ( s , 0 , 0 ) p d s 5 p 1 0 t t q ( α 1 ) E α , α ( μ ( t s ) α ) q d s p q · 2 p 1 E 0 t g ( s , x ( s ) , x ( s σ ) ) g ( s , 0 , 0 ) p d s + 0 t g ( s , 0 , 0 ) p d s 10 p 1 E α , α ( γ T α ) p T q α q + 1 q α q + 1 p q E 0 t C 1 p ( x ( s ) p + x ( s σ ) p ) d s + 0 t g ( s , 0 , 0 ) p d s 10 p 1 E α , α ( γ T α ) p T q α q + 1 q α q + 1 p q T C 1 p ( 2 x H p p + ϕ q ) + g L p ,
since
sup 0 t T E x ( t τ ) q max sup τ t 0 E ϕ ( t ) q , sup 0 t T E x ( t ) q = max ϕ q , x H q q ϕ q + x H q q .
For I 4 , by using (H1), (H2), Cauchy–Schwarz inequality, Ito’s isometry, Lemma 1, and Jensen inequality, we have
I 4 = 5 p 1 E 0 t X σ , α , α A , B ( t s ) κ ( s , x ( s ) , x ( s σ ) ) d W ( s ) 2 p 2 5 p 1 E 0 t X σ , α , α A , B ( t s ) 2 κ ( s , x ( s ) , x ( s σ ) ) 2 d s p 2 5 p 1 E 0 t 1 p p 2 d s p 2 p · 0 t X σ , α , α A , B ( t s ) p κ ( s , x ( s ) , x ( s σ ) ) p d s 2 p p 2 5 p 1 T p 2 1 E 0 t X σ , α , α A , B ( t s ) p κ ( s , x ( s ) , x ( s σ ) ) p d s 5 p 1 T p 2 1 2 p 1 E α , α ( γ T α ) p E 0 t ( t s ) p α p [ C 1 p ( x ( s ) p + x ( s σ ) p ) + κ ( s , 0 , 0 ) p ] d s 10 p 1 T p α p 2 p α p + 1 E α , α ( γ T α ) p ( C 1 p ( 2 x H p p + ϕ p ) + κ ( · , 0 , 0 ) p ) .
For I 5 , using (H1), (H2), Lemmas 1, 3 and Jensen inequality, we obtain
I 5 = 5 p 1 E 0 t V X σ , α , α A , B ( t s ) f ( s , x ( s ) , x ( s σ ) , v ) N ¯ ( d s , d v ) p 5 p 1 D p E 0 t V | X σ , α , α A , B ( t s ) | 2 f 2 ( s , x ( s ) , x ( s σ ) , v ) λ ( d v ) d s p 2 + 5 p 1 D p E 0 t V | X σ , α , α A , B ( t s ) | p f p ( s , x ( s ) , x ( s σ ) , v ) λ ( d v ) d s 5 p 1 D p ( T p 2 1 + 1 ) E 0 t | X σ , α , α A , B ( t s ) | p V f p ( s , x ( s ) , x ( s σ ) , v ) λ ( d v ) d s 5 p 1 D p ( T p 2 1 + 1 ) 2 p 1 E α , α ( γ T α ) p · E 0 t ( t s ) p α p [ C 1 p ( x ( s ) p + x ( s σ ) p ) + f ( s , 0 , 0 , 0 ) p ] d s 10 p 1 D p ( T p 2 1 + 1 ) T p α p + 1 p α p + 1 E α , α ( γ T α ) p ( C 1 p ( 2 x H p p + ϕ p ) + f ( · , 0 , 0 , 0 ) p ) .
Submitting (13)–(17) into (12) implies that T x H p < . Thus, the operator T is well-defined.  □
Theorem 1. 
Let 1 1 p < α < 1 . Assume that ( H 1 ) and ( H 2 ) hold, then (1) has a unique solution x H p ( [ 0 , T ] ) .
Proof. 
For T > 0 , we choose and fix a constant μ > 0 such that
μ > 2 · 3 p 1 C 1 p E α , α ( γ T α ) p ( T p q + ( D p + 1 ) T p 2 1 + 1 ) Γ ( p α p + 1 ) .
On the space H p ( [ 0 , T ] ) , we define a weighted norm · μ as below
x μ : = sup t [ 0 , T ] E ( x ( t ) p ) E p α p + 1 , 1 ( μ t p α p + 1 ) 1 p , x H p ( [ 0 , T ] ) .
Similarly to Theorem 1 in [18], it is easy to know that the norms · H p and · μ are equivalent. Hence, ( H p ( [ 0 , T ] ) , · μ ) is a Banach space. We can easily prove that T : H p ( [ 0 , T ] ) H p ( [ 0 , T ] ) defined in (10) is uniformly bounded operator by Lemma 5. Next, we only check that T is a contraction operator.
Firstly, by using Hölder inequality (H1) and Lemma 1, we obtain
0 t X σ , α , α A , B ( t s ) ( g ( s , x ( s ) , x ( s σ ) ) g ( s , y ( s ) , y ( s σ ) ) ) d s p 0 t 1 q d s p q · 0 t X σ , α , α A , B ( t s ) p g ( s , x ( s ) , x ( s σ ) ) g ( s , y ( s ) , y ( s σ ) ) p d s t p q 0 t ( t s ) p ( α 1 ) E α , α ( γ ( t s ) α ) p g ( s , x ( s ) , x ( s σ ) ) g ( s , y ( s ) , y ( s σ ) ) p d s T p q E α , α ( γ T α ) p C 1 p 0 t ( t s ) p ( α 1 ) ( x ( s ) y ( s ) p + x ( s σ ) y ( s σ ) p ) d s .
Secondly, similarly to the Proof of (16), one has
0 t X σ , α , α A , B ( t s ) ( κ ( s , x ( s ) , x ( s σ ) ) κ ( s , y ( s ) , y ( s σ ) ) ) d W ( s ) p T p 2 1 0 t X σ , α , α A , B ( t s ) p κ ( s , x ( s ) , x ( s σ ) ) κ ( s , y ( s ) , y ( s σ ) ) p d s T p 2 1 E α , α ( γ T α ) p 0 t ( t s ) p α p [ C 1 p ( x ( s ) y ( s ) p + x ( s σ ) y ( s σ ) p ) d s T p 2 1 C 1 p E α , α ( γ T α ) p 0 t ( t s ) p α p ( x ( s ) y ( s ) p + x ( s σ ) y ( s σ ) p ) d s .
Thirdly, similarly to the Proof of (17), we obtain
0 t X σ , α , α A , B ( t s ) V ( f ( s , x ( s ) , x ( s σ ) , v ) f ( s , y ( s ) , y ( s σ ) , v ) ) N ¯ ( d s , d v ) p D p 0 t X σ , α , α A , B ( t s ) 2 V f ( s , x ( s ) , x ( s σ ) , v ) f ( s , y ( s ) , y ( s σ ) , v ) 2 λ ( d v ) d s p 2 + D p 0 t | X σ , α , α A , B ( t s ) | p V | f ( s , x ( s ) , x ( s σ ) , v ) f ( s , y ( s ) , y ( s σ ) , v ) | p λ ( d v ) d s D p ( T p 2 1 + 1 ) 0 t X σ , α , α A , B ( T s ) p V f ( s , x ( s ) , x ( s σ ) , v ) f ( s , y ( s ) , y ( s σ ) , v ) p λ ( d v ) d s D p C 1 p ( T p 2 1 + 1 ) E α , α ( γ T α ) p 0 t ( t s ) p α p ( x ( s ) y ( s ) p + x ( s σ ) y ( s σ p ) d s .
For each x , y H p ( [ 0 , T ] ) , from (10), (11), and (19)–(21), we have
E ( T x ( t ) T y ( t ) p ) 3 p 1 E 0 t X σ , α , α A , B ( t s ) ( g ( s , x ( s ) , x ( s σ ) ) g ( s , y ( s ) , y ( s σ ) ) ) d s p + 3 p 1 E 0 t X σ , α , α A , B ( t s ) ( κ ( s , x ( s ) , x ( s σ ) ) κ ( s , y ( s ) , y ( s σ ) ) ) d W ( s ) p + 3 p 1 E 0 t X σ , α , α A , B ( t s ) V ( f ( s , x ( s ) , x ( s σ ) , v ) f ( s , y ( s ) , y ( s σ ) , v ) ) N ¯ ( d s , d v ) p 3 p 1 T p q C 1 p E α , α ( γ T α ) p 0 t ( t s ) p ( α 1 ) E ( x ( s ) y ( s ) p + x ( s σ ) y ( s σ p ) d s + 3 p 1 T p 2 1 C 1 p E α , α ( γ T α ) p 0 t ( t s ) p α p E ( x ( s ) y ( s ) p + x ( s σ ) y ( s σ p ) d s + 3 p 1 C 1 p D p ( T p 2 1 + 1 ) E α , α ( γ T α ) p 0 t ( t s ) p α p E ( x ( s ) y ( s ) p + x ( s σ ) y ( s σ p ) d s = ω 0 t ( t s ) p α p E ( x ( s ) y ( s ) p + x ( s σ ) y ( s σ p ) d s ,
where
ω : = 3 p 1 C 1 p E α , α ( γ T α ) p ( T p q + ( D p + 1 ) T p 2 1 + D p ) .
For t > σ , one has
0 t ( t s ) p ( α 1 ) x ( s σ ) y ( s σ ) p d s = 0 σ + σ t ( t s ) p ( α 1 ) x ( s σ ) y ( s σ ) p d s = σ t ( t s ) p ( α 1 ) x ( s σ ) y ( s σ ) p d s = 0 t σ ( t σ u ) p ( α 1 ) x ( u ) y ( u ) p d u .
From Lemma 2, combining (22) and (23) for each t [ 0 , T ] , we obtain
E ( T x ( t ) T y ( t ) p ) E p α p + 1 , 1 ( μ t p α p + 1 ) ω E p α p + 1 , 1 ( μ t p α p + 1 ) 0 t ( t s ) p α p E p α p + 1 , 1 ( μ s p α p + 1 ) d s x y μ p + ω E p α p + 1 , 1 ( μ ( t σ ) p α p + 1 ) 0 t σ ( t σ u ) p α p E p α p + 1 , 1 ( μ u p α p + 1 ) d u x y μ p 2 ω Γ ( p α p + 1 ) μ x y μ p ,
which implies that
T x T y μ ρ x y μ ,
where ρ = 2 ω Γ ( p α p + 1 ) μ 1 p .
Based on (18), one can obtain ρ < 1 and the operator T is a contractive. Thus, (1) has a unique solution using the Banach fixed-point theorem. This completes the proof of Theorem 1.  □

4. An Averaging Principle

To show the averaging principle for FSDDEs (1), let us consider the following standard form of (1)
x ϵ ( t ) = X σ , α , 1 A , B ( t + σ ) ϕ ( σ ) + σ 0 X σ , α , α A , B ( t s ) [ C D σ + α ϕ ) ( s ) A ϕ ( s ) ] d s + ϵ 0 t X σ , α , α A , B ( t s ) g ( s , x ϵ ( s ) , x ϵ ( s σ ) ) d s + ϵ 0 t X σ , α , α A , B ( t s ) κ ( s , x ϵ ( s ) , x ϵ ( s σ ) ) d W ( s ) + ϵ 0 t X σ , α , α A , B ( t s ) V f ( s , x ϵ ( s ) , x ϵ ( s σ ) , v ) N ¯ ( d s , d v ) ,
where ϵ ( 0 , ϵ 0 ] is a positive small parameter with ϵ 0 being a fixed number.
Consider the averaged form which corresponds to the standard form (24) as follows:
y ϵ ( t ) = X σ , α , 1 A , B ( t + σ ) ϕ ( σ ) + σ 0 X σ , α , α A , B ( t s ) [ C D σ + α ϕ ) ( s ) A ϕ ( s ) ] d s + ϵ 0 t X σ , α , α A , B ( t s ) g ^ ( s , y ϵ ( s ) , y ϵ ( s σ ) ) d s + ϵ 0 t X σ , α , α A , B ( t s ) κ ^ ( s , y ϵ ( s ) , y ϵ ( s σ ) ) d W ( s ) + ϵ 0 t X σ , α , α A , B ( t s ) V f ^ ( s , y ϵ ( s ) , y ϵ ( s σ ) , v ) N ¯ ( d s , d v ) ,
where g ^ : R n × R n R n , κ ^ : R n × R n R n × m , and f ^ : R n × R n × V R n satisfying the following averaging condition:
(H3) For each t J , x , y R n , and p 2 , there exists a positive bounded function φ i ( · ) , i = 1 , 2 , 3 such that
1 t 0 t g ( s , x , y ) g ^ ( x , y ) p d s φ 1 ( t ) ( 1 + x p + y p ) ,
1 t 0 t ( t s ) α 1 ( κ ( s , x , y ) κ ^ ( x , y ) ) p d s φ 2 ( t ) ( 1 + x p + y p ) ,
1 t 0 t V ( t s ) α 1 ( f ( s , x , y , v ) f ^ ( x , y , v ) ) p λ ( d v ) d s φ 3 ( t ) ( 1 + x p + y p ) , where lim t φ i ( t ) = 0 , i = 1 , 2 , 3 .
Theorem 2. 
Assume that (H1)–(H3) are satisfied. Then, for a given arbitrary small number δ > 0 ,   p 2 with 1 1 p < α < 1 , there exist L > 0 , ϵ 1 ( 0 , ϵ 0 ] and β ( 0 , 1 ) such that
E sup t [ τ , L ϵ β ] | x ϵ ( t ) y ϵ ( t ) | p δ ,
for all ϵ ( 0 , ϵ 1 ] .
Proof. 
If p = 2 , it is easy to prove that (26) holds using method similarly to that in [20]. In the following, we will only consider the case of p > 2 . From Equations (25), (26), and inequality (11), we obtain
x ϵ ( t ) y ϵ ( t ) p 3 p 1 ϵ p 0 t X σ , α , α A , B ( t s ) [ g ( s , x ϵ ( s ) , x ϵ ( s σ ) ) g ^ ( y ϵ ( s ) , y ϵ ( s σ ) ) ] d s p + 3 p 1 ϵ p 2 0 t X σ , α , α A , B ( t s ) [ κ ( s , x ϵ ( s ) , x ϵ ( s σ ) ) κ ^ ( y ϵ ( s ) , y ϵ ( s σ ) ) ] d W ( s ) p + 3 p 1 ϵ p 2 0 t X σ , α , α A , B ( t s ) V [ f ( s , x ϵ ( s ) , x ϵ ( s σ ) , v ) ) f ^ ( x ϵ ( s ) , x ϵ ( s σ ) , v ) ) ] N ¯ ( d s , d v ) p .
For any t [ 0 , u ] [ 0 , T ] , taking the expectation on both sides of Equation (27), we have
E sup 0 t u x ϵ ( t ) y ϵ ( t ) p 3 p 1 ϵ p E sup 0 t u 0 t X σ , α , α A , B ( t s ) [ g ( s , x ϵ ( s ) , x ϵ ( s σ ) ) g ^ ( y ϵ ( s ) , y ϵ ( s σ ) ) ] d s p + 3 p 1 ϵ p 2 E sup 0 t u 0 t X σ , α , α A , B ( t s ) [ κ ( s , x ϵ ( s ) , x ϵ ( s σ ) ) κ ^ ( y ϵ ( s ) , y ϵ ( s σ ) ) ] d W ( s ) p + 3 p 1 ϵ p 2 E sup 0 t u 0 t X σ , α , α A , B ( t s ) V [ f ( s , x ϵ ( s ) , x ϵ ( s σ ) , v ) ) f ^ ( x ϵ ( s ) , x ϵ ( s σ ) , v ) ) ] N ¯ ( d s , d v ) p = I 1 + I 2 + I 3 .
Applying Jensen’s inequality, we obtain
I 1 6 p 1 ϵ p · E sup 0 t u 0 t X σ , α , α A , B ( t s ) [ g ( s , x ϵ ( s ) , x ϵ ( s σ ) ) g ( s , y ϵ ( s ) , y ϵ ( s σ ) ) ] d s p + 6 p 1 ϵ p · E sup 0 t u 0 t X σ , α , α A , B ( t s ) [ g ( s , y ϵ ( s ) , y ϵ ( s σ ) ) g ^ ( y ϵ ( s ) , y ϵ ( s σ ) ) ] d s p = I 11 + I 12 .
Thanks to Hölder inequality and (H2), we obtain
I 11 6 p 1 ϵ p 0 u 1 q d s p q · E sup 0 t u 0 t X σ , α , α A , B ( t s ) p g ( s , x ϵ ( s ) , x ϵ ( s σ ) ) g ( s , y ϵ ( s ) , y ϵ ( s σ ) ) p d s 6 p 1 ϵ p u p 1 C 1 p E α , α ( γ u α ) p · E sup 0 t u 0 t ( t s ) p ( α 1 ) [ x ϵ ( s ) y ϵ ( s ) p + x ϵ ( s σ ) ) y ϵ ( s σ ) ) p ] d s 6 p 1 ϵ p u p 1 C 1 p E α , α ( γ u α ) p · 0 u ( u s ) p ( α 1 ) E sup 0 θ s x ϵ ( θ ) y ϵ ( θ ) p + E sup 0 θ s x ϵ ( θ σ ) y ϵ ( θ σ ) p d s 2 · 6 p 1 ϵ p u p 1 C 1 p E α , α ( γ u α ) p 0 u ( t s ) p ( α 1 ) E sup 0 θ s x ϵ ( θ ) y ϵ ( θ ) p d s ,
since
sup 0 θ s x ϵ ( θ τ ) y ϵ ( θ τ ) p sup 0 θ s x ϵ ( θ ) y ϵ ( θ ) p .
Applying Hölder inequality, we obtain
I 12 6 p 1 ϵ p · E sup 0 t u 0 t X σ , α , α A , B ( t s ) q d s p q · 0 t g ( s , y ϵ ( s ) , y ϵ ( s σ ) ) g ^ ( y ϵ ( s ) , y ϵ ( s σ ) ) p d s 6 p 1 ϵ p E α , α ( γ u α ) p u q α q + 1 q α q + 1 p q · u φ 1 1 + E sup 0 t u y ϵ ( t ) p + E sup 0 t u y ϵ ( t σ ) p = 6 p 1 φ 1 M 1 ( q α q + 1 ) ( p 1 ) ϵ p E α , α ( γ u α ) p u p α ,
here φ 1 = sup t [ 0 , u ] | φ 1 ( t ) | , M 1 = 1 + E sup 0 t u y ϵ ( t ) p + E sup 0 t u y ϵ ( t τ ) p .
For the second term I 2 , we have
I 2 6 p 1 ϵ p 2 E sup 0 t u 0 t X σ , α , α A , B ( t s ) [ κ ( s , x ϵ ( s ) , x ϵ ( s σ ) ) κ ( s , y ϵ ( s ) , y ϵ ( s σ ) ) ] d W ( s ) p + 6 p 1 ϵ p 2 E sup 0 t u 0 t X σ , α , α A , B ( t s ) [ κ ( s , y ϵ ( s ) , y ϵ ( s σ ) ) κ ^ ( y ϵ ( s ) , y ϵ ( s σ ) ) ] d W ( s ) p = I 21 + I 22 .
In view of the Burkholder–Davis–Gundy’s inequality, Hölder’s inequality and Doob’s martingale inequality, and (H1), one has
I 21 6 p 1 ϵ p 2 E sup 0 t u 0 t X σ , α , α A , B ( t s ) 2 κ ( s , x ϵ ( s ) , x ϵ ( s σ ) ) κ ( s , y ϵ ( s ) , y ϵ ( s σ ) ) 2 d s p 2 6 p 1 ϵ p 2 E sup 0 t u 0 t 1 p p 2 d s p 2 p · p 2 · 0 t X σ , α , α A , B ( t s ) p κ ( s , x ϵ ( s ) , x ϵ ( s σ ) ) κ ( s , y ϵ ( s ) , y ϵ ( s σ ) ) p d s 6 p 1 ϵ p 2 C 1 p u p 2 1 E α , α ( γ u α ) p · 0 u ( u s ) p α p E sup 0 θ s x ϵ ( θ ) y ϵ ( θ ) p + E sup 0 θ s x ϵ ( θ σ ) y ϵ ( θ σ ) p d s 2 · 6 p 1 ϵ p 2 C 1 p u p 2 1 E α , α ( γ u α ) p 0 u ( u s ) p α p E sup 0 θ s x ϵ ( θ ) y ϵ ( θ ) p d s .
Applying (H3) and an estimation method similar to Equation (33), we obtain
I 22 6 p 1 ϵ p 2 u p 2 1 E sup 0 θ u 0 t X σ , α , α A , B ( t s ) p κ ( s , y ϵ ( s ) , y ϵ ( s σ ) ) κ ^ ( y ϵ ( s ) , y ϵ ( s σ ) ) p d s 6 p 1 ϵ p 2 u p 2 1 E α , α ( γ u α ) p u φ 2 1 + E sup 0 t u | y ϵ ( t ) | p + E sup 0 t u | y ϵ ( t σ ) | p = 6 p 1 M 1 φ 2 ϵ p 2 E α , α ( γ u α ) p u p 2 .
For the third term I 3 , we have
I 3 3 p 1 ϵ p 2 E sup 0 t u 0 t V X σ , α , α A , B ( t s ) [ f ( s , x ϵ ( s ) , x ϵ ( s σ ) , v ) f ( s , y ϵ ( s ) , y ϵ ( s σ ) , v ) ] N ¯ ( d s , d v ) p + 3 p 1 ϵ p 2 E sup 0 t u 0 t V X σ , α , α A , B ( t s ) [ f ( s , y ϵ ( s ) , y ϵ ( s σ ) , v ) f ^ ( y ϵ ( s ) , y ϵ ( s σ ) , v ) ] N ¯ ( d s , d v ) p = I 31 + I 32 .
From Lemma 3, similarly to the Proof of (17), one has
I 31 3 p 1 ϵ p 2 D p E sup 0 t u 0 t V X σ , α , α A , B ( t s ) 2 f ( s , x ϵ ( s ) , x ϵ ( s σ ) , v ) f ( s , y ϵ ( s ) , y ϵ ( s σ ) , v ) 2 λ ( d v ) d s p 2 + 3 p 1 ϵ p 2 E sup 0 t u 0 t V X σ , α , α A , B ( t s ) p f ( s , x ϵ ( s ) , x ϵ ( s σ ) , v ) f ( s , y ϵ ( s ) , y ϵ ( s σ ) , v ) p λ ( d v ) d s 3 p 1 ϵ p 2 ( D p u p 2 1 + 1 ) · E sup 0 t u 0 t V X σ , α , α A , B ( t s ) p f ( s , x ϵ ( s ) , x ϵ ( s σ ) , v ) f ( s , y ϵ ( s ) , y ϵ ( s σ ) , v ) p λ ( d v ) d s 3 p 1 ϵ p 2 ( D p u p 2 1 + 1 ) E α , α ( γ u α ) p C 1 p · 0 u ( u s ) p α p E sup 0 θ s x ϵ ( θ ) y ϵ ( θ ) p + E sup 0 θ s x ϵ ( θ σ ) y ϵ ( θ σ ) p d s 2 · 3 p 1 C 1 p ϵ p 2 ( D p u p 2 1 + 1 ) E α , α ( γ u α ) p 0 u ( t s ) p α p E sup 0 θ s x ϵ ( θ ) y ϵ ( θ ) p d s .
Moreover, by (H3), we also have
I 32 3 p 1 ϵ p 2 D p E sup 0 t u 0 t V X σ , α , α A , B ( t s ) 2 f ( s , y ϵ ( s ) , y ϵ ( s σ ) , v ) f ^ ( y ϵ ( s ) , y ϵ ( s σ ) , v ) 2 λ ( d v ) d s p 2 + 3 p 1 ϵ p 2 E sup 0 t u 0 t V X σ , α , α A , B ( t s ) p f ( s , y ϵ ( s ) , y ϵ ( s σ ) , v ) f ^ ( y ϵ ( s ) , y ϵ ( s σ ) , v ) p λ ( d v ) d s 3 p 1 ϵ p 2 ( D p u p 2 1 + 1 ) · E sup 0 t u 0 t V X σ , α , α A , B ( t s ) p f ( s , y ϵ ( s ) , y ϵ ( s σ ) , v ) f ^ ( y ϵ ( s ) , y ϵ ( s σ ) , v ) p λ ( d v ) d s 3 p 1 ϵ p 2 ( D p u p 2 1 + 1 ) E α , α ( γ u α ) p u φ 3 1 + E sup 0 t u y ϵ ( t ) p + E sup 0 t u y ϵ ( t σ ) p 3 p 1 M 1 φ 3 ϵ p 2 E α , α ( γ u α ) p ( D p u p 2 + u ) .
From (28)–(37), for u ( 0 , T ] , we obtain
E sup 0 t u x ϵ ( t ) y ϵ ( t ) p A ( u ) + B ( u ) 0 u ( u s ) p α p E sup 0 θ s x ϵ ( θ ) y ϵ ( θ ) p d s ,
where
A ( u ) = 6 p 1 φ 1 M 1 ( q α q + 1 ) ( p 1 ) ϵ p E α , α ( γ u α ) p u p α + 6 p 1 M 1 φ 2 ϵ p 2 E α , α ( γ u α ) p u p 2 + 3 p 1 M 1 φ 3 ϵ p 2 E α , α ( γ u α ) p ( D p u p 2 + u ) ,
and
B ( u ) = 2 · 6 p 1 C 1 p ϵ p E α , α ( γ u α ) p u p 1 + 2 · 6 p 1 C 1 p ϵ p 2 E α , α ( γ u α ) p u p 2 1 + 2 · 3 p 1 C 1 p ϵ p 2 E α , α ( γ u α ) p ( D p u p 2 1 + 1 ) .
By using of Lemma 4, we obtain
E sup 0 t u x ϵ ( t ) y ϵ ( t ) p A ( u ) E p ( α 1 ) + 1 B ( u ) Γ ( p ( α 1 ) + 1 ) u p ( α 1 ) + 1 .
Choose L > 0 and β ( 0 , 1 ) such that, for all t ( 0 , L ϵ β ] ( 0 , T ] satisfies the following
E sup 0 < t L ϵ β x ϵ ( t ) y ϵ ( t ) p A ¯ ( ϵ ) E p ( α 1 ) + 1 ( B ¯ ( ϵ ) Γ ( p ( α 1 ) + 1 ) ) ϵ 1 β ,
where
A ¯ ( ϵ ) = 6 p 1 E α , α ( γ T α ) p φ 1 M 1 ( q α q + 1 ) ( p 1 ) L p α ϵ p ( 1 α β ) + 6 p 1 M 1 φ 2 E α , α ( γ T α ) p L p 2 ϵ p 2 ( 1 β ) + 3 p 1 M 1 φ 3 E α , α ( γ T α ) p ( D p L p 2 ϵ p 2 ( 1 β ) + L ϵ p 2 β ) ,
and
B ¯ ( ϵ ) = 2 · 6 p 1 C 1 p E α , α ( γ T α ) p L p 1 ϵ p ( p 1 ) β + 2 · 6 p 1 C 1 p E α , α ( γ T α ) p L p 2 1 ϵ p 2 ( p 2 1 ) β + 2 · 3 p 1 C 1 p E α , α ( γ T α ) p ( D p L p 2 1 ϵ p 2 ( p 2 1 ) β + ϵ p 2 ) ,
are two constants. Thus, for any given number δ > 0 , there exists ϵ 1 ( 0 , ϵ 0 ] such that, for each ϵ ( 0 , ϵ 1 ] and t [ 0 , L ϵ β ] J ,
E sup t [ 0 , L ϵ β ] x ϵ ( t ) y ϵ ( t ) p δ .
Remark 1. 
If p = 2 and f 0 , then FSDDEs (1) reduces to FSDDSs (1.1) in [16]. Therefore, Theorems 1 and 2 generalize the main results of [16].
By using Theorem 2 and Chebyshev–Markov inequality, we can obtain the following corollary.
Corollary 1. 
Assume that (H1)–(H3) are satisfied. Then, for a given arbitrary small number δ > 0 , p 2 with 1 1 p < α < 1 , then for arbitrarily number δ ¯ > 0 such that for L > 0 , ϵ 1 ( 0 , ϵ 0 ] and β ( 0 , 1 ) satisfying for all ϵ ( 0 , ϵ 1 ]
lim ϵ 0 P sup t [ 0 , L ϵ β ] x ϵ ( t ) y ϵ p > δ ¯ = 0 .

5. Applications

In this section, we will provide two examples to illustrate the application of our main results.
Example 1. 
Consider the following Caputo-type FSDDSs with Poisson jumps:
( C D 0 0.9 x ) ( t ) = A x ( t ) + B x ( t 0.4 ) + g ( t , x ( t ) , x ( t 0.4 ) ) + κ ( t , x ( t ) , x ( t 0.4 ) ) d W ( t ) d t + V f ( t , x ( t ) , x ( t 0.4 ) , v ) N ¯ ( d t , d v ) , t J , x ( t ) = ϕ ( t ) , 0.4 t 0 ,
where α = 0.9 , σ = 0.4 , J = [ 0 , 4 ] , x ( t ) = ( x 1 ( t ) , x 2 ( t ) ) T , and
A = 0.3 0.1 0.15 0.2 , B = 0.2 0.1 0.15 0.25 , ϕ ( t ) = 1 1 2 ,
and
g ( t , x ( t ) , x ( t 0.4 ) ) = 1 3 e 2 t sin ( x 1 ( t ) ) + 1 4 e t sin 3 t arctan ( x 1 ( t 0.4 ) ) + 1 7 1 3 e 2 t cos ( x 2 ( t ) ) + 1 4 e t cos 3 t arctan ( x 2 ( t 0.4 ) ) + 1 6 ,
and
κ ( t , x ( t ) , x ( t 0.4 ) ) = 1 4 e t arctan ( x 1 ( t ) ) + 1 3 e 2 t cos 2 t sin ( x 1 ( t 0.4 ) ) + 1 3 1 4 e t sin ( x 2 ( t ) ) + 1 3 e t sin 2 t arctan ( x 2 ( t 0.4 ) ) + 1 6 ,
and
f ( t , x ( t ) , x ( t 0.4 ) , v ) = 1 2 1 3 .
For each x ( t ) ,   y ( t ) Υ and t [ 0 , T ] , we have
g ( t , x ( t ) , x ( t 0.4 ) ) g ( t , y ( t ) , y ( t 0.4 ) ) 1 3 | x 1 ( t ) y 1 ( t ) | + 1 4 | x 1 ( t 0.4 ) y 1 ( t 0.4 ) | + 1 3 | x 2 ( t ) y 2 ( t ) | + 1 4 | x 2 ( t 0.4 ) y 2 ( t 0.4 ) | 1 3 ( x ( t ) y ( t ) + x ( t 0.4 ) y ( t 0.4 ) ) .
Thus
g ( t , x ( t ) , x ( t 0.4 ) ) g ( t , y ( t ) , y ( t 0.4 ) ) 3 2 2 3 3 ( x ( t ) y ( t ) 3 + x ( t 0.4 ) y ( t 0.4 ) 3 ) ,
which implies that the function g satisfies the assumption (H1). Similarly, we can obtain that the functions κ and f satisfy the assumptions (H1) and (H2).
Let p = 3 . By calculation, we have γ =   A + B   = 0.8 , g L p = 0 4 g ( t , 0 , 0 ) 3 d t = 0.0651 , κ ( · , 0 , 0 ) = 1 2 , f ( · , 0 , 0 , 0 ) = 5 6 , C 1 = 4 27 and
Ξ = 0.4 0 ( C D 0 . 4 + 0.9 ϕ ) ( s ) A ϕ ( s ) 3 2 d s 2 2 0.4 0 ( ( C D 0 . 4 + 0.9 ϕ ) ( s ) 3 2 + A ϕ ( s ) 3 2 ) d s 2 2 0.4 0 1 Γ ( 0.1 ) 0.4 s ( s t ) 0.9 d t 1 2 Γ ( 0.1 ) 0.4 s ( s t ) 0.9 d t 3 2 + 0.35 0.25 3 2 2 2 0.4 0 3 2 Γ ( 1.1 ) 3 2 ( s + 0.4 ) 0.15 + 0 . 6 1.5 d s 2 = 1.5722 .
Hence, we may choose a suitable value μ > 0 such that
2 · 3 2 C 1 3 E 0.9 , 0.9 ( 0.8 · 4 0.9 ) 3 ( 4 2 + ( D 3 + 1 ) 2 + 1 ) Γ ( 0.7 ) < μ .
By Theorem 1, FSDDEs (39) have a unique solution x H 3 ( [ 0 , 4 ] ) .
Example 2. 
In the following, we consider the standard form of (39) as follows
( C D 0 0.9 x ϵ ) t ) = A x ϵ ( t ) + B x ϵ ( t 0.4 ) + ϵ g ( t , x ϵ ( t ) , x ϵ ( t 0.4 ) ) + ϵ κ ( t , x ϵ ( t ) , x ϵ ( t 0.4 ) ) d W ( t ) d t + ϵ V f ( t , x ϵ ( t ) , x ϵ ( t 0.4 ) , v ) N ¯ ( d t , d v ) , t J , x ϵ ( t ) = ϕ ( t ) , 0.4 t 0 ,
where x ϵ ( t ) = ( x 1 , ϵ ( t ) , x 2 , ϵ ( t ) ) T , and
g ( t , x ϵ ( t ) , x ϵ ( t 0.4 ) ) = 1 3 e 2 t sin ( x 1 , ϵ ( t ) ) + 1 4 e t sin 3 t arctan ( x 1 , ϵ ( t 0.4 ) ) + 1 7 1 3 e 2 t cos ( x 2 , ϵ ( t ) ) + 1 4 e t cos 3 t arctan ( x 2 , ϵ ( t 0.4 ) ) + 1 6 ,
and
κ ( t , x ϵ ( t ) , x ϵ ( t 0.4 ) ) = 1 4 e t arctan ( x 1 , ϵ ( t ) ) + 1 3 e 2 t cos 2 t sin ( x 1 , ϵ ( t 0.4 ) ) + 1 3 1 4 e t sin ( x 2 , ϵ ( t ) ) + 1 3 e 2 t sin 2 t arctan ( x 2 , ϵ ( t 0.4 ) ) + 1 6 ,
and
f ( t , x ϵ ( t ) , x ϵ ( t 0.4 ) , v ) = 1 2 1 3 .
We can easily check that the conditions (H1) and (H2) hold, and according to Theorem 1, FSDDSs (40) have a unique solution x ϵ given by
x ϵ ( t ) = X 0.4 , 0.9 , 1 A , B ( t + 0.4 ) ϕ ( 0.4 ) + 0.4 0 X 0.4 , 0.9 , 0.9 A , B ( t s ) [ C D 0 . 4 + 0.9 ϕ ) ( s ) A ϕ ( s ) ] d s + ϵ 0 t X 0.4 , 0.9 , 0.9 A , B ( t s ) g ( s , x ϵ ( s ) , x ϵ ( s 0.4 ) ) d s + ϵ 0 t X 0.4 , 0.9 , 0.9 A , B ( t s ) κ ( s , x ϵ ( s ) , x ϵ ( s 0.4 ) ) d W ( s ) + ϵ 0 t X 0.4 , 0.9 , 0.9 A , B ( t s ) V f ( s , x ϵ ( s ) , x ϵ ( s 0.4 ) , v ) N ¯ ( d s , d v ) .
By calculation, one has
g ^ ( x ε ( t ) , x ε ( t σ ) ) = lim t 1 t 0 t f ( s , x ϵ ( s ) , x ϵ ( s σ ) ) d s = 1 7 1 6 ,
κ ^ ( x ε ( t ) , x ε ( t σ ) ) = lim t 1 t 0 t σ ( s , x ε ( s ) , x ε ( s σ ) ) d s = 1 3 1 6 ,
f ^ ( x ε ( t ) , x ε ( t σ ) , v ) = lim t 1 t 0 t g ( s , x ε ( s ) , x ε ( s σ ) , v ) d s = 1 2 1 3 .
We are now checking that condition (H3) holds. In fact, one has
1 t 0 t g ( s , x ε ( s ) , x ε ( s σ ) ) g ^ ( x ε ( s ) , x ε ( s σ ) ) p d s = 1 t 0 t 1 3 e 2 s sin ( x 1 , ϵ ( s ) ) + 1 4 e s sin 3 s arctan ( x 1 , ϵ ( s 0.4 ) ) 1 3 e 2 s cos ( x 2 , ϵ ( s ) ) + 1 4 e s cos 3 s arctan ( x 2 , ϵ ( s 0.4 ) ) p d s = 1 3 p t 0 t ( e s ( | x 1 , ϵ ( s ) | + | x 2 , ϵ ( s ) | ) + e s ( | x 1 , ϵ ( s 0.4 ) | + | x 2 , ϵ ( s 0.4 ) | ) ) p d s = 1 3 p t 0 t ( e p s ( x ϵ ( s ) + x ϵ ( s 0.4 ) ) p d s 2 p 1 3 p t ( x ϵ ( s ) p + x ϵ ( s 0.4 ) p ) 0 t e p s d s = 2 p 1 ( 1 e p t ) 3 p p t ( 1 + x ϵ ( s ) p + x ϵ ( s 0.4 ) p ) ,
1 t 0 t ( t s ) α 1 ( κ ( s , x ϵ ( s ) , x ϵ ( s σ ) ) κ ^ ( x ϵ ( s ) , x ϵ ( s σ ) ) ) p d s = 1 t 0 t 1 4 ( t s ) α 1 e s arctan ( x 1 , ϵ ( s ) ) + 1 3 ( t s ) α 1 e 2 s cos 2 s sin ( x 1 , ϵ ( s 0.4 ) ) 1 4 ( t s ) α 1 e s sin ( x 2 , ϵ ( s ) ) + 1 3 ( t s ) α 1 e s sin 2 s arctan ( x 2 , ϵ ( s 0.4 ) ) p d s 2 p 1 3 p t ( x ϵ ( s ) p + x ϵ ( s 0.4 ) p ) 0 t ( t s ) p α p d s = 2 p 1 ( p α p + 1 ) 3 p t p α p ( 1 + x ϵ ( s ) p + x ϵ ( s 0.4 ) p ) ,
and
1 t 0 t V ( t s ) α 1 ( f ( s , x ϵ ( s ) , x ϵ ( s 0.4 ) , v ) f ^ ( x ϵ ( s ) , x ϵ ( s 0.4 ) , v ) p λ ( d v ) d s = 1 t 0 t V 1 2 ( t s ) α 1 1 3 ( t s ) α 1 p λ ( d v ) d s = 1 t 5 6 p λ ( V ) 0 t ( t s ) p ( α 1 ) d s 5 p λ ( V ) ( p α p + 1 ) 6 p t p α p ( 1 + x ϵ ( s ) p + x ϵ ( s 0.4 ) p ) .
Thus, (H3) is satisfied with
φ 1 ( t ) = 2 p 1 ( 1 e p t ) 3 p p t , φ 2 ( t ) = 2 p 1 ( p α p + 1 ) 3 p t p α p , and φ 3 ( t ) = 5 p λ ( V ) ( p α p + 1 ) 6 p t p α p .
Therefore, the conditions of Theorem 2 and Corollary 1 are satisfied. So, as ϵ 0 , the original solution x ϵ ( · ) y ϵ ( · ) in the sense of p square ( p = 3 ) and in the probability, where
y ϵ ( t ) = X 0.4 , 0.9 , 1 A , B ( t + 0.4 ) ϕ ( 0.4 ) + 0.4 0 X 0.4 , 0.9 , 0.9 A , B ( t s ) [ C D 0 . 4 + 0.9 ϕ ) ( s ) A ϕ ( s ) ] d s + ϵ 0 t X 0.4 , 0.9 , 0.9 A , B ( t s ) g ^ ( y ϵ ( s ) , y ϵ ( s 0.4 ) ) d s + ϵ 0 t X 0.4 , 0.9 , 0.9 A , B ( t s ) κ ^ ( y ϵ ( s ) , y ϵ ( s 0.4 ) ) d W ( s ) + ϵ 0 t X 0.4 , 0.9 , 0.9 A , B ( t s ) V f ^ ( y ϵ ( s ) , y ϵ ( s 0.4 ) , v ) N ¯ ( d s , d v ) .

6. Conclusions

In this article, we established and proved the existence and uniqueness theorem for solutions of Caputo-type fractional stochastic delay differential systems (FSDDSs) with Poisson jumps. By utilizing Burkholder–Davis–Gundy’s inequality, Doob’s martingale inequality, fractional Gronwall’s inequality, Hölder’s inequality, and Jensen’s inequality, we proved the averaging principle for FSDDSs in the sense of L p . This provides an effective stochastic approximation of the solutions of FSDDSs. Our method for fractional averaging will be beneficial for the study of the dynamics behavior of FSDDSs. Our results enrich the research field of fractional-order stochastic delay differential equations. Finally, we provided two examples to show the usefulness of our results.

Author Contributions

Conceptualization, Z.B. and C.B.; formal analysis, Z.B.; investigation, Z.B. and C.B.; and writing—review and editing, Z.B. and C.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Natural Science Foundation of China (11571136).

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors are grateful to the reviewers for their valuable comments and suggestions.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zhu, Q. Stabilization of stochastic nonlinear delay systems with exogenous disturbances and the event-triggered feedback control. IEEE Trans. Autom. Control. 2019, 64, 3764–3771. [Google Scholar] [CrossRef]
  2. Zhu, Q.; Huang, T. Stability analysis for a class of stochastic delay nonlinear systems driven by G-Brownian motion. Syst. Control Lett. 2020, 140, 104699. [Google Scholar] [CrossRef]
  3. Arthi, G.; Suganya, K. Controllability of higher order stochastic fractional control delay systems involving damping behavior. Appl. Math. Comput. 2021, 410, 126439. [Google Scholar] [CrossRef]
  4. Chen, Z.; Wang, B. Existence, exponential mixing and convergence of periodic measures of fractional stochastic delay reaction-diffusion equations on Rn. J. Differ. Equ. 2022, 336, 505–564. [Google Scholar] [CrossRef]
  5. Xu, L.; Li, Z. Stochastic fractional evolution equations with fractional brownian motion and infinite delay. Appl. Math. Comput. 2018, 336, 36–46. [Google Scholar] [CrossRef]
  6. Li, M.; Niu, Y.; Zou, J. A result regarding finite-time stability for Hilfer fractional stochastic differential equations with delay. Fractal Fract. 2023, 7, 622. [Google Scholar] [CrossRef]
  7. Khasminskii, R.Z. On the principle of averaging the Itö stochastic differential equations. Kibernetika 1968, 4, 260–279. [Google Scholar]
  8. Xu, Y.; Duan, J.Q.; Xu, W. An averaging principle for stochastic dynamical systems with Lévy noise. Phys. D 2011, 240, 1395–1401. [Google Scholar] [CrossRef]
  9. Xu, W.; Xu, W.; Zhang, S. The averaging principle for stochastic differential equations with Caputo fractional derivative. Appl. Math. Lett. 2019, 93, 79–84. [Google Scholar] [CrossRef]
  10. Luo, D.; Zhu, Q.; Luo, Z. An averaging principle for stochastic fractional differential equations with time-delays. Appl. Math. Lett. 2020, 105, 106290. [Google Scholar] [CrossRef]
  11. Ahmed, H.M.; Zhu, Q. The averaging principle of Hilfer fractional stochastic delay differential equations with Poisson jumps. Appl. Math. Lett. 2021, 112, 106755. [Google Scholar] [CrossRef]
  12. Ahmed, H.M. Impulsive conformable fractional stochastic differential equations with poisson jumps. Evol. Equ. Control. Theory 2022, 11, 2073–2080. [Google Scholar] [CrossRef]
  13. Wang, Z.; Lin, P. Averaging principle for fractional stochastic differential equations with Lp convergence. Appl. Math. Lett. 2022, 130, 108024. [Google Scholar] [CrossRef]
  14. Xu, W.; Xu, W.; Lu, K. An averaging principle for stochastic differential equations of fractional order 0 < α < 1. Fract. Calc. Appl. Anal. 2020, 23, 908–919. [Google Scholar]
  15. Yang, D.; Wang, J.; Bai, C. Averaging principle for ψ-Capuo fractional stochastic delay differential equations with Poisson jumps. Symmetry 2023, 15, 1346. [Google Scholar] [CrossRef]
  16. Li, M.; Wang, J. The existence and averaging principle for Caputo fractional stochastic delay differential systems. Fract. Calc. Appl. Anal. 2023, 26, 893–912. [Google Scholar] [CrossRef]
  17. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  18. Mahmudov, N.I. Delayed perturbation of Mittag–Leffler functions and their applications to fractional linear delay differential equations. Math. Meth. Appl. Sci. 2019, 42, 5489–5497. [Google Scholar] [CrossRef]
  19. You, Z.; Fečkan, M.; Wang, J. Relative controllability of fractional delay differential equations via delayed perturbation of Mittag–Leffler functions. J. Comput. Appl. Math. 2020, 378, 112939. [Google Scholar] [CrossRef]
  20. Son, D.; Huong, P.; Kloeden, P.; Tuan, H. Asymptotic separation between solutions of Caputo fractional stochastic differential equations. Stoch. Anal. Appl. 2018, 36, 654–664. [Google Scholar] [CrossRef]
  21. Applebaum, D. Lévy Process and Stochastic Calculus; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  22. Kunita, H. Stochastic differential equations based on Lévy processes and stochastic flows of diffeomorphisms. In Real and Stochastic Analysis, New Perspectives; Birkhauser: Basel, Switzerland, 2004; pp. 305–373. [Google Scholar]
  23. Ye, H.; Gao, J.; Ding, Y. A generalized Gronwall inequality and its application to a fractional differential equation. J. Math. Anal. Appl. 2007, 328, 1075–1081. [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.

Share and Cite

MDPI and ACS Style

Bai, Z.; Bai, C. The Existence and Averaging Principle for Caputo Fractional Stochastic Delay Differential Systems with Poisson Jumps. Axioms 2024, 13, 68. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms13010068

AMA Style

Bai Z, Bai C. The Existence and Averaging Principle for Caputo Fractional Stochastic Delay Differential Systems with Poisson Jumps. Axioms. 2024; 13(1):68. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms13010068

Chicago/Turabian Style

Bai, Zhenyu, and Chuanzhi Bai. 2024. "The Existence and Averaging Principle for Caputo Fractional Stochastic Delay Differential Systems with Poisson Jumps" Axioms 13, no. 1: 68. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms13010068

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