Next Article in Journal
Stability Analysis for a Fractional-Order Coupled FitzHugh–Nagumo-Type Neuronal Model
Previous Article in Journal
New Results Involving Riemann Zeta Function Using Its Distributional Representation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Averaging Principle for Stochastic Fractional Differential Equations Driven by fBm Involving Impulses

1
School of Applied Science, Taiyuan University of Science and Technology, Taiyuan 030024, China
2
School of Mathematics and Statistics, Northwestern Polytechnical University, Xi’an 710129, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(5), 256; https://doi.org/10.3390/fractalfract6050256
Submission received: 2 April 2022 / Revised: 1 May 2022 / Accepted: 3 May 2022 / Published: 7 May 2022

Abstract

:
In contrast to previous research on periodic averaging principles for various types of impulsive stochastic differential equations (ISDEs), we establish an averaging principle without periodic assumptions of coefficients and impulses for impulsive stochastic fractional differential equations (ISFDEs) excited by fractional Brownian motion (fBm). Under appropriate conditions, we demonstrate that the mild solution of the original equation is approximately equivalent to that of the reduced averaged equation without impulses. The obtained convergence result guarantees that one can study the complex system through the simplified system. Better yet, our techniques dealing with multi-time scales and impulsive terms can be applied to improve some existing results. As for application, three examples are worked out to explain the procedure and validity of the proposed averaging principles.

1. Introduction

In the past few decades, fractional calculus and fractional dynamical systems have gained widespread concern from scholars in various research fields. Different researchers explored fractional dynamical systems from both theoretical and practical aspects due to their applications and prospects in many scientific areas [1,2,3]. Specifically, taking into account the inevitable uncertainties and random factors, massive research achievements on stochastic fractional differential equations have emerged one after another [4,5,6]. Besides, some real-world systems are affected by instantaneous disturbances or undergo sudden changes, and impulsive differential equations are naturally used to characterize them. A wide range of mathematical models in the study of control theory, telecommunications, biology, ecology, epidemiology, finance and economics can be described by integer or fractional-order impulsive stochastic differential equations (ISDEs) [7,8].
As we know, random perturbations with LRD (long-range dependence) exist in various application fields; for instance, hydrology, meteorology, physical chemistry, bioengineering, mathematical finance, etc. [9]. Fractional Brownian motion has a long correlation time and can be used to model these stochastic excitations. Therefore, various types of stochastic fractional differential equations (SFDEs) excited by fBm are being emphasized and some preliminary study results [10,11,12] have been acquired. Meanwhile, some recent research [13,14] considered the well-posedness and controllability of different classes of ISFDEs. In particular, Pedjeu and Ladde [15] introduce a stochastic model under multi-time scales, which can model the complex dynamical processes in sciences, engineering, ecology and epidemiology. Correspondingly, considering that the correlated noises can not be ignored, Abouagwa, Cheng and Li [16] establish the existence and uniqueness results of SFDEs under fBm involving impulses based on Carathéodory successive approximation method.
On the other hand, the averaging method is an effective approach to explore different types of non-linear dynamical systems. This is because you can concentrate on the reduced autonomous averaged equation instead of the original complex time-varying one, providing an approximate way to remove some of the complexity. The core of demonstrating an averaging principle is examining the conditions when the averaged system is equivalent to the original system in some sense. Since Khasminskii’s contribution [17], stochastic averaging methods for SDEs have been widely developed and applied [18,19]. However, most of the noises they had considered are uncorrelated. Such results ruled out these stochastic differential systems under the random perturbations with long-term dependence. Fortunately, Pei, Xu and Guo [20,21] have recently obtained the averaging principles for several different kinds of SDEs under fBm with Hurst index H ( 1 2 , 1 ) . With the development of research on ISDEs, in two recent years, a few study achievements on periodic averaging methods of ISDEs have appeared [22,23,24,25]. There are two disadvantages of this research. One is that the periodic assumptions of the coefficient functions and the impulses are strict; the other is the estimate of the difference between the impulsive term of the original equation and the impulsive integral term of the averaged equation is flawed. Here, it is noted that Liu and Xu [26] proved an averaging theorem for impulsive stochastic partial differential equations without periodic assumptions. In light of complex characteristics of stochastic fractional differential systems, Xu et al. [27] developed stochastic averaging for several classes of SFDEs driven by different noises. Abouagwa and Li [28] considered an averaging method for SFDEs under a non-Lipschitz condition. Luo and Zhu [29] discussed an averaging principle for SFDEs with delay. Shen et al. [30] showed an averaging principle for SFDEs with Lévy noise and Markovian switching. Guo et al. [31] examined an averaging theorem under a weaker condition for a kind of Caputo SFDEs. Liu and Xu [32] proved an averaging principle for neutral ISFDEs. These works offer posible ways to reduce various types of SFDEs. However, these papers do not involve the impulsive effects and the random perturbations with long-range dependence with which we are concerned. Another gap is that they showed no uniformity of time scales. To be precise, the aforementioned results cannot solve the problem if the SFDEs are excited by fBm and involve impulses. In this work, our main objective is to fill these two gaps. We will establish an averaging principle for the SFDEs driven by fBm with impulses in the following form:
d x t = f t , x t d t + g t , x t ( d t ) β + h t , x t d B H ( t ) , t [ 0 , T ] , t t k , Δ x t k = x ( t k + ) x ( t k ) = I k x t k , t = t k , k N , x 0 = x 0 R d ,
where f , g : [ 0 , T ] × R d R d ,   h : [ 0 , T ] × R d R d × R m are measurable functions; the processes B H ( t ) stand for the independent fBms, in which the Hurst parameter H ( 1 2 , 1 ) ,   0 < β < 1 ,   0 < T < ,   I k : R d R d indicate the pulse size of x at time t k , with x t k + = lim τ 0 + x ( t k + τ ) , x t k = lim τ 0 x ( t k + τ ) , x ( t k ) = x t k + ; and the impulsive time sequence satisfies that 0 t 1 < t 2 < < t k < t k + 1 < with lim k t k = . The initial value x 0 R d is a given random variable with E x 0 2 < .
Since Equation (1) is introduced under multi-time scales, the time scale of the standard form of system (1) should be uniform in the process of deriving the averaging principle. However, this point was overlooked in the present literature [28,29,30]. Additionally, the involvement of impulses makes the proof of the main results more technical. The techniques shown in the aforementioned research on periodic averaging methods for ISDEs are not applicable to our case. The innovation of this paper is not only that an averaging theorem for a new class of SFDEs is discussed, but also that some published results are improved. Naturally, we highlight the contributions of this article here.
  • We propose an effective approximation for the original system (1); hence, the complexity of ISFDEs under fBm can be reduced.
  • The problem of no match on each of the time scales of the standard stochastic fractional differential equations is pointed out and corrected.
  • The obtained averaging principle is valid for stochastic fractional differential equations driven by fBm; that is, our results are new even for non-impulsive SFDEs with fBm.
  • We present a method to estimate the impulsive terms, which is helpful to develop averaging principles for different types of ISDEs.
