Next Article in Journal
Extremal Bipartite Graphs with Given Parameters on the Resistance–Harary Index
Next Article in Special Issue
Iteration Process for Fixed Point Problems and Zeros of Maximal Monotone Operators
Previous Article in Journal
Identities of Symmetry for Type 2 Bernoulli and Euler Polynomials
Previous Article in Special Issue
Non-Unique Fixed Point Theorems in Modular Metric Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Basic Concepts of Riemann–Liouville Fractional Differential Equations with Non-Instantaneous Impulses

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”, Plovdiv 4000, Bulgaria
4
School of Mathematics, Statistics and Applied Mathematics National University of Ireland, H91 TK33 Galway, Ireland
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(5), 614; https://doi.org/10.3390/sym11050614
Submission received: 27 March 2019 / Revised: 18 April 2019 / Accepted: 22 April 2019 / Published: 2 May 2019
(This article belongs to the Special Issue Fixed Point Theory and Fractional Calculus with Applications)

Abstract

:
In this paper a nonlinear system of Riemann–Liouville (RL) fractional differential equations with non-instantaneous impulses is studied. The presence of non-instantaneous impulses require appropriate definitions of impulsive conditions and initial conditions. In the paper several types of initial value problems are considered and their mild solutions are given via integral representations. In the linear case the equivalence of the solution and mild solutions is established. Conditions for existence and uniqueness of initial value problems are presented. Several examples are provided to illustrate the influence of impulsive functions and the interpretation of impulses in the RL fractional case.

1. Introduction

Fractional differential equations model nonlocal phenomena in time. One of the basic fractional derivatives is the Riemann–Liouville one which arises naturally in real world phenomena. For example, Heymans and Podlubny [1] provide several examples from the field of viscoelasticity. Several applications of fractional calculus to control theory, electrical circuits, fractional-order multipoles in electromagnetism, electrochemistry, and the neurons in biology are provided in [2,3,4,5]. In particular, several applications of fractional derivatives in physics are given in the book [6].
Many physical phenomena have short-term perturbations at some points caused by external interventions during their evolution. Adequate models for this kind of phenomena are impulsive differential equations. Two types of impulses are popular in the literature: instantaneous impulses (whose duration is negligible small) and non-instantaneous impulses (these changes start impulsively and remain active on finite initially given time intervals). There are mainly two approaches for the interpretation of the solutions of impulsive fractional differential equations: one by keeping the lower bound of the fractional derivative at the fixed initial time and the other by switching the lower limit of the fractional derivative at the impulsive points. The statement of the problem depends significantly on the type of fractional derivative. Caputo fractional derivatives have some properties similar to ordinary derivatives (such as the derivative of a constant) which lead to similar initial value problems as well as similar impulsive conditions (instantaneous and non-instantaneous). In the literature many types of initial value problems and boundary value problems for Caputo fractional differential equations with instantaneous and non-instantaneous impulses are studied (see, for example, [7,8,9]). For Riemann–Liouville (RL) fractional differential equations with instantaneous impulses several results are obtained in [8,10,11,12]. However for the RL fractional derivative, the physical interpretation of the initial condition (see [1]) requires different impulsive conditions, so as a result the statement of the problem considered is crucial. To the best of our knowledge this is the first paper concerning Riemann–Liouville fractional differential equations with non-instantaneous impulses.
In this paper the basic ideas in introducing non-instantaneous impulses for RL fractional differential equations are presented. The appropriate definition of both the initial conditions and the impulsive conditions for RL fractional differential equations is an important starting point for their qualitative investigation. We set up and discuss several types of initial conditions and impulsive conditions which are deeply connected with the RL fractional derivative. We use both the RL integral and the weighted limit to present the initial condition and the impulsive conditions. We consider both approaches in the literature in the fractional case with the presence of impulses in the equations: when the lower bound of the fractional derivative is fixed at the initial time and when the lower bound of the fractional derivative is changed at any point of impulse. In all cases mild solutions are defined by appropriate Volterra integro–algebraic representations and some conditions for existence and uniqueness are given. Note since instantaneous impulses are a special case of non-instantaneous ones, a brief overview of impulsive conditions in the instantaneous case is given.

2. Some Preliminary Results From Fractional Calculus

Let t 0 R + = [ 0 , ) be the initial time. Let L 1 loc ( J , R n ) be the linear space of all locally Lebesgue integrable functions m : J R m , J R . Let | | . | | be a norm in R n .
In this paper we will use the following definitions for fractional derivatives and integrals:
-
Riemann–Liouville fractional integral of order q ( 0 , 1 ) ([13])
t 0 I t q m ( t ) = 1 Γ ( q ) t 0 t m ( s ) ( t s ) 1 q d s , t t 0 ,
where m L 1 loc ( [ t 0 , ) , R ) and Γ ( . ) is the Gamma function.
This is called by some authors the left Riemann–Liouville fractional integral of order q.
Note sometimes the notation t 0 D t q m ( t ) = t 0 I t q m ( t ) is used.
-
Riemann–Liouville fractional derivative of order q ( 0 , 1 ) ([13])
t 0 R L D t q m ( t ) = d d t t 0 I t 1 q m ( t ) = 1 Γ 1 q d d t t 0 t t s q m ( s ) d s , t t 0
where m L 1 loc ( [ t 0 , ) , R ) .
This is also called the left Riemann–Liouville fractional derivative.
-
Caputo fractional derivative of order q ( 0 , 1 ) ([13])
t 0 C D t q m ( t ) = 1 Γ 1 q t 0 t t s q m ( s ) d s , t t 0
where m A C 1 ( [ t 0 , ) , R ) .
Proposition 1.
([13]) For q ( 0 , 1 ) , β > 0 the following hold:
t 0 I t q ( t t 0 ) β 1 = Γ ( β ) Γ ( β + q ) ( t t 0 ) β + q 1 ,
t 0 R L D t q ( t t 0 ) β 1 = Γ ( β ) Γ ( β q ) ( t t 0 ) β q 1 .
From Proposition 1 we have:
Corollary 1.
([13]) For q ( 0 , 1 ) :
t 0 I t q ( t t 0 ) q = Γ ( 1 q ) , t 0 R L D t q 1 = 1 Γ ( 1 q ) ( t t 0 ) q , t 0 R L D t q ( t t 0 ) q 1 = 0 .
Proposition 2.
([13]) For any m A C 1 ( [ t 0 , ) , R ) , m ( t 0 ) 0 we get
t 0 C D t q m ( t ) t 0 R L D t q m ( t ) .
Corollary 2.
([13]) The unique solution of 0 R L D t q u ( t ) = 0 is u ( t ) = c t q 1 .
The definitions of the initial condition of fractional differential equations with RL-derivatives are based on the following result:
Lemma 1.
(Lemma 3.2 [14]). Let q ( 0 , 1 ) and b R : t 0 < b , m : [ t 0 , b ] R be a Lebesgue measurable function.
(a) 
If there exists a limit lim t t 0 + [ ( t t 0 ) q 1 m ( t ) ] = c R , then there also exists a limit
t 0 I t 1 q m ( t ) | t = t 0 : = lim t t 0 + t 0 I t 1 q m ( t ) = c Γ ( q ) .
(b) 
If there exists a limit t 0 I t 1 q m ( t ) | t = t 0 = b R , and if the limit lim t t 0 + [ ( t t 0 ) 1 q m ( t ) ] exists, then
lim t t 0 + [ ( t t 0 ) 1 q m ( t ) ] = b Γ ( q ) .
Let a < b be real numbers and consider the nonlinear RL fractional differential equation (RLFrDE)
a R L D t q x ( t ) = F ( t , x ( t ) ) , t ( a , b ] .
Note that according to [14] the initial conditions to (1) could be:
-
the integral form (see (3.1.6) [14])
a I t 1 q x ( t ) | t = a = B R
-
a weighted Cauchy type problem (see (3.1.7) [14])
lim t a ( t a ) 1 q x ( t ) = C R
-
the initial condition at the inner point of a finite interval (see (3.4.71) [14])
x ( ξ ) = B , ξ ( a , b ) .
Remark 1.
According to Lemma 1 if the function x ( t ) satisfies the initial condition (3), then, x ( t ) also satisfies the condition (2) with B = C Γ ( q ) .
Any of the above initial value problem (IVP)’s of RLFrDE (1) has an equivalent integral representation.
Lemma 2.
Let G be an open set in R and F ( t , x ) : ( a , b ] × G R be such that F L ( a , b ) for any x G .
Then
-
the IVP for RLFrDE (1), (2) is equivalent to (see Corollary 3.1 [14] )
x ( t ) = B Γ ( q ) ( t a ) q 1 + 1 Γ ( q ) a t F ( s , x ( s ) ) ( t s ) 1 q d s , t ( a , b ] .
-
the weighted IVP for RLFrDE (1), (3) is equivalent to (see Corollary 3.1 [14] and Remark 1)
x ( t ) = C ( t a ) q 1 + 1 Γ ( q ) a t F ( s , x ( s ) ) ( t s ) 1 q d s , t ( a , b ] .
-
the IVP for RLFrDE (1), (4) is equivalent to (see (3.4.27) [14])
x ( t ) = ( ξ a t a ) 1 q B 1 Γ ( q ) a ξ F ( s , x ( s ) ) ( ξ s ) 1 q d s + 1 Γ ( q ) a t F ( s , x ( s ) ) ( t s ) 1 q d s , t ( ξ , b ] .
Lemma 3.
(Corollary 3.12 [14]). Let G be an open set in R and let F : ( a , b ] × G R ] F ( t , x ) C 1 q [ a , b ] (see Section 3) and for all t ( a , b ] and x , y G the inequality | F ( t , x ) F ( t , y ) | L | x y | , L > 0 holds.
Then there exists a unique solution to the integral type initial value problem for RLFrDE (1), (2) in the space C 1 q [ a , b ] .
Remark 2.
Note that the global existence result is true for the weighted Cauchy type problem for RLFrDE (1), (3) (see Theorem 3.12 [14]).
Lemma 4.
([14])
(i) 
If g C ( t 0 , T ] , then for any point t ( t 0 , T ]
t 0 R L D t q ( t 0 I t q g ( t ) ) = g ( t ) .
(ii) 
If g C ( t 0 , T ] and I t 0 1 q g ( t ) C ( t 0 , T ] , then for any point t ( t 0 , T ]
t 0 I t q ( t 0 R L D t q g ( t ) ) = g ( t ) t 0 I t 1 q g ( t ) | t = t 0 Γ ( q ) ( t t 0 ) q 1 .
Remark 3.
In the vector case, the fractional derivatives with the same fractional order are taken for all components.

