Next Article in Journal
The Darboux Transformation and N-Soliton Solutions of Coupled Cubic-Quintic Nonlinear Schrödinger Equation on a Time-Space Scale
Previous Article in Journal
Numerical Simulation of Fractional Delay Differential Equations Using the Operational Matrix of Fractional Integration for Fractional-Order Taylor Basis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Existence Results for Hilfer Fractional Differential Equations with Variable Coefficient

School of Mathematics, Yunnan Normal University, Kunming 650092, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(1), 11; https://doi.org/10.3390/fractalfract6010011
Submission received: 16 November 2021 / Revised: 11 December 2021 / Accepted: 22 December 2021 / Published: 28 December 2021

Abstract

:
The aim of this paper is to establish the existence and uniqueness results for differential equations of Hilfer-type fractional order with variable coefficient. Firstly, we establish the equivalent Volterra integral equation to an initial value problem for a class of nonlinear fractional differential equations involving Hilfer fractional derivative. Secondly, we obtain the existence and uniqueness results for a class of Hilfer fractional differential equations with variable coefficient. We verify our results by providing two examples.

1. Introduction

Fractional differential equation is an interesting research field. The reason is that it can be used to solve practical problems from the fields of science and engineering, such as physics, chemistry, electrodynamics of complex media, polymer rheology and so on ([1,2,3,4,5,6]). Recently the research of fractional differential equations has made great progress. In the literature, there are several definitions of fractional integrals and derivatives, the most popular definitions are in the sense of the Riemann–Liouville and Caputo derivatives. In 2000, Hilfer introduced a generalized Riemann–Liouville fractional derivative, the socalled Hilfer fractional derivative. Many authors studied the existence of solutions for fractional differential equations involving Hilfer fractional derivative (see [7,8,9,10,11]).
However, there are few studies on fractional differential equations involving Hilfer fractional derivative with variable coefficient. Due to the variable coefficient function, it is difficult to obtain the analytical solution of such equations ([12,13,14,15,16,17,18]). The representation of explicit solutions has become a problem, even for the case involving Riemann–Liouville or Caputo derivative, the problem is not completely solved. In [15], the authors considered the following linear Caputo fractional differential equation with constant coefficients and obtained solutions by Adomian decomposition method:
i = 0 n α i L D β i y i ( t ) = f ( t ) , t [ a , b ] .
In [16], using Neumann series for the corresponding Volterra integral equations and the generalized Mittag–Leffler functions, the authors obtained solutions of initial value problems of fractional differential equations with constant coefficients:
L D β n y n ( t ) + i = 0 n 1 α i L D β i y i ( t ) = f ( t ) , t [ a , b ] .
In [17], the authors investigated solutions around an ordinary point for linear homogeneous sequential Caputo fractional differential equations with variable coefficients:
L D β n y n ( t ) + i = 0 n 1 α i ( t ) L D β i y i ( t ) = f ( t ) , t [ a , b ] .
In [18], the authors study the existence, uniqueness and stability of solutions of the implicit fractional differential equations as follows:
H D 0 + α , β x ( t ) = f ( t , H D 0 + α , β x ( t ) ) , t J : = ( 0 , T ] , ( I 0 + n γ x ) ( n j ) ( 0 + ) = i = 0 m c i x ( τ i ) , α < γ = α + β α β < 1 , j = 1 , 2 , , n ,
where H D 0 + α , β is the Hilfer fractional derivative, α , β ( 0 , 1 ) , f : J × R × R R is a given continuous function.
In [19], the authors study the following fractional differential equation with continuous variable coefficients and Hilfer fractional derivatives with respect to another function:
H D 0 + β 0 , μ 0 , φ y ( t ) + i = 1 m σ i ( t ) H D 0 + β i , μ i , φ y ( t ) = g ( t ) , t [ 0 , t 0 ] , m N , 1 φ ( t ) d d t k ( I 0 + ( 1 μ 0 ) ( n β 0 ) ; φ y ) ( 0 + ) = a k R , k = 0 , 1 , , n 0 1 ,
where β 0 > β 1 > > β m 0 , 0 μ i 1 , σ i , g C [ 0 , t 0 ] and n i are non-negative integers satisfying n i 1 < β i < n i , i = 0 , 1 , , m .
In this paper, we consider the following initial value problem for Hilfer fractional differential equations with variable coefficient:
H D 0 + α , β x ( t ) λ ( t ) x ( t ) = f ( t , x ( t ) ) , t J : = ( 0 , T ] , ( I 0 + n γ x ) ( n j ) ( 0 + ) = c j , α < γ = α + n β α β < n , j = 1 , 2 , , n ,
where α ( n 1 , n ) , β ( 0 , 1 ) , H D 0 + α , β is the Hilfer fractional derivative, f : J × R R is a function to be specified later. I 0 + n γ is the left-sided Riemann–Liouville fractional integral of order n γ .
Using new techniques, we prove the existence result under the weak assumptions for the variable coefficient λ L 1 p and nonlinearity f ( t , x ( t ) ) (Theorem 6), the main advantage of our techniques is that we are able to consider a wide range of function spaces, and in which the Hilfer fractional derivative of f exists.
In many cases, a representation of the solution of variable coefficient fractional differential equations is still an open question. In this paper, we give the structure of solutions for variable coefficient fractional differential equations with the Hilfer-type fractional derivative and provide a better understanding of the structure of solutions of Hilfer-type fractional differential equations with variable coefficients. The results are new even for the special case: β = 0 or β = 1 . Moreover, we find the explicit solutions for the linear Hilfer fractional differential equations with variable coefficient as follows:
H D 0 + α , β x ( t ) λ ( t ) x ( t ) = f ( t ) , t J : = ( 0 , T ] , ( I 0 + n γ x ) ( n j ) ( 0 + ) = c j , α < γ = α + n β α β < n , j = 1 , 2 , , n .
Under appropriate assumptions, we can prove the solution x is given by
x ( t ) = j = 1 n k = 0 c j ( I 0 + α λ ( · ) ) k t γ j Γ ( γ j + 1 ) + k = 1 ( I 0 + α λ ( · ) ) k 1 ( I 0 + α f ) ( t ) .
Furthermore, we obtain the existence of solutions for the problem (1) under the new assumptions for f, and we can find that the linear nonhomogeneous Hilfer fractional differential equations and the Hilfer fractional differential equations with constant coefficient are the special cases of our conclusion.
The article is partitioned as follows: In Section 2, the basic definitions and conclusions are presented. In Section 3, we present the equivalent Volterra integral equation. The main results are obtained in Section 4. The applications are shown in Section 5.

