Next Article in Journal
Improvement of Some Hayashi–Ostrowski Type Inequalities with Applications in a Probability Setting
Next Article in Special Issue
On Some Characterizations for Uniform Dichotomy of Evolution Operators in Banach Spaces
Previous Article in Journal
Multiple q-Integral and Records from Geometrically Distributed Sequences
Previous Article in Special Issue
Generalized Proportional Caputo Fractional Differential Equations with Delay and Practical Stability by the Razumikhin Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generalized Proportional Caputo Fractional Differential Equations with Noninstantaneous Impulses: Concepts, Integral Representations, and Ulam-Type Stability

1
Department of Mathematics, Texas A&M University-Kingsville, Kingsville, TX 78363, USA
2
Faculty of Mathematics and Informatics, University of Plovdiv Paisii Hilendarski, 4000 Plovdiv, Bulgaria
3
School of Mathematical and Statistical Sciences, National University of Ireland, H91 TK33 Galway, Ireland
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(13), 2315; https://doi.org/10.3390/math10132315
Submission received: 28 May 2022 / Revised: 26 June 2022 / Accepted: 27 June 2022 / Published: 1 July 2022

Abstract

:
The generalized proportional Caputo fractional derivative is a comparatively new type of derivative that is a generalization of the classical Caputo fractional derivative, and it gives more opportunities to adequately model complex phenomena in physics, chemistry, biology, etc. In this paper, the presence of noninstantaneous impulses in differential equations with generalized proportional Caputo fractional derivatives is discussed. Generalized proportional Caputo fractional derivatives with fixed lower limits at the initial time as well as generalized proportional Caputo fractional derivatives with changeable lower limits at each impulsive time are considered. The statements of the problems in both cases are set up and the integral representation of the solution of the defined problem in each case is presented. Ulam-type stability is also investigated and some examples are given illustrating these concepts.

1. Introduction

There are different real life processes and phenomena that are characterized by rapid changes in their state at certain times, and the duration of these changes is not negligibly short. The dynamic of such processes can be adequately modeled with the help of differential equations with noninstantaneous impulses, and examples of such processes can be found in physics, biology, population dynamics, ecology, pharmacokinetics, and other areas.
Impulsive differential equations have become important in recent years as appropriate mathematical models, and there has been a significant development in impulsive theory, especially in the area of impulsive differential equations of the integer order (see for instance, the monographs [1,2,3] and the references therein).
Fractional calculus is the theory of integrals and derivatives of arbitrary non-integer order, which unifies and generalizes the concepts of ordinary differentiation and integration; for more details on geometric and physical interpretations of fractional derivatives and for a general historical perspective, we refer the reader to the monographs [4,5,6,7] and the cited references therein. There are also several papers concerning impulses and fractional derivatives [8]. However, to the best of our knowledge, this is the first paper that discusses and studies noninstantaneous impulses in differential equations with generalized proportional Caputo fractional derivatives.
The main goal of this paper is to present some basic points when introducing impulses in generalized proportional Caputo fractional differential equations. We will present two approaches in the interpretation of impulses in generalized proportional Caputo fractional differential equations. We note that in the case of differential equations with fractional order derivatives, these two approaches differ and the presence of impulses has a deep influence on the applied fractional derivatives. Both approaches are compared and their advantages/disadvantages are illustrated with examples. Three types of Ulam stability are defined for these problems with noninstantaneous impulses, and some sufficient conditions are obtained. Examples are given illustrating these types of stabilities.

2. Preliminary Notes on Generalized Proportional Fractional Derivatives

We recall that the generalized proportional fractional integral and the generalized Caputo proportional fractional derivative of a function u : [ a , b ] R , ( b ) are defined respectively by (as long as all integrals are well defined, see [9,10])
( a I α , ρ u ) ( t ) = 1 ρ α Γ ( α ) a t e ρ 1 ρ t s t s α 1 u ( s ) d s , t ( a , b ] , α 0 , ρ ( 0 , 1 ] ,
and
( a C D α , ρ u ) ( t ) = ( a I 1 α , ρ ( D 1 , ρ u ) ) ( t ) = 1 ρ 1 α Γ ( 1 α ) a t e ρ 1 ρ t s t s α ( D 1 , ρ u ) ( s ) d s , for t ( a , b ] , α ( 0 , 1 ) , ρ ( 0 , 1 ] ,
where ( D 1 , ρ u ) ( t ) = ( D ρ u ) ( t ) = ( 1 ρ ) u ( t ) + ρ u ( t ) .
Remark 1.
If ρ = 1 , then the generalized Caputo proportional fractional derivative is reduced to the classical Caputo fractional derivative.
We introduce the following classes of functions
C a α , ρ [ a , b ] = { u : [ a , b ] R : ( a D α , ρ u ) ( t ) e x i s t s o n ( a , b ] } , I a α , ρ [ a , b ] = { u : [ a , b ] R : ( a I α , ρ u ) ( t ) e x i s t s o n ( a , b ] } .
Note that if u C a α , ρ [ a , b ] then D 1 , ρ u ( . ) I a 1 α , ρ [ a , b ] .
Remark 2.
Note that the generalized proportional Caputo fractional derivative could easily be generalized for u : [ a , b ] R n component-wise.
In the case u C a α , ρ [ a , b ] , a C D α , ρ u ( . ) I a α , ρ [ a , b ] , we have the following result:
Lemma 1
(Theorem 5.3 [10]). Let ρ ( 0 , 1 ] and α ( 0 , 1 ) . Then we have
( a I α , ρ ( a C D α , ρ u ) ) ( t ) = u ( t ) u ( a ) e ρ 1 ρ ( t a ) , t ( a , b ] .
In the case u I a α , ρ [ a , b ] , a I α , ρ u ( . ) C a α , ρ [ a , b ] , ρ ( 0 , 1 ] , we obtain the following result:
Corollary 1
([10]). Let α ( 0 , 1 ) . Then
( a C D α , ρ ( a I α , ρ u ) ) ( t ) = u ( t ) , t ( a , b ] .
Lemma 2
(Theorem 5.2 [10]). For ρ ( 0 , 1 ] and α ( 0 , 1 ) , we have
( a I α , ρ e ρ 1 ρ t ( t a ) β 1 ) ( τ ) = Γ ( β ) ρ α Γ ( β + α ) e ρ 1 ρ τ ( τ a ) β 1 + α , β > 0 .
Remark 3.
The generalized proportional Caputo fractional derivative of a constant is not zero for ρ ( 0 , 1 ) (compare with the Caputo fractional derivative of a constant).
Corollary 2
(Remark 3.2 [10]). For ρ ( 0 , 1 ] and α ( 0 , 1 ) , the equality
( a C D α , ρ e ρ 1 ρ ( . ) ) ( t ) = 0 f o r t > a
holds.
From Lemma 1, we have the following result for the initial value problem for the generalized proportional Caputo fractional differential equation
( a C D α , ρ u ) ( t ) = f ( t , u ( t ) ) , t ( a , b ] , u ( a ) = u 0 , α ( 0 , 1 ) , ρ ( 0 , 1 ] .
Remark 4.
In our paper we will assume that the solution u of (1) satisfies u C a α , ρ [ a , b ] and a C D α , ρ u ( . ) I a α , ρ [ a , b ] .
Lemma 3.
For ρ ( 0 , 1 ] and α ( 0 , 1 ) , the solution u ( t ) of (1) satisfies the integral equation
u ( t ) = u 0 e ρ 1 ρ ( t a ) + 1 ρ α Γ ( α ) a t e ρ 1 ρ t s t s α 1 f ( s , u ( s ) ) d s , t ( a , b ] .

3. Noninstantaneous Impulses in Generalized Proportional Caputo Fractional Differential Equations

The fractional derivative leads to considering and studying two main cases: the case of a fixed lower limit of the fractional derivative and the case of a changeable lower limit of the fractional derivative. In both cases, we will set up the appropriate initial value problems.
Denote by Z 0 the set of all non-negative integers, by Z + the set of all natural numbers, and by Z [ a , b ] the set of all integers k : a k b where a , b Z 0 .
We will consider the case of long lasting impulses when the time of action is not negligibly small.
Let both sequences { t k } k = 1 , { s k } k = 0 : 0 s k < t k + 1 s k + 1 , lim k t k = be given. Let t 0 [ 0 , s 0 ) k = 1 [ t k , s k ) be the given initial time. Without loss of generality, we can assume t 0 [ 0 , s 0 ) .
Then, the impulsive differential equations are expressed by the following:
-
differential equation (continuous part) given for t ( t k , s k ] , k Z 0 ;
-
impulsive part (jump condition) given for t ( s k , t k + 1 ] , k Z 0 .
Remark 5.
The intervals ( s k , t k + 1 ] , k Z 0 , are called intervals of noninstantaneous impulses (they are connected with the time of action of the impulses).
Remark 6.
In the paper, we will use the notation x ( s k 0 ) = lim t s k , t < s k x ( t ) = x ( s k ) and x ( s k + 0 ) = lim t s k , t > s k x ( t ) , k Z 0 .

3.1. Noninstantaneous Impulses in Ordinary Differential Equations