3. Non-Instantaneous Impulses in RL Fractional Differential Equations

The appropriate definition of both the initial conditions and the impulsive conditions for RL fractional differential equations is an important starting point. We will discuss the initial value problem for nonlinear Riemann–Liouville fractional differential equations with non-instantaneous impulses (NIRLFrDE). We will set up in an appropriate way both the initial conditions and the impulsive conditions.
Note RL fractional functional differential equations with non-instantaneous impulses were studied in [15] but the impulsive conditions as well as the initial condition do not depend on the fractional order.
Let t 0 , T R + : t 0 < T < and points { t i } i = 1 p and { s i } i = 0 p be given such that 0 t 0 < s 0 < t i < s i < t i + 1 , i = 1 , 2 , p 1 , T = s p , and t 0 R + with p a natural number.
Definition 1.
The intervals ( s k 1 , t k ] , k = 1 , 2 , , p are called intervals of non-instantaneous impulses.
Remark 4.
If t k = s k 1 , k = 1 , 2 , , p , then the intervals of non-instantaneous impulses are reduced to points of instantaneous impulses.
We consider the following sets:
C 1 q [ t 0 , T ] = { y : [ t 0 , T ] R n : y ( t ) C [ t 0 , T ] and there exists t 0 I t q y ( t ) for t [ t 0 , T ] } , P C [ t 0 , T ] = { x : [ t 0 , T ] R n : x C ( [ t 0 , s 0 ) k = 1 p 1 ( s k , s k + 1 ) ) and x ( s k + 0 ) , x ( s k 0 ) exist for k = 0 , 1 , , p 1 } , P C 1 q [ t 0 , T ] = y ( t ) P C [ t 0 , T ] : ( t t 0 ) 1 q y ( t ) P C [ t 0 , T ] for t k = 0 p ( t k , s k ] , PC 1 q [ t 0 , T ] = { y ( t ) k = 0 p C ( t k , s k ) k = 0 p C ( s k , t k + 1 ) : ( t t k ) 1 q y ( t ) C ( t k , s k ) for k = 0 , 1 , , p and lim t t k + t k I t 1 q y ( t ) = lim t t k y ( t ) , k = 1 , 2 , , p }
with the norms
| | y | | P C 1 q [ t 0 , T ] = max { max k = 0 , 1 , , p sup t ( t k , s k ] | | ( t t 0 ) 1 q y ( t ) | | , max k = 0 , 1 , , p 1 sup t ( s k , t k + 1 ] | | y ( t ) | | }
and
| | y | | PC 1 q [ t 0 , T ] = max { max k = 0 , 1 , , p sup t ( t k , s k ] | | ( t t k ) 1 q y ( t ) | | , max k = 0 , 1 , , p 1 sup t ( s k , t k + 1 ] | | y ( t ) | | } .
When impulses are involved in fractional differential equations there are two main approaches for interpretation of the solutions. We will consider both approaches.

3.1. Fixed Lower Bound of the RL Fractional Derivative at the Given Initial Time

For our considerations in this case we need the function f ( t , x ) to be defined for all t [ t 0 , T ] , x R n .
As is mentioned in Section 2 there are two types of initial conditions to RL fractional differential equations: the integral form (2) and the weighted form (3). Following this idea we will define two different types of initial value problems for NIRLFrDE.

3.1.1. Integral Form of the Initial Value Problem

Consider the nonlinear non-instantaneous impulsive Riemann–Liouville fractional differential equation (NIRLFrDE)
t 0 R L D t q x ( t ) = f ( t , x ) f o r t ( t k , s k ] , k = 0 , 1 , , p ,
with impulsive conditions
x ( t ) = ϕ k ( t , x ( t ) , x ( s k 1 0 ) ) f o r t ( s k 1 , t k ] , k = 1 , 2 , , p ,
and the initial condition in integral form
t 0 I t 1 q x ( t ) | t = t 0 = x 0 ,
where x 0 R n , f : [ 0 , T ] × R n R n , ϕ k : [ s k 1 , t k ] × R n × R n R n , ( k = 1 , 2 , 3 , , p ) .
The impulses in problem (8), (9) start abruptly at the points s k , k = 0 , 1 , 2 , , p 1 and their action continues on the interval ( s k , t k + 1 ] . The function x takes an impulse at s k , k = 0 , 1 , 2 , , p 1 and it follows different rules in the two consecutive intervals ( s k , t k + 1 ] and ( t k + 1 , s k + 1 ] . At the point t k + 1 , k = 0 , 1 , 2 , , p 1 , the function x is continuous.
We will define a mild solution of the initial value problem (8), (9) by properly handling the RL fractional derivative and both the initial condition and the impulsive conditions.
Next we prove some auxiliary results.
Lemma 5.
Let h : [ t 0 , T ] R : h C 1 q [ t 0 , T ] , for all k = 1 , 2 , , p the functions ψ k C ( [ t k , s k ] × R ) and x : [ t 0 , T ] R be a function defined by
x ( t ) = x 0 Γ ( q ) ( t t 0 ) q 1 + 1 Γ ( q ) t 0 t ( t s ) q 1 h ( s ) d s , t ( t 0 , s 0 ] ψ k ( t , x ( s k 1 0 ) ) , t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p t k t 0 t t 0 1 q ψ k ( t k , x ( s k 1 0 ) ) t k t 0 t t 0 1 q 1 Γ ( q ) t 0 t k ( t k s ) q 1 h ( s ) d s + 1 Γ ( q ) t 0 t ( t s ) q 1 h ( s ) d s , t ( t k , s k ] , k = 1 , 2 , , p .
Then:
(i) 
the function x satisfies the linear problem
t 0 R L D t q x ( t ) = h ( t ) f o r t ( t k , s k ] , k = 0 , 1 , , p , x ( t ) = ψ k ( t , x ( s k 1 0 ) ) f o r t ( s k 1 , t k ] , k = 1 , 2 , , p , t 0 I t 1 q x ( t ) | t = t 0 = x 0 .
(ii) 
The function x is continuous at t i , i = 1 , 2 , , p .
Proof. 
(i). Let t ( t k , s k ] , k = 1 , 2 , , p (the proof in the case t ( t 0 , s 0 ] is obvious and we omit it). According to Lemma 4i and Proposition 1 the derivative t 0 R L D t q x ( t ) exists. Apply the operator t 0 R L D t q to both sides of (11), use Proposition 1, Lemma 4i and we obtain
t 0 R L D t q x ( t ) = ( t k t 0 ) 1 q ψ k ( t k , x ( s k 1 0 ) ) t 0 R L D t q ( t t 0 ) q 1 1 Γ ( q ) ( t k t 0 ) 1 q t 0 t k ( t k s ) q 1 h ( s ) d s t 0 R L D t q ( t t 0 ) q 1 + 1 Γ ( q ) t 0 R L D t q t 0 t ( t s ) q 1 h ( s ) d s = h ( t ) .
Therefore the function x ( t ) satisfies the first equation in (12).
Applying the operator t 0 I t 1 q x ( t ) to the first equation of (11), and using Proposition 1 and the equality t 0 I t 1 q t 0 I t q = t 0 I t 1 (see p. 10, (1.10) [16]), we get
t 0 I t 1 q x ( t ) = x 0 Γ ( q ) t 0 I t 1 q ( t t 0 ) q 1 + t 0 I t 1 q t 0 I t q h ( t ) = x 0 + 1 Γ ( 1 ) t 0 t h ( s ) d s .
Therefore, lim t t 0 t 0 I t 1 q x ( t ) = x 0 and the function x ( t ) satisfies the initial condition in (12).
(ii) Let k = 1 , 2 , , p . Taking the limit in (11), we obtain
lim t t k + x ( t ) = lim t t k + t k t 0 t t 0 1 q ψ k ( t k , x ( s k 1 0 ) ) lim t t k + t k t 0 t t 0 1 q 1 Γ ( q ) t 0 t k ( t k s ) q 1 h ( s ) d s + lim t t k + 1 Γ ( q ) t 0 t ( t s ) q 1 h ( s ) d s = ψ k ( t k , x ( s k 1 0 ) ) = x ( t k ) .
 □
