Next Article in Journal
High-Dimensional Chaotic Lorenz System: Numerical Treatment Using Changhee Polynomials of the Appell Type
Next Article in Special Issue
Derivative-Free Conformable Iterative Methods for Solving Nonlinear Equations
Previous Article in Journal
A New Class of Generalized Fractal and Fractal-Fractional Derivatives with Non-Singular Kernels
Previous Article in Special Issue
Fractional Scale Calculus: Hadamard vs. Liouville
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Boundary Value Problem for Impulsive Delay Fractional Differential Equations with Several Generalized Proportional Caputo Fractional Derivatives

by
Ravi P. Agarwal
1,* and
Snezhana Hristova
2
1
Department of Mathematics, Texas A&M University-Kingsville, Kingsville, TX 78363, USA
2
Faculty of Mathematics and Informatics, Plovdiv University “P. Hilendarski”, 4000 Plovdiv, Bulgaria
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(5), 396; https://doi.org/10.3390/fractalfract7050396
Submission received: 27 April 2023 / Revised: 5 May 2023 / Accepted: 9 May 2023 / Published: 12 May 2023
(This article belongs to the Special Issue Feature Papers for the 'General Mathematics, Analysis' Section)

Abstract

:
A scalar nonlinear impulsive differential equation with a delay and generalized proportional Caputo fractional derivatives (IDGFDE) is investigated. The linear boundary value problem (BVP) for the given fractional differential equation is set up. The explicit form of the unique solution of BVP in the special linear case is obtained. This formula is a generalization of the explicit solution of the case without any delay as well as the case of Caputo fractional derivatives. Furthermore, this integral form of the solution is used to define a special proportional fractional integral operator applied to the determination of a mild solution of the studied BVP for IDGFDE. The relation between the defined mild solution and the solution of the BVP for the IDGFDE is discussed. The existence and uniqueness results for BVP for IDGFDE are proven. The obtained results in this paper are a generalization of several known results.

1. Introduction

Fractional differential equations play an important role in the modeling and development of many fields in science and engineering. In the literature, several different types of fractional derivatives are defined and applied (see, for example, Refs. [1,2] for Hilfer operators, Refs. [3,4,5] for derivatives depending on another function, and Refs. [6,7] for derivatives involving arbitrary kernels). In the last decade, the generalized proportional fractional derivative is defined [8,9] and studied (see, for example, for stability properties [10,11,12] and for stochastic fractional differential equations [13]).
In this paper, the boundary value problem (BVP) for nonlinear impulsive fractional differential equations with a delay and generalized proportional Caputo-type derivatives (GPCFD) is investigated. Initially, the linear case of a scalar impulsive differential equation with several GPCFDs is considered. The explicit proportional fractional integral form of the solution is obtained. As partial cases of the new formula, the BVP for a linear impulsive differential equation with one GPCFD with several Caputo fractional derivatives and one Caputo fractional derivative is considered and the integral form of the exact solutions is presented. These forms are generalizations of the solutions for the initial value problem of a linear fractional equation with and without impulses. Some of the partial cases coincide with the known results in the literature. The new integral form of the exact solution is applied to define a proportional fractional integral operator and the mild solution of the studied nonlinear problem. Existence and uniqueness results are provided. We use Banach’s fixed point theorem and the defined proportional fractional integral operator to prove the existence of a mild solution.
Note that the BVP for a couple of multi-term delay-generalized proportional Caputo fractional differential equations is studied in [14]. In this paper, impulses are involved in the differential equation. These impulses totally change the differential equation as well as the behavior of the solution. It is connected with the memory property of the fractional derivatives. We consider the case when the lower limit of the GPCFD is changed at any impulsive time (see [14]).

2. Basic Notes on Fractional Calculus

For a better understanding of the main results, we will provide the basic definitions and some known in the literature results which will be used in the proofs.
Let τ < T be given numbers and y : [ τ , T ] R .
The generalized proportional fractional integral (GPFI) is defined by (see [8,9])
( τ I α , ρ y ) ( t ) = 1 ρ α Γ ( α ) τ t e ρ 1 ρ t σ t σ α 1 y ( σ ) d σ , t ( τ , T ] , α 0 , ρ ( 0 , 1 ] ,
and the generalized Caputo proportional fractional derivative (GPCFD) is defined by
( τ C D α , ρ y ) ( t ) = 1 ρ ρ 1 α Γ ( 1 α ) τ t e ρ 1 ρ t s t σ α y ( σ ) d σ + ρ α Γ ( 1 α ) τ t e ρ 1 ρ t σ t σ α y ( σ ) d σ for t ( τ , T ] , α ( 0 , 1 ) , ρ ( 0 , 1 ] .
Remark 1.
The GPFI and the GPCFD are generalizations of the fractional integral τ I α y ( t ) and the Caputo fractional derivative τ C D α y ( t ) , respectively, ( ρ = 1 in (1)) (see, for example, [15]).
We will provide some known results from the literature about GPFI and GPCFD which will be used later.
Lemma 1 
(Proposition 3.7 [8]).  For ρ ( 0 , 1 ] , α , γ > 0 we have
τ I α , ρ e ρ 1 ρ . ( . τ ) γ 1 ( t ) = Γ ( γ ) ρ α Γ ( γ + α ) e ρ 1 ρ t ( t τ ) α + γ 1 , t > τ .
Corollary 1.
Let ρ ( 0 , 1 ] , α > 0 . Then τ I α , ρ e ρ 1 ρ . ( t ) = 1 ρ α Γ ( 1 + α ) e ρ 1 ρ t ( t τ ) α , t > τ , holds.
We will use the following partial case of proposition 5.2 [8].
Lemma 2.
For ρ ( 0 , 1 ] , α ( 0 , 1 ) , γ > 1 we have
τ C D α , ρ e ρ 1 ρ . ( . τ ) γ 1 ( t ) = ρ α Γ ( γ ) Γ ( γ α ) e ρ 1 ρ t ( t τ ) γ 1 α , t > τ ,
and
τ C D α , ρ e ρ 1 ρ . ( t ) = 0 , t > τ .
Lemma 3 
(Theorem 3.8 [8]).  Let ρ ( 0 , 1 ] , β , γ > 0 and ν C ( [ τ , T ] , R ) . Then
τ I β , ρ τ I γ , ρ ν ( t ) = ( τ I β + γ , ρ ν ) ( t ) , t ( τ , T ] .
Furthermore, consider the following partial case of Theorem 5.3 [8] for α ( 0 , 1 ) .
Lemma 4.
For ρ ( 0 , 1 ] , α ( 0 , 1 ) and ν C 1 [ τ , T ] , we have
τ I α , ρ τ C D α , ρ ν ( t ) = ν ( t ) ν ( τ ) e ρ 1 ρ ( t τ ) , t ( τ , T ] .
Corollary 2 
([8]).  Let α ( 0 , 1 ) , ρ ( 0 , 1 ] and ν C [ τ , T ] . Then ( τ C D α , ρ ( τ I α , ρ ν ) ) ( t ) = ν ( t ) for t ( τ , T ] .

