Next Article in Journal
Efficient Generators of the Generalized Fractional Gaussian Noise and Cauchy Processes
Next Article in Special Issue
Local Stability, Global Stability, and Simulations in a Fractional Discrete Glycolysis Reaction–Diffusion Model
Previous Article in Journal
Local Error Estimate of an L1-Finite Difference Scheme for the Multiterm Two-Dimensional Time-Fractional Reaction–Diffusion Equation with Robin Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lyapunov-Type Inequalities for Systems of Riemann-Liouville Fractional Differential Equations with Multi-Point Coupled Boundary Conditions

College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(6), 454; https://doi.org/10.3390/fractalfract7060454
Submission received: 8 May 2023 / Revised: 27 May 2023 / Accepted: 29 May 2023 / Published: 1 June 2023
(This article belongs to the Special Issue Stability Analysis for Fractional-Order Equations)

Abstract

:
We consider a system of Riemann–Liouville fractional differential equations with multi-point coupled boundary conditions. Using some techniques from matrix analysis and the properties of the integral operator defined on two Banach spaces, we establish some Lyapunov-type inequalities for the problem considered. Moreover, the comparison between two Lyapunov-type inequalities is given under certain special conditions. The inequalities obtained compliment the existing results in the literature.

1. Introduction

We are interested in the systems of nonlinear fractional differential equations with multi-point coupled boundary conditions
D a + β 1 u ( t ) + f ( t , u ( t ) , v ( t ) ) = 0 , t ( a , b ) , D a + β 2 v ( t ) + g ( t , u ( t ) , v ( t ) ) = 0 , t ( a , b ) , u ( a ) = 0 , u ( b ) = i = 1 n a 1 i u ( ξ i ) + j = 1 n a 2 j v ( η j ) , v ( a ) = 0 , v ( b ) = i = 1 n a 3 i u ( ξ i ) + j = 1 n a 4 j v ( η j ) ,
where a , b R , 0 < a < b , a < ξ 1 < ξ 2 < < ξ n < b , a < η 1 < η 2 < < η n < b , a i j 0 ( i = 1 , 2 , 3 , 4 ; j = 1 , 2 , , n ) , 1 < β i 2 ( i = 1 , 2 ) , D a + β i ( i = 1 , 2 ) is the Riemann–Liouville fractional derivative, and f , g : [ a , b ] × R 2 R are given functions. Using the spectral radius of the matrix, we establish Lyapunov-type inequality for (1). The well-known Lyapunov inequality [1] shows that a necessary condition for the second-order linear differential equation
y ( t ) + q ( t ) y ( t ) = 0 , a < t < b , u ( a ) = u ( b ) = 0
to have nontrivial solutions is that
a b | q ( s ) | d s > 4 b a .
Lyapunov inequality has found many practical applications such as estimates for intervals of disconjugacy [2] and eigenvalue problems [3] in investigating the qualitative properties of solutions of differential equations, marking the difference between Equations (4) and integral Equations (5). Since then, inequality (2) was rediscovered and generalized many times; see, for example [2,3,4,5].
The first Lyapunov-type inequalities for fractional boundary value problems is due to Ferreira [6], where he established the following result:
If a nontrivial solution to the Riemann–Liouville fractional boundary value problem
D a + α y ( t ) + q ( t ) y ( t ) = 0 , a < t < b , 1 < α < 2 , u ( a ) = u ( b ) = 0
exists, where q is a real and continuous function, then
a b | q ( s ) | d s > Γ ( α ) 4 b a α 1 .
By substituting the Riemann–Liouville fractional derivative D a + α in (3) by the Caputo fractional derivative C D a + α , a Lyapunov-type inequality [7] was obtained as follows:
If q C [ a , b ] , then
a b | q ( s ) | d s > Γ ( α ) α α [ ( α 1 ) ( b a ) ] α 1
holds if there is a nontrivial solution for the following Caputo boundary value problem
C D a + α y ( t ) + q ( t ) y ( t ) = 0 , a < t < b , 1 < α < 2 , u ( a ) = u ( b ) = 0 .
In [8], Ferreira addresses the issue of further research directions for Lyapunov-type inequality. After that, some results related to the study of Lyapunov-type inequalities for various types of fractional differential equations were obtained; see the survey article of Ntouyas et al. [9] and its complemented survey [10,11], and the papers of [12,13,14,15] and references therein. For example, in [12], a Lyapunov-type inequality was obtained for a higher order Riemann–Liouville fractional differential equation with fractional integral boundary conditions, and the lower bound for the eigenvalues of nonlocal boundary value problems was also presented. In [13], the authors proved a Lyapunov-type inequality for a class of Riemann–Liouville fractional boundary value problems with fractional boundary conditions. In [14], a Lyapunov-type inequality was obtained for Riemann–Liouville-type fractional boundary value problems with fractional boundary conditions
D a + α y ( t ) + q ( t ) y ( t ) = 0 , t ( a , b ) , y ( a ) = D a + β u ( b ) = 0 ,
where 1 < α 2 , 0 β 1 , q C [ a , b ] . It was proved that if (4) has a nontrivial solution, then
a b ( b s ) α β 1 | q ( s ) | d s > Γ ( α ) ( b a ) β 1 , if 0 < α β < 1 ,
and
a b | q ( s ) | d s > Γ ( α ) 2 2 α β 2 ( b a ) α 1 , if 1 α β < 2 .
In [15], Wang et al. considered multi-point boundary value problem of the form
D a + α u ( t ) + q ( t ) u ( t ) = 0 , t ( a , b ) , 2 < α < 3 , u ( a ) = u ( a ) = 0 , D a + β + 1 u ( b ) = i = 1 m 2 b i D a + β u ( ξ i ) ,
where D a + α is the Riemann–Liouville fractional derivative, a < ξ 1 < ξ 2 < < ξ m 2 < b , b i 0 , 0 δ = i = 1 m 2 b i ( ξ i a ) α β 1 < ( α β 2 ) ( b a ) α β 2 , and obtained the following Lyapunov-type inequalities:
a b ( b s ) α β 2 ( b a ) β + 1 ( b s ) β + 1 + ( b a ) α 1 δ i = 1 m 2 b i ( s a ) | q ( s ) | d s Γ ( α ) .
Recently, Jleli, O’Regan and Same [16] studied a coupled system of Caputo fractional differential equations
( C D a + α u ) ( t ) = f ( t , u , v ) , t ( a , b ) , ( C D a + β v ) ( t ) = g ( t , u , v ) , t ( a , b ) , u ( a ) = u ( b ) = 0 , v ( a ) = v ( b ) = 0 ,
where 1 < α , β < 2 . Let I a b ( θ , h ) = ( θ 1 ) θ 1 ( b a ) θ 1 θ θ Γ ( θ ) a b h ( s ) d s for h C [ a , b ] . They proved the following Lyapunov-type inequalities:
I a b ( α , p 11 ) + I a b ( β , p 22 ) + I a b ( α , p 11 ) I a b ( β , p 22 ) 2 + 4 I a b ( α , p 12 ) I a b ( β , p 21 ) 1 2 2 ,
if (5) has a nontrivial solution and there exist positive functions p i j C [ a , b ] ( i , j = 1 , 2 ) such that
| f ( t , u , v ) f ( t , w , z ) | p 11 ( t ) | u w | + p 12 ( t ) | v z | , t [ a , b ] , u , v , w , z R ,
and
| g ( t , u , v ) g ( t , w , z ) | p 21 ( t ) | u w | + p 22 ( t ) | v z | , t [ a , b ] , u , v , w , z R .
It is worth mentioning that, for the above fractional differential equation, the method used to obtain the fractional Lyapunov-type inequalities is the Green’s function approach that derives the Green’s function of the equivalent integral form of the boundary value problem being considered and then finding the maximum or an upper bound of its Green’s function. In addition, compared with a large number of references devoted to the study of Lyapunov-type inequalities for fractional differential equations, there is not much undertaken for systems of fractional differential equations. To the best of our knowledge, there is no paper to study the Lyapunov-type inequality for systems of Riemann–Liouville fractional differential equations with coupled boundary conditions. The objective of the present paper is to fill the gap in this area and, more extensively, to study the Lyapunov-type inequalities for the systems of nonlinear fractional differential equations with multi-point coupled boundary conditions (1). Coupled boundary conditions appear in the study of Sturm–Liouville problems and reaction-diffusion Equations [17], and have applications in many fields of sciences and engineering, such as thermal conduction [18] and mathematical biology [19]. The reader may consult the paper [20] for the initial study of differential equations under coupled boundary conditions and the paper [21,22] for fractional differential equations equipped with multi-point boundary conditions.
This paper has three significant features: the boundary conditions contain coupled multi-point boundary conditions and uncoupled multi-point boundary conditions; the establishment of three Lyapunov-type inequalities for (1) is mainly based on matrix analysis and the properties of the operator T defined on two Banach spaces E × E and E 1 × E 2 ; the comparison between two Lyapunov-type inequalities for (1) is given under certain special conditions.
The organization of this paper is as follows. In Section 2, we present preliminary definitions and properties of fractional calculus, several essential lemmas associated with this work. The main results and their proof are given in Section 3.
We make the following assumptions on the systems of nonlinear fractional differential equations (1):
(H0) a i j 0 ( i = 1 , 2 , 3 , 4 ; j = 1 , 2 , , n ) , κ i j 0 ( i , j = 1 , 2 ) and κ = Δ κ 11 κ 22 κ 12 κ 21 > 0 , where
κ 11 = 1 i = 1 n a 1 i ( ξ i a ) β 1 1 ( b a ) β 1 1 , κ 12 = j = 1 n a 2 j ( η j a ) β 2 1 ( b a ) β 2 1 ,
κ 21 = i = 1 n a 3 i ( ξ i a ) β 1 1 ( b a ) β 1 1 , κ 22 = 1 j = 1 n a 4 j ( η j a ) β 2 1 ( b a ) β 2 1 .
(H1) f , g : [ a , b ] × R 2 R are continuous.
(H2) There exist positive functions p 11 , p 12 C [ a , b ] such that
| f ( t , x , y ) | p 11 ( t ) | x | + p 12 ( t ) | y | , t [ a , b ] , x , y R .
(H3) There exist positive functions p 21 , p 22 C [ a , b ] such that
| g ( t , x , y ) | p 21 ( t ) | x | + p 22 ( t ) | y | , t [ a , b ] , x , y R .