Lemma 6.
Let h C 1 q [ t 0 , T ] . If x PC q [ t 0 , T ] and satisfies (12), then x satisfies (11).
Proof. 
Let t ( t 0 , s 0 ] and t 0 R L D t q x ( t ) = h ( t ) . Then according to Lemma 4ii we get
t 0 I t q h ( t ) = t 0 I t q t 0 R L D t q x ( t ) = x ( t ) t 0 I t 1 q x ( t ) | t = t 0 Γ ( q ) ( t t 0 ) q 1 ,
i.e.,
x ( t ) = t 0 I t q h ( t ) + t 0 I t 1 q x ( t ) | t = t 0 Γ ( q ) ( t t 0 ) q 1 = x 0 Γ ( q ) ( t t 0 ) q 1 + t 0 I t q h ( t ) ,
i.e., the function x ( t ) satisfies the first equation of (11).
Let t ( t 1 , s 1 ] and t 0 R L D t q x ( t ) = h ( t ) . According to Lemma 2 with a = t 0 , ξ = t 1 , b = s 1 , B = ψ 1 ( t 1 , x ( s 0 0 ) ) we get
x ( t ) = ( t 1 t 0 t t 0 ) 1 q B 1 Γ ( q ) t 0 t 1 h ( s ) ( t 1 s ) 1 q d s + 1 Γ ( q ) t 0 t h ( s ) ( t s ) 1 q d s ,
i.e., the last equality in (11) is satisfied.
Similarly, we can prove that x satisfies (11) for the other subintervals. □
Now, we introduce the concept of a mild solution for IVP for NIRLFrDE (8)–(10).
Definition 2.
A function x : [ t 0 , T ] R n is called a mild solution of the IVP for NIRLFrDE (8)–(10) if it satisfies the following Volterra integral–algebraic equation
x ( t ) = x 0 Γ ( q ) ( t t 0 ) q 1 + 1 Γ ( q ) t 0 t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t 0 , s 0 ] ϕ k ( t , x ( t ) , x ( s k 1 0 ) ) , t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p , t k t 0 t t 0 1 q ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) t k t 0 t t 0 1 q 1 Γ ( q ) t 0 t k ( t k s ) q 1 f ( s , x ( s ) ) d s + 1 Γ ( q ) t 0 t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t k , s k ] , k = 1 , 2 , , p .
Remark 5.
Note in formula (16) the function f ( t , x ) has to be defined on the whole interval [ t 0 , T ] .
Now we will establish existence results for mild solutions to the integral form of IVP for NIRLFrDE (8)–(10) on a finite interval.
Theorem 1.
(Existence and uniqueness). Let the following assumptions be satisfied:
1. 
The function f ( t , x ) C 1 q [ t 0 , T ] for any x R n and the inequality | | f ( t , x ) f ( t , y ) | | L | | x y | | holds for all t [ 0 , T ] and x , y R n , where L > 0 .
2. 
For all k = 1 , 2 , , p the functions ϕ k ( t , x , y ) C [ s k 1 , t k ] for any x , y R n and the inequality | | ϕ k ( t , x 1 , y 1 ) ϕ k ( t , x 2 , y 2 ) | | l k ( | | x 1 x 2 | | + | | y 1 y 2 | | ) holds for t [ s k 1 , t k ] and x 1 , x 2 , y 1 , y 2 R n where l k > 0 , k = 1 , 2 , , p .
3. 
The inequality K < 1 holds where
K = max { max k = 1 , 2 , p l k , 2 1 2 q L π ( s 0 t 0 ) q Γ ( 0 . 5 + q ) , l k ( t k t 0 ) 1 q + t k t 0 s k 1 t 0 1 q + L P k Γ ( q ) + L ( s k t k ) q Γ ( q + 1 ) s k t 0 t k t 0 1 q } , P k = max t [ t k , s k ] ( ( t k t 0 ) 1 q q i = 0 , 1 , , k ( t k t i ) q ( t k s i ) q ( t i t 0 ) 1 q ( t t 0 ) 1 q q i = 0 , 1 , , k ( t t i ) q ( t s i ) q ( t i t 0 ) 1 q + ( t k t 0 ) 1 q q i = 0 , 1 , , k 1 ( t k s i ) q ( t k t i + 1 ) q ( t t 0 ) 1 q q i = 0 , 1 , , k 1 ( t s i ) q ( t t i + 1 ) q ) , k = 1 , 2 , , p .
Then there exists a unique mild solution to the integral form of IVP for NIRLFrDE (8)–(10) in the space P C 1 q [ t 0 , T ] .
We will apply the Banach contraction principle. For any function x P C 1 q [ t 0 , T ] we define the operator
T x ( t ) = x 0 Γ ( q ) ( t t 0 ) q 1 + 1 Γ ( q ) t 0 t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t 0 , s 0 ] ϕ k ( t , x ( t ) , x ( s k 1 0 ) ) , t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p , t k t 0 t t 0 1 q ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) t k t 0 t t 0 1 q 1 Γ ( q ) t 0 t k ( t k s ) q 1 f ( s , x ( s ) ) d s + 1 Γ ( q ) t 0 t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t k , s k ] , k = 1 , 2 , , p .
From condition 1 it follows that the operator T is well defined.
Step 1. We prove that T x ( t ) P C 1 q [ t 0 , T ] for x P C 1 q [ t 0 , T ] .
From (17) and lim t t k + 0 T x ( t ) = lim t t k 0 T x ( t ) it follows that T x ( t ) C ( t 0 , T ] / { s k } .
Let t ( t 0 , s 0 ] . Then from condition 1 we get
( t t 0 ) 1 q T x ( t ) = x 0 Γ ( q ) + 1 Γ ( q ) t 0 t t t 0 t s 1 q f ( s , x ( s ) ) d s C ( t 0 , s 0 ] .
Let t ( t k , s k ] , k = 1 , 2 , , p . Then from (17) and condition 1 we get
( t t 0 ) 1 q T x ( t ) C [ t k , s k ) .
Step 2. The operator T is a contraction in P C 1 q [ t 0 , T ] .
Let x 1 , x 2 P C 1 q [ t 0 , T ] . First let t ( t 0 , s 0 ] . Then we obtain
sup t [ t 0 , s 0 ] | | ( t t 0 ) 1 q T x 1 ( t ) ( t t 0 ) 1 q T x 2 ( t ) | | sup t [ t 0 , s 0 ] L Γ ( q ) t 0 t t t 0 t s 1 q | | x 1 ( s ) x 2 ( s ) ) | | d s = sup t [ t 0 , s 0 ] L Γ ( q ) t 0 t t t 0 ( t s ) ( s t 0 ) 1 q | | ( s t 0 ) 1 q x 1 ( s ) ( s t 0 ) 1 q x 2 ( s ) ) | | d s sup t [ t 0 , s 0 ] L Γ ( q ) t 0 t t t 0 ( t s ) ( s t 0 ) 1 q | | ( s t 0 ) 1 q x 1 ( s ) ( s t 0 ) 1 q x 2 ( s ) ) | | d s | | x 1 x 2 | | P C 1 q [ t 0 , T ] L Γ ( q ) t 0 t t t 0 ( t s ) ( s t 0 ) 1 q d s = | | x 1 x 2 | | P C 1 q [ t 0 , T ] 2 1 2 q L π ( t t 0 ) q Γ ( 0 . 5 + q ) K | | x 1 x 2 | | P C 1 q [ t 0 , T ] .
Next let t ( s 0 , t 1 ] . Then
sup t ( s 0 , t 1 ] | | T x 1 ( t ) T x 2 ( t ) | | = sup t ( s 0 , t 1 ] | | ϕ 1 ( t , x 1 ( s 0 0 ) ) ϕ 1 ( t , x 2 ( s 0 0 ) ) | | l 1 | | x 1 x 2 | | P C 1 q [ t 0 , T ] K | | x 1 x 2 | | P C 1 q [ t 0 , T ] .
Continuing this procedure. For example let t ( t k , s k ] . Then we get
sup t ( t k , s k ] | | ( t t 0 ) 1 q T x 1 ( t ) ( t t 0 ) 1 q T x 2 ( t ) | | | | ( t k t 0 ) 1 q ϕ k ( t k , x 1 ( t k ) , x 1 ( s k 1 0 ) ) ϕ k ( t k , x 2 ( t k ) , x 2 ( s k 1 0 ) ) | | | | 1 Γ ( q ) t 0 t k t k t 0 t k s 1 q t t 0 t s 1 q ( f ( s , x 1 ( s ) ) f ( s , x 2 ( s ) ) ) d s | | + 1 Γ ( q ) t k t t t 0 t s 1 q | | f ( s , x 1 ( s ) ) f ( s , x 2 ( s ) ) | | d s l k ( t k t 0 ) 1 q | | x 1 ( t k ) x 2 ( t k ) | | + l k t k t 0 s k 1 t 0 1 q ( s k 1 t 0 ) 1 q | | x 1 ( s k 1 0 ) x 2 ( s k 1 0 ) | | + L Γ ( q ) t 0 t k | t k t 0 t k s 1 q t t 0 t s 1 q | | | x 1 ( s ) ) x 2 ( s ) | | d s + L Γ ( q ) t k t t t 0 ( t s ) ( t k t 0 ) 1 q | | ( s t 0 ) 1 q x 1 ( s ) x 2 ( s ) ) | | d s K | | x 1 x 2 | | P C 1 q [ t 0 , T ]
where the inequalities
| | x 1 ( t k ) x 2 ( t k ) | | sup t ( s k 1 , t k ] | | x 1 ( t ) x 2 ( t ) | | | | x 1 x 2 | | P C 1 q [ t 0 , T ] , | | ( s k 1 t 0 ) 1 q ( x 1 ( s k 1 0 ) x 2 ( s k 1 0 ) ) | | sup t ( t k 1 , s k 1 ] | | ( t t 0 ) 1 q ( x 1 ( t ) x 2 ( t ) ) | | | | x 1 x 2 | | P C 1 q [ t 0 , T ] , s i t i + 1 ( t s ) q 1 d s = ( t s i ) q ( t t i + 1 ) q q , s i t i + 1 ( t k s ) q 1 d s = ( t k s i ) q ( t k t i + 1 ) q q , t i s i ( t s ) q 1 d s = ( t t i ) q ( t s i ) q q , t i s i ( t k s ) q 1 d s = ( t k t i ) q ( t k s i ) q q , t k t ( t s ) q 1 d s = ( t t k ) q q ,
t k t t t 0 ( t s ) ( t k t 0 ) 1 q d s = t t 0 t k t 0 1 q ( t t k ) q q s k t 0 t k t 0 1 q ( s k t k ) q q
and
t 0 t k t k t 0 t k s 1 q t t 0 t s 1 q | | x 1 ( s ) x 2 ( s ) | | d s i = 0 , 1 , , k t i s i t k t 0 ( t k s ) ( t i t 0 ) 1 q t t 0 ( t s ) ( t i t 0 ) 1 q × × | | ( s t 0 ) 1 q ( x 1 ( s ) x 2 ( s ) ) | | d s + i = 0 , 1 , , k 1 s i t i + 1 t k t 0 t k s 1 q t t 0 t s 1 q | | x 1 ( s ) x 2 ( s ) ) | | d s ( ( t k t 0 ) 1 q q ( i = 0 , 1 , , k ( t k t i ) q ( t k s i ) q ( t i t 0 ) 1 q + i = 0 , 1 , , k 1 ( t k s i ) q ( t k t i + 1 ) q ( t t 0 ) 1 q q i = 0 , 1 , , k ( t t i ) q ( t s i ) q ( t i t 0 ) 1 q + i = 0 , 1 , , k 1 ( t s i ) q ( t t i + 1 ) q × × | | x 1 x 2 | | P C 1 q [ t 0 , T ] P k | | x 1 x 2 | | P C 1 q [ t 0 , T ]
are used.
From the inequalities (18)–(20) it follows that | | T x 1 T x 2 | | P C 1 q [ t 0 , T ] K | | x 1 x 2 | | P C 1 q [ t 0 , T ] which proves the Theorem.

3.1.2. Weighted Initial Value Problem

Consider the NIRLFrDE (8) with impulsive conditions (9) and the weighted Cauchy type condition
lim t t 0 ( t t 0 ) 1 q x ( t ) = x 0 R .
Then the following auxiliary results hold.
Lemma 7.
Let h : [ t 0 , T ] R : h C 1 q [ t 0 , T ] , for all k = 1 , 2 , , p the functions ψ k C ( [ t k , s k ] × R ) and x : [ 0 , T ] R be a function defined by
x ( t ) = x 0 ( t t 0 ) q 1 + 1 Γ ( q ) t 0 t ( t s ) q 1 h ( s ) d s , t ( t 0 , s 0 ] ψ k ( t , x ( s k 1 0 ) ) , t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p , t k t 0 t t 0 1 q ψ k ( t k , x ( s k 1 0 ) ) t k t 0 t t 0 1 q 1 Γ ( q ) t 0 t k ( t k s ) q 1 h ( s ) d s + 1 Γ ( q ) t 0 t ( t s ) q 1 h ( s ) d s , t ( t k , s k ] , k = 1 , 2 , , p .
Then:
(i) 
the function x satisfies the linear problem
t 0 R L D t q x ( t ) = h ( t ) f o r t ( t k , s k ] , k = 0 , 1 , , p , x ( t ) = ψ k ( t , x ( s k 1 0 ) ) f o r t ( s k 1 , t k ] , k = 1 , 2 , , p , lim t t 0 ( t t 0 ) 1 q x ( t ) = x 0 .
(ii) 
The function x is continuous at t i , i = 1 , 2 , , p .
The proof follows from Lemma 5 and Remark 1 so we omit it.
Lemma 8.
Let h C 1 q [ t 0 , T ] . If x PC q [ t 0 , T ] and satisfies (23), then x satisfies (22).
Now, we introduce the concept of a mild solution for IVP for NIRLFrDE (8), (9), (21).
Definition 3.
A function x : [ t 0 , T ] R n is called a mild solution of the IVP for NIRLFrDE (8), (9), (21) if it satisfies the following Volterra integral–algebraic equation (compare with (16)):
x ( t ) = x 0 ( t t 0 ) q 1 + 1 Γ ( q ) t 0 t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t 0 , s 0 ] ϕ k ( t , x ( t ) , x ( s k 1 0 ) ) , t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p , t k t 0 t t 0 1 q ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) t k t 0 t t 0 1 q 1 Γ ( q ) t 0 t k ( t k s ) q 1 f ( s , x ( s ) ) d s + 1 Γ ( q ) t 0 t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t k , s k ] , k = 1 , 2 , , p .
Remark 6.
The existence and uniqueness Theorem 1 is true for the mild solution to the weighted form of IVP for NIRLFrDE (8), (9), (21) defined by (24).