2. Preliminaries

In this paper, let Ω = [ a , b ] be a finite interval on R . We denote by L p ( Ω , R ) the Banach space of all Lebesgue measurable functions l : Ω R with the norm l L p = Ω | l ( t ) | p d t 1 p < and by A C ( [ a , b ] , R ) the space of all the absolutely continuous functions defined on [ a , b ] . Moreover, we use the following notation:
A C n ( [ a , b ] , R ) = { f ; f C n 1 ( [ a , b ] , R ) and f ( n 1 ) A C ( [ a , b ] , R ) } .
In particular, A C 1 ( [ a , b ] , R ) = A C ( [ a , b ] , R ) .
For 0 ν < 1 , we denote the weighted spaces
C ν [ a , b ] : = { x : ( a , b ] R ; ( t a ) ν x ( t ) C ( [ a , b ] , R ) } with the norm x C ν = max t [ a , b ] ( t a ) ν | x ( t ) | .
First, we recall some basic concepts and results which will be used in the sequel.
Definition 1
([3,4]). The left-sided fractional integral of order q for a function x ( t ) L 1 is defined by
( I a + q x ) ( t ) = 1 Γ ( q ) a t ( t s ) q 1 x ( s ) d s , t > a , q > 0 ,
where Γ ( · ) is the Gamma function.
Definition 2
([3,4]). If k ( t ) A C n ( [ a , b ] , R ) , then the left-sided Riemann–Liouville fractional derivative ( L D a + q k ) ( t ) of order q exists almost everywhere on [ a , b ] and can be written as
( L D a + q k ) ( t ) = 1 Γ ( n q ) d n d t n a t ( t s ) n q 1 k ( s ) d s , t > a , n 1 < q < n .
Lemma 1
([4]). If μ > 0 and θ > 0 , then
[ I a + μ ( s a ) θ 1 ] ( t ) = Γ ( θ ) Γ ( θ + μ ) ( t a ) θ + μ 1 , L D a + μ ( s a ) θ 1 ( t ) = Γ ( θ ) Γ ( θ μ ) ( t a ) θ μ 1 , L D a + μ ( s a ) μ j ( t ) = 0 , j = 1 , 2 , , [ μ ] + 1 .
Lemma 2
([4]). If k ( t ) L p ( a , b ) ( 1 p ) and θ 1 , θ 2 > 0 , then the following relations hold:
(i) 
( L D a + θ 1 I a + θ 1 k ) ( t ) = k ( t ) a.e. t [ a , b ] ;
(ii) 
For θ 1 > n , ( d n d t n I a + θ 1 ) k ( t ) = I a + θ 1 n k ( t ) .
(iii) 
For almost every point t [ a , b ] , ( I a + θ 1 I a + θ 2 k ) ( t ) = ( I a + θ 1 + θ 2 k ) ( t ) .
Lemma 3
([4]). Let ν ( n 1 , n ) ( n N ) . If y L 1 ( a , b ) and I a + n ν y A C n [ a , b ] , then
I a + ν L D a + ν y ( t ) = y ( t ) j = 1 n ( I a + n ν y ) ( n j ) ( a + ) Γ ( ν j + 1 ) ( t a ) ν j
holds almost everywhere on [ a , b ] .
In [7], R. Hilfer studied a generalized fractional operator having the Riemann–Liouville and Caputo derivatives as specific cases.
Definition 3
([4]). The left-sided Hilfer fractional derivative of order α ( n 1 , n ) , β [ 0 , 1 ] of y ( t ) is defined by:
H D a + α , β y ( t ) = ( I a + β ( n α ) D n I a + ( 1 β ) ( n α ) y ) ( t ) = ( I a + γ α D n I a + n γ y ) ( t ) , t [ a , b ] ,
where D : = d d t , γ = α + n β α β .
Remark 1. 
(i) 
The operator H D a + α , β can be written as:
H D a + α , β = I a + β ( n α ) D n ( I a + ( 1 β ) ( n α ) ) = I a + β ( n α ) L D a + γ = I a + γ α L D a + γ .
(ii) 
when β = 0 , the left-sided Riemann–Liouville fractional derivative can be presented as L D a + α : = H D a + α , 0 .
(iii) 
when β = 1 , the left-sided Caputo fractional derivative can be presented as C D a + α : = H D a + α , 1 .
From Lemma 2 (iii), we can obtain the following result.
Lemma 4.
Let y L 1 ( a , b ) and L D a + γ α y L 1 ( a , b ) exists, then
H D a + α , β I a + α y ( t ) = I a + γ α L D a + γ α y ( t ) .
Proof. 
H D a + α , β I a + α y ( t ) = I a + γ α D n I a + n γ I a + α y ( t ) = I a + γ α D n I a + n + α γ y ( t ) = I a + γ α L D a + γ α y ( t ) .
 □
