Next Article in Journal / Special Issue
Unique Existence Result of Approximate Solution to Initial Value Problem for Fractional Differential Equation of Variable Order Involving the Derivative Arguments on the Half-Axis
Previous Article in Journal
Means as Improper Integrals
Previous Article in Special Issue
On the Solvability of a Mixed Problem for a High-Order Partial Differential Equation with Fractional Derivatives with Respect to Time, with Laplace Operators with Spatial Variables and Nonlocal Boundary Conditions in Sobolev Classes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Random Coupled Hilfer and Hadamard Fractional Differential Systems in Generalized Banach Spaces

1
Laboratory of Mathematics, Geometry, Analysis, Control and Applications, Tahar Moulay University of Saïda, P.O. Box 138, EN-Nasr, Saïda 20000, Algeria
2
College of Science, Geology and Geophysics Department, King Saud University, Riyadh 11451, Saudi Arabia
3
Laboratory of Mathematics, Djillali Liabes University of Sidi Bel-Abbès, P.O. Box 89, Sidi Bel-Abbès 22000, Algeria
4
Department of Mathematics, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
5
Faculty of Mathematics and Computational Science, Xiangtan University, Xiangtan 411105, China
6
Nonlinear Analysis and Applied Mathematics (NAAM) Research Group, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Submission received: 30 January 2019 / Revised: 12 March 2019 / Accepted: 15 March 2019 / Published: 20 March 2019

Abstract

:
This article deals with some existence and uniqueness result of random solutions for some coupled systems of Hilfer and Hilfer–Hadamard fractional differential equations with random effects. Some applications are made of generalizations of classical random fixed point theorems on generalized Banach spaces.

1. Introduction

Fractional calculus is an extension of the ordinary differentiation and integration to arbitrary non-integer order. In recent years, this theory has become an important object of investigations due to its demonstrated applications in different areas of physics and engineering (see, for example, [1,2] and the references therein). In particular, time fractional differential equations are used when attempting to describe transport processes with long memory. Recently, the study of time fractional ordinary and partial differential equations has received great attention from many researchers, both in theory and in applications; we refer the reader to the monographs of Abbas et al. [3,4,5], Samko et al. [6], and Kilbas et al. [7], and the papers [8,9,10,11,12,13,14] and the references therein. On the other hand, the existence of solutions of initial and boundary value problems for fractional differential equations with the Hilfer fractional derivative have started to draw attention. For the related works, see for example [1,15,16,17,18,19,20] and the references therein.
Functional differential equations with random effects are differential equations with a stochastic process in their vector field [21,22,23,24,25]. They play a fundamental role in the theory of random dynamical systems.
Consider the following coupled system of Hilfer fractional differential equations:
( D 0 α 1 , β 1 u ) ( t , w ) = f 1 ( t , u ( t , w ) , v ( t , w ) , w ) ( D 0 α 2 , β 2 v ) ( t , w ) = f 2 ( t , u ( t , w ) , v ( t , w ) , w ) ; t I : = [ 0 , T ] , w Ω ,
with the following initial conditions:
( I 0 1 γ 1 u ) ( 0 , w ) = ϕ 1 ( w ) ( I 0 1 γ 2 v ) ( 0 , w ) = ϕ 2 ( w ) ; w Ω ,
where T > 0 , α i ( 0 , 1 ) , β i [ 0 , 1 ] , ( Ω , A ) is a measurable space, γ i = α i + β i α i β i , ϕ i : Ω R m , f i : I × R m × R m × Ω R m ; i = 1 , 2 , are given functions, I 0 1 γ i is the left-sided mixed Riemann–Liouville integral of order 1 γ i , and D 0 α i , β i is the generalized Riemann–Liouville derivative (Hilfer) operator of order α i and type β i : i = 1 , 2 . Next, we discuss the following coupled system of Hilfer–Hadamard fractional differential equations:
( H D 1 α 1 , β 1 u ) ( t , w ) = g 1 ( t , u ( t , w ) , v ( t , w ) , w ) ( H D 1 α 2 , β 2 v ) ( t , w ) = g 2 ( t , u ( t , w ) , v ( t , w ) , w ) ; t [ 1 , T ] , w Ω ,
with the following initial conditions:
( H I 1 1 γ 1 u ) ( 1 , w ) = ψ 1 ( w ) ( H I 1 1 γ 2 v ) ( 1 , w ) = ψ 2 ( w ) ; w Ω ,
where T > 1 , α i ( 0 , 1 ) , β i [ 0 , 1 ] , γ i = α i + β i α i β i , ψ i : Ω R m , g i : [ 1 , T ] × R m × R m × Ω R m ; i = 1 , 2 are given functions, R m ; m N * , H I 1 1 γ i is the left-sided mixed Hadamard integral of order 1 γ i , and H D 1 α i , β i is the Hilfer–Hadamard fractional derivative of order α i and type β i ; i = 1 , 2 .

2. Preliminaries