3.2. Changed Lower Bounds of the RL Fractional Derivative at the Impulsive Points

Here our function f ( t , x ) is defined only for t k = 0 p [ t k , s k ] and x R n .
In this case, we will consider the RL fractional derivative with a changeable lower bound at any point of jump, i.e., instead of t 0 R L D t q for all t > t 0 we consider t k R L D t q on ( t k , s k ] , k = 0 , 1 , . It does not seem to be possible to consider the usual Cauchy conditions at the point t k to the RLFrDE t k R L D t q x ( t ) = f ( t , x ( t ) ) , because the solution to this problem, in general, has a singularity at t k and therefore it is not bounded and continuous at the point t k (see Section 3.4.2 [14]). There are two types of initial conditions to RL fractional differential equations: the integral form (2) and the weighted form (3). Following this idea we will define four different types of initial value problems for NIRLFrDE.

3.2.1. Integral form of the Initial Conditions and Impulses

Consider the IVP for NIRLFrDE
t k R L D t q x ( t ) = f ( t , x ) f o r t ( t k , s k ] , k = 0 , 1 , , p ,
with impulsive conditions in the form of RL integrals
x ( t ) = ϕ k ( t , x ( t ) , x ( s k 1 0 ) ) f o r t ( s k 1 , t k ) , k = 1 , 2 , , p , t k I t 1 q x ( t ) | t = t k = ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) , k = 1 , 2 , , p ,
and initial conditions
t 0 I t 1 q x ( t ) | t = t 0 = x 0 ,
where x 0 R n , f : k = 0 p [ t k , s k ] × R n R n , ϕ k : [ s k 1 , t k ] × R n × R n R n , ( k = 1 , 2 , , p ) .
Lemma 9.
Let h : k = 0 p [ t k , s k ] R : h k = 0 p C 1 q [ t k , s k ] , for all k = 1 , 2 , , p the functions ψ k C ( [ t k , s k ] × R ) and x : [ t 0 , T ] R be a function defined by
x ( t ) = x 0 Γ ( q ) ( t t 0 ) q 1 + 1 Γ ( q ) t 0 t ( t s ) q 1 h ( s ) d s , t ( t 0 , s 0 ] ψ k ( t , x ( s k 1 0 ) ) , t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p ψ k ( t k , x ( s k 1 0 ) ) Γ ( q ) ( t t k ) q 1 + 1 Γ ( q ) t k t ( t s ) q 1 h ( s ) d s , t ( t k , s k ] , k = 1 , 2 , , p .
Then:
(i) 
the function x satisfies the linear problem
t k R L D t q x ( t ) = h ( t ) f o r t ( t k , s k ] , k = 0 , 1 , , p , x ( t ) = ψ k ( t , x ( s k 1 0 ) ) f o r t ( s k 1 , t k ) , k = 1 , 2 , , p , t k I t 1 q x ( t ) | t = t k = ψ k ( t k , x ( s k 1 0 ) ) , k = 1 , 2 , , p t 0 I t 1 q x ( t ) | t = t 0 = x 0 .
(ii) 
The function x satisfies lim t t k + t k I t 1 q x ( t ) = lim t t k x ( t ) for k = 1 , 2 , , p .
Proof. 
(i). Let t ( t k , s k ) , k = 1 , 2 , , p . According to Lemma 4i and Proposition 1 the derivative t k R L D t q x ( t ) exists. Applying the operator t k R L D t q to both sides of the last equality of (28), and useing Corollary 1 and Lemma 4i, we obtain
t k R L D t q x ( t ) = ψ k ( t k , x ( s k 1 0 ) ) Γ ( q ) t k R L D t q ( t t k ) q 1 + 1 Γ ( q ) t k R L D t q t k t ( t s ) q 1 h ( s ) d s = h ( t ) .
Therefore the function x ( t ) satisfies the first equation in (29). A similar argument is needed for t ( t 0 , s 0 ] .
Let t ( t k , s k ) , k = 1 , 2 , , p and apply the operator t k I t 1 q x ( t ) to both sides of the last equality of (28), use Corollary 1 and the equality t k I t 1 q t k I t q = t k I t 1 (see p. 10, (1.10) [16]), to get
t k I t 1 q x ( t ) = ψ k ( t k , x ( s k 1 0 ) ) Γ ( q ) t k I t 1 q ( t t k ) q 1 + t k I t 1 q t k I t q h ( t ) = ψ k ( t k , x ( s k 1 0 ) ) + 1 Γ ( 1 ) t k t h ( s ) d s .
Therefore, lim t t k t k I t 1 q x ( t ) = ψ k ( t k , x ( s k 1 0 ) ) and the function x ( t ) satisfies the impulsive condition in (29).
Apply the operator t 0 I t 1 q x ( t ) to the first equation of (28), use Corollary 1 and the equality t k I t 1 q t k I t q = t 0 I t 1 (see p. 10, (1.10) [16]) and we get
t k I t 1 q x ( t ) = x 0 Γ ( q ) t k I t 1 q ( t t k ) q 1 + t k I t 1 q t 0 I t q h ( t ) = x 0 + 1 Γ ( 1 ) t 0 t h ( s ) d s .
Therefore, lim t t 0 t 0 I t 1 q x ( t ) = x 0 and the function x ( t ) satisfies the initial condition in (29).
(ii) Let k = 1 , 2 , , p . According to (31) we obtain
lim t t k + x ( t ) = lim t t k + t k I t 1 q x ( t ) = ψ k ( t k , x ( s k 1 0 ) = lim t t k ψ k ( t , x ( s k 1 0 ) ) = lim t t k x ( t ) .
 □
Lemma 10.
Let h k = 0 p C 1 q [ t k , s k ] . If x PC q [ t 0 , T ] and satisfies (29), then x satisfies (28).
Proof. 
Let t ( t 0 , s 0 ] and t 0 R L D t q x ( t ) = h ( t ) . Then according to Lemma 4ii we get
t 0 I t q h ( t ) = t 0 I t q t 0 R L D t q x ( t ) = x ( t ) t 0 I t 1 q x ( t ) | t = t 0 Γ ( q ) ( t t 0 ) q 1 ,
i.e.,
x ( t ) = t 0 I t q h ( t ) + t 0 I t 1 q x ( t ) | t = t 0 Γ ( q ) ( t t 0 ) q 1 = x 0 Γ ( q ) ( t t 0 ) q 1 + t 0 I t q h ( t ) ,
i.e., the function x ( t ) satisfies the first equation of (28).
Let t ( t 1 , s 1 ] and t 1 R L D t q x ( t ) = h ( t ) . According to Lemma 2 and (6) with a = t k , b = s k , B = ψ 1 ( t 1 , x ( s 0 0 ) ) we get
x ( t ) = ψ 1 ( t 1 , x ( s 0 0 ) ) Γ ( q ) ( t t 1 ) q 1 + 1 Γ ( q ) t 1 t f ( s , x ( s ) ) ( t s ) 1 q d s ,
i.e., the last equality in (28) is satisfied for k = 1 .
Similarly, we can prove that x satisfies (28) for the other subintervals. □
Now, we define a mild solution for IVP for NIRLFrDE (25)–(27) following the idea with the presence of impulses in differential equations ([17]) and Lemma 2 with t 0 = t k , k = 1 , 2 , and b = ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) .
Definition 4.
A function x : [ t 0 , T ] R n is called a mild solution of the IVP for NIRLFrDE (25)–(27) if it satisfies the following Volterra integral–algebraic equations:
x ( t ) = x 0 Γ ( q ) ( t t 0 ) q 1 + 1 Γ ( q ) t 0 t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t 0 , s 0 ] ϕ k ( t , x ( t ) , x ( s k 1 0 ) ) , t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p , ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) Γ ( q ) ( t t k ) q 1 + 1 Γ ( q ) t k t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t k , s k ] , k = 1 , 2 , , p .
We give conditions for the existence of a mild solution to the integral type IVP for NIRLFrDE (25)–(27).
Theorem 2.
(Existence and uniqueness). Let the following assumptions be satisfied:
1. 
The function f ( t , x ) : k = 0 p [ t k , s k ] × R n R n , f ( t , x ) k = 0 p C 1 q [ t k , s k ] for any x R n and the inequality | | f ( t , x ) f ( t , y ) | | L | | x y | | holds for all t [ t k , s k ] , k = 0 , 1 , , p and x , y R n , where L > 0 .
2. 
For all k = 1 , 2 , , p the functions ϕ k ( t , x , y ) C [ s k 1 , t k ] for any x , y R n and the inequality | | ϕ k ( t , x 1 , y 1 ) ϕ k ( t , x 2 , y 2 ) | | l k ( | | x 1 x 2 | | + | | y 1 y 2 | | ) holds for t [ s k 1 , t k ] and x 1 , x 2 , y 1 , y 2 R n where l k > 0 , k = 1 , 2 , , p .
3. 
The inequality K < 1 holds where
K = max { max k = 1 , 2 , p l k , 2 1 2 q L π ( s 0 t 0 ) q Γ ( 0 . 5 + q ) , l k Γ ( q ) + l k Γ ( q ) ( s k 1 t k 1 ) 1 q + L 2 1 2 q π ( t t k ) q Γ ( 0 . 5 + q ) } .
Then there exists a unique mild solution to the integral form of IVP for NIRLFrDE (8)–(10) in the space PC 1 q [ t 0 , T ] .
We will apply the Banach contraction principle. For any function x PC 1 q [ t 0 , T ] we define the operator
T x ( t ) = x 0 Γ ( q ) ( t t 0 ) q 1 + t 0 I t q f ( t , x ( t ) ) , t ( t 0 , s 0 ] ϕ k ( t , x ( t ) , x ( s k 1 0 ) ) , t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p , ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) Γ ( q ) ( t t k ) q 1 + t k I t q f ( t , x ( t ) ) , t ( t k , s k ] , k = 1 , 2 , , p .
From condition 1 it follows that the operator T is well defined.
Step 1. We prove that T x ( t ) PC 1 q [ t 0 , T ] for x PC 1 q [ t 0 , T ] .
From (34) it follows that T x ( t ) k = 0 p C ( t k , s k ) k = 0 p C ( s k , t k + 1 ) .
Let t ( t 0 , s 0 ] . Then from condition 1 we get
( t t 0 ) 1 q T x ( t ) = x 0 Γ ( q ) + 1 Γ ( q ) t 0 t t t 0 t s 1 q f ( s , x ( s ) ) d s C ( t 0 , s 0 ] .
Let t ( t k , s k ] , k = 1 , 2 , , p . Then from (34) and condition 1 we get
( t t k ) 1 q T x ( t ) C ( t k , s k ) .
From (34), condition 1, Corollary 1 and t k I t 1 q t k I t q = t k I t 1 we get
t k I t 1 q T x ( t ) = ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) Γ ( q ) t k I t 1 q ( t t k ) q 1 + t k I t 1 q t k I q f ( t , x ( t ) ) = ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) + t k I t 1 f ( t , x ( t ) ) .
Therefore,
lim t t k + t k I t 1 q x ( t ) = lim t t k ϕ k ( t , x ( t ) , x ( s k 1 0 ) = lim t t k x ( t ) .
Step 2. The operator T is a contraction in PC 1 q [ t 0 , T ] .
Let x 1 , x 2 PC 1 q [ t 0 , T ] . First let t ( t 0 , s 0 ] . Then similar to (18) we get
sup t [ t 0 , s 0 ] | | ( t t 0 ) 1 q T x 1 ( t ) ( t t 0 ) 1 q T x 2 ( t ) | | | | x 1 x 2 | | PC 1 q [ t 0 , T ] 2 1 2 q L π ( t t 0 ) q Γ ( 0 . 5 + q ) K | | x 1 x 2 | | PC 1 q [ t 0 , T ]
where the equality
a t 1 ( t s ) ( s a ) 1 q d s = 2 1 2 q π ( t a ) 2 q 1 Γ ( q ) Γ ( 0 . 5 + q )
is applied. Next let t ( s 0 , t 1 ] . As in (19) we get
sup t ( s 0 , t 1 ] | | T x 1 ( t ) T x 2 ( t ) | | l 1 | | x 1 x 2 | | PC 1 q [ t 0 , T ] K | | x 1 x 2 | | PC 1 q [ t 0 , T ] .
Continuing this procedure. For example let t ( t k , s k ] . Then from (34) we get
sup t ( t k , s k ] | | ( t t k ) 1 q T x 1 ( t ) ( t t k ) 1 q T x 2 ( t ) | | | | ϕ k ( t k , x 1 ( t k ) , x 1 ( s k 1 0 ) ) ϕ k ( t k , x 2 ( t k ) , x 2 ( s k 1 0 ) ) Γ ( q ) | | + 1 Γ ( q ) t k t t t k t s 1 q | | f ( s , x 1 ( s ) ) f ( s , x 2 ( s ) ) | | d s l k Γ ( q ) | | x 1 ( t k ) x 2 ( t k ) | | + l k Γ ( q ) ( s k 1 t k 1 ) 1 q ( s k 1 t k 1 ) 1 q | | x 1 ( s k 1 0 ) x 2 ( s k 1 0 ) | | + L Γ ( q ) t k t t t k ( t s ) ( s t k ) 1 q | | ( s t k ) 1 q ( x 1 ( s ) x 2 ( s ) ) | | d s l k Γ ( q ) + l k Γ ( q ) ( s k 1 t k 1 ) 1 q + L 2 1 2 q π ( t t k ) q Γ ( 0 . 5 + q ) | | x 1 x 2 | | PC 1 q [ t 0 , T ] K | | x 1 x 2 | | PC 1 q [ t 0 , T ]
where the inequality (36) with a = t k and
| | x 1 ( t k ) x 2 ( t k ) | | sup t ( s k 1 , t k ] | | x 1 ( t ) x 2 ( t ) | | | | x 1 x 2 | | PC 1 q [ t 0 , T ] , | | ( s k 1 t 0 ) 1 q ( x 1 ( s k 1 0 ) x 2 ( s k 1 0 ) ) | | sup t ( t k 1 , s k 1 ] | | ( t t 0 ) 1 q ( x 1 ( t ) x 2 ( t ) ) | | | | x 1 x 2 | | PC 1 q [ t 0 , T ]
are used.
From the inequalities (35), (37), (38) it follows that | | T x 1 T x 2 | | P C 1 q [ t 0 , T ] K | | x 1 x 2 | | P C 1 q [ t 0 , T ] which proves the Theorem.

3.2.2. Weighted Form of the Initial Conditions and Impulses

Consider the NIRLFrDE (25) with impulsive conditions in weighted form
x ( t ) = ϕ k ( t , x ( t ) , x ( s k 1 0 ) ) f o r t ( s k 1 , t k ) , k = 1 , 2 , , p , lim t t k ( t t k ) 1 q x ( t ) = ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) , k = 1 , 2 , , p
and initial conditions in the weighted form:
lim t t 0 ( t t 0 ) 1 q x ( t ) = x 0 .
Applying Lemmas 2, 9 and 10 and (6) we could define a mild solution of IVP for NIRLFrDE (25), (39), (40):
Definition 5.
A function x : [ t 0 , T ] R n is called a mild solution of the IVP for NIRLFrDE (25), (39), (40) if it satisfies the following Volterra integral–algebraic equations:
x ( t ) = x 0 ( t t 0 ) q 1 + 1 Γ ( q ) t 0 t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t 0 , s 0 ] ϕ k ( t , x ( t ) , x ( s k 1 0 ) ) , t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p , ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) ( t t k ) q 1 + 1 Γ ( q ) t k t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t k , s k ] , k = 1 , 2 , , p .
Note Theorem 2 for existence and uniqueness could be easily converted to the mild solution of the IVP for NIRLFrDE (25), (39), (40) given by (41).