Initially, we will discuss the case of an ordinary differential equation with noninstantaneous impulses.
Let us consider the following problem for the nonlinear differential equation
x ( t ) = f ( t , x ( t ) ) f o r t ( t k , s k ] , k Z 0 ,
with impulsive conditions
x ( t ) = ψ k ( t , x ( s k ) ) f o r t ( s k , t k + 1 ] , k Z 0 ,
and initial condition
x ( t 0 ) = x 0 ,
where x 0 R n , f : k = 0 ( t k , s k ] × R n R n , ψ k : [ s k , t k ] × R n R n .
The problem (3)–(5) is equivalent to the following three equivalent integral representations (for more detailed explanations see the book [2]):
x ( t ) = ψ k 1 ( t k , x ( s k 1 ) ) + t 0 t f ( s , x ( s ) ) d s t 0 t k f ( s , x ( s ) ) d s , t ( t k , s k ] , k Z 0 ψ k ( t , x ( s k ) ) , for t ( s k , t k + 1 ] , k Z 0 .
x ( t ) = x 0 + i = 0 k 1 J i ( t i + 1 , x ( s i ) ) s i t i + 1 f ( s , x ( s ) ) d s + t 0 t f ( s , x ( s ) ) d s , t ( t k , s k ] , k Z 0 ψ k ( t , x ( s k ) ) , for t ( s k , t k + 1 ] , k Z 0 ,
and
x ( t ) = ψ k 1 ( t k , x ( s k 1 ) + t k t f ( s , x ( s ) ) d s , t ( t k , s k ] , k Z 0 ψ k ( t , x ( s k ; t 0 , x 0 ) ) , for t ( s k , t k + 1 ] , k Z 0 ,
where s 1 = t 0 , ψ 1 ( t 0 , x 0 ) x 0 , J k ( t , x ) = ψ k ( t , x ) x , x R n .
In the case s k = t k + 1 , Equations (3)–(5) are reduced to well-known and studied impulsive differential equations with impulses at point t k , k Z + , and Formula (7) is reduced to the known integral representation of the solution of impulsive equations
x ( t ) = x 0 + i = 1 k J i ( t i , x ( t i ) ) + t 0 t f ( s , x ( s ) ) d s , t ( t k , t k + 1 ] , k Z 0 .

3.2. Generalized Proportional Caputo Fractional Derivative with Fixed Lower Limit and Non-Instantaneous Impulses

Based on the idea of the noninstantaneous impulses in ordinary differential equations, we will consider the noninstantaneous impulses in generalized proportional Caputo fractional differential equations. In the case of any type of fractional derivative in differential equations, the lower limit of the fractional derivative is very important. We will consider two cases.
We will first consider the case of a fixed lower limit at the initial time.
Consider the initial value problem (IVP) for the nonlinear differential equation with generalized proportional Caputo fractional derivative and noninstantaneous impulses (NIFrDE)
( t 0 C D α , ρ x ) ( t ) = f ( t , x ( t ) ) f o r t ( t k , s k ] , k Z 0 ,
with impulsive conditions
x ( t ) = ψ k ( t , x ( s k ) ) f o r t ( s k , t k + 1 ] , k Z 0 ,
and initial condition
x ( t 0 ) = x 0 ,
where x 0 R n , f : [ t 0 , ) × R n R n , ψ k : [ s k , t k + 1 ] × R n R n , ( k Z 0 ) .
Denote by PC [ t 0 , T ] = u : [ a , b ] R n : u C ( [ t 0 , T ] / { s k } k = 0 m 1 ) , where T .
We introduce the following assumptions:
S1.
The function f C ( [ t 0 , ) × R n ) and for any function u PC [ t 0 , ) , the function f ( . , u ( . ) ) I t 0 α , ρ [ t 0 , ) ;
S2.
The functions ψ k C ( [ s k , t k + 1 ] × R n ) , k Z 0 .
Remark 7.
Note that when Assumption (S2) is satisfied, then the equality x ( s k + 0 ) = ψ k ( s k , x ( s k ) ) follows from the equality x ( t ) = ψ k ( t , x ( s k ) ) for t ( s k , t k + 1 ] .
Following the integral representation (6) in the case of the ordinary derivative, we will define a mild solution.
Definition 1.
The function x ( t ) = x ( t ; t 0 , x 0 ) is called a mild solution of the IVP for NIFrDE (9)–(11) if it satisfies the integral equation
x ( t ) = ψ k 1 ( t k , x ( s k 1 ) ) e ρ 1 ρ ( t t k ) + 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s 1 ρ α Γ ( α ) t 0 t k e ρ 1 ρ t k s t k s α 1 f ( s , x ( s ) ) d s e ρ 1 ρ ( t t k ) t ( t k , s k ] , k Z 0 ψ k ( t , x ( s k ) ) for t ( s k , t k + 1 ] , k Z 0 ,
where s 1 = t 0 , ψ 1 ( t 0 , x 0 ) x 0 .
Remark 8.
In (12) for t ( t k , s k ] , we could also write equivalently
t 0 t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s = i = 0 k 1 t i s i e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s + i = 0 k 1 s i t i + 1 e ρ 1 ρ t s t s α 1 f ( s , ψ i ( s , x ( s i ) ) ) d s + t k t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s .
We will now give an equivalent integral presentation of the mild solution similar to (7) in the case of ordinary derivatives.
Lemma 4.
The mild solution x ( t ) of the IVP for NIFrDE (9)–(11) defined in Definition 1 satisfies the integral equation (and also conversely)
x ( t ) = x 0 e ρ 1 ρ i = 0 k 1 ( s i t i ) + ( t t k ) + 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s 1 ρ α Γ ( α ) t 0 t k e ρ 1 ρ t k s t k s α 1 f ( s , x ( s ) ) d s e ρ 1 ρ ( t t k ) + i = 0 k 1 ( J i ( t i + 1 , x ( s i ) ) + 1 ρ α Γ ( α ) t 0 s i e ρ 1 ρ s i s s i s α 1 f ( s , x ( s ) ) d s ) e ρ 1 ρ j = i + 1 k 1 ( s j t j ) + ( t t k ) t ( t k , s k ] , k Z 0 ψ k ( t , x ( s k ) ) , for t ( s k , t k + 1 ] , k Z 0 ,
where J k ( t , x ) = ψ k ( t , x ) x , k Z 0 , x R n .
The equivalence of (12) and (13) follows an induction argument.
Remark 9.
In the case ρ = 1 , the generalized proportional Caputo fractional derivative is reduced to the Caputo fractional derivative and the integral representation (12) is reduced to the one given in [11].
The connection between the mild solution defined by Definition 1 and the solution of the IVP for NIFrDE (9)–(11) is given in the following theorem:
Theorem 1.
Let Assumptions (S1) and (S2) hold. Then the mild solution of the IVP for NIFrDE (9)–(11) is a solution of the same problem.
Proof. 
It is easy to see that x ( t 0 ; t 0 , x 0 ) = x 0 .
For an arbitrary t ( t k , s k ] , k Z 0 , we take the generalized proportional Caputo fractional derivative of both sides of the first equation in (12), apply Corollary 1 and Corollary 2 with a = t 0 , and u ( t ) = f ( t , x ( t ; t 0 , x 0 ) ) , and obtain
( t 0 C D α , ρ x ) ( t ) = ψ k 1 ( t k , x ( s k 1 ) ) e ρ 1 ρ ( t 0 t k ) ( t 0 C D α , ρ e ρ 1 ρ ( t t 0 ) ) + ( t 0 C D α , ρ 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s ) 1 ρ α Γ ( α ) t 0 t k e ρ 1 ρ t k s t k s α 1 f ( s , x ( s ) ) d s e ρ 1 ρ ( t 0 t k ) ( t 0 C D α , ρ e ρ 1 ρ ( t t 0 ) ) = f ( t , x ( t ) ) .
For any t ( s k , t k + 1 ] , k Z 0 , it is clear that x ( t ) = ψ k ( t , x ( s k ) ) . □
Thus, from the above, the mild solution x ( . ) of the IVP for NIFrDE (9)–(11), defined by Equality (12), has the following properties:
P1.
In (12), the equality x ( t k ) = ψ k 1 ( t k , x ( s k 1 ) ) = lim t t k , t > t k x ( t ) , k Z + holds, i.e., the mild solution is continuous at t k and discontinuous at any t = s k , k Z 0 ;
P2.
In the case when there are no impulses, i.e., the points t k , s k do not have an influence on fractional Equation (9) so they do not exist here or t k = s k = t 0 , k Z 0 , then Equation (9) is defined for all t > t 0 and integral Equation (12) is reduced to integral Equation (2) with a = t 0 , b = ;
P3.
The mild solution has an equivalent form (13);
P4.
The mild solution is a solution of the IVP for NIFrDE (9)–(11).
P5.
In the case ρ = 1 , i.e., the case of the Caputo fractional derivative, Formula (12) is reduced to (2.35) [11].
Remark 10.
Note the mild solution of the IVP for NIFrDE (9)–(11) could be defined differently than in (12), but we choose this definition because of Properties 1–5.
Remark 11.
We will now discuss solutions of the IVP for NIFrDE (9)–(11). For simplicity in writing, let k = 1 ; we describe the solution x ( t ) as follows:
( A )
For t [ t 0 , s 0 ] our x ( t ) will be a solution of the IVP
( t 0 C D α , ρ x ) ( t ) = f ( t , x ( t ) ) , t ( t 0 , s 0 ] , x ( t 0 ) = x 0 .
Note from Lemma 3 the function x ( t ) has the integral representation in (12) for t ( t 0 , s 0 ] ;
( B )
For t ( s 0 , t 1 ] our x ( t ) will be x ( t ) = ψ 0 ( t , x ( s 0 ) ) where x ( s 0 ) is given in (A);
( C )
For t ( t 1 , s 1 ] our x ( t ) will be x ( t ) = y ( t ) , i.e., x = y | ( t 1 , s 1 ] , where y ( s ) is the solution of the IVP
( t 0 C D α , ρ y ) ( s ) = F ( s , y ( s ) ) , s ( t 0 , s 1 ] , y ( t 0 ) = y 0 ,
where
y 0 = e ρ 1 ρ ( t 0 t 1 ) [ ψ 0 ( t 1 , x ( s 0 ) ) 1 ρ α Γ ( α ) t 0 s 0 e ρ 1 ρ t 1 s t 1 s α 1 f ( s , x ( s ) ) d s 1 ρ α Γ ( α ) s 0 t 1 e ρ 1 ρ t 1 s t 1 s α 1 f ( s , ψ 0 ( s , x ( s 0 ) ) d s ]
and
F ( t , y ( s ) ) = f ( s , x ( s ) ) , s ( t 0 , s 0 ] f ( s , ψ 0 ( s , x ( s 0 ) ) , s ( s 0 , t 1 ] , f ( s , y ( s ) ) , s ( t 1 , s 1 ] ;
here x ( s ) for s ( t 0 , s 0 ] is given in (A). From Lemma 3, note y ( s ) for s ( t 0 , s 1 ] has the integral representation
y ( s ) = y 0 e ρ 1 ρ ( s t 0 ) + 1 ρ α Γ ( α ) t 0 s e ρ 1 ρ s z s z α 1 F ( z , y ( z ) ) d z .
For t ( t 1 , s 1 ] , note y ( t ) has the representation
y ( t ) = y 0 e ρ 1 ρ ( t t 0 ) + 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 F ( s , y ( s ) ) d s = y 0 e ρ 1 ρ ( t t 0 ) + 1 ρ α Γ ( α ) t 0 s 0 e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s + 1 ρ α Γ ( α ) s 0 t 1 e ρ 1 ρ t s t s α 1 f ( s , ψ 0 ( s , x ( s 0 ) ) d s + 1 ρ α Γ ( α ) t 1 t e ρ 1 ρ t s t s α 1 f ( s , y ( s ) ) d s ,
and since y ( s ) = x ( s ) for s ( t 1 , s 1 ] , we have
x ( t ) = y 0 e ρ 1 ρ ( t t 0 ) + 1 ρ α Γ ( α ) t 0 s 0 e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s + 1 ρ α Γ ( α ) s 0 t 1 e ρ 1 ρ t s t s α 1 f ( s , ψ 0 ( s , x ( s 0 ) ) d s + 1 ρ α Γ ( α ) t 1 t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s = e ρ 1 ρ ( t t 1 ) [ ψ 0 ( t 1 , x ( s 0 ) ) 1 ρ α Γ ( α ) t 0 s 0 e ρ 1 ρ t 1 s t 1 s α 1 f ( s , x ( s ) ) d s 1 ρ α Γ ( α ) s 0 t 1 e ρ 1 ρ t 1 s t 1 s α 1 f ( s , ψ 0 ( s , x ( s 0 ) ) d s ] + 1 ρ α Γ ( α ) t 0 s 0 e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s + 1 ρ α Γ ( α ) s 0 t 1 e ρ 1 ρ t s t s α 1 f ( s , ψ 0 ( s , x ( s 0 ) ) d s + 1 ρ α Γ ( α ) t 1 t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s
and we have the integral representation in (12) for t ( t 1 , s 1 ] .
Summarizing so far, note that for t ( t 0 , s 0 ] , our x ( t ) is given in (A); for t ( s 0 , t 1 ] , our x ( t ) is given in (B); and for t ( t 1 , s 1 ] , our x ( t ) is given in (C). Note that our solution x ( t ) is such that x P C [ t 0 , s 1 ] since ψ 0 ( t 1 , x ( s 0 ) ) = x ( t 1 ) and
x ( t 1 + 0 ) = y 0 e ρ 1 ρ ( t 1 t 0 ) + 1 ρ α Γ ( α ) t 0 s 0 e ρ 1 ρ t 1 s t 1 s α 1 f ( s , x ( s ) ) d s + 1 ρ α Γ ( α ) s 0 t 1 e ρ 1 ρ t 1 s t 1 s α 1 f ( s , ψ 0 ( s , x ( s 0 ) ) ) d s = [ ψ 0 ( t 1 , x ( s 0 ) ) 1 ρ α Γ ( α ) t 0 s 0 e ρ 1 ρ t 1 s t 1 s α 1 f ( s , x ( s ) ) d s 1 ρ α Γ ( α ) s 0 t 1 e ρ 1 ρ t 1 s t 1 s α 1 f ( s , ψ 0 ( s , x ( s 0 ) ) ) d s ] + 1 ρ α Γ ( α ) t 0 s 0 e ρ 1 ρ t 1 s t 1 s α 1 f ( s , x ( s ) ) d s + 1 ρ α Γ ( α ) s 0 t 1 e ρ 1 ρ t 1 s t 1 s α 1 f ( s , ψ 0 ( s , x ( s 0 ) ) ) d s = ψ 0 ( t 1 , x ( s 0 ) ) .
Now, we will study the existence of the mild solution of the IVP for NIFrDE (9)–(11) on a finite interval using the Banach contraction principle; i.e., let T ( t 0 , ) be given and let both finite sequences { t k } k = 1 m , { s k } k = 0 m : 0 t 0 s k < t k + 1 s k + 1 < t m < s m = T , k Z [ 0 , m 2 ] be given.
We use the following assumptions, which are similar to (S1) and (S2) but restricted on the finite interval [ 0 , T ] :
A1. 
The function f C ( [ t 0 , T ] × R n ) and for any u PC [ t 0 , T ] , the function f ( . , u ( . ) ) I t 0 α , ρ [ t 0 , T ] ;
A2. 
The functions ψ k C ( [ s k , t k + 1 ] × R n ) , k Z [ 0 , m 1 ] .
For u PC [ t 0 , T ] , we introduce the norm u PC = sup t [ t 0 , T ] | | u ( t ) | | .
Theorem 2.
(Existence result). Let the following conditions be satisfied:
1 .
Assumption A1 is satisfied and there exists a constant L f > 0 , such that
| | f ( t , u ) f ( t , v ) | | L f | | u v | | , f o r t [ t 0 , T ] , a n d u , v R n ;
2 .
Assumption A2 is satisfied and there exist constants L ψ k > 0 such that
| | ψ k ( t , x ) ψ k ( t , y ) | | L ψ k | | x y | | , t [ s k , t k + 1 ] , x , y R n , k Z [ 0 , m 1 ] ;
3 .
The inequality
2 L f ( T t 0 ) α ρ α Γ ( 1 + α ) + max { L ψ i , i = 1 , 2 , , m } < 1
holds.
Then the IVP for NIFrDE (9)–(11) has a unique mild solution x PC [ t 0 , T ] .
Proof. 
Define the operator Ω : PC [ t 0 , T ] PC [ t 0 , T ] by
Ω ( x ( t ) ) = ψ k 1 ( t k , x ( s k 1 ) ) e ρ 1 ρ ( t t k ) + 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s 1 ρ α Γ ( α ) t 0 t k e ρ 1 ρ t k s t k s α 1 f ( s , x ( s ) ) d s e ρ 1 ρ ( t t k ) for t ( t k , s k ] , k Z [ 0 , m ] ψ k ( t , x ( s k ) ) , for t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] .
The fixed point of the operator Ω (if any) is a mild solution of the IVP for NIFrDE (9)–(11).
Let y , y * PC [ t 0 , T ] . For t [ t 0 , s 0 ] we obtain
| | Ω ( y ( t ) ) Ω ( y * ( t ) ) | | y ( t 0 ) y * ( t 0 ) e ρ 1 ρ ( t t 0 ) + 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 L f | | y ( s ) y * ( s ) | | d s L f ρ α Γ ( 1 + α ) ( s 0 t 0 ) α | | y y * | | PC = C 0 | | y y * | | PC .
Note from (14), we have that C 0 < 1 .
For t ( s 0 , t 1 ] from Condition 2 and (15) for t = s 0 , we obtain
| | Ω ( y ( t ) ) Ω ( y * ( t ) ) | | = ψ 0 ( t , y ( s 0 ) ) ψ 0 ( t , y * ( s 0 ) ) L ψ 0 C 0 | | y y * | | PC .
from (14) it follows L ψ 0 C 0 < 1 .
For t ( t 1 , s 1 ] , we apply (15) for t = s 0 and obtain
| | Ω ( y ( t ) ) Ω ( y * ( t ) ) | | L ψ 0 y ( s 0 ) y * ( s 0 ) e ρ 1 ρ ( t t 1 ) + 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 L f | | y ( s ) y * ( s ) | | d s , + 1 ρ α Γ ( α ) t 0 t 1 e ρ 1 ρ t 1 s t 1 s α 1 L f | | y ( s ) y * ( s ) | | d s e ρ 1 ρ t t 1 2 L f ( t 1 t 0 ) α ρ α Γ ( 1 + α ) + L ψ 0 C 0 | | y y * | | PC = C 1 | | y y * | | PC .
Note that from (14) and C 0 < 1 , we have C 1 < 1 . For t ( t 2 , s 2 ] , we apply (15) for t = s 1 and obtain
| | Ω ( y ( t ) ) Ω ( y * ( t ) ) | | L ψ 1 y ( s 1 ) y * ( s 1 ) e ρ 1 ρ ( t t 2 ) + 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 L f | | y ( s ) y * ( s ) | | d s , + 1 ρ α Γ ( α ) t 0 t 2 e ρ 1 ρ t 2 s t 2 s α 1 L f | | y ( s ) y * ( s ) | | d s e ρ 1 ρ t t 2 2 L f ( t 2 t 0 ) α ρ α Γ ( 1 + α ) + L ψ 1 C 1 | | y y * | | PC = C 2 | | y y * | | PC ,
with C 2 < 1 .
From Inequalities (16)–(19), it follows that the operator Ω ( . ) is a contraction and according to the Banach theorem, it has a fixed point. This proves the claim. □

3.3. Generalized Proportional Caputo Fractional Derivative with Changeable Lower Limit at Each End Point of the Action of an Impulse

Initially, we will start with an interpretation of the solution of the corresponding fractional differential equation. Let the point ( t , x ) be on the integral curve of the IVP for the given fractional differential equation. The point ( t , x ) starts its motion from the initial point ( t 0 , x 0 ) t 0 is the initial time, x 0 is the initial value—and it continues to move along the integral curve ( t , x ( t ) ) of the solution of the corresponding fractional differential equation with a lower limit of the fractional derivative at the point t 0 . At the time s 0 > t 0 (the beginning time of the action of the impulse), the point instantaneously moves from position ( s 0 , x ( s 0 ) ) to position ( s 0 , ψ 0 ( s 0 , x ( s 0 ) ) . Then, the point continues its motion on the integral curve of the solution determined by the impulsive condition up to the point t 1 (the end point of the action of the impulse). Then, on the interval ( t 1 , s 1 ] , the point continues its motion on the integral curve of the solution determined by the corresponding fractional differential equation with a new lower limit at the time t 1 and initial value ( t 1 , x ( t 1 ) ) up to the time s 1 > t 1 , at which the impulse start to act and the point ( t , x ) jumps from the position ( s 1 , x ( s 1 ) ) to position ( s 1 , ψ 1 ( s 1 , x ( s 1 ) ) ) and so on. Note that the solution has a discontinuity/jump at any point s k , k Z 0 , with the amount of jump equal to ψ k ( s k , x ( s k ) ) x ( s k ) , k Z 0 .
According to the above description, the IVP for the nonlinear differential equation with generalized proportional Caputo fractional derivative and non-instantaneous impulses (NIFrDE) could be written in the form
( t k C D α , ρ x ) ( t ) = f ( t , x ( t ) ) f o r t ( t k , s k ] , k Z 0 , x ( t ) = ψ k ( t , x ( s k ) ) f o r t ( s k , t k + 1 ] , k Z 0 , x ( t 0 ) = x 0 ,
where x 0 R n , f : k = 0 [ t k , s k ] × R n R n , ψ k : [ s k , t k + 1 ] × R n R n , ( k Z 0 ) .
We introduce the following assumption:
S3. 
The function f C ( i = 0 [ t k , s k ] × R n ) .
Based on the integral presentation (8) in the case of ordinary derivative, we will introduce a mild solution in the case of the fractional derivative.
Definition 2.
The function x ( t ) = x ( t ; t 0 , x 0 ) is called a mild solution of the IVP for NIFrDE (20) if it satisfies the integral equation
x ( t ) = ψ k 1 ( t k , x ( s k 1 ) ) e ρ 1 ρ ( t t k ) + 1 ρ α Γ ( α ) t k t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s , for t ( t k , s k ] , k Z 0 ψ k ( t , x ( s k ) ) , for t ( s k , t k + 1 ] , k Z 0 .
where s 1 = t 0 , ψ 1 ( t 0 , x 0 ) x 0 .
Remark 12.
Note that in the case of ordinary derivatives and integrals, the equality
t 0 t f ( s , x ( s ) ) d s t 0 t k f ( s , x ( s ) ) d s = t k t f ( s , x ( s ) ) d s
holds, and in the case of integer derivatives, both integral representations (6) and (8) are equivalent.
In the case of fractional derivatives and fractional integrals, note that
t 0 t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s t 0 t k e ρ 1 ρ t k s t k s α 1 f ( s , x ( s ) ) d s t k t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s
holds, and the integral representations (12) and (21) are not equivalent.
Remark 13.
Note that in the description below (9), the right hand side part of the fractional equation f ( t , x ) is defined for all t t 0 , including the impulsive intervals (in spite of the fact that the fractional equation does not have an influence on the impulsive intervals). In (20), the right hand side part of the fractional equation f ( t , x ) is defined only on the intervals without impulses.
We will prove the equivalence of (20) and (21).
Theorem 3.
The mild solution of the IVP for NIFrDE (20) defined by (21) satisfies the IVP for NIFrDE (20) and conversely as well.
Proof. 
Let x ( t ) = x ( t ; t 0 , x 0 ) be a solution of the IVP for NIFrDE (20). Then, for any k Z + and t ( t k , s k ] from Lemma 1 with a = t k and b = s k , we have
x ( t ) x ( t k + 0 ) e ρ 1 ρ ( t t k ) = x ( t ) + ψ k 1 ( t k , x ( s k 1 ) ) e ρ 1 ρ ( t t k ) = ( t k I α , ρ ( t k C D α , ρ x ) ) ( t ) = ( t k I α , ρ f ( . , x ( . ) ) ) ( t ) .
Now, (22) establishes (21).
Let x ( t ) = x ( t ; t 0 , x 0 ) be a mild solution of the IVP for NIFrDE (20) defined by (21).
It is easy to check that x ( t 0 ) = x 0 .
Let k Z 0 and t ( t k , s k ] . Take a generalized proportional Caputo fractional derivative with a lower limit at the point t k of both sides of the first and the second equation in (21) and obtain
t k C D α , ρ x ( t ) = ψ k ( s k , x ( s k ) ) t k C D α , ρ e ρ 1 ρ ( t t k ) + t k C D α , ρ ( t k I α , ρ f ( . , x ( . ) ) ) ( t ) .
From Corollary 1 and Corollary 2 with a = t k and u ( t ) = f ( t , x ( t ) ) , we obtain the first equation of (20).
It is clear that x ( t ) = ψ k ( s k , x ( s k ) ) , t ( s k , t k + 1 ] , k Z 0 . □
Thus, the mild solution x ( . ) of the IVP for NIFrDE (20), defined by equality (21) has the following properties:
P1.
In (21), the equality x ( t k ) = ψ k 1 ( t k , x ( s k 1 ) ) = lim t t k , t > t k x ( t ) , k Z + holds; i.e., the mild solution is continuous at t k and discontinuous at any t = s k , k Z 0 ;
P2.
In the case when there are no impulses, i.e., the points t k , s k do not have an influence on the fractional equation (9) so they do not exist here or t k = s k = t 0 , k Z 0 , then equation (20) is reduced to the equation defined for all t > t 0 with a lower limit at t 0 and the integral equation (21) is reduced to the integral equation (2) with a = t 0 , b = ;
P3.
The integral equation (21) is equivalent to the fractional differential equation (20);
P4.
In the case ρ = 1 , i.e., the case of Caputo fractional derivative, the formula (12) is reduced to (2.38) [11].
We will prove the existence of a mild solution on a finite time interval.
We introduce the following assumption:
A3.
The function f C ( k = 0 m [ t k , s k ] × R n ) , and for any function u C ( [ t k , s k ] ) , k Z [ 0 , m ] , the function f ( . , u ( . ) ) I t k α , ρ [ t k , s k ] .
Theorem 4.
(Existence result). Let the following conditions be satisfied:
1.
Assumption A3 is satisfied and there exists a constant L f > 0 , such that
| | f ( t , u ) f ( t , v ) | | L f | | u v | | f o r t k = 0 m [ t k , s k ] , a n d u , v R n ;
2.
Assumption A2 is satisfied and there exist constants L ψ k > 0 , k Z [ 0 , m 1 ] , such that
| | ψ k ( t , x ) ψ k ( t , y ) | | L ψ k | | x y | | , t [ s k , t k + 1 ] , x , y R n , k Z [ 0 , m 1 ] ;
3.
For k Z [ 0 , m ] , the inequalities
L f ( s k t k ) α ρ α Γ ( 1 + α ) + L ψ k 1 < 1
hold where L ψ 1 = 0 .
Then, the IVP for NIFrDE (20) has a unique mild solution x PC [ t 0 , T ] .
Proof. 
Define the operator Ω : PC [ t 0 , T ] PC [ t 0 , T ] by
Ω ( x ( t ) ) = x 0 e ρ 1 ρ ( t t 0 ) + 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s , for t ( t 0 , s 0 ] ψ k 1 ( t k , x ( s k 1 ) ) e ρ 1 ρ ( t t k ) + 1 ρ α Γ ( α ) t k t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) d s , for t ( t k , s k ] , k Z [ 1 , m ] ψ k ( t , x ( s k ) ) , for t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] .
The fixed point of the operator Ω (if any) is a mild solution of the IVP for NIFrDE (20).
Let y , y * PC [ t 0 , T ] . For t [ t 0 , s 0 ] , we obtain
| | Ω ( y ( t ) ) Ω ( y * ( t ) ) | | 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 L f | | y ( s ) y * ( s ) | | d s L f ρ α Γ ( 1 + α ) ( s 0 t 0 ) α | | y y * | | PC .
For t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] , we obtain
| | Ω ( y ( t ) ) Ω ( y * ( t ) ) | | L ψ k | | y y * | | PC .
For t ( t k , s k ] , k Z [ 1 , m ] , we obtain
| | Ω ( y ( t ) ) Ω ( y * ( t ) ) | | 1 ρ α Γ ( α ) t k t e ρ 1 ρ t s t s α 1 L f | | y ( s ) y * ( s ) | | d s + ψ k 1 ( t k , y ( s k 1 ) ) ψ k 1 ( t k , y * ( s k 1 ) ) e ρ 1 ρ ( t t k ) L f ρ α Γ ( 1 + α ) ( t t k ) α + L ψ k 1 e ρ 1 ρ t t k | | y y * | | PC L f ( s k t k ) α ρ α Γ ( 1 + α ) + L ψ k 1 | | y y * | | PC .
From Inequalities (25)–(27), it follows that | | Ω ( y ) Ω ( y * ) | | PC γ | | y y * | | PC , i.e., the operator Ω is a contraction. □
Remark 14.
In the case of a changeable lower limit in the fractional derivative, i.e., the IVP for NIFrDE (20), the right hand side function f ( t , x ) and Condition (23) are defined only on the subintervals k = 0 m [ t k , s k ] . In the case of a fixed lower limit at the initial time, in spite of the fact that we do not use the function f ( t , x ) on the intervals of impulses, this function as well as Condition (14) have to be defined on the whole interval [ t 0 , T ] .

4. Ulam-Type Stability of Non-Instantaneous Impulsive Generalized Proportional Caputo Fractional Differential Equations

In this section, we will consider a finite time interval; i.e., let T ( t 0 , ) be given and let both finite sequences { t k } k = 1 m , { s k } k = 0 m : 0 t 0 s k < t k + 1 s k + 1 < t m < s m = T , k Z [ 0 , m 2 ] be given.

4.1. Fixed Lower Limit at the Initial Time of the Generalized Proportional Caputo Fractional Derivative

Let ε > 0 and Φ C ( [ t 0 , T ] / { s k } ) , Φ ( t ) 0 for t [ t 0 , T ] , Φ ( . ) be nondecreasing. We consider the following inequalities:
y ( t ) ψ k 1 ( t k , y ( s k 1 ) ) e ρ 1 ρ ( t t k ) + 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 f ( s , y ( s ) ) d s 1 ρ α Γ ( α ) t 0 t k e ρ 1 ρ t k s t k s α 1 f ( s , y ( s ) ) d s e ρ 1 ρ ( t t k ) ε f o r t ( t k , s k ] , k Z [ 0 , m ] , y ( t ) ψ k ( t , y ( s k ) ) ε , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] ,
and
y ( t ) ψ k 1 ( t k , y ( s k 1 ) ) e ρ 1 ρ ( t t k ) + 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 f ( s , y ( s ) ) d s 1 ρ α Γ ( α ) t 0 t k e ρ 1 ρ t k s t k s α 1 f ( s , y ( s ) ) d s e ρ 1 ρ ( t t k ) Φ ( t ) f o r t ( t k , s k ] , k Z [ 0 , m ] , y ( t ) ψ k ( t , y ( s k ) ) Φ ( t ) , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] ,
and
y ( t ) ψ k 1 ( t k , y ( s k 1 ) ) e ρ 1 ρ ( t t k ) + 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 f ( s , y ( s ) ) d s 1 ρ α Γ ( α ) t 0 t k e ρ 1 ρ t k s t k s α 1 f ( s , y ( s ) ) d s e ρ 1 ρ ( t t k ) ε Φ ( t ) f o r t ( t k , s k ] , k Z [ 0 , m ] , y ( t ) ψ k ( t , y ( s k ) ) ε Φ ( t ) , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] ,
We will define Ulam-types of stability for a differential equation with the generalized proportional Caputo fractional derivative with the fixed lower limits (9) and (10).
Definition 3.
The problem (9)–(11) is Ulam–Hyers stable if there exists a real number C > 0 such that for each ε > 0 and for each solution y PC [ t 0 , T ] of the inequality (28), there exists a mild solution x PC [ t 0 , T ] of the problem (9)–(11) with x 0 = y ( t 0 ) such that
| | y ( t ) x ( t ) | | C ε , t [ t 0 , T ] .
Definition 4.
The problem (9)–(11) is Ulam–Hyers-Rassias stable w.r.t. Φ if there exists a positive real number C such that for each ε > 0 and for each solution y PC [ t 0 , T ] of the inequality (30), there exists a mild solution x PC [ t 0 , T ] of the problem (9)–(11) with x 0 = y ( t 0 ) such that
x ( t ) y ( t ) C ε Φ ( t ) , t [ t 0 , T ] .
Definition 5.
The problem (9)–(11) is generalized Ulam–Hyers-Rassias stable with respect to Φ if there exists a positive real number C such that for each solution y PC [ t 0 , T ] of the inequality (29), there exists a mild solution x PC of the problem (9)–(11) with x 0 = y ( t 0 ) such that
x ( t ) y ( t ) C Φ ( t ) , t [ t 0 , T ] .
Theorem 5.
(Stability results). Assume the conditions of Theorem 2 are satisfied.
( i )
Suppose for any ε > 0 , Inequalities (28) have at least one solution. Then, Problem (9)–(11) is Ulam–Hyers stable;
( i i )
Suppose for any ε > 0 , Inequalities (30) have at least one solution. Then, Problem (9)–(11) is Ulam–Hyers–Rassias stable with respect to Φ;
( i i i )
Suppose Inequalities (29) have at least one solution. Then, Problem (9)–(11) is generalized Ulam–Hyers–Rassias stable with respect to Φ.
Proof. 
(i). Let ε > 0 be an arbitrary number and y PC [ t 0 , T ] be a solution of Inequalities (28). According to Theorem 2, there exists a mild solution x ( t ) = x ( t ; t 0 , x 0 ) of the IVP for NIFDE (9)–(11) on [ t 0 , T ] with x 0 = y ( t 0 ) .
We use induction to prove Inequality (31) with
C = max i Z [ 0 , k 1 ] { C i , C i L ψ i } f o r t ( t 0 , t k ]
and
C 0 = E α L f ρ α Γ ( 1 + α ) ( s 0 t 0 ) α ,
C i = ( t i t 0 ) α ρ α Γ ( α + 1 ) max j Z [ 0 , i 1 ] { C j , 1 + L ψ j C j } + 1 + L ψ i 1 C i 1 E α L f ρ α Γ ( 1 + α ) ( s i t 0 ) α ,
i Z [ 1 , m ] .
The inequality (31) holds for t = t 0 . □
-
Let t ( t 0 , s 0 ] be an arbitrary fixed point. Then, we obtain
x ( t ) y ( t ) 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ ( t * s ) ( t s ) α 1 f ( s , x ( s ) ) f ( s , y ( s ) ) d s + y ( t ) y ( t 0 ) e ρ 1 ρ ( t t 0 ) 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , y ( s ) ) d s L f ρ α Γ ( α ) t 0 t ( t s ) α 1 x ( s ) y ( s ) d s + ε .
According to the fractional generalization of the Gronwall inequality (Corollary 2 [12]), we obtain
x ( t ) y ( t ) ε E α L f ρ α Γ ( 1 + α ) ( t t 0 ) α .
Therefore, inequality (31) holds for t [ t 0 , s 0 ] with
C = C 0 = E α L f ρ α Γ ( 1 + α ) ( s 0 t 0 ) α .
-
Let t ( s 0 , t 1 ] be a fixed point. Then, from (31) for t = s 0 , we obtain x ( t ) y ( t ) = | ψ 0 ( t , x ( s 0 ) ) ψ 0 ( t , y ( s 0 ) ) + y ( t ) ψ 0 ( t , y ( s 0 ) ) L ψ 0 x ( s 0 ) y ( s 0 ) + ε ε ( 1 + L ψ 0 C 0 ) , i.e., Inequality (31) holds for t [ s 0 , t 1 ] with C = 1 + C 0 L ψ 0 .
Therefore, inequality (31) holds for t ( t 0 , t 1 ] with C = K 0 = max { C 0 , 1 + L ψ 0 C 0 } .
-
Let t ( t 1 , s 1 ] be an arbitrary fixed point.
Then
x ( t ) y ( t ) ψ 0 ( t 1 , x ( s 0 ) ) e ρ 1 ρ ( t t 1 ) ψ 0 ( t 1 , y ( s 0 ) ) e ρ 1 ρ ( t t 1 ) + 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 f ( s , x ( s ) ) f ( s , y ( s ) ) d s + 1 ρ α Γ ( α ) t 0 t 1 e ρ 1 ρ t 1 s t 1 s α 1 f ( s , x ( s ) ) f ( s , y ( s ) ) d s e ρ 1 ρ ( t t 1 ) + y ( t ) ψ 0 ( t 1 , y ( s 0 ) ) e ρ 1 ρ ( t t 1 ) 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , y ( s ) ) d s + 1 ρ α Γ ( α ) t 0 t 1 e ρ 1 ρ t 1 s t 1 s α 1 f ( s , y ( s ) ) d s e ρ 1 ρ ( t t 1 ) L ψ 0 x ( s 0 ) ) y ( s 0 ) e ρ 1 ρ ( t t 1 ) + ε + 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 L f x ( s ) y ( s ) d s + 1 ρ α Γ ( α ) t 0 t 1 e ρ 1 ρ t 1 s t 1 s α 1 L f x ( s ) y ( s ) d s e ρ 1 ρ ( t t 1 ) L f ρ α Γ ( α ) t 0 t * e ρ 1 ρ ( t s ) ( t s ) α 1 x ( s ) y ( s ) d s + ( t 1 t 0 ) α ρ α Γ ( α + 1 ) K 0 + 1 + L ψ 0 C 0 ε .
According to the fractional generalization of the Gronwall inequality (Corollary 2 [12]), we obtain
x ( t ) y ( t ) ( t 1 t 0 ) α ρ α Γ ( α + 1 ) K 0 + 1 + L ψ 0 C 0 E α L f ρ α Γ ( 1 + α ) ( t t 0 ) α ε .
Therefore, inequality (31) holds for t [ t 0 , s 1 ] with C = max { K 0 , C 1 } where
C = C 1 = ( t 1 t 0 ) α ρ α Γ ( α + 1 ) K 0 + 1 + L ψ 0 C 0 E α L f ρ α Γ ( 1 + α ) ( s 1 t 0 ) α .
-
Let t ( s 1 , t 2 ] be a fixed point. Then, x ( t ) y ( t ) = ψ 1 ( t , x ( s 1 ) ) ψ 1 ( t , y ( s 1 ) ) + y ( t ) ψ 0 ( t , y ( s 0 ) ) L ψ 1 x ( s 1 ) y ( s 1 ) + y ( t ) ψ 0 ( t , y ( s 0 ) ) ε L ψ 1 C 1 + ε , i.e., Inequality (31) holds for t [ s 1 , t 2 ] with C = 1 + C 1 L ψ 1 and on [ t 0 , t 2 ] with C = K 1 = max { K 0 , C 1 , 1 + C 1 L ψ 1 } = max { C 0 , 1 + C 0 L ψ 0 , C 1 , 1 + C 1 L ψ 1 } .
Following this procedure, we prove Inequality (31) and claim (i).
(iii).
The proof is similar to that of (i) with replacing ε by Φ ( t ) and using the monotonicity property of Φ ( . ) ;
(ii).
The proof is similar to that in (i) with replacement of ε by ε Φ ( t ) .