The rest of this work is as follows: Some essential notations and preliminaries are presented in Section 2. In Section 3, an averaging principle for system (1) is established under a Non-Lipschitz condition, and some useful remarks are also presented. In the end, Section 4 shows three examples to explain the procedure and validity of the averaging theorem for SFDEs and ISFDEs.

2. Preliminary

Let ϕ ( r , t ) = H ( 2 H 1 ) r t 2 H 2 , in which H ( 1 12 2 , 1 ) is a constant. Defining
L ϕ 2 ( R + ) = p : p ϕ 2 = R + R + p ( r ) p ( t ) ϕ ( r , t ) d t d r < ,
where p : R + R is a Borel measurable function; L ϕ 2 ( R + ) comes to a separable Hilbert space when it is endowed with the following scalar product:
p 1 , p 2 ϕ = R + R + p 1 ( r ) p 2 ( t ) ϕ ( r , t ) d t d r , p 1 , p 2 L ϕ 2 ( R + ) .
Denote S as the family of cylindrical and smooth random variables:
F ( ω ) = p B H ( ψ 1 ) , B H ( ψ 2 ) , , B H ( ψ n ) ,
in which p C b ( R n ) , n 1 and ψ i H , where H represents the complete set of all measurable mappings with ψ φ 2 < with ψ i , ψ j φ = δ i j . Now, we can define the space H of p:
p H 2 = 0 T 0 T p ( t ) p ( s ) φ ( t , s ) d s d t < .
The Malliavin derivative D t H of F S is determined as follows:
D t H F = i = 1 n p x i B H ( ψ 1 ) , B H ( ψ 2 ) , , B H ( ψ n ) ψ i ( t ) .
Moreover, for any positive integer k, D t H , k symbolizes the iteration of the above formula. Then D k , q represents the closure of S and is normed by:
F k , q q = E F q + E j = 1 k D t H , j F H j q ,
in which q 1 , stands for the tensor product. Likewise, let us introduce D 1 , q ( H ) as the subspace of D 1 , q ( H ) . Then, the Malliavin ϕ -derivative of F can be given: D t ϕ F = R + ϕ ( t , s ) D s H F d s .
There are three types of pathwise integrals definitions for fBm (i.e., forward pathwise integral, symmetric pathwise integral and backward pathwise integral) [33]. In this paper, we use the forward pathwise integral of fBm. For more details on fBm and stochastic integrals of fBm, one can refer [20,21,33,34].
Definition 1
([33]). Assume that y ( t ) is a stochastic process on [ 0 , T ] with integrable trajectories, H ( 0 , 1 ) . If the limit of
lim δ 0 0 T y ( t ) B H ( t + δ ) B H ( t ) δ d t
exists in probability, then the above expression is called as the forward integral of y ( t ) for B H ( t ) , denoted as 0 T y ( t ) d B H ( t ) .
Remark 1.
Defining L ϕ [ 0 , T ] as the space of stochastic processes y ( t ) satisfying E y ( t ) ϕ 2 < , in which y ( t ) is ϕ-differentiable, and D s ϕ y ( t ) meet that E 0 T 0 T [ D t φ y ( s ) ] 2 d s d t < . More details about L φ [ 0 , T ] are given in [21].
Lemma 1
([20]). Let y ( t ) L ϕ [ 0 , T ] D 1 , 2 ( H ) , and B H ( t ) is a fBm with H ( 1 2 , 1 ) , then for any 0 < T < , there exists a positive constant C, such that
E 0 T y ( s ) d B H ( s ) 2 2 H T 2 H 1 E 0 T y ( s ) 2 d s + 4 T E 0 T [ D s φ y ( s ) ] 2 d s 2 H T 2 H 1 E 0 T y ( s ) 2 d s + 4 C T 2 .
Definition 2
([15]). Given β ( 0 , 1 ] , g ( t ) is a continuous function, then the integral of g ( t ) with ( d t ) β is determined as follows:
0 t g ( τ ) ( d τ ) β = β 0 t ( t τ ) β 1 g ( τ ) d τ .
Following Refs. [16,35], we show the definition of the mild solution for system (1).
Definition 3.
A stochastic process x ( t ) defined on [ 0 , T ] is called a mild solution of (1) with initial value x ( 0 ) = x 0 , if
(i)
x ( t ) R d is F t -adapted, and has càdlàg path a.e. on [ 0 , T ] ;
(ii)
for all t [ 0 , T ] , x ( t ) meets the integral equation below:
x ( t ) = x ( 0 ) + 0 t f ( s , x ( s ) ) d s + β 0 t g ( s , x ( s ) ) ( t s ) 1 β d s + 0 t h ( s , x ( s ) ) d B H ( s ) + 0 < t k < t I k ( x ( t k ) ) .
To achieve the approximate theoretical result, the following hypotheses are supposed throughout this work.
Condition 1.
f ( t , · ) , g ( t , · ) , h ( t , · ) L ϕ [ 0 , T ] D 1 , 2 ( H ) satisfy that
f ( t , x ) f ( t , y ) 2 + g ( t , x ) g ( t , y ) 2 + h ( t , x ) h ( t , y ) 2 + D t ϕ ( h ( t , x ) h ( t , y ) ) 2 T ( t , x y 2 ) , t [ 0 , T ] ,
where T : 0 + × R + R + is locally integrable with t, given any fixed t 0 , T is a continuous, non-decreasing, concave function with respect to x, T ( t , 0 ) = 0 and 0 + 1 T ( t , x ) d x = + . Additionally, for any t R + and positive constant κ , if the non-negative continuous function Y ( t ) fulfills that Y ( t ) κ 0 t T ( s , Y ( s ) ) d s , we have Y ( t ) 0 .
Condition 2.
f ( t , · ) , g ( t , · ) , h ( t , · ) L ϕ [ 0 , T ] D 1 , 2 ( H ) satisfy that
f ( t , x ) 2 + g ( t , x ) 2 + h ( t , x ) 2 + D t ϕ h ( t , x ) 2 M ( t , x 2 ) , t [ 0 , T ] ,
where M : 0 + × R + R + is locally integrable for t and given every fixed t, M is continuous, monotone non-decreasing, concave with x. Moreover, given any positive constant κ and initial data x 0 , the integral formula x ( t ) = x 0 + κ 0 t M ( s , x ( s ) ) d s admits a global solution.
Condition 3.
The impulsive functions I k satisfy the following Lipschitz condition and boundedness condition: there are positive constants c k , d k for all x , y R d , such that
I k ( x ) I k ( y ) 2 c k x y 2 , I k ( x ) 2 d k .
Remark 2.
Following Theorem 3.1 in Ref. [16], the system (1) has a unique mild solution under Conditions 1–3. Additionally, the unique solution of (1) has the property that
E sup 0 t T x ( t ) 2 C ,
where C is a positive constant.