3.2.3. Mixed Forms of the Initial Conditions and Impulses

We have the following cases:
Case 1. Consider the IVP for NIRLFrDE (25) with mixed impulsive and initial conditions of the form:
x ( t ) = ϕ k ( t , x ( t ) , x ( s k 1 0 ) ) f o r t ( s k 1 , t k ) , k = 1 , 2 , , p , t k I t 1 q x ( t ) | t = t k = ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) , k = 1 , 2 , , p
and initial conditions in the weighted form:
lim t t 0 ( t t 0 ) 1 q x ( t ) = x 0 ,
Then the mild solution to the IVP (25), (42), (43) will be:
Definition 6.
A function x : [ t 0 , T ] R n is called a mild solution of the IVP for NIRLFrDE (25), (42), (43) if it satisfies the following Volterra integral–algebraic equations:
x ( t ) = x 0 ( t t 0 ) q 1 + 1 Γ ( q ) t 0 t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t 0 , s 0 ] ϕ k ( t , x ( t ) , x ( s k 1 0 ) ) , t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p , ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) Γ ( q ) ( t t k ) q 1 + 1 Γ ( q ) t k t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t k , s k ] , k = 1 , 2 , , p .
Case 2. Consider the IVP for NIRLFrDE (25) with impulsive conditions in the weighted form:
x ( t ) = ϕ k ( t , x ( t ) , x ( s k 1 0 ) ) f o r t ( s k 1 , t k ) , k = 1 , 2 , , p , lim t t k ( t t k ) 1 q x ( t ) = ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) , k = 1 , 2 , , p
and initial conditions in integral form:
t 0 I t 1 q x ( t ) | t = t 0 = x 0 .
Then the mild solution to the IVP (25), (45), (46) will be:
Definition 7.
A function x : [ t 0 , T ] R n is called a mild solution of the IVP for NIRLFrDE (25), (45), (46) if it satisfies the following Volterra integral–algebraic equations:
x ( t ) = x 0 Γ ( q ) ( t t 0 ) q 1 + 1 Γ ( q ) t 0 t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t 0 , s 0 ] ϕ k ( t , x ( t ) , x ( s k 1 0 ) ) , t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p , ϕ k ( t k , x ( t k ) , x ( s k 1 0 ) ) ( t t k ) q 1 + 1 Γ ( q ) t k t ( t s ) q 1 f ( s , x ( s ) ) d s , t ( t k , s k ] , k = 1 , 2 , , p .
Note Theorem 2 could be easily converted to mild solutions of the IVP for NIRLFrDE with mixed initial and impulsive conditions.

