Next Article in Journal
Research on the Stability of Anti-Slip Pile Support Structures for Railway Pile Slopes
Next Article in Special Issue
Some Existence and Uniqueness Results for a Class of Fractional Stochastic Differential Equations
Previous Article in Journal
Mathematical Analysis of an SIVRWS Model for Pertussis with Waning and Naturally Boosted Immunity
Previous Article in Special Issue
On the Generalized Liouville–Caputo Type Fractional Differential Equations Supplemented with Katugampola Integral Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mittag-Leffler Type Stability of Delay Generalized Proportional Caputo Fractional Differential Equations: Cases of Non-Instantaneous Impulses, Instantaneous Impulses and without Impulses

1
Department of Mathematics, Texas A&M University-Kingsville, Kingsville, TX 78363, USA
2
Faculty of Mathematics and Informatics, Plovdiv University “P. Hilendarski”, 4000 Plovdiv, Bulgaria
3
School of Mathematical and Statistical Sciences, National University of Ireland, H91 TK33 Galway, Ireland
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(11), 2290; https://doi.org/10.3390/sym14112290
Submission received: 10 October 2022 / Revised: 22 October 2022 / Accepted: 24 October 2022 / Published: 1 November 2022
(This article belongs to the Special Issue Applied Mathematics and Fractional Calculus II)

Abstract

:
In this paper, nonlinear differential equations with a generalized proportional Caputo fractional derivative and finite delay are studied in this paper. The eventual presence of impulses in the equations is considered, and the statement of initial value problems in three cases is defined: namely non-instantaneous impulses, instantaneous impulses and no impulses. The relations between these three cases are discussed. Additionally, some stability properties are investigated. We apply the Mittag–Leffler function which plays a vital role and which gives well-known bounds on the norm of the solutions. The symmetry of this function about a line and the bounds is a property that plays an important role in stability. Several sufficient conditions are presented via appropriate new comparison results and the modified Razumikhin method. The results generalize several known results in the literature.

1. Introduction

Fractional calculus in real world phenomena is very applicable because of some typical properties such as memory. Various types of kernels in fractional integrals and fractional derivatives are applied (for example, in [1,2] the fourth-order time-fractional integro-differential equation with various types of kernels are studied numerically). A very general type of kernel was studied in [3] and called a general fractional integral/derivative. These general fractional integrals and derivatives were systematically studied by Y. Luchko [4,5] in appropriate function spaces in the framework of fractional calculus. Luchko also studied some qualitative properties of solutions of various types of differential equations with general fractional derivatives (see, [5]). In this paper, we focus on stability for a particular kernel (to be described in Section 3). Stability properties for fractional differential equations were studied by many authors (see, for example, [6,7]). As mentioned in [8], the generalized energy of a system does not have to decay exponentially for the system to be stable in the sense of Lyapunov, and recently the Mittag–Leffler stability and the fractional Lyapunov direct method were introduced for various types of fractional differential equations (see, for example, [9,10,11,12]) and applied in fractional models ([13,14,15,16,17]).
Many real processes are characterized by rapid changes in their state, and they are adequately modeled by differential equations with impulses. The acting time of these changes could be short relative to the duration of the whole process and they could be modeled as instantaneous impulses (see, for example, the classical book for ordinary differential equations [18] and the cited references therein). In some processes, the duration of changes might not be negligible, i.e., they start at arbitrary fixed points and remain active on finite time intervals. These types of changes could be modeled by non-instantaneous impulses (see, the overview given in the book [19]).
Even though fractional derivatives have memory, often various types of delays are involved in the fractional differential equations to represent some dynamics of the corresponding processes. When one studies fractional differential equations with delays and any type of impulse, there are a number of technical and theoretical difficulties.
In this paper, we study nonlinear differential equations with finite delay and with a generalized proportional Caputo fractional derivative. We consider three main cases: the case when there are non-instantaneous impulses in the equation, the case when there are instantaneous impulses in the equation and the case without any impulses. In all of these cases, we set up the initial value problem and we discuss the relation between them. The appropriate Mittag–Leffler type stability is defined, and several sufficient conditions are obtained. Our study is based on the Razumikhin method and its appropriate modifications. Some of the obtained results are generalizations of results known in the literature for the case of Caputo fractional differential equations.
Our contributions in this paper include:
  • The statement of the initial value problem for nonlinear systems of generalized proportional Caputo fractional differential equations with finite delays, and we consider three cases:
    -
    With non-instantaneous impulses;
    -
    With instantaneous impulses;
    -
    Without impulses.
  • An appropriate interpretation and connection between the three cases are provided.
  • Generalized proportional Mittag–Leffler stability of the three types of systems is defined.
  • The appropriate modifications of the Razumikhin method are applied in the three cases.
  • Some extensions of the comparison principle are provided.
  • Sufficient conditions for the Mittag–Leffler-type stability are obtained.
The paper is organized as follows. In Section 2, we recall some basic definitions about generalized proportional fractional integrals and Caputo-type derivatives, and some basic results are presented. In Section 3, we discuss the statements of fractional order delay systems in our three cases, and the relationships between them is provided. In Section 4, in the three cases, the generalized proportional Mittag–Leffler stability is defined, some comparison results are proved and several sufficient conditions are obtained with the help of appropriate modifications of the Razumikhin method.

2. Preliminary Notes and Results

