Next Article in Journal
Acknowledgement to Reviewers of Axioms in 2018
Next Article in Special Issue
Note on Limit-Periodic Solutions of the Difference Equation xt + 1 − [h(xt) + λ]xt = rt, λ > 1
Previous Article in Journal
Separability of Topological Groups: A Survey with Open Problems
Previous Article in Special Issue
A New Efficient Method for the Numerical Solution of Linear Time-Dependent Partial Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lipschitz Stability for Non-Instantaneous Impulsive Caputo Fractional Differential Equations with State Dependent Delays

1
Department of Mathematics, Texas A&M University-Kingsville, Kingsville, TX 78363, USA
2
Distinguished University Professor of Mathematics, Florida Institute of Technology, Melbourne, FL 32901, USA
3
Department of Applied Mathematics and Modeling, University of Plovdiv “Paisii Hilendarski”, 4000 Plovdiv, Bulgaria
4
School of Mathematics, Statistics and Applied Mathematics, National University of Ireland, H91 CF50 Galway, Ireland
*
Author to whom correspondence should be addressed.
Submission received: 21 November 2018 / Revised: 24 December 2018 / Accepted: 25 December 2018 / Published: 29 December 2018
(This article belongs to the Special Issue New Trends in Differential and Difference Equations and Applications)

Abstract

:
In this paper, we study Lipschitz stability of Caputo fractional differential equations with non-instantaneous impulses and state dependent delays. The study is based on Lyapunov functions and the Razumikhin technique. Our equations in particular include constant delays, time variable delay, distributed delay, etc. We consider the case of impulses that start abruptly at some points and their actions continue on given finite intervals. The study of Lipschitz stability by Lyapunov functions requires appropriate derivatives among fractional differential equations. A brief overview of different types of derivative known in the literature is given. Some sufficient conditions for uniform Lipschitz stability and uniform global Lipschitz stability are obtained by an application of several types of derivatives of Lyapunov functions. Examples are given to illustrate the results.

1. Introduction

Many papers in the literature study stability of solutions of differential equations via Lyapunov functions. One type of stability, useful in real world problems, is the so-called Lipschitz stability and Dannan and Elaydi [1] introduced the notion of Lipschitz stability for ordinary differential equations. As noted in [1], this type of stability is important only for nonlinear problems since it coincides with uniform stability in linear systems. Based on theoretical results for Lipschitz stability in [1], the dynamic behavior of a spacecraft when a single magnetic torque-rod is used for achieving a pure spin condition is studied in [2]. Recently, stability properties of delay fractional differential equations without any type of impulse are considered and we refer the reader to [3] and the references therein.
In this paper, we study the Lipschitz stability for a nonlinear system of non-instantaneous impulsive fractional differential equations with state dependent delay (NIFrDDE). The impulses start abruptly at some points and their actions continue on given finite intervals. Non-instantaneous impulsive differential equations were introduced by Hernandez and O’Regan in 2013 (see, for example, [4]). The systematic description of solutions of both ordinary and Caputo fractional differential equations with non-instantaneous impulse and without delays is given in the monograph [5]. In addition, some results for non-instantaneous fractional equations without any type of delay are presented in [6,7,8]. In [9], Caputo fractional differential equations with time varying delays is considered (we note that the model had no impulses). However, in this paper, for the first time, we consider together
  • Lipschitz stability;
  • state dependent delays (note a special case is time varying delays); and
  • models with non-instantaneous impulses.
There are two different approaches in the literature for the interpretation of the solution of fractional differential equations with impulses (for more details, see [6] and Chapter 2 of the book [5]). In the first interpretation, the lower limit of the fractional derivative is one and the same on the whole interval of study and at each point of jump we consider a boundary value problem defined by the impulsive function. In the second interpretation, the lower limit of the fractional derivative changes at each time of jump with the idea of considering an initial value problem at each jump point.
In this paper, we use the second approach to study Lipschitz stability properties of nonlinear non-instantaneous impulsive delay differential equations. The delays are bounded and depend on both the time and the state. Note several stability properties are studied in the literature for Caputo fractional differential equations (for example, see [10] (without delays), [3] (with delays and no impulses), and [11] (with multiple discrete delays without impulses)). Our study is based on Lyapunov functions and the Razumikhin technique. A brief overview in the literature of different types of derivatives of Lyapunov functions among the studied fractional differential equation is given. Several sufficient conditions for uniform Lipschitz stability and global uniform Lipschitz stability are obtained by an application of these derivatives. Some examples illustrating the results are given.

2. Notes on Fractional Calculus

We give the main definition of fractional derivatives used in the literature (see, for example, [12,13,14]). We give these definitions for scalar functions. Throughout the paper, we assume q ( 0 , 1 ) .
-
Riemann–Liouville (RL) fractional derivative:
D t q t 0     R L m ( t ) = 1 Γ 1 q d d t t 0 t t s q m ( s ) d s , t t 0
where Γ ( . ) denotes the Gamma function.
-
Caputo fractional derivative
D t q t 0 C m ( t ) = 1 Γ 1 q t 0 t t s q m ( s ) d s , t t 0 .
Note that for a constant m the equality D t q t 0 C m = 0 holds. However, for any given t * , we denote D t q t 0 C m ( t * ) = D t q t 0 C m ( t ) | t = t * .
-
The Grünwald–Letnikov fractional derivative is given by
D t q t 0     G L m ( t ) = lim h 0 1 h q r = 0 [ t t 0 h ] 1 r q C r m t r h , t t 0
and the Grünwald–Letnikov fractional Dini derivative by
D + q t 0     G L m ( t ) = lim   sup h 0 + 1 h q r = 0 [ t t 0 h ] ( 1 ) r q C r m ( t r h ) , t t 0 ,
where C r q = q ( q 1 ) ( q r + 1 ) r ! and [ t t 0 h ] denotes the integer part of the fraction t t 0 h .
From the relation between the Caputo fractional derivative and the Grünwald–Letnikov fractional derivative using Equation (1), we define the Caputo fractional Dini derivative of a function as
D + q t 0 C m ( t ) = t 0 G L D + q [ m ( t ) m ( t 0 ) ] ,
i.e.,
D + q t 0     G L m ( t ) = lim sup h 0 + 1 h q [ m ( t ) m ( t 0 ) r = 1 [ t t 0 h ] ( 1 ) r + 1 q r ( m ( t r h ) m ( t 0 ) ) ] .
The fractional derivatives for scalar functions could be easily generalized to the vector case by taking fractional derivatives with the same fractional order for all components.

3. Statement of the Problem and Basic Definitions