3.3. Examples

To illustrate the application of the above formulas for the mild solutions of non-instantaneous impulsive RL fractional differential equations and the types of impulsive functions we consider the following examples in the scalar case.
Example 1. Let t k = 2 k , s k = 2 k + 1 , k = 0 , 1 , 2 , , p , T = 2 p + 1 , p is a natural number, and f ( t , x ) = 1 t 1.5 for t [ 0 , T ] .
Consider the IVP for the scalar NIRLFrDE (8)–(10) with t 0 = 0 . The function f ( t , x ) is defined on the whole interval [ 0 , T ] . However (16) cannot be applied because the integral 0 t ( t s ) q 1 1 s 1.5 d s is not convergent for t > 1.5 .
Now, consider the IVP for the scalar NIRLFrDE (25), (27). The application of formula (33) for the solution causes no problem since we use the integral 2 k t ( t s ) q 1 1 s 1.5 d s for t ( 2 k , 2 k + 1 ] , k = 0 , 1 , 2 , , p which is convergent. Example 2. Consider the partial case of (8)–(10) with n = 1 , t 0 = 0 , q = 0.8 and f ( t , x ) 1 for t [ 0 , T ] .
We will consider several types of impulsive functions ϕ k ( t , x , y ) .
Case 1.
Let ϕ k ( t , x , y ) = g k ( t ) for t [ s k 1 , t k ] , ( k = 1 , 2 , 3 , , p ) .
Case 1.1.
Consider (8)–(10) and (16) to obtain the solution:
x ( t ) = x 0 t 0 . 2 Γ ( 0 . 8 ) + 1 . 25 t 0 . 8 Γ ( 0 . 8 ) t ( 0 , s 0 ] g k ( t ) t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p t k t 0 . 2 g k ( t k ) 1 . 25 t k Γ ( 0 . 8 ) t 0 . 2 + 1 . 25 t 0 . 8 Γ ( 0 . 8 ) t ( t k , s k ] , k = 1 , 2 , , p .
The solution depends on the initial value x 0 only on the interval ( 0 , s 0 ] . Therefore, two solutions with different initial values x 0 x 0 ˜ will coincide for all t > s 0 .
Case 1.2.
Consider (25), (27) and apply (33) to obtain the solution:
x ( t ) = x 0 t 0 . 2 Γ ( 0 . 8 ) + 1 . 25 t 0 . 8 Γ ( 0 . 8 ) t ( 0 , s 0 ] g k ( t ) t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p g k ( t k ) Γ ( 0 . 8 ) ( t t k ) 0 . 2 + 1 . 25 ( t t k ) 0 . 8 Γ ( 0 . 8 ) t ( t k , s k ] , k = 1 , 2 , , p .
Similarly to Case 1.1 we obtain that two solutions with different initial values x 0 x 0 ˜ coincide for all t > s 0 .
Case 2.
Let ϕ k ( t , x , y ) = a k x + b k for t [ s k 1 , t k ] , ( k = 1 , 2 , 3 , , p ) where a k , b k are constants.
-
If a k = 1 , b k = 0 , then the impulsive condition (9) is reduced to the condition x ( t ) = x ( t ) for t [ s k 1 , t k ] , ( k = 1 , 2 , 3 , , p ) and obviously the IVP for NIRLFrDE (8)–(10), respectively (25), (27), will have an infinite number of solutions.
-
If a k = 1 , b k 0 , then the impulsive condition (9) is reduced to the condition x ( t ) = x ( t ) + b for t [ s k 1 , t k ] , ( k = 1 , 2 , 3 , , p ) which has no solution and the IVP for NIRLFrDE (8)–(10), respectively (25), (27), will have no solution.
-
If a k 1 , b k = 0 , then the impulsive condition (9) is reduced to the condition x ( t ) = a x ( t ) for t [ s k 1 , t k ] , ( k = 1 , 2 , 3 , , p ) which has only the zero solution, and therefore any solution of IVP for NIRLFrDE (8)–(10), respectively (25), (27) will be zero on ( s k 1 , t k ] , ( k = 1 , 2 , 3 , , p ) . In this case from (16) we obtain the solution for (8)–(10):
x ( t ) = x 0 t 0 . 2 Γ ( 0 . 8 ) + 1 . 25 t 0 . 8 Γ ( 0 . 8 ) t ( 0 , s 0 ] 0 t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p , 1 . 25 t k Γ ( 0 . 8 ) t 0 . 2 + 1 . 25 t 0 . 8 Γ ( 0 . 8 ) t ( t k , s k ] , k = 1 , 2 , , p .
and from (33) we obtain the solution for (25), (27):
x ( t ) = x 0 t 0 . 2 Γ ( 0 . 8 ) + 1 . 25 t 0 . 8 Γ ( 0 . 8 ) t ( 0 , s 0 ] 0 t ( s k 1 , t k ] , k = 1 , 2 , 3 , , p , 1 . 25 ( t t k ) 0 . 8 Γ ( 0 . 8 ) t ( t k , s k ] , k = 1 , 2 , , p .
-
If a k 1 , b k 0 , then the impulsive condition (9) is reduced to the condition x ( t ) = a x ( t ) + b for t [ s k 1 , t k ] , ( k = 1 , 2 , , p ) which has an unique solution and we can talk about uniqueness of the solution IVP for NIRLFrDE (8)–(10), respectively (25), (27).
Case 3.
Let ϕ k ( t , x , y ) = arctan ( x ) + cos ( x ) + y for t [ s k 1 , t k ] , ( k = 1 , 2 , 3 , , p ) . Then the impulsive condition (9) is reduced to the algebraic equation x = arctan ( x ) + cos ( x ) + y which could have more than one solution (for example if y = 1 , then there are 5 constant solutions), i.e., we do not have uniqueness for the IVP for NIRLFrDE (8)–(10), respectively (25), (27).
Example 3. Consider the partial case of (8)–(10) with t 0 = 0 , q = 0.5 , f ( t , x ) = a x for t [ 0 , T ] , x R , a = c o n s t 0 and ϕ k ( t , x , y ) = b k y , k = 1 , 2 , , p . In this case we can consider both IVP (8), (10) and IVP (25), (27).
Case 1. When we apply (16) we obtain the solution of (8), (10):
x ( t ) = x 0 t E 0 . 5 , 0 . 5 ( a t ) for t [ 0 , s 0 ] , b k x ( s k 1 ) for t ( s k 1 , t k ] , t k t b k x ( s k 1 ) a Γ ( 0 . 5 ) 0 t k x ( s ) t k s d s + a Γ ( 0 . 5 ) 0 t x ( s ) t s d s for t ( t k , s k ] .
Case 2. When we apply (33) and Proposition 1 we obtain the solution of (25), (27):
x ( t ) = x 0 E 0 . 5 , 0 . 5 ( a t ) t for t [ 0 , s 0 ] , x 0 j = 1 k b j E 0 . 5 , 0 . 5 ( a s j 1 t j 1 ) s j 1 t j 1 for t ( s k 1 , t k ] , x 0 E 0 . 5 , 0 . 5 ( a t t k ) t t k j = 1 k b j E 0 . 5 , 0 . 5 ( a s j 1 t j 1 ) s j 1 t j 1 for t [ t k , s k ] ,