We denote by C ; the Banach space of all continuous functions from I into R m with the supremum (uniform) norm · . As usual, A C ( I ) denotes the space of absolutely continuous functions from I into R m . By L 1 ( I ) , we denote the space of Lebesgue-integrable functions v : I R m with the norm:
v 1 = 0 T v ( t ) d t .
By C γ ( I ) and C γ 1 ( I ) , we denote the weighted spaces of continuous functions defined by:
C γ ( I ) = { w : ( 0 , T ] R m : t 1 γ w ( t ) C } ,
with the norm:
w C γ : = sup t I t 1 γ w ( t ) ,
and:
C γ 1 ( I ) = { w C : d w d t C γ } ,
with the norm:
w C γ 1 : = w + w C γ .
Furthermore, by C : = C γ 1 × C γ 2 , we denote the product weighted space with the norm:
( u , v ) C = u C γ 1 + v C γ 2 .
Now, we give some definitions and properties of fractional calculus.
Definition 1.
[4,6,7] The left-sided mixed Riemann–Liouville integral of order r > 0 of a function w L 1 ( I ) is defined by:
( I 0 r w ) ( t ) = 1 Γ ( r ) 0 t ( t s ) r 1 w ( s ) d s ; f o r a . e . t I ,
where Γ ( · ) is the (Euler’s) Gamma function.
Notice that for all r , r 1 , r 2 > 0 and each w C , we have I 0 r w C , and:
( I 0 r 1 I 0 r 2 w ) ( t ) = ( I 0 r 1 + r 2 w ) ( t ) ; f o r a . e . t I .
Definition 2.
[4,6,7] The Riemann–Liouville fractional derivative of order r ( 0 , 1 ] of a function w L 1 ( I ) is defined by:
( D 0 r w ) ( t ) = d d t I 0 1 r w ( t ) = 1 Γ ( 1 r ) d d t 0 t ( t s ) r w ( s ) d s ; f o r a . e . t I .
Let r ( 0 , 1 ] , γ [ 0 , 1 ) and w C 1 γ ( I ) . Then, the following expression leads to the left inverse operator as follows.
( D 0 r I 0 r w ) ( t ) = w ( t ) ; f o r a l l t ( 0 , T ] .
Moreover, if I 0 1 r w C 1 γ 1 ( I ) , then the following composition is proven in [6]:
( I 0 r D 0 r w ) ( t ) = w ( t ) ( I 0 1 r w ) ( 0 + ) Γ ( r ) t r 1 ; f o r a l l t ( 0 , T ] .
Definition 3.
[4,6,7] The Caputo fractional derivative of order r ( 0 , 1 ] of a function w L 1 ( I ) is defined by:
( c D 0 r w ) ( t ) = I 0 1 r d d t w ( t ) = 1 Γ ( 1 r ) 0 t ( t s ) r d d s w ( s ) d s ; f o r a . e . t I .
In [1], R.Hilfer studied applications of a generalized fractional operator having the Riemann–Liouville and the Caputo derivatives as specific cases (see also [17,19]).
Definition 4.
(Hilfer derivative). Let α ( 0 , 1 ) , β [ 0 , 1 ] , w L 1 ( I ) , and I 0 ( 1 α ) ( 1 β ) w A C ( I ) . The Hilfer fractional derivative of order α and type β of w is defined as:
( D 0 α , β w ) ( t ) = I 0 β ( 1 α ) d d t I 0 ( 1 α ) ( 1 β ) w ( t ) ; f o r a . e . t I .
Property 1.
Let α ( 0 , 1 ) , β [ 0 , 1 ] , γ = α + β α β , and w L 1 ( I ) .
  • The operator ( D 0 α , β w ) ( t ) can be written as:
    ( D 0 α , β w ) ( t ) = I 0 β ( 1 α ) d d t I 0 1 γ w ( t ) = I 0 β ( 1 α ) D 0 γ w ( t ) ; f o r a . e . t I .
    Moreover, the parameter γ satisfies:
    γ ( 0 , 1 ] , γ α , γ > β , 1 γ < 1 β ( 1 α ) .
  • The generalization (5) for β = 0 , coincides with the Riemann–Liouville derivative and for β = 1 with the Caputo derivative.
    D 0 α , 0 = D 0 α , a n d D 0 α , 1 = c D 0 α .
  • If D 0 β ( 1 α ) w exists and is in L 1 ( I ) , then:
    ( D 0 α , β I 0 α w ) ( t ) = ( I 0 β ( 1 α ) D 0 β ( 1 α ) w ) ( t ) ; f o r a . e . t I .
    Furthermore, if w C γ ( I ) and I 0 1 β ( 1 α ) w C γ 1 ( I ) , then:
    ( D 0 α , β I 0 α w ) ( t ) = w ( t ) ; f o r a . e . t I .
  • If D 0 γ w exists and is in L 1 ( I ) , then:
    ( I 0 α D 0 α , β w ) ( t ) = ( I 0 γ D 0 γ w ) ( t ) = w ( t ) I 0 1 γ ( 0 + ) Γ ( γ ) t γ 1 ; f o r a . e . t I .
