Next Article in Journal / Special Issue
Generalized Fractional Integrals Involving Product of a Generalized Mittag–Leffler Function and Two H-Functions
Previous Article in Journal
Fundamental Spacetime Representations of Quantum Antenna Systems
Previous Article in Special Issue
On Fractional Lyapunov Functions of Nonlinear Dynamic Systems and Mittag-Leffler Stability Thereof
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Note on a Coupled System of Hilfer Fractional Differential Inclusions

by
Aurelian Cernea
1,2
1
Faculty of Mathematics and Computer Science, University of Bucharest, Academiei 14, 010014 Bucharest, Romania
2
Academy of Romanian Scientists, Ilfov 3, 050044 Bucharest, Romania
Submission received: 3 February 2022 / Revised: 21 February 2022 / Accepted: 28 February 2022 / Published: 3 March 2022
(This article belongs to the Special Issue Recent Advances in Fractional Differential Equations and Inclusions)

Abstract

:
A coupled system of Hilfer fractional differential inclusions with nonlocal integral boundary conditions is considered. An existence result is established when the set-valued maps have non-convex values. We treat the case when the set-valued maps are Lipschitz in the state variables and we avoid the applications of fixed point theorems as usual. An illustration of the results is given by a suitable example.

1. Introduction

The fractional derivative introduced by Hilfer in [1] was recently used in the study of many boundary value problems concerned with fractional derivatives. This fractional derivative generalizes both Riemann–Liouville and Caputo derivatives; in fact, this derivative is an interpolation between Riemann–Liouville and Caputo derivatives. Several properties and applications of the Hilfer fractional derivative may be found in [2]. Additionally, we recall that the literature is full of explanations and motivations for considering systems defined by fractional order derivatives (e.g., [3,4,5,6,7] etc.).
Recently, many papers in the literature have been devoted to the study of fractional differential equations and inclusions defined by the Hilfer fractional derivative (e.g., [8,9,10,11] etc.). We point out that a complete survey on this field of study may be found in [8]. Taking into account this new trend in research, our intention is to contribute to the development of this topic by establishing new results for a particular class of problems.
The present note is devoted to the following boundary value problem
D H α 1 , β 1 x 1 ( t ) F 1 ( t , x 1 ( t ) , x 2 ( t ) ) , a . e . t [ a , b ] , D H α 2 , β 2 x 2 ( t ) F 2 ( t , x 1 ( t ) , x 2 ( t ) ) , a . e . t [ a , b ] ,
x 1 ( a ) = 0 , x 1 ( b ) = i = 1 m θ i I φ i x 2 ( ξ i ) , x 2 ( a ) = 0 , x 2 ( b ) = j = 1 n ζ j I ψ j x 1 ( z j ) ,
where F 1 ( . , . , . ) : [ a , b ] × R 2 P ( R ) , F 2 ( . , . , . ) : [ a , b ] × R 2 P ( R ) are given set-valued maps, D H α , β is the Hilfer fractional derivative of order α and type β , I φ is the Riemann–Liouville fractional integral of order φ > 0 , α 1 , α 2 ( 1 , 2 ) , β 1 , β 2 [ 0 , 1 ] , φ i , ψ j > 0 , ξ i , z j [ a , b ] , a > 0 , θ i , ζ j R , i = 1 , m ¯ , j = 1 , n ¯ .
Our study is motivated by two recent papers [10,11]. In [10], problems (1) and (2) is studied in the single-valued case; namely, the right-hand side in (1) is given by single-valued maps. Existence and uniqueness results are provided by using well-known fixed point theorems: Banach, Leray–Schauder and Krasnoselskii. In [11], a “simple” (not coupled) set-valued problem as in (1)–(2) is studied, and existence results are also obtained by applying known set-valued fixed point results as Leray–Schauder and Nadler.
The purpose of our note is twofold. On one hand, we extend the study in [10] to the set-valued framework, and on the other hand, we generalize the study in [11] to the coupled case. The approach we present here avoids the applications of fixed point theorems and takes into account the case when the values of F 1 and F 2 are not convex; but these set-valued maps are assumed to be Lipschitz in the second and third variable. In this case, we establish an existence result for problems (1) and (2). Our result use Filippov’s technique [12]; more exactly, the existence of solutions is obtained by starting from a pair of given “quasi”solutions. In addition, the result provides an estimate between the “quasi” solutions and the solutions obtained.
Even if the technique used here may be seen at other classes of coupled systems of fractional differential inclusions [13,14,15], to the best of our knowledge, the present paper is the first in literature which contains an existence result of Filippov type for coupled systems of Hilfer fractional differential inclusions.