Let the positive constant r be given and the points { t i } 1 , { s i } 1 be such that 0 < s i < t i < s i + 1 , i = 1 , 2 , . Let t 0 0 be the given initial time. Without loss of generality, we can assume t 0 [ 0 , s 1 ) .
Consider the space P C 0 of all functions y : [ r , 0 ] R n , which are piecewise continuous endowed with the norm | | y | | P C 0 = sup t [ r , 0 ] { | | y ( t ) | | : y P C 0 } where | | . | | is a norm in R n .
The intervals ( t i , s i + 1 ) , i = 0 , 1 , 2 , are the intervals on which the fractional differential equations are given and on the intervals ( s i , t i ) , i = 1 , 2 , the impulsive conditions are given.
The Caputo fractional derivative has a memory and it depends significantly on its lower derivative. This property as well as the meaning of impulses in the differential equation lead to two basic approaches to Caputo fractional differential equations with non-instantaneous impulses:
-
Unchangeable lower limit of the Caputo fractional derivative: the lower limit of the fractional derivative is equal to the initial time t 0 on the whole interval of consideration.
-
Changeable lower limit of the Caputo fractional derivative: the lower limit of the fractional derivative is equal to the left end t i on the interval ( t i , s i + 1 ) , i = 0 , 1 , 2 , without impulses.
In this paper, we study the case of changeable lower limit of the Caputo fractional derivative.
Consider the initial value problem (IVP) for a nonlinear system of non-instantaneous impulsive fractional differential equations with state dependent delay (NIFrDDE) with q ( 0 , 1 ) :
D t q t i C x ( t ) = f ( t , x ( t ) , x ρ ( t , x t ) ) for t ( t i , s i + 1 ] , i = 0 , 1 , 2 , , x ( t ) = ϕ i ( t , x ( s i ) ) , t ( s i , t i ] , i = 1 , 2 , , x ( t + t 0 ) = φ ( t ) for t [ r , 0 ] ,
where x R n , D t q t i C x ( t ) denotes the Caputo fractional derivative with lower limit t i for the state x ( t ) , the functions f : [ 0 , s 1 ] i = 1 [ t i , s i + 1 ] × R n × P C 0 R n ; ρ : [ 0 , s 1 ] i = 1 [ t i , s i + 1 ] × P C 0 R , φ P C 0 ; ϕ i : [ s i , t i ] × R n R n , i = 1 , 2 , . Here, x t ( s ) = x ( t + s ) , s [ r , 0 ] , i.e., represents the history of the state from time t r up to the present time t. Note that for any t 0 we let x ρ ( t , x t ) = x ( ρ ( t , x ( t + s ) ) ) , s [ r , 0 ] , i.e., the function ρ determines the state-dependent delay. Note, the integer order differential equations with non-instantaneous impulses and state dependent delay are studied in [15].
Let P C [ t 0 , ) be the space of all functions y : [ t 0 r , ) R n which are piecewise continuous on [ t 0 r , ) with points of discontinuity s i , i = 1 , 2 , , the limits y ( s i 0 ) = lim t s i , t < s i y ( t ) = y ( s i ) and y ( s i + ) = lim t s i , t > s i y ( t ) exist, for any t ( t i , s i ] the Caputo fractional derivative D t q t i C y ( t ) , i = 0 , 1 , , exists and it is endowed with the norm | | y | | P C = sup t [ t 0 r , ) { | | y ( t ) | | : y P C [ t 0 , ) } where | | . | | is a norm in R n .
Define the set P C q [ t 0 , ) = { y C ( i = 0 ( t i , s i ] , R n ) such   that   for   any t ( t i , s i ] : t i t ( t s ) q 1 y ( s ) d s < , i = 1 , 2 , } .
We introduce the assumptions:
A1. 
The function f C ( [ 0 , s 1 ] i = 1 [ t i , s i + 1 ] × R n × P C 0 , R n ) is such that for any y R n , u P C 0 the inclusion f ( . , y , u ) P C q [ 0 , ) holds.
A2. 
The function ρ C ( [ 0 , s 1 ] i = 1 [ t i , s i + 1 ] × P C 0 , [ r , ) ) and for any ( t , y ) i = 0 k [ t i , s i + 1 ] × P C 0 the inequalities t r ρ ( t , y ) t holds.
A3. 
The functions ϕ i C ( [ s i , t i ] × R n , R n ) , i = 1 , 2 , .
A4. 
The function φ P C 0 .
A5. 
The function f ( t , 0 ) = 0 for t [ 0 , s 1 ] i = 1 [ t i , s i + 1 ] and ϕ i ( t , 0 ) = 0 for t [ s i , t i ] , i = 1 , 2 , .
Remark 1.
Assumption A5 guarantees the existence of the zero solution of IVP for NIFrDDE (Equation (1)) with the zero initial function φ 0 .
Remark 2.
Assumption A2 guarantees the delay of the argument in Equation (1).
Definition 1.
Let the conditions A1–A4 be satisfied. The function x P C [ t 0 , ) is a solution of the IVP in Equation (1) iff it satisfies the following integral-algebraic equation
x ( t ) = φ ( t ) , t [ r , 0 ] , φ ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 f ( s , x ( s ) , x ρ ( s , x s ) ) d s , t ( 0 , s 1 ] , ϕ i ( t , x ( s i ) ) , t ( s i , t i ] , i = 1 , 2 , , ϕ i ( t i , x ( s i ) ) + 1 Γ ( q ) t i t ( t s ) q 1 f ( s , x ( s ) , x ρ ( s , x s ) ) d s , t ( t i , s i + 1 ] , i = 1 , 2 , .
Definition 2.
The functions f , ρ are defined only on the intervals without impulses on which the differential equation is given.
We generalize Lipschitz stability ([1]) for ordinary differential equations to systems of Caputo fractional non-instantaneous impulsive differential equations with state dependent delay.
Definition 3.
The zero solution of NIFrDDE (Equation (1)) is said to be:
-
Uniformly Lipschitz stable if there exists M 1 and δ > 0 such that, for any for any initial time t 0 [ 0 , s 1 ] k = 1 [ t k , s k ] and any initial function φ P C 0 , the inequality | | φ | | P C 0 < δ implies | | x ( t ; t 0 , φ ) | | M | | φ | | P C 0 for t t 0 where x ( t ; t 0 , φ ) is a solution of Equation (1).
-
Globally uniformly Lipschitz stable if there exists M 1 such that, for any initial time t 0 [ 0 , s 1 ] k = 1 [ t k , s k ] and any initial function φ P C 0 , the inequality | | φ | | P C 0 < implies | | x ( t ; t 0 , φ ) | | M | | φ | | P C 0 for t t 0 .
Let J R + , 0 J , ρ > 0 . Consider the following sets:
M ( J ) = { a C [ J , R + ] : a ( 0 ) = 0 , a ( r ) is strictly increasing in J , and a 1 ( α r ) r q a ( α ) for some function q a : q a ( α ) 1 , if α 1 } , K ( J ) = { a C [ J , R + ] : a ( 0 ) = 0 , a ( r ) is strictly increasing in J , and a ( r ) K a r for some constant K a > 0 } , S ρ = { x R n : | | x | | ρ } .
Remark 3.
The function a ( u ) = K 1 u , K 1 > 0 is from the class K ( R + ) with K a = K 1 . The function a ( u ) = K 2 u 2 , K 2 ( 0 , 1 ] is from the class M ( [ 1 , ) ) with q ( u ) = u K 2 1 for u 1 .

4. Lyapunov Functions and Their Derivatives among Nonlinear Non-Instantaneous Caputo Delay Fractional Differential Equations