From (2), we can derive the following result.
Lemma 5.
For γ = α + n β α β ( n 1 , n ) , a general solution of the fractional differential equation H D a + α , β x ( t ) = 0 is given by
x ( t ) = i = 1 n c i ( t a ) γ i , t > a ,
where c i R , i = 1 , 2 , , n are arbitrary constants.
Next, we present the following lemmas.
Lemma 6.
If ω > 0 , 1 < τ 0 , then for ψ L 1 p ( 0 < p < 1 ) ,
(i) 
a t ( t s ) ω 1 ψ ( s ) d s 1 p ω p 1 p ( t a ) ω p ψ L 1 p , for ω > p , t > a ,
(ii) 
0 t ( t s ) ω 1 s τ ψ ( s ) d s B ω p 1 p , τ + 1 p 1 p 1 p t ω + τ p ψ L 1 p , for p < min { ω , 1 + τ } , t > 0 ,
where B ( · , · ) is the Beta function.
Proof. 
With the help of Hölder’s inequality, we get
a t ( t s ) ω 1 ψ ( s ) d s a t ( t s ) ω 1 1 p d s 1 p ψ L 1 p = 1 p ω p 1 p ( t a ) ω p ψ L 1 p , for ω > p , t > a . 0 t ( t s ) ω 1 s τ ψ ( s ) d s 0 t ( t s ) ω 1 1 p s τ 1 p d s 1 p ψ L 1 p = B ω p 1 p , 1 + τ p 1 p 1 p t ω + τ p ψ L 1 p , for p < min { ω , 1 + τ } , t > 0 .
Lemma 7.
If ω ( 0 , 1 ) , τ ( 1 , 0 ] and ω + τ > 0 , then for 0 < p < ω + τ 2 and y L 1 p , ( I 0 + ω s τ y ( s ) ) ( t ) A C ( [ 0 , T ] , R ) .
Proof. 
For every finite collection { ( a i , b i ) } 1 i n on J with i = 1 n ( b i a i ) 0 , it follows from Hölder’s inequality that
i = 1 n | ( I 0 + ω s τ y ( s ) ) ( b i ) ( I 0 + ω s τ y ( s ) ) ( a i ) | 1 Γ ( ω ) i = 1 n 0 a i [ ( b i s ) ω 1 ( a i s ) ω 1 ] s τ y ( s ) d s + i = 1 n a i b i ( b i s ) ω 1 s τ | y ( s ) | d s y L 1 p Γ ( ω ) i = 1 n 0 a i | ( b i s ) ω 1 ( a i s ) ω 1 | 1 1 p s τ 1 p d s 1 p + i = 1 n a i b i ( b i s ) ω 1 1 p s τ 1 p d s 1 p y L 1 p Γ ( ω ) { ( 1 ω ) i = 1 n 0 a i a i b i ( ζ s ) ω 2 d ζ 1 1 p s τ 1 p d s 1 p + B ω p 1 p , τ + 1 p 1 p 1 p i = 1 n ( b i a i ) ω + τ p } M ¯ i = 1 n 0 a i ( ( a i s ) θ ( b i s ) θ ) s τ 1 p d s 1 p + y L 1 p Γ ( ω ) B ω p 1 p , τ + 1 p 1 p 1 p i = 1 n ( b i a i ) ω + τ p = M ¯ i = 1 n { ( a i θ + τ + 1 p 1 p b i θ + τ + 1 p 1 p ) B θ + 1 , τ + 1 p 1 p + b i θ + τ + 1 p 1 p B θ + 1 , τ + 1 p 1 p 0 a i b i ( 1 μ ) θ μ τ 1 p d μ } 1 p + y L 1 p Γ ( ω ) B ω p 1 p , τ + 1 p 1 p 1 p i = 1 n ( b i a i ) ω + τ p 0 ,
where M ¯ > 0 is a constant and θ = ω 1 p 1 p ( 1 , 0 ) . Now, ( I 0 + ω s τ y ( s ) ) ( t ) A C ( [ 0 , T ] , R ) . □
The following theorem will be used to prove the main result.
Theorem 1
([20] Schauder Fixed Point Theorem). If U is a nonempty, closed, bounded convex subset of a Banach space X and T : U U is completely continuous, then T has a fixed point in U.

3. Equivalent Volterra Integral Equation

We consider the following linear Hilfer fractional differential equation with initial value conditions
H D 0 + α , β x ( t ) = f ( t ) , t J : = ( 0 , T ] , ( I 0 + n γ x ) ( n j ) ( 0 + ) = c j , α < γ = α + n β α β < n , j = 1 , 2 , , n ,
where f L 1 p ( J ) ( 0 < p < 1 γ + α 2 ) .
Definition 4.
A function x C n γ satisfying (3) is called a solution of (3).
Theorem 2.
A function x C n γ satisfies a.e. (3) if and only if x satisfies a.e. the following integral equation
x ( t ) = j = 1 n c j t γ j Γ ( γ j + 1 ) + ( I 0 + α f ) ( t ) .
Proof. 
(Necessity) For t J , γ α ( 0 , 1 ) , it follows from Lemma 7 ( τ = 0 ) that I 0 + 1 ( γ α ) f A C [ 0 , T ] and L D 0 + γ α f L 1 ( J ) . From Lemmas 3 and 4, we arrive at
( H D 0 + α , β I 0 + α f ) ( t ) = I 0 + γ α L D 0 + γ α f ( t ) = f ( t ) ( I 0 + 1 γ + α f ) ( 0 + ) Γ ( γ α ) t γ α 1 , a . e . t [ 0 , T ] ,
then
H D 0 + α , β x ( t ) ( I 0 + α f ) ( t ) ( I 0 + 1 γ + α f ) ( 0 + ) Γ ( γ α ) t γ α 1 = 0 .
Since ( I 0 + 1 γ + α f ) ( 0 + ) = 0 , then
H D 0 + α , β x ( t ) ( I 0 + α f ) ( t ) = 0 .
Using Lemma 5, we find
x ( t ) = j = 1 n d j t γ j + ( I 0 + α f ) ( t ) .
Noting that the initial value conditions ( I 0 + n γ x ) ( n j ) ( 0 + ) = c j ( j = 1 , 2 , , n ) , we can obtain d j = Γ ( γ j + 1 ) ( I 0 + n γ x ) ( n j ) ( 0 + ) = c j ( j = 1 , 2 , , n ) .
(Sufficiency) Let x ( t ) satisfy (4). It follows from Lemmas 1, 3 and 7 that H D 0 + α , β x ( t ) exists and x ( t ) satisfies (3) by direct computation. This completes the proof. □
Next, we consider the following nonlinear Hilfer fractional differential equation with initial value conditions
H D 0 + α , β x ( t ) = f ( t , x ( t ) ) , t J : = ( 0 , T ] , ( I 0 + n γ x ) ( n j ) ( 0 + ) = c j , α < γ = α + n β α β < n , j = 1 , 2 , , n .
Definition 5.
A function x C n γ satisfying (5) is called a solution of (5).
Theorem 3.
Let f : ( 0 , T ] × R R be a function such that f ( · , x ( · ) ) L 1 p ( J , R ) ( 0 < p < 1 γ + α 2 ) for any x C n γ ( [ 0 , T ] , R ) . If x ( t ) C n γ ( [ 0 , T ] , R ) , then x satisfies a.e. (5) if and only if x satisfies a.e. the following equation
x ( t ) = j = 1 n c j t γ j Γ ( γ j + 1 ) + 0 t ( t s ) α 1 f ( s , x ( s ) ) Γ ( α ) d s .
The proof of Theorem 3 is nearly same as Theorem 2, we omit it here.

