Next Article in Journal
Three-Dimensional Hydro-Magnetic Flow Arising in a Long Porous Slider and a Circular Porous Slider with Velocity Slip
Next Article in Special Issue
Difference Mappings Associated with Nonsymmetric Monotone Types of Fejér’s Inequality
Previous Article in Journal
A Modified Fletcher–Reeves Conjugate Gradient Method for Monotone Nonlinear Equations with Some Applications
Previous Article in Special Issue
Refinements of Majorization Inequality Involving Convex Functions via Taylor’s Theorem with Mean Value form of the Remainder
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Gronwall–Bellman Inequality in Frame of Generalized Proportional Fractional Derivative

1
Department of Mathematics and General Sciences, Prince Sultan University, P.O. Box 66833, 11586 Riyadh, Saudi Arabia
2
Department of General Education, Faculty of Science and Health Technology, Navamindradhiraj University, Bangkok 10300, Thailand
3
Department of Mathematics, Van Yuzuncu Yil University, 65080 Van, Turkey
4
Department of Mathematics, Faculty of Mathematics, University of Sistan and Baluchestan, Zahedan, Iran
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(8), 747; https://doi.org/10.3390/math7080747
Submission received: 19 July 2019 / Revised: 10 August 2019 / Accepted: 12 August 2019 / Published: 15 August 2019
(This article belongs to the Special Issue Inequalities)

Abstract

:
New versions of a Gronwall–Bellman inequality in the frame of the generalized (Riemann–Liouville and Caputo) proportional fractional derivative are provided. Before proceeding to the main results, we define the generalized Riemann–Liouville and Caputo proportional fractional derivatives and integrals and expose some of their features. We prove our main result in light of some efficient comparison analyses. The Gronwall–Bellman inequality in the case of weighted function is also obtained. By the help of the new proposed inequalities, examples of Riemann–Liouville and Caputo proportional fractional initial value problems are presented to emphasize the solution dependence on the initial data and on the right-hand side.

1. Introduction

Integral inequalities have been used as fabulous instruments to explore the qualitative properties of differential equations [1]. Over the years, there have appeared many inequalities which have been established by many authors such as Ostrowski type inequality, Hardy type inequality, Olsen type inequality, Gagliardo–Nirenberg type inequality, Lyapunove type inequality, Opial type inequality and Hermite–Hadamard type inequality [2,3]. However, the most common and significant inequality is the Gronwall–Bellman inequality, which they introduced in [4,5]. The Gronwall–Bellman inequality allows one to provide an estimate for a function that is known to satisfy a certain integral inequality by the solution of the corresponding integral equation. In particular, it has been employed to provide a comparison that can be used to prove uniqueness of a solution to an initial value problem (see some recent relevant papers [6,7,8,9]).
Fractional differential equations (FDEs) is a rich area of research that has widespread applications in science and engineering. Indeed, it describes a large number of nonlinear phenomena in different fields such as physics, chemistry, biology, viscoelasticity, control hypothesis, speculation, fluid dynamics, hydrodynamics, aerodynamics, information processing system networking, notable and picture processing, control theory, etc. FDEs also provide marvellous tools for the depiction of memory and inherited properties of many materials and processes. In view of recent developments, one can consequently conclude that FDEs have emerged significant achievements in the last couple of decades [10,11,12,13,14,15,16]. The study of integral equations in the scope of non-integer-order equations has been in the spotlight in the recent years. Many mathematicians in the field of applied and pure mathematics have dedicated their efforts to extend, generalize and refine the integral inequalities carried over from integer order equations to the non-integer order equations. Meanwhile, different definitions of fractional derivatives have been recently introduced [17,18]. The Gronwall–Bellman inequality, which is our concern herein, has been under investigation and different versions of it have been established for different types of fractional operators [19,20,21,22,23,24,25].
In this paper, new versions for a Gronwall–Bellman inequality in the frame of the newly defined generalized (Riemann–Liouville and Caputo) proportional fractional derivative are provided. Before proceeding to the main results, we define the generalized Riemann–Liouville and Caputo proportional fractional derivatives and integrals and expose some of their features [26]. We prove our main result in light of some efficient comparison analysis. The Gronwall–Bellman inequality in the case of a weighted function is also obtained. By the help of the new proposed inequalities, examples of Riemann–Liouville and Caputo generalized proportional fractional initial value problems are presented to emphasize the solution dependence on the initial data and on the right-hand side. It worth mentioning that the new proposed derivative is well-behaved. Indeed, it has nonlocal character and satisfies the semigroup or the so-called index property. Besides, the resulting inequalities converge to the classical ones upon considering particular cases of the derivative. That is, our results not only extend the classical inequalities but also generalize the existing ones for non-integer-order equations.

2. The GPF Derivatives and Integrals