2. Preliminaries

In this part, we first give some basic definitions, lemmas and theorems.
Definition 1 
([23,24]). The derivative with fractional order α > 0 of Riemann–Liouville type is defined for the function σ defined on [ a , b ] as
D a + α σ ( t ) = 1 Γ ( n α ) d d t ( n ) a t σ ( s ) ( t s ) α + 1 n d s , t [ a , b ] , n = [ α ] + 1 ,
where Γ ( n α ) is Euler gamma function.
Definition 2 
([23,24]). The integral with fractional order α > 0 of Riemann–Liouville type is defined for the function σ as
I a + α σ ( t ) = 1 Γ ( α ) a t ( t s ) α 1 σ ( s ) d s , t [ a , b ] ,
where Γ ( α ) is Euler gamma function.
Lemma 1 
([23,24]). Suppose that φ C [ a , b ] , n N + , and n 1 < α n . Thus, the general solution of D a + α u ( t ) = φ ( t ) is
u ( t ) = I a + α φ ( t ) + c 1 ( t a ) α 1 + c 2 ( t a ) α 2 + + c n ( t a ) α n ,
such that c j R , j = 1 , , n .
Let E = C [ a , b ] be the Banach space equipped with norm x = max t [ a , b ] | x ( t ) | for x E .
Consider the system of linear fractional differential equations
D a + β 1 u ( t ) + φ 1 ( t ) = 0 , t ( a , b ) , u ( a ) = 0 , u ( b ) = i = 1 n a 1 i u ( ξ i ) + j = 1 n a 2 j v ( η j ) , D a + β 2 v ( t ) + φ 2 ( t ) = 0 , t ( a , b ) , v ( a ) = 0 , v ( b ) = i = 1 n a 3 i u ( ξ i ) + j = 1 n a 4 j v ( η j ) ,
then one has the following lemma.
Lemma 2. 
Let φ 1 , φ 2 C [ a , b ] then ( u , v ) is a solution of (6) if and only if ( u , v ) is a solution of the integral equation
u ( t ) = a b G 11 ( t , s ) φ 1 ( s ) d s + a b G 12 ( t , s ) φ 2 ( s ) d s , v ( t ) = a b G 21 ( t , s ) φ 1 ( s ) d s + a b G 22 ( t , s ) φ 2 ( s ) d s ,
where
G 11 ( t , s ) = G β 1 ( t , s ) + ( t a ) β 1 1 κ ( b a ) β 1 1 i = 1 n ( κ 22 a 1 i + κ 12 a 3 i ) G β 1 ( ξ i , s ) ,
G 12 ( t , s ) = ( t a ) β 1 1 κ ( b a ) β 1 1 j = 1 n ( κ 22 a 2 j + κ 12 a 4 j ) G β 2 ( η j , s ) ,
G 21 ( t , s ) = ( t a ) β 2 1 κ ( b a ) β 2 1 i = 1 n ( κ 21 a 1 i + κ 11 a 3 i ) G β 1 ( ξ i , s ) ,
G 22 ( t , s ) = G β 2 ( t , s ) + ( t a ) β 2 1 κ ( b a ) β 2 1 j = 1 n ( κ 21 a 2 j + κ 11 a 4 j ) G β 2 ( η j , s ) ,
and
G β i ( t , s ) = 1 Γ ( β i ) ( t a ) β i 1 ( b a ) β i 1 ( b s ) β i 1 ( t s ) β i 1 , a s t b ; ( t a ) β i 1 ( b a ) β i 1 ( b s ) β i 1 , a t s b .
Proof. 
Deduced from Lemma 1, we obtain
u ( t ) = I a + β 1 φ 1 ( t ) + c 1 ( t a ) β 1 1 + c 2 ( t a ) β 1 2 , v ( t ) = I a + β 2 φ 2 ( t ) + d 1 ( t a ) β 2 1 + d 2 ( t a ) β 2 2 .
Therefore, the general solution of (6) is
u ( t ) = a t ( t s ) β 1 1 Γ ( β 1 ) φ 1 ( s ) d s + c 1 ( t a ) β 1 1 + c 2 ( t a ) β 1 2 ,
v ( t ) = a t ( t s ) β 2 1 Γ ( β 2 ) φ 2 ( s ) d s + d 1 ( t a ) β 2 1 + d 2 ( t a ) β 2 2 .
Since u ( a ) = v ( a ) = 0 , it is clear that c 2 = d 2 = 0 . Let u ( b ) = c and v ( b ) = d , so we conclude that
c = u ( b ) = a b ( b s ) β 1 1 Γ ( β 1 ) φ 1 ( s ) d s + c 1 ( b a ) β 1 1 ,
d = v ( b ) = a b ( b s ) β 2 1 Γ ( β 2 ) φ 2 ( s ) d s + d 1 ( b a ) β 2 1 .
Hence, (8) implies
u ( t ) = a b G β 1 ( t , s ) φ 1 ( s ) d s + c ( t a ) β 1 1 ( b a ) β 1 1 , v ( t ) = a b G β 2 ( t , s ) φ 2 ( s ) d s + d ( t a ) β 2 1 ( b a ) β 2 1 ,
where G β i ( t , s ) is given by (7).
It is worth noting that u ( 1 ) = c and v ( 1 ) = d . In order to determine c , d , we require that the function in (9) should satisfy multi-point boundary conditions in (6), i.e.,
c = u ( 1 ) = i = 1 n a 1 i a b G β 1 ( ξ i , s ) φ 1 ( s ) d s + j = 1 n a 2 j a b G β 2 ( η j , s ) φ 2 ( s ) d s + c i = 1 n a 1 i ( ξ i a ) β 1 1 ( b a ) β 1 1 + d j = 1 n a 2 j ( η j a ) β 2 1 ( b a ) β 2 1 ,
d = v ( 1 ) = i = 1 n a 3 i a b G β 1 ( ξ i , s ) φ 1 ( s ) d s + j = 1 n a 4 j a b G β 2 ( η j , s ) φ 2 ( s ) d s + c ( ξ i a ) β 1 1 ( b a ) β 1 1 i = 1 n a 3 i + d j = 1 n a 4 j ( η j a ) β 2 1 ( b a ) β 2 1 .
The above two equations are written in matrix form
κ 11 κ 12 κ 21 κ 22 c d = i = 1 n a 1 i a b G β 1 ( ξ i , s ) φ 1 ( s ) d s + j = 1 n a 2 j a b G β 2 ( η j , s ) φ 2 ( s ) d s i = 1 n a 3 i a b G β 1 ( ξ i , s ) φ 1 ( s ) d s + j = 1 n a 4 j a b G β 2 ( η j , s ) φ 2 ( s ) d s .
Since κ = κ 11 κ 22 κ 12 κ 21 0 , there is
c = κ 22 κ i = 1 n a 1 i a b G β 1 ( ξ i , s ) φ 1 ( s ) d s + κ 12 κ i = 1 n a 3 i a b G β 1 ( ξ i , s ) φ 1 ( s ) d s + κ 22 κ j = 1 n a 2 j a b G β 2 ( η j , s ) φ 2 ( s ) d s + κ 12 κ j = 1 n a 4 j a b G β 2 ( η j , s ) φ 2 ( s ) d s = 1 κ i = 1 n ( κ 22 a 1 i + κ 12 a 3 i ) a b G β 1 ( ξ i , s ) φ 1 ( s ) d s + 1 κ j = 1 n ( κ 22 a 2 j + κ 12 a 4 j ) a b G β 2 ( η j , s ) φ 2 ( s ) d s ,
and
d = κ 21 κ i = 1 n a 1 i a b G β 1 ( ξ i , s ) φ 1 ( s ) d s + κ 11 κ i = 1 n a 3 i a b G β 1 ( ξ i , s ) φ 1 ( s ) d s + κ 21 κ j = 1 n a 2 j a b G β 2 ( η j , s ) φ 2 ( s ) d s + κ 11 κ j = 1 n a 4 j a b G β 2 ( η j , s ) φ 2 ( s ) d s = 1 κ i = 1 n ( κ 21 a 1 i + κ 11 a 3 i ) a b G β 1 ( ξ i , s ) φ 1 ( s ) d s + 1 κ j = 1 n ( κ 21 a 2 j + κ 11 a 4 j ) a b G β 2 ( η j , s ) φ 2 ( s ) d s .
Hence,
u ( t ) = a b G β 1 ( t , s ) φ 1 ( s ) d s + ( t a ) β 1 1 κ ( b a ) β 1 1 i = 1 n ( κ 22 a 1 i + κ 12 a 3 i ) a b G β 1 ( ξ i , s ) φ 1 ( s ) d s + ( t a ) β 1 1 κ ( b a ) β 1 1 j = 1 n ( κ 22 a 2 j + κ 12 a 4 j ) a b G β 2 ( η j , s ) φ 2 ( s ) d s = a b G 11 ( t , s ) φ 1 ( s ) d s + a b G 12 ( t , s ) φ 2 ( s ) d s ,
v ( t ) = a b G β 2 ( t , s ) φ 2 ( s ) d s + ( t a ) β 2 1 κ ( b a ) β 2 1 i = 1 n ( κ 21 a 1 i + κ 11 a 3 i ) a b G β 1 ( ξ i , s ) φ 1 ( s ) d s + ( t a ) β 2 1 κ ( b a ) β 2 1 j = 1 n ( κ 21 a 2 j + κ 11 a 4 j ) a b G β 2 ( η j , s ) φ 2 ( s ) d s = a b G 21 ( t , s ) φ 1 ( s ) d s + a b G 22 ( t , s ) φ 2 ( s ) d s
and the proof is complete. □
Lemmas 3–5 below give some important properties of G β i ( t , s ) and G i j ( t , s ) . Parts (1)–(3) of Lemma 3 are taken from [6] and part (4) of Lemma 3 follows from the expression of G β i ( t , s ) .
Lemma 3 
([6]). The Green function G β i ( t , s ) defined above satisfies the following conditions:
(1) G β i ( t , s ) 0 for all t , s [ a , b ] .
(2) max t [ a , b ] G β i ( t , s ) = G β i ( s , s ) , s [ a , b ] .
(3) G β i ( s , s ) has a unique maximum, given by
max s [ a , b ] G β i ( s , s ) = G β i a + b 2 , a + b 2 = 1 Γ ( β i ) b a 4 β i 1 .
(4) G β i ( t , s ) 1 Γ ( β i ) ( t a ) β i 1 ( b a ) β i 1 ( b s ) β i 1 for all t , s [ a , b ] .
Lemma 4. 
For a < ξ i < b , we have
max s [ a , b ] G β i ( ξ i , s ) = G β i ( ξ i , ξ i ) = 1 Γ ( β i ) ( ξ i a ) β i 1 ( b a ) β i 1 ( b ξ i ) β i 1 .
Proof. 
Let us start to define two functions
g 1 ( s ) = ( ξ i a ) β i 1 ( b a ) β i 1 ( b s ) β i 1 ( ξ i s ) β i 1 , a s ξ i < b ,
and
g 2 ( s ) = ( ξ i a ) β i 1 ( b a ) β i 1 ( b s ) β i 1 , a < ξ i s b .
It is easy to see that g 2 ( s ) is a decreasing function on [ ξ i , b ] , and we have
g 2 ( s ) g 2 ( ξ i ) = G β i ( ξ i , ξ i ) = ( ξ i a ) β i 1 ( b a ) β i 1 ( b s ) β i 1 , s [ ξ i , b ] .
On the other hand, since 1 < β i 2 and s [ a , ξ i ] , we have
d g 1 ( s ) d s = ( β i 1 ) ( ξ i a ) β i 1 ( b a ) β i 1 ( b s ) β i 2 ( ξ i s ) β i 2 = ( β i 1 ) ( b s ) β i 2 ( ξ i a ) β i 1 ( b a ) β i 1 ( ξ i s ) β i 2 ( b s ) β i 2 0 .
Therefore, g 1 ( s ) is an increasing function on [ a , ξ i ] and g 1 ( s ) g 1 ( ξ i ) = G β i ( ξ i , ξ i ) . □
Lemma 5. 
Four functions G i j ( i , j = 1 , 2 ) defined in Lemma 2 satisfy the following conditions:
(i) G i j ( t , s ) λ i j for all a t , s b ,
(ii) G i j ( t , s ) μ i j ( t a ) β i 1 ( b s ) β j 1 for all a t , s b ,
where λ i j , μ i j ( i , j = 1 , 2 ) are given by
λ 11 = 1 Γ ( β 1 ) b a 4 β 1 1 + 1 κ Γ ( β 1 ) i = 1 n ( κ 22 a 1 i + κ 12 a 3 i ) ( ξ i a ) β 1 1 ( b a ) β 1 1 ( b ξ i ) β 1 1 ,
λ 12 = 1 κ Γ ( β 2 ) ( b a ) β 1 1 j = 1 n ( κ 22 a 2 j + κ 12 a 4 j ) ( η j a ) β 2 1 ( b η j ) β 2 1 ,
λ 21 = 1 κ Γ ( β 1 ) ( b a ) β 2 1 i = 1 n ( κ 21 a 1 i + κ 11 a 3 i ) ( ξ i a ) β 1 1 ( b ξ i ) β 1 1 ,
λ 22 = 1 Γ ( β 2 ) b a 4 β 2 1 + 1 κ Γ ( β 2 ) j = 1 n ( κ 21 a 2 j + κ 11 a 4 j ) ( η j a ) β 2 1 ( b a ) β 2 1 ( b η j ) β 2 1 ,
μ 11 = κ 22 κ Γ ( β 1 ) ( b a ) β 1 1 , μ 12 = 1 κ Γ ( β 2 ) ( b a ) β 1 1 ,
and
μ 21 = = 1 κ Γ ( β 1 ) ( b a ) β 2 1 , μ 22 = κ 11 κ Γ ( β 2 ) ( b a ) β 2 1 .
Proof. 
For sake of simplicity, we only prove Lemma 5 for function G 11 ( t , s ) . Similar arguments apply for the other function.
Using Lemmas 3 and 4, we obtain
G 11 ( t , s ) = G β 1 ( t , s ) + ( t a ) β 1 1 κ ( b a ) β 1 1 i = 1 n ( κ 22 a 1 i + κ 12 a 3 i ) G β 1 ( ξ i , s ) G β 1 ( s , s ) + 1 κ i = 1 n ( κ 22 a 1 i + κ 12 a 3 i ) G β 1 ( ξ i , ξ i ) 1 Γ ( β 1 ) b a 4 β 1 1 + 1 κ Γ ( β 1 ) i = 1 n ( κ 22 a 1 i + κ 12 a 3 i ) ( ξ i a ) β 1 1 ( b a ) β 1 1 ( b ξ i ) β 1 1 = λ 11 ,
and
G 11 ( t , s ) = G β 1 ( t , s ) + ( t a ) β 1 1 κ ( b a ) β 1 1 i = 1 n ( κ 22 a 1 i + κ 12 a 3 i ) G β 1 ( ξ i , s ) 1 Γ ( β 1 ) ( t a ) β 1 1 ( b a ) β 1 1 ( b s ) β 1 1 + ( t a ) β 1 1 κ ( b a ) β 1 1 i = 1 n ( κ 22 a 1 i + κ 12 a 3 i ) 1 Γ ( β 1 ) ( ξ i a ) β 1 1 ( b a ) β 1 1 ( b s ) β 1 1 = 1 Γ ( β 1 ) ( t a ) β 1 1 ( b a ) β 1 1 ( b s ) β 1 1 1 + κ 22 i = 1 n a 1 i + κ 12 i = 1 n a 3 i κ ( b a ) β 1 1 ( ξ i a ) β 1 1 = 1 Γ ( β 1 ) ( t a ) β 1 1 ( b a ) β 1 1 ( b s ) β 1 1 1 + κ 22 ( 1 κ 11 ) + κ 12 κ 21 κ = κ 22 κ Γ ( β 1 ) ( t a ) β 1 1 ( b a ) β 1 1 ( b s ) β 1 1 = μ 11 ( t a ) β 1 1 ( b s ) β 1 1 .
In the next section, nonnegative square matrices will be used in order to present Lyapunov-type inequalities for systems of nonlinear fractional differential Equation (1). In this proof, a key role will be played by the so-called convergent to zero matrices. A nonnegative square matrix M is said to be convergent to zero if
M k 0 , a s k .
Let us recall that a real square matrix M = ( b i j ) 2 × 2 is said to be nonnegative and write M 0 if b i j 0 for i , j = 1 , 2 . For square matrix M 1 , M 2 , we say M 1 M 2 if M 1 M 2 0 . Similar definitions and notation apply for vectors.
We denote by M 2 + the set of square nonnegative matrices. For a matrix M M 2 + , let Trace ( M ) , det ( M ) and ρ ( M ) denote the trace, the determinant and the spectral radius of M, respectively.
Lemma 6 
([25]). Let C M 2 + . If ρ ( C ) < 1 , then
lim n C n = 0 .
Lemma 7 
([16]). Let C M 2 + . Then
ρ ( C ) = T r a c e ( C ) + [ T r a c e ( C ) ] 2 4 d e t ( C ) 2 .
Lemma 8 
([26]). If M 1 M 2 0 , then ρ ( M 2 ) ρ ( M 1 ) .