4.2. Changeable Lower Limit at Each Ending Time of the Impulses of the Generalized Proportional Caputo Fractional Derivative

Let ε > 0 , Ψ k : [ s k , t k + 1 ] R + , k Z [ 0 , m 1 ] and the function Φ C α , ρ ( k = 0 m [ t k , s k ] ) be non-negative nondecreasing on each interval [ t k , s k ] , k Z [ 0 , m ] . We consider the following inequalities:
( t k C D α , ρ y ) ( t ) f ( t , y ( t ) ) ε f o r t ( t k , s k ] , k Z [ 0 , m ] , y ( t ) ψ k ( t , y ( s k ) ) ε , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] ,
and
( t k C D α , ρ y ) ( t ) f ( t , y ( t ) ) Φ ( t ) f o r t ( t k , s k ] , k Z [ 0 , m ] , y ( t ) ψ k ( t , y ( s k ) ) Ψ k ( t ) , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] ,
and
( t k C D α , ρ y ) ( t ) f ( t , y ( t ) ) ε Φ ( t ) f o r t ( t k , s k ] , k Z [ 0 , m ] , y ( t ) ψ k ( t , y ( s k ) ) ε Ψ k ( t ) , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] ,
We well define Ulam-types of stability for the IVP for NIFrDE (20) for k Z [ 0 , m ] .
Definition 6.
The problem (20) is Ulam–Hyers stable if there exists a real number C > 0 such that for each ε > 0 and for each solution y PC [ t 0 , T ] of the inequality (36), there exists a solution x PC [ t 0 , T ] of the problem (20) with x 0 = y ( t 0 ) such that
| | y ( t ) x ( t ) | | C ε , t [ t 0 , T ] .
Definition 7.
The problem (20) is Ulam–Hyers–Rassias stable w.r.t. Φ and Ψ k , k Z [ 0 , m 1 ] if there exist positive real numbers C k , k Z [ 0 , m ] and K k , k Z [ 0 , m 1 ] , such that for each ε > 0 and for each solution y PC [ t 0 , T ] of the inequality (38), there exists a solution x PC [ t 0 , T ] of the problem (20) with x 0 = y ( t 0 ) such that
x ( t ) y ( t ) ε ( C k Φ ( t ) + K k ) , t ( t k , s k ] , k Z [ 0 , m ] ε ( C k Φ ( t ) + K k + Ψ k ( t ) ) , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] .
Definition 8.
The problem (20) is generalized Ulam–Hyers–Rassias stable with respect to Φ and Ψ k , k Z [ 0 , m 1 ] , if there exist positive real numbers C k , k Z [ 0 , m ] and K k , k Z [ 0 , m 1 ] , such that for each solution y PC [ t 0 , T ] of the inequality (37), there exists a solution x PC of the problem (20) with x 0 = y ( t 0 ) such that
x ( t ) y ( t ) C k Φ ( t ) + K k , t ( t k , s k ] , k Z [ 0 , m ] C k Φ ( t ) + K k + Ψ k ( t ) , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] .
Remark 15.
If Assumptions A2, A3 are satisfied, then the function y PC [ t 0 , T ] is a solution of Inequality (36) if there exist a function G C α , ρ ( k = 0 m [ t k , s k ] ) : G ( t ) ε for t k = 0 m [ t k , s k ] and a sequence of functions G k C ( [ s k , t k + 1 ] ) : G k ( t ) ε for t [ s k , t k + 1 ] , k Z [ 0 , m 1 ] , such that
( i )
t k C D t q y ( t ) = f ( t , y ( t ) ) + G ( t ) f o r t ( t k , s k ] , k Z [ 0 , m ] ;
( i i )
y ( t ) = ψ k ( t , y ( s k ) ) + G k ( t ) , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] .
Remark 16.
If Assumptions A2, A3 are satisfied, then the function y PC [ t 0 , T ] is a solution of Inequality (37) if there exist a function G C α , ρ ( k = 0 m [ t k , s k ] ) : G ( t ) Φ ( t ) for t k = 0 m [ t k , s k ] and a sequence of functions G k C ( [ s k , t k + 1 ] ) : | | G k ( t ) | | Ψ k ( t ) for t [ s k , t k + 1 ] , k Z [ 0 , m 1 ] ,
( i )
t k C D t q y ( t ) = f ( t , y ( t ) ) + G ( t ) f o r t ( t k , s k ] , k Z [ 0 , m ] ;
( i i )
y ( t ) = ψ k ( t , y ( s k ) ) + G k ( t ) , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] .
Remark 17.
If assumptions A2, A3 are satisfied then the function y PC [ t 0 , T ] is a solution of the inequality (38) if there exist a function G C α , ρ ( k = 0 m [ t k , s k ] ) : G ( t ) ε Φ ( t ) for t k = 0 m [ t k , s k ] and a sequence of functions G k C ( [ s k , t k + 1 ] ) : | | G k ( t ) | | ε Ψ k ( t ) for t [ s k , t k + 1 ] , k Z [ 0 , m 1 ] ,
( i )
t k C D t q y ( t ) = f ( t , y ( t ) ) + G ( t ) f o r t ( t k , s k ] , k Z [ 0 , m ] ;
( i i )
y ( t ) = ψ k ( t , y ( s k ) ) + G k ( t ) , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] .
Lemma 5.
Let Assumptions A2, A3 be satisfied. If y PC [ t 0 , T ] is a solution of Inequalities (36), then it satisfies the following integral-algebraic inequalities:
y ( t ) ) y ( t 0 ) e ρ 1 ρ ( t t 0 ) 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 f ( s , y ( s ) ) d s ( s 0 t 0 ) α ρ α Γ ( α + 1 ) + e ρ 1 ρ ( t t 0 ) ε t [ t 0 , s 0 ] , y ( t ) ψ k 1 ( t k , y ( s k 1 ) ) e ρ 1 ρ ( t t k ) 1 ρ α Γ ( α ) t k t e ρ 1 ρ t s t s α 1 f ( s , y ( s ) ) d s ( s k t k ) α ρ α Γ ( α + 1 ) + e ρ 1 ρ ( t t k ) ε , t ( t k , s k ] , k Z [ 1 , m ] , y ( t ) ψ k ( t , y ( s k ) ) ε , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] .
Proof. 
For any k Z [ 0 , m ] and t ( t k , s k ] , applying Remark 15 and Definition 2, we have
y ( t ) = ψ k 1 t k , y ( s k 1 ) + G k 1 ( t k ) e ρ 1 ρ ( t t k ) + 1 ρ α Γ ( α ) t k t e ρ 1 ρ t s t s α 1 f ( s , y ( s ) ) + G ( s ) d s
with ψ 1 ( t 0 , x ) y ( t 0 ) , x R n and G 1 ( t 0 ) = 0 .
This proves the claim of Lemma 5. □
Similarly, we have the following results.
Lemma 6.
Let Assumptions A2, A3 be satisfied. If y PC [ t 0 , T ] is a solution of Inequalities (37), then it satisfies the following integral-algebraic inequalities:
y ( t ) y ( t 0 ) e ρ 1 ρ ( t t 0 ) 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 f ( s , y ( s ) ) d s ( s 0 t 0 ) α ρ α Γ ( α + 1 ) Φ ( t ) t [ t 0 , s 0 ] , y ( t ) ψ k 1 ( t k , y ( s k 1 ) ) e ρ 1 ρ ( t t k ) 1 ρ α Γ ( α ) t k t e ρ 1 ρ t s t s α 1 f ( s , y ( s ) ) d s ( s k t k ) α ρ α Γ ( α + 1 ) Φ ( t ) + Ψ k 1 ( t k ) e ρ 1 ρ ( t t k ) , t ( t k , s k ] , k Z [ 1 , m ] , y ( t ) ψ k ( t , y ( s k ) ) Ψ k ( t ) , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] .
Proof. 
The proof is similar to the one of Lemma 5 applying Remark 16, the monotonicity of the function Φ ( . ) and the inequalities t k t e ρ 1 ρ ( t s ) ( t s ) α 1 Φ ( s ) d s t k t e ρ 1 ρ ( t s ) ( t s ) α 1 d s Φ ( t ) ( t t k ) α α Φ ( t ) and G k 1 ( t k ) e ρ 1 ρ ( t t k ) Ψ k 1 ( t k ) e ρ 1 ρ ( t t k ) , k Z [ 1 , m ] . □
From Remark 17 and Definition 2, we have:
Lemma 7.
Let Assumptions A2, A3 be satisfied. If y PC [ t 0 , T ] is a solution of Inequalities (38), then it satisfies the following integral-algebraic inequalities:
y ( t ) y ( t 0 ) e ρ 1 ρ ( t t 0 ) 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ t s t s α 1 f ( s , y ( s ) ) d s ( s 0 t 0 ) α ρ α Γ ( α + 1 ) Φ ( t ) ε t [ t 0 , s 0 ] , y ( t ) ψ k 1 ( t k , y ( s k 1 ) ) e ρ 1 ρ ( t t k ) 1 ρ α Γ ( α ) t k t e ρ 1 ρ t s t s α 1 f ( s , y ( s ) ) d s ( s k t k ) α ρ α Γ ( α + 1 ) Φ ( t ) + Ψ k 1 ( t k ) e ρ 1 ρ ( t t k ) ε , t ( t k , s k ] , k Z [ 1 , m ] , y ( t ) ψ k ( t , y ( s k ) ) Ψ k ( t ) ε , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] .
Now, we will study Ulam-type stability of Problem (20) on the finite interval [ 0 , T ] .
Theorem 6.
(Stability results). Assume the conditions of Theorem 4 are satisfied.
( i )
Suppose for any ε > 0 , Inequalities (36) have at least one solution. Then, Problem (20) is Ulam–Hyers stable.
( i i )
Suppose for any ε > 0 , Inequalities (38) have at least one solution. Then, Problem (20) is Ulam–Hyers–Rassias stable with respect to Φ and Ψ k , k Z [ 0 , m 1 ] .
( i i i )
Suppose Inequalities (37) have at least one solution. Then, Problem (20) is generalized Ulam–Hyers–Rassias stable with respect to Φ and ψ k , k Z [ 0 , m 1 ] .
Proof. 
(i). Let ε > 0 be an arbitrary number and y PC [ t 0 , T ] be a solution of Inequalities (36). Therefore, the integral-algebraic inequalities (42) hold. According to Theorem 4, there exists a mild solution x ( t ) = x ( t ; t 0 , x 0 ) of the IVP for NIFDE (20) on [ t 0 , T ] with x 0 = y ( t 0 ) .
For any t [ t 0 , T ] , we define the function γ ( t ) = sup s [ t 0 , t ] | | x ( s ) y ( s ) | | .
We use induction to prove
γ ( t ) C ε , t [ t 0 , T ] ,
where
C = C k = ( s k t k ) α ρ α Γ ( α + 1 ) + 1 + L ψ k 1 C k 1 E α L f ρ α Γ ( 1 + α ) ( s k t k ) α t ( t k , s k ] , k Z [ 0 , m ] , max { C k , 1 + L ψ k C k } , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ]
with C 0 = ( s 0 t 0 ) α ρ α Γ ( α + 1 ) .
Inequality (45) holds for t = t 0 . □
-
Let t ( t 0 , s 0 ] be an arbitrary fixed point. Since y C ( [ t 0 , s 0 ] ) , there exists a point t * [ t 0 , t ] such that γ ( t ) = | | x ( t * ) y ( t * ) | | .
If t * = t 0 , then Inequality (45) holds.
If t * ( t 0 , t ] , then according to Definition 2 and Lemma 5, we obtain
x ( t * ) y ( t * ) 1 ρ α Γ ( α ) t 0 t * e ρ 1 ρ ( t * s ) ( t * s ) α 1 f ( s , x ( s ) ) f ( s , y ( s ) ) d s + y ( t * ) y ( t 0 ) e ρ 1 ρ ( t * t 0 ) 1 ρ α Γ ( α ) t 0 t * e ρ 1 ρ ( t * s ) ( t * s ) α 1 f ( s , y ( s ) ) d s L f ρ α Γ ( α ) t 0 t * ( t * s ) α 1 x ( s ) y ( s ) d s + ( s 0 t 0 ) α ρ α Γ ( α + 1 ) ε .
According to the fractional generalization of the Gronwall inequality (Corollary 2 [12]), we obtain
γ ( t ) = x ( t * ) y ( t * ) ( s 0 t 0 ) α ρ α Γ ( α + 1 ) ε E α L f ρ α Γ ( 1 + α ) ( t t 0 ) α .
Therefore, Inequality (45) holds for t [ t 0 , s 0 ] with
C = C 0 = ( s 0 t 0 ) α ρ α Γ ( α + 1 ) E α L f ρ α Γ ( 1 + α ) ( s 0 t 0 ) α .
with M 0 = ( s 0 t 0 ) α ρ α Γ ( α + 1 ) .
-
Let t ( s 0 , t 1 ] be a fixed point.
If there exists t * ( s 0 , t 1 ] such that γ ( t ) = x ( t * ) y ( t * ) , then γ ( t ) = ψ 0 ( t * , x ( s 0 ) ) ψ 0 ( t * , y ( s 0 ) ) G 0 ( t * ) L ψ 0 x ( s 0 ) y ( s 0 ) + ε ε 1 + L ψ 0 C 0 , i.e., Inequality (45) holds for t [ t 0 , t 1 ] with C = 1 + C 0 L ψ 0 .
If there exists t * [ t 0 , s 0 ] such that γ ( t ) = x ( t * ) y ( t * ) , then according to the above, (45) holds with C = C 0 .
Therefore, Inequality (45) holds for t ( s 0 , t 1 ] with
C = max { C 0 , 1 + L ψ 0 C 0 } .
-
Let t ( t 1 , s 1 ] be an arbitrary fixed point.
If there exist a point t * [ t 0 , t 1 ] such that γ ( t ) = x ( t * ) y ( t * ) , then Inequality (45) holds with C = max { C 0 , 1 + L ψ 0 C 0 } .
If there exist a point t * ( t 1 , s 1 ] such that γ ( t ) = x ( t * ) y ( t * ) , then
x ( t * ) y ( t * ) 1 ρ α Γ ( α ) t 1 t * e ρ 1 ρ ( t * s ) ( t * s ) α 1 f ( s , x ( s ) ) f ( s , y ( s ) ) d s | | + | | y ( t * ) ψ 0 ( t 1 , y ( s 0 ) ) e ρ 1 ρ ( t * t 1 ) 1 ρ α Γ ( α ) t 1 t * e ρ 1 ρ t * s t * s α 1 f ( s , y ( s ) ) d s + ψ 0 ( t 1 , x ( s 0 ) ) ψ 0 ( t 1 , y ( s 0 ) ) e ρ 1 ρ ( t * t 1 ) L f ρ α Γ ( α ) t 1 t * e ρ 1 ρ ( t * s ) ( t * s ) α 1 x ( s ) y ( s ) d s + ( s 1 t 1 ) α ρ α Γ ( α + 1 ) + e ρ 1 ρ ( t t 1 ) ε + L ψ 0 C 0 ε L f ρ α Γ ( α ) t 1 t * e ρ 1 ρ ( t * s ) ( t * s ) α 1 x ( s ) y ( s ) d s + ( s 1 t 1 ) α ρ α Γ ( α + 1 ) + 1 + L ψ 0 C 0 ε .
According to the fractional generalization of the Gronwall inequality (Corollary 2 [12]), we obtain
γ ( t ) = x ( t * ) y ( t * ) ( s 1 t 1 ) α ρ α Γ ( α + 1 ) + 1 + L ψ 0 C 0 E α L f ρ α Γ ( 1 + α ) ( t * t 1 ) α ε .
Therefore, Inequality (45) holds for t [ t 1 , s 1 ] with
C = C 1 = ( s 1 t 1 ) α ρ α Γ ( α + 1 ) + 1 + L ψ 0 C 0 E α L f ρ α Γ ( 1 + α ) ( s 1 t 1 ) α .
Following this procedure, we prove Inequality (45).
This proves claim (i).
(iii).
Let y PC [ t 0 , T ] be a solution of Inequality (37). Therefore, the integral-algebraic inequalities (43) hold.
According to Theorem 4, there exists a mild solution x ( t ) = x ( t ; t 0 , x 0 ) of the IVP for NIFDE (20) on [ t 0 , T ] with x 0 = y ( t 0 ) .
We use induction to prove the inequality
x ( t ) y ( t ) C k Φ ( t ) + K k , t ( t k , s k ] , k Z [ 0 , m ] C k L ψ k Φ ( t ) + L ψ k K k + Ψ k ( t ) , t ( s k , t k + 1 ] , k Z [ 0 , m 1 ] ,
where
C k = ( s k t 0 ) α ρ α Γ ( 1 + α ) + L ψ k 1 C k 1 E α L f ρ α Γ ( 1 + α ) ( s k t k ) α , k Z [ 1 , m ]
and
K k = L ψ k 1 K k 1 + Ψ k 1 ( t k ) E α L f ρ α Γ ( 1 + α ) ( s k t k ) α , k Z [ 1 , m ]
C 0 = ( s 0 t 0 ) α ρ α Γ ( 1 + α ) E α L f ρ α Γ ( 1 + α ) ( s 0 t 0 ) α , K 0 = 0 .
-
Let t ( t 0 , s 0 ] be an arbitrary fixed point. According to Lemma 6, we obtain
x ( t ) y ( t ) 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , x ( s ) ) f ( s , y ( s ) ) d s + y ( t ) y ( t 0 ) e ρ 1 ρ ( t t 0 ) 1 ρ α Γ ( α ) t 0 t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , y ( s ) ) d s L f ρ α Γ ( α ) t 0 t ( t s ) α 1 x ( s ) y ( s ) d s + ( s 0 t 0 ) α ρ α Γ ( 1 + α ) Φ ( t ) .
According to the fractional generalization of the Gronwall inequality (Corollary 2 [12]), we obtain
x ( t ) y ( t ) ( s 0 t 0 ) α ρ α Γ ( 1 + α ) Φ ( t ) E α L f ρ α Γ ( 1 + α ) ( t t 0 ) α C 0 Φ ( t ) .
Therefore, Inequality (48) holds for t [ t 0 , s 0 ] with
C = C 0 = ( s 0 t 0 ) α ρ α Γ ( 1 + α ) E α L f ρ α Γ ( 1 + α ) ( s 0 t 0 ) α .
-
Let t ( s 0 , t 1 ] be a fixed point. Then, applying the monotonicity of the function Φ ( . ) , we obtain x ( t ) y ( t ) = ψ 0 ( t , x ( s 0 ) ) ψ 0 ( t , y ( s 0 ) ) G 0 ( t ) L ψ 0 x ( s 0 ) y ( s 0 ) + Ψ 0 ( t ) L ψ 0 C 0 Φ ( s 0 ) + Ψ 0 ( t ) L ψ 0 C 0 Φ ( t ) + Ψ 0 ( t ) ; i.e., Inequality (48) holds for t ( s 0 , t 1 ] with C = C 0 L ψ 0 .
-
Let t ( t 1 , s 1 ] be an arbitrary fixed point.
Then, according to Lemma 6 and Inequality (50), for t = s 0 we obtain
x ( t ) y ( t ) 1 ρ α Γ ( α ) t 1 t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , x ( s ) ) f ( s , y ( s ) ) d s | | + | | y ( t ) ψ 0 ( t 1 , y ( s 0 ) ) e ρ 1 ρ ( t t 1 ) 1 ρ α Γ ( α ) t 1 t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , y ( s ) ) d s + ψ 0 ( t 1 , x ( s 0 ) ) ψ 0 ( t 1 , y ( s 0 ) ) e ρ 1 ρ ( t t 1 ) L f ρ α Γ ( α ) t 1 t e ρ 1 ρ ( t s ) ( t s ) α 1 x ( s ) y ( s ) d s + ( s 1 t 1 ) α ρ α Γ ( 1 + α ) Φ ( t ) + Ψ 0 ( t 1 ) + L ψ 0 C 0 Φ ( s 0 ) L f ρ α Γ ( α ) t 1 t e ρ 1 ρ ( t s ) ( t s ) α 1 x ( s ) y ( s ) d s + ( s 1 t 1 ) α ρ α Γ ( 1 + α ) + L ψ 0 C 0 Φ ( t ) + Ψ 0 ( t 1 ) .
According to the fractional generalization of the Gronwall inequality (Corollary 2 [12]), we obtain
x ( t ) y ( t ) ( s 1 t 1 ) α ρ α Γ ( 1 + α ) + L ψ 0 C 0 E α L f ρ α Γ ( 1 + α ) ( t t 1 ) α Φ ( t ) + E α L f ρ α Γ ( 1 + α ) ( t t 1 ) α Ψ 0 ( t 1 ) ( s 1 t 1 ) α ρ α Γ ( 1 + α ) + L ψ 0 C 0 Φ ( t ) + Ψ 0 ( t 1 ) E α L f ρ α Γ ( 1 + α ) ( s 1 t 1 ) α = C 1 Φ ( t ) + K 1 .
where
C 1 = ( s 1 t 1 ) α ρ α Γ ( 1 + α ) + L ψ 0 C 0 E α L f ρ α Γ ( 1 + α ) ( s 1 t 1 ) α ,
and
K 1 = Ψ 0 ( t 1 ) E α L f ρ α Γ ( 1 + α ) ( s 1 t 1 ) α .
-
Let t ( s 1 , t 2 ] be a fixed point.
Then, applying the monotonicity of the function Φ ( . ) and Inequality (51) for t = s 1 , we obtain x ( t ) y ( t ) = ψ 1 ( t , x ( s 1 ) ) ψ 1 ( t , y ( s 1 ) ) G 1 ( t ) L ψ 1 x ( s 1 ) y ( s 1 ) + Ψ 1 ( t ) C 1 L ψ 1 Φ ( t ) + K 1 L ψ 1 + Ψ 1 ( t ) .
-
Let t ( t 2 , s 2 ] be an arbitrary fixed point.
Then, according to Lemma 6 and Inequality (51) for t = s 1 , we obtain
x ( t ) y ( t ) 1 ρ α Γ ( α ) t 2 t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , x ( s ) ) f ( s , y ( s ) ) d s + | | y ( t ) ψ 1 ( t 2 , y ( s 1 ) ) e ρ 1 ρ ( t t 2 ) 1 ρ α Γ ( α ) t 2 t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , y ( s ) ) d s + ψ 1 ( t 2 , x ( s 1 ) ) ψ 1 ( t 2 , y ( s 1 ) ) e ρ 1 ρ ( t t 2 ) L f ρ α Γ ( α ) t 2 t e ρ 1 ρ ( t s ) ( t s ) α 1 x ( s ) y ( s ) d s + ( s 2 t 2 ) α ρ α Γ ( 1 + α ) Φ ( t ) + Ψ 1 ( t 2 ) + L ψ 1 ( C 1 Φ ( s 1 ) + K 1 ) L f ρ α Γ ( α ) t 2 t e ρ 1 ρ ( t s ) ( t s ) α 1 x ( s ) y ( s ) d s + ( s 2 t 2 ) α ρ α Γ ( 1 + α ) + L ψ 1 C 1 Φ ( t ) + L ψ 1 K 1 + Ψ 1 ( t 2 ) .
According to the fractional generalization of the Gronwall inequality (Corollary 2 [12]), we obtain
x ( t ) y ( t ) ( s 2 t 2 ) α ρ α Γ ( 1 + α ) + L ψ 1 C 1 E α L f ρ α Γ ( 1 + α ) ( t t 2 ) α Φ ( t ) + L ψ 1 K 1 + Ψ 1 ( t 2 ) E α L f ρ α Γ ( 1 + α ) ( t t 2 ) α C 2 Φ ( t ) + K 2
where
C 2 = ( s 2 t 2 ) α ρ α Γ ( 1 + α ) + L ψ 1 C 1 E α L f ρ α Γ ( 1 + α ) ( s 2 t 2 ) α
and
K 2 = L ψ 1 K 1 + Ψ 1 ( t 2 ) E α L f ρ α Γ ( 1 + α ) ( s 2 t 2 ) α .
Following this procedure, we prove Inequality (48). From Inequality (48), the claim (iii).
(ii).
The proof is similar to the one in (iii) with the application of Lemma 7 instead of Lemma 6.