3. Main Results

In this part, we devote to develop an averaging principle for ISFDEs excited by fBm. That is, an approximate way of simplifying the system is presented.
The standard form of ISFDEs driven by fBm is defined as
x ε ( t ) = x ( 0 ) + ε 0 t f ( s , x ε ( s ) ) d s + ε β β 0 t g ( s , x ε ( s ) ) ( t s ) 1 β d s + ε H 0 t h ( s , x ε ( s ) ) d B H ( s ) + ε 0 < t k < t I k ( x ε ( t k ) ) ,
where t [ 0 , T ] , ε 0 , ε 1 ( 0 < ε 1 1 ) represents a time-scale parameter; the coefficients f , g , h , I k meet Conditions 1–3.
Remark 3.
As discussed in Ref. [32], β 0 t g ( s , x ε ( s ) ) ( t s ) 1 β d s is gained from the integral of g ( t , x ε ( t ) ) with ( d t ) β , so the time scale of this term should be ε β .
Next, we will demonstrate that when ε tends to zero, the solution x ε ( t ) for Equation (1) converges to the solution z ε ( t ) of the following averaged equation:
z ε ( t ) = x ( 0 ) + ε 0 t f ¯ ( z ε ( s ) ) + I ¯ ( z ε ( s ) ) d s + ε β β 0 t g ¯ ( z ε ( s ) ) ( t s ) 1 β d s + ε H 0 t h ¯ ( z ε ( s ) ) d B H ( s ) ,
the coefficients f ¯ ( · ) , g ¯ ( · ) , h ¯ ( · ) D 1 , 2 ( H ) L φ [ 0 , T ] and I ¯ : R d R d are all measurable, while meeting Conditions 1–3 and the hypotheses below:
Hypothesis 1 (H1):
1 S 0 S f ( s , x ) f ¯ ( x ) 2 d s γ 1 ( S ) ϑ ( x 2 ) ,
Hypothesis 2 (H2):
1 S 0 S g ( s , x ) g ¯ ( x ) 2 d s γ 2 ( S ) ϑ ( x 2 ) ,
Hypothesis 3 (H3):
1 S 0 S h ( s , x ) h ¯ ( x ) 2 d s γ 3 ( S ) ϑ ( x 2 ) ,
Hypothesis 4 (H4):
1 S 0 < t j < S I j ( x ) S I ¯ ( x ) γ 4 ( S ) ( 1 + x ) ,
where S [ 0 , T ] , γ i ( S ) are positive bounded functions with lim S γ i ( S ) = 0 , i = 1 , 2 , 3 , 4 , and the function ϑ : R + R + is non-decreasing, continuous and concave.
Remark 4.
As in Remark 2, under Hypotheses H 1 H 4 , it is not difficult to derive that E sup 0 t T z ε ( t ) 2 < .
Now, the main theorem of this paper is presented as follows.
Theorem 1.
Let Conditions 1–3 and Hypotheses H 1 H 4 hold. For any given small positive number δ, there are P > 0 ,   ε 2 0 , ε 1 and γ 0 , 1 such that, for all 0 < ε ε 2 ,   0 t P ε H γ , we have
E x ε ( t ) z ε ( t ) 2 δ .
For the sake of showing Theorem 1, we state the following lemma in advance.
Lemma 2.
Let the conditions of Theorem 1 hold; then, for any 0 s t T , we have
E z ε ( t ) z ε ( s ) 2 ε 2 Q 1 ( t s ) 2 + ε 2 β Q 2 ( t s ) 2 β + ε 2 H Q 3 ( t s ) 2 H + ε 2 H Q 4 ( t s ) 2 ,
in which Q i ( i = 1 , 2 , 3 , 4 ) are four positive constants.
The proof of Lemma 2 is analogous to that of Theorem 3.1 in [16], so we just outline it here.
From Equation (4), with basic inequality, Cauchy–Schwarz inequality and Lemma 1, one can get
E z ε ( t ) z ε ( s ) 2 3 ε 2 ( t s ) s t f ¯ ( z ε ( θ ) ) + I ¯ ( z ε ( θ ) ) 2 d θ + 3 ε 2 β β 2 ( t s ) 2 β 1 2 β 1 s t g ¯ ( z ε ( θ ) ) 2 d θ + 3 ε 2 H 2 H ( t s ) 2 H 1 E s t h ¯ ( z ε ( θ ) 2 d θ + 4 C ( t s ) 2 ,
then, with the aid of Condition 2, H 1 H 4 and Remark 4, one can conclude the statement of Lemma 2 without difficulty.
We are now turning back to the proof of Theorem 1.
Proof. 
Employing the elementary inequality, we obtain from (3) and (4):
E x ε ( t ) z ε ( t ) 2 4 ε 2 E 0 t f ( s , x ε ( s ) ) f ¯ ( z ε ( s ) ) d s 2 + 4 β 2 ε 2 β E 0 t ( g ( s , x ε ( s ) ) g ¯ ( z ε ( s ) ) ) ( t s ) 1 β d s 2 + 4 ε 2 H E 0 t h ( s , x ε ( s ) ) h ¯ ( z ε ( s ) ) d B H ( s ) 2 + 4 ε 2 E 0 < t k < t I k ( x ε ( t k ) ) 0 t I ¯ ( z ε ( s ) ) d s 2 : = J 1 + J 2 + J 3 + J 4 .
Denote by 0 , t 0 , u ˜ 0 , T . Next, we estimate each term separately. For the first term, we have
J 1 8 ε 2 E 0 u ˜ f ( s , x ε ( s ) ) f ( s , z ε ( s ) ) d s 2 + 8 ε 2 E 0 u ˜ f ( s , z ε ( s ) ) f ¯ ( z ε ( s ) ) d s 2 : = J 11 + J 12 .
With the virtue of Cauchy–Schwarz inequality and Condition 1, it follows that
J 11 8 ε 2 u ˜ E 0 u ˜ T ( s , x ε ( s ) z ε ( s ) 2 ) d s .
On the basis of Remark 4, the elementary inequality, Condition 1 and H 1 , one has
J 12 8 ε 2 u ˜ 2 E 1 u ˜ 0 u ˜ f ( s , z ε ( s ) ) f ¯ ( z ε ( s ) ) 2 d s 8 ε 2 u ˜ 2 C 12 ,
where C 12 is a positive constant. Then, we arrive at
J 1 8 ε 2 u ˜ E 0 u ˜ T ( s , x ε ( s ) z ε ( s ) 2 ) d s + 8 ε 2 u ˜ 2 C 12 .
Adopting Cauchy–Schwarz inequality to J 2 , we get
J 2 4 β 2 ε 2 β 0 t ( t s ) 2 β 1 d s E 0 t g ( s , x ε ( s ) ) g ¯ ( z ε ( s ) ) 2 d s 4 ε 2 β β 2 t 2 β 1 2 β 1 E 0 t g ( s , x ε ( s ) ) g ¯ ( z ε ( s ) ) 2 d s ,
for E 0 t g ( s , x ( s ) ) g ¯ ( z ε ( s ) ) 2 d s , estimating analogously as above, one can easily obtain
J 2 8 β 2 ε 2 β u ˜ 2 β 1 2 β 1 E 0 t T ( s , x ε ( s ) z ε ( s ) 2 ) d s + 8 β 2 ε 2 β u ˜ 2 β 2 β 1 sup 0 t u ˜ γ 2 ( t ) E ϑ sup 0 t u ˜ z ε ( t ) 2 8 β 2 ε 2 β u ˜ 2 β 1 2 β 1 E 0 u ˜ T ( s , x ε ( s ) z ε ( s ) 2 ) d s + 8 β 2 ε 2 β u ˜ 2 β 2 β 1 C 22 ,
where 1 2 < β < 1 , C 22 is a positive constant.
To proceed, thanks to Lemma 1, the elementary inequality, Condition 1, Hypothesis H 2 and Remark 4, one can derive that
J 3 8 ε 2 H H t 2 H 1 E 0 u ˜ h ( s , x ε ( s ) ) h ¯ ( z ε ( s ) ) 2 d s + 16 ε 2 H C t 2 16 ε 2 H H t 2 H 1 E 0 u ˜ h ( s , x ε ( s ) ) h ( s , z ε ( s ) ) 2 d s + 16 ε 2 H H t 2 H 1 E 0 u ˜ h ( s , z ε ( s ) ) h ¯ ( z ε ( s ) ) 2 d s + 16 ε 2 H C t 2 16 ε 2 H u ˜ 2 H 1 H E 0 u ˜ T ( s , x ε ( s ) z ε ( s ) 2 ) d s + 16 ε 2 H H u ˜ 2 H sup 0 t u ˜ γ 3 ( t ) E ϑ sup 0 t u ˜ z ε ( t ) 2 + 16 ε 2 H C u ˜ 2 16 ε 2 H u ˜ 2 H 1 H E 0 u ˜ T ( s , x ε ( s ) z ε ( s ) 2 ) d s + 16 ε 2 H u ˜ 2 H 1 H C 32 + 16 ε 2 H C u ˜ 2 ,
where C 32 is a positive constant.
For the last term J 4 , we have
J 4 16 ε 2 E 0 < t k < t I k ( x ε ( t k ) ) 0 t I ¯ ( x ε ( t k ) ) d s 2 + 16 ε 2 E 0 t I ¯ ( x ε ( t k ) ) d s 0 t I ¯ ( z ε ( t k ) ) d s 2 + 16 ε 2 E 0 t I ¯ ( z ε ( t k ) ) d s 0 t I ¯ ( z ε ( s ) ) d s 2 : = J 41 + J 42 + J 43 .
Notice that
0 < t k < t I k ( x ε ( t k ) ) = i = 1 m 0 < t k < t i + 1 I k ( x ε ( t k ) ) 0 < t k < t i I k ( x ε ( t k ) )
and
0 t I ¯ ( x ε ( t k ) ) d s = i = 0 m ( t i + 1 t i ) I ¯ ( x ε ( t k ) ) .
Estimating J 41 in the following manner:
J 41 16 ε 2 E i = 0 m t i + 1 1 t i + 1 0 < t k < t i + 1 I k ( x ε ( t k ) ) I ¯ ( x ε ( t k ) ) + t i 1 t i 0 < t k < t i I k ( x ε ( t k ) ) I ¯ ( x ε ( t k ) ) 2 ,
then Hypothesis H 4 yields that
J 41 32 ε 2 ( m + 1 ) u ˜ 2 sup 0 t u ˜ γ 4 ( t ) 2 2 E 1 + x ε ( t k ) 2 64 ε 2 ( m + 1 ) u ˜ 2 sup 0 t u ˜ γ 4 ( t ) 2 1 + E sup 0 t u ˜ x ε ( t ) 2 64 ε 2 ( m + 1 ) u ˜ 2 C 41 ,
where C 41 is a positive constant.
As I ¯ meets the Lipschitz condition, there is a positive constant N involving with c k such that
I ¯ ( x ) I ¯ ( y ) 2 N x y 2 , x , y R d .
Using Cauchy-Schwarz inequality, we have
J 42 16 ε 2 u ˜ N 0 u ˜ E x ε ( t k ) z ε ( t k ) 2 d s 16 ε 2 u ˜ N 0 u ˜ sup 0 θ s E x ε ( θ ) z ε ( θ ) 2 d s .
With the aid of Lemma 2, one can reach
J 43 16 ε 2 u ˜ N 0 t E z ε ( t k ) z ε ( s ) 2 d s 16 ε 4 N Q 1 u ˜ 4 + 16 ε 2 β + 2 N Q 2 u ˜ 2 β + 2 + 16 ε 2 H + 2 N Q 3 u ˜ 2 H + 2 + 16 ε 2 H + 2 N Q 4 u ˜ 4 .
Substituting (11), (12) and (13) together into (10), we get
J 4 64 ε 2 ( m + 1 ) u ˜ 2 C 41 + 16 ε 2 u ˜ N 0 u ˜ sup 0 θ s E x ε ( θ ) z ε ( θ ) 2 d s + 16 ε 4 N Q 1 u ˜ 4 + 16 ε 2 β + 2 N Q 2 u ˜ 2 β + 2 + 16 ε 2 H + 2 N Q 3 u ˜ 2 H + 2 + 16 ε 2 H + 2 N Q 4 u ˜ 4 .
Thanks to the fact that T is concave in x, there is w ( t ) > 0 and e ( t ) > 0 such that
T ( t , x ) w ( t ) + e ( t ) x , 0 T w ( t ) d t < , 0 T e ( t ) d t < ,
then, combining the estimates (7), (8), (9) and (14), we obtain
sup 0 t u ˜ E x ε ( t ) z ε ( t ) 2 8 ε 2 u ˜ + 16 ε 2 H u ˜ 2 H 1 + 8 ε 2 β u ˜ 2 β 1 2 β 1 C 1 0 u ˜ sup 0 θ s E x ε ( θ ) z ε ( θ ) 2 d s + 16 ε 2 N u ˜ 0 u ˜ sup 0 θ s E x ε ( θ ) z ε ( θ ) 2 d s + 80 ε 2 u ˜ 2 C 2 + 32 ε 2 H u ˜ 2 H C 3 + 64 ε 2 H C 4 u ˜ 2 + 16 ε 2 β u ˜ 2 β C 5 2 β 1 + 16 ε 4 C 6 u ˜ 4 + 16 ε 2 β + 2 C 7 u ˜ 2 β + 2 + 16 ε 2 H + 2 C 8 u ˜ 2 H + 2 + 16 ε 2 H + 2 C 9 u ˜ 4 ,
where C i , i = 1 , 2 , , 9 are positive constants. Gronwall inequality shows that there are P > 0 and γ 0 , 1 fulfilling
sup t 0 , P ε H γ E x ε ( t ) z ε ( t ) 2 Q ε ( 1 H γ ) ,
where
Q = 80 C 2 P 2 ε 1 H γ + 32 C 3 P 2 H ε ( 2 H 1 ) ( 1 H γ ) + 64 C 4 P 2 ε ( 2 H 1 ) ( 1 γ ) + 16 C 5 P 2 β ε ( 2 β 1 ) ( 1 H γ ) 16 C 5 P 2 β ε ( 2 β 1 ) ( 1 H γ ) ( 2 β 1 ) ( 2 β 1 ) + 16 C 6 P 4 ε 3 ( 1 H γ ) + 16 C 7 P 2 β + 2 ε ( 2 β + 1 ) ( 1 H γ ) + 16 C 8 P 2 H + 2 ε ( 2 H + 1 ) ( 1 H γ ) + 16 C 9 P 4 ε 2 H ( 1 γ ) + ( 1 H γ ) exp 8 P 2 C 1 ε 2 ( 1 H γ ) + 16 P 2 H C 1 ε 2 H ( 1 H γ ) + 8 P 2 β C 1 ε 2 β ( 1 H γ ) 8 P 2 β C 1 ε 2 β ( 1 H γ ) ( 2 β 1 ) ( 2 β 1 ) + 16 N P 2 ε 2 ( 1 H γ )
is a positive constant.
Therefore, for arbitrarily given small positive number δ , there is ε 2 0 , ε 1 such that
E x ε ( t ) z ε ( t ) 2 δ ,
for ε 0 , ε 2 , t 0 , P ε H γ 0 , T . □
Remark 5.
Noted that the order of convergence interval is ε H γ in this study. It is not stronger as ε γ of Bm case and Lévy case [27,28,29,30] because of the weakness that stochastic integral of fBm is not a martingale. However, the proposed result is still certainly a good theoretical method to simplify non-autonomous ISFDEs driven by fBm with long-term dependence.
Remark 6.
When I k = 0 , system (1) degenerates into stochastic fractional differential equation driven by fBm. Evidently, our obtained result is applicable for SFDEs with fBm.

4. Example

Example 1.
Consider the SFDE under fBm below:
d x ε ( t ) = ε sin 2 ( t ) x ε ( t ) d t + ε β ( d t ) β + ε H d B H ( t ) , t [ 0 , T ] x ε ( 0 ) = x 0 ,
where E x 0 2 < , f ( t , x ε ( t ) ) = sin 2 ( t ) x ε ( t ) , g ( t , x ε ( t ) ) = 1 , h ( t , x ε ( t ) ) = 1 . Let
f ¯ ( z ε ( t ) ) = 1 2 z ε ( t ) , g ¯ ( z ε ( t ) ) = 1 , h ¯ ( z ε ( t ) ) = 1 .
Obviously, Conditions 1, 2 and Hypotheses H 1 H 3 are satisfied. The simplified averaged SFDE is given by
d z ε ( t ) = 1 2 ε z ε ( t ) d t + ε β ( d t ) β + ε H d B H ( t ) .
Remark 6 illustrates that z ε ( t ) is approximately equivalent to x ε ( t ) on 0 , T . Then, a numerical comparison between the solution x ε ( t ) of original Equation (15) and the solution z ε ( t ) of the averaged Equation (16) is presented. Not surprisingly, good agreement can be observed in Figure 1.
Example 2.
Consider the SFDE with fBm below:
d x ε ( t ) = 2 ε sin 2 ( t ) x ε ( t ) d t + 2 ε β cos 2 ( t ) ( d t ) β + ε H d B H ( t ) , t [ 0 , T ] x ε ( 0 ) = x 0 ,
where E x 0 2 < , f ( t , x ε ( t ) ) = 2 sin 2 ( t ) x ε ( t ) , g ( t , x ε ( t ) ) = 2 cos 2 ( t ) , h ( t , x ε ( t ) ) = 1 . Let
f ¯ ( z ε ( t ) ) = z ε ( t ) , g ¯ ( z ε ( t ) ) = 1 , h ¯ ( z ε ( t ) ) = 1 .
One can easily verify that Conditions 1, 2 and Hypotheses H 1 H 3 are satisfied. Meanwhile, a simplified averaged SFDE is presented as
d z ε ( t ) = ε z ε ( t ) d t + ε β ( d t ) β + ε H d B H ( t ) .
Theorem 1 guarantees that z ε ( t ) can approximate x ε ( t ) on 0 , P ε H γ 0 , T . In succession, we show a numerical comparison between the solution x ε ( t ) to the original Equation (17) and the solution z ε ( t ) to the averaged Equation (18). It can be found that the two are consistent in a certain range in Figure 2.
Example 3.
Consider the ISFDE with fBm below:
d x ε ( t ) = ε d t + 2 ε β sin 2 ( t ) ( d t ) β + ε H cos 2 ( t ) x ε ( t ) d B H ( t ) , t t k , Δ x ε ( t k ) = x ε ( t k + ) x ε ( t k ) = ε arctan ( x ε ( t k ) ) , t = t k , x ε ( 0 ) = x 0 ,
where E x 0 2 < , f ( t , x ε ( t ) ) = 1 , g ( t , x ε ( t ) ) = 2 sin 2 ( t ) , h ( t , x ε ( t ) ) = cos 2 ( t ) x ε ( t ) . Let
f ¯ ( z ε ( t ) ) = 1 , g ¯ ( z ε ( t ) ) = 1 , h ¯ ( z ε ( t ) ) = 1 2 z ε ( t ) , I ¯ ( z ε ( t ) ) = 1 5 arctan ( z ε ( t ) ) ,
then, a new simplified averaged SFDE without impulses is given as
d z ε ( t ) = ε 1 + 1 5 arctan ( z ε ( t ) ) d t + ε β ( d t ) β + 1 2 ε H z ε ( t ) d B H ( t ) .
As we can see, the impulses are averaged out and the original system is simplified. Even better, only a small difference is generated in this proximate process. In the end, we carry out a numerical comparison between the solution process x ε ( t ) of original Equation (19) and the solution process z ε ( t ) of the simplified Equation (20). As expected, good agreement can be seen in Figure 3.

Author Contributions

Formal analysis, J.L.; methodology, J.L. and W.W.; writing—original draft preparation, J.L.; writing—review and editing, J.L. and W.W.; supervision, W.X. All authors have read and agreed to the published version of the manuscript.

Funding

We thank the support of the National Natural Science Foundation of China (Grant Nos. 12072261, 11872305) for our work, and J. Liu also thank the partial support of Fundamental Research Program of Shanxi Province (NO. 202103021223274) and TYUST SRIF (No. 20212074).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ragusa, M.A. On weak solutions of ultraparabolic equations. Nonlinear Anal. Theor. 2001, 47, 503–511. [Google Scholar] [CrossRef]
  2. Machado, J.T.; Mainardi, F.; Kiryakova, V. Fractional calculus: Quo vadimus? (Where are we going?). Fract. Calc. Appl. Anal. 2015, 18, 495–526. [Google Scholar] [CrossRef]
  3. Khan, I.; Ullah, H.; AlSalman, H. Fractional analysis of MHD boundary layer flow over a stretching sheet in porous medium: A new stochastic Method. J. Func. Spaces 2021, 2021, 5844741. [Google Scholar] [CrossRef]
  4. Yang, M. (Weighted pseudo) almost automorphic solutions in distribution for fractional stochastic differential equations driven by levy noise. Filomat 2021, 35, 2403–2424. [Google Scholar] [CrossRef]
  5. Sakthivel, R.; Revathib, P.; Ren, Y. Existence of solutions for nonlinear fractional stochastic differential equations. Nonlinear Anal. Theor. 2013, 81, 70–86. [Google Scholar] [CrossRef]
  6. Kamrani, M. Numerical solution of stochastic fractional differential equations. Numer. Algor. 2015, 68, 81–93. [Google Scholar] [CrossRef]
  7. Lakshmikantham, V.; Bainov, D.D.; Simeonov, P.S. Theory of Impulsive Differential Equations; World Scientific: Singapore, 1989. [Google Scholar]
  8. Perestyuk, N.A.; Plotnikov, V.A.; Samoilenko, A.M.; Skripnik, N.V. Differential Equations with Impulse Effects: Multivalued Right-Hand Sides with Discontinuities; Walter de Gruyter: Berlin, Germany, 2011. [Google Scholar]
  9. Hu, Y.; Øksendal, B. Fractional white noise calculus and applications to finance. Infin. Dimens. Anal. Quantum Probab. Relat. Top. 2003, 6, 1–32. [Google Scholar] [CrossRef]
  10. Liu, J.; Yan, L.; Cang, Y. On a jump-type stochastic fractional partial differential equation with fractional noises. Nonlinear Anal. Theor. 2012, 75, 6060–6070. [Google Scholar] [CrossRef]
  11. Li, K. Stochastic delay fractional evolution equations driven by fractional Brownian motion. Math. Methods Appl. Sci. 2015, 38, 1582–1591. [Google Scholar] [CrossRef] [Green Version]
  12. 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]
  13. Chadha, A.; Pandey, D.N. Existence results for an impulsive neutral stochastic fractional integro-differential equation with infinite delay. Nonlinear. Anal. Theor. 2015, 128, 149–175. [Google Scholar] [CrossRef]
  14. Dhayal, R.; Malik, M.; Abbas, S. Approximate controllability for a class of non-instantaneous impulsive stochastic fractional differential equation driven by fractional Brownian motion. Differ. Equ. Dyn. Syst. 2021, 29, 175–191. [Google Scholar] [CrossRef]
  15. Pedjeu, J.C.; Ladde, G.S. Stochastic fractional differential equations: Modeling, method and analysis. Chaos Solitons Fract. 2012, 45, 279–293. [Google Scholar] [CrossRef]
  16. Abouagwa, M.; Cheng, F.; Li, J. Impulsive stochastic fractional differential equations driven by fractional Brownian motion. Adv. Differ. Equ. 2020, 2020, 57. [Google Scholar] [CrossRef]
  17. Khasminskii, R.Z. A limit theorem for the solution of differential equations with random right-hand sides. Theory Probab. Appl. 1963, 11, 390–405. [Google Scholar] [CrossRef]
  18. Roberts, J.B.; Spanos, P.D. Stochastic averaging: An approximate method of solving random vibration problems. Int. J. Nonlin. Mech. 1986, 21, 111–134. [Google Scholar] [CrossRef]
  19. Zhu, W.Q. Stochastic Averaging Methods in Random Vibration. Appl. Mech. Rev. 1988, 41, 189–199. [Google Scholar] [CrossRef]
  20. Xu, Y.; Pei, B.; Guo, R. Stochastic averaging for slow-fast dynamical systems with fractional Brownian motion. Discret. Contin. Dyn. Syst. B 2015, 20, 2257–2267. [Google Scholar] [CrossRef]
  21. Xu, Y.; Pei, B.; Wu, J.L. Stochastic averaging principle for differential equations with non-Lipschitz coefficients driven by fractional Brownian motion. Stoch. Dyn. 2017, 17, 1750013. [Google Scholar] [CrossRef] [Green Version]
  22. Ma, S.; Kang, Y.M. Periodic averaging method for impulsive stochastic differential equations with Lévy noise. Appl. Math. Lett. 2019, 93, 91–97. [Google Scholar] [CrossRef]
  23. Khalaf, A.D.; Abouagwa, M.; Wang, X. Periodic averaging method for impulsive stochastic dynamical systems driven by fractional Brownian motion under non-Lipschitz condition. Adv. Differ. Equ. 2019, 526. [Google Scholar] [CrossRef] [Green Version]
  24. Cui, J.; Bi, N. Averaging principle for neutral stochastic functional differential equations with impulses and non-Lipschitz coefficients. Stat. Probabil. Lett. 2020, 163, 108775. [Google Scholar] [CrossRef]
  25. Wang, P.; Xu, Y. Periodic averaging principle for neutral stochastic delay differential equations with impulses. Complexity 2020, 2020, 6731091. [Google Scholar] [CrossRef]
  26. Liu, J.K.; Xu, W.; Guo, Q. Averaging principle for impulsive stochastic partial differential equations. Stoch. Dynam. 2021, 21, 2150014. [Google Scholar] [CrossRef]
  27. Xu, W.J.; Duan, J.Q.; Xu, W. An averaging principle for fractional stochastic differential equations with lévy noise. Chaos Interdiscip. J. Nonlinear Sci. 2020, 30, 083126. [Google Scholar] [CrossRef] [PubMed]
  28. Abouagwa, M.; Li, J. Approximation properties for solutions to Itô–Doob stochastic fractional differential equations with non-Lipschitz coefficients. Stoch. Dynam. 2019, 19, 1950029. [Google Scholar] [CrossRef]
  29. Luo, D.F.; Zhu, Q.X.; Luo, Z.G. An averaging principle for stochastic fractional differential equations with time-delays. Appl. Math. Lett. 2020, 105, 106290. [Google Scholar] [CrossRef]
  30. Shen, G.J.; Xiao, R.D.; Yin, X.W. Averaging principle and stability of hybrid stochastic fractional differential equations driven by Lévy noise. Int. J. Syst. Sci. 2020, 51, 2115–2133. [Google Scholar] [CrossRef]
  31. Guo, Z.K.; Hu, J.H.; Yuan, C.G. Averaging principle for a type of Caputo fractional stochastic differential equations. Chaos Interdiscip. J. Nonlinear Sci. 2021, 31, 053123. [Google Scholar] [CrossRef]
  32. Liu, J.K.; Xu, W. An averaging result for impulsive fractional neutral stochastic differential equations. Appl. Math. Lett. 2021, 114, 106892. [Google Scholar] [CrossRef]
  33. Russo, F.; Vallois, P. Forward, backward and symmetric stochastic integration. Probab. Theory Relat. Fields 1993, 97, 403–421. [Google Scholar] [CrossRef]
  34. Biagini, F.; Hu, Y.Z.; Øksendal, B.; Zhang, T.S. Stochastic Calculus for Fractional Brownian Motion and Applications; Springer: London, UK, 2008. [Google Scholar]
  35. Shen, L.J.; Sun, J.T. Existence and uniqueness of mild solutions for nonlinear stochastic impulsive differential equation. Abstr. Appl. Anal. 2011, 2011, 439724. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Comparison of x ε ( t ) and z ε ( t ) . (a) x 0 = 0 , ε = 0.045 , H = 0.65 , β = 0.8 , (b) x 0 = 0.5 , ε = 0.01 , H = 0.85 , β = 0.6 .