3. Main Results

For p i j C [ a , b ] ( i , j = 1 , 2 ) , let
J i j ( p 1 j , p 2 j ) = λ i 1 a b p 1 j ( s ) d s + λ i 2 a b p 2 j ( s ) d s , i , j = 1 , 2 .
Theorem 1. 
Suppose that (H0)-(H3) are satisfied. If (1) have a nontrivial solution, then
J 11 ( p 11 , p 21 ) + J 22 ( p 12 , p 22 ) + [ J 11 ( p 11 , p 21 ) J 22 ( p 12 , p 22 ) ] 2 + 4 J 12 ( p 12 , p 22 ) J 21 ( p 11 , p 21 ) 2 .
Proof. 
Let ( u * , v * ) E × E be a nontrivial solution of (1), and suppose that
J 11 ( p 11 , p 21 ) + J 22 ( p 12 , p 22 ) + [ J 11 ( p 11 , p 21 ) J 22 ( p 12 , p 22 ) ] 2 + 4 J 12 ( p 12 , p 22 ) J 21 ( p 11 , p 21 ) < 2 .
Let us introduce the operator T : E × E E × E given by
T ( u , v ) = ( T 1 ( u , v ) , T 2 ( u , v ) ) , u , v E ,
where
T 1 ( u , v ) ( t ) = a b G 11 ( t , s ) f ( s , u ( s ) , v ( s ) ) d s + a b G 12 ( t , s ) g ( s , u ( s ) , v ( s ) ) d s ,
T 2 ( u , v ) ( t ) = a b G 21 ( t , s ) f ( s , u ( s ) , v ( s ) ) d s + a b G 22 ( t , s ) g ( s , u ( s ) , v ( s ) ) d s ,
and G i j ( i , j = 1 , 2 ) is defined in Lemma 2. By Lemma 2, ( u * , v * ) is a nontrivial fixed point of T.
Using (H2) and Lemma 5, for all t [ 0 , 1 ] , we obtain
| u * ( t ) | = | T 1 ( u * , v * ) ( t ) | a b G 11 ( t , s ) | f ( s , u * ( s ) , v * ( s ) ) | d s + a b G 12 ( t , s ) | g ( s , u * ( s ) , v * ( s ) ) | d s a b G 11 ( t , s ) ( p 11 ( s ) | u * ( s ) | + p 12 ( s ) | v * ( s ) | ) d s + a b G 12 ( t , s ) ( p 21 ( s ) | u * ( s ) | + p 22 ( s ) | v * ( s ) | ) d s λ 11 a b p 11 ( s ) d s + λ 12 a b p 21 ( s ) d s u * + λ 11 a b p 12 ( s ) d s + λ 12 a b p 22 ( s ) d s v * .
Therefore, we obtain
u * J 11 ( p 11 , p 21 ) u * + J 12 ( p 12 , p 22 ) v * .
Similarly, using Lemma 5 and (H3), we have
v * J 21 ( p 11 , p 21 ) u * + J 22 ( p 12 , p 22 ) v * .
Combining (12) with (13), we deduce that
u * v * J 11 ( p 11 , p 21 ) J 12 ( p 12 , p 22 ) J 21 ( p 11 , p 21 ) J 22 ( p 12 , p 22 ) u * v * .
By induction, for n N , we have
u * v * J 11 ( p 11 , p 21 ) J 12 ( p 12 , p 22 ) J 21 ( p 11 , p 21 ) J 22 ( p 12 , p 22 ) n u * v * .
Next, using Lemma 6, Lemma 7 and (11), we deduce that
u * = v * = 0 ,
which contradicts the nontriviality of ( u * , v * ) . This proves (10). □
For p i j C [ a , b ] ( i , j = 1 , 2 ) , let
I i j ( p 1 j , p 2 j ) = μ i 1 a b p 1 j ( s ) ( b s ) β 1 1 ( s a ) β j 1 d s + μ i 2 a b p 2 j ( s ) ( b s ) β 2 1 ( s a ) β j 1 d s , i , j = 1 , 2 .
Theorem 2. 
Suppose that (H0)-(H3) are satisfied. If (1) have a nontrivial solution, then
I 11 ( p 11 , p 21 ) + I 22 ( p 12 , p 22 ) + [ I 11 ( p 11 , p 21 ) I 22 ( p 12 , p 22 ) ] 2 + 4 I 12 ( p 12 , p 22 ) I 21 ( p 11 , p 21 ) 2 .
Proof. 
Let ( u * , v * ) E × E be a nontrivial solution of (1), and suppose that
I 11 ( p 11 , p 21 ) + I 22 ( p 12 , p 22 ) + [ I 11 ( p 11 , p 21 ) I 22 ( p 12 , p 22 ) ] 2 + 4 I 12 ( p 12 , p 22 ) I 21 ( p 11 , p 21 ) < 2 .
For u , v E , by the bounded property of continuous functions, there is K > 0 such that
| f ( s , u ( s ) , v ( s ) ) | K , | g ( s , u ( s ) , v ( s ) ) | K , s [ a , b ] .
With the use of Lemma 5, we have
| T 1 ( u , v ) ( t ) | a b G 11 ( t , s ) | f ( s , u ( s ) , v ( s ) ) | d s + a b G 12 ( t , s ) | g ( s , u ( s ) , v ( s ) ) | d s K μ 11 ( t a ) β 1 1 a b ( b s ) β 1 1 d s + K μ 12 ( t a ) β 1 1 a a ( b s ) β 2 1 d s = K μ 11 ( b a ) β 1 β 1 + K μ 12 ( b a ) β 2 β 2 ( t a ) β 1 1 ,
and
| T 2 ( u , v ) ( t ) | K μ 21 ( b a ) β 1 β 1 + K μ 22 ( b a ) β 2 β 2 ( t a ) β 2 1 .
This implies that T i maps all of E × E into the vector subspace E i of E, where E i is given by
E i = { u E : there is M > 0 such that | u ( t ) | M ( t a ) β i 1 , t [ a , b ] } .
Evidently, E i ( i = 1 , 2 ) are Banach spaces with the norm
u i = inf { M > 0 : | u ( t ) | M ( t a ) β i 1 , t [ a , b ] } .
Therefore, ( u * , v * ) is a nontrivial fixed point of T in E 1 × E 2 .
Using (H2) and Lemma 5, for all t [ 0 , 1 ] , we obtain
| u ( t ) | a b G 11 ( t , s ) | f ( s , u ( s ) , v ( s ) ) | d s + a b G 12 ( t , s ) | g ( s , u ( s ) , v ( s ) ) | d s a b G 11 ( t , s ) ( p 11 ( s ) | u ( s ) | + p 12 ( s ) | v ( s ) | ) d s a b G 12 ( t , s ) ( p 21 ( s ) | u ( s ) | + p 22 ( s ) | v ( s ) | ) d s μ 11 ( t a ) β 1 1 a b ( b s ) β 1 1 ( p 11 ( s ) u 1 ( s a ) β 1 1 + p 12 ( s ) v 2 ( s a ) β 2 1 ) d s + μ 12 ( t a ) β 1 1 a b ( b s ) β 2 1 ( p 21 ( s ) u 1 ( s a ) β 1 1 + p 22 ( s ) v 2 ( s a ) β 2 1 ) d s μ 11 a b p 11 ( s ) ( b s ) β 1 1 ( s a ) β 1 1 d s + μ 12 a b p 21 ( s ) ( b s ) β 2 1 ( s a ) β 1 1 d s u 1 ( t a ) β 1 1 + μ 11 a b p 12 ( s ) ( b s ) β 1 1 ( s a ) β 2 1 d s + μ 12 a b p 22 ( s ) ( b s ) β 2 1 ( s a ) β 2 1 d s v 2 ( t a ) β 1 1 .
Therefore, we obtain
u 1 I 11 ( p 11 , p 21 ) u 1 + I 12 ( p 12 , p 22 ) v 2 .
Similarly, using Lemma 5 and (H3), we conclude
v 1 I 21 ( p 11 , p 21 ) u 1 + I 22 ( p 12 , p 22 ) v 2 .
Combining (16) with (18), we deduce that
u 1 v 2 I 11 ( p 11 , p 21 ) I 12 ( p 12 , p 22 ) I 21 ( p 11 , p 21 ) I 22 ( p 12 , p 22 ) u 1 v 2 .
With the consideration of Lemma 6, Lemma 7 and (15), we deduce that
u * 1 = v * 2 = 0 ,
which contradicts the nontriviality of ( u * , v * ) . This proves (14). □
Let
g i j ( s ) = ( b s ) β i 1 ( s a ) β j 1 , s [ a , b ] , i , j = 1 , 2 .
Now, we differentiate g i j ( s ) on ( a , b ) and we obtain
g i j ( s ) = ( b s ) β i 2 ( s a ) β j 2 [ a ( β i 1 ) + b ( β j 1 ) s ( β i + β j 2 ) ]
which implies that g 3 ( s ) = 0 only at s i j * = b ( β j 1 ) + a ( β i 1 ) β i + β j 2 . Note that g i j ( a ) = g i j ( b ) = 0 . By the continuity of g i j , we conclude that
g i j ( s ) g i j ( s i j * ) = ( b a ) β i + β j 2 ( β i 1 ) β i 1 ( β j 1 ) β j 1 ( β i + β j 2 ) β i + β j 2 .
For p i j C [ a , b ] ( i , j = 1 , 2 ) , let
ζ i j ( p 1 j , p 2 j ) = μ i 1 g 1 j ( s 1 j * ) a b p 1 j ( s ) d s + μ i 2 g 2 j ( s 2 j * ) a b p 2 j ( s ) d s , i , j = 1 , 2 .
The following theorems are immediate.
Theorem 3. 
Suppose that (H0)-(H3) are satisfied. If (1) have a nontrivial solution, then
ζ 11 ( p 11 , p 21 ) + ζ 22 ( p 12 , p 22 ) + [ ζ 11 ( p 11 , p 21 ) ζ 22 ( p 12 , p 22 ) ] 2 + 4 ζ 12 ( p 12 , p 22 ) ζ 21 ( p 11 , p 21 ) 2 .
Finally, we will compare two inequalities (10) and (19) in the case that β 1 = β 2 . If β 1 = β 2 , by (18), we obtain
g ( s * ) g i j ( s i j * ) = b a 2 2 β 1 2 , i , j = 1 , 2 ,
and square nonnegative matrices M 1 = ( I i j ( p 1 j , p 2 j ) ) 2 × 2 and M 2 = ( ζ i j ( p 1 j , p 2 j ) ) 2 × 2 become
M 1 = λ 11 λ 12 λ 21 λ 22 a b p 11 ( s ) d s a b p 12 ( s ) d s a b p 21 ( s ) d s a b p 22 ( s ) d s ,
and
M 2 = g ( s * ) μ 11 μ 12 μ 21 μ 22 a b p 11 ( s ) d s a b p 12 ( s ) d s a b p 21 ( s ) d s a b p 22 ( s ) d s .
It is worth noting that
λ 11 = 1 Γ ( β 1 ) b a 4 β 1 1 + 1 κ Γ ( β 1 ) i = 1 n ( κ 22 a 1 i + κ 12 a 3 i ) ( ξ i a ) β 1 1 ( b ξ i ) β 1 1 ( b a ) β 1 1 ,
and
μ 11 g ( s * ) = 1 Γ ( β 1 ) b a 4 β 1 1 + 1 κ Γ ( β 1 ) i = 1 n ( κ 22 a 1 i + κ 12 a 3 i ) ( ξ i a ) β 1 1 1 4 β 1 1 .
Then, it is easily seen that λ 11 μ 11 g ( s * ) if a + 3 b 4 ξ i < b for all i = 1 , 2 , , n , and λ 11 μ 11 g ( s * ) if a < ξ i a + 3 b 4 for all i = 1 , 2 , , n . In the same way, for i , j = 1 , 2 , we can prove
(i)
λ i j μ i j g ( s * ) if a + 3 b 4 ξ k , η k < b for all k = 1 , 2 , , n ,
(ii)
λ j 1 μ j 1 g ( s * ) if a < ξ k , η k a + 3 b 4 for all k = 1 , 2 , , n .
Therefore, if a + 3 b 4 ξ k , η k < b for all k = 1 , 2 , , n , we obtain that 0 M 1 M 2 . This together with Lemma 8 show that inequality (10) is the improvement of inequality (19). Similarly, inequality (19) is the improvement of inequality (10) if a < ξ k , η k a + 3 b 4 for all k = 1 , 2 , , n .