We assemble in this section some fundamental preliminaries that are used throughout the remaining part of the paper. For their justifications and proofs, the reader can consult the work in [26].
In control theory, a proportional derivative controller (PDC) for controller output u at time t with two tuning parameters has the algorithm
u ( t ) = κ p E ( t ) + κ d d d t E ( t ) ,
where κ p is the proportional gain, κ d is the derivative gain, and E is the input deviation or the error between the state variable and the process variable. Recent investigations have shown that PDC has direct incorporation in the control of complex networks models (see [27] for more details).
For ρ [ 0 , 1 ] , let the functions κ 0 , κ 1 : [ 0 , 1 ] × R [ 0 , ) be continuous such that for all t R we have
lim ρ 0 + κ 1 ( ρ , t ) = 1 , lim ρ 0 + κ 0 ( ρ , t ) = 0 , lim ρ 1 κ 1 ( ρ , t ) = 0 , lim ρ 1 κ 0 ( ρ , t ) = 1 ,
and κ 1 ( ρ , t ) 0 , ρ [ 0 , 1 ) , κ 0 ( ρ , t ) 0 , ρ ( 0 , 1 ] . Then, Anderson et al. [28] defined the proportional derivative of order ρ by
D ρ ξ ( t ) = κ 1 ( ρ , t ) ξ ( t ) + κ 0 ( ρ , t ) ξ ( t )
provided that the right-hand side exists at t R and ξ : = d d t ξ . For the operator given in Equation (1), κ 1 is a type of proportional gain κ p , κ 0 is a type of derivative gain κ d , ξ is the error and u = D ρ ξ is the controller output. The reader can consult the work in [29] for more details about the control theory of the proportional derivative and its component functions. We only consider here the case when κ 1 ( ρ , t ) = 1 ρ and κ 0 ( ρ , t ) = ρ . Therefore, Equation (1) becomes
D ρ ξ ( t ) = ( 1 ρ ) ξ ( t ) + ρ ξ ( t ) .
It is easy to find that lim ρ 0 + D ρ ξ ( t ) = ξ ( t ) and lim ρ 1 D ρ ξ ( t ) = ξ ( t ) . Thus, the derivative in Equation (2) is somehow more general than the conformable derivative, which certainly does not converge to the original functions as ρ tends to 0.
In what follows, we define the generalized proportional fractional (GPF) integral and derivative:
Definition 1
([26]). For 0 < ρ 1 , α C and R e ( α ) > 0 , the GPF integral of ξ of order α is
( a I α , ρ ξ ) ( t ) = 1 ρ α Γ ( α ) a t e ρ 1 ρ ( t τ ) ( t τ ) α 1 ξ ( τ ) d τ = ρ α e ρ 1 ρ t a I α e 1 ρ ρ t ξ ( t ) .
Definition 2
([26]). For 0 < ρ 1 , α C , R e ( α ) 0 and n = [ R e ( α ) ] + 1 . Then, the Riemann–Liouville type GPF derivative of f of order α is
( a D α , ρ ξ ) ( t ) = D n , ρ a I n α , ρ ξ ( t ) = D t n , ρ ρ n α Γ ( n α ) a t e ρ 1 ρ ( t τ ) ( t τ ) n α 1 ξ ( τ ) d τ .
Remark 1.
If we let ρ = 1 in Definition 2, then one can obtain the left Riemann–Liouville fractional derivative [12,14,15]. Moreover, it is obvious that
lim α 0 ( D α , ρ ξ ) ( t ) = ξ ( t ) a n d lim α 1 ( D α , ρ ξ ) ( t ) = ( D ρ ξ ) ( t ) .
Proposition 1
([26]). Let α , β C be such that R e ( α ) 0 and R e ( β ) > 0 . Then, for any 0 < ρ 1 , we have
(1)
a I α , ρ e ρ 1 ρ t ( t a ) β 1 ( x ) = Γ ( β ) Γ ( β + α ) ρ α e ρ 1 ρ x ( x a ) α + β 1 , R e ( α ) > 0 .
(2)
a D α , ρ e ρ 1 ρ t ( t a ) β 1 ( x ) = ρ α Γ ( β ) Γ ( β α ) e ρ 1 ρ x ( x a ) β 1 α , R e ( α ) 0 .
In the following lemmas, we expose some features of Riemann–Liouville type GPF operator. The first result concerns with the index property of GPF which is of great significance.
Lemma 1
([26]). If 0 < ρ 1 , R e ( α ) > 0 and R e ( β ) > 0 . For a continuous function ξ defined on [ a , ) , we have
a I α , ρ ( a I β , ρ ξ ) ( t ) = a I β , ρ ( a I α , ρ ξ ) ( t ) = ( a I α + β , ρ ξ ) ( t ) .
The action of the operator a D α , ρ on the integral operator is demonstrated in the following results.
Lemma 2
([26]). Let 0 < ρ 1 , 0 m < [ R e ( α ) ] + 1 and ξ be integrable in each interval [ a , t ] , t > a . Then,
a D m , ρ ( a I α , ρ ξ ) ( t ) = ( a I α m , ρ ξ ) ( t ) .
Corollary 1
([26]). Let 0 < ρ 1 , 0 < R e ( β ) < R e ( α ) and m 1 < R e ( β ) m . Then, we have
a D β , ρ ( a I α , ρ ξ ) ( t ) = ( a I α β , ρ ξ ) ( t ) .
Lemma 3
([26]). Let f be integrable on t a and R e [ α ] > 0 , 0 < ρ 1 , n = [ R e ( α ) ] + 1 . Then, we have
a D α , ρ ( a I α , ρ ξ ) ( t ) = ξ ( t ) .
Lemma 4
([26]). Let 0 < ρ 1 , R e ( α ) > 0 , n = [ R e ( α ) ] + 1 , ξ L 1 ( a , b ) and ( a I α , ρ ξ ) ( t ) A C n [ a , b ] . Then,
a I α , ρ ( a D α , ρ ξ ) ( t ) = ξ ( t ) e ρ 1 ρ ( t a ) j = 1 n ( a I j α , ρ ξ ) ( a + ) ( t a ) α j ρ α j Γ ( α + 1 j ) .
The GPF derivative of Caputo type is defined as follows:
Definition 3
([26]). For 0 < ρ 1 , α C , R e ( α ) 0 and n = [ R e ( α ) ] + 1 . Then, the GPF derivative of Caputo type of ξ of order α is
( a C D α , ρ ξ ) ( t ) = a I n α , ρ ( D n , ρ ξ ) ( t ) = 1 ρ n α Γ ( n α ) a t e ρ 1 ρ ( t τ ) ( t τ ) n α 1 ( D n , ρ ξ ) ( τ ) d τ .
Proposition 2
([26]). Let α , β C be such that R e ( α ) > 0 and R e ( β ) > 0 . Then, for any 0 < ρ 1 and n = [ R e ( α ) ] + 1 , we have
a C D α , ρ e ρ 1 ρ t ( t a ) β 1 ( x ) = ρ α Γ ( β ) Γ ( β α ) e ρ 1 ρ x ( x a ) β 1 α , R e ( β ) > n .
For k = 0 , 1 , , n 1 , we have a C D α , ρ e ρ 1 ρ t ( t a ) k ( x ) = 0 .
Lemma 5
([26]). For ρ ( 0 , 1 ] , R e ( α ) > 0 and n = [ R e ( α ) ] + 1 . Then, we have
a I α , ρ ( a C D α , ρ ξ ) ( t ) = ξ ( t ) e ρ 1 ρ ( t a ) k = 0 n 1 ( a D k , ρ ξ ) ( a ) ρ k k ! ( t a ) k .