One approach to study Lipschitz stability of solutions of Equation (1) is based on using Lyapunov-like functions. The first step is to define a Lyapunov function. The second step is to define its derivative among the fractional equation.
We use the class Λ of Lyapunov-like functions, defined and used for impulsive differential equations in [16].
Definition 4.
Let J R + be a given interval, and Δ R n be a given set. We say that the function V ( t , x ) : J × Δ R + , belongs to the class Λ ( J , Δ ) if
-
The function V ( t , x ) is continuous on J / { s k J } × Δ and it is locally Lipschitz with respect to its second argument.
-
For each s k J and x Δ , there exist finite limits
V ( s k , x ) = 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 ) .
In connection with the Caputo fractional derivative, it is necessary to define in an appropriate way the derivative of Lyapunov functions among the studied equation. We give a brief overview of the derivatives of Lyapunov functions among solutions of fractional differential equations known and used in the literature. There are mainly three types of derivatives of Lyapunov functions from the class Λ ( J , Δ ) used in the literature to study stability properties of solutions of Caputo fractional differential in Equation (1):
-
First type: the Caputo fractional derivative of the function V ( t , x ( t ) ) Λ ( [ a , b ) , Δ ) defined by
D q t k c   V ( t , x ( t ) ) = 1 Γ 1 q t k t t s q d d s V ( s , x ( s ) ) d s , t [ t k , s k + 1 )
where x ( t ) is a solution of Equation (1).
-
Second type: Dini fractional derivative of the Lyapunov function V Λ ( [ t 0 , ) , R n ) among Equation (1): Let ϕ P C 0 and t ( t k , s k + 1 ) for a non-negative integer k. Then,
D ( 1 ) + V ( t , ϕ ( 0 ) , t k , ϕ ) = lim   sup h 0 1 h q [ V ( t , ϕ ( 0 ) ) r = 1 [ t t k h ] ( 1 ) r + 1 q C r V ( t r h , ϕ ( 0 ) h q f ( t , ϕ ( 0 ) , ϕ ( ρ ( t , ϕ 0 ) t ) ) ) ]
where ϕ 0 ( s ) = ψ ( s ) and ϕ ( ρ ( t , ϕ 0 ) t ) = ϕ ( ρ ( t , ϕ ( s ) ) t ) for any s [ r , 0 ] . We note that, because of Assumption A2, the inequality t r < ρ ( t , ϕ ( s ) ) < t holds ( r < ρ ( t , ϕ ( s ) ) t < 0 ), i.e ϕ ( ρ ( t , ϕ ( s ) ) t ) is well defined.
The derivative of Equation (4) keeps the concept of fractional derivatives because it has a memory.
-
Third type: Caputo fractional Dini derivative of a Lyapunov function V Λ ( [ t 0 , ) , R n ) among Equation (1): Let the initial function φ P C 0 be given and the function ϕ P C 0 and t ( t k , s k + 1 ) for a non-negative integer k. Then,
D + q ( 1 ) c     V ( t , ϕ ; t k , φ ( 0 ) ) = lim   sup h 0 + 1 h q { V ( t , ϕ ( 0 ) ) V ( t k , φ ( 0 ) ) r = 1 [ t t k h ] ( 1 ) r + 1 q C r V ( t r h , ϕ ( 0 ) h q f ( t , ϕ ( 0 ) , ϕ ( ρ ( t , ϕ 0 ) t ) ) ) V ( t k , φ ( 0 ) ) } ,
or its equivalence
D + q ( 1 ) c     V ( t , ϕ ; t k , φ ( 0 ) ) = lim   sup h 0 + 1 h q { V ( t , ϕ ( 0 ) ) + r = 1 [ t t k h ] ( 1 ) r q C r V ( t r h , ϕ ( 0 ) h q f ( t , ϕ ( 0 ) , ϕ ( ρ ( t , ϕ 0 ) t ) ) ) } V ( t k , φ ( 0 ) ) ( t t k ) q Γ ( 1 q ) .
The derivative D + q ( 1 ) c     V ( t , ϕ ; t k , φ ( 0 ) ) given by Equation (6) depends significantly on both the fractional order q and the initial data ( t k , φ ) of IVP for FrDDE (Equation (1)) and it makes this type of derivative close to the idea of the Caputo fractional derivative of a function.
Remark 4.
For any initial data ( t k , φ ) R + × P C 0 of the IVP for NIFrDDE (Equation (1)) and any function ϕ P C 0 and any point t ( t k , s i + 1 ) for a non-negative integer k the relations
D + q ( 1 ) c     V ( t , ϕ ; t k , φ ( 0 ) ) = D ( 1 ) + V ( t , ϕ ( 0 ) , t k , ϕ ) D q t k    R L V ( t k , φ ( 0 ) ) ,
D + q ( 1 ) c     V ( t , ϕ ; t k , φ ( 0 ) ) = D ( 1 ) + V ( t , ϕ ( 0 ) , t k ϕ ) , i f V ( t k , φ ( 0 ) ) = 0
D + q ( 1 ) c     V ( t , ϕ ; t k , φ ( 0 ) ) < D ( 1 ) + V ( t , ϕ ( 0 ) , t k , ϕ ) , i f V ( t k , φ ( 0 ) ) > 0 .
are satisfied.
Remark 5.
A derivative of V ( t , x ) Λ ( J , Δ ) among a system of Caputo fractional differential equations without delays was introduced by V. Lakshmikantham et al. [17] in 2009. Later, it was generalized for fractional equations with delays ([18,19,20]):
D ( 1 ) + V ( t , ϕ ( 0 ) , ϕ ) = lim   sup h 0 1 h q V ( t , ϕ ( 0 ) ) V ( t h , ϕ ( 0 ) h q f ( t , ϕ ) ) , t t 0
where ϕ C ( [ τ , 0 ] , Δ ) .
This definition is a direct generalization of the well known Dini derivative among differential equations with ordinary derivatives. However, for equations with fractional derivatives, it seems strange. It does not depend on the order q of the fractional derivative nor on the initial time t 0 . The operator defined by Equation (9) has no memory, which is typical for the fractional derivative.
The derivative D ( 1 ) + V ( t , ϕ ( 0 ) , ϕ ) defined by Equation (9) is applied in [18] to study stability of fractional delay differential equations where in the proof of the main comparison result (Theorem 4.3 [18]) the derivative D ( 1 ) + V ( t , ϕ ( 0 ) , ϕ ) is incorrectly substituted by the Caputo fractional derivative (see Equations (20) and (30) in [18]). A similar situation occurs with the application of the derivative of Equation (9) in [20] for studying stability of impulsive fractional differential equations.
In the next example to simplify the calculations and to emphasize the derivatives and their properties, we consider the scalar case, i.e., n = 1 .
Example 1. 
(Lyapunov function depending directly on the time variable). Let V ( t , x ) = m ( t ) x 2 where m C 1 ( R + , R + ) .
Case 1. Caputo fractional derivative. Let x be a solution of NIFrDDE (Equation (1)). Then, the fractional derivative
D q t 0 c   V ( t , x ( t ) ) = D q t 0 c   m ( t ) x 2 ( t ) = 1 Γ ( 1 q ) t 0 t m ( s ) x 2 ( s ) + 2 m ( s ) x ( s ) x ( s ) ( t s ) q d s
is difficult to obtain in the general case for any solution of Equation (1). In addition, the solution x ( t ) might not be differentiable on the intervals of impulses.
Case 2. Dini fractional derivative. Let ϕ P C 0 and t ( t k , s k + 1 ) for a non-negative integer k. Then, applying Equation (4), we obtain
D ( 1 ) + V ( t , ϕ ( 0 ) , t k , ϕ ) = lim sup h 0 1 h q [ m ( t ) ( ϕ ( 0 ) ) 2 r = 1 [ t t k h ] ( 1 ) r + 1 q C r m ( t r h ) ( ϕ ( 0 ) h q f ( t , ϕ ( 0 ) , ϕ ( ρ ( t , ϕ 0 ) t ) ) ) ) 2 ] = lim sup h 0 1 h q [ m ( t ) ( ϕ ( 0 ) ) 2 ( ϕ ( 0 ) h q f ( t , ϕ ( 0 ) , ϕ ( ρ ( t , ϕ 0 ) t ) ) ) ) 2 + ( ϕ ( 0 ) h q f ( t , ϕ ( 0 ) , ϕ 0 ) ) 2 r = 0 [ t t k h ] ( 1 ) r q C r m ( t r h ) ] = ϕ ( 0 ) m ( t ) f ( t , ϕ ( 0 ) , ϕ ( ρ ( t , ϕ 0 ) t ) ) ) + ( ϕ ( 0 ) ) 2 D q t k   R L m ( t ) .
Case 2. Caputo fractional Dini derivative. Let φ , ϕ P C 0 and t ( t k , s k + 1 ) for a non-negative integer k. Then, we use Equation (6) and obtain
D + q ( 1 ) c     V ( t , ϕ ; t k , φ ( 0 ) ) = lim sup h 0 + 1 h q { ϕ ( 0 ) 2 m ( t ) r = 1 [ t t k h ] ( 1 ) r + 1 q C r m ( t r h ) ( ϕ ( 0 ) h q f ( t , ϕ ( 0 ) , ϕ ( ρ ( t , ϕ 0 ) t ) ) ) ) 2 } ( φ ( 0 ) ) 2 m ( t k ) ( t t k ) q Γ ( 1 q ) = 2 ϕ ( 0 ) m ( t ) f ( t , ϕ ( 0 ) , ϕ ( ρ ( t , ϕ 0 ) t ) ) ) + ( ϕ ( 0 ) ) 2 D q t k   R L t k R L D q m ( t ) ( φ ( 0 ) ) 2 m ( t k ) ( t t k ) q Γ ( 1 q ) = D ( 1 ) + V ( t , ϕ ( 0 ) , t k , ϕ ) V ( t k , φ ( 0 ) ) ( t t k ) q Γ ( 1 q ) .
 □

5. Comparison Results

