Next Article in Journal
Limiting Behaviors of Stochastic Spread Models Using Branching Processes
Next Article in Special Issue
Synchronization of Fractional Partial Difference Equations via Linear Methods
Previous Article in Journal
Higher-Order Blended Compact Difference Scheme on Nonuniform Grids for the 3D Steady Convection-Diffusion Equation
Previous Article in Special Issue
The Dynamics of a General Model of the Nonlinear Difference Equation and Its Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solvability of a Boundary Value Problem Involving Fractional Difference Equations

1
School of Sciences and Arts, Suqian University, Suqian 223800, China
2
School of Mathematical Sciences, Chongqing Normal University, Chongqing 401331, China
3
School of Mathematical and Statistical Sciences, University of Galway, H91 TK33 Galway, Ireland
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(7), 650; https://doi.org/10.3390/axioms12070650
Submission received: 8 June 2023 / Revised: 25 June 2023 / Accepted: 28 June 2023 / Published: 29 June 2023
(This article belongs to the Special Issue Advances in Difference Equations)

Abstract

:
In this current work, we apply the topological degree and fixed point theorems to investigate the existence, uniqueness, and multiplicity of solutions for a boundary value problem associated with a fractional-order difference equation. Moreover, we provide some appropriate examples to verify our main conclusions.
MSC:
34B10; 34B18; 34A34; 45G15; 45M20

1. Introduction