3. Main Results

This section is devoted to provide our main results of this paper. We formulate new versions of the Gronwall–Bellman inequality within GPF operators in Riemann–Liouville and Caputo settings.

3.1. Gronwall–Bellman Inequality via the GPF Derivative of Riemann–Liouville Type

Consider the following generalized proportional Riemann–Liouville fractional initial value problem
a D α , ρ y ( t ) = f ( t , y ( t ) ) , 0 < α 1 , t [ a , b ] , lim t a + a I 1 α , ρ y ( t ) = y ( a ) = y a .
Applying the operator a I α , ρ to both sides of Equation (10), we obtain
y ( t ) = e ρ 1 ρ ( t a ) ( t a ) α 1 y ( a ) + a I α , ρ f ( t , y ( t ) ) ,
In the following, we present a comparison result for the GPF integral operator.
Theorem 1.
Let η and ζ be nonnegative continuous functions defined on [ a , b ] and satisfying
η ( t ) e ρ 1 ρ ( t a ) ( t a ) α 1 η ( a ) + a I α , ρ f ( t , η ( t ) ) ,
and
ζ ( t ) e ρ 1 ρ ( t a ) ( t a ) α 1 ζ ( a ) + a I α , ρ f ( t , ζ ( t ) ) ,
respectively. Suppose further that f satisfies a one-sided Lipschitz condition of the form
f ( t , x ) f ( t , y ) L e ρ 1 ρ ( a t ) e ρ 1 ρ ( t a ) ( t a ) α 1 + ( t a ) α α + 1 ( x y ) , f o r x y , L > 0 ,
and f ( t , y ) is nondecreasing in y. Then, η ( a ) ζ ( a ) and L < 1 + α ( t a ) α Γ ( α ) ρ α e ρ 1 ρ ( t a ) imply that η ( t ) ζ ( t ) for all t [ a , b ] .
Proof. 
We start by setting
η ε ( t ) = η ( t ) + ε e ρ 1 ρ ( t a ) ( t a ) α 1 + ( t a ) α α + 1 , for   small ε > 0 ,
so that we have
η ε ( a ) = η ( a ) + ε > η ( a ) and η ε ( t ) > η ( t ) , t [ a , b ] .
It follows that
η ε ( t ) e ρ 1 ρ ( t a ) ( t a ) α 1 η ( a ) + 1 Γ ( α ) ρ α a t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , η ( s ) ) d s + ε e ρ 1 ρ ( t a ) ( t a ) α 1 + ( t a ) α α + 1
or
η ε ( t ) e ρ 1 ρ ( t a ) ( t a ) α 1 η ( a ) + 1 Γ ( α ) ρ α a t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , η ( s ) ) d s 1 Γ ( α ) ρ α a t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , η ε ( s ) ) d s + 1 Γ ( α ) ρ α a t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , η ε ( s ) ) d s + ε e ρ 1 ρ ( t a ) ( t a ) α 1 + ( t a ) α α + 1 .
Using the Lipschitz condition in Equation (14) and the relations in Equations (15) and (16), we obtain
η ε ( t ) e ρ 1 ρ ( t a ) ( t a ) α 1 η ε ( a ) ε L Γ ( α ) ρ α a t ( t s ) α 1 d s + 1 Γ ( α ) ρ α a t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , η ε ( s ) ) d s + ε ( t a ) α α + 1
Since a t ( t s ) α 1 d s = ( t a ) α α and L < 1 + α ( t a ) α Γ ( α ) ρ α e ρ 1 ρ ( t + a ) , we arrive at
η ε ( t ) > e ρ 1 ρ ( t a ) ( t a ) α 1 η ε ( a ) + 1 Γ ( α ) ρ α a t e ρ 1 ρ ( t s ) ( t s ) α 1 f ( s , η ε ( s ) ) d s .
The remaining part of the proof can be completed by adopting the same steps followed in the proof of Theorem 2.1 in [30,31] to get η ε ( t ) ζ ( t ) , t [ a , b ] . However, and since ε is arbitrary, we conclude that η ( t ) ζ ( t ) , t [ a , b ] holds true. □
Remark 2.
The Lipschitz condition in Equation (14) can be relaxed by relaxing the upper bound for the constant L.
For our purpose, we replace f ( t , y ( t ) ) in Equation (11) by x ( t ) y ( t ) where | x ( t ) | < 1 , t [ a , b ] . Define the following operator
Ω x ϕ = a I α , ρ x ( t ) ϕ ( t ) .
The following results are important in the proof of the main theorem. We only state these lemmas as their proofs are straightforward.
Lemma 6.
For any constant λ, one has
| Ω λ e ρ 1 ρ ( t a ) ( t a ) α 1 | Ω | λ | e ρ 1 ρ ( t a ) ( t a ) α 1 .
Lemma 7.
For any constant λ, one has
| Ω λ n e ρ 1 ρ ( t a ) ( t a ) α 1 | = | λ | n ( t a ) ( n + 1 ) α 1 Γ ( α ) ρ n α Γ ( ( n + 1 ) α ) e ρ 1 ρ ( t a ) , n = 0 , 1 , 2 , .
Lemma 8.
Let λ > 0 be such that | y ( t ) | λ for t [ a , b ] . Then,
| Ω y n e ρ 1 ρ ( t a ) ( t a ) α 1 | Ω λ n e ρ 1 ρ ( t a ) ( t a ) α 1 , n = 0 , 1 , 2 , .
Theorems 1 and 2 together give us the desired proportional Riemann–Liouville fractional Gronwall–Bellman-type inequality.
Theorem 2.
Let y be a nonnegative function on [ a , b ] . Then, the GPF integral equation
y ( t ) = e ρ 1 ρ ( t a ) ( t a ) α 1 y ( a ) + a I α , ρ x ( t ) y ( t ) , t [ a , b ] ,
has a solution
y ( t ) = y ( a ) k = 0 Ω x k e ρ 1 ρ ( t a ) ( t a ) α 1 .
Proof. 
The proof is accomplished by applying the successive approximation method. Set
y 0 ( t ) = e ρ 1 ρ ( t a ) ( t a ) α 1 y ( a )
and
y n ( t ) = e ρ 1 ρ ( t a ) ( t a ) α 1 y ( a ) + a I α , ρ x ( t ) y n 1 ( t ) , n 1 .
We observe that
y 1 ( t ) = e ρ 1 ρ ( t a ) ( t a ) α 1 y ( a ) + a I α , ρ x ( t ) y 0 ( t ) = y ( a ) Ω x 0 e ρ 1 ρ ( t a ) ( t a ) α 1 + y ( a ) Ω x 1 e ρ 1 ρ ( t a ) ( t a ) α 1 ,
and
y 2 ( t ) = e ρ 1 ρ ( t a ) ( t a ) α 1 y ( a ) + a I α , ρ x ( t ) y 1 ( t ) = y ( a ) Ω x 0 e ρ 1 ρ ( t a ) ( t a ) α 1 + Ω x 1 y ( a ) Ω x 0 e ρ 1 ρ ( t a ) ( t a ) α 1 + y ( a ) Ω x 1 e ρ 1 ρ ( t a ) ( t a ) α 1 = y ( a ) Ω x 0 e ρ 1 ρ ( t a ) ( t a ) α 1 + y ( a ) Ω x 1 e ρ 1 ρ ( t a ) ( t a ) α 1 + y ( a ) Ω x 2 e ρ 1 ρ ( t a ) ( t a ) α 1 .
It follows inductively that
y n ( t ) = y ( a ) k = 0 n Ω x k e ρ 1 ρ ( t a ) ( t a ) α 1 , n 0 .
Formally, taking the limit as n to obtain
y ( t ) = y ( a ) k = 0 Ω x k e ρ 1 ρ ( t a ) ( t a ) α 1 .
We use Lemmas 6–8, the comparison test and the d’Alembert ratio test to show the absolute convergence of the series in Equation (24). Indeed, the infinite series
n = 0 λ n ( t a ) ( n + 1 ) α 1 Γ ( α ) ρ n α Γ ( ( n + 1 ) α ) e ρ 1 ρ ( t a ) ,
is convergent for all t [ a , b ] and for all 0 < λ , ρ 1 . Let a n be defined as
a n = λ n ( t a ) ( n + 1 ) α 1 Γ ( α ) ρ n α Γ ( ( n + 1 ) α ) e ρ 1 ρ ( t a ) .
Then, we have
lim n a n + 1 a n = λ ( t a ) α ρ α lim n Γ ( ( n + 1 ) α ) Γ ( ( n + 2 ) α ) .
Next, we use Stirling approximation formula for the Gamma function x Γ ( x ) 2 π x x e x , where x is large enough. It is a straightforward computation using this formula to show that
lim x x Γ ( x ) 2 π x x e x = 1 and lim x x x + 1 x = 1 e ,
which are all we need. Hence, we have
lim n ( n + 1 ) α Γ ( ( n + 1 ) α ) 2 π ( n + 1 ) α ( n + 1 ) α e ( n + 1 ) α = 1 and lim n ( n + 2 ) α Γ ( ( n + 2 ) α ) 2 π ( n + 2 ) α ( n + 2 ) α e ( n + 2 ) α = 1 .
Thus,
lim n a n + 1 a n = λ ( t a ) α ρ α lim n Γ ( ( n + 1 ) α ) Γ ( ( n + 2 ) α ) = λ ( t a ) α ρ α lim n n + 2 n + 1 n + 1 n + 2 α e α α n + 1 n + 2 n α 1 n + 2 α = λ ( t a ) α ρ α e α α 1 e α 0 = 0 < 1 .
Hence, convergence is guaranteed. Besides, one can easily show that Equation (22) solves Equation (21). □
Remark 3.
Note that Equation (22) solves the inequality
ζ ( t ) e ρ 1 ρ ( t a ) ( t a ) α 1 ζ ( a ) + a I α , ρ ζ ( t ) y ( t ) , t [ a , b ] ,
where ζ and y are nonnegative real valued functions such that 0 y ( t ) < λ < 1 .
Now, we are in a position to state the main theorem, which is a new version of the Gronwall–Bellman inequality within the generalized proportional fractional Riemann–Liouville settings.
Corollary 2.
Let ζ and y be nonnegative real valued functions such that 0 y ( t ) < λ < 1 and
ζ ( t ) e ρ 1 ρ ( t a ) ( t a ) α 1 ζ ( a ) + a I α , ρ ζ ( t ) y ( t ) , t [ a , b ] .
Then,
ζ ( t ) ζ ( a ) k = 0 Ω y k e ρ 1 ρ ( t a ) ( t a ) α 1 .
The proof of the corollary is a straightforward implementation of Theorems 1 and 2. Indeed, it is immediately obtained by setting η ( t ) = ζ ( a ) k = 0 Ω y k e ρ 1 ρ ( t a ) ( t a ) α 1 .