3. Impulsive Linear Differential Equations with Caputo-Type Fractional Derivatives

Let the points 0 = t 0 < t 1 < t 2 < < t p < t p + 1 = T be given, and T < .
Let α ( 0 , 1 ) , ρ ( 0 , 1 ] . Consider the classes of functions
PC [ 0 , T ] = { y C ( k = 0 p ( t k , t k + 1 ] ] , R ) : lim t t k , t > t k y ( t ) = y ( t k + 0 ) < , lim t t k , t < t k y ( t ) = y ( t k ) } ; PC 1 [ 0 , T ] = y PC [ 0 , T ] : y C 1 ( ( t k , t k + 1 ] , R ) , k = 0 , 1 , 2 , , p , ; PC α , ρ [ 0 , T ] = y PC 1 [ 0 , T ] : ( t k D α , ρ y ) ( t ) exists   for t ( t k , t k + 1 ] , k = 0 , 1 , , p ; PI α , ρ [ 0 , T ] = y PC [ 0 , T ] : ( t k I α , ρ y ) ( t ) exists   for t ( t k , t k + 1 ] , k = 0 , 1 , , p .

3.1. Generalized Proportional Caputo Fractional Derivatives

Let the sequence of numbers 1 > α 1 > α 2 > > α N > 0 be given.
Consider the linear impulsive fractional differential equation with several generalized proportional Caputo fractional derivatives (LIGPCDE)
i = 1 N C i ( t k C D α i , ρ z ) ( t ) = F ( t ) , for t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , p , z ( t k + 0 ) = B k z ( t k ) , k = 1 , 2 , 3 , , p ,
with the boundary condition
z ( 0 ) + μ z ( T ) = Ψ ,
where ρ ( 0 , 1 ] , B k R ,   ( k = 1 , 2 , , p ) , C i R ,   ( i = 2 , N ) , C 1 = 1 , μ R , F : [ 0 , T ] R , Ψ R and
1 + μ e ρ 1 ρ T m = 0 p B m j = 1 N C j ρ α 1 α j Γ ( 1 + α 1 α j ) ( t m + 1 t m ) α 1 α j 0 .
Theorem 1.
Let the function F P I α 1 , ρ [ 0 , T ] and inequality (5) be fulfilled. Then the BVP for the LIGPCDE (3) (4) has a unique solution satisfying the integral presentation
z ( t ) = Ψ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ M p ( T ) m = 0 p 1 M m ( t m + 1 ) + μ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ M p ( T ) m = 0 p 1 M m ( t m + 1 ) × j = 2 N C j [ M p ( T ) m = 1 p l = 1 m 1 M p l ( t p l + 1 ) ( t p m I α 1 α j , ρ z ) ( t p m + 1 ) + ( t p I α 1 α j , ρ z ) ( T ) ]
μ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ M p ( T ) m = 0 p 1 M m ( t m + 1 ) M p ( T ) × m = 1 p l = 1 m 1 M p l ( t p l + 1 ) ( t p m I α 1 , ρ F ) ( t p m + 1 ) μ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ M p ( T ) m = 0 p 1 M m ( t m + 1 ) ( t p I α 1 , ρ F ) ( T ) j = 2 N C j [ M k ( t ) m = 1 k l = 1 m 1 M k l ( t k l + 1 ) ( t k m I α 1 α j , ρ z ) ( t k m + 1 ) + ( t k I α 1 α j , ρ z ) ( t ) ] + M k ( t ) m = 1 k l = 1 m 1 M k l ( t k l + 1 ) ( t k m I α 1 , ρ F ) ( t k m + 1 ) + ( t k I α 1 , ρ F ) ( t ) for t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , ,
where
A k ( t ) = j = 1 N C j ρ α 1 α j Γ ( 1 + α 1 α j ) ( t t k ) α 1 α j , k = 0 , 1 , 2 , , p ,
M k ( t ) = B k A k ( t ) e ρ 1 ρ ( t t k ) , k = 0 , 1 , 2 , , p , B 0 = 1 .
Proof. 
Let t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , p . Since α 1 > α j , j = 2 , 3 , , N , we take a GPFI ( t k I α 1 , ρ z ) ( t ) from both sides of (3), use Lemma 1 with τ = t k , α = α 1 α j , Lemma 3, Lemma 4 and get for any j = 2 , 3 , , N the equalities
t k I α 1 , ρ t k C D α j , ρ z ) ( t ) = t k I α 1 α j , ρ t k I α j , ρ t k C D α j , ρ z ) ( t ) = t k I α 1 α j , ρ z ( t ) z ( t k + 0 ) e ρ 1 ρ ( t t k ) = ( t k I α 1 α j , ρ z ) ( t ) z ( t k + 0 ) ( t k I α 1 α j , ρ e ρ 1 ρ ( t t k ) ) = ( t k I α 1 α j , ρ z ) ( t ) z ( t k + 0 ) e ρ 1 ρ ( t t k ) ρ α 1 α j Γ ( 1 + α 1 α j ) ( t t k ) α 1 α j .
From (3) on ( t k , t k + 1 ] , Lemma 4, Corollary 1 with τ = t k , α = α 1 α j and Equation (7) we get
z ( t ) = z ( t k + 0 ) e ρ 1 ρ ( t t k ) j = 2 N C j ( ( t k I α 1 α j , ρ z ) ( t ) z ( t k + 0 ) ( t k I α 1 α j , ρ e ρ 1 ρ ( t t k ) ) + ( t k I α 1 , ρ F ) ( t ) = z ( t k + 0 ) e ρ 1 ρ ( t t k ) A k ( t ) j = 2 N C j ( t k I α 1 α j , ρ z ) ( t ) + ( t k I α 1 , ρ F ) ( t ) , t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , p .
From (8) with k = 1 we obtain
z ( t 1 ) = z ( 0 ) M 0 ( t 1 ) j = 2 N C j ( t 0 I α 1 α j , ρ z ) ( t 1 ) + ( t 0 I α 1 , ρ F ) ( t 1 ) .
Apply Equation (9) to (8) with k = 1 , use the impulsive condition for k = 1 in (3) to get
z ( t ) = M 1 ( t ) z ( t 1 ) j = 2 N C j ( t 1 I α 1 α j , ρ z ) ( t ) + ( t 1 I α 1 , ρ F ) ( t ) = M 0 ( t 1 ) M 1 ( t ) z ( 0 ) j = 2 N C j M 1 ( t ) ( t 0 I α 1 α j , ρ z ) ( t 1 ) + ( t 1 I α 1 α j , ρ z ) ( t ) + M 1 ( t ) ( t 0 I α 1 , ρ F ) ( t 1 ) + ( t 1 I α 1 , ρ F ) ( t ) , t ( t 1 , t 2 ] .
From (8) with k = 2 , (10) for t = t 2 and the impulsive condition for k = 2 in (3) we obtain
z ( t ) = M 2 ( t ) z ( t 2 ) j = 2 N C j ( t 2 I α 1 α j , ρ z ) ( t ) + ( t 2 I α 1 , ρ F ) ( t ) = M 0 ( t 1 ) M 1 ( t 2 ) M 2 ( t ) z ( 0 ) j = 2 N C j [ M 2 ( t ) M 1 ( t 2 ) ( t 0 I α 1 α j , ρ z ) ( t 1 ) + M 2 ( t ) ( t 1 I α 1 α j , ρ z ) ( t 2 ) + ( t 2 I α 1 α j , ρ z ) ( t ) ] + M 1 ( t 2 ) M 2 ( t ) ( t 0 I α 1 , ρ F ) ( t 1 ) + M 2 ( t ) ( t 1 I α 1 , ρ F ) ( t 2 ) + ( t 2 I α 1 , ρ F ) ( t ) , t ( t 2 , t 3 ] .
By induction with respect to the intervals between two consecutive impulses, we obtain
z ( t ) = z ( 0 ) M k ( t ) m = 0 k 1 M m ( t m + 1 ) j = 2 N C j { m = 1 k l = 1 m 1 M k l ( t k l + 1 ) M k ( t ) ( t k m I α 1 α j , ρ z ) ( t k m + 1 ) + ( t k I α 1 α j , ρ z ) ( t ) } + m = 1 k l = 1 m 1 M k l ( t k l + 1 ) M k ( t ) ( t k m I α 1 , ρ F ) ( t k m + 1 ) + ( t k I α 1 , ρ F ) ( t ) , t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , p .
Therefore,
z ( T ) = z ( t p + 1 ) = z ( 0 ) M p ( T ) m = 0 p 1 M m ( t m + 1 ) j = 2 N C j [ m = 1 p l = 1 m 1 M p l ( t p l + 1 ) M p ( T ) ( t p m I α 1 α j , ρ z ) ( t p m + 1 ) + ( t p I α 1 α j , ρ z ) ( T ) ] + m = 1 p l = 1 m 1 M p l ( t p l + 1 ) M p ( T ) ( t p m I α 1 , ρ F ) ( t p m + 1 ) + ( t p I α 1 , ρ F ) ( T ) .
From Equation (13) and the boundary condition (3) we get
z ( 0 ) = Ψ 1 + μ M p ( T ) m = 0 p 1 M m ( t m + 1 ) + μ 1 + μ M p ( T ) m = 0 p 1 M m ( t m + 1 ) × j = 2 N C j [ m = 1 p l = 1 m 1 M p l ( t p l + 1 ) M p ( T ) ( t p m I α 1 α j , ρ z ) ( t p m + 1 ) + ( t p I α 1 α j , ρ z ) ( T ) ] μ 1 + μ M p ( T ) m = 0 p 1 M m ( t m + 1 ) m = 1 p ( l = 1 m 1 M p l ( t p l + 1 ) × M p ( T ) ( t p m I α 1 , ρ F ) ( t p m + 1 ) ) μ 1 + μ M p ( T ) m = 0 p 1 M m ( t m + 1 ) ( t p I α 1 , ρ F ) ( T ) .
Substitute equality (14) in (12) and obtain (7). □
Case 3.1.1. Generalized proportional Caputo fractional derivatives and no impulses.
Let us consider the partial case of (3) without impulses, i.e., t 0 = t 1 = t 2 = = t p . Thus, consider the BVP for the scalar linear differential equation with generalized proportional Caputo fractional derivatives (LGPCDE)
i = 1 N C i ( 0 C D α i , ρ z ) ( t ) = F ( t ) , f o r t ( 0 , T ] ,
with the boundary condition (4). Then as a partial case of Equation (7) it follows the solution of BVP for LGPCDE (15), (4) is
z ( t ) = Ψ M ( t ) 1 + μ M ( T ) + μ M ( t ) 1 + μ M ( T ) j = 2 N C j ( 0 I α 1 α j , ρ z ) ( T ) μ M ( t ) 1 + μ M ( T ) ( 0 I α 1 , ρ F ) ( T ) j = 2 N C j ( 0 I α 1 α j , ρ z ) ( t ) + ( 0 I α 1 , ρ F ) ( t ) , t ( 0 , T ] ,
where M ( t ) = e ρ 1 ρ t j = 1 N C j ρ α 1 α j Γ ( 1 + α 1 α j ) t α 1 α j , with
1 + μ e ρ 1 ρ T j = 1 N C j ρ α 1 α j Γ ( 1 + α 1 α j ) T α 1 α j 0 .
(Compare with [14] with γ = 1 , β = 0 , Ψ ( ξ ) c o n s t ).

3.2. Caputo Fractional Derivatives

In this section, we will use the fractional integral and the fractional derivative defined in Remark 1.
Let α ( 0 , 1 ) . Define the classes of functions
PCC α [ 0 , T ] = ν PC 1 [ 0 , T ] ) : ( t k D α ν ) ( t ) exists   for t ( t k , t k + 1 ] , k = 0 , 1 , , p , PCI α [ 0 , T ] = ν PC [ 0 , T ] : ( t k I α ν ) ( t ) exists   for t ( t k , t k + 1 ] , k = 0 , 1 , , p .
We will study the following linear impulsive fractional differential equation with several Caputo fractional derivatives (LICDE)
i = 1 N C i t k C D α i , ρ z ( t ) = F ( t ) , for t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , p , z ( t k + 0 ) = B k z ( t k ) , k = 1 , 2 , 3 , , p ,
where t k C D α i , ρ u ( t ) is the Caputo fractional derivative (see Remark 1), B k R ,   ( k = 1 , 2 , , p ) , C i R , ( i = 2 , N ) , C 1 = 1 , μ R , F : [ 0 , T ] R and
1 + μ m = 0 p B m j = 1 N C j ρ α 1 α j Γ ( 1 + α 1 α j ) ( t m + 1 t m ) α 1 α j 0 .
Theorem 2.
Let F P C I α 1 [ 0 , T ] and inequality (18) holds. Then the BVP for the LICDE (17) (4) has a unique solution
z ( t ) = Ψ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ M p ( T ) m = 0 p 1 M m ( t m + 1 ) + μ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ M p ( T ) m = 0 p 1 M m ( t m + 1 ) × j = 2 N C j [ M p ( T ) m = 1 p l = 1 m 1 M p l ( t p l + 1 ) t p m I α 1 α j z ( t p m + 1 ) + t p I α 1 α j z ( T ) ] μ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ M p ( T ) m = 0 p 1 M m ( t m + 1 ) M p ( T ) × m = 1 p l = 1 m 1 M p l ( t p l + 1 ) t p m I α 1 F ( t p m + 1 ) μ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ M p ( T ) m = 0 p 1 M m ( t m + 1 ) t p I α 1 F ( T ) j = 2 N C j [ M k ( t ) m = 1 k l = 1 m 1 M k l ( t k l + 1 ) t k m I α 1 α j z ( t k m + 1 ) + t k I α 1 α j z ( t ) ] + M k ( t ) m = 1 k l = 1 m 1 M k l ( t k l + 1 ) t k m I α 1 F ( t k m + 1 ) + t k I α 1 F ( t ) , t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , p ,
where M k ( t ) = B k j = 1 N C j Γ ( 1 + α 1 α j ) ( t t k ) α 1 α j , k = 0 , 1 , 2 , , p , B 0 = 1 .
Case 3.2.1. Several Caputo fractional derivatives and no impulses
Consider the partial case of (17) without impulses, i.e., t 0 = t 1 = t 2 = = t p . Thus, consider the scalar linear differential equation with Caputo fractional derivatives (LCDE)
i = 1 N C i 0 C D α i , ρ z ( t ) = F ( t ) for t ( 0 , T ] ,
where C i R , ( i = 2 , N ) , C 1 = 1 , μ R , F : [ 0 , T ] R .
As a partial case of Equation (19) the solution of BVP for LCDE (20) (4) is
z ( t ) = Ψ M ( t ) 1 + μ M ( T ) + μ M ( t ) 1 + μ M ( T ) j = 2 N C j 0 I α 1 α j z ( T ) μ M ( t ) 1 + μ M ( T ) 0 I α 1 F ( T ) j = 2 N C j 0 I α 1 α j z ( t ) + 0 I α 1 F ( t ) , t ( 0 , T ] ,
where M ( t ) = j = 1 N C j Γ ( 1 + α 1 α j ) t α 1 α j with 1 + μ j = 1 N C j Γ ( 1 + α 1 α j ) T α 1 α j 0 .
Case 3.2.2. One Caputo fractional derivative and no impulses
Consider the partial case of Equation (20) without impulses, i.e., t 0 = t 1 = t 2 = = t p and C k = 0 , k = 2 , 3 , , N . We will study the scalar linear Caputo fractional differential equation
0 C D α z ( t ) = F ( t ) for t ( 0 , T ] ,
where α ( 0 , 1 ) , F : [ 0 , T ] R
As a partial case of Equation (21), the solution of the BVP (22) (4) is
z ( t ) = Ψ 1 + μ μ 1 + μ 1 Γ ( α ) 0 T F ( s ) T s 1 α ds + 1 Γ ( α ) 0 t F ( s ) t s 1 α ds , t ( 0 , T ] ,
where 1 + μ 0 .
Case 3.2.3. One Caputo fractional derivative and no impulses—Initial value problem
Consider the scalar linear Caputo fractional differential Equation (22) with boundary condition (4) with μ = 0 . Then the solution of (22) with the initial condition z ( 0 ) = Ψ is given by
z ( t ) = Ψ + 1 Γ ( α ) 0 t F ( s ) t s 1 α ds .
(see, for example, [15]).