4. Existence and Uniqueness Result

Theorem 4.
Let f : ( 0 , T ] × R R be a function such that f ( · , x ( · ) ) L 1 p ( J , R ) ( n γ < p < 1 + α n 2 ) for any x C n γ ( [ 0 , T ] , R ) and there exists a function μ L 1 p ( J , R + ) ( n γ < p < 1 + α n 2 ) such that
| f ( t , u ( t ) ) f ( t , v ( t ) ) | μ ( t ) | u ( t ) v ( t ) |
and M f : = sup t J | f ( t , 0 ) | < . Then the problem (5) has a unique solution x C n γ .
Proof. 
Clearly, if x C n γ , then
| f ( t , x ( t ) ) | | f ( t , x ( t ) ) f ( t , 0 ) | + | f ( t , 0 ) | μ ( t ) | x ( t ) | + M f , t n γ | f ( t , x ( t ) ) | μ ( t ) t γ n x C n γ + t n γ M f .
Applying with Lemma 7 (with ω and τ replaced by 1 γ + α and γ n , respectively), we can verify that ( I 0 + 1 ( γ α ) f ) ( t ) A C [ 0 , T ] and ( L D 0 + γ α f ) ( t ) L 1 ( J ) for any x C n γ .
Define an operator F : C n γ C n γ as
( F x ) ( t ) = j = 1 n c j t γ j Γ ( γ j + 1 ) + ( I 0 + α F x ) ( t ) ,
where F x ( t ) = f ( t , x ( t ) ) . Clearly, F is well defined and the fixed point of F is the solution of the problem (5).
It is easy to check that
t n γ | ( F x ) ( t ) ( F y ) ( t ) | t n γ Γ ( α ) 0 t ( t s ) α 1 | μ ( s ) | | x ( s ) y ( s ) | d s t n γ Γ ( α ) 0 t ( t s ) α 1 1 p s γ n 1 p d s 1 p · μ L 1 p x y C n γ t α p Γ ( α ) B ( α p 1 p , γ n + 1 p 1 p ) 1 p · μ L 1 p x y C n γ
and
t n γ | ( F 2 x ) ( t ) ( F 2 y ) ( t ) | t n γ I 0 + α | μ ( s ) | s γ n s n γ | ( F x ) ( s ) ( F y ) ( s ) | ( t ) t n γ Γ 2 ( α ) 0 t ( t s ) α 1 | μ ( s ) | s γ n s α p d s B ( α p 1 p , γ n + 1 p 1 p ) 1 p · μ L 1 p x y C n γ t 2 ( α p ) Γ 2 ( α ) B ( α p 1 p , γ n + 1 p + α p 1 p ) 1 p B ( α p 1 p , γ n + 1 p 1 p ) 1 p · μ L 1 p 2 x y C n γ = t 2 ( α p ) Γ 2 ( α ) Γ 2 ( α p 1 p ) Γ ( γ n 1 p + 1 ) Γ ( γ n + 2 ( α p ) 1 p + 1 ) 1 p · μ L 1 p 2 x y C n γ .
By induction we deduce that
t n γ | ( F k x ) ( t ) ( F k y ) ( t ) | T k ( α p ) Γ k ( α ) Γ k ( α p 1 p ) Γ ( γ n 1 p + 1 ) Γ ( γ n + k ( α p ) 1 p + 1 ) 1 p · μ L 1 p k x y C n γ ,
for k large enough one has T k ( α p ) Γ k ( α ) Γ k ( α p 1 p ) Γ ( γ n 1 p + 1 ) Γ ( γ n + k ( α p ) 1 p + 1 ) 1 p · μ L 1 p k < 1 and by Banach contraction principle, F has a unique fixed point x C n γ . □
Remark 2.
Obviously, if n γ < p < 1 α + n 2 , then the relation C n γ L 1 p holds, but the converse may not be true, for example, n = 1 , α = 0.9 , β = 0.1 , γ = α + β α β = 0.91 , t 1 10 L 9 ( 0 , T ) but t 1 10 C 0.09 [ 0 . T ] .
Special Case I. If we take f ( t , x ( t ) ) = λ ( t ) x ( t ) + f ( t ) , obviously, | f ( t , x ( t ) ) f ( t , y ( t ) ) | | λ ( t ) | | x ( t ) y ( t ) | , then the following result holds.
Corollary 1.
Let λ ( t ) C ( [ 0 , T ] ) and f ( t ) L 1 p ( J ) ( n γ < p < 1 + α n 2 ) and sup t J | f ( t ) | < , then the linear Hilfer fractional differential equations with variable coefficient
H D 0 + α , β x ( t ) λ ( t ) x ( t ) = f ( t ) , t J : = ( 0 , T ] , ( I 0 + n γ x ) ( n j ) ( 0 + ) = c j , α < γ = α + n β α β < n , j = 1 , 2 , , n
has a unique solution x C n γ and x is given by
x ( t ) = j = 1 n k = 0 c j ( I 0 + α λ ( · ) ) k t γ j Γ ( γ j + 1 ) + k = 1 ( I 0 + α λ ( · ) ) k 1 ( I 0 + α f ) ( t ) .
Proof. 
By Theorem 4, we can see that the solution of the problem (7) is given by
x ( t ) = j = 1 n c j t γ j Γ ( γ j + 1 ) + ( I 0 + α λ ( · ) x ) ( t ) + ( I 0 + α f ) ( t ) .
It follows from Theorem 4 that the operator F : C n γ C n γ defined as
( F x ) ( t ) = j = 1 n c j t γ j Γ ( γ j + 1 ) + ( I 0 + α λ ( · ) x ) ( t ) + ( I 0 + α f ) ( t )
has a unique fixed point which is the solution of the problem (7).
Let
x 0 ( t ) = j = 1 n c j t γ j Γ ( γ j + 1 ) , x m ( t ) = x 0 ( t ) + ( I 0 + α λ ( · ) x m 1 ) ( t ) + ( I 0 + α f ) ( t ) ,
iterating, for m N , we can write
x m ( t ) = k = 0 m ( I 0 + α λ ( · ) ) k x 0 ( t ) + k = 1 m ( I 0 + α λ ( · ) ) k 1 ( I 0 + α f ) ( t ) .
Taking the limit as m , we can get
x ( t ) = j = 1 n k = 0 c j ( I 0 + α λ ( · ) ) k t γ j Γ ( γ j + 1 ) + k = 1 ( I 0 + α λ ( · ) ) k 1 ( I 0 + α f ) ( t ) .
Special Case II. In particular, when λ ( t ) λ R is a constant, we derive the following result.
Corollary 2.
Let f ( t ) L 1 p ( J ) ( n γ < p < 1 + α n 2 ) and sup t J | f ( t ) | < . Then the initial value problem
H D 0 + α , β x ( t ) λ x ( t ) = f ( t ) , t J : = ( 0 , T ] , λ R ( I 0 + n γ x ) ( n j ) ( 0 + ) = c j , α < γ = α + n β α β < n , j = 1 , 2 , , n
has a unique solution x C n γ and x is given by
x ( t ) = j = 1 n c j t γ j E α , γ j + 1 ( λ t α ) + 0 t ( t s ) α 1 E α , α ( λ ( t s ) α ) f ( s ) d s ,
where
E μ , ν ( z ) = k = 0 z k Γ ( μ k + ν ) , μ , ν , z C , R e ( μ ) > 0
is a Mittag–Leffler function ([4]).
Proof. 
It follows from Theorem 4 that the problem (8) has a unique solution x C n γ and the solution x is given by the limit x ( t ) = lim m + x m ( t ) of the sequence
x 0 ( t ) = j = 1 n c j t γ j Γ ( γ j + 1 ) , x m ( t ) = x 0 ( t ) + λ ( I 0 + α x m 1 ) ( t ) + ( I 0 + α f ) ( t ) .
Iterating, for m N , we arrive at
x m ( t ) = k = 0 m λ k I 0 + α k x 0 ( t ) + k = 1 m λ k 1 I 0 + α k f ( t ) = j = 1 n k = 0 m λ k c j Γ ( γ j + 1 ) I 0 + α k t γ j + 0 t k = 1 m λ k 1 ( t s ) α k 1 Γ ( α k ) f ( s ) d s = j = 1 n k = 0 m λ k c j t γ j + α k + 1 Γ ( γ j + α k ) + 0 t k = 1 m λ k 1 ( t s ) α k 1 Γ ( α k ) f ( s ) d s .
Taking the limit as m , we have
x ( t ) = lim m x m ( t ) = j = 1 n k = 0 λ k c j t γ j + α k Γ ( γ j + α k + 1 ) + 0 t k = 1 λ k 1 ( t s ) α k 1 Γ ( α k ) f ( s ) d s = j = 1 n c j t γ j E α , γ j + 1 ( λ t α ) + 0 t ( t s ) α 1 E α , α ( λ ( t s ) α ) f ( s ) d s .
This yields the explicit solution to the problem (8). □
Next, we deal with the existence of solutions to the problem (1). We consider the following assumptions.
(H1)
f : ( 0 , T ] × R R is a function such that f ( · , x ( · ) ) L 1 p ( J , R ) ( n γ < p < 1 + α n 2 ) for any x C n γ ( [ 0 , T ] , R ) and f ( t , · ) : R R is continuous for a.e. t J ;
(H2)
there exists a function μ L 1 p 1 ( J , R + ) ( n γ < p 1 < 1 + α n 2 ) such that
| f ( t , u ) | μ ( t ) | u ( t ) | σ , σ ( 0 , 1 ) .
Theorem 5.
Assume that λ L 1 p 2 ( n γ < p 2 < 1 + α n 2 ) and (H1) and (H2) are satisfied, a function x C n γ is a solution of (1) if and only if x is a solution of the following equation
x ( t ) = j = 1 n c j t γ j Γ ( γ j + 1 ) + ( I 0 + α λ ( · ) x ) ( t ) + ( I 0 + α F x ) ( t ) .
Proof. 
Clearly, if x C n γ , then
t n γ | f ( t , x ( t ) ) | μ ( t ) t σ ( γ n ) x C n γ σ .
According to Lemma 7 (with ω and τ replaced by 1 γ + α and σ ( γ n ) , respectively), we can see that I 0 + 1 ( γ α ) ( λ ( · ) x + F x ) ( t ) A C [ 0 , T ] and L D a + γ α ( λ ( · ) x + F x ) ( t ) L 1 ( J ) for any x C n γ . Then the above result can be proved in a similar way as in the proof of Theorem 2. □
For the sake of convenience, we set two constants
Λ 1 : = B α p 2 1 p 2 , γ n + 1 p 2 1 p 2 1 p 2 λ L 1 p 2 Γ ( α ) , Λ 2 : = B α p 1 1 p 1 , σ ( γ n ) + 1 p 1 1 p 1 1 p 1 μ L 1 p 1 Γ ( α ) .
Theorem 6.
Assume that λ L 1 p 2 ( n γ < p 2 < 1 + α n 2 ) and (H1) and (H2) are satisfied, then the problem (1) has at least a solution x C n γ if Λ 1 T α p 2 < 1 .
Proof. 
We define an operator F : C n γ C n γ as
( F x ) ( t ) = j = 1 n c j t γ j Γ ( γ j + 1 ) + ( I 0 + α λ ( · ) x ) ( t ) + ( I 0 + α F x ) ( t ) .
Clearly, by Lemma 5, F is well defined, and the fixed point of F is the solution of the problem (1).
For x C n γ , from Lemma 6 (ii) and (9), we have
| I 0 + α λ ( s ) x ( s ) ( t ) | 1 Γ ( α ) 0 t ( t s ) α 1 s γ n | λ ( s ) | d s · x C n γ Λ 1 t α n + γ p 2 · x C n γ , | ( I 0 + α F x ) ( t ) | 1 Γ ( α ) 0 t ( t s ) α 1 s σ ( γ n ) μ ( s ) d s · x C n γ σ Λ 2 t α + σ ( γ n ) p 1 · x C n γ σ , t n γ | ( F x ) ( t ) | t n γ j = 1 n c j t γ j Γ ( γ j + 1 ) + | ( I 0 + α λ ( · ) x ) ( t ) | + | ( I 0 + α F x ) ( t ) | j = 1 n c j t n j Γ ( γ j + 1 ) + Λ 1 t α p 2 · x C n γ + Λ 2 t α p 1 + ( 1 σ ) ( n γ ) · x C n γ σ .
Step I. For an r > 0 , we set the ball B r C n γ as B r = { x C n γ ; x C n γ r } . We claim that there exists an r 0 > 0 such that F B r 0 B r 0 . If this is not true, then for each positive number r, there exists a function x ˜ ( · ) B r , for some t 0 J ,
( F x ˜ ) C n γ : = t 0 n γ | ( F x ˜ ) ( t 0 ) | > r .
It follows from (10) that
r < t 0 n γ | ( F x ˜ ) ( t 0 ) | j = 1 n c j t 0 n j Γ ( γ j + 1 ) + Λ 1 T α p 2 r + Λ 2 T α p 1 + ( 1 σ ) ( n γ ) r σ .
Dividing both sides of the above inequality by r and taking the limit as r , we deduce that Λ 1 T α p 2 1 . This is a contradiction. Hence, F B r B r .
Step II. We show that F is continuous. Let { x m } be a sequence such that x m x in B r , then there exists ε > 0 such that x m x C n γ < ε for n sufficiently large. By (9), we have
| f ( t , x m ( t ) ) f ( t , x ( t ) ) | t σ ( γ n ) μ ( t ) ( ε σ + 2 r σ ) .
Moreover, by (H1), for almost every t J , we see that | f ( t , x m ( t ) ) f ( t , x ( t ) ) | 0 . From the Lebesgue’s dominated convergence theorem, it follows that | ( F x m ) ( t ) ( F x ) ( t ) | 0 as m . Now we see that F is continuous.
Step III. We prove that F maps bounded sets into equicontinuous sets of B r . Let t 1 , t 2 J , t 1 < t 2 and x B r . Similar to Lemma 7, we derive
t 1 n γ | ( I 0 + α λ ( · ) x ) ( t 2 ) ( I 0 + α λ ( · ) x ) ( t 1 ) | + | ( I 0 + α F x ) ( t 2 ) ( I 0 + α F x ) ( t 1 ) | t 1 n γ Γ ( α ) [ 0 t 1 | ( t 1 s ) α 1 ( t 2 s ) α 1 | s γ n | λ ( s ) | r + s σ ( n γ ) | μ ( s ) | r σ d s + t 1 t 2 ( t 2 s ) α 1 s γ n | λ ( s ) | r + s σ ( n γ ) | μ ( s ) | r σ d s ] 0 , as t 1 t 2 .
Then by (10) and (11),
| t 2 n γ ( F x ) ( t 2 ) t 1 n γ ( F x ) ( t 1 ) | | t 2 n γ t 1 n γ | | ( F x ) ( t 2 ) | + t 1 n γ | ( F x ) ( t 2 ) ( F x ) ( t 1 ) | | t 2 n γ t 1 n γ | | ( F x ) ( t 2 ) | + t 1 n γ j = 1 n c j | t 2 γ j t 1 γ j | Γ ( γ j + 1 ) + t 1 n γ [ | ( I 0 + α λ ( · ) x ) ( t 2 ) ( I 0 + α λ ( · ) x ) ( t 1 ) | + | ( I 0 + α F x ) ( t 2 ) ( I 0 + α F x ) ( t 1 ) | ] 0 , as t 1 t 2 .
From the above steps we deduce that F is completely continuous. It follows from Schauder’s fixed point theorem that F has a fixed point x B r . □