Corollary 1.
Let h C γ ( I ) . Then, the Cauchy problem:
( D 0 α , β u ) ( t ) = h ( t ) ; t I , ( I 0 1 γ u ) ( t ) | t = 0 = ϕ ,
has the following unique solution:
u ( t ) = ϕ Γ ( γ ) t γ 1 + ( I 0 α h ) ( t ) .
Let β R m be the σ -algebra of Borel subsets of R m . A mapping v : Ω R m is said to be measurable if for any B β R m ; one has:
v 1 ( B ) = { w Ω : v ( w ) B } A .
Definition 5.
Let A × β R m be the direct product of the σ-algebras A and β R m those defined in Ω and R m , respectively. A mapping T : Ω × R m R m is called jointly measurable if for any B β R m , one has:
T 1 ( B ) = { ( w , v ) Ω × E : T ( w , v ) B } A × β R m .
Definition 6.
A function T : Ω × R m R m is called jointly measurable if T ( · , u ) is measurable for all u R m and T ( w , · ) is continuous for all: w Ω .
A random operator is a mapping T : Ω × R m R m such that T ( w , u ) is measurable in w for all u R m , and it expressed as T ( w ) u = T ( w , u ) ; we also say that T ( w ) is a random operator on R m . The random operator T ( w ) on E is called continuous (resp. compact, totally bounded, and completely continuous) if T ( w , u ) is continuous (resp. compact, totally bounded, and completely continuous) in u for all w Ω . The details of completely continuous random operators in Banach spaces and their properties appear in Itoh [26].
Definition 7.
[27] Let P ( Y ) be the family of all nonempty subsets of Y and C be a mapping from Ω into P ( Y ) . A mapping T : { ( w , y ) : w Ω , y C ( w ) } Y is called a random operator with stochastic domain C if C is measurable (i.e., for all closed A Y , { w Ω , C ( w ) A } is measurable), and for all open D Y and all y Y , { w Ω : y C ( w ) , T ( w , y ) D } is measurable. T will be called continuous if every T ( w ) is continuous. For a random operator T , a mapping y : Ω Y is called a random (stochastic) fixed point of T if for P-almost all w Ω , y ( w ) C ( w ) and T ( w ) y ( w ) = y ( w ) , and for all open D Y , { w Ω : y ( w ) D } is measurable.
Definition 8.
A function f : I × R m × Ω R m is called random Carathéodory if the following conditions are satisfied:
(i) 
The map ( t , w ) f ( x , y , u , w ) is jointly measurable for all u R m and
(ii) 
The map u f ( t , u , w ) is continuous for all t I and w Ω .
Let x , , y R m with x = ( x 1 , x 2 , , x m ) , y = ( y 1 , y 2 , , y m ) .
By x y , we mean x i y i ; i = 1 , , m . Also | x | = ( | x 1 | , | x 2 | , , | x m | ) , max ( x , y ) = ( max ( x 1 , y 1 ) , max ( x 2 , y 2 ) , , max ( x m , y m ) ) , and R + m = { x R m : x i R + , i = 1 , , m } . If c R , then x c means x i c ; i = 1 , , m .
Definition 9.
Let X be a nonempty set. By a vector-valued metric on X, we mean a map d : X × X R m with the following properties:
(i) 
d ( x , y ) 0 for all x , y X , and if d ( x , y ) = 0 , then x = y ;
(ii) 
d ( x , y ) = d ( y , x ) for all x , y X ;
(iii) 
d ( x , z ) d ( x , y ) + d ( y , z ) for all x , y , z X .
We call the pair ( X , d ) a generalized metric space with d ( x , y ) : = d 1 ( x , y ) d 2 ( x , y ) · · · d m ( x , y ) .
Notice that d is a generalized metric space on X if and only if d i ; i = 1 , , m are metrics on X . For r = ( r 1 , , r m ) R m and x 0 X , we will denote by:
B r ( x 0 ) : = { x X : d ( x 0 , x ) < r } = { x X : d i ( x 0 , x ) < r i ; i = 1 , , m }
the open ball centered in x 0 with radius r and:
B ¯ r ( x 0 ) : = { x X : d ( x 0 , x ) r } = { x X : d i ( x 0 , x ) r i ; i = 1 , , m }
the closed ball centered in x 0 with radius r . We mention that for generalized metric spaces, the notations of open, closed, compact, convex sets, convergence, and Cauchy sequence are similar to those in usual metric spaces.
Definition 10.
[28,29] A square matrix of real numbers is said to be convergent to zero if and only if its spectral radius ρ ( M ) is strictly less than one. In other words, this means that all the eigenvalues of M are in the open unit disc, i.e., | λ | < 1 ; for every λ C with d e t ( M λ I ) = 0 ; where I denotes the unit matrix of M m × m ( R ) .
Example 1.
The matrix A M 2 × 2 ( R ) defined by:
A = a b c d ,
converges to zero in the following cases:
(1) 
b = c = 0 , a , d > 0 , and max { a , d } < 1 .
(2) 
c = 0 , a , d > 0 , a + d < 1 , and 1 < b < 0 .
(3) 
a + b = c + d = 0 , a > 1 , c > 0 , and | a c | < 1 .
In the sequel, we will make use of the following random fixed point theorems:
Theorem 1.
[23,24,25] Let ( Ω , F ) be a measurable space, X a real separable generalized Banach space, and F : Ω × X X a continuous random operator, and let M ( w ) M n × n ( R + ) be a random variable matrix such that for every w Ω , the matrix M ( w ) converges to zero and:
d ( F ( w , x 1 ) , F ( w , x 2 ) ) M ( w ) d ( x 1 , x 2 ) ; f o r e a c h x 1 , x 2 X a n d w Ω ,
then there exists a random variable x : Ω X that is the unique random fixed point of F .
Theorem 2.
[23,24,25] Let ( Ω , F ) be a measurable space, X be a real separable generalized Banach space, and F : Ω × X X be a completely continuous random operator. Then, either:
(i) 
the random equation F ( w , x ) = x has a random solution, i.e., there is a measurable function x : Ω X such that F ( w , x ( w ) ) = x ( w ) for all w Ω or
(ii) 
the set M = { x : Ω X i s m e a s u r a b l e : λ ( w ) F ( w , x ) = x } is unbounded for some measurable function λ : Ω X with 0 < λ ( w ) < 1 on Ω .
Furthermore, we will use the following Gronwall lemma:
Lemma 1.
[23] Let u : I [ 0 , ) be a real function and u ( · ) a nonnegative, locally-integrable function on I . Assume that there exist constants c > 0 and r < 1 such that:
u ( t ) v ( t ) + c 0 t u ( s ) ( t s ) r d s ,
then, there exists a constant K : = K ( r ) such that:
u ( t ) v ( t ) + c K 0 t v ( s ) ( t s ) r d s ,
for every t I .

3. Coupled Hilfer Fractional Differential Systems