3.3. One GPCFD and Impulses

Consider the following impulsive linear generalized proportional Caputo fractional differential equation (IGPFDE)
( t k C D α , ρ z ) ( t ) = F ( t ) , for t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , p , z ( t k + 0 ) = B k z ( t k ) , k = 1 , 2 , 3 , , p ,
where α ( 0 , 1 ) , ρ ( 0 , 1 ] , B k R , ( k = 1 , 2 , , p ) , μ R , F : [ 0 , T ] R and
1 + μ e ρ 1 ρ ( T t 0 ) m = 0 p B m 0 .
As a partial case of Theorem 1 and Formula (7) with C i = 0 , i = 2 , 3 , , N , we obtain
Theorem 3.
Let F P I α , ρ [ 0 , T ] and inequality (25) be satisfied. Then the BVP for the IGPFDE (25) (4) has an unique solution given by
z ( t ) = Ψ e ρ 1 ρ ( t t 0 ) m = 0 k B m 1 + μ e ρ 1 ρ ( T t 0 ) m = 0 p B m + μ e ρ 1 ρ ( t t 0 ) m = 0 k B m 1 + μ e ρ 1 ρ ( T t 0 ) m = 0 p B m ( t p I α , ρ F ) ( T ) μ e ρ 1 ρ ( t t 0 ) m = 0 k B m 1 + μ e ρ 1 ρ ( T t 0 ) m = 0 p B m
× m = 1 p l = 0 m 1 B p l e ρ 1 ρ ( T t p m + 1 ) ( t p m I α , ρ F ) ( t p m + 1 ) + m = 1 k l = 0 m 1 B k l e ρ 1 ρ ( t t k m + 1 ) ( t k m I α , ρ F ) ( t k m + 1 ) + ( t k I α , ρ F ) ( t ) , t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , p .
Case 3.3.1. GPCFD and no impulse
We will study the partial case of Equation (24) without impulses, i.e., t 0 = t 1 = t 2 = = t p . Thus, consider the scalar linear generalized proportional Caputo fractional differential equation (GPFDE)
( 0 C D α i , ρ z ) ( t ) = F ( t ) , t ( 0 , T ] ,
where μ R , F : [ 0 , T ] R .
As a partial case of Equation (27) the solution of BVP for GPFDE (27), (4) is given by
z ( t ) = Ψ e ρ 1 ρ t 1 + μ e ρ 1 ρ T μ e ρ 1 ρ t 1 + μ e ρ 1 ρ T 1 ρ α Γ ( α ) 0 T e ρ 1 ρ ( T s ) F ( s ) T s 1 α ds + 1 ρ α Γ ( α ) 0 t e ρ 1 ρ ( t s ) F ( s ) t s 1 α ds , t ( 0 , T ] ,
where 1 + μ e ρ 1 ρ T 0 .
Case 3.3.2. GPCFD and no impulse—Initial value problem
Consider GPFDE (27) with the boundary condition (4) with μ = 0 . Then as a partial case of (28) the GPFDE (27) with the initial condition z ( 0 ) = Ψ has a solution
z ( t ) = Ψ e ρ 1 ρ t + 1 ρ α Γ ( α ) 0 t e ρ 1 ρ t s F ( s ) t s 1 α ds ,
(see, Example 5.7 [8] with λ = 0 ).