Lemma 1.
[17]. Let v C ( [ a , b ] , R ) be such that ( t a ) 1 q v C ( [ a , b ] , R ) and there exists a point τ ( a , b ] : v ( τ ) = 0 and v ( t ) 0 for t [ a , τ ] . Then, D τ q a C v ( τ ) 0 .
We use the following comparison scalar fractional differential equation with non-instantaneous impulses:
D t q t i C u ( t ) = g ( t , u ( t ) ) f o r t ( t i , s i + 1 ] , i = 0 , 1 , 2 , , u ( t ) = ψ i ( t , u ( s i 0 ) ) , t ( s i , t i ] , i = 1 , 2 , , u ( t 0 ) = u 0 ,
where u , u 0 R , g : [ 0 , s 1 ] k = 1 [ t k , s k ] × R R , ψ k : [ s k , t k + 1 ] × R R ( k = 1 , 2 , 3 , ) .
We obtain some comparison results. Note some comparison results for fractional time delay differential equations are obtained in [18] by applying the derivative defined by Equation (9) and substituting it incorrectly as a Caputo fractional derivative (see Remark 5).
We introduce the following conditions:
A6. 
The function g ( t , u ) C ( [ 0 , s 1 ] k = 1 [ t k , s k + 1 ] × R + , R ) is strictly decreasing with respect to its second argument, and for any k = 1 , 2 , the functions ψ k : [ s k , t k ] × R + R + are nondecreasing with respect to their second argument.
A7. 
The function g ( t , 0 ) = 0 for t [ 0 , s 1 ] k = 1 [ t k , s k + 1 ] and for any k = 1 , 2 , the function ψ k ( t , 0 ) = 0 for t [ s k , t k ] .
A8. 
For all k = 1 , 2 , , the functions ψ k satisfies ψ k ( t , u ) u , t [ s k , t k ] , u R .
In our main results, we use the Lipschitz stability of the zero solution of the scalar comparison non-instantaneous impulsive fractional differential in Equation (10).
Example 2.
Let t k = 2 k , k = 0 , 1 , and s k = 2 k 1 , k = 1 , 2 , . Consider the scalar non-instantaneous impulsive fractional differential equation
D t 0.25 t i C u ( t ) = u ( t ) f o r t ( t i , s i + 1 ] , i = 0 , 1 , 2 , , u ( t ) = ψ k ( t , u ( s k 0 ) ) , t ( s i , t i ] , i = 1 , 2 , , u ( 0 ) = u 0 ,
where u , u 0 R .
Case 1. Suppose for all natural numbers k = 1 , 2 , the equality ψ k ( t , u ) = u 2 t , u R , t [ s k , t k ] holds. Then, the solution of Equation (11) is given by
u ( t ) = u 0 E 0.25 ( t 0.25 ) , t ( 0 , 1 ] , u 0 ( E 0.25 ( 1 ) ) k 2 t i = 1 k 1 ( 4 i ) , t ( 2 k 1 , 2 k ] , k = 1 , 2 , , u 0 ( E 0.25 ( 1 ) ) k i = 1 k ( 4 i ) E 0.25 ( ( t 2 k ) 0.25 ) , t ( 2 k , 2 k + 1 ] , k = 1 , 2 , .
The solution of Equation (11) is uniformly Lipschitz stable with M 1 = 30 (see Figure 1 for the graph of the solutions with various initial values).
Case 2. Suppose for all natural numbers k = 1 , 2 , the equality ψ k ( t , u ) = t u , u R , t [ s k , t k ] holds. Then, the solution of Equation (11) is given by
u ( t ) = u 0 E 0.25 ( t 0.25 ) , t ( 0 , 1 ] , u 0 ( E 0.25 ( 1 ) ) k i = 1 k 1 ( 2 i ) t , t ( 2 k 1 , 2 k ] , k = 1 , 2 , , u 0 ( E 0.25 ( 1 ) ) k i = 1 k ( 2 i ) E 0.25 ( ( t 2 k ) 0.25 ) , t ( 2 k , 2 k + 1 ] , k = 1 , 2 , .
The solution of Equation (11) is unbounded (see Figure 2 for the graph of the solution).
Therefore, for ψ k ( t , u ) = u 2 t u the solution is Lipschitz stable but for ψ k ( t , u ) = u 2 t u it is not (compare with condition (A8)).
In our study, we use some comparison results. When the Caputo fractional derivative is used, then the comparison result is:
Lemma 2.
(Caputo fractional derivative). Assume the following conditions are satisfied:
1. 
Assumptions A1–A4 and A6 are satisfied.
2. 
The function x * ( t ) = x ( t ; t 0 , φ ) : [ t 0 , T ) Δ , x * P C q ( [ t 0 , T ) ) is a solution of Equation (1) where Δ R n , 0 Δ , T .
3. 
The function V Λ ( [ t 0 , T ) , Δ ) is such that
(i) 
For any i = 0 , 1 , 2 , : ( t i , s i + 1 ) [ t 0 , T ) and for t ( t i , s i + 1 ) [ t 0 , T ) , the inequality
D t q t i C V ( t , x * ( t ) ) g ( t , V ( t , x * ( t ) ) )
holds.
(ii) 
For all i = 1 , 2 , 3 , : ( s i , t i ) [ t 0 , T ) the inequality
V ( t , ϕ i ( t , x * ( s i 0 ) ) ) ψ i ( t , V ( s i 0 , x * ( s i 0 ) ) ) f o r t ( s i , t i ] [ t 0 , T )
holds.
If s u p s [ r , 0 ] V ( t 0 , φ ( s ) ) u 0 , then the inequality V ( t , x * ( t ) ) r ( t ) for t [ t 0 , T ) holds, where r ( t ) = r ( t ; t 0 , u 0 ) is the maximal solution on [ t 0 , T ) of Equation (10) with u 0 0 .
Proof. 
We use induction with respect to the intervals to prove Lemma 2. Let m ( t ) = V ( t , x * ( t ) ) , t t 0 . We prove
m ( t ) u ( t ) , t t 0 .
Let t [ t 0 , s 1 ] . Let ε > 0 be an arbitrary number. We prove
m ( t ) < u ( t ) + ε , t [ t 0 , s 1 ] .
Note m ( t 0 ) = V ( t 0 , φ ( 0 ) ) sup s [ r , 0 ] V ( t 0 , φ ( s ) ) u 0 , i.e., the inequality in Equation (15) holds for t = t 0 . If the inequality in Equation (15) is not true, then there exists a point t * ( t 0 , s 1 ] such that m ( t * ) = u ( t * ) + ε , m ( t ) < u ( t ) + ε , t [ t 0 , t * ) .
From Lemma 1 with a = t 0 , b = s 1 , τ = t * and v ( t ) = m ( t ) u ( t ) ε the inequality D t q t 0 C m ( t * ) t 0 C D t q u ( t * ) = g ( t * , u ( t * ) ) holds.
From Assumption A6 and Condition 3(i), the inequality D t q t 0 C m ( t * ) g ( t * , m ( t * ) ) = g ( t * , u ( t * ) + ε ) < g ( t * , u ( t * ) ) holds. The contradiction proves the validity of Equation (15). Since ε is an arbitrary positive number, we obtain the inequality in Equation (14) for t [ t 0 , s 1 ] .
Let t ( s 1 , t 1 ] . Then, from the impulsive equality in Equation (1), Condition 3(ii), Assumption A6 and the inequality in Equation (14) for t = s 1 0 , we obtain m ( t ) = V ( t , x * ( t ) ) = V ( t , ϕ 1 ( t , x * ( s 1 0 ) ) ) ψ 1 ( t , V ( s 1 0 , x * ( s 1 0 ) ) ) = ψ 1 ( t , m ( s 1 0 ) ) ψ 1 ( t , u ( s 1 0 ) ) = u ( t ) , i.e., Equation (14) holds on ( s 1 , t 1 ] .
Let t ( t 1 , s 2 ] . Let ε > 0 be an arbitrary number. We prove Equation (15) for t [ t 1 , s 2 ] . Note that Equation (15) is true for t = t 1 . If the inequality in Equation (15) is not true, then there exists a point t * ( t 1 , s 2 ] such that m ( t * ) = u ( t * ) + ε , m ( t ) < u ( t ) + ε , t [ t 1 , t * ) .
From Lemma 1 with a = t 1 , b = s 2 , τ = t * and v ( t ) = m ( t ) u ( t ) ε , the inequality D t q t 0 C m ( t * ) t 1 C D t q u ( t * ) = g ( t * , u ( t * ) ) holds.
From Assumption A6 and Condition 3(i), the inequality D t q t 0 C m ( t * ) g ( t * , m ( t * ) ) = g ( t * , u ( t * ) + ε ) < g ( t * , u ( t * ) ) holds. The contradiction proves the validity of Equation (15) and the inequality in Equation (14) for t ( t 1 , s 2 ] . Continuing this process and an induction argument prove Equation (14) and Lemma 2. □
Lemma 3.
[10] Let m C ( [ t 0 , T ] , R ) and there exists τ ( t 0 , T ] , such that m ( τ ) = 0 and m ( t ) < 0 for t [ t 0 , τ ) . Then, the inequality D + q t 0     G L m ( τ ) > 0 holds.
When the Dini fractional derivative defined by Equation (4) or Caputo fractional Dini derivative defined by Equation (5) is used then the comparison result is:
Lemma 4.
(Dini fractional derivative/Caputo fractional Dini derivative). Assume:
1. 
Assumptions A1–A4 and A6 are satisfied.
2. 
The function x * ( t ) = x ( t ; t 0 , φ ) : [ t 0 , T ) Δ , x * P C q ( [ t 0 , T ) ) is a solution of Equation (1) where Δ R n , 0 Δ , T .
3. 
The function V Λ ( [ t 0 , T ) , Δ ) is such that
(i) 
For any i = 0 , 1 , 2 , : ( t i , s i + 1 ) [ t 0 , T ) and for t ( t i , s i + 1 ) [ t 0 , T ) , the inequality
D ( 1 ) V ( t , ϕ , t i ) g ( t , V ( t , ϕ ( 0 ) )
holds where ϕ ( Θ ) = x * ( t + Θ ) , Θ [ r , 0 ] , and D ( 1 ) V ( t , ϕ , t i ) is one of the following two derivatives: the Dini fractional derivative D ( 1 ) + V ( t , ϕ ( 0 ) , t i , ϕ ) defined by Equation (4) or the Caputo fractional Dini derivative D + q ( 1 ) c     V ( t , ϕ ; t i , φ ( 0 ) ) defined by Equation (5).
(ii) 
For all i = 1 , 2 , 3 , : ( s i , t i ) [ t 0 , T ) , the inequality
V ( t , ϕ i ( t , x * ( s i 0 ) ) ) ψ i ( t , V ( s i 0 , x * ( s i 0 ) ) ) f o r t ( s i , t i ] [ t 0 , T )
holds.
If s u p s [ r , 0 ] V ( t 0 , φ ( s ) ) u 0 , then the inequality V ( t , x * ( t ) ) r ( t ) for t [ t 0 , T ) holds, where r ( t ) = r ( t ; t 0 , u 0 ) is the maximal solution on [ t 0 , T ) of Equation (10) with u 0 0 .
Proof. 
The proof is similar to the one in Lemma 2 where instead of the Caputo fractional derivative of the Lyapunov function, we use the Dini fractional derivative or the Caputo fractional Dini derivative which are less restrictive with respect to the properties of Lyapunov functions (for example, differentiability is not required). We sketch the proof emphasizing the differences with Lemma 2.
Case 1. Let D ( 1 ) V ( t , ϕ , t i ) = D + q ( 1 ) c     V ( t , ϕ ; t i , φ ( 0 ) ) , i = 1 , 2 , in Condition 3(i) of Lemma 4.
We use induction with respect to the intervals to prove Lemma 4. We prove the inequality in Equation (14).
Case 1.1. Let t [ t 0 , s 1 ] . We prove Equation (15) with ε > 0 an arbitrary number. Note that Equation (15) holds for t = t 0 . If the inequality in Equation (15) is not true, then there exists a point t * ( t 0 , s 1 ] such that p ( t * ) = 0 and p ( t * ) < 0 for t [ t 0 , t * ) where p ( t ) = m ( t ) u ( t ) ε . From Lemma 3 with τ = t * we get the inequality
D + q t 0     G L m ( t * ) > t 0 G L D + q u ( t * ) + t 0 G L D + q ε .
Thus
D t q t 0 C m ( t * ) = t 0 G L D + q ( m ( t * ) m ( t 0 ) ) = t 0 G L D + q m ( t * ) t 0 G L D + q m ( t 0 )       > D t q t 0 C ( u ( t * ) u 0 ) = t 0 C D t q u ( t * ) .
Following the proof of Lemma 3 [3] from the choice of the point t * , the definition of the function m ( t ) , the definition of the derivative D + q ( 1 ) c     V ( τ , ϕ ( 0 ) ; t 0 , φ ( 0 ) ) , Assumption A2 and x ( t + s ) = ϕ ( s ) , x ρ ( t , x t ) = x ρ ( t , ϕ 0 ) = x ( ρ ( t , ϕ 0 ) ) = x ( t + ( ρ ( t , ϕ 0 ) t ) ) = ϕ ( ρ ( t , ϕ 0 ) t ) , Assumption A6 and Condition 3(i) of Lemma 4, we obtain the inequality
D t q t 0 C m ( τ ) = D + q t 0 G L ( m ( t * ) m ( t 0 ) ) D + q ( 1 ) c     V ( t * , ϕ ( 0 ) ; t 0 , φ ( 0 ) )       g ( t * , V ( t * , ϕ ( 0 ) ) = g ( t * , m ( t * ) ) = g ( t * , u ( t * ) + ε ) g ( t * , u ( t * ) )       = t 0 C D t q u ( t * )
with ϕ ( Θ ) = x ( τ + Θ ) , Θ [ τ , 0 ] .
The inequality in Equation (17) contradicts the inequality in Equation (16). The contradiction proves the validity of Equation (15) and, therefore, the validity of Equation (14) on [ t 0 , s 1 ] .
Case 1.2. Let t ( s 1 , t 1 ] . From the impulsive equality in Equation (1), Condition 3(ii) of Lemma 4, Assumption A6 and the inequality in Equation (14) for t = s 1 0 , we obtain for t ( s 1 , t 1 ] the inequalities m ( t ) = V ( t , x * ( t ) ) = V ( t , ϕ 1 ( t , x * ( s 1 0 ) ) ) ψ 1 ( t , V ( s 1 0 , x * ( s 1 0 ) ) ) = ψ 1 ( t , m ( s 1 0 ) ) ψ 1 ( t , u ( s 1 0 ) ) = u ( t ) , i.e., Equation (14) holds on ( s 1 , t 1 ) .
Case 1.3. Let t ( t 1 , s 2 ] . The proof of the inequality in Equation (15) for t ( t 1 , s 2 ] is similar to the one in Case 1.1 by replacing t 0 with t 1 .
Case 2. Let D ( 1 ) V ( t , ϕ , t i ) in Condition 3(i) of Lemma 4 be the Dini fractional derivative D ( 1 ) + V ( t , ϕ ( 0 ) , t i , φ ) defined by Equation (4). Then, based on the proof in Case 1 and Remark 4, we establish Lemma 4. □

6. Main Results

Theorem 1.
(Caputo fractional derivative) Let the following conditions be satisfied:
1. 
Assumptions A1–A8 are fulfilled.
2. 
There exist a function V Λ ( R + , R n ) and
(i) 
The inequalities
b ( | | x | | ) V ( t , x ) a ( | | x | | ) , x R n , t R +
holds, where a K ( [ 0 , ρ ] ) , b M ( [ 0 , ρ ] ) , ρ > 0 ;
(ii) 
For any initial data and any solution x ( t ) of Equation (1) defined on [ t 0 , ) such that for any τ ( t k , s k + 1 ) , k is a non-negative integer, such that x ( t ) S ρ , t [ t 0 , τ ] and V ( τ , x ( τ ) ) V ( s , x ( s ) ) for s [ t 0 , τ ] the inequality
D t q t i C V ( t , x ( t ) ) g ( t , V ( t , x ( t ) ) ) , t ( t i , s i + 1 ] [ t 0 , τ ] , i = 0 , 1 , 2 , , k
holds.
(iii) 
For any k = 0 , 1 , 2 , and t ( s k , t k + 1 ] , y S ρ the inequality
V ( t , ϕ k ( t , y ) ) ψ k ( t , V ( s k 0 , y ) )
holds.
3. 
The zero solution of Equation (10) is uniformly Lipschitz stable (uniformly globally Lipschitz stable).
Then, the zero solution of Equation (1) is uniformly Lipschitz stable (uniformly globally Lipschitz stable).
Proof. 
Let the zero solution of Equation (10) be uniformly Lipschitz stable. Let t 0 0 be an arbitrary. Without loss of generality, we assume t 0 [ 0 , s 1 ) . From Condition 3, there exist M 1 , δ 1 > 0 such that for any u 0 R : | u 0 | < δ 1 the inequality
| u ( t ; t 0 , u 0 ) | M | u 0 | for t t 0
holds, where u ( t ; t 0 , u 0 ) is a solution of Equation (10) with the initial data ( t 0 , u 0 ) .
From the inclusions a K ( [ 0 , ρ ] ) and b M ( [ 0 , ρ ] ) , there exist a function q b ( u ) and a positive constant K a . Without loss of generality, we can assume K a 1 . Choose the constant M 1 such that M 1 > max { 1 , q b ( K a ) , q b ( M ) K a } and δ 2 ρ 2 M 1 . Therefore, 2 M 1 δ 2 ρ .
Let δ = min δ 1 , δ 2 , δ 1 K a . Choose the initial function φ P C 0 ( [ r , 0 ] ) such that | | φ | | P C 0 < δ . Therefore, | | φ | | P C 0 < δ δ 2 ρ , i.e., φ ( s ) S ρ for s [ r , 0 ] . Consider the solution y ( t ) = y ( t ; t 0 , φ ) of the system in Equation (1) for the chosen initial data ( t 0 , φ ) .
Let u 0 * = sup s [ r , 0 ] V ( t 0 , φ ( s ) ) . From the choice of φ and the properties of the function a ( u ) applying condition 2 ( i ) we get u 0 * = V ( t 0 , φ ( ξ ) ) a ( | | φ ( ξ ) | | ) a ( | | φ | | P C 0 ) K a | | φ | | P C 0 < K a δ δ 1 . Therefore, the function u * ( t ) satisfies Equation (18) for t t 0 with u 0 = u 0 * , where u * ( t ) = u ( t ; t 0 , u 0 * ) is a solution of Equation (10) with initial data ( t 0 , u 0 * ) .
Let ε ( 0 , M 1 δ ] be an arbitrary number. We prove
V ( t , y ( t ) ) < b ( M 1 | | φ | | P C 0 + ε ) , t t 0 .
For t = t 0 , we get V ( t 0 , y ( t 0 ) ) = V ( t 0 , φ ( 0 ) ) a ( | | φ ( 0 ) | | ) a ( | | φ | | P C 0 K a | | φ | | P C 0 b ( q b ( K a ) | | φ | | P C 0 ) b ( M 1 | | φ | | P C 0 ) < b ( M 1 | | φ | | P C 0 + ε ) , i.e., the inequality in Equation (19) holds.
Assume Equation (19) is not true.
Case 1. There exists a point T > t 0 , T k = 0 ( t k , s k + 1 ] such that V ( t , y ( t ) ) < b ( M 1 | | φ | | P C 0 + ε ) for t [ t 0 , T ) , V ( T , y ( T ) ) = b ( M 1 | | φ | | P C 0 + ε ) , i.e., V ( s , y ( s ) ) V ( T , y ( T ) ) for s [ t 0 , T ] . Then, from Condition 2(i), we obtain the inequalities | | y ( t ) | | b 1 ( V ( t , y ( t ) ) M 1 | | φ | | P C 0 + ε < 2 M 1 δ 2 M 1 δ 2 ρ for t [ t 0 , T ] , i.e., y ( t ) S ρ for t [ t 0 , T ] and, according to Condition 2(ii) of Theorem 1 with τ = T , it follows that Condition 3(i) of Lemma 2 is satisfied for the solution y ( t ) on the interval [ t 0 , T ] and Δ = S ρ .
According to Lemma 2, we get
V ( t , y ( t ) ) u * ( t ) for t [ t 0 , T ] .
From the inequality in Equation (20) and Condition 2 ( i ) , we obtain
M 1 | | φ | | P C 0 = b 1 ( V ( T , y ( T ) ) ) b 1 ( u * ( T ) )        b 1 ( M | u 0 * | ) = b 1 ( M V ( t 0 , φ ( ξ ) ) ) q b ( M ) V ( t 0 , φ ( ξ ) )        q b ( M ) a ( | | φ ( ξ ) | | ) q b ( M ) a ( | | φ | | 0 ) q b ( M ) K a | | φ ) | | P C 0 < M 1 | | φ | | P C 0 .
The contradiction proves the validity of Equation (19). From the inequality in Equation (19) and Condition 2(i), we have Theorem 1.
Case 2. There exists a point T > t 0 , T k = 1 ( s k , t k ) such that V ( t , y ( t ) ) < b ( M 1 | | φ | | P C 0 + ε ) for t [ t 0 , T ) , V ( T , y ( T ) ) = b ( M 1 | | φ | | P C 0 + ε ) . Then, as in Case 1 we get y ( t ) S ρ for t [ t 0 , T ] . Let T ( s j , t j + 1 ) for a natural number j. According to Condition 2(iii) of Theorem 1, we obtain b ( M 1 | | φ | | P C 0 + ε ) = V ( T , y ( T ) ) = V ( T , ϕ j ( T , y ( s j 0 ) ) ) ψ j ( T , V ( s j 0 , y ( s j 0 ) ) ) ψ j ( T , b ( M 1 | | φ | | P C 0 ) ) < ψ j ( T , b ( M 1 | | φ | | P C 0 ) + ε ) . The contradiction proves this case is not possible.
Case 3. There exists a natural number k such that V ( t , y ( t ) ) < b ( M 1 | | φ | | 0 + ε ) for t [ t 0 , s k ] and V ( s k + 0 , y ( s k + 0 ) ) > b ( M 1 | | φ | | 0 + ε ) . Therefore, ψ k ( s k , b ( M 1 | | φ | | 0 ) ) ψ k ( s k , V ( s k , y ( s k ) ) ) V ( s k + 0 , ϕ k ( s k , y ( s k 0 ) ) ) = V ( s k + 0 , y ( s k + 0 ) ) > b ( M 1 | | φ | | 0 ) . The contradiction proves this case is not possible.
The proof of globally uniformly Lipschitz stability is analogous so we omit it. □
Theorem 2.
Let the conditions of Theorem 1 be satisfied where Condition 2(i) is replaced by:
2 * ( i ) the inequalities λ 1 ( t ) | | x | | 2 V λ 2 ( t ) | | x | | 2 , x S ρ , t R + holds, where λ 1 , λ 2 C ( R + , ( 0 , ) ) and there exists positive constant A 1 , A 2 : A 1 < A 2 such that λ 1 ( t ) A 1 , λ 2 ( t ) A 2 for t 0 , and ρ > 0 .
If the zero solution of Equation (10) is uniformly Lipschitz stable (uniformly globally Lipschitz stable), then the zero solution of Equation (1) is uniformly Lipschitz stable (uniformly globally Lipschitz stable).
Proof. 
The proof is similar to the one in Theorem 1 where M 1 = M A 2 A 1 .
Theorem 3.
(Dini fractional derivative/ Caputo fractional Dini derivative) Let the following conditions be satisfied:
1. 
Assumptions A1–A8 are fulfilled.
2. 
There exist a function V ( t , x ) Λ ( R + , R n ) , ρ > 0 and
(i) 
The inequalities
b ( | | x | | ) V ( t , x ) a ( | | x | | ) , x R n , t R +
holds, where a K ( [ 0 , ρ ] ) , b M ( [ 0 , ρ ] ) .
(ii) 
For any function ϕ P C 0 : ϕ ( s ) S ρ for s [ r , 0 ] such that for any t : t ( t k , s k + 1 ) , k is a non-negative integer, such that V ( t + s , ϕ ( s ) ) V ( t , ϕ ( 0 ) ) , s [ r , 0 ] the inequality
D ( 1 ) V ( t , ϕ , t k ) g ( t , V ( t , ϕ ( 0 ) ) )
holds where D ( 1 ) V ( t , ϕ , t k ) is one of the following two derivatives: the Dini fractional derivative D ( 1 ) + V ( t , ϕ ( 0 ) , t k , ϕ ) defined by Equation (4) or the Caputo fractional Dini derivative D + q ( 1 ) c     V ( t , ϕ ; t k , φ ( 0 ) ) defined by Equation (5) and ρ > 0 .
(iii) 
For any k = 0 , 1 , 2 , and t ( s k , t k + 1 ] , y S ρ the inequality
V ( t , ϕ k ( t , y ) ) ψ k ( t , V ( s k 0 , y ) )
holds.
3. 
The zero solution of Equation (10) is uniformly Lipschitz stable (uniformly globally Lipschitz stable).
Then, the zero solution of Equation (1) is uniformly Lipschitz stable (uniformly globally Lipschitz stable).
The proof of Theorem 3 is similar to the one in Theorem 1 where Lemma 4 is applied instead of Lemma 2.
Example 3.
Let t k = 2 k , k = 0 , 1 , 2 and s k = 2 k 1 , k = 1 , 2 , . Consider the non-instantaneous impulsive fractional differential equations
D t 0.25 t i C x 1 ( t ) = 0.25 x 1 ( t ) x 2 ( t ) + 0.25 x 1 ( t ) ( x ρ ( t , x t ) ) 2 2 , D t 0.25 t i C x 2 ( t ) = 0.25 x 2 ( t ) + x 1 ( t ) + 0.25 x 2 ( t ) ( x ρ ( t , x t ) ) 1 2 f o r t ( t i , s i + 1 ] , i = 0 , 1 , 2 , , x 1 ( t ) = x 1 ( s i 0 ) 2 i t , x 2 ( t ) = x 2 ( s i 0 ) 2 i t , t ( s i , t i ] , i = 1 , 2 , ,
where x = ( x 1 , x 2 ) , ρ ( t , u ) = t sin 2 ( u ) : t 1 ρ ( t , u ) t , x ρ ( t , x t ) = ( ( x ρ ( t , x t ) ) 1 , ( x ρ ( t , x t ) ) 2 ) and ( x ρ ( t , x t ) ) i = x i ( t sin 2 ( x i ( t + s ) ) ) , s [ 1 , 0 ] , i = 1 , 2 .
Let V ( t , x ) = x 1 2 + x 2 2 , x = ( x 1 , x 2 ) .
Let x ( t ) be a solution of Equation (22). Let the point τ ( t k , s k + 1 ] , k is a non-negative integer, be such that x ( t ) S 1 , t [ 0 , τ ] and x 1 ( τ ) 2 + x 2 ( τ ) 2 x 1 ( s ) 2 + x 2 ( s ) 2 , s [ 0 , τ ] . Using the notation x ρ ( τ , x τ ) and Assumption A2, it follows that ρ ( τ , x j ( τ + Θ ) ) [ τ r , τ ] , j = 1 , 2 , Θ [ r , 0 ] and therefore ( x ρ ( τ , x τ ) ) 1 2 + ( x ρ ( τ , x τ ) ) 2 2 x 1 ( τ ) 2 + x 2 ( τ ) 2 or
x 1 2 ( t ) ( x ρ ( t , x t ) ) 2 2 x 1 2 ( t ) ( x ρ ( t , x t ) ) 1 2 + x 1 2 ( t ) ( x ρ ( t , x t ) ) 2 2 x 1 2 ( t ) ( x 1 ( τ ) 2 + x 2 ( τ ) 2 ) x 1 ( t ) 2
and
x 2 2 ( t ) ( x ρ ( t , x t ) ) 1 2 x 2 2 ( t ) ( x ρ ( t , x t ) ) 1 2 + x 2 2 ( t ) ( x ρ ( t , x t ) ) 2 2 x 2 2 ( t ) ( x 1 ( τ ) 2 + x 2 ( τ ) 2 ) x 2 ( t ) 2 .
Then, for all i = 0 , 1 , 2 , , k ] and t ( t i , s i + 1 ] [ 0 , τ ] , we get the inequality
D t 0.25 t i C V ( t , x ( t ) ) = t i C D t 0.25 x 1 2 ( t ) + t i C D t 0.25 x 2 2 ( t )         2 x 1 ( t ) D t 0.25 t i C x 1 ( t ) + 2 x 2 ( t ) D t 0.25 t i C x 2 ( t )         = 2 x 1 ( t ) 0.25 x 1 ( t ) x 2 ( t ) + 0.25 x 1 ( t ) ( x ρ ( t , x t ) ) 2 2         + 2 x 2 ( t ) 0 . 25 x 2 ( t ) + x 1 ( t ) + 0.25 x 2 ( t ) ( x ρ ( t , x t ) ) 1 2         V ( t , x ( t ) ) .
In addition, for any natural number i, x S 1 R 2 and t [ s i , t i ] = [ 2 i 1 , 2 i , we get V ( t , x 2 i t ) = x 1 2 i t 2 + x 2 2 i t 2 = 1 2 t x 1 2 i + x 2 2 i x 1 2 + x 2 2 2 t = V ( s i , x ) 2 t = ψ i ( t , V ( s i , x ) ) with ψ i ( t , u ) = u 2 t .
According to Example 2, Case 1 and Theorem 1, the zero solution of Equation (22) is uniformly Lipschitz stable. □
Example 4.
Let t k = 2 k , k = 0 , 1 , 2 and s k = 2 k 1 , k = 1 , 2 , . Consider the non-instantaneous impulsive fractional differential equations
D t 0.25 t i C x 1 ( t ) = 0.5 x 1 ( t ) x 2 ( t ) + 0.5 x 1 ( t ) ( x ρ ( t , x t ) ) 2 2 x 1 ( t ) D q t i    R L ( cos ( 0.5 π ( t + t i + 1 ) ) + 1.1 ) cos ( 0.5 π ( t + t i + 1 ) ) + 1.1 , D t 0.25 t i C x 2 ( t ) = 0.5 x 2 ( t ) + x 1 ( t ) + 0.5 x 2 ( t ) ( x ρ ( t , x t ) ) 1 2 x 2 ( t ) D q t i    R L ( cos ( 0.5 π ( t + t i + 1 ) ) + 1.1 ) cos ( 0.5 π ( t + t i + 1 ) ) + 1.1 f o r t ( t i , s i + 1 ] , i = 0 , 1 , 2 , , x 1 ( t ) = x 1 ( s i 0 ) 2 i t , x 2 ( t ) = x 2 ( s i 0 ) 2 i t , t ( s i , t i ] , i = 1 , 2 , ,
where x = ( x 1 , x 2 ) , ρ ( t , u ) = t sin 2 ( u ) , t 0 . 5 ρ ( t , u ) t , x ρ ( t , x t ) = ( ( x ρ ( t , x t ) ) 1 , ( x ρ ( t , x t ) ) 2 ) and ( x ρ ( t , x t ) ) i = x i ( t 0.5 sin 2 ( x i ( t + s ) ) ) , s [ 0.5 , 0 ] , i = 1 , 2 , p ( t ) = cos ( 0.5 π ( t + t k + 1 ) ) + 1.1 for t [ t k , s k + 1 ] .
Note that, for any t [ t k , s k + 1 ] , the inequality p ( t ) p ( t + s ) , s [ 0 . 5 , 0 ] holds.
In this case, the quadratic function and Theorem 1 does not work (as it did in Example 3) because D t 0.25 t i C V ( t , x ( t ) ) 2 V ( t , x ( t ) ) ( 1 D q t i    R L p ( t ) p ( t ) ) ) 2 V ( t , x ( t ) ) ( 1 10 11 D q t i    R L ( p ( t ) ) and the solution of the comparison Equation (10) with g ( t , u ) = 2 u ( 1 10 21 D q t i    R L ( cos ( 0.5 π ( t + t k + 1 ) ) + 1.1 ) is difficult to obtain.
Consider the Lyapunov function V ( t , x ) = p ( t ) ( x 1 2 + x 2 2 ) , x = ( x 1 , x 2 ) .
Let the function ϕ P C 0 , r = 0.5 be such that ϕ ( s ) S 1 for s [ 0.5 , 0 ] . Let t : t ( t k , s k + 1 ) , k is a non-negative integer, be such that p ( t + s ) ( ϕ 1 ( s ) 2 + ϕ 2 ( s ) 2 ) p ( t ) ( ϕ 1 ( 0 ) 2 + ϕ 2 ( 0 ) 2 ) , s [ 1 , 0 ] . From the definition of the function ρ , it follows that ρ ( t , ϕ j ( s ) ) t = 0.5 sin 2 ( ϕ j ( s ) ) [ 0.5 , 0 ] for s [ 1 , 0 ] , j = 1 , 2 and therefore p ( t + s ) ( ( ϕ 1 ( ρ ( t , ϕ 1 ( s ) ) t ) ) 2 + ( ϕ 2 ( ρ ( t , ϕ 2 ( s ) ) t ) ) 2 ) p ( t ) ( ϕ 1 ( 0 ) 2 + ϕ 2 ( 0 ) 2 ) , s [ 0 . 5 , 0 ] . Then
p ( t ) ( ϕ 2 ( ρ ( t , ϕ 0 ) t ) ) 2 p ( t + s ) ( ϕ 2 ( ρ ( t , ϕ 2 ( s ) ) t ) ) 2    p ( t + s ) ( ϕ 1 ( ρ ( t , ϕ 1 ( s ) ) t ) ) 2 + ( ϕ 2 ( ρ ( t , ϕ 2 ( s ) ) t ) ) 2     p ( t ) ϕ 1 ( 0 ) 2 + ϕ 2 ( 0 ) 2 , s [ 0.5 , 0 ] .
Similarly, we get p ( t ) ( ϕ 1 ( ρ ( t , ϕ 0 ) t ) ) 2 p ( t ) ϕ 1 ( 0 ) 2 + ϕ 2 ( 0 ) 2 .
Then, using Example 1, Case 2 and the notations x t = ϕ 0 , x ρ ( t , x t ) = ϕ ( ρ ( t , ϕ 0 ) t ) , i.e., f ( t , x ( t ) , x ρ ( t , x t ) ) = f ( t , ϕ ( 0 ) , ϕ ( ρ ( t , ϕ 0 ) t ) ) , we get the inequality
D ( 24 ) + V ( t , ϕ ( 0 ) , t k , ϕ ) = ϕ 1 ( 0 ) p ( t ) f 1 ( t , ϕ ( 0 ) , ϕ ( ρ ( t , ϕ 0 ) t ) ) + ϕ 2 ( 0 ) p ( t ) f 2 ( t , ϕ ( 0 ) , ϕ ( ρ ( t , ϕ 0 ) t ) ) + ϕ 1 2 ( 0 ) + ϕ 2 2 ( 0 ) t k R L D q p ( t ) = ϕ 1 ( 0 ) p ( t ) ( 0.5 ϕ 1 ( 0 ) ϕ 2 ( 0 ) + 0.5 ϕ 1 ( 0 ) ( ϕ 2 ( ρ ( t , ϕ 0 ) t ) ) 2 ϕ 1 ( 0 ) D q t i    R L p ( t ) p ( t ) ) + ϕ 2 ( 0 ) p ( t ) ( 0.5 ϕ 2 ( 0 ) + ϕ 1 ( 0 ) + 0.5 ϕ 2 ( 0 ) ( ϕ 1 ( ρ ( t , ϕ 0 ) t ) ) 2 ϕ 2 ( 0 ) D q t i    R L p ( t ) p ( t ) ) + ϕ 1 2 ( 0 ) + ϕ 2 2 ( 0 ) t k R L D q p ( t ) V ( t , ϕ ( 0 ) ) .
In addition, for any natural number i, x S 1 R 2 and t [ s i , t i ] = [ 2 i 1 , 2 i , we get V ( t , x 2 i t ) = x 1 2 i t 2 + x 2 2 i t 2 = 1 2 t x 1 2 i + x 2 2 i x 1 2 + x 2 2 2 t = V ( s i , x ) 2 t = ψ i ( t , V ( s i , x ) ) with ψ i ( t , u ) = u 2 t .
According to Example 2, Case 1 and Theorem 3, the zero solution of Equation (24) is uniformly Lipschitz stable.

Author Contributions

Conceptualization, R.A., S.H. and D.O.; methodology, R.A., S.H. and D.O.; formal analysis, R.A., S.H. and D.O.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dannan, F.M.; Elaydi, S. Lipschitz stability of nonlinear systems of differential equations. J. Math. Anal. Appl. 1986, 113, 562–577. [Google Scholar] [CrossRef]
  2. Zavoli, A.; Giulietti, F.; Avanzini, G.; De Matteis, G. Spacecraft dynamics under the action of Y-dot magnetic control. Acta Astronaut. 2016, 122, 146–158. [Google Scholar] [CrossRef]
  3. Agarwal, R.; Hristova, S.; O’Regan, D. Lyapunov Functions and Stability of Caputo Fractional Differential Equations with Delays. Differ. Equ. Dyn. Syst. 2018, 1–22. [Google Scholar] [CrossRef]
  4. Hernandez, E.; O’Regan, D. On a new class of abstract impulsive differential equations. Proc. Am. Math. Soc. 2013, 141, 1641–1649. [Google Scholar] [CrossRef]
  5. Agarwal, R.; Hristova, S.; O’Regan, D. Non-Instantaneous Impulses in Differential Equations; Springer: Berlin, Germany, 2017. [Google Scholar]
  6. Agarwal, R.; Hristova, S.; O’Regan, D. Non-instantaneous impulses in Caputo fractional differential equations. Fract. Calc. Appl. Anal. 2017, 20, 595–622. [Google Scholar] [CrossRef] [Green Version]
  7. Bai, L.; Nieto, J.J. Variational approach to differential equations with not instantaneous impulses. Appl. Math. Lett. 2017, 73, 44–48. [Google Scholar] [CrossRef]
  8. Nieto, J.J.; Uzal, J.M. Pulse positive periodic solutions for some classes of singular nonlinearities. Appl. Math. Lett. 2018, 86, 134–140. [Google Scholar] [CrossRef]
  9. Agarwal, R.; Hristova, S.; O’Regan, D. Lyapunov Functions to Caputo Fractional Neural Networks with Time-Varying Delays. Axioms 2018, 7, 30. [Google Scholar] [CrossRef]
  10. Agarwal, R.P.; O’Regan, D.; Hristova, S. Stability of Caputo fractional differential equations by Lyapunov functions. Appl. Math. 2015, 60, 653–676. [Google Scholar] [CrossRef] [Green Version]
  11. Zhang, H.; Wu, D.; Cao, J. Asymptotic Stability of Caputo Type Fractional Neutral Dynamical Systems with Multiple Discrete Delays. Abstr. Appl. Analys. 2014, 2014, 138124. [Google Scholar] [CrossRef]
  12. Das, S. Functional Fractional Calculus; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  13. Diethelm, K. The Analysis of Fractional Differential Equations; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  14. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  15. Pandey, D.N.; Das, S.; Sukavanam, N. Existence of solutions for a second order neutral differential equation with state dependent delay and not instantaneous impulses. Intern. J. Nonlinear Sci. 2014, 18, 145–155. [Google Scholar]
  16. Lakshmikantham, V.; Bainov, D.; Simeonov, P. Theory of Impulsive Differential Equations; World Scientific: Singapore, 1989. [Google Scholar]
  17. Lakshmikantham, V.; Leela, S.; Devi, J.V. Theory of Fractional Dynamic Systems; CSP: Cambridge, UK, 2009; 170p. [Google Scholar]
  18. Sadati, S.J.; Ghaderi, R.; Ranjbar, A. Some fractional comparison results and stability theorem for fractional time delay systems. Rom. Rep. Phy. 2013, 65, 94–102. [Google Scholar]
  19. Stamova, I.; Stamov, G. Lipschitz stability criteria for functional differential systems of fractional order. J. Math. Phys. 2013, 54, 043502. [Google Scholar] [CrossRef]
  20. Stamova, I.; Stamov, G. Functional and Impulsive Differential Equations of Fractional Order: Qualitative Analysis and Applications; CRC Press: Boca Raton, FL, USA, 2016; 276p. [Google Scholar]
Figure 1. Example 2. Graph of the solution of Equation (11) with ψ k ( t , u ) = u 2 t for various initial values.
Figure 1. Example 2. Graph of the solution of Equation (11) with ψ k ( t , u ) = u 2 t for various initial values.
Axioms 08 00004 g001
Figure 2. Example 2. Graph of the solution of Equation (11) with ψ k ( t , u ) = t u .
Figure 2. Example 2. Graph of the solution of Equation (11) with ψ k ( t , u ) = t u .
Axioms 08 00004 g002

Share and Cite

MDPI and ACS Style

Agarwal, R.; Hristova, S.; O’Regan, D. Lipschitz Stability for Non-Instantaneous Impulsive Caputo Fractional Differential Equations with State Dependent Delays. Axioms 2019, 8, 4. https://doi.org/10.3390/axioms8010004

AMA Style

Agarwal R, Hristova S, O’Regan D. Lipschitz Stability for Non-Instantaneous Impulsive Caputo Fractional Differential Equations with State Dependent Delays. Axioms. 2019; 8(1):4. https://doi.org/10.3390/axioms8010004

Chicago/Turabian Style

Agarwal, Ravi, Snezhana Hristova, and Donal O’Regan. 2019. "Lipschitz Stability for Non-Instantaneous Impulsive Caputo Fractional Differential Equations with State Dependent Delays" Axioms 8, no. 1: 4. https://doi.org/10.3390/axioms8010004

APA Style

Agarwal, R., Hristova, S., & O’Regan, D. (2019). Lipschitz Stability for Non-Instantaneous Impulsive Caputo Fractional Differential Equations with State Dependent Delays. Axioms, 8(1), 4. https://doi.org/10.3390/axioms8010004

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