In this section, we are concerned with the existence and uniqueness results of the system (1) and (2).
Definition 11.
By a solution of the problem (1) and (2), we mean coupled measurable functions ( u , v ) C γ 1 × C γ 2 , which satisfy the Equation (1) on I , and the conditions ( I 0 1 γ 1 u ) ( 0 + ) = ϕ 1 , and ( I 0 1 γ 2 v ) ( 0 + ) = ϕ 2 .
The following hypotheses will be used in the sequel.
( H 1 )
The functions f i ; i = 1 , 2 are Carathéodory.
( H 2 )
There exist measurable functions p i , q i : Ω ( 0 , ) ; i = 1 , 2 such that:
f i ( t , u 1 , v 1 ) f i ( t , u 2 , v 2 ) p i ( w ) u 1 u 2 + q i ( w ) v 1 v 2 ;
f o r a . e . t I , a n d e a c h u i , v i R m , i = 1 , 2 .
( H 3 )
There exist measurable functions a i , b i : Ω ( 0 , ) ; i = 1 , 2 such that:
f i ( t , u , v ) a i ( w ) u + b i ( w ) v ; f o r a . e . t I , a n d e a c h u , v R m .
First, we prove an existence and uniqueness result for the coupled system (1)–(2) by using Banach’s random fixed point theorem in generalized Banach spaces.
Theorem 3.
Assume that the hypotheses ( H 1 ) and ( H 2 ) hold. If for every w Ω , the matrix:
M ( w ) : = T α 1 Γ ( 1 + α 1 ) p 1 ( w ) T α 1 Γ ( 1 + α 1 ) q 1 ( w ) T α 2 Γ ( 1 + α 2 ) p 2 ( w ) T α 2 Γ ( 1 + α 2 ) q 2 ( w )
converges to zero, then the coupled system (1) and (2) has a unique random solution.
Proof. 
Define the operators N 1 : C × Ω C γ 1 and N 2 : C × Ω C γ 2 by:
( N 1 ( u , v ) ) ( t , w ) = ϕ 1 ( w ) Γ ( γ 1 ) t γ 1 1 + 0 t ( t s ) α 1 1 f ( s , u ( s , w ) , v ( s , w ) , w ) Γ ( α 1 ) d s ,
and:
( N 2 ( u , v ) ) ( t , w ) = ϕ 2 ( w ) Γ ( γ 2 ) t γ 2 1 + 0 t ( t s ) α 2 1 f ( s , u ( s , w ) , v ( s , w ) , w ) Γ ( α 2 ) d s .
Consider the operator N : C × Ω C defined by:
( N ( u , v ) ) ( t , w ) = ( ( N 1 ( u , v ) ) ( t , w ) , ( N 2 ( u , v ) ) ( t , w ) ) .
Clearly, the fixed points of the operator N are random solutions of the system (1) and (2).
Let us show that N is a random operator on C . Since f i ; i = 1 , 2 are Carathéodory functions, then w f i ( t , u , v , w ) are measurable maps. We concluded that the maps:
w ( N 1 ( u , v ) ) ( t , w ) a n d w ( N 2 ( u , v ) ) ( t , w ) ,
are measurable. As a result, N is a random operator on C × Ω into C . We show that N satisfies all conditions of Theorem 1.
For any w Ω and each ( u 1 , v 1 ) , ( u 2 , v 2 ) C , and t I , we have:
t 1 γ 1 ( N 1 ( u 1 , v 1 ) ) ( t , w ) t 1 γ 1 ( N 1 ( u 2 , v 2 ) ) ( t , w ) t 1 γ 1 Γ ( α 1 ) 0 t ( t s ) α 1 1 f 1 ( s , u 1 ( s , w ) , v 1 ( s , w ) , w ) f 1 ( s , u 2 ( s , w ) , v 2 ( s , w ) , w ) d s t 1 γ 1 Γ ( α 1 ) 0 t ( t s ) α 1 1 ( p 1 ( w ) u 1 ( s , w ) v 1 ( s , w ) + q 1 ( w ) u 2 ( s , w ) v 2 ( s , w ) ) d s 1 Γ ( α 1 ) 0 t ( t s ) α 1 1 ( p 1 ( w ) s 1 γ 1 u 1 ( s , w ) v 1 ( s , w ) + q 1 ( w ) s 1 γ 1 u 2 ( s , w ) v 2 ( s , w ) ) d s p 1 ( w ) u 1 ( · , w ) v 1 ( · , w ) C γ 1 + q 1 ( w ) u 2 ( · , w ) v 2 ( · , w ) C γ 2 Γ ( α 1 ) × 0 t ( t s ) α 1 1 d s T α 1 Γ ( 1 + α 1 ) ( p 1 ( w ) u 1 ( · , w ) v 1 ( · , w ) C γ 1 + q 1 ( w ) u 2 ( · , w ) v 2 ( · , w ) C γ 2 ) .
Then,
( N 1 ( u 1 , v 1 ) ) ( · , w ) ( N 1 ( u 2 , v 2 ) ) ( · , w ) C γ 1 T α 1 Γ ( 1 + α 1 ) ( p 1 ( w ) u 1 ( · , w ) v 1 ( · , w ) C γ 1 + q 1 ( w ) u 2 ( · , w ) v 2 ( · , w ) C γ 2 ) .
Furthermore, for any w Ω and each ( u 1 , v 1 ) , ( u 2 , v 2 ) C , and t I , we get:
( N 2 ( u 1 , v 1 ) ) ( · , w ) ( N 2 ( u 2 , v 2 ) ) ( · , w ) C γ 2 T α 2 Γ ( 1 + α 2 ) ( p 2 ( w ) u 1 ( · , w ) v 1 ( · , w ) C γ 1 + q 2 ( w ) u 2 ( · , w ) v 2 ( · , w ) C γ 2 ) .
Thus,
d ( ( N ( u 1 , v 1 ) ) ( · , w ) , ( N ( u 2 , v 2 ) ) ( · , w ) ) M ( w ) d ( ( u 1 ( · , w ) , v 1 ( · , w ) ) , ( u 2 ( · , w ) , v 2 ( · , w ) ) ) ,
where:
d ( ( u 1 ( · , w ) , v 1 ( · , w ) ) , ( u 2 ( · , w ) , v 2 ( · , w ) ) ) = u 1 ( · , w ) v 1 ( · , w ) C γ 1 u 2 ( · , w ) v 2 ( · , w ) C γ 2 .
Since for every w Ω , the matrix M ( w ) converges to zero, then Theorem 1 implies that the operator N has a unique fixed point, which is a random solution of system (1) and (2). □
Now, we prove an existence result for the coupled system (1) and (2) by using the random nonlinear alternative of the Leray–Schauder type in generalized Banach space.
Theorem 4.
Assume that the hypotheses ( H 1 ) and ( H 3 ) hold. Then, the coupled system (1) and (2) has at least one random solution.
Proof. 
We show that the operator N : C × Ω C defined in (8) satisfies all conditions of Theorem 2. The proof will be given in four steps.
Step 1. N ( · , · , w ) is continuous.
Let ( u n , v n ) n be a sequence such that ( u n , v n ) ( u , v ) C as n . For any w Ω and each t I , we have:
t 1 γ 1 ( N 1 ( u n , v n ) ) ( t , w ) t 1 γ 1 ( N 1 ( u , v ) ) ( t , w ) t 1 γ 1 Γ ( α 1 ) 0 t ( t s ) α 1 1 f 1 ( s , u n ( s , w ) , v n ( s , w ) , w ) f 1 ( s , u ( s , w ) , v ( s , w ) , w ) d s 1 Γ ( α 1 ) 0 t ( t s ) α 1 1 s 1 γ 1 f 1 ( s , u n ( s , w ) , v n ( s , w ) , w ) f 1 ( s , u ( s , w ) , v ( s , w ) , w ) d s T α 1 Γ ( 1 + α 1 ) f 1 ( · , u n ( · , w ) , v n ( · , w ) , w ) f 1 ( · , u ( · , w ) , v ( · , w ) , w ) C γ 1 .
Since f 1 is Carathéodory, we have:
( N 1 ( u n , v n ) ) ( · , w ) ( N 1 ( u , v ) ) ( · , w ) C γ 1 0 a s n .
On the other hand, for any w Ω and each t I , we obtain:
t 1 γ 2 ( N 2 ( u n , v n ) ) ( t , w ) t 1 γ 2 ( N 2 ( u , v ) ) ( t , w ) T α 2 Γ ( 1 + α 2 ) f 2 ( · , u n ( · , w ) , v n ( · , w ) , w ) f 2 ( · , u ( · , w ) , v ( · , w ) , w ) C γ 2 .
Furthermore, from the fact that f 2 is Carathéodory, we get:
( N 2 ( u n , v n ) ) ( · , w ) ( N 2 ( u , v ) ) ( · , w ) C γ 2 0 a s n .
Hence, N ( · , · , w ) is continuous.
Step 2. N ( · , · , w ) maps bounded sets into bounded sets in C .
Let R > 0 , and set:
B R : = { ( μ , ν ) C : μ C γ 1 R , ν C γ 2 R } .
For any w Ω and each ( u , v ) B R and t I , we have:
t 1 γ 1 ( N 1 ( u , v ) ) ( t , w ) ϕ 1 ( w ) Γ ( γ 1 ) + t 1 γ 1 Γ ( α 1 ) 0 t ( t s ) α 1 1 f 1 ( s , u ( s , w ) , v ( s , w ) , w ) d s ϕ 1 ( w ) Γ ( γ 1 ) + 1 Γ ( α 1 ) 0 t ( t s ) α 1 1 s 1 γ 1 ( a 1 ( w ) u ( s , w + b 1 ( w ) v ( s , w ) d s ϕ 1 ( w ) Γ ( γ 1 ) + R Γ ( α 1 ) 0 t ( t s ) α 1 1 s 1 γ 1 ( a 1 ( w ) + b 1 ( w ) ) d s ϕ 1 ( w ) Γ ( γ 1 ) + ( a 1 ( w ) + b 1 ( w ) T α 1 Γ ( 1 + α 1 ) : = 1 .
Thus,
( N 1 ( u , v ) ) ( · , w ) C γ 1 1 .
Furthermore, for any w Ω and each ( u , v ) B R and t I , we get:
( N 2 ( u , v ) ) ( · , w ) C γ 2 ϕ 2 ( w ) Γ ( γ 2 ) + ( a 2 ( w ) + b 2 ( w ) T α 2 Γ ( 1 + α ) : = 2 .
Hence,
( N ( u , v ) ) ( · , w ) C ( 1 , 2 ) : = .
Step 3. N ( · , · , w ) maps bounded sets into equicontinuous sets in C .
Let B R be the ball defined in Step 2. For each t 1 , t 2 I with t 1 t 2 and any ( u , v ) B R and w Ω , we have:
t 1 1 γ 1 ( N 1 ( u , v ) ) ( t 1 , w ) t 2 1 γ 1 ( N 1 ( u , v ) ) ( t 2 , w ) t 2 1 γ 1 Γ ( α 1 ) t 1 t 2 ( t 2 s ) α 1 1 f 1 ( s , u ( s , w ) , v ( s , w ) , w ) d s T α 1 Γ ( 1 + α 1 ) ( t 2 t 1 ) α 1 ( a 1 ( w ) u ( · , w ) C γ 1 + b 1 ( w ) v ( · , w ) C γ 2 ) R T α 1 ( a 1 ( w ) + b 1 ( w ) ) Γ ( 1 + α 1 ) ( t 2 t 1 ) α 1 0 a s t 1 t 2 .
Furthermore, we get:
t 1 1 γ 2 ( N 2 ( u , v ) ) ( t 1 , w ) t 2 1 γ 2 ( N 2 ( u , v ) ) ( t 2 , w ) R T α 1 2 ( a 2 ( w ) + b 2 ( w ) ) Γ ( 1 + α 2 ) ( t 2 t 1 ) α 2 0 a s t 1 t 2 .
As a consequence of Steps 1–3, with the Arzela–Ascoli theorem, we conclude that N ( · , · , w ) maps B R into a precompact set in C .
Step 4. The set E ( w ) consisting of ( u ( · , w ) , v ( · , w ) ) C such that ( u ( · , w ) , v ( · , w ) ) = λ ( w ) ( N ( ( u , v ) ) ( · , w ) for some measurable function λ : Ω ( 0 , 1 ) is bounded in C .
Let ( u ( · , w ) , v ( · , w ) ) C such that ( u ( · , w ) , v ( · , w ) ) = λ ( w ) ( N ( ( u , v ) ) ( · , w ) . Then, u ( · , w ) = λ ( w ) ( N 1 ( ( u , v ) ) ( · , w ) and v ( · , w ) = λ ( w ) ( N 2 ( ( u , v ) ) ( · , w ) . Thus, for any w Ω and each t I , we have:
t 1 γ 1 u ( t , w ) ϕ 1 ( w ) Γ ( γ 1 ) + t 1 γ 1 Γ ( α 1 ) 0 t ( t s ) α 1 1 f 1 ( s , u ( s , w ) , v ( s , w ) , w ) d s ϕ 1 ( w ) Γ ( γ 1 ) + 1 Γ ( α 1 ) 0 t ( t s ) α 1 1 s 1 γ 1 ( a 1 ( w ) u ( s , w + b 1 ( w ) v ( s , w ) d s .
Furthermore, we get:
t 1 γ 2 v ( t , w ) ϕ 2 ( w ) Γ ( γ 2 ) + 1 Γ ( α 2 ) 0 t ( t s ) α 2 1 s 1 γ 2 ( a 2 ( w ) u ( s , w + b 2 ( w ) v ( s , w ) d s .
Hence, we obtain:
t 1 γ 1 u ( t , w ) + t 1 γ 2 v ( t , w ) a + b c 0 t ( t s ) α 1 ( s 1 γ 1 u ( s , w + s 1 γ 2 v ( s , w ) d s ,
where:
a : = ϕ 1 ( w ) Γ ( γ 1 ) + ϕ 2 ( w ) Γ ( γ 2 ) , b : = 1 Γ ( α 1 ) + 1 Γ ( α 2 ) ,
c : = max { a 1 ( w ) + a 2 ( w ) , b 1 ( w ) + b 2 ( w ) } , α : = max { α 1 , α 2 } .
Lemma 1 implies that there exists ρ : = ρ ( α ) > 0 such that:
t 1 γ 1 u ( t , w ) + t 1 γ 2 v ( t , w ) a + a b c ρ 0 t ( t s ) α 1 d s a + a b c ρ T α α = L .
This gives:
u ( · , w ) C γ 1 + v ( · , w ) C γ 2 L .
Hence:
( u ( · , w ) , v ( · , w ) ) C L .
This shows that the set E ( w ) is bounded. As a consequence of Steps 1–4 together with Theorem 2, we can conclude that N has at least one fixed point in B R , which is a solution for the system (1) and (2). □

4. Coupled Hilfer–Hadamard Fractional Differential Systems

Now, we are concerned with the coupled system (3) and (4). Set C : = C ( [ 1 , T ] ) , and denote the weighted space of continuous functions defined by:
C γ , ln ( [ 1 , T ] ) = { w ( t ) : ( ln t ) 1 γ w ( t ) C } ,
with the norm:
w C γ , ln : = sup t [ 1 , T ] | ( ln t ) 1 r w ( t ) | .
Furthermore, by C γ 1 , γ 2 , ln ( [ 1 , T ] ) : = C γ 1 , ln ( [ 1 , T ] ) × C γ 2 , ln ( [ 1 , T ] ) , we denote the product weighted space with the norm:
( u , v ) C γ 1 , γ 2 , ln ( [ 1 , T ] ) = u C γ 1 , ln + v C γ 2 , ln .
Let us recall some definitions and properties of Hadamard fractional integration and differentiation. We refer to [7] for a more detailed analysis.
Definition 12.
[7] (Hadamard fractional integral) The Hadamard fractional integral of order q > 0 for a function g L 1 ( [ 1 , T ] ) is defined as:
( H I 1 q g ) ( x ) = 1 Γ ( q ) 1 x ln x s q 1 g ( s ) s d s ,
provided the integral exists.
Example 2.
Let 0 < q < 1 . Let g ( x ) = ln x , x [ 0 , e ] . Then:
( H I 1 q g ) ( x ) = 1 Γ ( 2 + q ) ( ln x ) 1 + q ; f o r a . e . x [ 0 , e ] .
Set:
δ = x d d x , q > 0 , n = [ q ] + 1 ,
and:
A C δ n : = { u : [ 1 , T ] E : δ n 1 [ u ( x ) ] A C ( I ) } .
Definition 13.
[7] The Hadamard fractional derivative of order q > 0 applied to the function w A C δ n is defined as:
( H D 1 q w ) ( x ) = δ n ( H I 1 n q w ) ( x ) .
In particular, if q ( 0 , 1 ] , then:
( H D 1 q w ) ( x ) = δ ( H I 1 1 q w ) ( x ) .
Example 3.
Let 0 < q < 1 . Let w ( x ) = ln x , x [ 0 , e ] . Then:
( H D 1 q w ) ( x ) = 1 Γ ( 2 q ) ( ln x ) 1 q , f o r a . e . x [ 0 , e ] .
It has been proven (see, e.g., Kilbas [30], Theorem 4.8) that in the space L 1 ( I ) , the Hadamard fractional derivative is the left-inverse operator to the Hadamard fractional integral, i.e.:
( H D 1 q ) ( H I 1 q w ) ( x ) = w ( x ) .
From [7], we have:
( H I 1 q ) ( H D 1 q w ) ( x ) = w ( x ) ( H I 1 1 q w ) ( 1 ) Γ ( q ) ( ln x ) q 1 .
The Caputo–Hadamard fractional derivative is defined in the following way:
Definition 14.
The Caputo–Hadamard fractional derivative of order q > 0 applied to the function w A C δ n is defined as:
( H c D 1 q w ) ( x ) = ( H I 1 n q δ n w ) ( x ) .
In particular, if q ( 0 , 1 ] , then:
( H c D 1 q w ) ( x ) = ( H I 1 1 q δ w ) ( x ) .
Definition 15.
Let α ( 0 , 1 ) , β [ 0 , 1 ] , γ = α + β α β , w L 1 ( I ) , and H I 1 ( 1 α ) ( 1 β ) w A C ( I ) . The Hilfer–Hadamard fractional derivative of order α and type β applied to the function w is defined as:
( H D 1 α , β w ) ( t ) = H I 1 β ( 1 α ) ( H D 1 γ w ) ( t ) = H I 1 β ( 1 α ) δ ( H I 1 1 γ w ) ( t ) ; f o r a . e . t [ 1 , T ] .
This new fractional derivative (9) may be viewed as interpolating the Hadamard fractional derivative and the Caputo–Hadamard fractional derivative. Indeed, for β = 0 , this derivative reduces to the Hadamard fractional derivative, and when β = 1 , we recover the Caputo–Hadamard fractional derivative.
H D 1 α , 0 = H D 1 α , a n d H D 1 α , 1 = H c D 1 α .
From [31], we conclude the following lemma.
Lemma 2.
Let g : [ 1 , T ] × E E be such that g ( · , u ( · ) ) C γ , ln ( [ 1 , T ] ) for any u C γ , ln ( [ 1 , T ] ) . Then, Problem (3) is equivalent to the following Volterra integral equation:
u ( t ) = ϕ 0 Γ ( γ ) ( ln t ) γ 1 + ( H I 1 α g ( · , u ( · ) ) ) ( t ) .
Definition 16.
By a random solution of the coupled system (3) and (4), we mean a coupled measurable function ( u , v ) C γ 1 , ln × C γ 2 , ln that satisfies the conditions (4) and Equation (3) on [ 1 , T ] .
Now, we give (without proof) similar existence and uniqueness results for the system (3) and (4). Let us introduce the following hypotheses:
( H 1 )
The functions g i ; i = 1 , 2 are Carathéodory.
( H 2 )
There exist measurable functions p i , q i : Ω ( 0 , ) ; i = 1 , 2 such that:
g i ( t , u 1 , v 1 ) g i ( t , u 2 , v 2 ) p i ( w ) u 1 u 2 + q i ( w ) v 1 v 2 ;
f o r a . e . t [ 1 , T ] , a n d e a c h u i , v i R m , i = 1 , 2 .
( H 3 )
There exist measurable functions a i , b i : Ω ( 0 , ) ; i = 1 , 2 such that:
g i ( t , u , v ) a i ( w ) u + b i ( w ) v ; f o r a . e . t [ 1 , T ] , a n d e a c h u , v R m .
Theorem 5.
Assume that the hypotheses ( H 1 ) and ( H 2 ) hold. If for every w Ω , the matrix:
( ln T ) α 1 Γ ( 1 + α 1 ) p 1 ( w ) ( ln T ) α 1 Γ ( 1 + α 1 ) q 1 ( w ) ( ln T ) α 2 Γ ( 1 + α 2 ) p 2 ( w ) ( ln T ) α 2 Γ ( 1 + α 2 ) q 2 ( w )
converges to zero, then the coupled system (3) and (4) has a unique random solution.
Theorem 6.
Assume that the hypotheses ( H 1 ) and ( H 3 ) hold. Then, the coupled system (3) and (4) has at at least a random solution.

5. An Example

We equip the space R * : = ( , 0 ) with the usual σ -algebra consisting of Lebesgue measurable subsets of R * . Consider the following random coupled Hilfer fractional differential system:
( D 0 1 2 , 1 2 u ) ( t , w ) = f ( t , u ( t , w ) , v ( t , w ) , w ) ; ( D 0 1 2 , 1 2 v ) ( t ) = g ( t , u ( t , w ) , v ( t , w ) , w ) ; ( I 0 1 4 u ) ( 0 , w ) = cos w , ( I 0 1 4 v n ) ( 0 , w ) = sin w , ; w R * , t [ 0 , 1 ] ,
where:
f ( t , u , v , w ) = t 1 4 w 2 ( u ( t ) + v ( t ) ) sin t 64 ( 1 + w 2 + t ) ( 1 + | u | + | v | ) ; t [ 0 , 1 ] ,
g ( t , u , v ) = w 2 ( u ( t ) + v ( t ) ) cos t 64 ( 1 + | u | + | v | ) ; w R * , t [ 0 , 1 ] .
Set α i = β i = 1 2 ; i = 1 , 2 , then γ i = 3 4 ; i = 1 , 2 . The hypothesis ( H 2 ) is satisfied with:
p 1 ( w ) = p 2 ( w ) = q 1 ( w ) = q 2 ( w ) = w 2 64 ( 1 + w 2 ) .
Furthermore, if for every w R * , the matrix:
w 2 64 ( 1 + w 2 ) Γ ( 1 2 ) 1 1 1 1
converges to zero, hence, Theorem 3 implies that the system (10) has a unique random solution defined on [ 0 , 1 ] .

Author Contributions

All the authors contributed in obtaining the results and writing the paper. All authors have read and approved the final manuscript.

Funding

This research was funded by National Natural Science Foundation of China (11671339).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  2. Tarasov, V.E. Fractional Dynamics: Application of Fractional Calculus to Dynamics of Particles, Fields and Media; Springer: Heidelberg, Germany; Higher Education Press: Beijing, China, 2010. [Google Scholar]
  3. Abbas, S.; Benchohra, M.; Graef, J.R.; Henderson, J. Implicit Fractional Differential and Integral Equations: Existence and Stability; De Gruyter: Berlin, Germany, 2018. [Google Scholar]
  4. Abbas, S.; Benchohra, M.; N’Guérékata, G.M. Topics in Fractional Differential Equations; Springer: New York, NY, USA, 2012. [Google Scholar]
  5. Abbas, S.; Benchohra, M.; N’Guérékata, G.M. Advanced Fractional Differential and Integral Equations; Nova Science Publishers: New York, NY, USA, 2015. [Google Scholar]
  6. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives. Theory and Applications; Engl. Trans. from the Russian; Gordon and Breach: Amsterdam, The Netherlands, 1987. [Google Scholar]
  7. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  8. Abbas, S.; Benchohra, M.; Lazreg, J.E.; Zhou, Y. A Survey on Hadamard and Hilfer fractional differential equations: Analysis and Stability. Chaos Solitons Fractals 2017, 102, 47–71. [Google Scholar] [CrossRef]
  9. Benchohra, M.; Henderson, J.; Ntouyas, S.K.; Ouahab, A. Existence results for functional differential equations of fractional order. J. Math. Anal. Appl. 2008, 338, 1340–1350. [Google Scholar] [CrossRef]
  10. Zhou, Y. Attractivity for fractional differential equations in Banach space. Appl. Math. Lett. 2018, 75, 1–6. [Google Scholar] [CrossRef]
  11. Zhou, Y. Attractivity for fractional evolution equations with almost sectorial operators. Fract. Calc. Appl. Anal. 2018, 21, 786–800. [Google Scholar] [CrossRef]
  12. Zhou, Y.; Shangerganesh, L.; Manimaran, J.; Debbouche, A. A class of time-fractional reaction-diffusion equation with nonlocal boundary condition. Math. Meth. Appl. Sci. 2018, 41, 2987–2999. [Google Scholar] [CrossRef]
  13. Zhou, Y.; Peng, L.; Huang, Y.Q. Duhamel’s formula for time-fractional Schrödinger equations. Math. Meth. Appl. Sci. 2018, 41, 8345–8349. [Google Scholar] [CrossRef]
  14. Zhou, Y.; Peng, L.; Huang, Y.Q. Existence and Hölder continuity of solutions for time-fractional Navier-Stokes equations. Math. Meth. Appl. Sci. 2018, 41, 7830–7838. [Google Scholar] [CrossRef]
  15. Furati, K.M.; Kassim, M.D. Non-existence of global solutions for a differential equation involving Hilfer fractional derivative. Electron. J. Differ. Equ. 2013, 2013, 1–10. [Google Scholar]
  16. Furati, K.M.; Kassim, M.D.; Tatar, N.E. Existence and uniqueness for a problem involving Hilfer fractional derivative. Comput. Math. Appl. 2012, 64, 1616–1626. [Google Scholar] [CrossRef] [Green Version]
  17. Kamocki, R.; Obczńnski, C. On fractional Cauchy-type problems containing Hilfer’s derivative. Electron. J. Qual. Theory Differ. Equ. 2016, 2016, 1–12. [Google Scholar] [CrossRef]
  18. Qassim, M.D.; Furati, K.M.; Tatar, N.-E. On a differential equation involving Hilfer-Hadamard fractional derivative. Abstr. Appl. Anal. 2012, 2012, 391062. [Google Scholar] [CrossRef]
  19. Tomovski, Ž.; Hilfer, R.; Srivastava, H.M. Fractional and operational calculus with generalized fractional derivative operators and Mittag-Leffler type functions. Integral Transform. Spec. Funct. 2010, 21, 797–814. [Google Scholar] [CrossRef]
  20. Wang, J.-R.; Zhang, Y. Nonlocal initial value problems for differential equations with Hilfer fractional derivative. Appl. Math. Comput. 2015, 266, 850–859. [Google Scholar] [CrossRef]
  21. Abbas, S.; Albarakati, W.; Benchohra, M.; Zhou, Y. Weak solutions for partial pettis Hadamard fractional integral equations with random effects. J. Integral Equ. Appl. 2017, 29, 473–491. [Google Scholar] [CrossRef]
  22. Abbas, S.; Benchohra, M.; Zhou, Y. Coupled Hilfer fractional differential systems with random effects. Adv. Differ. Equ. 2018, 2018, 369. [Google Scholar] [CrossRef]
  23. Graef, J.R.; Henderson, J.; Ouahab, A. Some Krasnosel’skii type random fixed point theorems. J. Nonlinear Funct. Anal. 2017, 2017, 46. [Google Scholar]
  24. Petre, I.R.; Petrusel, A. Krasnoselskii’s theorem in generalized Banach spaces and applications. Electron. J. Qual. Theory Differ. Equ. 2012, 85, 1–20. [Google Scholar] [CrossRef]
  25. Sinacer, M.L.; Nieto, J.J.; Ouahab, A. Random fixed point theorems in generalized Banach spaces and applications. Random Oper. Stoch. Equ. 2016, 24, 93–112. [Google Scholar] [CrossRef]
  26. Itoh, S. Random fixed point theorems with applications to random differential equations in Banach spaces. J. Math. Anal. Appl. 1979, 67, 261–273. [Google Scholar] [CrossRef]
  27. Engl, H.W. A general stochastic fixed-point theorem for continuous random operators on stochastic domains. J. Math. Anal. Appl. 1978, 66, 220–231. [Google Scholar] [CrossRef]
  28. Allaire, G.; Kaber, S.M. Numerical Linear Algebra; Ser. Texts in Applied Mathematics; Springer: New York, NY, USA, 2008. [Google Scholar]
  29. Varga, R.S. Matrix Iterative Analysis, Second Revised and Expanded Edition; Springer Series in Computational Mathematics 27; Springer: Berlin, Germany, 2000. [Google Scholar]
  30. Kilbas, A.A. Hadamard-type fractional calculus. J. Korean Math. Soc. 2001, 38, 1191–1204. [Google Scholar]
  31. Qassim, M.D.; Tatar, N.-E. Well-posedness and stability for a differential problem with Hilfer-Hadamard fractional derivative. Abstr. Appl. Anal. 2013, 2013, 605029. [Google Scholar]

Share and Cite

MDPI and ACS Style

Abbas, S.; Arifi, N.A.; Benchohra, M.; Zhou, Y. Random Coupled Hilfer and Hadamard Fractional Differential Systems in Generalized Banach Spaces. Mathematics 2019, 7, 285. https://0-doi-org.brum.beds.ac.uk/10.3390/math7030285

AMA Style

Abbas S, Arifi NA, Benchohra M, Zhou Y. Random Coupled Hilfer and Hadamard Fractional Differential Systems in Generalized Banach Spaces. Mathematics. 2019; 7(3):285. https://0-doi-org.brum.beds.ac.uk/10.3390/math7030285

Chicago/Turabian Style

Abbas, Saïd, Nassir Al Arifi, Mouffak Benchohra, and Yong Zhou. 2019. "Random Coupled Hilfer and Hadamard Fractional Differential Systems in Generalized Banach Spaces" Mathematics 7, no. 3: 285. https://0-doi-org.brum.beds.ac.uk/10.3390/math7030285

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