3.2. Gronwall–Bellman Inequality via the GPF Derivative of Caputo Type

Consider the following generalized proportional Caputo fractional initial value problem
a C D α , ρ y ( t ) = f ( t , y ( t ) ) , 0 < α 1 , t [ a , b ] , y ( a ) = y a .
Applying the operator a I α , ρ to both sides of Equation (29), we obtain
y ( t ) = e ρ 1 ρ ( t a ) y ( a ) + a I α , ρ f ( t , y ( t ) ) ,
The results of this subsection resemble the ones proved in Section 3.1. To avoid redundancy, therefore, we skip some steps of the proofs. We start by the following comparison result for the generalized proportional Caputo fractional integral operator.
Theorem 3.
Let η and ζ be nonnegative continuous functions defined on [ a , b ] and satisfy
η ( t ) e ρ 1 ρ ( t a ) η ( a ) + a I α , ρ f ( t , η ( t ) ) ,
and
ζ ( t ) e ρ 1 ρ ( t a ) ζ ( a ) + a I α , ρ f ( t , ζ ( t ) ) ,
respectively. Suppose further that f satisfies one-sided Lipschitz condition of the form
f ( t , x ) f ( t , y ) L e ρ 1 ρ ( a t ) e ρ 1 ρ ( t a ) + ( t a ) α α ( x y ) , f o r x y , L > 0 ,
and f ( t , y ) is nondecreasing in y. Then, η ( a ) ζ ( a ) and L < Γ ( α ) ρ α e ρ 1 ρ ( t a ) imply that η ( t ) ζ ( t ) for all t [ a , b ] .
The proof of the above theorem can be completed by setting η ε ( t ) = η ( t ) + ε e ρ 1 ρ ( t a ) + ( t a ) α α , for small ε > 0 , and following similar steps as the proof of Theorem 1.
In the sequel, we replace f ( t , y ( t ) ) in Equation (30) by x ( t ) y ( t ) , where | x ( t ) | < 1 , t [ a , b ] . Define the following operator
Φ x ϕ = a I α , ρ x ( t ) ϕ ( t ) .
In similar manner, the following lemmas are formulated for Caputo type operator.
Lemma 9.
For any constant λ, one has
| Φ λ e ρ 1 ρ ( t a ) | Φ | λ | e ρ 1 ρ ( t a ) .
Lemma 10.
For any constant λ, one has
| Φ λ n e ρ 1 ρ ( t a ) | = | λ | n ( t a ) n α ρ n α Γ ( n α + 1 ) e ρ 1 ρ ( t a ) , n = 0 , 1 , 2 , .
Lemma 11.
Let λ > 0 be such that | y ( t ) | λ for t [ a , b ] . Then,
| Φ y n e ρ 1 ρ ( t a ) | = Φ λ n e ρ 1 ρ ( t a ) , n = 0 , 1 , 2 , .
Theorem 4.
Let y be a nonnegative function on [ a , b ] . Then, the generalized proportional fractional integral equation
y ( t ) = e ρ 1 ρ ( t a ) y ( a ) + a I α , ρ x ( t ) y ( t ) , t [ a , b ] ,
has a solution
y ( t ) = y ( a ) k = 0 Φ x k e ρ 1 ρ ( t a ) .
Proof. 
We employ the successive approximation method to complete the proof. Set
y 0 ( t ) = e ρ 1 ρ ( t a ) y ( a ) y n ( t ) = e ρ 1 ρ ( t a ) y ( a ) + a I α , ρ x ( t ) y n 1 ( t ) , n 1 .
We observe that
y 1 ( t ) = y ( a ) Φ x 0 e ρ 1 ρ ( t a ) + y ( a ) Φ x 1 e ρ 1 ρ ( t a )
and
y 2 ( t ) = e ρ 1 ρ ( t a ) y ( a ) + a I α , ρ x ( t ) y 1 ( t ) = y ( a ) Φ x 0 e ρ 1 ρ ( t a ) + y ( a ) Φ x 1 e ρ 1 ρ ( t a ) + y ( a ) Φ x 2 e ρ 1 ρ ( t a ) .
It follows inductively that y n ( t ) = y ( a ) k = 0 n Φ x k e ρ 1 ρ ( t a ) . Taking the limit as n to obtain
y ( t ) = y ( a ) k = 0 Φ x k e ρ 1 ρ ( t a ) .
Following the same arguments as in the proof of Theorem 2, we use Lemmas 9–11, the comparison test and the d’Alembert ratio test to show the absolute convergence of the series in Equation (40). Moreover, it is clear to verify that Equation (39) solves Equation (38). The proof is finished. □
Remark 4.
Note that Equation (39) solves the inequality
ζ ( t ) e ρ 1 ρ ( t a ) ζ ( a ) + a I α , ρ ζ ( t ) y ( t ) , t [ a , b ] ,
where ζ and y are nonnegative functions on [ a , b ] such that 0 y ( t ) < λ < 1 .
The Gronwall–Bellman inequality in generalized proportional Caputo fractional is stated as follows.
Corollary 3.
Let ζ and y be nonnegative real valued functions such that 0 y ( t ) < λ < 1 and
ζ ( t ) e ρ 1 ρ ( t a ) ζ ( a ) + a I α , ρ ζ ( t ) y ( t ) , t [ a , b ] .
Then,
ζ ( t ) ζ ( a ) k = 0 Φ y k e ρ 1 ρ ( t a ) .
To prove Equation (43), we set η ( t ) = ζ ( a ) k = 0 Φ y k e ρ 1 ρ ( t a ) and the rest follows as a direct application of Theorems 3 and 4.