Figure 1. Comparison of x ε ( t ) and z ε ( t ) . (a) x 0 = 0 , ε = 0.045 , H = 0.65 , β = 0.8 , (b) x 0 = 0.5 , ε = 0.01 , H = 0.85 , β = 0.6 .
Fractalfract 06 00256 g001
Figure 2. Comparison of x ε ( t ) and z ε ( t ) . (a) x 0 = 0 , ε = 0.0045 , H = 0.65 , β = 0.7 , (b) x 0 = 1 , ε = 0.01 , H = 0.55 , β = 0.65 .
Figure 2. Comparison of x ε ( t ) and z ε ( t ) . (a) x 0 = 0 , ε = 0.0045 , H = 0.65 , β = 0.7 , (b) x 0 = 1 , ε = 0.01 , H = 0.55 , β = 0.65 .
Fractalfract 06 00256 g002
Figure 3. Comparison of x ε ( t ) and z ε ( t ) with t k = 4.8 k . (a) x 0 = 0 , ε = 0.005 , H = 0.65 , β = 0.7 , (b) x 0 = 1 , ε = 0.01 , H = 0.55 , β = 0.6 .
Figure 3. Comparison of x ε ( t ) and z ε ( t ) with t k = 4.8 k . (a) x 0 = 0 , ε = 0.005 , H = 0.65 , β = 0.7 , (b) x 0 = 1 , ε = 0.01 , H = 0.55 , β = 0.6 .
Fractalfract 06 00256 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, J.; Wei, W.; Xu, W. An Averaging Principle for Stochastic Fractional Differential Equations Driven by fBm Involving Impulses. Fractal Fract. 2022, 6, 256. https://doi.org/10.3390/fractalfract6050256

AMA Style

Liu J, Wei W, Xu W. An Averaging Principle for Stochastic Fractional Differential Equations Driven by fBm Involving Impulses. Fractal and Fractional. 2022; 6(5):256. https://doi.org/10.3390/fractalfract6050256

Chicago/Turabian Style

Liu, Jiankang, Wei Wei, and Wei Xu. 2022. "An Averaging Principle for Stochastic Fractional Differential Equations Driven by fBm Involving Impulses" Fractal and Fractional 6, no. 5: 256. https://doi.org/10.3390/fractalfract6050256

APA Style

Liu, J., Wei, W., & Xu, W. (2022). An Averaging Principle for Stochastic Fractional Differential Equations Driven by fBm Involving Impulses. Fractal and Fractional, 6(5), 256. https://doi.org/10.3390/fractalfract6050256

Article Metrics

Back to TopTop