5. Examples

We will give some examples to illustrate the main results in this paper: the integral representations and the Ulam-type stability.

5.1. Integral Representation of the Solution

We will illustrate the application of both types of fractional differential equations and their suggested integral representations.
Example 1.
Let t k = 2 k < s k = 2 k + 1 , k Z 0 .
Fixed lower limit of the generalized proportional Caputo fractional derivative.
Consider the scalar IVP for IFrDE
0 c D α , ρ x ( t ) = 0 f o r t ( 2 k , 2 k + 1 ] , k Z 0 , x ( t ) = b k t x ( 2 k + 1 ) f o r t ( 2 k + 1 , 2 k + 2 ] , k Z 0 , x ( 0 ) = x 0 ,
where b k R .
From the integral representation (12), we obtain the mild solution
x ( t ) = x 0 e ρ 1 ρ t , t ( 0 , 1 ] , 2 k b k 1 x ( 2 k 1 ) e ρ 1 ρ ( t 2 k ) , t ( 2 k , 2 k + 1 ] , k 1 b k t x ( 2 k + 1 ) for t ( 2 k + 1 , 2 k + 2 ] , k Z 0 .
Note that the mild solution is piecewise continuous only at points 2 k + 1 , k Z 0 .
Changeable lower limit of the generalized proportional Caputo fractional derivative.
Consider the scalar IVP for IFrDE
2 k c D α , ρ x ( t ) = 0 f o r t ( 2 k , 2 k + 1 ] , k Z 0 , x ( t ) = b k t x ( 2 k ) f o r t ( 2 k + 1 , 2 k + 2 ] k Z 0 , x ( 0 ) = x 0 ,
where b k R .
From the integral representation (21), we obtain
x ( t ) = x 0 e ρ 1 ρ t , t ( 0 , 1 ] , 2 k b k 1 x ( 2 k 1 ) e ρ 1 ρ ( t k ) , t ( 2 k , 2 k + 1 ] , k 1 b k t x ( 2 k + 1 ) for t ( 2 k + 1 , 2 k + 2 ] , k Z 0 .
Note that the solutions obtained by both types of lower limits of generalized proportional Caputo fractional derivatives coincide.
Example 2.
Let t k = 2 k < s k = 2 k + 1 , k Z 0 .
Fixed lower limit of the fractional derivative.
Consider the scalar differential equation with the generalized proportional Caputo fractional derivative with a fixed lower limit at the initial time t 0 = 0
0 c D q , ρ x ( t ) = e ρ 1 ρ t f o r t ( 2 k , 2 k + 1 ] , k Z 0 , x ( t ) = b k t x ( 2 k + 1 ) f o r t ( 2 k + 1 , 2 k + 2 ] k Z 0 x ( 0 ) = x 0 ,
where b k R .
According to Definition 1, the mild solution of the IVP for NIFrDE (56) is given by
x ( t ) = x 0 + t α ρ α Γ ( 1 + α ) e ρ 1 ρ t , t [ 0 , 1 ] , 2 k b k 1 x ( 2 k 1 ) e ρ 1 ρ ( t 2 k ) + t α ρ α Γ ( 1 + α ) e ρ 1 ρ t ( 2 k ) α ρ α Γ ( 1 + α ) e ρ 1 ρ t , t ( 2 k , 2 k + 1 ] , k 1 , b k t x ( 2 k + 1 ) for t ( 2 k + 1 , 2 k + 2 ] , k Z 0 ,
Changeable lower limit of the fractional derivative.
Consider the scalar differential equation with the generalized proportional Caputo fractional derivative with a changeable lower limit at any point t k = 2 k at which the impulse stops its action
2 k c D q , ρ x ( t ) = e ρ 1 ρ t f o r t ( 2 k , 2 k + 1 ] , k Z 0 , x ( t ) = b k t x ( 2 k + 1 ) f o r t ( 2 k + 1 , 2 k + 2 ] k Z 0 x ( 0 ) = x 0 .
According to Definition 2, the mild solution of the IVP for IFrDE (58) is given by
x ( t ; t 0 , x 0 ) = x 0 + t α ρ α Γ ( 1 + α ) e ρ 1 ρ t , for t ( 0 , 1 ] 2 k b k x ( 2 k 1 ) + ( t 2 k ) α ρ α Γ ( 1 + α ) e ρ 1 ρ t , t ( 2 k , 2 k + 1 ] , k 1 b k t x ( 2 k + 1 ) , for t ( 2 k + 1 , 2 k + 2 ] , k Z 0 .
From Equalities (57) and (59), it can be seen that the lower limit of the fractional derivative in the differential equation has a huge influence on the type of solution. Therefore, it is very important in studying fractional differential equations with non-instantaneous impulses to start with the correct statement of the problem and to use the correct integral representation of the solution.