4. Gronwall–Bellman Inequality via Weighted Function

In this section, we extend the results obtained in Section 3 to the case of weighted function. The analysis can be carried out for the Riemann–Liouville and Caputo operators. However, we only present the results for the case of Riemann–Liouville proportional fractional operator. Unlike previous relevant results in the literature [32], the weighted function w in the following first two theorems requires no monotonic restriction.
Theorem 5.
Let η , ζ , w be nonnegative continuous functions on [ a , b ] where η and ζ satisfy
η ( t ) e ρ 1 ρ ( t a ) ( t a ) α 1 η ( a ) + w ( t ) a I α , ρ f ( t , η ( t ) ) ,
and
ζ ( t ) e ρ 1 ρ ( t a ) ( t a ) α 1 ζ ( a ) + w ( t ) a I α , ρ f ( t , ζ ( t ) ) ,
respectively. Suppose further that f satisfies one-sided Lipschitz condition of the form
f ( t , x ) f ( t , y ) L e ρ 1 ρ ( a t ) e ρ 1 ρ ( t a ) ( t a ) α 1 + w ( t ) ( t a ) α α + 1 ( x y ) , f o r x y , L > 0 ,
and f ( t , y ) is nondecreasing in y. Then, η ( a ) ζ ( a ) and L < 1 + α w ( t ) ( t a ) α Γ ( α ) ρ α e ρ 1 ρ ( t a ) imply that η ( t ) ζ ( t ) for all t [ a , b ] .
To prove the above theorem, we set η ε ( t ) = η ( t ) + ε e ρ 1 ρ ( t a ) + w ( t ) ( t a ) α α + 1 , f o r s m a l l ε > 0 , and follow similar steps as the proof of Theorem 1.
Remark 5.
The Lipschitz condition in Equation (46) can be relaxed by relaxing the upper bound for the constant L.
Theorem 6.
Let x , y be nonnegative functions on [ a , b ] and w be a nonnegative continuous function defined on [ a , b ] . Further, assume that | x ( t ) | < 1 for t [ a , b ] and max t [ a , b ] w ( t ) = M . Then, the generalized proportional fractional integral equation
y ( t ) = e ρ 1 ρ ( t a ) ( t a ) α 1 y ( a ) + w ( t ) a I α , ρ x ( t ) y ( t ) , t [ a , b ] ,
has a solution
y ( t ) = y ( a ) Ω x 0 e ρ 1 ρ ( t a ) ( t a ) α 1 + y ( a ) w ( t ) k = 1 M k 1 Ω x k e ρ 1 ρ ( t a ) ( t a ) α 1 .
Remark 6.
Note that Equation (48) solves the inequality
ζ ( t ) e ρ 1 ρ ( t a ) ( t a ) α 1 ζ ( a ) + w ( t ) a I α , ρ ζ ( t ) y ( t ) , t [ a , b ] ,
where ζ , y are nonnegative functions on [ a , b ] and w is a nonnegative continuous function defined on [ a , b ] and 0 y ( t ) < λ < 1 and max t [ a , b ] w ( t ) = M .
The Gronwall–Bellman inequality in case of weighted function w is stated as follows.
Theorem 7.
Let ζ , y be nonnegative functions on [ a , b ] and w be a nonnegative continuous function defined on [ a , b ] . Further, assume that 0 y ( t ) < λ < 1 for t [ a , b ] and max t [ a , b ] w ( t ) = M and
ζ ( t ) e ρ 1 ρ ( t a ) ( t a ) α 1 ζ ( a ) + w ( t ) a I α , ρ ζ ( t ) y ( t ) , t [ a , b ] .
Then,
ζ ( t ) ζ ( a ) Ω y 0 e ρ 1 ρ ( t a ) ( t a ) α 1 + ζ ( a ) w ( t ) k = 1 M k 1 Ω y k e ρ 1 ρ ( t a ) ( t a ) α 1 .
If the weighted function w possesses a monotonic behavior, then Theorem 6 and Theorem 7 can be reformulated, respectively, in the following forms.
Theorem 8.
Let y , x be nonnegative functions on [ a , b ] and w be a nonnegative continuous function defined on [ a , b ] . Further, assume that | x ( t ) | < 1 for t [ a , b ] and w is a nondecreasing function. Then, the generalized proportional fractional integral equation
y ( t ) = e ρ 1 ρ ( t a ) ( t a ) α 1 y ( a ) + w ( t ) a I α , ρ x ( t ) y ( t ) , t [ a , b ] ,
has a solution
y ( t ) = y ( a ) k = 0 w k ( t ) Ω x k e ρ 1 ρ ( t a ) ( t a ) α 1 .
Theorem 9.
Let ζ , y be nonnegative functions on [ a , b ] and w be a nonnegative continuous function defined on [ a , b ] . Assume that 0 y ( t ) < λ < 1 for t [ a , b ] and w is a nondecreasing function and
ζ ( t ) e ρ 1 ρ ( t a ) ( t a ) α 1 ζ ( a ) + w ( t ) a I α , ρ ζ ( t ) y ( t ) , t [ a , b ] .
Then,
ζ ( t ) ζ ( a ) k = 0 w k ( t ) Ω y k e ρ 1 ρ ( t a ) ( t a ) α 1 .