4. Brief Overview of RL Fractional Equations with Instantaneous Impulses

From Remark 4 note the case of instantaneous impulses is a special case of non-instantaneous impulses.
Several authors set up and studied various types of instantaneous impulsive differential equations for RL fractional differential equations and here we will give a brief overview of the types of instantaneous impulsive conditions in RL fractional differential equations. We emphasize only the case of the fractional order q ( 0 , 1 ) . Also, we will assume that the points { t k } are the points of impulses.
In [18] Kosmatov studied the RL fractional differential equations of fractional order q with the following impulsive conditions:
0 R L D t β x ( t ) | t = t i + 0 0 R L D t q x ( t ) | t = t i 0 = J i ( x ( t i ) ) , i = 1 , 2 , , p ,
where β ( 0 , q ) . As is shown in [12] the impulse functions 0 R L D t β x ( t ) are singular at t = t i , i = 1 , 2 , , p , for β ( 0 , q ) and the impulsive conditions (48) are unsuitable (see Remarks 4.1 [12]). Note the existence of a mild solution for Riemann–Liouville fractional differential equations with fractional order q and impulsive conditions of the form (48) is studied in [10]. In [19] Zhao studied a nonlinear Riemann–Liouville fractional differential equation with instantaneous impulsive conditions of the type
Δ x ( t i ) = J i ( x ( t i ) ) , i = 1 , 2 , , p ,
where Δ x ( t i ) = x ( t i + 0 ) x ( t i 0 ) . Note similar impulsive conditions are used in [20,21]. As is shown in [12] the integral representation of the solution with impulsive condition (49) is not correct (see Remark 4.3 [12]). In [8] the initial value problems of nonlinear impulsive RL fractional differential equations are studied with instantaneous impulsive conditions of the type
lim t t i + 0 ( t t i ) 1 q x ( t ) = J i ( t i , x ( t i ) , 0 R L D t β x ( t ) | t = t i 0 ) , i = 1 , 2 , , p ,
with β ( 0 , q ) . Similar impulsive conditions are studied in [11,22] but written in the form:
Δ I t i 1 q x ( t i ) = J i ( t i , x ( t i ) ) , i = 1 , 2 , , p ,
where Δ I t i 1 q x ( t i ) = I t i + 1 q x ( t i + ) I t i 1 q x ( t i ) = lim t t i + 0 ( t t i ) 1 q x ( t ) lim t t i ( t t i ) 1 q x ( t ) . In [23] the instantaneous impulsive conditions are given by
0 + R L I t 1 q x ( t i + 0 ) 0 + R L I t 1 q x ( t i 0 ) = J i ( x ( t i ) ) , i = 1 , 2 , , p .
The integral representation of the solution of instantaneous impulsive RL fractional derivative of order q ( 0 , 1 ) is given in [8] in the case when the impulsive conditions are
lim t t i + 0 ( t t i ) 1 + β q x ( t ) = x i , i = 1 , 2 , , p ,
and the initial condition is
0 R L I t 1 q x ( t ) | t = 0 = x 0
where β ( 0 , q ) .