4. Conclusions

In this article, we investigate a system of Riemann–Liouville fractional differential equations with multi-point coupled boundary conditions. The first Lyapunov-type inequality is obtained by matrix analysis and the Green’s function approach, where the Green’s function approach implies deriving the Green’s function of the equivalent integral form of the boundary value problem being considered and then finding an upper bound of its Green’s function (note that the maximum value of the Green’s function of (1) exists, but the analysis is somewhat complicated). The second Lyapunov-type inequality is obtained by matrix analysis and the properties of the operator T defined on Banach spaces E 1 × E 2 . The third Lyapunov-type inequality is the corollary of the second Lyapunov-type inequality. Finally, the comparison between the first and the third Lyapunov-type inequalities is given under certain special conditions. We expect that the second approach used in this paper can be applied to study a system of various fractional boundary value problems, and we will continue to discuss the optimal version of the constants appearing in the Lyapunov-type inequality and seek other ways to obtain fractional Lyapunov-type inequality in future papers.

Author Contributions

Investigation and formal analysis, Y.Z. and Y.C.; writing—original draft, Y.C.; writing—review and editing, Y.Z. and Y.C. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by the National Natural Science Foundation of China (11571207), the Shandong Natural Science Foundation (ZR2018MA011).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lyapunov, A.M. Probleme général de la stabilité du mouvement, (French Translation of a Russian paper dated 1893). Ann. Fac. Sci. Univ. Toulouse 1907, 2, 27–247. [Google Scholar]
  2. Brown, R.C.; Hinton, D.B. Lyapunov inequalities and their applications, in Survey on Classical Inequalities. In Mathematics and Its Applications; Springer: Dordrecht, The Netherlands, 2000; Volume 517. [Google Scholar]
  3. Pinasco, J.P. Lyapunov-Type Inequalities with Applications to Eigenvalue Problems; Springer: New York, NY, USA, 2013. [Google Scholar]
  4. Cheng, S.S. Lyapunov inequalities for differential and difference equations. Fasc. Math. 1991, 23, 25–41. [Google Scholar]
  5. Tiryaki, A. Recent developments of Lyapunov-type inequalities. Adv. Dyn. Syst. Appl. 2010, 5, 231–248. [Google Scholar]
  6. Ferreira Rui, A.C. A Lyapunov-type inequality for a fractional boundary value problem. Fract. Calc. Appl. Anal. 2013, 16, 978–984. [Google Scholar] [CrossRef]
  7. Ferreira Rui, A.C. On a Lyapunov-type inequality and the zeros of a certain MittagLeffler function. J. Math. Anal. Appl. 2014, 412, 1058–1063. [Google Scholar] [CrossRef]
  8. Ferreira Rui, A.C. Lyapunov-type inequality for an anti-periodic fractional boundary value problem. Fract. Calc. Appl. Anal. 2017, 20, 284–291. [Google Scholar] [CrossRef]
  9. Ntouyas, S.K.; Ahmad, B.; Horikis, T.P. Recent developments of Lyapunov-type inequalities for fractional differential equations. In Differential and Integral Inequalities; Andrica, D., Rassias, T.M., Eds.; Springer Optimization and Its Applications; Springer: Berlin/Heidelberg, Germany, 2019; pp. 619–686. [Google Scholar]
  10. Ntouyas, S.K.; Ahmad, B. Lyapunov-type inequalities for fractional differential equations: A survey. Surv. Math. Appl. 2021, 16, 43–93. [Google Scholar]
  11. Ntouyas, S.K.; Ahmad, B.; Tariboon, J. A Survey on Recent Results on Lyapunov-Type Inequalities for Fractional Differential Equations. Fractal Fract. 2022, 6, 273. [Google Scholar] [CrossRef]
  12. Jleli, M.; Nieto, J.J.; Samet, B. Lyapunov-type inequalities for a higher order fractional differential equation with fractional integral boundary conditions. Electron. J. Qual. Theory Differ. Equations 2017, 16, 1–17. [Google Scholar] [CrossRef]
  13. Bachar, I.; Eltayeb, H. Hartman-type and Lyapunov-type inequalities for a fractional differential equation with fractional boundary conditions. Discret. Dyn. Nat. Soc. 2020, 2020, 8234892. [Google Scholar] [CrossRef]
  14. Basua, D.; Jonnalagadda, J.M.; Satpathi, D.K. Lyapunov-type inequalities for Riemann–Liouville type fractional boundary value problems with fractional boundary conditions. Adv. Theory Nonlinear Anal. Its Appl. 2019, 3, 53–63. [Google Scholar]
  15. Wang, Y.; Wang, Q. Lyapunov-type inequalities for fractional differential equations under multi-point boundary conditions. J. Comput. Anal. Appl. 2019, 26, 707–716. [Google Scholar] [CrossRef]
  16. Jleli, M.; O’Regan, D.; Samet, B. Lyapunov-type inequalities for coupled systems of nonlinear fractional differential equations via a fixed point approach. J. Fixed Point Theory Appl. 2019, 21, 1–15. [Google Scholar] [CrossRef]
  17. Amann, H. Parabolic evolution equations and nonlinear boundary conditions. J. Differ. Equations 1988, 72, 201–269. [Google Scholar] [CrossRef]
  18. Deng, K. Global existence and blow-up for a system of heat equations with nonlinear boundary condition. Math. Methods Appl. Sci. 1995, 18, 307–315. [Google Scholar] [CrossRef]
  19. Aronson, D.G. A comparison method for stability analysis of nonlinear parabolic problems. SIAM Rev. 1978, 20, 245–264. [Google Scholar] [CrossRef]
  20. Asif, N.A.; Khan, R.A. Positive solutions to singular system with four-point coupled boundary conditions. J. Math. Anal. Appl. 2012, 386, 848–861. [Google Scholar] [CrossRef]
  21. Agarwal, R.P.; Ahmad, B.; Garout, D.; Alsaedi, A. Existence results for coupled nonlinear fractional differential equations equipped with nonlocal coupled flux and multi-point boundary conditions. Chaos Solitons Fractals 2017, 102, 149–161. [Google Scholar] [CrossRef]
  22. Ahmad, B.; Alghamdi, N.; Alsaedi, A.; Ntouyas, S.K. A system of coupled multi-term fractional differential equations with three-point coupled boundary conditions. Fract. Calc. Appl. Anal. 2019, 22, 601–618. [Google Scholar] [CrossRef]
  23. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematics Studies; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  24. Diethelm, K. The Analysis of Fractional Differential Equations; Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 2010; Volume 2004. [Google Scholar]
  25. Perov, A.I. On the Cauchy problem for a system of ordinary differential equations. Pviblizhen. Met. Reshen. Differ. Uvavn. 1964, 2, 115–134. [Google Scholar]
  26. Bapat, R.B.; Raghavan, T. Nonnegative Matrices and Applications; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zou, Y.; Cui, Y. Lyapunov-Type Inequalities for Systems of Riemann-Liouville Fractional Differential Equations with Multi-Point Coupled Boundary Conditions. Fractal Fract. 2023, 7, 454. https://doi.org/10.3390/fractalfract7060454

AMA Style

Zou Y, Cui Y. Lyapunov-Type Inequalities for Systems of Riemann-Liouville Fractional Differential Equations with Multi-Point Coupled Boundary Conditions. Fractal and Fractional. 2023; 7(6):454. https://doi.org/10.3390/fractalfract7060454

Chicago/Turabian Style

Zou, Yumei, and Yujun Cui. 2023. "Lyapunov-Type Inequalities for Systems of Riemann-Liouville Fractional Differential Equations with Multi-Point Coupled Boundary Conditions" Fractal and Fractional 7, no. 6: 454. https://doi.org/10.3390/fractalfract7060454

APA Style

Zou, Y., & Cui, Y. (2023). Lyapunov-Type Inequalities for Systems of Riemann-Liouville Fractional Differential Equations with Multi-Point Coupled Boundary Conditions. Fractal and Fractional, 7(6), 454. https://doi.org/10.3390/fractalfract7060454

Article Metrics

Back to TopTop