5.2. Ulam-Type Stability

We will provide some examples about Ulam-type stability for both defined types of fractional differential equations.
A. Fixed lower limit of the fractional derivative.
Example 3.
Let t = 0 , s 0 = 1.5 , t 1 = 2.5 , s 1 = T = 4 , α = 0.1 , ρ = 0.9 . Consider the following system of nonlinear differential equations with generalized proportional Caputo fractional derivative and noninstantaneous impulses
( 0 C D 0.1 , 0.9 x ) ( t ) = 0.1 sin ( y ( t ) ) , ( 0 C D 0.1 , 0.9 y ) ( t ) = 0.1 cos ( x ( t ) ) , f o r t ( 0 , 1.5 ] ( 2.5 , 4 ]
with impulsive conditions
x ( t ) = 0.01 x ( s 0 0 ) , y ( t ) = 0.01 y ( s 0 0 ) f o r t ( 1.5 , 2.5 ] ,
where x , y : [ 0 , 4 ] R . For functions f = ( f 1 , f 2 ) , f 1 ( t , x , y ) = 0.1 sin ( y ) , f 2 ( t , x , y ) = 0.1 cos ( x ) and ψ = ( p s i 1 , ψ 2 ) : ψ 1 ( x , y ) = 0.01 x , ψ 2 ( x , y ) = 0.01 y we obtain L f = 0.1 and L ψ = 0.01 > 0 . Then 2 L f ( T t 0 ) α ρ α Γ ( 1 + α ) + max { L ψ i } = 0.2 ( 4 0.1 ) 0 . 9 0.1 Γ ( 1.1 ) + 0.01 0.25405 < 1 .
According to Theorem 2, the IVP for NIFrDE (60) and (61) has a unique mild solution ( x , y ) PC [ 0 , 4 ] .
Now, consider the fractional inequalities with impulses
| u ( t ) 0.01 u ( 1.5 ) e ρ 1 ρ ( t t k ) + 0.1 0 . 9 0.1 Γ ( 0.1 ) 0 t e ρ 1 ρ t s t s 0.9 sin ( v ( s ) ) d s 0.1 0 . 9 0.1 Γ ( 0.1 ) 0 t k e ρ 1 ρ t k s t k s 0.9 sin ( v ( s ) ) d s e 1 9 ( t t k ) | 0.1 | v ( t ) 0.01 v ( 1.5 ) e ρ 1 ρ ( t t k ) + 0.1 0 . 9 0.1 Γ ( 0.1 ) 0 t e ρ 1 ρ t s t s 0.9 cos ( u ( s ) ) d s 0.1 0 . 9 0.1 Γ ( 0.1 ) 0 t k e ρ 1 ρ t k s t k s 0.9 cos ( u ( s ) ) d s e 1 9 ( t t k ) | 0.25 + t , f o r t ( t k , s k ] , k = 0 , 1 , | u ( t ) 0.01 u ( 1.5 ) | 0.1 , | v ( t ) 0.01 v ( 1.5 ) | 0.25 + t f o r t ( 1.5 , 2.5 ]
where u , v : [ 0 , 4 ] R .
The functions u ( t ) 0 and v ( t ) = t satisfy Inequalities (62) because | v ( t ) 0 + 0.1 t 0.1 0 . 9 0.1 Γ ( 1.1 ) 0.25 + t , t [ 0 , 1.5 ] , | t 0.01 ( 1.5 ) | 0.25 + t , t ( 1 , 5 , 2.5 ] and | t 0.01 ( 1.5 ) e 1 9 ( t 2.5 ) | + 0.1 t 0.1 0 . 9 0.1 Γ ( 1.1 ) + ( 0.1 ) 2 . 5 0.1 0 . 9 0.1 Γ ( 1.1 ) 0.25 + t , t ( 2.5 , 4 ] (see Figure 1).
Therefore, Inequalities (29) are satisfied with Φ ( t ) = ( 0.1 , 0.25 + t ) . According to Theorem 5, Problem (60) and (61) is generalized Ulam–Hyers–Rassias stable w.r.t. to Φ ( t ) .
B. Changeable lower limit of the fractional derivative.
Example 4.
Let t = 0 , s 0 = 1.5 , t 1 = 2.5 , s 2 = T = 4 , α = 0.1 , ρ = 0.9 . Consider the following system of nonlinear differential equations with generalized proportional Caputo fractional derivative and noninstantaneous impulses
( 0 C D 0.1 , 0.9 x ) ( t ) = 0.1 sin ( y ( t ) ) , ( 0 C D 0.1 , 0.9 y ) ( t ) = 0.1 cos ( x ( t ) ) , f o r t ( 0 , 1.5 ] ( 2.5 C D 0.1 , 0.9 x ) ( t ) = 0.1 sin ( y ( t ) ) , ( 2.5 C D 0.1 , 0.9 y ) ( t ) = 0.1 cos ( x ( t ) ) , f o r t ( 2.5 , 4 ]
with impulsive conditions
x ( t ) = 0.01 x ( s 0 0 ) , y ( t ) = 0.01 y ( s 0 0 ) f o r t ( 1.5 , 2.5 ] .
In this case, L f = 0.1 and L ψ = 0.01 .
Thus, the inequality 0.1 ( 1.5 0 ) 0.1 0 . 9 0.1 Γ ( 1.1 ) 0.100121 < 1 holds.
Then according to Theorem 4, the IVP for NIFrDE (63) and (64) has a unique mild solution ( x , y ) PC [ 0 , 4 ] .
Now, consider the fractional inequalities with impulses
| ( 0 C D 0.1 , 0.9 u ) ( t ) 0.1 sin 2 ( v ( t ) ) | 0.05 , | ( 0 C D 0.1 , 0.9 v ) ( t ) 0.1 cos ( u ( t ) ) | 0.05 , f o r t ( 0 , 1.5 ] , | ( 2.5 C D 0.1 , 0.9 u ) ( t ) 0.1 sin 2 ( v ( t ) ) | 0.05 , | ( 2.5 C D 0.1 , 0.9 v ) ( t ) 0.1 cos ( u ( t ) ) | 0.05 , f o r ( 2.5 , 4 ] , | u ( t ) 0.01 u ( s 0 0 ) | 0.85 , | v ( t ) 0.01 v ( s 0 0 ) | 0.85 f o r t ( 1.5 , 2.5 ]
The functions
u ( t ) = e 1 9 t , for t ( 0 , 1.5 ] 0.01 e 1 9 1.5 , for t ( 1.5 , 2.5 ] e 1 9 ( t 2.5 ) , for t ( 2.5 , 4 ] , v ( t ) = 0.05 e 1 9 t t 0.2 , for t ( 0 , 1.5 ] 0.0005 e 1 9 1.5 1 . 5 0.2 , for t ( 1.5 , 2.5 ] 0.05 e 1 9 ( t 2.5 ) ( t 2.5 ) 0.2 , for t ( 2.5 , 4 ] ,
satisfy Inequalities (65). Therefore, Inequalities (37) are satisfied with Φ ( t ) 0.05 and Ψ ( t ) 0.85 . According to Theorem 6, Problem (63) and (64) is generalized Ulam–Hyers–Rassias stable w.r.t. to 0.05 and 0.85 .

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.; writing—original draft preparation, R.A., S.H. and D.O.; writing—review and editing, R.A., S.H. and D.O.; supervision, R.A., S.H., and D.O.; funding acquisition, S.H. All authors have read and agreed to the published version of the manuscript.