4. Nonlinear Impulsive Delay Differential Equations with Several GPCFDE

Let the numbers 1 > α 1 > α 2 > > α N > 0 and the points 0 = t 0 < t 1 < t 2 < < t p < t p + 1 = T be given.
Consider the nonlinear impulsive delay differential equation with several GPCDE (NIGPDE)
i = 1 N C i ( t k C D α i , ρ x ) ( t ) = f ( t , x ( t ) , x ( λ t ) ) , f o r t ( t k , t k + 1 ] , k = 0 , 1 , , p , x ( t k + 0 ) = B k x ( t k ) , k = 1 , 2 , , p ,
with the boundary value condition (4), where λ ( 0 , 1 ) , C i , i = 2 , N , C 1 = 1 , B i R , i = 1 , 2 , 3 , , p , f : [ 0 , T ] × R 2 R .

4.1. Mild Solution of the BVP for NIGPDE

Based on the integral form (7) of the solution of BVP for LIGPCDE (3) (4) we define the fractional integral operator Ω : k = 2 N P I α 1 α k , ρ [ 0 , T ] P C [ 0 , T ] by the equality
( Ω x ) ( t ) = Ψ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) + μ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) j = 2 N C j [ M p ( T ) m = 1 p l = 1 m 1 M p l ( t p l + 1 ) × 1 ρ α 1 α j Γ ( α 1 α j ) t p m t p m + 1 e ρ 1 ρ t p m + 1 s x ( s ) t p m + 1 s 1 α 1 + α j ds + 1 ρ α 1 α j Γ ( α 1 α j ) t p T e ρ 1 ρ T s x ( s ) T s 1 α 1 + α j ds ] μ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) M p ( T ) × m = 1 p l = 1 m 1 M p l ( t p l + 1 ) 1 ρ α 1 Γ ( α 1 ) t p m t p m + 1 e ρ 1 ρ t p m + 1 s f ( s , x ( s ) , x ( λ s ) ) t p m + 1 s 1 α 1 ds μ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) 1 ρ α 1 Γ ( α 1 ) t p T e ρ 1 ρ T s f ( s , x ( s ) , x ( λ s ) ) T s 1 α 1 ds j = 2 N C j [ M k ( t ) m = 1 k l = 1 m 1 M k l ( t k l + 1 ) × 1 ρ α 1 α j Γ ( α 1 α j ) t k m t k m + 1 e ρ 1 ρ t k m + 1 s x ( s ) t k m + 1 s 1 α 1 + α j ds + 1 ρ α 1 α j Γ ( α 1 α j ) t k t e ρ 1 ρ t s x ( s ) t s 1 α 1 + α j ds ] + M k ( t ) m = 1 k l = 1 m 1 M k l ( t k l + 1 ) 1 ρ α 1 Γ ( α 1 ) t k m t k m + 1 e ρ 1 ρ t k m + 1 s f ( s , x ( s ) , x ( λ s ) ) t k m + 1 s 1 α 1 ds + 1 ρ α 1 Γ ( α 1 ) t k t e ρ 1 ρ t s f ( s , x ( s ) , x ( λ s ) ) t s 1 α 1 ds , t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , p ,
where
M 0 ( t ) = e ρ 1 ρ t j = 1 N C j ρ α 1 α j Γ ( 1 + α 1 α j ) t α 1 α j ,
and
M k ( t ) = B k e ρ 1 ρ ( t t k ) j = 1 N C j ρ α 1 α j Γ ( 1 + α 1 α j ) ( t t k ) α 1 α j , k = 1 , 2 , , p .
Definition 1.
The function x PI α 1 α k , ρ [ 0 , T ] , k = 2 , 3 , , N , is a mild solution of the BVP for NIGPDE (29) (4) if it is a fixed point of the operator Ω, defined by (30).
Theorem 4.
Let the inequality (5) hold and x ( t ) , t [ 0 , T ] , be a mild solution of the BVP for NIGPDE (29) (4) such that x P C α k , ρ [ 0 , T ] for k = 1 , 2 , , N . Then the function x ( t ) is a solution to the same problem.
Proof. 
From Equation (30) applying M 0 ( 0 ) = 1 and the following two equalities
x ( T ) = Ψ m = 0 p M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) + μ m = 0 p M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) 1 j = 2 N C j M p ( T ) m = 1 p l = 1 m 1 M p l ( t p l + 1 ) × 1 ρ α 1 α j Γ ( α 1 α j ) t p m t p m + 1 e ρ 1 ρ t p m + 1 s x ( s ) t p m + 1 s 1 α 1 + α j ds + μ m = 0 p M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) 1 j = 2 N C j 1 ρ α 1 α j Γ ( α 1 α j ) × t p T e ρ 1 ρ T s x ( s ) T s 1 α 1 + α j ds μ m = 0 p M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) 1 M p ( T ) m = 1 p l = 1 m 1 M p l ( t p l + 1 ) 1 ρ α 1 Γ ( α 1 ) × t p m t p m + 1 e ρ 1 ρ t p m + 1 s f ( s , x ( s ) , x ( λ s ) ) t p m + 1 s 1 α 1 ds μ m = 0 p M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) 1 1 ρ α 1 Γ ( α 1 ) t p T e ρ 1 ρ T s f ( s , x ( s ) , x ( λ s ) ) T s 1 α 1 ds ,
and
x ( 0 ) = Ψ 1 + μ m = 0 p M m ( t m + 1 ) + μ 1 + μ m = 0 p M m ( t m + 1 ) j = 2 N C j M p ( T ) m = 1 p l = 1 m 1 M p l ( t p l + 1 ) × 1 ρ α 1 α j Γ ( α 1 α j ) t p m t p m + 1 e ρ 1 ρ t p m + 1 s x ( s ) t p m + 1 s 1 α 1 + α j ds + μ 1 + μ m = 0 p M m ( t m + 1 ) j = 2 N C j 1 ρ α 1 α j Γ ( α 1 α j ) t p T e ρ 1 ρ T s x ( s ) T s 1 α 1 + α j ds μ 1 + μ m = 0 p M m ( t m + 1 ) M p ( T ) m = 1 p l = 1 m 1 M p l ( t p l + 1 ) 1 ρ α 1 Γ ( α 1 ) × t p m t p m + 1 e ρ 1 ρ t p m + 1 s f ( s , x ( s ) , x ( λ s ) ) t p m + 1 s 1 α 1 ds μ 1 + μ m = 0 p M m ( t m + 1 ) 1 ρ α 1 Γ ( α 1 ) t p T e ρ 1 ρ T s f ( s , x ( s ) , x ( λ s ) ) T s 1 α 1 ds ,
we obtain
x ( 0 ) + μ x ( T ) = Ψ 1 + μ m = 0 p M m ( t m + 1 ) + μ Ψ m = 0 p M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) + μ 1 + μ m = 0 p M m ( t m + 1 ) j = 2 N C j M p ( T ) m = 1 p l = 1 m 1 M p l ( t p l + 1 ) × 1 ρ α 1 α j Γ ( α 1 α j ) t p m t p m + 1 e ρ 1 ρ t p m + 1 s x ( s ) t p m + 1 s 1 α 1 + α j ds + μ μ m = 0 p M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) 1 j = 2 N C j M p ( T ) m = 1 p l = 1 m 1 M p l ( t p l + 1 ) × 1 ρ α 1 α j Γ ( α 1 α j ) t p m t p m + 1 e ρ 1 ρ t p m + 1 s x ( s ) t p m + 1 s 1 α 1 + α j ds + μ 1 + μ m = 0 p M m ( t m + 1 ) j = 2 N C j 1 ρ α 1 α j Γ ( α 1 α j ) t p T e ρ 1 ρ T s x ( s ) T s 1 α 1 + α j ds + μ μ m = 0 p M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) 1 j = 2 N C j 1 ρ α 1 α j Γ ( α 1 α j ) t p T e ρ 1 ρ T s x ( s ) T s 1 α 1 + α j ds μ 1 + μ m = 0 p M m ( t m + 1 ) M p ( T ) × m = 1 p l = 1 m 1 M p l ( t p l + 1 ) 1 ρ α 1 Γ ( α 1 ) t p m t p m + 1 e ρ 1 ρ t p m + 1 s f ( s , x ( s ) , x ( λ s ) ) t p m + 1 s 1 α 1 ds μ μ m = 0 p M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) 1 M p ( T ) × m = 1 p l = 1 m 1 M p l ( t p l + 1 ) 1 ρ α 1 Γ ( α 1 ) t p m t p m + 1 e ρ 1 ρ t p m + 1 s f ( s , x ( s ) , x ( λ s ) ) t p m + 1 s 1 α 1 ds μ 1 + μ m = 0 p M m ( t m + 1 ) 1 ρ α 1 Γ ( α 1 ) t p T e ρ 1 ρ T s f ( s , x ( s ) , x ( λ s ) ) T s 1 α 1 ds μ μ m = 0 p M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) 1 1 ρ α 1 Γ ( α 1 ) t p T e ρ 1 ρ T s f ( s , x ( s ) , x ( λ s ) ) T s 1 α 1 ds = Ψ .
Equalities (35) proves the boundary condition (4) holds for x ( t ) .
Define w j ( t ) = ( t k I α 1 α j , ρ x ) ( t ) , j = 2 , 3 , , N , t ( t k , t k 1 ] , k = 1 , 2 , , p . According to Corollary 2, the equalities
t k C D α 1 , ρ t k I α 1 α j , ρ x ( t ) = t k C D α 1 , ρ w j ( t ) = t k C D α 1 , ρ t k C D α j , ρ t k I α j , ρ w j ( t )
hold. According to Lemma 2 we get
t k C D α 1 , ρ t k I α 1 α j , ρ x ( t ) = 0 C D α j , ρ t k C D α 1 , ρ t k I α 1 , ρ x ( t ) = t k C D α j , ρ x ( t ) .
According to Lemma 1 with α = α 1 α j , we have the equality
1 ρ α 1 α j Γ ( 1 + α 1 α j ) e ρ 1 ρ ( t t k ) ( t t k ) α 1 α j = t k I α 1 α j , ρ e ρ 1 ρ t
and in view of Lemma 3 with α = α j , y = x we get
x ( t ) x ( t k + 0 ) e ρ 1 ρ ( t t k ) = t k I α j , ρ t k C D α j , ρ x ( t ) .
Then using (30) for t k + 0 ( t k , t k + 1 ] ) , M k ( t k + 0 ) = B k , k = 1 , 2 , , p , we obtain x ( t k + 0 ) = B k x ( t k ) , i.e., the mild solution is satisfying the impulsive condition in (29).
Thus, equality x ( t ) = ( Ω x ) ( t ) for the mild solution, with operator Ω defined by (30). Thus,
x ( t ) = M k ( t ) Ξ k + j = 2 N C j ( t k I α 1 α j , ρ x ) ( t ) + ( t k I α 1 , ρ f ( . , x ( . ) , x ( λ . ) ) ( t ) t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , p .
where
Ξ k = Ψ m = 0 k 1 M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) + μ m = 0 k 1 M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) j = 2 N C j [ M p ( T ) m = 1 p l = 1 m 1 M p l ( t p l + 1 ) × ( t p m I α 1 α j , ρ x ) ( t p m + 1 ) + ( t p I α 1 α j , ρ x ) ( T ) ] μ m = 0 k 1 M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) M p ( T ) m = 1 p l = 1 m 1 M p l ( t p l + 1 ) × ( t p m I α 1 , ρ f ( . , x ( . ) , x ( λ 1 . ) ) ) ( t p m + 1 ) μ m = 0 k 1 M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) ( t p I α 1 , ρ f ( . , x ( . ) , x ( λ . ) ) ( T ) j = 2 N C j m = 1 k l = 1 m 1 M k l ( t k l + 1 ) ( t k m I α 1 α j , ρ x ) ( t k m + 1 ) + M k ( t ) m = 1 k l = 1 m 1 M k l ( t k l + 1 ) ( t k m I α 1 , ρ f ( . , x ( . ) , x ( λ . ) ) ( t k m + 1 ) .
We take the GPCFD t k C D α 1 , ρ of both sides of (37) apply Equation (36), Lemma 3, Corollary 2, the equality
t k C D α 1 , ρ t k I α 1 α j , ρ e ρ 1 ρ ( . t k ) ( t ) = t k C D α j , ρ e ρ 1 ρ ( , t k ) ( t ) = 0
and obtain
t k C D α 1 , ρ x ( t ) = j = 2 N C j t k C D α 1 , ρ t k I α 1 α j , ρ x ( t ) + t k C D α 1 , ρ t k I α 1 , ρ f . , x ( . ) , x ( λ . ) ( t ) = j = 2 N C j t k C D α j , ρ x ( t ) + f ( t , x ( t ) , x ( λ t ) ) , t ( t k , t k + 1 ] , k = 0 , 1 , , p .
It proves the claim. □
Theorem 5.
Let the inequality (5) hold, x ( t ) be a solution of the BVP for NIGPDE (29) (4) and F P I α 1 , ρ [ 0 , T ] , where F ( t ) = f ( t , x ( t ) , x ( λ t ) . Then the function x ( t ) is a mild solution of the same problem.
The proof of Theorem 5 is similar to the one of Theorem 1 and we omit it.
We will obtain some sufficient conditions for the existence of the mild solutions of (29) (4). The proofs are based on fixed point theorems and the application of the operator Ω .
Theorem 6.
Let:
1.
The inequality (44) holds.
2.
There exist constants L , M > 0 :
| f ( t , x 1 , z 1 ) f ( t , x 2 , z 2 ) | L | x 1 x 2 | + M | z 1 z 2 | , t [ 0 , T ] , x i , z i R , i = 1 , 2 .
3.
For  k = 0 , 1 , 2 , , p ,  the inequality
( j = 2 N | C j | ξ α 1 α j Γ ( 1 + α 1 α j ) + ( L + M ) ξ α 1 Γ ( 1 + α 1 ) ) Φ k < 1 ,
holds where
ξ k = t k + 1 t k ρ , k = 0 , 1 , 2 , , p , ξ = max k = 0 , 1 , 2 , , p ξ k , σ m = j = 1 N | C j | ξ m α j Γ ( 1 + α 1 α j ) , m = 0 , 1 , 2 , , p , β m = j = 1 N C j ξ m α j Γ ( 1 + α 1 α j ) , m = 0 , 1 , 2 , , p , η = 1 + μ m = 0 p B m ξ m α 1 e ( ρ 1 ) ξ m β m , Φ k = | μ | m = 0 k | B m | ξ m α 1 σ m η m = 0 p l = 0 m 1 | B p l | ξ p l α 1 σ p l + m = 0 k l = 0 m 1 | B k l | ξ k l α 1 σ k l .
Then the BVP for NIGPDE (29) (4) has a unique mild solution. 
Proof. 
We use the norm | | x | | = sup t [ 0 , T ] | x ( t ) | for any x P C [ 0 , T ] .
Then | M k ( t ) m = 0 k 1 M m ( t m + 1 ) | m = 0 k | B m ξ m α 1 | j = 1 N | C j | ξ m α j Γ ( 1 + α 1 α j ) = m = 0 k | B m | ξ m α 1 σ m , where M k ( t ) , k = 0 , 1 , 2 , , p are defined by (31) and (32).
Consider the integral fractional operator Ω defined by (30). Let x , y P C ( [ 0 , T ] ) . Then t k t e ρ 1 ρ t s t s 1 α ds ( t t k ) α α . Thus, we obtain for t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , p ,
| Ω x ( t ) Ω y ( t ) | { μ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) j = 2 N C j [ M p ( T ) × m = 1 p l = 1 m 1 M p l ( t p l + 1 ) ( t p m + 1 t p m ) α 1 α j ρ α 1 α j Γ ( 1 + α 1 α j ) + ( T t p ) α 1 α j ) ρ α 1 α j Γ ( 1 + α 1 α j ) ] + ( L + M ) μ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) M p ( T ) × m = 1 p l = 1 m 1 M p l ( t p l + 1 ) ( t p m + 1 t p m ) α 1 ρ α 1 Γ ( 1 + α 1 )
+ ( L + M ) μ M k ( t ) m = 0 k 1 M m ( t m + 1 ) 1 + μ m = 0 p M m ( t m + 1 ) ( T t p ) α 1 ρ α 1 Γ ( 1 + α 1 ) + j = 2 N C j [ M k ( t ) m = 1 k l = 1 m 1 M k l ( t k l + 1 ) ( t k m + 1 t k m ) α 1 α j ρ α 1 α j Γ ( 1 + α 1 α j ) + ( t t k ) α 1 α j ρ α 1 α j Γ ( 1 + α 1 α j ) ] + ( L + M ) M k ( t ) m = 1 k l = 1 m 1 M k l ( t k l + 1 ) ( t k m + 1 t k m ) α 1 ρ α 1 Γ ( 1 + α 1 ) + ( L + M ) ( t t k ) α 1 ρ α 1 Γ ( 1 + α 1 ) } | | x y | |
or
| Ω x ( t ) Ω y ( t ) | { | μ | m = 0 k | B m | ξ m α 1 σ m η j = 2 N | C j | [ m = 1 p l = 0 m l | B p l | ξ p l α 1 σ p l × ξ p m α 1 α j Γ ( 1 + α 1 α j ) + ξ p α 1 α j Γ ( 1 + α 1 α j ) ] + j = 2 N | C j | m = 1 k l = 0 m 1 | B k l | ξ k l α 1 σ k l ξ k m α 1 α j Γ ( 1 + α 1 α j ) + ξ k α 1 α j Γ ( 1 + α 1 α j ) + ( L + M ) [ | μ | m = 0 k | B m | ξ m α 1 σ m η m = 1 p l = 0 m 1 | B p l | ξ p l α 1 σ p l + m = 1 k l = 0 m 1 | B k l | ξ k l α 1 σ k l + | μ | m = 0 k | B m | ξ m α 1 σ m η + 1 ] ξ α 1 Γ ( 1 + α 1 ) } | | x y | | { j = 2 N | C j | ξ α 1 α j Γ ( 1 + α 1 α j ) [ | μ | m = 0 k | B m | ξ m α 1 σ m η m = 0 p l = 0 m l | B p l | ξ p l α 1 σ p l + m = 0 k l = 0 m 1 | B k l | ξ k l α 1 σ k l ] + ( L + M ) [ | μ | m = 0 k | B m | ξ m α 1 σ m η m = 0 p l = 0 m 1 | B p l | ξ p l α 1 σ p l + m = 0 k l = 0 m 1 | B k l | ξ k l α 1 σ k l ] ξ α 1 Γ ( 1 + α 1 ) } | | x y | | = j = 2 N | C j | ξ α 1 α j Γ ( 1 + α 1 α j ) + ( L + M ) ξ α 1 Γ ( 1 + α 1 ) Φ k | | x y | | .
Inequality (42) and condition 3 proves the existence of a constant K ( 0 , 1 ) such that
| | Ω x Ω y | | K | | x 1 y | | .
Thus, the Banach contraction principle applied to the operator Ω proves that there exists a unique fixed point x * P C [ 0 , T ] of Ω , which is a mild solution of (29) and (4). □