5. Applications

In this section, two examples of Riemann–Liouville and Caputo generalized proportional fractional initial value problems are presented. By the help of the new proposed Gronwall–Bellman inequalities in Theorems 2 and 3, we show that the solution of the initial value problems depend on the initial data and on the right-hand side.
Consider the proportional Riemann–Liouville fractional initial value problem in Equation (10). In the remaining part of this section, we assume that the nonlinearity function f ( t , y ) satisfies a Lipschitz condition with a constant L [ 0 , 1 ) for all ( t , y ) .
Example 1.
Consider the following Riemann–Liouville proportional fractional initial value problems of the form
( a D α , ρ β ) ( t ) = f ( t , β ( t ) ) , lim t a + a I 1 α , ρ β ( t ) = β ( a ) = β 0 , 0 < α 1 , t [ a , b ] ,
and
( a D α , ρ γ ) ( t ) = f ( t , γ ( t ) ) , lim t a + a I 1 α , ρ γ ( t ) = γ ( a ) = γ 0 , 0 < α 1 , t [ a , b ] .
We claim that a small change in the initial condition implies a small change in the solution.
Proof. 
Applying the generalized proportional fractional integral operator in Equations (56) and (57), we have
β ( t ) = e ρ 1 ρ ( t a ) ( t a ) α 1 β 0 + a I α , ρ f ( t , β ( t ) ) ,
and
γ ( t ) = e ρ 1 ρ ( t a ) ( t a ) α 1 γ 0 + a I α , ρ f ( t , γ ( t ) ) .
It follows that
β ( t ) γ ( t ) = e ρ 1 ρ ( t a ) ( t a ) α 1 β 0 γ 0 + a I α , ρ [ f ( t , β ( t ) ) f ( t , γ ( t ) ) ] .
Taking the absolute value, we obtain
| β ( t ) γ ( t ) | e ρ 1 ρ ( t a ) ( t a ) α 1 | β 0 γ 0 | + a I α , ρ | f ( t , β ( t ) ) f ( t , γ ( t ) ) | e ρ 1 ρ ( t a ) ( t a ) α 1 | β 0 γ 0 | + L a I α , ρ | β ( t ) γ ( t ) | .
By employing Theorem 2, we get
| β ( t ) γ ( t ) | | β 0 γ 0 | k = 0 Ω L k e ρ 1 ρ ( t a ) ( t a ) α 1 = | β 0 γ 0 | k = 0 L k ( t a ) ( k + 1 ) α 1 Γ ( α ) ρ k α Γ ( ( k + 1 ) α ) e ρ 1 ρ ( t a ) .
Consider the initial value problem
( a D α , ρ ν ) ( t ) = f ( t , ν ( t ) ) , 0 < α 1 , t [ a , b ] lim t a + a I 1 α , ρ ν ( t ) = ν ( a ) = β n ,
where β n β 0 . If the solution of Equation (58) is denoted by ν n , then, for all t [ a , b ] , we have
| β ( t ) ν n ( t ) | | β 0 β n | k = 0 L k ( t a ) ( k + 1 ) α 1 Γ ( α ) ρ k α Γ ( ( k + 1 ) α ) e ρ 1 ρ ( t a ) .
Hence, | β ( t ) ν n ( t ) | 0 when β n β 0 as n . We conclude that a small change in the initial condition implies a small change in the solution. □
Example 2.
Consider the following Caputo generalized proportional fractional initial value problems of the form
( a C D α , ρ β ) ( t ) = f ( t , β ( t ) ) , β ( a ) = β 0 , α ( 0 , 1 ] , t [ a , b ] .
and
( a C D α , ρ σ ) ( t ) = f ( t , σ ( t ) ) + g ( t , σ ( t ) ) , σ ( a ) = σ 0 , α ( 0 , 1 ] , t [ a , b ] .
We claim that the solution of Equation (60) depends continuously on the right-hand side of Equation (60) if | g ( t , σ ) | K e ρ 1 ρ ( t a ) for all t [ a , b ] and for a positive number K.
Proof. 
If the solution of Equation (60) is denoted by σ , then, for all t [ a , b ] , we have
| β ( t ) σ ( t ) | e ρ 1 ρ ( t a ) | β 0 σ 0 | + a I α , ρ | f ( t , β ( t ) ) f ( t , σ ( t ) ) | + a I α , ρ | g ( t , σ ( t ) ) | e ρ 1 ρ ( t a ) | β 0 σ 0 | + L a I α , ρ | β ( t ) σ ( t ) | + a I α , ρ | g ( t , σ ( t ) ) | .
By the assumption, we have
| β ( t ) σ ( t ) | e ρ 1 ρ ( t a ) | β 0 σ 0 | + L a I α , ρ | β ( t ) σ ( t ) | + a I α , ρ K e ρ 1 ρ ( t a )
or
| β ( t ) σ ( t ) | + K L e ρ 1 ρ ( t a ) e ρ 1 ρ ( t a ) | β 0 σ 0 | + K L + L a I α , ρ | β ( t ) σ ( t ) | + K L e ρ 1 ρ ( t a ) .
Let r ( t ) = | β ( t ) σ ( t ) | + K L e ρ 1 ρ ( t a ) . Then, if we apply Theorem 3, we obtain
r ( t ) | β 0 σ 0 | + K L k = 0 Φ L k e ρ 1 ρ ( t a ) ,
or
| β ( t ) σ ( t ) | | β 0 σ 0 | + K L k = 0 L k ( t a ) ( k + 1 ) α 1 ρ k α Γ ( ( k + 1 ) α ) e ρ 1 ρ ( t a ) K L e ρ 1 ρ ( t a ) .
For a t b , letting K e ρ 1 ρ ( t a ) < δ implies that
| β ( t ) σ ( t ) | | β 0 σ 0 | k = 0 L k ( t a ) ( k + 1 ) α 1 K ρ k α Γ ( ( k + 1 ) α ) δ + δ L k = 0 L k ( t a ) ( k + 1 ) α 1 ρ k α Γ ( ( k + 1 ) α ) 1 δ | β 0 σ 0 | k = 0 L k b ( k + 1 ) α 1 K ρ k α Γ ( ( k + 1 ) α ) + 1 L k = 0 L k b ( k + 1 ) α 1 ρ k α Γ ( ( k + 1 ) α ) 1 = ϵ ,
which implies that a small change on the right-hand side of Equation (59) implies a small change in its solution. □