Funding

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

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Benchohra, M.; Henderson, J.; Ntouyas, S.K. Impulsive Differential Equations and Inclusions; Hindawi Publishing Corporation: New York, NY, USA, 2006; Volume 2. [Google Scholar]
  2. Lakshmikantham, V.; Bainov, D.D.; Simeonov, P.S. Theory of Impulsive Differential Equations; World Scientific: Singapore, 1989. [Google Scholar]
  3. Samoilenko, A.M.; Perestyuk, N.A. Impulsive Differential Equations; World Scientific Series on Nonlinear Science, Series A: Monographs and Treatises, 14; World Scientific Publishing Co., Inc.: River Edge, NJ, USA, 1995. (In Russian) [Google Scholar]
  4. Das, S. Functional Fractional Calculus; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  5. Diethelm, K. The Analysis of Fractional Differential Equations; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  6. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  7. Samko, G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach: Philadelphia, PA, USA, 1993. [Google Scholar]
  8. Hristova, S.; Abbas, M.I. Explicit solutions of initial value problems for fractional generalized proportional differential equations with and without impulses. Symmetry 2021, 13, 996. [Google Scholar] [CrossRef]
  9. Jarad, F.; Abdeljawad, T. Generalized fractional derivatives and Laplace transform. Discret. Contin. Dyn. Syst. Ser. S 2020, 13, 709–722. [Google Scholar] [CrossRef] [Green Version]
  10. Jarad, F.; Abdeljawad, T.; Alzabut, J. Generalized fractional derivatives generated by a class of local proportional derivatives. Eur. Phys. J. Spec. Top. 2017, 226, 3457–3471. [Google Scholar] [CrossRef]
  11. Agarwal, R.P.; Hristova, S.; O’Regan, D. Non-Instantaneous Impulses in Differential Equations; Springer: New York, NY, USA, 2017. [Google Scholar]
  12. Ye, N.; Gao, J.; Ding, Y. Gronwall inequality and its application to fractional differential equation. L. Math. Anal. Appl. 2007, 328, 1075–1081. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Graph of the bound 0.25 + t and the difference of the function v ( t ) = t and its integral presentation on [ 0 , 4 ] .