4.2. Mild Solution of BVP for the Impulsive Delay Fractional Differential Equation with GPCDE

We consider the nonlinear impulsive delay differential equation with one GPCDE (NIGDE) of the type
( t k C D α , ρ x ) ( t ) = f ( t , x ( t ) , x ( λ t ) ) ) , for t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , p z ( t k + 0 ) = B k z ( t k ) , k = 1 , 2 , 3 , , p ,
with the boundary value condition (4), where α ( 0 , 1 ) , ρ ( 0 , 1 ] , λ ( 0 , 1 ) , the numbers B i R , k = 1 , 2 , 3 , , p , the function f : [ 0 , T ] × R 2 R and
1 + μ e ρ 1 ρ T m = 0 p B m 0 .
Based on the integral presentation (24) we define the integral fractional operator W : PC [ 0 , T ] PC [ 0 , T ] by
( W x ) ( t ) = Ψ e ρ 1 ρ ( t t 0 ) m = 0 k B m 1 + μ e ρ 1 ρ ( T t 0 ) m = 0 p B m + μ e ρ 1 ρ ( t t 0 ) m = 0 k B m 1 + μ e ρ 1 ρ ( T t 0 ) m = 0 p B m ( t p I α 1 , ρ f ( . , x ( . ) , x ( λ . ) ) ) ( T ) μ e ρ 1 ρ ( t t 0 ) m = 0 k B m 1 + μ e ρ 1 ρ ( T t 0 ) m = 0 p B m × m = 1 p l = 0 m 1 B p l e ρ 1 ρ ( T t p m + 1 ) ( t p m I α 1 , ρ f ( . , x ( . ) , x ( λ . ) ) ) ( t p m + 1 ) + m = 1 k l = 0 m 1 B k l e ρ 1 ρ ( t t k m + 1 ) ( t k m I α 1 , ρ f ( . , x ( . ) , x ( λ . ) ) ) ( t k m + 1 ) + ( t k I α 1 , ρ f ( . , x ( . ) , x ( λ . ) ) ) ( t ) , t ( t k , t k + 1 ] , k = 0 , 1 , 2 , , p .
Definition 2.
The function x ( t ) is a mild solution of the BVP for NIGPDE (43) (4) if it is a fixed point of the integral fractional operator W , defined by (45).
Theorem 7.
Let:
1.
The inequality (5) holds.
2.
There exist constants L , M > 0 :
| f ( t , x 1 , z 1 ) f ( t , x 2 , z 2 ) | L | x 1 x 2 | + M | z 1 z 2 | , t [ 0 , T ] , x i , z i R , i = 1 , 2 .
3.
For k = 0 , 1 , 2 , , p , the inequality
( L + M ) ξ α Γ ( 1 + α ) Φ k < 1 ,
holds where
ξ = max k = 0 , 1 , 2 , , p t k + 1 t k ρ , η = 1 + μ e ρ 1 ρ T m = 0 p B m ξ m α 1 , Φ k = | μ | m = 0 k | B m | ξ m α 1 η m = 0 p l = 0 m 1 | B p l | ξ p l α 1 + m = 0 k l = 0 m 1 | B k l | ξ k l α 1 .
Then the BVP for NIGPDE (43) (4) has a unique mild solution. 
The claim of Theorem 7 follows from Theorem 6.
Remark 2.
From the above results for BVP for NIGPDE, we obtain some results for BVP for impulsive fractional differential equations with Caputo fractional derivatives with lower changeable limits at the impulsive points. The case of Caputo fractional derivative with a fixed lower limit at the initial time is studied in [16].