2. Preliminaries

We set by I the interval [ a , b ] . We denote by C ( I , R ) the Banach space of all continuous functions x ( . ) : I R endowed with the norm | x ( . ) | C = sup t I | x ( t ) | and by L 1 ( I , R ) the Banach space of all integrable functions x ( . ) : I R endowed with the norm | x ( . ) | 1 = a b | x ( t ) | d t .
The Pompeiu–Hausdorff distance of the closed subsets A , B R is defined by d H ( A , B ) = max { d * ( A , B ) , d * ( B , A ) } , where d * ( A , B ) = sup { d ( a , B ) ; a A } and d ( x , B ) = inf y B d ( x , y ) .
The fractional integral of order α > 0 of a Lebesgue integrable function f : ( 0 , ) R is defined by
I α f ( t ) = 0 t ( t s ) α 1 Γ ( α ) f ( s ) d s ,
provided the right-hand side is pointwise defined on ( 0 , ) and Γ ( . ) is the (Euler’s) Gamma function defined by Γ ( α ) = 0 t α 1 e t d t .
The Riemann–Liouville fractional derivative of order α > 0 of a Lebesgue integrable function f : ( 0 , ) R is defined by
D α f ( t ) = 1 Γ ( n α ) d n d t n 0 t ( t s ) α + n 1 f ( s ) d s ,
where n = [ α ] + 1 , provided the right-hand side is pointwise defined on ( 0 , ) .
The Caputo fractional derivative of order α > 0 of a function f : [ 0 , ) R is defined by
D C α f ( t ) = 1 Γ ( n α ) 0 t ( t s ) α + n 1 f ( n ) ( s ) d s ,
where n = [ α ] + 1 . It is assumed implicitly that f is n times differentiable whose n-th derivative is absolutely continuous.
The Hilfer fractional derivative of order α ( n 1 , n ) and type β [ 0 , 1 ] of a function f : ( 0 , ) R is defined by
D H α , β f ( t ) = I β ( n α ) d n d t n I ( 1 β ) ( n α ) f ( t )
As it was already recalled, this derivative interpolates between Riemann–Liouville and Caputo derivatives. When β = 0 the Hilfer fractional derivative gives Riemann–Liouville fractional derivative D H α , 0 f ( t ) = d n d t n I n α f ( t ) and when β = 1 the Hilfer fractional derivative gives Caputo fractional derivative D H α , 1 f ( t ) = I n α d n d t n f ( t ) .
The next technical result proved in [10] considers a linear version of problems (1) and (2), for which an integral representation of the solution is provided.
Lemma 1.
Let f 1 ( . ) : [ a , b ] R , f 2 ( . ) : [ a , b ] R be continuous mappings and α 1 , α 2 ( 1 , 2 ) , β 1 , β 2 [ 0 , 1 ) , φ i , ψ j > 0 , γ 1 = α 1 + 2 β 1 α 1 β 1 , γ 2 = α 2 + 2 β 2 α 2 β 2 , ξ i , z j [ a , b ] , a > 0 , θ i , ζ j R , i = 1 , m ¯ , j = 1 , n ¯ .
Then, the solution of the system
D H α 1 , β 1 x 1 = f 1 ( t ) t [ a , b ] , D H α 2 , β 2 x 2 = f 2 ( t ) t [ a , b ]
with boundary conditions (2) is given by
x 1 ( t ) = I α 1 f 1 ( t ) + ( t a ) γ 1 1 C Γ ( γ 1 ) [ ( b a ) γ 2 1 Γ ( γ 2 ) ( i = 1 m θ i I α 2 + φ i f 2 ( ξ i ) I α 1 f 1 ( b ) ) + ( i = 1 m θ i ( ξ i a ) γ 2 + φ i 1 Γ ( γ 2 + φ i ) ) ( j = 1 n ζ j I α 1 + ψ j f 1 ( z j ) I α 2 f 2 ( b ) ) ] , x 2 ( t ) = I α 2 f 2 ( t ) + ( t a ) γ 2 1 C Γ ( γ 2 ) [ ( b a ) γ 1 1 Γ ( γ 1 ) ( j = 1 n ζ j I α 1 + ψ j f 1 ( z j ) I α 2 f 2 ( b ) ) + ( j = 1 n ζ j ( z j a ) γ 1 + ψ j 1 Γ ( γ 1 + ψ j ) ) ( i = 1 m θ i I α 2 + φ i f 2 ( ξ i ) I α 1 f 1 ( b ) ) ] ,
where C = ( b a ) γ 1 + γ 2 2 Γ ( γ 1 ) Γ ( γ 2 ) ( i = 1 m θ i ( ξ i a ) γ 2 + φ i 1 Γ ( γ 2 + φ i ) ) ( j = 1 n ζ j ( z j a ) γ 1 + ψ j 1 Γ ( γ 1 + ψ j ) ) 0 .
Definition 1.
The mappings x 1 ( . ) , x 2 ( . ) C ( I , R ) are said to be solutions of problems (1) and (2) if there exists f 1 ( . ) , f 2 ( . ) L 1 ( I , R ) such that f 1 ( t ) F 1 ( t , x 1 ( t ) , x 2 ( t ) ) a . e . ( I ) , f 2 ( t ) F 2 ( t , x 1 ( t ) , x 2 ( t ) ) a . e . ( I ) and x 1 ( . ) and x 2 ( . ) are given by (3).
In what follows, χ A ( · ) denotes the characteristic function of the set A R .
Remark 1.
Let us introduce the following notations
K 1 ( t , τ ) = ( t τ ) α 1 1 χ [ a , t ] ( τ ) Γ ( α 1 ) ( t a ) γ 1 1 ( b a ) γ 2 1 ( b τ ) α 1 1 C Γ ( γ 1 ) Γ ( γ 2 ) Γ ( α 1 ) + ( i = 1 m θ i ( ξ i a ) γ 2 + φ i 1 Γ ( γ 2 + φ i ) ) ( j = 1 n ζ j ( z j a ) α 1 + ψ j 1 Γ ( α 1 + ψ j ) χ [ a , z j ] ( τ ) ) ,
K 2 ( t , τ ) = ( t a ) γ 1 1 ( b a ) γ 2 1 C Γ ( γ 1 ) Γ ( γ 2 ) i = 1 m θ i ( ξ i τ ) α 2 + φ i 1 χ [ a , ξ i ] ( τ ) Γ ( α 2 + φ i ) i = 1 m θ i ( ξ i a ) γ 2 + φ i 1 ( b a ) α 2 1 Γ ( γ 2 + φ i ) Γ ( α 2 ) ,
K 3 ( t , τ ) = ( t a ) γ 2 1 ( b a ) γ 1 1 C Γ ( γ 1 ) Γ ( γ 2 ) j = 1 n ζ j ( z j τ ) α 1 + ψ j 1 χ [ a , z j ] ( τ ) Γ ( α 1 + ψ j ) j = 1 n ζ j ( z j a ) γ 1 + ψ j 1 ( b a ) α 1 1 Γ ( γ 1 + ψ j ) Γ ( α 1 ) ,
K 4 ( t , τ ) = ( t τ ) α 2 1 χ [ a , t ] ( τ ) Γ ( α 2 ) ( t a ) γ 2 1 ( b a ) γ 1 1 ( b τ ) α 2 1 C Γ ( γ 1 ) Γ ( γ 2 ) Γ ( α 2 ) + ( j = 1 n ζ j ( z j a ) γ 1 + ψ j 1 Γ ( γ 1 + ψ j ) ) ( i = 1 m θ i ( ξ i a ) α 2 + φ i 1 Γ ( α 2 + φ i ) χ [ a , ξ i ] ( τ ) ) .
Then, the solutions ( x 1 ( . ) , x 2 ( . ) ) in Lemma 1 may be put as
x 1 ( t ) = a b K 1 ( t , τ ) f 1 ( τ ) d τ + a b K 2 ( t , τ ) f 2 ( τ ) d τ , t I x 2 ( t ) = a b K 3 ( t , τ ) f 1 ( τ ) d τ + a b K 4 ( t , τ ) f 2 ( τ ) d τ , t I .
Moreover, we have the following estimates:
| K 1 ( t , τ ) | ( b a ) α 1 1 Γ ( α 1 ) + ( b a ) γ 1 + γ 2 + α 1 3 | C | Γ ( γ 1 ) Γ ( γ 2 ) Γ ( α 1 ) + i = 1 m | θ i | ( ξ i a ) γ 2 + φ i 1 Γ ( γ 2 + φ i ) · j = 1 n | ζ j | ( z j a ) α 1 + ψ j 1 Γ ( α 1 + ψ j ) = : k 1 t , τ I ,
| K 2 ( t , τ ) | ( b a ) γ 1 + γ 2 2 | C | Γ ( γ 1 ) Γ ( γ 2 ) i = 1 m | θ i | ( ξ i a ) α 2 + φ i 1 Γ ( α 2 + φ i ) + i = 1 m | θ i | ( ξ i a ) γ 2 + φ i 1 ( b a ) α 2 1 Γ ( γ 2 + φ i ) Γ ( α 2 ) = : k 2 t , τ I ,
| K 3 ( t , τ ) | ( b a ) γ 2 + γ 1 2 | C | Γ ( γ 1 ) Γ ( γ 2 ) j = 1 n | ζ j | ( z j a ) α 1 + ψ j 1 Γ ( α 1 + ψ j ) + j = 1 n | ζ j | ( z j a ) γ 1 + ψ j 1 ( b a ) α 1 1 Γ ( γ 1 + ψ j ) Γ ( α 1 ) = : k 3 t , τ I ,
| K 4 ( t , τ ) | ( b a ) α 2 1 Γ ( α 2 ) ( b a ) γ 1 + γ 2 + α 2 3 | C | Γ ( γ 1 ) Γ ( γ 2 ) Γ ( α 2 ) + j = 1 n | ζ j | ( z j a ) γ 1 + ψ j 1 Γ ( γ 1 + ψ j ) · i = 1 m | θ i | ( ξ i a ) α 2 + φ i 1 Γ ( α 2 + φ i ) = : k 4 t , τ I .
Finally, in the proof of our main result we need the following selection result for set-valued maps (e.g., [16]). It is, in fact, a variant of the well known selection theorem due to Kuratowski and Ryll-Nardzewski which, briefly, states that a measurable set-valued map with nonempty closed values admits a measurable selection.
Lemma 2.
Let Z be a separable Banach space, B its closed unit ball, A : I P ( Z ) is a set-valued map whose values are nonempty closed and b : I Z , c : I R + are two measurable functions. If
A ( t ) ( b ( t ) + c ( t ) B ) a . e . ( I ) ,
then the set-valued map t A ( t ) ( b ( t ) + c ( t ) B ) admits a measurable selection.