Figure 1. Graph of the bound 0.25 + t and the difference of the function v ( t ) = t and its integral presentation on [ 0 , 4 ] .
Mathematics 10 02315 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Agarwal, R.; Hristova, S.; O’Regan, D. Generalized Proportional Caputo Fractional Differential Equations with Noninstantaneous Impulses: Concepts, Integral Representations, and Ulam-Type Stability. Mathematics 2022, 10, 2315. https://doi.org/10.3390/math10132315

AMA Style

Agarwal R, Hristova S, O’Regan D. Generalized Proportional Caputo Fractional Differential Equations with Noninstantaneous Impulses: Concepts, Integral Representations, and Ulam-Type Stability. Mathematics. 2022; 10(13):2315. https://doi.org/10.3390/math10132315

Chicago/Turabian Style

Agarwal, Ravi, Snezhana Hristova, and Donal O’Regan. 2022. "Generalized Proportional Caputo Fractional Differential Equations with Noninstantaneous Impulses: Concepts, Integral Representations, and Ulam-Type Stability" Mathematics 10, no. 13: 2315. https://doi.org/10.3390/math10132315

APA Style

Agarwal, R., Hristova, S., & O’Regan, D. (2022). Generalized Proportional Caputo Fractional Differential Equations with Noninstantaneous Impulses: Concepts, Integral Representations, and Ulam-Type Stability. Mathematics, 10(13), 2315. https://doi.org/10.3390/math10132315

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