Author Contributions

Conceptualization, R.A., S.H. and D.O.; methodology, R.A., S.H. and D.O.; validation, R.A., S.H. and D.O.; formal analysis, R.A., S.H. and D.O.; writing—original draft preparation, R.A., S.H. and D.O.; writing—review and editing, 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. Heymans, N.; Podlubny, I. Physical interpretation of initial conditions for fractional differential equations with Riemann–Liouville fractional derivatives. Rheol. Acta 2006, 45, 765–771. [Google Scholar] [CrossRef]
  2. Debnath, L. Recent applications of fractional calculus to science and engineering. IJMMS 2003, 54, 3413–3442. [Google Scholar] [CrossRef]
  3. Magin, R.L. Fractional Calculus in Bioengineering; Begell House Publishers: Danbury, CT, USA, 2006. [Google Scholar]
  4. Sabatier, J.; Agrawal, O.P.; Machado, J.A.T. (Eds.) Advances in Fractional Calculus: Theoretical Developments and Applications in Physics and Engineering; Springer: Dordrecht, The Netherlands, 2007. [Google Scholar]
  5. Vazquez, L.; Velasco, M.P.; Usero, D.; Jimenez, S. Fractional calculus as a modeling framework. Monografias Matematicas Garcia de Galdean 2018, 41, 187–197. [Google Scholar]
  6. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  7. Agarwal, R.; Hristova, S.; O’Regan, D. Non-Instantaneous Impulses in Differential Equations; Springer International Publishing: Cham, Switzerland, 2017. [Google Scholar]
  8. Liu, Y. Solvability of Anti-periodic BVPs for Impulsive Fractional Differential Systems Involving Caputo and Riemann–Liouville Fractional Derivatives. IJNSNS 2018, 19, 125–152. [Google Scholar] [CrossRef]
  9. Wang, J.R.; Feckan, M. Non-Instantaneous Impulsive Differential Equations. Basic Theory And Computation; IOP Publishing Ltd.: Bristol, UK, 2018; ISBN 978-0-7503-1704-7. [Google Scholar]
  10. Anguraj, A.; Kanjanadevi, S.; Nieto, J.J. Mild solutions of Riemann–Liouville fractional differential equations with fractional impulses. Nonlinear Anal. Model. Control. 2017, 22, 753–764. [Google Scholar] [CrossRef]
  11. Chaudhary, R.; Pandey, D.N. Monotone iterative technique for impulsive Riemann–Liouville fractional differential equations. Filomat 2018, 32, 3381–3395. [Google Scholar] [CrossRef]
  12. Liu, Y. Solvability of BVPs for impulsive fractional differential equations involving the Riemann–Liouville fractional derivatives. Stud. Univ. Babes-Bolyai Math. 2018, 63, 79–108. [Google Scholar] [CrossRef]
  13. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  14. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications Of Fractional Differential Equations; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  15. Yan, Z.M.; Lu, F.X. Existence of a new class of impulsive Riemann–Liouville fractional partial neutral functional differential equations with infinite delay. IAENG Intern. J. Appl. Math. 2015, 45, 300–312. [Google Scholar]
  16. Bajlekova, E. Fractional Evolution Equations in Banach Spaces. Ph.D. Thesis, Eindhoven University of Technology, Eindhoven, The Netherlands, 2001. [Google Scholar]
  17. Lakshmikantham, V.; Bainov, D.D.; Simeonov, P.S. Theory of Impulsive Differential Equations; World Scientific: Singapore, 1989. [Google Scholar]
  18. Kosmatov, N. Initial value problems of fractional order with fractional impulsive conditions. Results. Math. 2013, 63, 1289–1310. [Google Scholar] [CrossRef]
  19. Zhao, K. Impulsive fractional differential equation higher order problems of the higher- order fractional differential equation with eigenvalue arguments. Adv. Differ. Equ. 2015, 382, 16. [Google Scholar]
  20. Liu, Y. Bifurcation techniques for a class of boundary value problems of fractional impulsive differential equations. J. Nonlinear Sci. Appl. 2015, 8, 340–353. [Google Scholar] [CrossRef]
  21. Zhao, J.; Meng, F. Some Properties for Solutions of Riemann–Liouville Fractional Differential Systems with a Delay. J. Funct. Spaces 2018. [Google Scholar] [CrossRef]
  22. Liu, S.; Wu, X.; Zhou, X.F.; Jiang, W. Asymptotical stability of Riemann–Liouville fractional nonlinear systems. Nonlinear Dyn. 2016. [Google Scholar] [CrossRef]
  23. Wang, H.; Lin, X. Anti-periodic BVP of fractional order with fractional impulsive conditions and variable parameter. J. Appl. Math. Comput. 2017, 53, 285–301. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Agarwal, R.; Hristova, S.; O’Regan, D. Basic Concepts of Riemann–Liouville Fractional Differential Equations with Non-Instantaneous Impulses. Symmetry 2019, 11, 614. https://doi.org/10.3390/sym11050614

AMA Style

Agarwal R, Hristova S, O’Regan D. Basic Concepts of Riemann–Liouville Fractional Differential Equations with Non-Instantaneous Impulses. Symmetry. 2019; 11(5):614. https://doi.org/10.3390/sym11050614

Chicago/Turabian Style

Agarwal, Ravi, Snezhana Hristova, and Donal O’Regan. 2019. "Basic Concepts of Riemann–Liouville Fractional Differential Equations with Non-Instantaneous Impulses" Symmetry 11, no. 5: 614. https://doi.org/10.3390/sym11050614

APA Style

Agarwal, R., Hristova, S., & O’Regan, D. (2019). Basic Concepts of Riemann–Liouville Fractional Differential Equations with Non-Instantaneous Impulses. Symmetry, 11(5), 614. https://doi.org/10.3390/sym11050614

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