5. Applications

In this section, we present two examples to illustrate our results.
Example 1.
We consider the following initial value problem for the variable parameter fractional differential equation:
H D 0 + 9 10 , 1 5 x ( t ) sin t 2 t 5 x ( t ) = cos t t 5 | x ( t ) | 1 2 , t J : = ( 0 , 1 ] , ( I 0 + 2 25 x ) ( 0 ) = c 2 , ( I 0 + 2 25 x ) ( 0 ) = c 1 .
Set α = 9 10 , β = 1 5 , n = 1 , γ = α + β α β = 23 25 , n γ = 2 25 , λ ( t ) = sin t 2 t 5 , f ( t , x ( t ) ) = cos t t 5 | x ( t ) | 1 2 . Notice that this problem is a particular case of (1).
We choose p = 1 4 ( 2 25 , 9 20 ) , then for x ( t ) C 2 25 ,
0 1 | cos t t 5 | x ( t ) | 1 2 | 4 d t 1 4 0 1 1 t 24 25 d t 1 4 · x C 2 25 1 2 < + ,
i.e., f satisfies (H1). Moreover, | f ( t , x ( t ) ) | μ ( t ) | x ( t ) | 1 2 where μ ( t ) = 1 t 5 L 1 p 1 ( p 1 = 1 3 ) , μ L 3 = ( 5 2 ) 1 3 and λ L 1 p 2 ( p 2 = 1 4 ) , λ L 4 5 1 4 2 . By simple calculation one can obtain
Λ 1 : = B α p 2 1 p 2 , γ n + 1 p 2 1 p 2 1 p 2 λ L 1 p 2 Γ ( α ) B 13 15 , 67 75 3 4 5 1 4 2 Γ ( 9 10 ) 0.84 < 1 .
By Theorem 6, the problem (12) has at least one solution x C 2 25 .
Example 2.
We consider the following initial value problem for the variable parameter fractional differential equation:
H D 0 + 19 10 , 1 10 x ( t ) λ t ν x ( t ) = b t ω , t J : = ( 0 , 1 ] , ( I 0 + 9 100 x ) ( 0 + ) = c 2 , ( I 0 + 9 100 x ) ( 0 + ) = c 1 ,
where λ , b R and ν > 1 100 , ω > 0 .
Set α = 19 10 , β = 1 10 , n = 2 , γ = α + 2 β α β = 191 100 , n γ = 9 100 , λ ( t ) = λ t ν , f ( t , x ( t ) ) = λ t ν x ( t ) + b t ω . Notice that this problem is a particular case of (5).
We choose p = 1 10 ( 9 100 , 9 20 ) , then for x ( t ) C 9 100 ,
0 1 | λ t ν x ( t ) | 10 d t 1 10 | λ | 0 1 t 10 ν 9 10 d t 1 10 x C 9 100 = | λ | 10 ν + 1 10 1 10 x C 9 100 < + ,
i.e., f ( · , x ( · ) ) L 10 . Moreover,
| f ( t , x ( t ) ) f ( t , y ( t ) ) | | λ | t ν | x ( t ) y ( t ) | : = μ ( t ) | x ( t ) y ( t ) | , μ ( t ) L 10
and M f : = sup t J | f ( t , 0 ) | = | b | . Thus, by Theorem 4, the problem (13) has a unique solution x C 9 100 .
Let
x 0 ( t ) = j = 1 2 c j t γ j Γ ( γ j + 1 ) , x m ( t ) = x 0 ( t ) + ( I 0 + α λ ( · ) x m 1 ) ( t ) + b I 0 + α t ω ,
iterating, for m N , we can write
x ( t ) = j = 1 2 k = 0 c j ( I 0 + α λ ( · ) ) k t γ j Γ ( γ j + 1 ) + k = 1 ( I 0 + α λ ( · ) ) k 1 ( b I 0 + α t ω ) .
Noting that
( I 0 + α λ ( · ) ) t θ = λ I 0 + α t ν + θ = λ Γ ( ν + θ + 1 ) Γ ( α + ν + θ + 1 ) t α + ν + θ ,
( I 0 + α λ ( · ) ) 2 t θ = λ 2 Γ ( ν + θ + 1 ) Γ ( α + ν + θ + 1 ) I 0 + α t α + 2 ν + θ = λ 2 Γ ( ν + θ + 1 ) Γ ( α + 2 ν + θ + 1 ) Γ ( α + ν + θ + 1 ) Γ ( 2 ( α + ν ) + θ + 1 ) t 2 ( α + ν ) + θ ,
we obtain
( I 0 + α λ ( · ) ) k t θ = i = 0 k 1 Γ [ i α + ( i + 1 ) ν + θ + 1 ] Γ [ ( i + 1 ) ( α + ν ) + θ + 1 ] ( λ t α + ν ) k · t θ = i = 0 k 1 Γ α i 1 + ν α + ν + θ α + 1 Γ α i 1 + ν α + ν + θ α + 1 + 1 ( λ t α + ν ) k · t θ , k N .
Now we have
j = 1 2 k = 0 c j ( I 0 + α λ ( · ) ) k t γ j Γ ( γ j + 1 ) = j = 1 2 c j t γ j Γ ( γ j + 1 ) E α , 1 + ν α , ν + γ j α ( λ t α + ν ) ,
where
E α , ζ , l ( z ) = k = 0 d k z k , z R , α > 0 , ζ > 0 , α ( i ζ + l ) 1 , 2 ,
is a generalized Mittag–Leffler function ([4])and d 0 = 1 ,
d k = i = 0 k 1 Γ α ( i ζ + l ) + 1 Γ α ( i ζ + l + 1 ) + 1 .
Similarly, we have
k = 1 ( I 0 + α λ ( · ) ) k 1 t α + ω = t α + ω E α , 1 + ν α , ν + α + ω α ( λ t α + ν ) .
Finally, we get the following explicit solution for (13)
x ( t ) = j = 1 2 c j t γ j Γ ( γ j + 1 ) E α , 1 + ν α , ν + γ j α ( λ t α + ν ) + b t α + ω Γ ( ω + 1 ) Γ ( α + ω + 1 ) E α , 1 + ν α , ν + α + ω α ( λ t α + ν ) .