Let [ κ , δ ] N κ : = { κ , κ + 1 , κ + 2 , , δ } ( δ κ   N 1 , where N κ : = { κ , κ + 1 , κ + 2 , } . In the current work, we shall discuss the solvability of the fractional difference boundary value problem
Δ v 3 v ψ ( t ) = g ( t + v 1 , ψ ( t + v 1 ) ) , t [ 0 , b + 2 ] N 0 , ψ ( v 3 ) = Δ v 3 α ψ ( t ) t = v α 2 = Δ v 3 β ψ ( t ) t = v + b + 2 β = 0 ,
where v ( 2 , 3 ] , β ( 1 , 2 ) , v β ( 1 , + ) , α ( 0 , 1 ) , b ( 3 , + ) ( b N ) , and Δ v 3 v is a discrete fractional-order operator defined by
Δ a ν ψ ( t ) : = 1 Γ ( ν ) s = a t + ν ( t s 1 ) ν 1 ̲ ψ ( s ) , N 1 < ν < N , Δ N ψ ( t ) , ν = N ,
where N N with 0 N 1 < v N . As in [1], this definition is equivalent to (2) in Section 2.
The theory of fractional calculus has been widely used in modern mathematics for more than 300 years, and the study of solutions of fractional differential (difference) equations arises in real-world problems in the field of physics, mechanics, chemistry, and engineering. For example, in [2], the authors extended the variational approach to the fractional discrete case and introduced the Gompertz fractional difference equation
Δ 0 α ln G ( t α + 1 ) = ( b 1 ) ln G ( t ) + a ,
which can be used to describe tumor growth, a special relationship between tumor size and time, and is of special interest since growth estimation is very critical in clinical practice. Here, a , b are parameters and α ( 0 , 1 ] . One can also find some other applications for the Gompertz fractional difference equation in [1]. In [3], the authors introduced the following discrete logistic map and investigated the chaotic behavior:
C Δ a v ψ ( t ) = μ ψ ( t + v 1 ) ( 1 ψ ( t + v 1 ) ) , t N a + 1 v , 0 < v 1 , ψ ( a ) = c ,
where C Δ a v is the left Caputo-like delta difference defined by
C Δ a v ψ ( t ) = 1 Γ ( m v ) s = a t ( m v ) ( t s 1 ) m v 1 ̲ Δ s m ψ ( s ) ,
where t N a + m v , m = [ v ] + 1 .
We note that in [4], the author mentioned that discretization is inevitable for fractional differential equations. To date, they are only used as the starting point for approximate solution calculations, and there is no special research on fractional difference equations. Therefore, from the perspective of theory and application, this is a big gap. Many developments in the theory are now taking place, and two books [5,6] are sources for mathematicians who are interested in this area. However, we still note that most works focus on fractional-order differential equations, while the research on fractional-order difference equations is quite small (we refer the reader to [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26]). In [7], the authors investigated positive solutions for the discrete fractional boundary value problems
Δ v 2 v χ ( t ) = F ( t + v 1 , χ ( t + v 1 ) ) , 1 < v 2 , χ ( v 2 ) = Δ v 1 v 1 χ ( v + N ) = 0 ,
where t [ 0 , N + 1 ] N 0 and F : [ v 1 , v + N ] N v 1 × R R + satisfies some superlinear or sublinear conditions. In [8], the authors utilized fixed-point methods to investigate the solvability of a fractional difference equation with a p-Laplacian operator
Δ β [ ϕ p ( Δ α χ ) ] ( t ) + F ( α + β + t 1 , χ ( α + β + t 1 ) ) = 0 , t [ 0 , b ] N 0 , Δ α χ ( β 2 ) = Δ α χ ( β + b ) = 0 , χ ( α + β 4 ) = χ ( α + β + b ) = 0 ,
where F : [ α + β 4 , α + β + b ] N α + β 4 × R R satisfies a Lipschitz condition, and ϕ p ( z ) = | z | p 2 z , p > 1 , z R . In [9], the authors utilized the fixed point index to consider the solvability of the system of fractional-order difference boundary value problems
Δ v χ 1 ( t ) = F 1 ( t + v 1 , χ 1 ( t + v 1 ) , χ 2 ( t + v 1 ) ) , t [ 0 , T ] Z , Δ v χ 2 ( t ) = F 2 ( t + v 1 , χ 1 ( t + v 1 ) , χ 2 ( t + v 1 ) ) , t [ 0 , T ] Z , χ 1 ( v 1 ) = χ 1 ( v + T ) , χ 2 ( v 1 ) = χ 2 ( v + T ) ,
where F i ( i = 1 , 2 ) are semipositone nonlinearities.
We note that usually one expresses the solutions of fractional-order equations by a Green’s function. However, not all fractional-order difference equations can be obtained in this way, for example, in [10], the authors studied the problem
Δ α χ ( t ) = F ( t + α 1 , χ ( t + α 1 ) ) , t [ 0 , T ] N 0 , α ( 1 , 2 ] , χ ( α 2 ) = 0 , χ ( α + T ) = Δ β χ ( ζ + β ) , ζ N α 2 , α + T 1 , β > 0 ,
and showed it is equivalent to
χ ( t ) = t α 1 ̲ Θ Γ ( α ) 1 Γ ( β ) s = α ζ ξ = 0 s α ( ζ + β σ ( s ) ) β 1 ( s σ ( ξ ) ) α 1 F ( ξ + α 1 , χ ( ξ + α 1 ) ) s = 0 T ( T + α σ ( s ) ) α 1 F ( s + α 1 , χ ( s + α 1 ) ) + 1 Γ ( α ) s = 0 t α ( t σ ( s ) ) α 1 F ( s + α 1 , χ ( s + α 1 ) ) ,
where
Θ = s = 1 ζ α + 2 ( ζ + β s α + 1 ) β 1 Γ ( s + α 1 ) Γ ( β ) Γ ( s ) Γ ( α + T + 1 ) Γ ( T + 2 ) .
Clearly, the integral form is very complicated and cannot be formulated via some suitable Green’s function.
Inspired by the aforementioned works, in this paper, via a Green’s function, we use the topological degree and fixed point theorems to consider the existence, uniqueness, and multiplicity of solutions to (1). Furthermore, we present some examples to illustrate our main results.

2. Preliminary

In this section, we first offer some basic materials for discrete fractional calculus; see [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26] and the references therein.
Definition 1. 
Let
t v ̲ : = Γ ( t + 1 ) Γ ( t + 1 v ) , t , v R . If t + 1 v is a pole of Γ ( · ) and t + 1 is not a pole, then t v ̲ = 0 .
Definition 2. 
For v > 0 , a function F ’s v-th fractional sum is defined by
Δ a v F ( t ) = 1 Γ ( v ) s = a t v ( t s 1 ) v 1 ̲ F ( s ) , t N a + v .
F ’s v-th fractional difference is defined by
Δ a v F ( t ) = Δ N Δ a v N F ( t ) , t N a + N v ,
where N N with 0 N 1 < v N .
Let χ : [ v 1 , b + v + 1 ] N v 1 R be a given function. Then, we consider the problem
Δ v 3 v ψ ( t ) = χ ( t + v 1 ) , t [ 0 , b + 2 ] N 0 , ψ ( v 3 ) = Δ v 3 α ψ ( t ) t = v α 2 = Δ v 3 β ψ ( t ) t = v + b + 2 β = 0 ,
where v , α , β , b can be founded in (1).
Lemma 1 
(see [11]). Problem (3) has a unique solution
ψ ( t ) = s = 0 b + 2 G ( t , s ) χ ( s + v 1 ) , t [ v 1 , b + v + 1 ] N v 1 ,
where G is the Green’s function given by
G ( t , s ) = 1 Γ ( v ) t v 1 ̲ ( v + b β s + 1 ) v β 1 ̲ ( v + b β + 2 ) v β 1 ̲ ( t s 1 ) v 1 ̲ , 0 s < t v + 1 b + 2 , t v 1 ̲ ( v + b β s + 1 ) v β 1 ̲ ( v + b β + 2 ) v β 1 ̲ , 0 t v + 1 s b + 2 .
Lemma 2 
(see [11]). The Green’s function (4) has the properties
(G1) G ( t , s ) > 0 , ( t , s ) [ v 1 , b + v + 1 ] N v 1 × [ 0 , b + 2 ] N 0 ,
(G2) t v 1 ̲ G ( b + v + 1 , s ) ( b + v + 1 ) v 1 ̲ G ( t , s ) t v 1 ̲ ( v + b β s + 1 ) v β 1 ̲ Γ ( v ) ( v + b β + 2 ) v β 1 ̲ , ( t , s ) [ v 1 , b + v + 1 ] N v 1 × [ 0 , b + 2 ] N 0 ,
(G3) q ( t ) G ( b + v + 1 , s ) G ( t , s ) G ( b + v + 1 , s ) , q ( t ) = t v 1 ̲ ( b + v + 1 ) v 1 ̲ , ( t , s ) [ v 1 , b + v + 1 ] N v 1 × [ 0 , b + 2 ] N 0 .
Lemma 3 
(see [11]). Let φ ( s + v 1 ) = G ( b + v + 1 , s ) , s [ 0 , b + 2 ] N 0 . Then, the following inequalities hold:
t = v 1 b + v + 1 G ( t , s ) φ ( t ) κ 2 φ ( s + v 1 ) , κ 2 = t = v 1 b + v + 1 φ ( t ) , s [ 0 , b + 2 ] N 0 ,
and
s = 0 b + 2 G ( t , s ) q ( s + v 1 ) κ 1 q ( t ) , κ 1 = s = 0 b + 2 G ( b + v + 1 , s ) q ( s + v 1 ) , t [ v 1 , b + v + 1 ] N v 1 .
Let E be a set of all maps from [ v 3 , b + v + 1 ] N v 3 to R , and ψ = max t [ v 3 , b + v + 1 ] N v 3 | ψ ( t ) | . Then, E is a Banach space. Moreover, define a set P = ψ E : ψ ( t ) 0 , t [ v 1 , b + v + 1 ] N v 1 . Then, P is a cone on E. Lemma 3 enables us to obtain that (1) is equivalent to the sum equation
ψ ( t ) = s = 0 b + 2 G ( t , s ) g ( s + v 1 , ψ ( s + v 1 ) ) : = ( B ψ ) ( t ) , t [ v 1 , b + v + 1 ] N v 1 ,
where G is defined in Lemma 3. Obviously, ψ E { 0 } is a solution for (1) when ψ E { 0 } is a fixed point of B .
Lemma 4. 
Let P 0 = ψ P : ψ ( t ) q ( t ) ψ , t [ v 1 , b + v + 1 ] N v 1 . Then, L ( P ) P 0 , where
( L ψ ) ( t ) = s = 0 b + 2 G ( t , s ) ψ ( s + v 1 ) , t [ v 1 , b + v + 1 ] N v 1 .
Lemma 5 
(see [27] Theorem A.3.3). Let E be a Banach space, Ω E a bounded open set, and T : Ω E be a continuous compact operator. If there is an x 0 E { 0 } such that
x T x μ x 0 , x Ω , μ 0 ,
then deg ( I T , Ω , 0 ) = 0 , where deg denotes the topological degree.
Lemma 6 
(see [27] Lemma 2.5.1). Let E be a Banach space, Ω E a bounded open set with 0 Ω , and T : Ω E be a continuous compact operator. If
T x μ x , x Ω , μ 1 ,
then deg ( I T , Ω , 0 ) = 1 .
Lemma 7 
(see [28,29]). Let X be a Banach space and P be a cone on X. Define functionals as follows: α , γ : P R + are continuous increasing and β : P R + is continuous. Moreover, there exists M > 0 , 0 < a ˜ < c ˜ such that
α ( 0 ) < a ˜ , γ ( x ) β ( x ) α ( x ) and x M γ ( x ) , x P ( γ , c ˜ ) ¯ : = { x P : γ ( x ) < c ˜ } ¯ .
Furthermore, there is a completely continuous operator T : P ( γ , c ˜ ) ¯ P and a constant b ˜ > 0 with 0 < a ˜ < b ˜ < c ˜ such that β ( λ x ) λ β ( x ) for λ ( 0 , 1 ] , x P ( β , b ˜ ) , and
(E1) γ ( T x ) < c ˜ , x P ( γ , c ˜ ) ;
(E2) β ( T x ) > b ˜ , x P ( β , b ˜ ) ;
(E3) α ( T x ) < a ˜ , x P ( α , a ˜ ) .
Then, T has at least three fixed points x 1 , x 2 , x 3 P ( γ , c ˜ ) ¯ such that
0 α x 1 < a ˜ < α x 2 , β x 2 < b ˜ < β x 3 , γ x 3 < c ˜ .
In the following, we present some lemmas involving the theory of mixed monotone operators. Let ( E , · ) be a real Banach space which is partially ordered by a cone P E , i.e., x y   y x P . If x y and x y , then we mean that x < y or y > x . Moreover, for a fixed h > 0 , we define P h = { x E x h } , in which ∼ is an equivalence relation, i.e., x y implies that there are λ , μ > 0 such that λ x y μ x , x , y E .
Definition 3 
(see [30,31]). If u i , v i ( i = 1 , 2 ) P , u 1 < u 2 , v 1 > v 2 imply A u 1 , v 1 A u 2 , v 2 , then A : P × P P is called a mixed monotone operator.
Definition 4 
(see [30,31]). If A ( t x ) t A x , t ( 0 , 1 ) , x P , then A : P P is said to be sub-homogeneous.
Lemma 8 
(see [30,31]). Let B : P P be an increasing sub-homogeneous operator, A : P × P P a mixed monotone operator and satisfy
A t x , t 1 y t α A ( x , y ) , t , α ( 0 , 1 ) , x , y P .
If
(C1) There is a h 0 P h such that A h 0 , h 0 P h and B h 0 P h ;
(C2) There is a constant δ 0 > 0 such that A ( x , y ) δ 0 B x , x , y P .
Then,
(D1) A : P h × P h P h , B : P h P h ;
(D2) There are u 0 , v 0 P h and r ( 0 , 1 ) such that r v 0 u 0 < v 0 , u 0 A u 0 , v 0 + B u 0 A v 0 , u 0 + B v 0 v 0 ;
(D3) A ( x , x ) + B x = x has a unique solution x * in P h ;
(D4) For any initial values x 0 , y 0 P h , the sequences x n = A x n 1 , y n 1 + B x n 1 , y n = A y n 1 , x n 1 + B y n 1 converge to x * as n .

3. Main Results

In the section, we will state our main theorems and give their proof. In the first theorem, we obtain an existence result on nontrivial solutions for (1) when the nonlinearity can change sign.
Theorem 1. 
Suppose that the following assumptions hold:
(H1) g ( t , ψ ) : [ v 1 , b + v + 1 ] N v 1 × R R is a continuous function;
(H2) There are nonnegative continuous functions γ 1 ( t ) , γ 2 ( t ) and M ( ψ ) with γ 2 ( t ) 0 , t [ v 1 , b + v + 1 ] N v 1 such that
g ( t , ψ ) γ 1 ( t ) γ 2 ( t ) M ( ψ ) , ( t , ψ ) [ v 1 , b + v + 1 ] N v 1 × R ;
(H3) lim | ψ | + M ( ψ ) | ψ | = 0 ;
(H4) lim inf | ψ | + g ( t , ψ ) | ψ | > κ 1 1 , uniformly in t [ v 1 , b + v + 1 ] N v 1 ;
(H5) lim inf | ψ | 0 + | g ( t , ψ ) | | ψ | < κ 2 1 , uniformly in t [ v 1 , b + v + 1 ] N v 1 .
Then, (1) has one nontrivial solution.
Proof. 
From (H3), for any given ε > 0 , there exists Y 0 > 0 such that M ( ψ ) ε | ψ | for | ψ | > Y 0 . Let M * = max | ψ | [ 0 , Y 0 ] M ( ψ ) . Then, we have
M ( ψ ) ε | ψ | + M * , ψ R .
By (H4), there exist δ 1 > 0 and Y 1 Y 0 such that g ( t , ψ ) ( κ 1 1 + δ 1 ) | ψ | for | ψ | > Y 1 and t [ v 1 , b + v + 1 ] N v 1 . Furthermore, let C ˜ g = max ( t , ψ ) [ v 1 , b + v + 1 ] N v 1 × [ 0 , Y 1 ] | g ( t , ψ ) | . Then, we obtain
g ( t , ψ ) ( κ 1 1 + δ 1 ) | ψ | C ˜ g , t [ v 1 , b + v + 1 ] N v 1 , ψ R .
Note that δ 1 can be greater than ε γ 2 ; then, from (H2) and (H3) and (8), we have
g ( t , ψ ) ( κ 1 1 + δ 1 ε γ 2 ) | ψ | γ 1 ( t ) C g , t [ v 1 , b + v + 1 ] N v 1 , ψ R ,
where C g = C ˜ g + γ 2 M * . Let
R > max κ 2 γ 1 + γ 2 M * + C g 1 ε κ 2 γ 2 , ( γ 1 + γ 2 M * + C g ) κ 2 ( δ 1 ε γ 2 ) ( b + v + 1 ) v 1 ̲ + ( κ 1 1 + δ 1 ε γ 2 ) N 1 δ 1 ε γ 2 ( b + v + 1 ) v 1 ̲ ( 1 ε κ 2 γ 2 ) ε γ 2 N 1 ( κ 1 1 + δ 1 ε γ 2 ) ,
where
N 1 = τ = 0 b + 2 ( v + b β τ + 1 ) v β 1 ̲ ( v + b β + 2 ) v β 1 ̲ Γ ( v ) .
We prove that
ψ B ψ μ q , ψ B R , μ 0 ,
where q is given in Lemma 4, and
B R = { ψ E : ψ < R } , B R = { ψ E : ψ = R } .
Proof by contradiction. Then, there are ψ B R , μ 0 such that
ψ B ψ = μ q .
Note that if μ = 0 and ψ B R is a nontrivial solution to (1), the theorem has been obtained. So, we only consider the case μ > 0 . Moreover, we also find that
q P 0 .
In order to prove our theorem, we need to define a function ψ ˜ as follows:
ψ ˜ ( t ) = s = 0 b + 2 G ( t , s ) γ 1 ( s + v 1 ) + γ 2 ( s + v 1 ) M ( ψ ( s + v 1 ) ) + C g , t [ v 1 , b + v + 1 ] N v 1 , ψ B R .
Then, we get the following claims:
Claim i. Note that γ 1 + γ 2 M ( ψ ) + C g P , and Lemma 6 implies that
ψ ˜ P 0 .
Claim ii. From (12), we find
ψ ( t ) + ψ ˜ ( t ) = ( B ψ ) ( t ) + ψ ˜ ( t ) + μ q ( t ) = s = 0 b + 2 G ( t , s ) g ( s + v 1 , ψ ( s + v 1 ) ) + γ 1 ( s + v 1 ) + γ 2 ( s + v 1 ) M ( ψ ( s + v 1 ) ) + C g + μ q ( t ) ,
for all t [ v 1 , b + v + 1 ] N v 1 . Note that g + γ 1 + γ 2 M + C g P and q P 0 , and we have
ψ + ψ ˜ P 0 .
Claim iii. From (8) and (10), we have
ψ ˜ s = 0 b + 2 G ( b + v + 1 , s ) γ 1 + γ 2 ( ε ψ + M * ) + C g = κ 2 γ 1 + γ 2 ( ε ψ + M * ) + C g < R .
From Claim ii and (9), we have
( B ψ ) ( t ) + ψ ˜ ( t ) = s = 0 b + 2 G ( t , s ) g ( s + v 1 , ψ ( s + v 1 ) ) + γ 1 ( s + v 1 ) + γ 2 ( s + v 1 ) M ( ψ ( s + v 1 ) ) + C g s = 0 b + 2 G ( t , s ) g ( s + v 1 , ψ ( s + v 1 ) ) + γ 1 ( s + v 1 ) + C g s = 0 b + 2 G ( t , s ) ( κ 1 1 + δ 1 ε γ 2 ) | ψ ( s + v 1 ) | γ 1 ( s + v 1 ) C g + γ 1 ( s + v 1 ) + C g ( κ 1 1 + δ 1 ε γ 2 ) s = 0 b + 2 G ( t , s ) ψ ( s + v 1 ) = ( κ 1 1 + δ 1 ε γ 2 ) s = 0 b + 2 G ( t , s ) ψ ( s + v 1 ) + ψ ˜ ( s + v 1 ) s = 0 b + 2 G ( t , s ) ψ ˜ ( s + v 1 ) κ 1 1 s = 0 b + 2 G ( t , s ) ψ ( s + v 1 ) + ψ ˜ ( s + v 1 ) .
The last inequality in (15) holds if
( δ 1 ε γ 2 ) s = 0 b + 2 G ( t , s ) ψ ( s + v 1 ) + ψ ˜ ( s + v 1 ) ( κ 1 1 + δ 1 ε γ 2 ) s = 0 b + 2 G ( t , s ) ψ ˜ ( s + v 1 ) 0 ,
for t [ v 1 , b + v + 1 ] N v 1 . In what follows, we prove (16). Indeed, from Claim ii we have ψ ( t ) + ψ ˜ ( t ) q ( t ) ψ + ψ ˜ q ( t ) ( ψ ψ ˜ ) , t [ v 1 , b + v + 1 ] N v 1 . Therefore, from (4) and (10), we obtain
( δ 1 ε γ 2 ) s = 0 b + 2 G ( t , s ) ψ ( s + v 1 ) + ψ ˜ ( s + v 1 ) ( κ 1 1 + δ 1 ε γ 2 ) s = 0 b + 2 G ( t , s ) ψ ˜ ( s + v 1 ) ( δ 1 ε γ 2 ) s = 0 b + 2 G ( t , s ) q ( s + v 1 ) ( ψ ψ ˜ )   ( κ 1 1 + δ 1 ε γ 2 ) s = 0 b + 2 G ( t , s ) τ = 0 b + 2 G ( s + v 1 , τ ) γ 1 ( τ + v 1 ) + γ 2 ( τ + v 1 ) M ( ψ ( τ + v 1 ) ) + C g ( δ 1 ε γ 2 ) s = 0 b + 2 G ( t , s ) ( s + v 1 ) v 1 ̲ ( b + v + 1 ) v 1 ̲ ( ψ ψ ˜ )   ( κ 1 1 + δ 1 ε γ 2 ) s = 0 b + 2 G ( t , s ) τ = 0 b + 2 ( s + v 1 ) v 1 ̲ ( v + b β τ + 1 ) v β 1 ̲ ( v + b β + 2 ) v β 1 ̲ Γ ( v ) γ 1 + γ 2 ( ε ψ + M * ) + C g s = 0 b + 2 G ( t , s ) ( s + v 1 ) v 1 ̲ [ δ 1 ε γ 2 ( b + v + 1 ) v 1 ̲ R κ 2 γ 1 + γ 2 ( ε R + M * ) + C g   ( κ 1 1 + δ 1 ε γ 2 ) γ 1 + γ 2 ( ε R + M * ) + C g τ = 0 b + 2 ( v + b β τ + 1 ) v β 1 ̲ ( v + b β + 2 ) v β 1 ̲ Γ ( v ) ] 0 .
This implies that (15) holds, as required. Consequently, we have
( B ψ ) ( t ) + ψ ˜ ( t ) κ 1 1 L ( ψ + ψ ˜ ) ( t ) , t [ v 1 , b + v + 1 ] N v 1 .
Using (12), we obtain
ψ + ψ ˜ = B ψ + ψ ˜ + μ q κ 1 1 L ( ψ + ψ ˜ ) + μ q μ q , ψ B R , μ > 0 .
Define
μ * = sup { μ > 0 : ψ + ψ ˜ μ q } .
Note that μ * μ and ψ + ψ ˜ μ * q , and from (6), we have
ψ + ψ ˜ κ 1 1 L ( μ * q ) + μ q = κ 1 1 μ * L q + μ q ( μ * + μ ) q ,
which contradicts the definition of μ * . Hence, (11) holds, and Lemma 7 enables us to find
deg I B , B R , 0 = 0 .
From (H5), there exist δ 2 ( 0 , κ 2 1 ) and r > 0 such that
| g ( t , ψ ) | ( κ 2 1 δ 2 ) | ψ | , | ψ | [ 0 , r ] , t [ v 1 , b + v + 1 ] N v 1 .
For this r, we prove that
B ψ μ ψ , ψ B r , μ 1 .
Proof by contradiction. Then, there are ψ B r , μ 1 such that
B ψ = μ ψ | ψ | = 1 μ | B ψ | | B ψ | .
This, together with (18), implies that
| ψ ( t ) | s = 0 b + 2 G ( t , s ) | g ( s + v 1 , ψ ( s + v 1 ) ) | ( κ 2 1 δ 2 ) s = 0 b + 2 G ( t , s ) | ψ ( s + v 1 ) | .
Multiplying by φ ( t ) on the both sides of (20) and summing over [ v 1 , b + v + 1 ] , then (5) implies that
t = v 1 b + v + 1 | ψ ( t ) | φ ( t ) ( κ 2 1 δ 2 ) t = v 1 b + v + 1 s = 0 b + 2 G ( t , s ) φ ( t ) | ψ ( s + v 1 ) | ( κ 2 1 δ 2 ) κ 2 s = 0 b + 2 | ψ ( s + v 1 ) | φ ( s + v 1 ) = ( κ 2 1 δ 2 ) κ 2 t = v 1 b + v + 1 | ψ ( t ) | φ ( t ) .
This implies that t = v 1 b + v + 1 | ψ ( t ) | φ ( t ) = 0 , and thus ψ ( t ) 0 , t [ v 1 , b + v + 1 ] N v 1 . Clearly, this is contradictory to ψ B r . Hence, Lemma 8 shows that
deg I B , B r , 0 = 1 .
Equations (17) and (21) enable us to obtain
deg I B , B R B ¯ r , 0 = deg I B , B R , 0 deg I B , B r , 0 = 1 .
This implies that B has a fixed point in B R B ¯ r , and (1) has a nontrivial solution. □
In the following theorem, using the generalized Avery–Henderson fixed point theorem (Lemma 9), we obtain triple positive solutions for (1) when the nonlinearity satisfies some bounded conditions.
Theorem 2. 
Suppose that there exist positive constants a ˜ , b ˜ , c ˜ with a ˜ < b ˜ < c ˜ , c ˜ b ˜ > s = 0 b + 2 G ( b + v + 1 , s ) t = r b + v + 1 G ( b + v + 1 , t v + 1 ) (r is a fixed point in ( v 1 , b + v + 1 ) N v 1 ) such that
(H6) g ( t , ψ ) : [ v 1 , b + v + 1 ] N v 1 × R + R + is a continuous function, and g ( t , 0 ) 0 , t [ v 1 , b + v + 1 ] N v 1 ;
(H7) g ( t , ψ ) < c ˜ q 0 s = 0 b + 2 G ( b + v + 1 , s ) , for t [ v 1 , b + v + 1 ] N v 1 , ψ [ 0 , c ˜ q 0 2 ] ;
(H8) g ( t , ψ ) > b ˜ q 0 t = r b + v + 1 G ( b + v + 1 , t v + 1 ) for t [ r , b + v + 1 ] N v 1 , ψ [ b ˜ , b ˜ q 0 2 ] ;
(H9) g ( t , ψ ) < a ˜ s = 0 b + 2 G ( b + v + 1 , s ) for t [ v 1 , b + v + 1 ] N v 1 , ψ [ 0 , a ˜ ] .
Then, (1) has at least three positive solutions y 1 , y 2 and y 3 satisfying
0 < α y 1 < a ˜ < α y 2 , β y 2 < b ˜ < β y 3 , γ y 3 < c ˜ .
Proof. 
Note that if q 0 = min t [ v 1 , b + v + 1 ] N v 1 q ( t ) > 0 , then from Lemma 6 and (H6) we have
B ( P ) P 0 .
Let α ( ψ ) = max t [ v 1 , b + v + 1 ] N v 1 ψ ( t ) , β ( ψ ) = min t [ r , b + v + 1 ] N v 1 ψ ( t ) and γ ( ψ ) = q 0 max t [ v 1 , r ] N v 1 ψ ( t ) . We easily know that α , γ : P R + are continuous, increasing functionals with α ( 0 ) = 0 , t [ v 1 , b + v + 1 ] N v 1 , ψ P and β ( λ ψ ) = λ β ( ψ ) . Moreover, for ψ P 0 , we have
γ ( ψ ) = q 0 max t [ v 1 , r ] N v 1 ψ ( t ) q 0 ψ min t [ v 1 , b + v + 1 ] N v 1 ψ ( t ) min t [ r , b + v + 1 ] N v 1 ψ ( t ) = β ( ψ ) α ( ψ ) ,
and
γ ( ψ ) q 0 min t [ v 1 , r ] N v 1 ψ ( t ) q 0 min t [ v 1 , b + v + 1 ] N v 1 ψ ( t ) q 0 2 ψ ,
i.e.,
ψ 1 q 0 2 γ ( ψ ) .
(i) For ψ P ( γ , c ˜ ) , we have
c ˜ = γ ( ψ ) q 0 2 ψ ,
which implies that
0 ψ ( t ) c ˜ q 0 2 , t [ v 1 , b + v + 1 ] N v 1 .
By (H7), we find
γ ( B ψ ) = q 0 max t [ v 1 , r ] N v 1 s = 0 b + 2 G ( t , s ) g ( s + v 1 , ψ ( s + v 1 ) ) q 0 max t [ v 1 , b + v + 1 ] N v 1 s = 0 b + 2 G ( t , s ) g ( s + v 1 , ψ ( s + v 1 ) ) < c ˜ q 0 s = 0 b + 2 G ( b + v + 1 , s ) q 0 s = 0 b + 2 G ( b + v + 1 , s ) = c ˜ .
(ii) For ψ P ( β , b ˜ ) , we have
b ˜ = β ( ψ ) = min t [ r , b + v + 1 ] N v 1 ψ ( t ) ψ 1 q 0 2 γ ( ψ ) 1 q 0 2 β ( ψ ) = b ˜ q 0 2 .
This implies that
b ˜ ψ ( t ) b ˜ q 0 2 , ψ P ( β , b ˜ ) , t [ r , b + v + 1 ] N v 1 .
This, combined with (H8), enables us to obtain
β ( B ψ ) = min t [ r , b + v + 1 ] N v 1 s = 0 b + 2 G ( t , s ) g ( s + v 1 , ψ ( s + v 1 ) ) s = 0 b + 2 min t [ r , b + v + 1 ] N v 1 q ( t ) G ( b + v + 1 , s ) g ( s + v 1 , ψ ( s + v 1 ) ) s = 0 b + 2 min t [ v 1 , b + v + 1 ] N v 1 q ( t ) G ( b + v + 1 , s ) g ( s + v 1 , ψ ( s + v 1 ) ) = q 0 s = 0 b + 2 G ( b + v + 1 , s ) g ( s + v 1 , ψ ( s + v 1 ) ) = q 0 t = v 1 b + v + 1 G ( b + v + 1 , t v + 1 ) g ( t , ψ ( t ) ) q 0 t = r b + v + 1 G ( b + v + 1 , t v + 1 ) g ( t , ψ ( t ) ) > b ˜ q 0 t = r b + v + 1 G ( b + v + 1 , t v + 1 ) q 0 t = r b + v + 1 G ( b + v + 1 , t v + 1 ) = b ˜ .
(iii) For ψ P ( α , a ˜ ) , we have
0 α ( ψ ) = max t [ v 1 , b + v + 1 ] N v 1 ψ ( t ) = a ˜ ,
and
0 ψ ( t ) a ˜ , ψ P ( α , a ˜ ) , t [ v 1 , b + v + 1 ] N v 1 .
This, together with (H9), implies that
α ( B ψ ) = max t [ v 1 , b + v + 1 ] N v 1 s = 0 b + 2 G ( t , s ) g ( s + v 1 , ψ ( s + v 1 ) ) < a ˜ s = 0 b + 2 G ( b + v + 1 , s ) s = 0 b + 2 G ( b + v + 1 , s ) = a ˜ .
Now, we have established that all the conditions in Lemma 9 hold, and note that B 0 0 , so we conclude that (1) has at least three positive solutions y i P { 0 } such that 0 < α y 1 < a ˜ < α y 2 , β y 2 < b ˜ < β y 3 , γ y 3 < c ˜ . □
In what follows, we study the problem
Δ v 3 v ψ ( t ) = f ( t + v 1 , ψ ( t + v 1 ) , ψ ( t + v 1 ) ) + g ( t + v 1 , ψ ( t + v 1 ) ) , t [ 0 , b + 2 ] N 0 , ψ ( v 3 ) = Δ v 3 α ψ ( t ) t = v α 2 = Δ v 3 β ψ ( t ) t = v + b + 2 β = 0 ,
where v , α , β , b are founded in (1). By Lemma 3, (22) is equivalent to the following equation
ψ ( t ) = s = 0 b + 2 G ( t , s ) [ f ( s + v 1 , ψ ( s + v 1 ) , ψ ( s + v 1 ) ) + g ( s + v 1 , ψ ( s + v 1 ) ) ] , t [ v 1 , b + v + 1 ] N v 1 ,
and let A : P × P P and B : P × P P be defined by
A ( y , x ) ( t ) = s = 0 b + 2 G ( t , s ) f ( s + v 1 , y ( s + v 1 ) , x ( s + v 1 ) ) , ( B y ) ( t ) = s = 0 b + 2 G ( t , s ) g ( s + v 1 , y ( s + v 1 ) ) .
Obviously, y * is a solution of (1) when y * = A ( y * , y * ) + B y * . In the following theorem, we study the operators A , B to help us to obtain the existence of solutions to (22). Moreover, the positive solution is unique, and it can be uniformly approximated by two appropriate iterative sequences.
Now, we list some assumptions for our nonlinearities f , g as follows:
(H10) f ( t , u , v ) : [ v 1 , b + v + 1 ] N v 1 × R + × R + R + , g ( t , u ) : [ v 1 , b + v + 1 ] N v 1 × R + R + are continuous functions;
(H11) f ( t , u , v ) is increasing about u R + for fixed t [ v 1 , b + v + 1 ] N v 1 and v R + and decreasing about v R + for fixed t [ v 1 , b + v + 1 ] N v 1 and u R + , and g ( t , u ) is increasing about u R + for fixed t [ v 1 , b + v + 1 ] N v 1 ;
(H12) For every t [ v 1 , b + v + 1 ] N v 1 , γ ( 0 , 1 ) , u , v R + , there is a constant ξ ( 0 , 1 ) such that f t , γ u , γ 1 v γ ξ f ( t , u , v ) and g ( t , γ u ) γ g ( t , u ) ;
(H13) For every t [ v 1 , b + v + 1 ] N v 1 and u , v R + , there is a constant δ 0 > 0 such that f ( t , u , v ) δ 0 g ( t , u ) .
Theorem 3. 
Suppose that (H10)–(H13) hold. Then, we get
(T1) There are y 0 , x 0 P h and r ( 0 , 1 ) such that r x 0 y 0 < x 0 ,
y 0 ( t ) s = 0 b + 2 G ( t , s ) [ f ( s + v 1 , y 0 ( s + v 1 ) , x 0 ( s + v 1 ) ) + g ( s + v 1 , y 0 ( s + v 1 ) ) ] ,
and
x 0 ( t ) s = 0 b + 2 G ( t , s ) [ f ( s + v 1 , x 0 ( s + v 1 ) , y 0 ( s + v 1 ) ) + g ( s + v 1 , x 0 ( s + v 1 ) ) ] ,
where h ( t ) = t v 1 ̲ ;
(T2) (22) has a unique positive solution y * P h ;
(T3) For each initial value x 0 , y 0 P h , the sequences
x n = s = 0 b + 2 G ( t , s ) [ f ( s + v 1 , x n 1 ( s + v 1 ) , y n 1 ( s + v 1 ) ) + g ( s + v 1 , x n 1 ( s + v 1 ) ) ] ,
y n = s = 0 b + 2 G ( t , s ) [ f ( s + v 1 , y n 1 ( s + v 1 ) , x n 1 ( s + v 1 ) ) + g ( s + v 1 , y n 1 ( s + v 1 ) ) ] ,
converge to y * as n .
Proof. 
From (H10) and (H11), we know that A : P × P P is a mixed monotone operator and B : P P is an increasing operator. Using (H12), for all γ ( 0 , 1 ) and x , y P , we obtain
A ( γ y , γ 1 x ) ( t ) = s = 0 b + 2 G ( t , s ) f ( s + v 1 , γ y ( s + v 1 ) , γ 1 x ( s + v 1 ) ) γ ξ s = 0 b + 2 G ( t , s ) f ( s + v 1 , y ( s + v 1 ) , x ( s + v 1 ) ) = γ ξ A ( y , x ) ( t ) ,
and hence A satisfies (7) in Lemma 12. In addition, for any y P and γ ( 0 , 1 ) we find
( B γ y ) ( t ) = s = 0 b + 2 G ( t , s ) g ( s + v 1 , γ y ( s + v 1 ) ) γ s = 0 b + 2 G ( t , s ) g ( s + v 1 , y ( s + v 1 ) ) = γ ( B y ) ( t ) .
Thus, B is a sub-homogeneous operator.
Let h 0 = h = t v 1 ̲ , and h 0 P h . From Lemma 4, we have
A ( h 0 , h 0 ) ( t ) = s = 0 b + 2 G ( t , s ) f ( s + v 1 , ( s + v 1 ) v 1 ̲ , ( s + v 1 ) v 1 ̲ ) s = 0 b + 2 t v 1 ̲ ( v + b β s + 1 ) v β 1 ̲ Γ ( v ) ( v + b β + 2 ) v β 1 ̲ f ( s + v 1 , ( s + v 1 ) v 1 ̲ , ( s + v 1 ) v 1 ̲ ) s = 0 b + 2 t v 1 ̲ ( v + b β s + 1 ) v β 1 ̲ Γ ( v ) ( v + b β + 2 ) v β 1 ̲ f ( s + v 1 , ( b + v + 1 ) v 1 ̲ , 0 ) = s = 0 b + 2 ( v + b β s + 1 ) v β 1 ̲ Γ ( v ) ( v + b β + 2 ) v β 1 ̲ f ( s + v 1 , ( b + v + 1 ) v 1 ̲ , 0 ) · h 0 ,
and
A ( h 0 , h 0 ) ( t ) = s = 0 b + 2 G ( t , s ) f ( s + v 1 , ( s + v 1 ) v 1 ̲ , ( s + v 1 ) v 1 ̲ ) s = 0 b + 2 t v 1 ̲ G ( b + v + 1 , s ) ( b + v + 1 ) v 1 ̲ f ( s + v 1 , ( s + v 1 ) v 1 ̲ , ( s + v 1 ) v 1 ̲ ) s = 0 b + 2 t v 1 ̲ G ( b + v + 1 , s ) ( b + v + 1 ) v 1 ̲ f ( s + v 1 , 0 , ( b + v + 1 ) v 1 ̲ ) = s = 0 b + 2 G ( b + v + 1 , s ) ( b + v + 1 ) v 1 ̲ f ( s + v 1 , 0 , ( b + v + 1 ) v 1 ̲ ) · h 0 .
Let l = s = 0 b + 2 G ( b + v + 1 , s ) ( b + v + 1 ) v 1 ̲ f ( s + v 1 , 0 , ( b + v + 1 ) v 1 ̲ ) , L = s = 0 b + 2 ( v + b β s + 1 ) v β 1 ̲ Γ ( v ) ( v + b β + 2 ) v β 1 ̲ f ( s + v 1 , ( b + v + 1 ) v 1 ̲ , 0 ) . Then, we have l h 0 A ( h 0 , h 0 ) L h 0 , i.e., A ( h 0 , h 0 ) P h 0 . Similarly, from (H11), we have
t v 1 ̲ s = 0 b + 2 G ( b + v + 1 , s ) ( b + v + 1 ) v 1 ̲ g ( s + v 1 , 0 ) ( B h 0 ) ( t ) = s = 0 b + 2 G ( t , s ) g ( s + v 1 , ( s + v 1 ) v 1 ̲ ) t v 1 ̲ s = 0 b + 2 ( v + b β s + 1 ) v β 1 ̲ Γ ( v ) ( v + b β + 2 ) v β 1 ̲ g ( s + v 1 , ( b + v + 1 ) v 1 ̲ ) .
Thus, we obtain B h 0 P h 0 . Therefore, (C1) in Lemma 12 holds.
Finally, for every x , y P , from (H13) we have
A ( y , x ) ( t ) = s = 0 b + 2 G ( t , s ) f ( s + v 1 , y ( s + v 1 ) , x ( s + v 1 ) ) δ 0 s = 0 b + 2 G ( t , s ) g ( s + v 1 , y ( s + v 1 ) ) = δ 0 ( B y ) ( t ) .
Thus, (C2) in Lemma 12 holds. Then, our conclusions are true from Lemma 12. □

4. Examples

In this section, we will provide some examples to verify our main results.
Example 1. 
Let g ( t , ψ ) = a | ψ | b M ( ψ ) , M ( ψ ) = ln ( | ψ | + 1 ) , ψ R , t [ v 1 , b + v + 1 ] N v 1 , where a κ 1 1 , + and b a , a + κ 2 1 . Then, lim | ψ | + M ( ψ ) | ψ | = 0 , and lim | ψ | + a | ψ | b M ( ψ ) | ψ | = a > κ 1 1 , lim | ψ | 0 + | a | ψ | b M ( ψ ) | | ψ | = | a b | < κ 2 1 . Therefore, (H1)–(H5) hold.
Example 2. 
Let b = 4 , v = 2.5 , α = 0.5 , β = 1.4 . Then, [ v 1 , b + v + 1 ] N v 1 = { 1.5 , 2.5 , 3.5 , 4.5 , 5.5 , 6.5 , 7.5 } ,   [ 0 , b + 2 ] N 0 = { 0 , 1 , 2 , 3 , 4 , 5 , 6 } , q 0 = min t [ v 1 , b + v + 1 ] N v 1 t 1.5 ̲ 7 . 5 1.5 ̲ = 0.068 , and if r = 6.5 we also obtain
s = 0 6 G ( 7.5 , s ) = s = 0 6 7 . 5 1.5 ̲ ( 6.1 s ) 0.1 ̲ 7 . 1 0.1 ̲ 5 . 5 1.5 ̲ = 112.26 ,
t = r b + v + 1 G ( b + v + 1 , t v + 1 ) = t = r + 1 v b + 2 G ( b + v + 1 , t ) = t = 5 6 G ( 7.5 , t ) = 19.86 .
Let a ˜ = 1 , b ˜ = 4 , c ˜ = 24 , and
g ( t , ψ ) = 0.008 , ψ [ 0 , 1 ] , ψ 0.992 , ψ [ 1 , 4 ] , 3.008 , ψ [ 4 , + ) .
Then, g satisfies
(I) g ( t , ψ ) < c ˜ q 0 s = 0 b + 2 G ( b + v + 1 , s ) = 3.144 , for t [ v 1 , b + v + 1 ] N v 1 , ψ [ 0 , 5190.3 ] ;
(II) g ( t , ψ ) > b ˜ q 0 t = r b + v + 1 G ( b + v + 1 , t v + 1 ) = 2.96 for t [ r , b + v + 1 ] N v 1 , ψ [ 4 , 865.1 ] ;
(III) g ( t , ψ ) < a ˜ s = 0 b + 2 G ( b + v + 1 , s ) = 0.009 for t [ v 1 , b + v + 1 ] N v 1 , ψ [ 0 , 1 ] .
Therefore, (H6)–(H9) hold.
Example 3. 
Let f ( t , u , v ) = ( b + v + 1 t ) 1 3 t 2 3 u 1 3 + v 1 5 , g ( t , u ) = ( b + v + 1 t ) 1 3 t 2 3 u 1 3 , ( t , u , v ) [ v 1 , b + v + 1 ] N v 1 × R + × R + . Then, f is increasing about u and decreasing about v , and g is increasing about u . For any γ ( 0 , 1 ) , u , v R + , taking ξ = 1 2 , then γ ξ ( γ , 1 ) and we obtain
f t , γ u , γ 1 v = ( b + v + 1 t ) 1 3 t 2 3 ( γ u ) 1 3 + γ 1 v 1 5 = γ 1 3 ( b + v + 1 t ) 1 3 t 2 3 u 1 3 + γ 1 5 v 1 5 γ 1 2 ( b + v + 1 t ) 1 3 t 2 3 u 1 3 + v 1 5 = γ ξ f ( t , u , v ) ,
and
g ( t , γ u ) = ( b + v + 1 t ) 1 3 t 2 3 ( γ u ) 1 3 = γ 1 3 ( b + v + 1 t ) 1 3 t 2 3 u 1 3 γ ( b + v + 1 t ) 1 3 t 2 3 u 1 3 = γ g ( t , u ) .
Moreover, it is easy to see that f ( t , u , v ) g ( t , u ) for ( t , u , v ) [ v 1 , b + v + 1 ] N v 1 × R + × R + . Therefore, (H10)–(H13) hold.
Example 4. 
In [31], the authors consider nonlinearities like:
f ( t , u , v ) = u 1 4 + [ v + 2 ] 1 3 + b ( t ) + d , g ( t , u ) = u 1 + u a ( t ) + c d , ( t , u , v ) [ v 1 , b + v + 1 ] N v 1 × R + × R + ,
where a , b : [ v 1 , b + v + 1 ] N v 1 R + with a 0 , and c , d are positive constants with c > d > 0 . Note that f is increasing about u and decreasing about v , and g is increasing about u . Moreover, for γ ( 0 , 1 ) , t [ v 1 , b + v + 1 ] N v 1 , u , v s . R + , we have
g ( t , γ u ) = γ u 1 + γ u a ( t ) + c d γ u 1 + u a ( t ) + γ ( c d ) = γ g ( t , u ) ,
and
f t , γ u , γ 1 v = γ 1 4 u 1 4 + γ 1 3 [ v + 2 γ ] 1 3 + b ( t ) + d γ 1 3 u 1 4 + [ v + 2 ] 1 3 + b ( t ) + d = γ ξ f ( t , u , v ) .
Furthermore, we note that
f ( t , u , v ) d = d max t [ v 1 , b + v + 1 ] N v 1 a ( t ) + c d × max t [ v 1 , b + v + 1 ] N v 1 a ( t ) + c d d max t [ v 1 , b + v + 1 ] N v 1 a ( t ) + c d × u 1 + u a ( t ) + c d : = δ 0 g ( t , u ) , ( t , u , v ) [ v 1 , b + v + 1 ] N v 1 × R + × R + .
Therefore, (H10)–(H13) hold.

5. Conclusions

Fractional-order difference equations are a new form of differential equation that have wider applications compared to traditional integer-order differential equations. They are generalized differential equations whose derivative index can be a decimal or a fraction, rather than an integer. This form of differential equation has wide applications in fields such as physics, engineering, and finance. Therefore, the importance of studying fractional difference equations is now becoming apparent. In this paper, we consider a boundary value problem with a fractional-order difference equation and use Green’s function to express its solution. Moreover, we obtain some existence theorems for the considered problem, i.e., when the nonlinearities satisfy some appropriate conditions, we study the existence, uniqueness, and multiplicity of solutions via the topological degree and fixed point theorems. Finally, we provide some examples to verify our main results.

Author Contributions

Conceptualization, Z.L.; methodology, Z.L., C.W., D.O. and J.X.; validation, Z.L., C.W., D.O. and J.X.; formal analysis, Z.L., C.W., D.O. and J.X.; investigation, D.O. and J.X. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Suqian Sci and Tech Program (grant no. K202134) anf the Natural Science Foundation of Chongqing (grant no. cstc2019jcyj-msxmX0390).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Atıcı, F.M.; Atıcı, M.; Belcher, M.; Marshall, D. A new approach for modeling with discrete fractional equations. Fund. Inform. 2017, 151, 313–324. [Google Scholar] [CrossRef]
  2. Atıcı, F.M.; Senguel, S. Modeling with fractional difference equations. J. Math. Anal. Appl. 2010, 369, 1–9. [Google Scholar] [CrossRef] [Green Version]
  3. Wu, G.C.; Baleanu, D. Discrete fractional logistic map and its chaos. Nonlinear Dyn. 2014, 75, 283–287. [Google Scholar] [CrossRef]
  4. Zheng, Z.X. On the developments and applications of fractional differential equations. J. Xuzhou Normal Univ. 2008, 26, 1–10. [Google Scholar]
  5. Goodrich, C.S.; Peterson, A.C. Discrete Fractional Calculus; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  6. Cheng, J. The Theory of Fractional-Order Difference Equations; Xiamen University Press: Xiamen, China, 2011. [Google Scholar]
  7. Rehman, M.; Lqbal, F.; Seemab, A. On existence of positive solutions for a class of discrete fractional boundary value problems. Positivity 2017, 21, 1173–1187. [Google Scholar] [CrossRef]
  8. Zhao, Y.; Sun, S.; Zhang, Y. Existence and uniqueness of solutions to a fractional difference equation with p-Laplacian operator. J. Appl. Math. Comput. 2017, 54, 183–197. [Google Scholar] [CrossRef]
  9. Xu, J.; Goodrich, C.S.; Cui, Y. Positive solutions for a system of first-order discrete fractional boundary value problems with semipositone nonlinearities. Rev. R. Acad. Cienc. Exactas Fís. Nat. Ser. A Mater. 2019, 113, 1343–1358. [Google Scholar] [CrossRef]
  10. Sitthiwirattham, T.; Tariboon, J.; Ntouyas, S.K. Boundary value problems for fractional difference equations with three-point fractional sum boundary conditions. Adv. Differ. Equ. 2013, 2013, 296. [Google Scholar] [CrossRef] [Green Version]
  11. Xu, J.; O’Regan, D.; Hou, C.; Cui, Y. Positive solutions for a class of fractional difference boundary value problems. Differ. Equ. Appl. 2017, 9, 479–493. [Google Scholar]
  12. Cheng, W.; Xu, J.; O’Regan, D.; Cui, Y. Positive solutions for a nonlinear discrete fractional boundary value problems with a p-Laplacian operator. J. Anal. Appl. Comput. 2019, 9, 1959–1972. [Google Scholar]
  13. Cheng, W.; Xu, J.; Cui, Y.; Ge, Q. Positive solutions for a class of fractional difference systems with coupled boundary conditions. Adv. Differ. Equ. 2019, 2019, 249. [Google Scholar] [CrossRef]
  14. Xu, J.; O’Regan, D. Existence and uniqueness of solutions for a first-order discrete fractional boundary value problem. Rev. R. Acad. Cienc. Exactas Fís. Nat. Ser. A Mat. 2018, 112, 1005–1016. [Google Scholar] [CrossRef]
  15. Bohner, M.; Fewster-Young, N. Discrete fractional boundary value problems and inequalities. Fract. Calc. Appl. Anal. 2021, 24, 1777–1796. [Google Scholar] [CrossRef]
  16. Eralp, B.; Topal, F.S. Monotone method for discrete fractional boundary value problems. Int. J. Nonlinear Anal. Appl. 2022, 13, 1989–1997. [Google Scholar]
  17. Bourguiba, R.; Cabada, A.; Kalthoum, W.O. Existence of solutions of discrete fractional problem coupled to mixed fractional boundary conditions. Rev. Real Acad. Cienc. Exactas Fis. Nat. Ser. A Mat. 2022, 116, 175. [Google Scholar] [CrossRef]
  18. Gopal, N.S.; Jonnalagadda, J.M. Positive solutions of nabla fractional boundary value problem. Cubo Math. J. 2022, 24, 467–484. [Google Scholar] [CrossRef]
  19. Dahal, R.; Duncan, D.; Goodrich, C.S. Systems of semipositone discrete fractional boundary value problems. J. Differ. Equ. Appl. 2014, 20, 473–491. [Google Scholar] [CrossRef]
  20. Goodrich, C.S. Systems of discrete fractional boundary value problems with nonlinearities satisfying no growth conditions. J. Differ. Equ. Appl. 2015, 21, 437–453. [Google Scholar] [CrossRef]
  21. Lv, W. Existence of solutions for discrete fractional boundary value problems with a p-Laplacian operator. Adv. Differ. Equ. 2012, 2012, 163. [Google Scholar] [CrossRef] [Green Version]
  22. Ferreira, R.A.C. Existence and uniqueness of solution to some discrete fractional boundary value problems of order less than one. J. Differ. Equ. Appl. 2013, 19, 712–718. [Google Scholar] [CrossRef]
  23. Goodrich, C.S. On a first-order semipositone discrete fractional boundary value problem. Arch. Math. 2012, 99, 509–518. [Google Scholar] [CrossRef]
  24. Goodrich, C.S. On discrete sequential fractional boundary value problems. J. Math. Anal. Appl. 2012, 385, 111–124. [Google Scholar] [CrossRef] [Green Version]
  25. Sitthiwirattham, T. Boundary value problem for p-Laplacian Caputo fractional difference equations with fractional sum boundary conditions. Math. Methods Appl. Sci. 2016, 39, 1522–1534. [Google Scholar] [CrossRef]
  26. Kunnawuttipreechachan, E.; Promsakon, C.; Sitthiwirattham, T. Nonlocal fractional sum boundary value problem for a coupled system of fractional sum-difference equations. Dyn. Syst. Appl. 2019, 28, 73–92. [Google Scholar] [CrossRef]
  27. Guo, D.; Lakshmikantham, V. Nonlinear Problems in Abstract Cones; Academic Press: Orlando, FL, USA, 1988. [Google Scholar]
  28. Zhang, W.; Ni, J. New multiple positive solutions for Hadamard-type fractional differential equations with nonlocal conditions on an infinite interval. Appl. Math. Lett. 2021, 118, 107165. [Google Scholar] [CrossRef]
  29. Xu, J.; Jiang, J.; O’Regan, D. Positive solutions for a class of p-Laplacian Hadamard fractional-order three-point boundary value problems. Mathematics 2020, 8, 308. [Google Scholar] [CrossRef] [Green Version]
  30. Zhai, C.; Hao, M. Fixed point theorems for mixed monotone operators with perturbation and applications to fractional differential equation boundary value problems. Nonlinear Anal. 2012, 75, 2542–2551. [Google Scholar] [CrossRef]
  31. Zhai, C.; Hao, M. Mixed monotone operator methods for the existence and uniqueness of positive solutions to Riemann–Liouville fractional differential equation boundary value problems. Bound. Value Probl. 2013, 2013, 85. [Google Scholar] [CrossRef] [Green Version]
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

Lv, Z.; Wu, C.; O’Regan, D.; Xu, J. Solvability of a Boundary Value Problem Involving Fractional Difference Equations. Axioms 2023, 12, 650. https://doi.org/10.3390/axioms12070650

AMA Style

Lv Z, Wu C, O’Regan D, Xu J. Solvability of a Boundary Value Problem Involving Fractional Difference Equations. Axioms. 2023; 12(7):650. https://doi.org/10.3390/axioms12070650

Chicago/Turabian Style

Lv, Zhiwei, Chun Wu, Donal O’Regan, and Jiafa Xu. 2023. "Solvability of a Boundary Value Problem Involving Fractional Difference Equations" Axioms 12, no. 7: 650. https://doi.org/10.3390/axioms12070650

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