We will give some basic notations used in this paper.
Let u : [ 0 , b ] R n , b > 0 , b and τ ( 0 , b ) . Then, we will use the following notations u ( τ ) = u ( τ 0 ) = lim t τ u ( t ) and u ( τ + 0 ) = lim t τ u ( t ) .
Let r > 0 be a given number and consider the set E = { ϕ : [ r , 0 ] R n is continuous everywhere except at a finite number of points τ j ( r , 0 ) : ϕ ( τ j 0 ) = ϕ ( τ j ) , ϕ ( τ j + 0 ) < } with a norm | | ϕ | | 0 = sup s [ r , 0 ] | | ϕ ( s ) | | , where | | . | | is a norm in R n .
Let two sequences of points { t i } i = 1 and { s i } i = 0 be given such that 0 < s i 1 t i < s i < t i + 1 , i = 1 , 2 , , and lim k s k = . Denote t 0 = 0 .
Let J [ 0 , ) be a given interval. Consider the following classes of functions:
N P C ( J , R n ) = { u : J R n : u C [ J k = 0 ( t k , s k ] , R n ] : u ( s k ) = u ( s k 0 ) = lim t s k u ( t ) < , u ( s k + 0 ) = lim t s k u ( t ) < , k : s k J } ,
and
P C ( J , R n ) = { v : J R n : v C [ J ( [ 0 , ) / { t k } k = 1 ) , R n ] : v ( t k ) = v ( t k 0 ) = lim t t k v ( t ) < , v ( t k + 0 ) = lim t t k v ( t ) < , k : t k J } ,
We will give a brief overview of the literature on fractional integrals and derivatives with general kernels. In [4], Luchko described what was known in the literature on general fractional integrals (GFI) and general fractional derivatives (GFD) and studied GFI and GFD with the Sonine kernel. In [5], Luchko studied some analytical properties of initial-value problems for single and multi-term fractional differential equations with GFD with a Sonine kernel that possess integrable singularities of power function-type at the point zero. Luchko introduced the set of Sonine kernels S 1 and he considered GFI with a kernel κ S 1 (Definition 3.2 [5]):
( I ( κ ) f ) ( t ) = 0 t κ ( t τ ) f ( τ ) d τ , t > 0 ,
GFD of Riemann–Liouville type (Definition 3.3 [5]):
( D ( κ ) f ) ( t ) = d d t 0 t κ ( t τ ) f ( τ ) d τ , t > 0 ,
and GFD of Caputo-type (Definition 3.3 [5]):
( * D ( κ ) f ) ( t ) = ( D ( κ ) f ) ( t ) f ( 0 ) κ ( t ) , t > 0 .
In [5], the first fundamental theorem of fractional calculus for the GFD (Theorem 3.1 [5]) and the second fundamental theorem of FC for the GFD (Theorem 3.2 [5]) are proved. Additionally, an explicit form of the solution of the initial value problem (IVP) for the linear fractional differential equation with Caputo type GFD is obtained. This formula significantly depends on the kernel κ S 1 . Since the main goal of this paper is the study of fractional generalization of exponential stability, i.e., so-called Mittag–Leffler-type of stability, we will use a spacial type of the kernel κ S 1 :
κ ( t ; α , ρ ) = ρ α 1 t α Γ ( 1 α ) e ρ 1 ρ t S 1 , α ( 0 , 1 ) , ρ ( 0 , 1 ] , t 0 .
Then, the definitions of GFI and GFD given by (1)–(3) are reduced:
( I 1 α , ρ f ) ( t ) = ( I ( κ ( t , 1 α , ρ ) ) f ) ( t ) = 0 t ρ α ( t s ) α 1 Γ ( α ) e ρ 1 ρ ( t s ) f ( s ) d s , α > 0 , ρ ( 0 , 1 ] , ( R L D α , ρ f ) ( t ) = ( D ( κ ( t ; α , ρ ) ) f ) ( t ) = 1 ρ 1 α Γ ( 1 α ) d d t 0 t ( t s ) α e ρ 1 ρ ( t s ) f ( s ) d s , α ( 0 , 1 ) , ρ ( 0 , 1 ] , ( C D α , ρ f ) ( t ) = ( * D ( κ ( t ; α , ρ ) ) f ) ( t ) = ρ α 1 Γ ( 1 α ) d d t 0 t ( t s ) α e ρ 1 ρ ( t s ) f ( s ) d s f ( 0 ) ρ α 1 t α Γ ( 1 α ) e ρ 1 ρ t , f o r t > 0 , α ( 0 , 1 ) , ρ ( 0 , 1 ] .
Remark 1.
The fractional integral ( I 1 α , ρ f ) ( t ) , the fractional derivatives ( R L D α , ρ f ) ( t ) and ( C D α , ρ f ) ( t ) are called generalized proportional fractional integral, generalized proportional Rieman–Liouville fractional integral and generalized proportional Caputo fractional derive, respectively, and they are studied in [20,21].
Remark 2.
(see Remark 3.2 [20]) If α ( 0 , 1 ) and ρ ( 0 , 1 ] then the relation ( a C D α , ρ e ρ 1 ρ ( . ) ) ( t ) = 0 for t > a holds. At the same time ( a C D α , ρ K ) ( t ) 0 for K R , K 0 .
We recall some results about generalized proportional Caputo fractional derivatives and their applications in differential equations, which will be applied in the main result in the paper.
Lemma 1.
(Proposition 5.2 [20]) For ρ ( 0 , 1 ] and α ( 0 , 1 ) we have
( a C D α , ρ ( e ρ 1 ρ t ( t a ) β 1 ) ( t ) = ρ α Γ ( β ) Γ ( β α ) e ρ 1 ρ t ( t a ) β 1 α , β > 0 .
Lemma 2.
(Lemma 3.2 [22]) Let u C 1 ( [ a , b ] , R ) with a , b R , b (if b = then the interval is half open), and q ( 0 , 1 ) , ρ ( 0 , 1 ] be two reals. Then,
( a C D α , ρ u 2 ) ( t ) 2 u ( t ) ( a C D α , ρ u ) ( t ) , t ( a , b ] .
Lemma 3.
(Lemma 5 [23]) Let u C ( [ t 0 , T , R ) , T > t 0 , and there exists a point t * ( t 0 , T ] such that u ( t * ) = 0 , and u ( t ) < 0 , for t 0 t < t * . Then, if the generalized proportional Caputo fractional derivative of u exists for t = t * , then the inequality ( t 0 c D α , ρ u ) ( t ) | t = t * > 0 holds.
Lemma 4.
(Example 5.7 [20]) The scalar linear generalized proportional Caputo fractional initial value problem
( a C D α , ρ u ) ( t ) = λ u ( t ) , u ( a ) = u 0 , α ( 0 , 1 ) , ρ ( 0 , 1 ]
has a solution
u ( t ) = u 0 e ρ 1 ρ ( t a ) E α ( λ ( t a ρ ) α ) , t > a ,
where λ R , E α ( z ) = i = 0 z i Γ ( i α + 1 ) is the Mittag–Leffler function of one parameter.
Lemma 5.
Let α ( 0 , 1 ) and ρ ( 0 , 1 ] . Then
( a C D α , ρ e ρ 1 ρ ( t a ) E α λ ( t a ) ρ α = λ e ρ 1 ρ ( t a ) E α λ ( t a ) ρ α .
Proof. 
From Lemma 1 and the definition of Mittag–Leffler function with one parameter, we obtain
( a C D α , ρ ( E α ( λ ( t a ρ α ) ) ) e ρ 1 ρ ( t a ) ) = i = 0 ( a C D α , ρ ( e ρ 1 ρ ( t a ) ) ( λ ( t a ρ ) α ) i Γ ( i α + 1 ) = i = 1 λ i ρ α Γ ( α i + 1 ) e = ρ 1 ρ ( t a ) ) ( t a ) α i α ρ α i Γ ( α i + 1 α ) Γ ( i α + 1 ) = λ e ρ 1 ρ ( t a ) i = 1 λ i 1 ( t a ) α ( i 1 ) ρ α ( i 1 ) Γ ( α ( i 1 ) + 1 ) = λ e ρ 1 ρ ( t a ) E α λ ( t a ) ρ α .

3. Statement of the Problems

In this paper, we will consider three cases: non-instantaneous impulses, instantaneous impulses and without impulses, and we give the relations between them.

3.1. Non-Instantaneous Impulses

Let two sequences of points { t i } i = 1 and { s i } i = 0 be given such that 0 < s i 1 t i < s i < t i + 1 , i = 1 , 2 , , and lim k s k = . Let t 0 0 be the given fixed initial time. Without loss of generality, we will assume 0 t 0 < s 0 < t 1 .
Remark 3.
The intervals ( s k , t k + 1 ] , k = 0 , 1 , 2 , are called intervals of non-instantaneous impulses.
Let J R be a given interval. Consider the following class of functions:
N P C α , ρ ( J , R n ) = { u : J R n : u N P C ( J , R n ) : for   any k = 0 , 1 , 2 , : t k J , ( t k C D α , ρ u ) ( t ) e x i s t s f o r t ( t k , s k ] J } ,
Consider the system of non-instantaneous impulsive delay differential equations (NIDDE) with the generalized proportional Caputo fractional derivative
( t k C D α , ρ x ) ( t ) = f ( t , x t ) for t ( t k , s k ] , k = 0 , 1 , 2 , x ( t ) = Φ k ( t , x ( s k 0 ) ) for t ( s k , t k + 1 ] , k = 0 , 1 , 2 , ,
with initial condition
x ( t + t 0 ) = ϕ ( t ) for t [ r , 0 ] ,
where f : [ t 0 , s 0 ] k = 1 [ t k , s k ] × R n R n , Φ i : [ s i , t i + 1 ] × R n R n , ( i = 0 , 1 , 2 , 3 , ) , r > 0 is a given number, ϕ : [ r , 0 ] R n and x t = x ( t + s ) , s [ r , 0 ] .
Remark 4.
The functions Φ k ( t , x ) , k = 1 , 2 , , are called non-instantaneous impulsive functions.
Remark 5.
For some detailed explanations about non-instantaneous impulses in generalized proportional Caputo fractional differential equations without delays, see [24].
We will introduce the following conditions:
(A 1.1.) The function f C ( k = 0 [ t k , s k ] × R n , R n ) .
(A 1.2.) For any natural number k the functions Φ k C ( [ s k , t k ] × R n , R n ) , k = 1 , 2 , .
Remark 6.
We will assume that for any initial function ϕ E the IVP for the system of NIDDE (6) and (7) has a solution x ( t ; t 0 , ϕ ) N P C α , ρ ( [ t 0 , ) , R n ) .
We now give a brief description of the solution of IVP for NIDDE (6) and (7). The solution x ( t ; t 0 , ϕ ) of (6) and (7) is given by
x ( t ; t 0 , ϕ ) = X k ( t ) , for t ( t k , s k ] , k = 0 , 1 , 2 , , Φ k ( t , X k ( s k 0 ) ) , for t ( s k , t k + 1 ] k = 1 , 2 ,
where
-
On the interval [ t 0 r , t 0 ] , the solution satisfies the initial condition (7);
-
On the interval [ t 0 , s 0 ] , the solution coincides with X 0 ( t ) which is the solution of ( t 0 C D α , ρ x ) ( t ) = f ( t , x t ) , t ( t 0 , s 0 ] with initial condition (7);
-
On the interval ( s 0 , t 1 ] , the solution x ( t ; t 0 , ϕ ) satisfies the equation
x ( t ; t 0 , ϕ ) = Φ 0 ( t , X 0 ( s 0 0 ) ) ;
-
On the interval ( t 1 , s 1 ] , the solution coincides with X 1 ( t ) which is the solution of ( t 1 C D α , ρ x ) ( t ) = f ( t , x t ) , t ( t 1 , s 1 ] and initial condition x ( t + t 1 ) = ϕ ˜ ( t ) , t [ r , 0 ] with
ϕ ˜ ( t ) = Φ 0 ( t 1 , X 0 ( s 0 0 ) ) t = 0 x ( t t 1 ; t 0 , ϕ ) t [ r , 0 ) ;
-
On the interval ( s 1 , t 2 ] , the solution x ( t ; t 0 , ϕ ) satisfies the equation
x ( t ; t 0 , x 0 ) = Φ 1 ( t , X 1 ( s 1 0 ) ) ;
and so on.
In connection with the study of the stability properties of zero solutions, we introduce the following assumption:
(A 1.3.) The equalities f ( t , 0 ) = 0 and Φ k ( t , 0 ) 0 , k = 0 , 1 , 2 , , hold.

3.2. Instantaneous Impulses

Let the sequence of points { t i } i = 1 be given such that 0 < t i t i + 1 , i = 1 , 2 , , and lim k t k = . Let t 0 0 be the given fixed initial time. Without loss of generality we will assume 0 t 0 < t 1 .
Remark 7.
The points t k , k = 0 , 1 , 2 , are called points of impulses.
Let J R be a given interval. Consider the following class of functions
P C α , ρ ( J , R n ) = { v : J R n : v P C ( J , R n ) : for   any t k J , k = 0 , 1 , 2 , : ( t k C D α , ρ v ) ( t ) exists   for t ( t k , t k + 1 ] J } .
Consider the system of instantaneous impulsive delay differential equations (IDDE) with the generalized proportional Caputo fractional derivative
( t k C D α , ρ x ) ( t ) = f ( t , x t ) for t ( t k , t k + 1 ] , k = 0 , 1 , 2 , x ( t k + 0 ) = Ψ k ( x ( t k 0 ) ) for k = 1 , 2 , ,
with initial condition (7), where f : [ t 0 , ) × R n R n , Ψ i : R n R n , ( i = 1 , 2 , 3 , ) .
Remark 8.
The functions Ψ k ( y ) , k = 1 , 2 , , are called impulsive functions.
Remark 9.
In the case in Section 3.1 that both sequences coincide, i.e., s i = t i + 1 , i = 0 , 1 , 2 , , the system (6) is reduced to the system (9) with Φ k ( t , u ) = Ψ k ( u ) , k = 0 , 1 , 2 , , i.e., the case of non-instantaneous impulses could be considered as a generalization of the case of instantaneous impulses.
We will introduce the following conditions:
(A 2.1.) The function f C ( [ t 0 , t 1 ] k = 1 ( t k , t k + 1 ] × R n , R n ) .
(A 2.2.) The functions Φ k C ( R n , R n ) , k = 1 , 2 , .
(A 2.3.) The function f ( t , 0 ) = 0 , t t 0 and the functions Ψ k ( 0 ) = 0 , k = 1 , 2 , .
If condition (A 2.3) is satisfied, then for the zero initial function, the IVP for IDDE (7) and (9) has a zero solution.
Remark 10.
We will assume that for any initial function ϕ E the IVP for the system of IDDE (7) and (9) has a solution x ( t ; t 0 , ϕ ) P C α , ρ ( [ t 0 , ) , R n )

3.3. No Impulses

Consider the system of delay differential equations (DDE) with the generalized proportional fractional derivative
( t 0 C D α , ρ x ) ( t ) = f ( t , x t ) for t > t 0
with initial condition (7), where f : [ t 0 , ) × R n R n .
Remark 11.
The system (10) could be considered as a partial case of (9) in the case when there are no impulses, i.e., in Section 3.2  t i = t 0 , i = 1 , 2 , , i.e., the case of instantaneous impulses could be considered as a generalization of the case of without impulses.
Let J R be a given interval. Consider the following classes of functions
C α , ρ ( J , R n ) = { u : J R n : u C ( J [ a , ) , R n ) : ( a C D α , ρ u ) ( t ) exists   for t [ a , ) J } .
We will introduce the following conditions:
(A 3.1.) The function f C ( [ t 0 , ) × R n , R n ) .
(A 3.2.) The function f ( t , 0 ) = 0 , t t 0 .
Remark 12.
We will assume that for any initial function ϕ E , the IVP for the system of DDE (7) and (10) has a solution x ( t ; t 0 , ϕ ) C α , ρ ( [ t 0 , ) , R n )

4. Mittag–Leffer-Type Stability Properties

We will study the Mittag–Leffler-type stability for NIDDE (6), IDDE (9) and DDE (10) by Lyapunov functions and an appropriate modification of the Razumikhin method.

4.1. Non-Instantaneous Impulses

Definition 1.
The zero solution of the system NIDDE (6) and (7) is said to be generalized proportional Mittag–Leffler stable if there exist constants β , γ , C , λ > 0 such that the inequality
| | x ( t ; t 0 , ϕ ) | | C | | ϕ | | 0 β i = 0 k 1 e ρ 1 ρ ( s i t i ) E α ( λ ( s i t i ρ ) α ) e ρ 1 ρ ( t t k ) E α ( λ ( t t k ρ ) α ) γ , t ( t k , s k ] , k = 0 , 1 , , C | | ϕ | | 0 β i = 0 k e ρ 1 ρ ( s i t i ) E α ( λ ( s i t i ρ ) α ) γ , t ( s k , t k + 1 ] , k = 0 , 1 , 2 ,
holds, where x ( t ; t 0 , ϕ ) is a solution of the IVP for NIDDE (6) and (7) (with an arbitrary initial function ϕ E ).
Remark 13.
The definition for generalized proportional Mittag–Leffler stability for NIDDE (6) and (7) depends significantly on the type of intervals—the intervals of differential equations and the intervals of non-instantaneous impulses (see, the first and the second line, respectively, in (11)).
We will use the following class of Lyapunov-like functions (for more details, see the book [19]):
Definition 2.
Let a < b be given numbers, Ω R n , 0 Ω . Then, the function V : [ a r , b ] × Ω [ 0 , ) is from the class N Λ ( [ a r , b ] , Ω ) if:
-
V C ( [ a , b ] / { s k } × Ω , [ 0 , ) ) and it is Lipschitz with respect to the second argument;
-
For any s k ( a , b ) , x Ω , there exist finite limits V ( s k 0 , x ) = lim t s k V ( t , x ) and V ( s k + 0 , x ) = lim t s k V ( t , x ) .
We will consider the following scalar non-instantaneous impulsive differential equation (NIDE) as a comparison equation
( t k C D α , ρ u ) ( t ) = λ u ( t ) , for t ( t k , s k ] , k = 0 , 1 , 2 , , u ( t ) = Ξ k ( t , u ( s k 0 ) ) for t ( s k , t k + 1 ] , k = 0 , 1 , 2 , , u ( t 0 ) = u 0 .
According to Lemma 4, the solution of the IVP for NIDE (12) is given by
u ( t ) = u 0 e ρ 1 ρ ( t t 0 ) E α ( λ ( t t 0 ρ ) α ) , t [ t 0 , s 0 ] Ξ k ( t , u ( s k 0 ) ) , t ( s k , t k + 1 ] , k = 0 , 1 , 2 , , Ξ k 1 ( t k , u ( s k 1 0 ) ) e ρ 1 ρ ( t t k ) E α ( λ ( t t k ρ ) α ) , t ( t k , s k ] , k = 1 , 2 , .
Applying the scalar NIDE (12) as a comparison equation, we will obtain the following comparison result for NIDDE (6).
Lemma 6.
Suppose:
1. 
The function x * ( t ) = x ( t ; t 0 , ϕ ) N P C α , ρ ( [ t 0 , ) , Δ ) is a solution of the NIDDE (6) and (7), where Δ R n .
2. 
The functions Ξ k C ( [ s k , t k + 1 ] × R , R ) and Ξ k ( t , u ) u for t [ s k , t k + 1 ] , u 0 , k = 0 , 1 , 2 , .
3. 
The function V N Λ ( [ t 0 r , ) , Δ ) and
(i) 
for any t ( t k , s k ] with k = 0 , 1 , such that
V ( t , x * ( t ) ) e 1 ρ ρ ( t t k ) E α λ ( t t k ) ρ α sup s [ t r , t ] [ t k , t ] e 1 ρ ρ ( s t k ) E α λ ( s t k ) ρ α V ( s , x * ( s ) )
the inequality
t k C D α , ρ V ( t , x * ( t ) ) λ V ( t , x * ( t ) )
holds where λ > 0 is a given number.
(ii) 
For any k = 0 , 1 , the inequalities
V ( t , Φ k ( t , x * ( s k 0 ) ) ) Ξ k ( t , V ( s k 0 , x * ( s k 0 ) ) ) f o r t ( s k , t k + 1 ] .
hold.
Then, the inequality
V ( t , x * ( t ) ) M i = 0 k 1 e ρ 1 ρ ( s i t i ) E α ( λ ( s i t i ρ ) α ) e ρ 1 ρ ( t t k ) E α ( λ ( t t k ρ ) α ) , t ( t k , s k ] , k = 0 , 1 , , M i = 0 k e ρ 1 ρ ( s i t i ) E α ( λ ( s i t i ρ ) α ) , t ( s k , t k + 1 ] , k = 0 , 1 , 2 , ,
holds where M = max s [ r , 0 ] V ( t 0 + s , ϕ ( s ) ) .
Proof. 
Case 1. Let t [ t 0 , s 0 ] . Define the function m ( t ) = V ( t , x * ( t ) ) for t [ t 0 r , s 0 ] . Then, the function m ( t ) C α , ρ ( [ t 0 , s 0 ] , R + ) and the inequality m ( t 0 ) = V ( t 0 , ϕ ( 0 ) ) sup s [ r , 0 ] V ( t 0 + s , ϕ ( s ) ) = M hold. We will prove that
m ( t ) < M e ρ 1 ρ ( t t 0 ) E α λ t t 0 ρ α + ε e ρ 1 ρ ( t t 0 ) , t [ t 0 , s 0 ] ,
where ε > 0 is a small enough number. Note for t = t 0 inequality (15) holds. Assume (15) is not true on ( t 0 , s 0 ] . Therefore, there exists t * ( t 0 , s 0 ] such that
m ( t ) < M e ρ 1 ρ ( t t 0 ) E α λ t t 0 ρ α + ε e ρ 1 ρ ( t t 0 ) , t [ t 0 , t * ) , m ( t * ) = e ρ 1 ρ ( t * t 0 ) E α λ t * t 0 ρ α + ε e ρ 1 ρ ( t * t 0 ) .
Consider the function ξ ( t ) = m ( t ) M e ρ 1 ρ ( t t 0 ) E α λ t t 0 ρ α ε e ρ 1 ρ ( t t 0 ) for t [ t 0 , s 0 ] . According to Lemma 3 with u ( t ) ξ ( t ) the inequality ( t 0 c D α , ρ ξ ) ( t ) | t = t * > 0 holds. Therefore, according to Lemma 5 and Remark 2, we obtain
( t 0 c D α , ρ m ) ( t ) | t = t * > λ M e ρ 1 ρ ( t * t 0 ) E α λ t * t 0 ρ α .
Case 1.1. Let r < t * t 0 . Then, t * r > t 0 and [ t * r , t * ] ( t 0 , t * ] , i.e., [ t * r , t * ] [ t 0 , t * ] = [ t * r , t * ] . Therefore, since the function E α ( λ t ) is decreasing for t ( t 0 , t * ] , i.e., 1 E α λ ( t t 0 ) ρ α 1 E α λ ( t * t 0 ) ρ α for t [ t * r , t * ] by (16), we obtain
m ( t ) e 1 ρ ρ ( t t 0 ) E α λ t t 0 ρ α < M + ε 1 E α λ t t 0 ρ α M + ε 1 E α λ t * t 0 ρ α = m ( t * ) e 1 ρ ρ ( t * t 0 ) E α λ t * t 0 ρ α , t [ t * r , t * ] ,
i.e., inequality (13) is satisfied for t = t * .
According to condition 3(i) the inequality
t 0 c D α , ρ m ( t ) | t = t * = t 0 c D α , ρ V ( t , x * ( t ) ) | t = t * λ V ( t * , x * ( t * ) ) = λ m ( t * ) = λ M e ρ 1 ρ ( t * t 0 ) E α λ t * t 0 ρ α λ ε e ρ 1 ρ ( t * t 0 )
holds.
From inequalities (17) and (19), it follows that λ ε e ρ 1 ρ ( t * t 0 ) > 0 . The obtained contradiction proves the inequality (15) on [ t 0 , s 0 ] .
Case 1.2. Let r t * t 0 . Then, t * r t 0 and [ t * r , t * ] [ t 0 , t * ] = [ t 0 , t * ] = { t 0 } ( t 0 , t * ] . Similar to the proof in Case 1.1, we obtain the inequality
m ( t ) e 1 ρ ρ ( t t 0 ) E α λ t t 0 ρ α m ( t * ) e 1 ρ ρ ( t * t 0 ) E α λ t * t 0 ρ α , t ( t 0 , t * ] .
For t = t 0 , apply (16), E α ( 0 ) = 1 and obtain m ( t * ) e 1 ρ ρ ( t * t 0 ) E α λ t * t 0 ρ α > M m ( t 0 ) .
Therefore, inequality (13) holds for t = t * .
Thus, condition 3(i) is applicable and as in Case 1.1 we obtain a contradiction.
The contradiction proves inequality (15). From inequality (15) as ε 0 follows the validity of (14) on [ t 0 , s 0 ] .
Case 2. Let t ( s 0 , t 1 ] . Then, x * ( t ) = Φ 1 ( t , x * ( s 0 0 ) ) . From conditions 2, 3(ii) for k = 0 and Case 1, we obtain
V ( t , x * ( t ) ) = V ( t , Φ 0 ( t , x * ( s 0 0 ) ) ) Ξ 0 ( t , V ( s 0 0 , x * ( s 0 0 ) ) ) V ( s 0 0 , x * ( s 0 0 ) ) M e ρ 1 ρ ( s 0 t 0 ) E α λ s 0 t 0 ρ α , t ( s 0 , t 1 ] .
Therefore, inequality (14) holds on ( s 0 , t 1 ] .
Case 3. Let t ( t 1 , s 2 ] . Define the function
m 1 ( t ) = V ( t 1 , x * ( t 1 ) ) for t [ t 1 r , t 1 ] , V ( t , x * ( t ) ) for t ( t 1 , s 1 ] .
Then, the function m 1 ( t ) C α , ρ ( [ t 1 , s 1 ] , R + ) . Denote M 1 = V ( t 1 , x * ( t 1 ) ) . Then,
max s [ r , 0 ] m 1 ( t 1 + s ) = M 1
and according to Case 2, the inequality
M 1 < M e ρ 1 ρ ( s 0 t 0 ) E α λ s 0 t 0 ρ α
holds.
Similar to the proof of inequality (15) in Case 1, we have the validity of the inequality
m 1 ( t ) < M 1 e ρ 1 ρ ( t t 1 ) E α λ t t 1 ρ α + ε e ρ 1 ρ ( t t 1 ) . t [ t 1 , s 1 ] .
Thus,
m 1 ( t ) < M e ρ 1 ρ ( s 0 t 0 ) E α λ s 0 t 0 ρ α e ρ 1 ρ ( t t 1 ) E α λ t t 1 ρ α + ε e ρ 1 ρ ( t t 1 ) , t ( t 1 , s 1 ] .
Taking the limit in (20) as ε 0 we obtain the claim of Lemma 6 on ( t 1 , s 1 ] .
Continue this process and an induction argument proves the claim in Lemma 6. □
Remark 14.
The condition (13) is a modified Razumikhin condition applied in connection with generalized proportional fractional derivatives.
Remark 15.
The inequality (13) in condition 3(i) of Lemma 6 could be replaced by
V ( t , x * ( t ) ) sup s [ t r , t ] [ t k , t ] e 1 ρ ρ ( s t k ) E α λ ( s t k ) ρ α V ( s , x * ( s ) )
Note that if (21) holds, then inequality (13) is also satisfied.
Remark 16.
If the condition (21) is satisfied, then the classical Razumikhin condition V ( t , x * ( t ) ) sup s [ t r , t ] [ t k , t ] V ( s , x * ( s ) ) holds.
Remark 17.
The condition 3(i) is satisfied only at some particular points of t from the studied interval.
We study the generalized Mittag–Leffler stability properties of the zero solution of NIDDE (6).
Theorem 1.
Suppose:
1. 
Conditions (A 1.1)–(A 1.3) are satisfied.
2. 
There exists a function V N Λ ( [ t 0 r , ) , R n ) such that
(i) 
There exist positive constants A , B , a , b such that the inequalities A | | x | | a V ( t , x ) B | | x | | a b , t t 0 , x R n hold.
(ii) 
For any point t ( t k , s k ] with k = 0 , 1 , 2 , and any function ψ C α , ρ ( t k , [ t r , t ] , R n ) such that t k C D α , ρ ψ ( t ) = f ( t , ψ t ) and
V ( t , ψ ( t ) ) e 1 ρ ρ ( t t k ) E α D ( t t k ) ρ α sup s [ t r , t ] [ t k , t ] e 1 ρ ρ ( s t k ) E α D ( s t k ) ρ α V ( s , ψ ( s ) )
the inequality
t k C D α , ρ V ( t , ψ ( t ) ) D V ( t , ψ ( t ) )
holds where D > 0 is a given number.
(iii) 
For any k = 0 , 1 , and u R n , the inequalities
V ( t , Φ k ( t , u ) ) C | | u | | a f o r t ( s k , t k + 1 ] .
hold where C ( 0 , A ] .
Then, the zero solution of NIDDE (6) with the zero initial function is generalized proportional Mittag–Leffler stable with C = B A a , β = b , λ = D , γ = 1 a .
Proof. 
Let ϕ E be an arbitrary initial function and now let x ( t ) = x ( t ; t 0 , ϕ ) N P C α , ρ ( [ t 0 , ) , R n ) be the solution of the IVP for NIDDE (6) and (7). Let t * ( t k , s k ] with k a non-negative integer, be such that the inequality (22) holds with ψ ( t ) = x ( t ) . Note that x C α , ρ ( t k , [ t * r , t * ] , R n ) and t k C D α , ρ x ( t ) | t = t * = f ( t * , x t * ) . Then, according to condition 2(ii) of Theorem 1, the inequality (23) holds, i.e., we have
t k C D α , ρ V ( t , x ( t ) ) | t = t * D V ( t * , x ( t * ) ) ,
i.e., the condition 3(i) of Lemma 6 is satisfied with λ = D .
Let k = 0 , 1 , be an arbitrary number. Then, from conditions 2(i) and 2(iii) of Theorem 1, we obtain V ( t , Φ k ( t , x ( s k 0 ) ) ) C | | x ( s k 0 ) | | a C A V ( s k 0 , x ( s k 0 ) ) , i.e., condition 3(ii) of Lemma 6 is satisfied with Ξ k ( t , u ) = C A u u according to the choice of the constants A , C .
According to Lemma 6, the inequality
V ( t , x ( t ) ) M i = 0 k 1 e ρ 1 ρ ( s i t i ) E α ( D ( s i t i ρ ) α ) e ρ 1 ρ ( t t k ) E α ( D ( t t k ρ ) α ) , t ( t k , s k ] , k = 0 , 1 , , M i = 0 k e ρ 1 ρ ( s i t i ) E α ( D ( s i t i ρ ) α ) , t ( s k , t k + 1 ] , k = 0 , 1 , 2 , .
holds where M B | | ϕ | | 0 a b .
Thus, from condition 2(i) of Theorem 1, we obtain
| | x ( t ) | | B A a | | ϕ | | 0 b i = 0 k 1 e ρ 1 ρ ( s i t i ) E α ( D ( s i t i ρ ) α ) e ρ 1 ρ ( t t k ) E α ( D ( t t k ρ ) α ) 1 a , t ( t k , s k ] , k = 0 , 1 , , B A a | | ϕ | | 0 b i = 0 k e ρ 1 ρ ( s i t i ) E α ( D ( s i t i ρ ) α ) 1 a , t ( s k , t k + 1 ] , k = 0 , 1 , 2 , .
Thus, the zero solution of (6) is generalized Mittag-Leffler stable with C = B A a , β = b , λ = D , γ = 1 a . □
Corollary 1.
Let the conditions of Theorem 1 be satisfied where the inequality (22) is replaced by
V ( t , ψ ( t ) ) sup s [ t r , t ] [ t k , t ] e 1 ρ ρ ( s t k ) E α D ( s t k ) ρ α V ( s , ψ ( s ) )
Then, the zero solution of NIDDE (6) with the zero initial function is generalized proportional Mittag–Leffler stable.
Proof. 
If the inequality (26) is satisfied for the point t, then we obtain
V ( t , ψ ( t ) ) e 1 ρ ρ ( t t k ) E α D ( t t k ) ρ α V ( t , ψ ( t ) ) ,
i.e., inequality (22) is satisfied. □
Corollary 2.
Let the conditions of Theorem 1 be satisfied where the condition 2(ii) is replaced by 2(ii)* for any point t ( t k , s k ] with k = 0 , 1 , 2 , and any function ψ C α , ρ ( t k , [ t r , t ] , R n ) such that t k C D α , ρ ψ ( t ) = f ( t , ψ t ) and
V ( t , ψ ( t ) ) sup s [ t r , t ] [ t k , t ] e 1 ρ ρ ( s t k ) E α D ( s t k ) ρ α V ( s , ψ ( s ) )
the inequality
t k C D α , ρ V ( t , ψ ( t ) ) D sup s [ t r , t ] [ t k , t ] | | ψ ( s ) | | a b
holds where D > 0 is a given number.
Then, the zero solution of NIDDE (6) with the zero initial function is generalized proportional Mittag–Leffler stable.
Proof. 
From condition 2(iii) of Theorem 1 and inequality (27), we have that | | ψ ( s ) | | a b V ( s , ψ ( s ) ) , s [ t r , t ] [ t k , t ] , i.e.,
D sup s [ t r , t ] [ t k , t ] | | ψ ( s ) | | a b D sup s [ t r , t ] [ t k , t ] V ( s , ψ ( s ) ) = D V ( t , ψ ( t ) ) .
Thus, from inequality (28) we have inequality (23). □
Corollary 3.
Let the conditions of Theorem1 be satisfied where the inequality (23) is replaced by
t k C D α , ρ V ( t , ψ ( t ) ) 0 ,
and condition 2(i) is changed by 2(i)*. There exist positive constants A , B such that the inequalities A | | x | | V ( t , x ) B | | x | | , t t 0 , x R n hold.
Then, the zero solution of NIDDE (6) with the zero initial function is stable.
Proof. 
Inequality (29) is a partial case of (23) with D = 0 , then use E α ( 0 ) = 1 and inequality (25) and we obtain | | x ( t ) | | B A | | ϕ | | 0 for t t 0 , which proves the stability of the solution. □
Example 1.
Consider the scalar IVP for NIDDE
t k C D α , ρ x ( t ) = 2 + t t + 1 ( x ( t ) 0.5 x t ( k ) ) , f o r t ( t k , s k ] , k = 0 , 1 , 2 , , x ( t ) = 0.5 sin t x ( s k 0 ) f o r t ( s k , t k + 1 ] , k = 0 , 1 , 2 , , x ( t 0 + s ) = ϕ ( s ) , s [ r , 0 ] ,
where for any t ( t k , s k ] we denote x t ( k ) ( s ) = x ( t + s ) , s [ m a x { r , t k t } , 0 ] .
The scalar IVP for NIDDE (30) with ϕ ( s ) 0 has a zero solution.
Consider the Lyapunov function V ( t , x ) = x 2 . Then, condition 2(i) of Theorem 1 is satisfied with A = 0.25 , B = 1 , a = 2 , b = 1 . Let k be a whole number and the point t ( t k , s k ] and the function ψ C α , ρ ( t k , [ t r , t ] , R ) be such that t k C D α , ρ ψ ( t ) = 2 + t t + 1 ( ψ ( t ) 0.5 ψ t ( k ) ) and
ψ 2 ( t ) sup s [ t r , t ] [ t k , t ] e 1 ρ ρ ( s t k ) E α ( s t k ) ρ α ψ 2 ( s ) .
Then applying sup s [ t r , t ] [ t k , t ] e 1 ρ ρ ( s t k ) E α ( s t k ) ρ α ψ 2 ( s ) sup s [ t r , t ] [ t k , t ] ψ 2 ( s ) we obtain
t k C D α , ρ ψ 2 ( t ) 2 ψ ( t ) t k C D α , ρ ψ ( t ) = 2 2 + t t + 1 ( ψ 2 ( t ) 0.5 ψ ( t ) ψ t ( k ) ) 2 + t t + 1 ( 2 ψ 2 ( t ) + 0.5 ψ 2 ( t ) + 0.5 ( ψ t ( k ) ) 2 ) 2 + t t + 1 ( 2 ψ 2 ( t ) + 0.5 ψ 2 ( t ) + 0.5 sup s [ t r , t ] [ t k , t ] ψ 2 ( s ) ) 2 + t t + 1 ( 1.5 ψ 2 ( t ) + 0.5 ψ 2 ( t ) ) = 2 + t t + 1 ψ 2 ( t ) < V ( t , ψ ( t ) ) .
Let t ( s k , t k + 1 ] where k = 0 , 1 , 2 , . Then, 0.5 sin t u 2 0.25 u 2 = 0.25 | u | 2 .
Therefore, the conditions of Corollary 1 are satisfied with D = 1 , C = A = 0.25 , B = 1 , a = 2 , b = 1 . According to Corollary 1 the zero solution of the scalar NIDDE (30) is generalized proportional Mittag–Leffler stable with C = 4 = 2 , β = 1 , λ = 1 , γ = 0.5 , i.e., the inequality
| | x ( t ) | | 2 | | ϕ | | 0 i = 0 k 1 e ρ 1 ρ ( s i t i ) E α ( ( s i t i ρ ) α ) e ρ 1 ρ ( t t k ) E α ( ( t t k ρ ) α ) , t ( t k , s k ] , k = 0 , 1 , , 2 | | ϕ | | 0 i = 0 k e ρ 1 ρ ( s i t i ) E α ( ( s i t i ρ ) α ) , t ( s k , t k + 1 ] , k = 0 , 1 , 2 , .
holds.
Remark 18.
The Mittag–Leffler type stability for the Caputo fractional differential equations (with ρ = 1 ) is studied in [25].

4.2. Instantaneous Impulses

As mentioned in Remark 9, the case of non-instantaneous impulses could be considered as a generalization of the case of instantaneous impulses. That is why we can translate the results from the previous section to instantaneous impulses.
Definition 3.
The zero solution of the system IDDE (7) and (9) (with ϕ 0 ) is said to be generalized proportional Mittag–Leffler stable if there exist constants β , γ , C , λ > 0 such that the inequality
| | x ( t ; t 0 , ϕ ) | | C | | ϕ | | 0 β e ρ 1 ρ ( t t k ) E α ( λ ( t t k ρ ) α ) γ , t ( t k , t k + 1 ] , k = 0 , 1 , ,
holds, where x ( t ; t 0 , ϕ ) is a solution on the IVP for IDDE (7) and (9) with an arbitrary initial function ϕ E .
We will use some comparison results for IDDE (9) by applying piecewise continuous Lyapunov functions and we introduce a class of Lyapunov-like functions:
Definition 4.
Let a < b be given numbers, Ω R n , 0 Ω . Then, the function V : [ a r , b ] × Ω [ 0 , ) is from the class P Λ ( [ a r , b ] , Ω ) if:
-
V C ( [ a , b ] / { t k } × Ω , [ 0 , ) ) and it is Lipschitz with respect to the second argument;
-
For any t k ( a , b ) , x Ω , there exist finite limits V ( t k 0 , x ) = lim t t k V ( t , x ) and V ( t k + 0 , x ) = lim t t k V ( t , x ) .
The comparison scalar equation (IDE) is
( t k C D α , ρ u ) ( t ) = λ u ( t ) , for t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , u ( t ) = Ξ k ( u ( t k 0 ) ) for k = 1 , 2 , , u ( t 0 ) = u 0 .
According to Lemma 4, the solution of the IVP for IDE (34) is given by
u ( t ) = u 0 e ρ 1 ρ ( t t 0 ) E α ( λ ( t t 0 ρ ) α ) t [ t 0 , t 1 ] Ξ k ( u ( t k 0 ) ) e ρ 1 ρ ( t t k ) E α ( λ ( t t k ρ ) α ) t ( t k , t k + 1 ] , k = 1 , 2 , .
The auxiliary Lemma, corresponding to Lemma 6, reduces to
Lemma 7.
Suppose:
1. 
The function x * ( t ) = x ( t ; t 0 , ϕ ) P C α , ρ ( [ t 0 , ) , Δ ) is a solution of the IDDE (7) and (9) where Δ R n .
2. 
The functions Ξ k C ( R , R ) and Ξ k ( u ) u for u 0 , k = 1 , 2 , .
3. 
The function V P Λ ( [ t 0 r , ) , Δ ) and
(i) 
For any t ( t k , t k + 1 ] with k = 0 , 1 , , such that
V ( t , x * ( t ) ) e 1 ρ ρ ( t t k ) E α λ ( t t k ) ρ α sup s [ t r , t ] [ t k , t ] e 1 ρ ρ ( s t k ) E α λ ( s t k ) ρ α V ( s , x * ( s ) )
the inequality
t k C D α , ρ V ( t , x * ( t ) ) λ V ( t , x * ( t ) )
holds where λ > 0 is a given number.
(ii) 
For any k = 1 , , the inequalities
V ( t k 0 , Ψ k ( x * ( t k 0 ) ) ) Ξ k ( V ( t k 0 , x * ( t k 0 ) ) ) ,
hold.
Then, the inequality
V ( t , x * ( t ) ) max s [ r , 0 ] V ( t 0 + s , ϕ ( s ) ) e ρ 1 ρ ( t t k ) E α ( λ ( t t k ρ ) α ) , t ( t k , t k + 1 ] , k = 0 , 1 , ,
holds.
Remark 19.
The comparison scalar Equation (34) is chosen such that its explicit solution is known and condition 3(i) will be satisfied for the Lyapunov function.
Theorem 2.
Suppose:
1. 
Conditions (A 2.1)–(A 2.3) are satisfied.
2. 
There exists a function V P Λ ( [ t 0 r , ) , R n ) such that
(i) 
There exist positive constants A , B , a , b such that the inequalities A | | x | | a V ( t , x ) B | | x | | a b , t t 0 , x R n hold.
(ii) 
For any point t ( t k , t k + 1 ] with k = 0 , 1 , 2 , and any function ψ C α , ρ ( t k , [ t r , t ] , R n ) such that t k C D α , ρ ψ ( t ) = f ( t , ψ t ) and
V ( t , ψ ( t ) ) e 1 ρ ρ ( t t k ) E α D ( t t k ) ρ α sup s [ t r , t ] [ t k , t ] e 1 ρ ρ ( s t k ) E α D ( s t k ) ρ α V ( s , ψ ( s ) )
the inequality
t k C D α , ρ V ( t , ψ ( t ) ) D V ( t , ψ ( t ) ) r
holds where D > 0 is a given number.
(iii) 
For any k = 1 , 2 , and u R n the inequalities
V ( t , Ψ k ( u ) ) C | | u | | a f o r t ( t k , t k + 1 ] .
hold where C ( 0 , A ] .
Then, the zero solution of IDDE (9) with the zero initial function is generalized proportional Mittag–Leffler stable with C = B A a , β = b , λ = D , γ = 1 a .
Now we will provide an example illustrating the application of the given above sufficient conditions. To be able to compare both cases about non-instantaneous impulses and instantaneous impulses we will consider the scalar IVP for NIDDE (30) with appropriate changes.
Example 2.
Consider the scalar IVP for IDDE
t k C D α , ρ x ( t ) = 2 + t t + 1 ( x ( t ) 0.5 x t ( k ) ) f o r t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , x ( t k + 0 ) = 0.5 sin t k x ( t k 0 ) f o r k = 1 , 2 , , x ( t 0 + s ) = ϕ ( s ) , s [ r , 0 ] .
The scalar IVP for IDDE (39) with ϕ ( s ) 0 has a zero solution.
Let V ( t , x ) = x 2 . Thus, the condition 2(i) of Theorem 2 is satisfied with A = 0.25 , B = 1 , a = 2 , b = 1 .
Let k be a given natural number and t ( t k , t k + 1 ) , and the function ψ C α , ρ ( t k , [ t r , t ] , R ) be such that
t k C D α , ρ ψ ( t ) = 2 + t t + 1 ( ψ ( t ) 0.5 ψ t ( k ) )
and
ψ 2 ( t ) sup s [ t r , t ] [ t k , t ] e 1 ρ ρ ( s t k ) E α ( s t k ) ρ α ψ 2 ( s ) .
Then, we obtain t k C D α , ρ ψ 2 ( t ) < V ( t , ψ ( t ) ) (see (32)), i.e., condition 2(ii) of Theorem 2 is satisfied with D = 1 .
For any k = 1 , 2 , we obtain 0.5 sin t k u 2 0.25 u 2 = 0.25 | u | 2 , i.e., the condition 2(iii) of Theorem 2 is satisfied with C = 0.25 .
According to Theorem 2, the zero solution of the scalar IDDE (39) is a generalized proportional Mittag–Leffler stable with C = 2 , β = 1 , λ = 1 , γ = 0.5 , i.e., the inequality
| | x ( t ; t 0 , ϕ ) | | 2 | | ϕ | | 0 e ρ 1 ρ ( t t i ) E α ( ( t t k ρ ) α ) , t ( t k , t k + 1 ] , k = 0 , 1 ,
holds (compare with the special case t k + 1 = s k , k = 0 , 1 , 2 , of Example 1).

4.3. No Impulses

As mentioned in Remark 11 the case of instantaneous impulses could be considered as a generalization of the case of no impulses, i.e., the system (10) could be considered as a partial case of (9) with t i = t 0 , i = 1 , 2 , . That is why we can translate the results from the previous section to the case without impulses.
Definition 5.
The zero solution of the system DDE (10) (with ϕ 0 ) is said to be generalized proportional Mittag-Leffler stable if there exist constants β , γ , C , λ > 0 such that the inequality
| | x ( t ; t 0 , ϕ ) | | C | | ϕ | | 0 β e ρ 1 ρ ( t t 0 ) E α ( λ ( t t 0 ρ ) α ) γ , t t 0 ,
holds, where x ( t ; t 0 , ϕ ) is a solution on the IVP for DDE (7) and (10).
Remark 20.
In the case ρ = 1 , Definition 5 is the same as in [26].
We will use some comparison results for DDE (10) by applying Lyapunov functions:
Definition 6.
Let a < b be given numbers, Ω R n , 0 Ω . Then, the function V : [ a r , b ] × Ω [ 0 , ) is from the class Λ ( [ a r , b ] , Ω ) if V C ( [ a , b ] / { t k } × Ω , [ 0 , ) ) and it is Lipschitz with respect to the second argument.
The comparison scalar equation (DE) is
( t 0 C D α , ρ u ) ( t ) = λ u ( t ) , for t > t 0 , u ( t 0 ) = u 0 .
According to Lemma 4, the solution of the IVP for DE (41) is given by u ( t ) = u 0 e ρ 1 ρ ( t t 0 ) E α ( λ ( t t 0 ρ ) α ) . t t 0 .
The auxiliary Lemma, corresponding to Lemma 6 reduces to
Lemma 8.
Suppose:
1. 
The function x * ( t ) = x ( t ; t 0 , ϕ ) C α , ρ ( [ t 0 , ) , Δ ) is a solution of the DDE (7) and (10), where Δ R n .
2. 
The function V C Λ ( [ t 0 r , ) , Δ ) and for any point t > t 0 such that
V ( t , x * ( t ) ) e 1 ρ ρ ( t t 0 ) E α λ ( t t 0 ) ρ α sup s [ t r , t ] [ t 0 , t ] e 1 ρ ρ ( s t 0 ) E α λ ( s t 0 ) ρ α V ( s , x * ( s ) )
the inequality
t 0 C D α , ρ V ( t , x * ( t ) ) λ V ( t , x * ( t ) )
holds where λ > 0 is a given number.
Then, the inequality
V ( t , x * ( t ) ) max s [ r , 0 ] V ( t 0 + s , ϕ ( s ) ) e ρ 1 ρ ( t t 0 ) E α ( λ ( t t 0 ρ ) α ) , t > t 0
holds.
Theorem 3.
Suppose:
1. 
Conditions (A 3.1), (A 3.2) are satisfied.
2. 
There exists a function V Λ ( [ t 0 r , ) , R n ) such that
(i) 
There exist positive constants A , B , a , b such that C A and the inequalities A | | x | | a V ( t , x ) B | | x | | a b , t t 0 , x R n hold.
(ii) 
For any point t > t 0 and any function ψ C α , ρ ( t 0 , [ t r , t ] , R n ) such that t 0 C D α , ρ ψ ( t ) = f ( t , ψ t ) and
V ( t , ψ ( t ) ) e 1 ρ ρ ( t t 0 ) E α λ ( t t 0 ) ρ α sup s [ t r , t ] [ t 0 , t ] e 1 ρ ρ ( s t k ) E α λ ( s t 0 ) ρ α V ( s , ψ ( s ) )
the inequality
t 0 C D α , ρ V ( t , ψ ( t ) ) D V ( t , ψ ( t ) )
holds where D > 0 is a given number.
Then, the zero solution of DDE (10) with the zero initial function is generalized proportional Mittag–Leffler stable with constants C = B A a , β = b , λ = D , γ = 1 a .
Example 3.
Consider the scalar IVP for DDE
t 0 C D α , ρ x ( t ) = 2 + t t + 1 ( x ( t ) 0.5 sup s [ r , 0 ] x ( t + s ) ) , t > t 0 , x ( t 0 + s ) = ϕ ( s ) , s [ r , 0 ] .
The scalar IVP for DDE (45) with ϕ ( s ) 0 has a zero solution.
Let V ( t , x ) = x 2 . Thus, the condition 2(i) of Theorem 3 is satisfied with A = 0.25 , B = 1 , a = 2 , b = 1 .
Let t > t 0 and the function ψ C α , ρ ( t 0 , [ t r , t ] , R ) be such that t 0 C D α , ρ ψ ( t ) = 2 + t t + 1 ( ψ ( t ) 0.5 sup s [ r , 0 ] ψ ( t + s ) and ψ 2 ( t ) sup s [ t r , t ] [ t 0 , t ] e 1 ρ ρ ( s t 0 ) E α ( s t 0 ) ρ α ψ 2 ( s ) . Then, we obtain
t 0 C D α , ρ ψ 2 ( t ) < V ( t , ψ ( t ) )
(see (32)), i.e., condition 2(ii) of Theorem 3 is satisfied with D = 1 .
According to Theorem 3, the zero solution of the scalar DDE (45) is generalized proportional Mittag–Leffler stable with C = 2 , β = 1 , λ = 1 , γ = 0 , i.e., the inequality
| | x ( t ; t 0 , ϕ ) | | 2 | | ϕ | | 0 e ρ 1 ρ ( t t 0 ) E α ( ( t t 0 ρ ) α ) , t t 0 .
holds (compare with the special case of t 0 = t k , k = 1 , 2 , of Example 2).

5. Conclusions

In this paper, a system of nonlinear differential equations with finite delay and with a generalized proportional Caputo fractional derivative is studied. The basic cases are presented: the case when there are non-instantaneous impulses in the equations, the case when there are instantaneous impulses in the equations, and the case without any impulses in all equations. The appropriate initial value problem is set up in all these cases, and the relation between them is discussed. It is shown that the case of non-instantaneous impulses is a generalization of the case of instantaneous impulses, and the case of instantaneous impulses could be considered as a generalization of the case without any impulses. These statements could be applied to study various qualitative properties of the solutions. In this paper, based on the application of Lyapunov functions and an appropriate modification of the Razumikhin method, the Mittag–Leffler type stability is investigated.

Author Contributions

Conceptualization, R.P.A., S.H. and D.O.; Formal analysis, R.P.A., S.H. and D.O.; Investigation, R.P.A., S.H. and D.O.; Methodology, R.P.A., S.H. and D.O.; Supervision, R.P.A., S.H. and D.O.; Writing—original draft, R.P.A., S.H. and D.O.; Writing—review and editing, R.P.A., S.H. and D.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Bulgarian National Science Fund under Project KP-06-N32/7.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Qiu, W.; Xu, D.; Guo, J. Numerical solution of the fourth-order partial integro-differential equation with multi-term kernels by the Sinc-collocation method based on the double exponential transformation. Appl. Math. Comput. 2021, 392, 125693. [Google Scholar] [CrossRef]
  2. Xu, D.; Qiu, W.; Guo, J. A compact finite difference scheme for the fourth-order time-fractional integro-differential equation with a weakly singular kernel. Numer. Meth. Part. Differ. Equ. 2020, 36, 439–458. [Google Scholar] [CrossRef]
  3. Kochubei, A.N. General fractional calculus, evolution equations, and renewal processes. Integr. Equ. Oper. Theory 2011, 71, 583–600. [Google Scholar] [CrossRef] [Green Version]
  4. Luchko, Y. General Fractional Integrals and Derivatives with the Sonine Kernels. Mathematics 2021, 9, 594. [Google Scholar] [CrossRef]
  5. Luchko, Y. Operational calculus for the general fractional derivatives with the Sonine kernels. Fract. Calc. Appl. Anal. 2021, 24, 338–375. [Google Scholar] [CrossRef]
  6. Momani, S.; Hadid, S. Lyapunov stability solutions of fractional integrodifferential equations. Intern. J. Math. Math. Sci. 2004, 47, 2503–2507. [Google Scholar] [CrossRef]
  7. Zhang, L.G.; Li, J.M.; Chen, G.P. Extension of Lyapunov second method by fractional calculus. Pure Appl. Math. 2005, 21, 291–294. [Google Scholar]
  8. Chen, Y.Q. Ubiquitous fractional order controls? IFAC Proc. Vol. 2006, 39, 481–492. [Google Scholar] [CrossRef] [Green Version]
  9. Li, Y.; Chen, Y.Q.; Podlubny, I. Mittag-Leffler stability of fractional order nonlinear dynamic systems. Automatica 2009, 45, 2965–2969. [Google Scholar] [CrossRef]
  10. Abed-Elhameed, T.M.; Aboelenen, T. Mittag–Leffler stability, control, and synchronization for chaotic generalized fractional-order systems. Adv. Cont. Discr. Mod. 2022, 50. [Google Scholar] [CrossRef]
  11. Stamova, I. Mittag-Leffler stability of impulsive differential equations of fractional order. Quart. Appl. Math. 2015, 73, 525–535. [Google Scholar] [CrossRef]
  12. Yang, X.; Li, C.; Huang, T.; Song, Q. Mittag–Leffler stability analysis of nonlinear fractional-order systems with impulses. Appl. Math. Comput. 2017, 293, 416–422. [Google Scholar] [CrossRef]
  13. Gao, Y. Mittag-Leffler stability for a new coupled system of fractional-order differential equations on network. Adv. Differ. Equ. 2018, 121. [Google Scholar] [CrossRef]
  14. Sene, N. Mittag-Leffler input stability of fractional differential equations and its applications. Discr. Cont. Dynam. Syst.-S 2020, 13, 867–880. [Google Scholar] [CrossRef] [Green Version]
  15. Tatar, N. Mittag—Leffler stability for a fractional Euler—Bernoulli problem. Chaos Solitons Fractals 2021, 149, 1110777. [Google Scholar] [CrossRef]
  16. Moharramnia, A.; Eghbali, N.; Rassias, J.M. Mittag-Leffler-Hyers-Ulam stability of Prabhakar fractional integral equation. Int. J. Nonlinear Anal. Appl. 2021, 12, 25–33. [Google Scholar]
  17. Wu, A.; Liu, L.; Huang, T.; Zeng, Z. Mittag-Leffler stability of fractional-order neural networks in the presence of generalized piecewise constant arguments. Neural Netw. 2017, 85, 118–127. [Google Scholar] [CrossRef]
  18. Lakshmikantham, V.; Bainov, D.D.; Simeonov, P.S. Theory of Impulsive Differential Equations; World Scientific: Singapore, 1989. [Google Scholar]
  19. Agarwal, R.; Hristova, S.; O’Regan, D. Non-Instantaneous Impulses in Differential Equations; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  20. Jarad, F.; Abdeljawad, T.; Alzabut, J. Generalized fractional derivatives generated by a class of local proportional derivatives. Eur. Phys. J. Spec. Top. 2017, 226, 3457–3471. [Google Scholar] [CrossRef]
  21. Jarad, F.; Abdeljawad, T. Generalized fractional derivatives and Laplace transform. Discret. Contin. Dyn. Syst. Ser. S 2020, 13, 709–722. [Google Scholar] [CrossRef] [Green Version]
  22. Bohner, M.; Hristova, S. Stability for generalized Caputo proportional fractional delay integro-differential equations. Bound Value Probl. 2022, 14, 14. [Google Scholar] [CrossRef]
  23. Almeida, R.; Agarwal, R.P.; Hristova, S.; O’Regan, D. Quadratic Lyapunov functions for stability of generalized proportional fractional differential equations with applications to neural networks. Axioms 2021, 10, 322. [Google Scholar] [CrossRef]
  24. Agarwal, R.; Hristova, S.; O’Regan, D. Generalized Proportional Caputo Fractional Differential Equations with Noninstantaneous Impulses: Concepts, Integral Representations, and Ulam-Type Stability. Mathematics 2022, 10, 2315. [Google Scholar] [CrossRef]
  25. Agarwal, R.; Hristova, S.; O’Regan, D. Mittag-Leffler stability for non-instantaneous impulsive Caputo fractional differential equations with delays. Math. Slovaca 2019, 69, 583–598. [Google Scholar] [CrossRef]
  26. Sadati, S.J.; Baleanu, D.; Ranjbar, A.; Ghaderi, R.; Abdeljawad, T. Mittag–Leffler stability theorem for fractional nonlinear systems with delay. Abstract Appli. Anal. 2010, 2010, 108651. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Agarwal, R.P.; Hristova, S.; O’Regan, D. Mittag-Leffler Type Stability of Delay Generalized Proportional Caputo Fractional Differential Equations: Cases of Non-Instantaneous Impulses, Instantaneous Impulses and without Impulses. Symmetry 2022, 14, 2290. https://doi.org/10.3390/sym14112290

AMA Style

Agarwal RP, Hristova S, O’Regan D. Mittag-Leffler Type Stability of Delay Generalized Proportional Caputo Fractional Differential Equations: Cases of Non-Instantaneous Impulses, Instantaneous Impulses and without Impulses. Symmetry. 2022; 14(11):2290. https://doi.org/10.3390/sym14112290

Chicago/Turabian Style

Agarwal, Ravi P., Snezhana Hristova, and Donal O’Regan. 2022. "Mittag-Leffler Type Stability of Delay Generalized Proportional Caputo Fractional Differential Equations: Cases of Non-Instantaneous Impulses, Instantaneous Impulses and without Impulses" Symmetry 14, no. 11: 2290. https://doi.org/10.3390/sym14112290

APA Style

Agarwal, R. P., Hristova, S., & O’Regan, D. (2022). Mittag-Leffler Type Stability of Delay Generalized Proportional Caputo Fractional Differential Equations: Cases of Non-Instantaneous Impulses, Instantaneous Impulses and without Impulses. Symmetry, 14(11), 2290. https://doi.org/10.3390/sym14112290

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