6. Conclusions

In this paper, using new techniques, the existence and uniqueness results for differential equations of Hilfer-type fractional order with variable coefficient are successfully obtained. The results are new even for the special case: β = 0 or β = 1 . The results hold for all α ( n 1 , n ) ( n N ) . The proposed techniques can be extended to other hybrid fractional differential equations.

7. Future Research

In this paper, the variable coefficient λ ( · ) L 1 p and is not necessarily continuous, therefore, the technique can be used to solve other types of equations. For example, inspired by the reference [19], as future work, we will study the fractional differential equation with variable coefficients with respect to another function under the weak assumptions for σ i , g :
H D 0 + β 0 , μ 0 , φ y ( t ) + i = 1 m σ i ( t ) H D 0 + β i , μ i , φ y ( t ) = g ( t ) , t [ 0 , t 0 ] , m N ,
where β 0 > β 1 > > β m 0 , 0 μ i 1 and n i are non-negative integers satisfying n i 1 < β i < n i , i = 0 , 1 , , m .

Author Contributions

Conceptualization, F.L. and H.W.; methodology, F.L.; software, C.W.; validation, F.L., C.W. and H.W.; formal analysis, F.L.; investigation, C.W.; resources, F.L.; writing—original draft preparation, F.L.; writing—review and editing, C.W. and H.W.; supervision, F.L.; project administration, F.L. and H.W.; funding acquisition, F.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge the financial support provided by the Natural Science Foundation of China (11971329,11561077).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley: New York, NY, USA, 1993. [Google Scholar]
  2. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach Science Publishers: Yverdon, Switzerland, 1993. [Google Scholar]
  3. Podlubny, I. Fractional Differential Equations, Vol. 198 of Mathematics in Science and Engineering; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  4. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematics Studies, 204; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  5. Li, F.; Liang, J.; Xu, H.K. Existence of mild solutions for fractional integrodifferential equations of Sobolev type with nonlocal conditions. J. Math. Anal. Appl. 2012, 391, 510–525. [Google Scholar] [CrossRef] [Green Version]
  6. Liang, J.; Mu, Y. Mild solutions to the Cauchy problem for some fractional differential equations with delay. Axioms 2017, 6, 30. [Google Scholar] [CrossRef] [Green Version]
  7. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  8. Hilfer, R. Threefold Introduction to Fractional Derivatives: In Anomalous Transport: Foundations and Applications; Wiley-VCH Verlag GmbH & Co. KGaA: Weinheim, Germany, 2008. [Google Scholar]
  9. Hilfer, R.; Luchko, Y.; Tomovski, Z. Operational method for the solution of fractional differential equations with generalized Riemann-Lioville fractional derivative. Fract. Calc. Appl. Anal. 2009, 12, 289–318. [Google Scholar]
  10. 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]
  11. Furati, K.M.; Kassim, M.D.; Tatar, N.E. Non-existence of global solutions for a differential equation involving Hilfer fractional derivative. Electron. J. Differ. Equ. 2013, 235, 1–10. [Google Scholar]
  12. Gu, H.; Trujillo, J.J. Existence of mild solution for evolution equation with Hilfer fractional derivative. Appl. Math. Comput. 2014, 257, 344–354. [Google Scholar] [CrossRef]
  13. Wang, J.; Ahmed, H.M. Null controllability of nonlocal Hilfer fractional stochastic differential equations. Miskolc Math. Notes 2017, 18, 1073–1083. [Google Scholar]
  14. Liang, J.; Mu, Y.; Xiao, T.-J. Initial-value/Nonlocal Cauchy problems for fractional differential equations involving ψ-Hilfer multivariable operator. Fract. Calc. Appl. Anal. 2020, 23, 1090–1124. [Google Scholar] [CrossRef]
  15. Hu, Y.; Luo, Y.; Lu, Z. Analytical Solution of the linear fractional differential equation by Adomian decomposition method. J. Comput. Appl. Math. 2008, 215, 220–229. [Google Scholar] [CrossRef] [Green Version]
  16. Morita, T.; Sato, K. Neumann-Series Solution of Fractional Differential Equation. Interdiscip. Inf. Sci. 2010, 16, 127–137. [Google Scholar] [CrossRef] [Green Version]
  17. Kilbas, A.A.; Rivero, M.; Rodrignez-Germa, L.; Trujillo, J.J. α-Analytic solutions of some linear fractional differential equations with variable coefficients. Appl. Math. Comput. 2007, 187, 239–249. [Google Scholar] [CrossRef]
  18. Vivek, D.; Kanagarajan, K.; Elsayed, E.M. Some Existence and Stability Results for Hilfer-fractional Implicit Differential Equations with Nonlocal Conditions. Mediterr. J. Math. 2018, 15, 1–21. [Google Scholar] [CrossRef]
  19. Restrepo, J.E.; Suragan, D. Hilfer-type fractional differential equations with variable coefficients. Chaos Solitons Fractals 2021, 150, 111146. [Google Scholar] [CrossRef]
  20. Henry, D. Geometric Theory of Semilinear Parabolic Partial Differential Equations; Springer: Berlin, Germany, 1989. [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

Li, F.; Wang, C.; Wang, H. Existence Results for Hilfer Fractional Differential Equations with Variable Coefficient. Fractal Fract. 2022, 6, 11. https://doi.org/10.3390/fractalfract6010011

AMA Style

Li F, Wang C, Wang H. Existence Results for Hilfer Fractional Differential Equations with Variable Coefficient. Fractal and Fractional. 2022; 6(1):11. https://doi.org/10.3390/fractalfract6010011

Chicago/Turabian Style

Li, Fang, Chenglong Wang, and Huiwen Wang. 2022. "Existence Results for Hilfer Fractional Differential Equations with Variable Coefficient" Fractal and Fractional 6, no. 1: 11. https://doi.org/10.3390/fractalfract6010011

APA Style

Li, F., Wang, C., & Wang, H. (2022). Existence Results for Hilfer Fractional Differential Equations with Variable Coefficient. Fractal and Fractional, 6(1), 11. https://doi.org/10.3390/fractalfract6010011

Article Metrics

Back to TopTop