5. Conclusions

In this paper, two main goals are reached. First, an integral form of the unique solutions of the BVP for various types of linear impulsive fractional differential equations is obtained. Second, based on these presentations, new types of integral fractional operators are defined. These operators are applied to obtain sufficient conditions for the uniqueness and existence of mild solutions of BVP for several types of nonlinear delay impulsive fractional differential equations.
The obtained integral presentations of the solutions of the BVP for the linear impulsive fractional differential equations could be applied in some approximate methods such as the monotone-iterative technique for constructing successive approximations. The newly defined integral fractional operator could be successfully used to study the Ulam-type stability of BVP for various types of nonlinear fractional differential equations with impulses.
The ideas in this paper are also applicable to studying various kinds of boundary conditions, such as nonlinear ones, nonlocal ones, and integral ones. Furthermore, the case of the fractional order α > 1 could be investigated.

Author Contributions

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

Funding

The work is partially supported by the Bulgarian National Science Fund under Project KP-06-PN62/1.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Haider, S.S.; Rehman, M.; Abdeljawad, T. On Hilfer fractional difference operator. Adv. Differ. Equ. 2020, 2020, 122. [Google Scholar] [CrossRef]
  2. Furati, K.M.; Kassim, M.D. Existence and uniqueness for a problem involving Hilfer fractional derivative. Comput. Math. Appl. 2012, 64, 1616–1626. [Google Scholar] [CrossRef]
  3. Almeida, R. A Caputo fractional derivative of a function with respect to another function. Commun. Nonlinear Sci. Numer. Simul. 2017, 44, 460–481. [Google Scholar] [CrossRef]
  4. Osler, T.J. Fractional derivatives of a composite function. SIAM J. Math. Anal. 1970, 1, 288–293. [Google Scholar] [CrossRef]
  5. Sousa, J.V.C.; Oliveira, E.C. On the g-Hilfer fractional derivative. Commun. Nonlinear Sci. Numer. Simul. 2018, 60, 72–91. [Google Scholar] [CrossRef]
  6. Odzijewicz, T.; Malinowska, A.B.; Torres, D.F.M. Generalized fractional calculus with applications to the calculus of variations. Comput. Math. Appl. 2012, 64, 3351–3366. [Google Scholar] [CrossRef]
  7. Odzijewicz, T.; Malinowska, A.B.; Torres, D.F.M. Fractional calculus of variations in terms of a generalized fractional integral with applications to physics. Abstr. Appl. Anal. 2012, 2012, 871912. [Google Scholar] [CrossRef]
  8. 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]
  9. Jarad, F.; Abdeljawad, T. Generalized fractional derivatives and Laplace transform. Discret. Contin. Dyn. Syst. Ser. S 2020, 13, 709–722. [Google Scholar] [CrossRef]
  10. Agarwal, R.P.; Hristova, S.; O’Regan, D. Stability of generalized proportional Caputo fractional differential equations by lyapunov functions. Fractal Fract. 2022, 6, 34. [Google Scholar] [CrossRef]
  11. Almeida, R.; Agarwal, R.P.; Hristova, S.; O’Regan, D. Stability of gene regulatory networks modeled by generalized proportional Caputo fractional differential equations. Entropy 2022, 24, 372. [Google Scholar] [CrossRef]
  12. Bohner, M.; Hristova, S. Stability for generalized Caputo proportional fractional delay integro-differential equations. Bound. Value Probl. 2022, 2022, 14. [Google Scholar] [CrossRef]
  13. Barakat, M.A.; Soliman, A.H.; Hyder, A. Langevin equations with generalized proportional Hadamard-Caputo fractional derivative. Comput. Intell. Neurosci. 2021, 2021, 6316477. [Google Scholar] [CrossRef] [PubMed]
  14. Agarwal, R.P.; Hristova, S. Boundary Value Problem for Multi-Term Nonlinear Delay Generalized Proportional Caputo Fractional Differential Equations. Fractal Fract. 2022, 6, 691. [Google Scholar] [CrossRef]
  15. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  16. Sun, Y.; Zeng, Z.; Song, J. Existence and Uniqueness for the Boundary Value Problems of Nonlinear Fractional Differential Equation. Appl. Math. 2017, 8, 3. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Agarwal, R.P.; Hristova, S. Boundary Value Problem for Impulsive Delay Fractional Differential Equations with Several Generalized Proportional Caputo Fractional Derivatives. Fractal Fract. 2023, 7, 396. https://doi.org/10.3390/fractalfract7050396

AMA Style

Agarwal RP, Hristova S. Boundary Value Problem for Impulsive Delay Fractional Differential Equations with Several Generalized Proportional Caputo Fractional Derivatives. Fractal and Fractional. 2023; 7(5):396. https://doi.org/10.3390/fractalfract7050396

Chicago/Turabian Style

Agarwal, Ravi P., and Snezhana Hristova. 2023. "Boundary Value Problem for Impulsive Delay Fractional Differential Equations with Several Generalized Proportional Caputo Fractional Derivatives" Fractal and Fractional 7, no. 5: 396. https://doi.org/10.3390/fractalfract7050396

APA Style

Agarwal, R. P., & Hristova, S. (2023). Boundary Value Problem for Impulsive Delay Fractional Differential Equations with Several Generalized Proportional Caputo Fractional Derivatives. Fractal and Fractional, 7(5), 396. https://doi.org/10.3390/fractalfract7050396

Article Metrics

Back to TopTop