3. Main Result

Our results are proved under the following hypotheses.
Hypothesis 1.
(i) F 1 : I × R 2 P ( R ) and F 2 : I × R 2 P ( R ) have nonempty closed values and the set-valued maps F 1 ( . , y 1 , y 2 ) , F 2 ( . , y 1 , y 2 ) are measurable for any y 1 , y 2 R .
(ii) There exist l 1 ( . ) , l 2 ( . ) L 1 ( I , ( 0 , ) ) such that, for almost all t I , F 1 ( t , . , . ) is l 1 ( t ) -Lipschitz and F 2 ( t , . , . ) is l 2 ( t ) -Lipschitz; i.e.,
d H ( F 1 ( t , y 1 , z 1 ) , F 1 ( t , y 2 , z 2 ) ) l 1 ( t ) ( | y 1 y 2 | + | z 1 z 2 | ) y 1 , y 2 , z 1 , z 2 R .
d H ( F 2 ( t , y 1 , z 1 ) , F 2 ( t , y 2 , z 2 ) ) l 2 ( t ) ( | y 1 y 2 | + | z 1 z 2 | ) y 1 , y 2 , z 1 , z 2 R .
In what follows l ( t ) = k 1 l 1 ( t ) + k 2 l 2 ( t ) + k 3 l 1 ( t ) + k 4 l 2 ( t ) , t I .
Theorem 1.
Assume that C 0 , Hypothesis 1 is satisfied and | l ( . ) | 1 < 1 . ( y 1 ( . ) , y 2 ( . ) ) C ( I , R ) 2 are considered such that there exist q 1 ( . ) , q 2 ( . ) L 1 ( I , R ) that verify d ( D H α 1 , β 1 y 1 ( t ) , F 1 ( t , y 1 ( t ) , y 2 ( t ) ) ) q 1 ( t ) a.e. t I , d ( D H α 2 , β 2 y 2 ( t ) , F 2 ( t , y 1 ( t ) , y 2 ( t ) ) ) q 2 ( t ) a.e. t I , y 1 ( a ) = y 2 ( a ) = 0 , y 1 ( b ) = i = 1 m θ i I φ i y 2 ( ξ i ) and y 2 ( b ) = j = 1 n ζ j I ψ j y 1 ( z j ) .
Then, there exists ( x 1 ( . ) , x 2 ( . ) ) C ( I , R ) 2 a solution of problem (1) and (2) satisfying for all t I
| x 1 ( t ) y 1 ( t ) | + | x 2 ( t ) y 2 ( t ) | ( k 1 + k 3 ) | q 1 ( . ) | 1 + ( k 2 + k 4 ) | q 2 ( . ) | 1 1 | l ( . ) | 1 .
Proof. 
From the assumptions of the theorem
F 1 ( t , y 1 ( t ) , y 2 ( t ) ) { D H α 1 , β 1 y 1 ( t ) + q 1 ( t ) [ 1 , 1 ] } a . e . ( I ) , F 2 ( t , y 1 ( t ) , y 2 ( t ) ) { D H α 2 , β 2 y 2 ( t ) + q 2 ( t ) [ 1 , 1 ] } a . e . ( I ) .
By Lemma 2, there exist mesurable selections f 1 1 ( t ) F 1 ( t , y 1 ( t ) , y 2 ( t ) ) , f 2 1 ( t ) F 2 ( t , y 1 ( t ) , y 2 ( t ) )   a . e . ( I ) such that
| f 1 1 ( t ) D H α 1 , β 1 y 1 ( t ) | q 1 ( t ) , | f 2 1 ( t ) D H α 2 , β 2 y 2 ( t ) | q 2 ( t ) a . e . ( I ) .
Define
x 1 1 ( t ) = a b K 1 ( t , τ ) f 1 1 ( τ ) d τ + a b K 2 ( t , τ ) f 2 1 ( τ ) d τ , t I x 2 1 ( t ) = a b K 3 ( t , τ ) f 1 1 ( τ ) d τ + a b K 4 ( t , τ ) f 2 1 ( τ ) d τ , t I .
We have the estimates
| x 1 1 ( t ) y 1 ( t ) | k 1 | q 1 ( . ) | 1 + k 2 | q 2 ( . ) | 1 t I , | x 2 1 ( t ) y 2 ( t ) | k 3 | q 1 ( . ) | 1 + k 4 | q 2 ( . ) | 1 t I ,
and so,
| x 1 1 ( t ) y 1 ( t ) | + | x 2 1 ( t ) y 2 ( t ) | ( k 1 + k 3 ) | q 1 ( . ) | 1 + ( k 2 + k 4 ) | q 2 ( . ) | 1 = : k .
In the next part of the proof, we construct, by induction, the sequences x n 1 ( . ) , x n 2 ( . ) C ( I , R ) and f n 1 ( . ) , f n 2 ( . ) L 1 ( I , R ) , n 1 with the following properties
x 1 n ( t ) = a b K 1 ( t , τ ) f 1 n ( τ ) d τ + a b K 2 ( t , τ ) f 2 n ( τ ) d τ , t I x 2 n ( t ) = a b K 3 ( t , τ ) f 1 n ( τ ) d τ + a b K 4 ( t , τ ) f 2 n ( τ ) d τ , t I .
f 1 n ( t ) F 1 ( t , x 1 n 1 ( t ) , x 2 n 1 ( t ) ) , f 2 n ( t ) F 2 ( t , x 1 n 1 ( t ) , x 2 n 1 ( t ) ) a . e . ( I ) ,
| f 1 n + 1 ( t ) f 1 n ( t ) | l 1 ( t ) ( | x 1 n ( t ) x 1 n 1 ( t ) | + | x 2 n ( t ) x 2 n 1 ( t ) | ) a . e . ( I ) , | f 2 n + 1 ( t ) f 2 n ( t ) | l 2 ( t ) ( | x 1 n ( t ) x 1 n 1 ( t ) | + | x 2 n ( t ) x 2 n 1 ( t ) | ) a . e . ( I ) .
We point out that from (5) to (7) it follows
| x 1 n + 1 ( t ) x 1 n ( t ) | + | x 2 n + 1 ( t ) x 2 n ( t ) | k ( | l ( . ) | 1 ) n a . e . ( I ) n N .
The case n = 0 is already proved. Now, we assume (8) valid for n 1 . For almost all t I ,
| x 1 n + 1 ( t ) x 1 n ( t ) | a b | K 1 ( t , τ ) | . | f 1 n + 1 ( τ ) f 1 n ( τ ) | d τ + a b | K 2 ( t , τ ) | . | f 2 n + 1 ( τ ) f 2 n ( τ ) | d τ k 1 a b | f 1 n + 1 ( τ ) f 1 n ( τ ) | d τ + k 2 a b | f 2 n + 1 ( τ ) f 2 n ( τ ) | d τ k 1 a b l 1 ( τ ) ( | x 1 n ( τ ) x 1 n 1 ( τ ) | + | x 2 n ( τ ) x 2 n 1 ( τ ) | ) d τ + k 2 a b l 2 ( τ ) ( | x 1 n ( τ ) x 1 n 1 ( τ ) | + | x 2 n ( τ ) x 2 n 1 ( τ ) | ) d τ k ( | l ( . ) | 1 ) n 1 ( k 1 a b l 1 ( τ ) d τ + k 2 a b l 2 ( τ ) d τ ) .
In a similar way, we obtain for almost all t I ,
| x 2 n + 1 ( t ) x 2 n ( t ) | k ( | l ( . ) | 1 ) n 1 ( k 3 a b l 1 ( τ ) d τ + k 4 a b l 2 ( τ ) d τ ) .
Therefore, (8) is true for n.
Inequality (8) shows that the sequences { x 1 n ( . ) } , { x 2 n ( . ) } are Cauchy in the space C ( I , R ) . Let x 1 ( . ) C ( I , R ) and x 2 ( . ) C ( I , R ) be their limits in C ( I , R ) . Additionally, from (7) we deduce that, for almost all t I , the sequences { f 1 n ( t ) } , { f 2 n ( t ) } are Cauchy in R . We consider f 1 ( . ) , f 2 ( . ) their pointwise limit.
At the same time, inequality (8) and Hypothesis 1 give
| x 1 n ( t ) y 1 ( t ) | + | x 2 n ( t ) y 2 ( t ) | | x 1 1 ( t ) y 1 ( t ) | + | x 2 1 ( t ) y 2 ( t ) | + i = 1 n 1 ( | x 1 i + 1 ( t ) x 1 i ( t ) | + | x 2 i + 1 ( t ) x 2 i ( t ) | ) k + i = 1 n k ( | l ( . ) | 1 ) i k 1 | l ( . ) | 1 .
and
| f 1 n ( t ) D H α 1 , β 1 y 1 ( t ) | + | f 2 n ( t ) D H α 2 , β 2 y 2 ( t ) | | f 1 1 ( t ) D H α 1 , β 1 y 1 ( t ) | + | f 2 1 ( t ) D H α 2 , β 2 y 2 ( t ) | + i = 1 n 1 ( | f 1 i + 1 ( t ) f 1 i ( t ) | + | f 2 i + 1 ( t ) f 2 i ( t ) | ) | f 1 1 ( t ) D H α 1 , β 1 y 1 ( t ) | + | f 2 1 ( t ) D H α 2 , β 2 y 2 ( t ) | + i = 1 n 1 ( l 1 ( t ) + l 2 ( t ) ) ( | x 1 i ( t ) x 1 i 1 ( t ) | + | x 2 i ( t ) x 2 i 1 ( t ) | ) q 1 ( t ) + q 2 ( t ) + ( l 1 ( t ) + l 2 ( t ) ) k 1 | l ( . ) | 1
for almost all t I .
This means that the sequences f 1 n ( . ) , f 2 n ( . ) are integrably bounded and, therefore, their limits f 1 ( . ) , f 2 ( . ) belong to L 1 ( I , R ) .
The next step of the proof contains the construction in (5)–(7). By induction, we suppose that for M 1 , x 1 m ( . ) , x 2 m ( . ) C ( I , R ) and f 1 m ( . ) , f 2 m ( . ) L 1 ( I , R ) , m = 1 , 2 , . . . M with (5) and (7) for m = 1 , 2 , . . . M and (6) for m = 1 , 2 , . . . M 1 are constructed.
Using again Hypothesis 1
F 1 ( t , x 1 M ( t ) , x 2 M ( t ) ) { f 1 M ( t ) + ( l 1 ( t ) | x 1 M ( t ) x 1 M 1 ( t ) | + l 1 ( t ) | x 2 M ( t ) x 2 M 1 ( t ) | ) [ 1 , 1 ] } , F 2 ( t , x 1 M ( t ) , x 2 M ( t ) ) { f 2 M ( t ) + ( l 2 ( t ) | x 1 M ( t ) x 1 M 1 ( t ) | + l 2 ( t ) | x 2 M ( t ) x 2 M 1 ( t ) | ) [ 1 , 1 ] }
for almost all t I . By Lemma 2, we obtain the existence of measurable selections f 1 M + 1 ( . ) of F 1 ( . , x M 1 ( . ) , x M 2 ( . ) ) and f 2 M + 1 ( . ) of F 2 ( . , x M 1 ( . ) , x M 2 ( . ) ) such that
| f 1 M + 1 ( t ) f 1 M ( t ) | l 1 ( t ) ( | x 1 M ( t ) x 1 M 1 ( t ) | + | x 2 M ( t ) x 2 M 1 ( t ) | ) a . e . ( I ) , | f 2 M + 1 ( t ) f 2 M ( t ) | l 2 ( t ) ( | x 1 M ( t ) x 1 M 1 ( t ) | + | x 2 M ( t ) x 2 M 1 ( t ) | ) a . e . ( I ) .
We define x 1 M + 1 ( . ) , x 2 M + 1 ( . ) as in (5) with n = M + 1 .
Finally, it remains to take n in (5) and (9) in order to finish the proof. □
Theorem 1 is the first in the literature that contains an existence result of Filippov type for coupled systems of Hilfer fractional differential inclusions. Due to the presence in its statement of the ”quasi” solutions, the formulation of Theorem 1 seems to look complicated. However, for a particular choice of the ”quasi” solutions (namely, y 1 ( . ) = y 2 ( . ) = 0 ), one may obtain a statement similar to a result that can be derived by using the set-valued contration principle. This may be seen in the following consequence of Theorem 1.
Corollary 1.
Assume that C 0 , Hypothesis 1 is satisfied, | l ( . ) | 1 < 1 , d ( 0 , F 1 ( t , 0 , 0 ) ) l 1 ( t ) a.e. t I and d ( 0 , F 2 ( t , 0 , 0 ) ) l 2 ( t ) a.e. t I .
Then, there exists ( x 1 ( . ) , x 2 ( . ) ) C ( I , R ) 2 a solution of problems (1) and (2) satisfying for all t I
| x 1 ( t ) | + | x 2 ( t ) | ( k 1 + k 3 ) | l 1 ( . ) | 1 + ( k 2 + k 4 ) | l 2 ( . ) | 1 1 | l ( . ) | 1 .
Proof. 
We apply Theorem 1 with y 1 ( . ) = y 2 ( . ) = 0 , q 1 ( . ) = l 1 ( . ) and q 2 ( . ) = l 2 ( . ) . □
Remark 2.
If in (1), F 1 and F 2 are single-valued maps, Corollary 1 provides a generalization to the set-valued framework of Theorem 1 in [10].
Example 1.
Let us consider the problem
D H 3 2 , 1 2 x 1 ( t ) [ 1 300 | x 1 ( t ) | 1 + | x 1 ( t ) | , 0 ] [ 0 , 1 300 | x 2 ( t ) | 1 + | x 2 ( t ) | ] a . e . [ 1 3 , 10 3 ] , D H 5 4 , 2 3 x 2 ( t ) [ 1 300 | sin ( x 1 ( t ) ) | 1 + | sin ( x 1 ( t ) ) | , 0 ] [ 0 , 1 300 | cos ( x 2 ( t ) ) | 1 + | cos ( x 2 ( t ) ) | ] a . e . [ 1 3 , 10 3 ] ,
with nonlocal integral boundary conditions as in [10]
x 1 ( 1 3 ) = 0 , x 1 ( 10 3 ) = 1 2 I 1 3 x 2 ( 2 3 ) + 2 3 I 1 2 x 2 ( 1 ) + 3 4 I 3 5 x 2 ( 4 3 ) + 4 5 I 2 3 x 2 ( 5 3 ) x 2 ( 1 3 ) = 0 , x 2 ( 10 3 ) = 3 2 I 3 2 x 1 ( 2 ) + 4 3 I 5 3 x 1 ( 7 3 ) + 5 4 I 7 4 x 1 ( 8 3 ) + 6 5 I 9 5 x 1 ( 3 ) .
Thus, F 1 ( t , ( x 1 , x 2 ) ) = [ 1 300 | x 1 | 1 + | x 1 | , 0 ] [ 0 , 1 300 | x 2 | 1 + | x 2 | ] , F 2 ( t , ( x 1 , x 2 ) ) = [ 1 300 | sin x 1 | 1 + | sin x 1 | , 0 ] [ 0 , 1 300 | cos x 2 | 1 + | cos x 2 | ] , α 1 = 3 2 , β 1 = 1 2 , γ 1 = 7 4 , α 2 = 5 4 , β 2 = 2 3 , γ 2 = 41 21 , a = 1 3 , b = 10 3 , m = n = 4 , θ i = i i + 1 , ζ i = 2 + i i + 1 , φ i = i i + 2 , ψ i = 2 i + 1 i + 1 , ξ i = i + 1 3 , z i = i + 5 3 , i = 1 , 2 , 3 , 4 and C = 6 , 371398411 .
For all t [ 1 3 , 10 3 ] and all x 1 , x 2 , y 1 , y 2 R we have
sup { | z | ; z F i ( t , ( x 1 , x 2 ) ) } 1 300 , i = 1 , 2 ,
d H ( F 1 ( t , ( x 1 , x 2 ) ) , F 2 ( t , ( y 1 , y 2 ) ) 1 300 | x 1 y 1 | + 1 300 | x 2 y 2 | i = 1 , 2 .
By standard computations (e.g., [10]) k 1 12.6 , k 2 3.26 , k 3 48.73 , k 4 24.1 and ( k 1 + k 3 ) 3 1 300 + ( k 2 + k 4 ) 3 1 300 0.89 < 1 . So, we may apply Corollary 1 and obtain the existence of a solution for problems (10) and (11).

4. Conclusions

In the present paper, we extended the research in [10] to multivalued problems and the research in [11] to the situation of coupled Hilfer fractional differential inclusions. We established an existence result for problems (1) and (2) when the set-valued maps are Lipschitz in the state variables without any assumptions concerning the convexity of the values of the set-valued maps. Our approach uses a technique due to Filippov ([12]) instead of an usual application of set-valued fixed point theorems. An illustration of our result is provided by a numerical example.
It is worth mentioning that Theorem 1 may be a basic tool in the study of optimal control problems defined by such kinds of coupled systems of Hilfer fractional differential inclusions.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  2. Hilfer, R.; Luchko, Y.; Tomovski, Z. Operational method for the solution of fractional differential equations with generalized Riemann-Liouville fractional derivative. Fract. Calc. Appl. Anal. 2009, 12, 299–318. [Google Scholar]
  3. Băleanu, D.; Diethelm, K.; Scalas, E.; Trujillo, J.J. Fractional Calculus Models and Numerical Methods; World Scientific: Singapore, 2012. [Google Scholar]
  4. Diethelm, K. The Analysis of Fractional Differential Equations; Springer: Berlin, Germany, 2010. [Google Scholar]
  5. Kilbas, A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  6. Miller, K.; Ross, B. An Introduction to the Fractional Calculus and Differential Equations; John Wiley: New York, NY, USA, 1993. [Google Scholar]
  7. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  8. Ntouyas, S.K. A survey on existence results for boundary value problems of Hilfer fractional derivative equations and inclusions. Foundations 2021, 1, 63–98. [Google Scholar] [CrossRef]
  9. Phuangthong, N.; Ntouyas, S.K.; Tariboon, J.; Nonlaopon, K. Nonlocal sequential boundary value problems fir Hilfer type fractional integro-differential equations and inclusions. Mathematics 2021, 9, 615. [Google Scholar] [CrossRef]
  10. Wongcharoen, A.; Ntouyas, S.K.; Tariboon, J. On coupled systems for Hilfer fractional differential equations with nonlocal integral boundary conditions. J. Math. 2020, 2020, 2875152. [Google Scholar] [CrossRef]
  11. Wongcharoen, A.; Ntouyas, S.K.; Tariboon, J. Boundary value problems for Hilfer fractional differential inclusions with nonlocal integral boundary conditions. Mathematics 2020, 8, 1905. [Google Scholar] [CrossRef]
  12. Filippov, A.F. Classical solutions of differential equations with multivalued right hand side. SIAM J. Control 1967, 5, 609–621. [Google Scholar] [CrossRef]
  13. Cernea, A. Existence of solutions for some coupled systems of fractional differential inclusions. Mathematics 2020, 8, 700. [Google Scholar] [CrossRef]
  14. Cernea, A. On some coupled systems of fractional differential inclusions. Fract. Differ. Calc. 2021, 11, 133–145. [Google Scholar] [CrossRef]
  15. Cernea, A. A note on a coupled system of Caputo-Fabrizio fractional differential inclusions. Ann. Commun. Math. 2021, 4, 190–195. [Google Scholar]
  16. Aubin, J.P.; Frankowska, H. Set-Valued Analysis; Birkhauser: Basel, Switzerland, 1990. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cernea, A. A Note on a Coupled System of Hilfer Fractional Differential Inclusions. Foundations 2022, 2, 290-297. https://0-doi-org.brum.beds.ac.uk/10.3390/foundations2010020

AMA Style

Cernea A. A Note on a Coupled System of Hilfer Fractional Differential Inclusions. Foundations. 2022; 2(1):290-297. https://0-doi-org.brum.beds.ac.uk/10.3390/foundations2010020

Chicago/Turabian Style

Cernea, Aurelian. 2022. "A Note on a Coupled System of Hilfer Fractional Differential Inclusions" Foundations 2, no. 1: 290-297. https://0-doi-org.brum.beds.ac.uk/10.3390/foundations2010020

Article Metrics

Back to TopTop