6. Conclusions

One of the most crucial issues in the theory of differential equations is to study qualitative properties for solutions of these equations. Integral inequalities are significant instruments that facilitate exploring such properties. In this paper, we accommodate a newly defined generalized proportional fractional (GPF) derivative to establish new versions for the well–known Gronwall–Bellman inequality. We prove our results in the frame of GPF operators within the Riemann–Liouville and Caputo settings. The main results are also extended to the weighted function case. One can easily figure out that the current results generalize the ones previously obtained in the literature. Indeed, the case ρ = 1 covers the results of classical Riemann–Liouville and Caputo fractional derivatives. As an application, we provide two efficient examples that demonstrate the solution dependence on the initial data and on the right-hand side of the initial value problems. The results of this paper have strong potential to be used for establishing new substantial investigations in the future for equations involving the GPF operators.

Author Contributions

All authors contributed equally and significantly to this paper. All authors have read and approved the final version of the manuscript.

Funding

The first author would like to thank Prince Sultan University for funding this work through research group Nonlinear Analysis Methods in Applied Mathematics (NAMAM) group number RG-DES-2017-01-17.

Acknowledgments

The authors would like to express their sincere thanks to the handling editor and the referees for their constructive comments and suggestions. We believe that their efforts have significantly helped in improving the contents of this paper.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Hardy, G.H.; Littlewood, J.E.; Pólya, G. Inequalities; Cambridge University Press: London, UK, 1952. [Google Scholar]
  2. Bainov, D.D.; Simeonov, P.S. Integral Inequalities and Applications (Mathematics and Its Applications); Springer: Berlin, Germany, 1992. [Google Scholar]
  3. Cloud, J.M.; Drachman, C.B.; Lebedev, P.L. Inequalities with Applications to Engineering; Springer: Berlin, Germany, 2014. [Google Scholar]
  4. Grönwall-Bellman, T.H. Note on the derivatives with respect to a parameter of the solutions of a system of differential equations. Ann. Math. 1919, 20, 292–296. [Google Scholar] [CrossRef]
  5. Bellman, R. The stability of solutions of linear differential equations. Duke Math. J. 1943, 10, 643–647. [Google Scholar] [CrossRef]
  6. Rasmussen, D.L. Gronwall’s inequality for functions of two independent variables. J. Math. Anal. Appl. 1976, 55, 407–417. [Google Scholar] [CrossRef]
  7. Dragomir, S.S. Some Gronwall Type Inequalities and Applications; Nova Science Pub Inc.: Hauppauge, NY, USA, 2003. [Google Scholar]
  8. Lin, X. A note on Gronwall’s inequality on time scales. Abstr. Appl. Anal. 2014, 2014, 623726. [Google Scholar] [CrossRef]
  9. Wang, W.; Feng, Y.; Wang, Y. Nonlinear Gronwall—Bellman type inequalities and their applications. Mathematics 2017, 5, 31. [Google Scholar] [CrossRef]
  10. Hilfer, R. Applications of Fractional Calculus in Physics; Word Scientific: Singapore, 2000. [Google Scholar]
  11. Debnath, L. Recent applications of fractional calculus to science and engineering. Int. J. Math. Math. Sci. 2003, 54, 3413–3442. [Google Scholar] [CrossRef]
  12. Kilbas, A.; Srivastava, H.M.; Trujillo, J.J. Theory and Application of Fractional Differential Equations. North-Holland Mathematical Studies; Elsevier (North-Holland) Science Publishers: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  13. Magin, R.L. Fractional Calculus in Bioengineering; Begell House: Redding, CT, USA, 2006. [Google Scholar]
  14. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach Science: Yverdon, Switzerland, 1993. [Google Scholar]
  15. Podlubny, I. Fractional Differential Equations. Mathematics in Science and Engineering 198; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  16. He, J.H.; Ji, F.Y. Two–scale mathematics and fractional calculus for thermodynamics. Therm. Sci. 2019. [Google Scholar] [CrossRef]
  17. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Progr. Fract. Differ. Appl. 2015, 1, 73–85. [Google Scholar]
  18. Atangana, A.; Baleanu, D. New fractional derivatives with non-local and non-singular kernels. Therm. Sci. 2016, 20, 757–763. [Google Scholar] [CrossRef]
  19. Ye, H.; Gao, J.; Ding, Y. A generalized Gronwall inequality and its application to a fractional differential equation. J. Math. Anal. Appl. 2007, 328, 1075–1081. [Google Scholar] [CrossRef] [Green Version]
  20. Ferreira, R.A.C. A Discrete fractional Gronwall inequality. Proc. Am. Math. Soc. 2012, 5, 1605–1612. [Google Scholar] [CrossRef]
  21. Abdeljawad, T.; Alzabut, J. The q-fractional analogue for Gronwall-type inequality. J. Funct. Spaces Appl. 2013. [Google Scholar] [CrossRef]
  22. Abdeljawad, T.; Alzabut, J.; Baleanu, D. A generalized q-fractional Gronwall inequality and its applications to nonlinear delay q-fractional difference systems. J. Inequal. Appl. 2016, 240, 1–13. [Google Scholar] [CrossRef]
  23. Zhang, Z.; Wei, Z. A generalized Gronwall inequality and its application to fractional neutral evolution inclusions. J. Inequal. Appl. 2016, 45, 1–18. [Google Scholar] [CrossRef]
  24. Sarikaya, M.Z. Gronwall type inequalities for conformable fractional integrals. Konuralp J. Math. 2016, 4, 217–222. [Google Scholar]
  25. Alzabut, J.; Abdeljawad, T. A generalized discrete fractional Gronwall inequality and its application on the uniqueness of solutions for nonlinear delay fractional difference system. Appl. Anal. Discret. Math. 2018, 12, 36–48. [Google Scholar] [CrossRef]
  26. Jarad, F.; Abdeljawad, T.; Alzabut, J. Generalized fractional derivatives generated by a class of local proportional derivatives. Eur. Phys. J. Spec. Top. 2017, 226, 3457–3471. [Google Scholar] [CrossRef]
  27. Ding, D.; Zhang, X.; Cao, J.; Wang, N.; Liang, D. Bifurcation control of complex networks model via PD controller. Neurocomputing 2016, 175, 1–9. [Google Scholar]
  28. Anderson, D.R.; Ulness, D.J. Newly defined conformable derivatives. Adv. Dyn. Syst. Appl. 2015, 10, 109–137. [Google Scholar]
  29. Anderson, D.R. Second–order self-adjoint differential equations using a proportional–derivative controller. Commun. Appl. Nonlinear Anal. 2017, 24, 17–48. [Google Scholar]
  30. Lakshmikantham, V.; Vatsala, A.S. Basic theory of fractional differential equations. Nonlinear Anal. 2008, 69, 2677–2682. [Google Scholar] [CrossRef]
  31. Denton, Z.; Vatsala, A.S. Fractional integral inequalities and applications. Comput. Math. Appl. 2010, 59, 1087–1094. [Google Scholar] [CrossRef] [Green Version]
  32. Alzabut, J.; Abdeljawad, T.; Jarad, F.; Sudsutad, W. A Gronwall inequality via the generalized proportional fractional derivative with applications. J. Inequal. Appl. 2019, 101, 1–12. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Alzabut, J.; Sudsutad, W.; Kayar, Z.; Baghani, H. A New Gronwall–Bellman Inequality in Frame of Generalized Proportional Fractional Derivative. Mathematics 2019, 7, 747. https://doi.org/10.3390/math7080747

AMA Style

Alzabut J, Sudsutad W, Kayar Z, Baghani H. A New Gronwall–Bellman Inequality in Frame of Generalized Proportional Fractional Derivative. Mathematics. 2019; 7(8):747. https://doi.org/10.3390/math7080747

Chicago/Turabian Style

Alzabut, Jehad, Weerawat Sudsutad, Zeynep Kayar, and Hamid Baghani. 2019. "A New Gronwall–Bellman Inequality in Frame of Generalized Proportional Fractional Derivative" Mathematics 7, no. 8: 747. https://doi.org/10.3390/math7080747

APA Style

Alzabut, J., Sudsutad, W., Kayar, Z., & Baghani, H. (2019). A New Gronwall–Bellman Inequality in Frame of Generalized Proportional Fractional Derivative. Mathematics, 7(8), 747. https://doi.org/10.3390/math7080747

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