Next Article in Journal
GIS Partial Discharge Pattern Recognition Based on Multi-Feature Information Fusion of PRPD Image
Next Article in Special Issue
A Method for Performing the Symmetric Anti-Difference Equations in Quantum Fractional Calculus
Previous Article in Journal
Interpolating Scaling Functions Tau Method for Solving Space–Time Fractional Partial Differential Equations
Previous Article in Special Issue
On q-Limaçon Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Positive Solutions for a Class of Integral Boundary Value Problem of Fractional q-Difference Equations

1
School of Mathematics and Statistics, Shanxi Datong University, Datong 037009, China
2
School of Mathematics and Computer Science, Shanxi Normal University, Taiyuan 030032, China
3
Department of Mathematics, Trent University, Peterborough, ON K9L 0G2, Canada
4
Department of Computer Science, Trent University, Peterborough, ON K9L 0G2, Canada
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(11), 2465; https://doi.org/10.3390/sym14112465
Submission received: 16 August 2022 / Revised: 7 November 2022 / Accepted: 7 November 2022 / Published: 21 November 2022
(This article belongs to the Special Issue Symmetry in Quantum Calculus)

Abstract

:
This paper studies a class of integral boundary value problem of fractional q-difference equations. We first give an explicit expression for the associated Green’s function and obtain an important property of the function. The new property allows us to prove sufficient conditions for the existence of positive solutions based on the associated parameter. The results are derived from the application of a fixed point theorem on order intervals.

1. Introduction

It is known that fractional differential equations and their discrete analogues, fractional difference equations, have broad applications in many areas such as neural computing [1], transportation modelling [2], and dynamical systems [3,4]. Compared to the traditional integer ordered calculus, fractional ordered derivatives have less requirements for the smoothness of the functions [5]. This property provides some advantages in modeling computational algorithms and other applications [6,7]. From the perspectives of numerical analysis, fractional integral equations can be solved via various neural networks [8].
The concept of quantum calculus (q-calculus) was first developed by Jackson [9,10]. The theory has advanced applications in mathematical physics, impulsive waves, and signal analysis [11,12,13,14,15]. Fractional q-difference calculus is the counterpart of fractional integral and difference operators [16,17]. For a comprehensive introduction to this subject, we refer to the book [12].
The study of positive solutions for Boundary Value Problems (BVPs) is paramount due to its significant applications. For example, positive solutions of a BVP arising from chemical reactor theory [18] represent the temperature during the reaction. In the literature, existence of positive solutions for differential equations (DE) and fractional differential equations (FDE) have been extensively studied [19,20]. However, results for q-difference fractional boundary value problems are relatively rare [21,22].
The following nonlinear q-fractional BVP was considered in [21]:
( D q α y ) ( t ) = f ( x , y ( x ) ) , 0 < x < 1 , y ( 0 ) = D q y ( 0 ) = 0 , D q y ( 1 ) = β 0 .
A common approach to prove the existence of positive solutions for BVPs is to convert the problem to an integral equation and then apply fixed point theorems. Results of [21] on sufficient conditions for the existence of positive solutions were obtained following this idea. The well-known Guo–Krasnoselskii fixed point theorem was applied.
Later, as a special case of the boundary conditions for BVP (1), when β = 0 , the following parameter dependent fractional-difference BVP was examined in [23]:
( D q α u ) ( x ) + λ h ( x ) f ( u ( x ) ) = 0 , 0 < x < 1 , u ( 0 ) = D q u ( 0 ) = D q u ( 1 ) = 0 .
BVP (2) is equivalent to a class of time-independent fractional-difference Schrödinger equations that simulate the evolution of the slowly varying amplitude of a nonlinear wave [23]. Again, applying the Guo–Krasnoselskii fixed point theorem, the existence of positive solutions depending on the values of the parameter λ were obtained.
Recently, the uniqueness of solutions for BVP involving a parameter and the integral operator was obtained in [24]:
( D q α u ) ( t ) + f ( t , u ( t ) ) = 0 , t ( 0 , 1 ) , u ( 0 ) = D q u ( 0 ) = 0 , u ( 1 ) = λ 0 1 u ( s ) d q s .
The approach of [24] relies on the iterative technique and a fixed point theorem that ensures uniqueness [25].
Motivated by the above studies, in this paper, we consider the fractional q-difference equation associated with a parameter and subject to integral boundary conditions:
( D q α u ) ( t ) + λ f ( u ( t ) ) = 0 , 0 < t < 1 , u ( 0 ) = D q u ( 0 ) = 0 , D q u ( 1 ) = [ α 1 ] q 0 1 u ( s ) d q s ,
where the parameters λ > 0 , 2 < α 3 , and f : [ 0 , + ) [ 0 , + ) is a continuous function. D q α denotes the q-fractional differential operator of order α .
The same as the work of [21,23,24], we will use the Riemann–Liouville type fractional derivative [26]. Different from methods applied in [21,23,24], we utilize a relatively new fixed point theorem on order intervals [27] that provides more information on the solutions.
In Section 2, we present some definitions and lemmas for fractional q-derivative and fractional q-integral. The Green’s function of the fractional integral BVP is obtained. Its properties are proved. In Section 3, new sufficient conditions for the existence of positive solutions are proved.

2. Preliminaries on q -Calculus

The following notations and definitions can be found in [21,28]. Suppose q ( 0 , 1 ) and define
[ a ] q = 1 q a 1 q , a R .
If α R , then
( a b ) ( α ) = a α n = 0 a b q n a b q α + n .
The q-gamma function is defined as follows:
Γ q ( x ) = ( 1 q ) ( x 1 ) ( 1 q ) x 1 , x R \ { 0 , 1 , 2 , } .
The q-derivative of a function f is defined by
( D q f ) ( x ) = f ( x ) f ( q x ) ( 1 q ) x .
The following four formulas will be used in the sequel:
[ a ( t s ) ] ( α ) = a α ( t s ) ( α ) ,
( a b ) ( α ) = ( a b q α 1 ) ( a b ) ( α 1 ) ,
t D q ( t s ) ( α ) = [ α ] q ( t s ) ( α 1 ) ,
t D q 0 t f ( t , s ) d q s ( t ) = 0 t t D q f ( t , s ) d q s + f ( q t , t ) ,
where t D q denotes the derivative with respect to variable t.
Remark 1
([22]). If α > 0 and a b t , then ( t a ) ( α ) ( t b ) ( α ) .
Definition 1
([17]). Assume that α 0 and f is a function defined on [ 0 , 1 ] . The fractional q-integral of the Riemann–Liouville type is defined by
( I q α f ) ( t ) = 1 Γ ( α ) 0 t ( t q s ) ( α 1 ) f ( s ) d q s , α > 0 , t [ 0 , 1 ] .
Definition 2
([14]). The fractional q-derivative of the Riemann–Liouville type of f is defined by
( D q α f ) ( t ) = ( D q m I q m α f ) ( t ) , α > 0 ,
where m is the smallest integer greater than or equal to α.
Lemma 1
([11,14]). Let α , β 0 and f be a function defined on [ 0 , 1 ] . We have the following formulas:
( I q β I q α f ) ( t ) = ( I q α + β f ) ( t ) ,
( D q α I q α f ) ( t ) = f ( t ) .
Lemma 2
([22]). Assuming that α > 0 and p is a positive integer, then
( I q α D q p f ) ( t ) = ( D q p I q α f ) ( t ) k = 0 p 1 x α p + k Γ q ( α + k p + 1 ) ( D q k f ) ( 0 ) .
Lemma 3.
Assume that parameter λ > 0 , 2 < α 3 , f : [ 0 , + ) [ 0 , + ) is a continuous function. Then, the solution of integral boundary value problem
( D q α u ) ( t ) + λ f ( u ( t ) ) = 0 , 0 < t < 1 , u ( 0 ) = D q u ( 0 ) = 0 , D q u ( 1 ) = [ α 1 ] q 0 1 u ( s ) d q s
is equivalent to the solution of
u ( t ) = λ 0 1 G ( t , q s ) f ( u ( s ) ) d q s , t [ 0 , 1 ] ,
where
G ( t , q s ) = t α 1 ( 1 q s ) ( α 2 ) [ α ] q ( 1 s q α ) ( 1 s q α 1 ) [ α ] q 1 ( t q s ) ( α 1 ) ( [ α ] q 1 ) Γ q ( α ) , 0 q s t 1 ; t α 1 ( 1 q s ) ( α 2 ) [ α ] q ( 1 s q α ) ( 1 s q α 1 ) ( [ α ] q 1 ) Γ q ( α ) , 0 t q s 1 .
Proof. 
In view of Definition 2 and Lemma 1, we obtain that
( D q α u ) ( t ) + λ f ( u ( t ) ) = 0 ( I q α D q 3 I q 3 α u ) ( t ) = λ I q α f ( u ( t ) ) .
It follows from Lemma 2 that the solution u of (9) is given by
u ( t ) = c 1 t α 1 + c 2 t α 2 + c 3 t α 3 λ Γ q ( α ) 0 t ( t q s ) ( α 1 ) f ( u ( s ) ) d q s ,
where c 1 , c 2 , c 3 R . Since u ( 0 ) = 0 , we obtain c 3 = 0 .
Taking the derivative of both sides of Equation (10), with the help of (7) and (8), we obtain
( D q u ) ( t ) = [ α 1 ] q c 1 t α 2 + [ α 2 ] q c 2 t α 3 λ Γ q ( α ) 0 t [ α 1 ] q ( t q s ) ( α 2 ) f ( u ( s ) ) d q s .
Using the boundary conditions of D q u ( 0 ) = 0 and D q u ( 1 ) = [ α 1 ] q 0 1 u ( s ) d q s , we have c 2 = 0 and
c 1 = λ Γ q ( α ) 0 1 ( 1 q s ) ( α 2 ) f ( u ( s ) ) d q s + 0 1 u ( s ) d q s .
Thus, u can be calculated as
u ( t ) = λ Γ q ( α ) 0 t ( t q s ) ( α 1 ) f ( u ( s ) ) d q s + λ t α 1 Γ q ( α ) 0 1 ( 1 q s ) ( α 2 ) f ( u ( s ) ) d q s + t α 1 0 1 u ( s ) d q s .
Let C = 0 1 u ( s ) d q s . Integrating (11) with respect to t from 0 to 1, we obtain
C = λ Γ q ( α ) 0 1 d q t 0 t ( t q s ) ( α 1 ) f ( u ( s ) ) d q s + λ Γ q ( α ) 0 1 t α 1 d q t 0 1 ( 1 q s ) ( α 2 ) f ( u ( s ) ) d q s + C 0 1 t α 1 d q t .
By using equations 0 1 t α 1 d q t = 1 [ α ] q and
0 1 d q t 0 t ( t q s ) ( α 1 ) f ( u ( s ) ) d q s = 1 [ α ] q 0 1 ( 1 q s ) ( α ) f ( u ( s ) ) d q s ,
We obtain
C = λ ( [ α ] q 1 ) Γ q ( α ) 0 1 ( 1 q s ) ( α ) f ( u ( s ) ) d q s + λ ( [ α ] q 1 ) Γ q ( α ) 0 1 ( 1 q s ) ( α 2 ) f ( u ( s ) ) d q s .
In view of (6), substituting C into the Formula (11), we deduce that
u ( t ) = λ Γ q ( α ) 0 t ( t q s ) ( α 1 ) f ( u ( s ) ) d q s + λ t α 1 Γ q ( α ) 0 1 ( 1 q s ) ( α 2 ) f ( u ( s ) ) d q s λ t α 1 ( [ α ] q 1 ) Γ q ( α ) 0 1 ( 1 q s ) ( α ) f ( u ( s ) ) d q s + λ t α 1 ( [ α ] q 1 ) Γ q ( α ) 0 1 ( 1 q s ) ( α 2 ) f ( u ( s ) ) d q s = λ Γ q ( α ) 0 t ( t q s ) ( α 1 ) f ( u ( s ) ) d q s + λ t α 1 ( [ α ] q 1 ) Γ q ( α ) 0 1 ( 1 q s ) ( α 2 ) [ α ] q ( 1 s q α ) ( 1 s q α 1 ) f ( u ( s ) ) d q s = λ 0 1 G ( t , q s ) f ( u ( s ) ) d q s .
Lemma 4.
The above function G satisfies
(i) 
G ( t , q s ) 0 , t , s [ 0 , 1 ] ;
(ii) 
t α 1 G ( 1 , q s ) G ( t , q s ) G ( 1 , q s ) , t , s [ 0 , 1 ] .
Proof. 
We first define two functions
g 1 ( t , q s ) = t α 1 ( 1 q s ) ( α 2 ) [ α ] q ( 1 s q α ) ( 1 s q α 1 ) ( [ α ] q 1 ) ( t q s ) ( α 1 ) ,
g 2 ( t , q s ) = t α 1 ( 1 q s ) ( α 2 ) [ α ] q ( 1 s q α ) ( 1 s q α 1 ) .
In view of Remark 1, we have
g 1 ( t , q s ) = t α 1 ( 1 q s ) ( α 2 ) ( [ α ] q 1 ) + [ 1 ( 1 s q α ) ( 1 s q α 1 ) ] ( [ α ] q 1 ) ( t q s ) ( α 1 ) = t α 1 ( 1 q s ) ( α 2 ) { ( [ α ] q 1 ) + [ 1 ( 1 s q α ) ( 1 s q α 1 ) ] } t α 1 ( [ α ] q 1 ) ( 1 q s t ) ( α 1 ) t α 1 ( 1 q s ) ( α 2 ) { ( [ α ] q 1 ) + [ 1 ( 1 s q α ) ( 1 s q α 1 ) ] } t α 1 ( [ α ] q 1 ) ( 1 q s ) ( α 1 ) = t α 1 ( [ α ] q 1 ) [ ( 1 q s ) ( α 2 ) ( 1 q s ) ( α 1 ) ] + t α 1 ( 1 q s ) ( α 2 ) [ 1 ( 1 s q α ) ( 1 s q α 1 ) ] 0 .
Notice that [ α ] q 1 > 0 , thus
g 2 ( t , q s ) = t α 1 ( 1 q s ) ( α 2 ) { ( [ α ] q 1 ) + [ 1 ( 1 s q α ) ( 1 s q α 1 ) ] } t α 1 ( 1 q s ) ( α 2 ) [ 1 ( 1 s q α ) ( 1 s q α 1 ) ] 0 .
Therefore, G ( t , q s ) 0 . In addition, for fixed t [ 0 , 1 ] ,
t D q g 1 ( t , q s ) = [ α 1 ] q t α 2 ( 1 q s ) ( α 2 ) { [ α ] q ( 1 s q α ) ( 1 s q α 1 ) } ( [ α ] q 1 ) [ α 1 ] q ( t q s ) ( α 2 ) = [ α 1 ] q t α 2 ( 1 q s ) ( α 2 ) { [ α ] q ( 1 s q α ) ( 1 s q α 1 ) } ( [ α ] q 1 ) [ α 1 ] q t α 2 ( 1 q s t ) ( α 2 ) [ α 1 ] q t α 2 ( 1 q s ) ( α 2 ) { ( [ α ] q 1 ) + [ 1 ( 1 s q α ) ( 1 s q α 1 ) ] } [ α 1 ] q t α 2 ( 1 q s ) ( α 2 ) ( [ α ] q 1 ) = [ α 1 ] q t α 2 ( 1 q s ) ( α 2 ) [ 1 ( 1 s q α ) ( 1 s q α 1 ) ] 0 .
t D q g 2 ( t , q s ) = [ α 1 ] q t α 2 ( 1 q s ) ( α 2 ) { ( [ α ] q 1 ) + [ 1 ( 1 s q α ) ( 1 s q α 1 ) ] } [ α 1 ] q t α 2 ( 1 q s ) ( α 2 ) [ 1 ( 1 s q α ) ( 1 s q α 1 ) ] 0 .
i.e., g 1 ( t , q s ) and g 2 ( t , q s ) are increasing functions of t, so G ( t , q s ) is an increasing function to t for fixed s [ 0 , 1 ] . If t q s , then
G ( t , q s ) G ( 1 , q s ) = t α 1 ( 1 q s ) ( α 2 ) { [ α ] q ( 1 s q α ) ( 1 s q α 1 ) } ( [ α ] q 1 ) ( t q s ) ( α 1 ) ( 1 q s ) ( α 2 ) { [ α ] q ( 1 s q α ) ( 1 s q α 1 ) } ( [ α ] q 1 ) ( 1 q s ) ( α 1 ) t α 1 ( 1 q s ) ( α 2 ) { [ α ] q ( 1 s q α ) ( 1 s q α 1 ) } t α 1 ( [ α ] q 1 ) ( 1 q s ) ( α 1 ) ( 1 q s ) ( α 2 ) { [ α ] q ( 1 s q α ) ( 1 s q α 1 ) } ( [ α ] q 1 ) ( 1 q s ) ( α 1 ) = t α 1 .
If t q s , then
G ( t , q s ) G ( 1 , q s ) = t α 1 ( 1 q s ) ( α 2 ) { [ α ] q ( 1 s q α ) ( 1 s q α 1 ) } ( 1 q s ) ( α 2 ) { [ α ] q ( 1 s q α ) ( 1 s q α 1 ) } ( [ α ] q 1 ) ( 1 q s ) ( α 1 ) t α 1 ( 1 q s ) ( α 2 ) { [ α ] q ( 1 s q α ) ( 1 s q α 1 ) } ( 1 q s ) ( α 2 ) { [ α ] q ( 1 s q α ) ( 1 s q α 1 ) } = t α 1 .
Thus, the proof is complete. □

3. Main Results

For convenience, we denote
F 0 = lim sup u 0 + f ( u ) u , f 0 = lim inf u 0 + f ( u ) u ,
F = lim sup u + f ( u ) u , f = lim inf u + f ( u ) u ,
M = 0 1 G ( 1 , q s ) d q s , N = 0 1 G ( 1 , q s ) h ( s ) d q s ,
where h ( s ) = s α 1 .
Let X = C [ 0 , 1 ] be an ordered Banach space with the cone X + = { u x , u ( t ) 0 , t [ 0 , 1 ] } ; suppose u 0 X + with u 1 , and define
P u 0 = { x X + , x x u 0 } ;
then, P u 0 is a cone of X + .
Let u 0 = t α 1 , define P = { x X + , x t α 1 x } , and define the operator H : P X + as following
( H u ) ( t ) = λ 0 1 G ( t , q s ) f ( u ( s ) ) d q s .
For any r > 0 , we denote Ω r = { x X : x < r } and Ω r = { x X : x = r } .
Theorem 1
([27]). Assume that X is a Banach space with the cone X + . Let 0 u 0 φ be such that u 0 1 , φ = 1 satisfying the condition:
If x X + , x 1 , then x φ .
If there exist positive numbers 0 < a < b such that T : P u 0 ( Ω ¯ b \ Ω a ) P u 0 is a completely continuous operator and the conditions:
T ( x ) x [ a u 0 , a φ ] a a n d T ( x ) x [ b u 0 , b φ ] b
or
T ( x ) x [ a u 0 , a φ ] a a n d T ( x ) x [ b u 0 , b φ ] b
are satisfied, then T has a fixed point x 0 [ a u 0 , b φ ] .
Lemma 5.
Assume f : [ 0 , + ) [ 0 , + ) is a continuous function; then, H : P P is a completely continuous operator.
Proof. 
First, it is easy to verify that H C ( P , P ) .
Next, from Lemma 4(ii), we have
( H u ) ( t ) = λ 0 1 G ( t , q s ) f ( s , u ( s ) ) d q s λ t α 1 0 1 G ( 1 , q s ) f ( s , u ( s ) ) d q s = t α 1 H u .
That is, H ( P ) P .
Let D P be bounded. There exists a positive constant M > 0 such that u M for each u D . Assume L = max u M | f ( u ) | + 1 ; then, in view of Lemma 4(ii), we have
| ( H u ) ( t ) | = | λ 0 1 G ( t , q s ) f ( u ( s ) ) d q s | λ 0 1 | G ( t , q s ) f ( u ( s ) ) | d q s λ L 0 1 G ( 1 , q s ) d q s [ α ] q λ L ( [ α ] q 1 ) Γ q ( α ) 0 1 ( 1 q s ) ( α 2 ) d q s [ α ] q λ L ( [ α ] q 1 ) Γ q ( α )
This shows that, for each parameter λ > 0 , T ( D ) is bounded.
On the other hand, for all u D , when 0 < t 1 < t 2 < 1 , there exists a sufficiently small δ , such that | t 2 t 1 | < δ . One has | H u ( t 2 ) H u ( t 1 ) | < ε , that is to say, H ( Ω ) is equicontinuous. In fact, when t 1 t 2 , we have
| H u ( t 2 ) H u ( t 1 ) | λ 0 1 | G ( t 2 , q s ) G ( t 1 , q s ) | f ( u ( s ) ) d q s λ L 0 1 | G ( t 2 , q s ) G ( t 1 , q s ) | d q s = λ L ( [ α ] q 1 ) Γ q ( α ) ( 0 t 1 | ( 1 q s ) ( α 2 ) ( [ α ] q ( 1 s q α ) ( 1 s q α 1 ) ) ( t 2 α 1 t 1 α 1 ) ( [ α ] q 1 ) ( ( t 2 q s ) ( α 1 ) ( t 1 q s ) ( α 1 ) ) | d q s + t 1 t 2 | t 2 α 1 ( 1 q s ) ( α 2 ) ( [ α ] q ( 1 s q α ) ( 1 s q α 1 ) ) ( [ α ] q 1 ) ( t 2 q s ) ( α 1 ) t 1 α 1 ( 1 q s ) ( α 2 ) ( [ α ] q ( 1 s q α ) ( 1 s q α 1 ) ) | d q s + t 2 1 | ( 1 q s ) ( α 2 ) ( [ α ] q ( 1 s q α ) ( 1 s q α 1 ) ) ( t 2 α 1 t 1 α 1 ) | d q s ) λ L ( [ α ] q 1 ) Γ q ( α ) 0 t 1 2 δ d q s + t 1 t 2 ( δ + 1 ) d q s + t 2 1 δ d q s 0 .
Therefore, T ( D ) is equicontinuous. By means of the Arzela-Ascoli theorem, T : P P is completely continuous. □
It is well-known that u X is a solution of (4) if and only if u is a fixed point of H.
Theorem 2.
Assume that M F 0 < N f . Then, BVP (4) has at least one positive solution for λ 1 N f , 1 M F 0 .
Proof. 
From the definition of the cone P, we have u 0 = h ( t ) , φ = 1 . For u X + , | | u | | = u ( 1 ) . Therefore, u 0 and φ satisfy the conditions of Theorem 1.
Since λ < 1 M F 0 , we can choose a small enough number of ε 1 > 0 , such that
λ ( F 0 + ε 1 ) M < 1 .
From F 0 = lim sup u 0 + f ( u ) u , there exists δ > 0 such that
f ( u ) u < F 0 + ε 1 f o r u ( 0 , δ ) .
Let a = δ 2 . Then, for u [ a t α 1 , a ] , we can obtain
H ( u ) = H ( u ) ( 1 ) = λ 0 1 G ( 1 , q s ) f ( u ( s ) ) d q s λ ( F 0 + ε 1 ) 0 1 G ( 1 , q s ) u ( s ) d q s a .
On the other hand, using the condition λ > 1 N f , there exist c > 0 and ε 2 > 0 such that
λ ( f ε 2 ) c 1 G ( 1 , q s ) h ( s ) d q s > 1 .
Let K > 0 such that
f ( u ) u > f ε 2 for u K .
Set K 0 = max { K , c α 1 δ } and b = K 0 c α 1 . Then, for u [ b t α 1 , b ] ( t [ c , 1 ] ) . We can obtain
H ( u ) = H ( u ) ( 1 ) = λ 0 1 G ( 1 , q s ) f ( u ( s ) ) d q s λ ( f ε 2 ) c 1 G ( 1 , q s ) u ( s ) d q s λ ( f ε 2 ) b c 1 G ( 1 , q s ) h ( s ) d q s b .
By Theorem 1, H has a fixed point u λ [ a t α 1 , b ] . It is a positive solution of (4). □
Theorem 3.
Assuming that M F < N f 0 holds, then BVP (4) has at least one positive solution for λ 1 N f 0 , 1 M F .
Proof. 
From the definition of the cone P, we have u 0 = h ( t ) , φ = 1 . For u X + , | | u | | = u ( 1 ) . Therefore, u 0 and φ satisfy the conditions of Theorem 1.
Since λ > 1 N f 0 , we can choose a small enough number of ε 1 > 0 , such that
λ ( f 0 ε 1 ) N > 1 .
Let δ > 0 such that
f ( u ) u > f 0 ε 1 f o r u ( 0 , δ ) ,
Let a = δ 2 . Then, for u [ a t α 1 , a ] , we can obtain
H ( u ) = H ( u ) ( 1 ) = λ 0 1 G ( 1 , q s ) f ( u ( s ) ) d q s λ ( f 0 ε 1 ) 0 1 G ( 1 , q s ) u ( s ) d q s λ ( f 0 ε 1 ) a 0 1 G ( 1 , q s ) h ( s ) d q s a .
On the other hand, using the condition λ < 1 M F , there exists ε 2 > 0 such that
λ ( F + ε 2 ) M < 1 .
Let K > 0 such that f ( u ) u < F + ε 2 for u K .
Set K 0 = max { K , t α 1 δ } = b t α 1 ( t [ 0 , 1 ] ) ; then, for u [ b t α 1 , b ] ( t [ 0 , 1 ] ) , we have
H ( u ) = H ( u ) ( 1 ) = λ 0 1 G ( 1 , q s ) f ( u ( s ) ) d q s λ ( F + ε 2 ) 0 1 G ( 1 , q s ) u ( s ) d q s b .
By Theorem 1, H has a fixed point u λ [ a t α 1 , b ] . It is a positive solution of (4).
Setting α = 5 2 , q = 1 2 , a simple calculation shows that M 0.3868 ,   N 0.1159 .
Example 1.
Consider the following fractional q-difference boundary value problem
( D q α u ) ( t ) + λ u 2 + 4 u 3 1 + u 2 = 0 , 0 < t < 1 ; u ( 0 ) = D q u ( 0 ) = 0 , D q u ( 1 ) = [ α 1 ] q 0 1 u ( s ) d q s ,
where F 0 = 0.5 , f = 4.5 . Applying Theorem 2, if λ ( 1.92 , 5.17 ) 1 N f , 1 M F 0 , then the above boundary value problem has at least one positive solution.
Example 2.
Consider the following fractional q-difference boundary value problem
( D q α u ) ( t ) + λ u 1 2 + 2 1 + u 2 = 0 , 0 < t < 1 ; u ( 0 ) = D q u ( 0 ) = 0 , D q u ( 1 ) = [ α 1 ] q 0 1 u ( s ) d q s ,
where f 0 = 2.5 , F = 0.5 . Applying Theorem 3, if λ ( 3.45 , 5.17 ) 1 N f 0 , 1 M F , then the above boundary value problem has at least one positive solution.

4. Conclusions

We study a new class of fractional q-difference boundary value problem associated with a parameter and subject to boundary conditions that involve the integral operator. After the explicit expression for the associated Green’s function is given in Lemma 3, two important properties are proved in Lemma 4. In particular, the property (ii) of Lemma 4 shows that the Green’s function does not satisfy the strong positivity condition that is necessary for some popular methods in studying nonlinear boundary value problems [20]. The obtained results on existence of positive solutions are new. Applications are shown by two examples.

Author Contributions

Methodology, S.K., Y.Z., H.C. and W.F.; software, S.K.; validation, H.C.; formal analysis, S.K., Y.Z., H.C. and W.F.; resources, S.K. and W.F.; writing—original draft preparation, S.K. and Y.Z.; writing—review and editing, H.C. and W.F.; visualization, S.K. and H.C.; supervision, S.K. and W.F.; project administration, S.K. and W.F.; funding acquisition, S.K. and W.F.; All authors have read and agreed to the published version of the manuscript.

Funding

The project was supported by the Natural Sciences Foundation of China (No. 11871314) and the Natural Sciences and Engineering Research Council of Canada (NSERC).

Data Availability Statement

Data sharing is not applicable to this paper as no datasets were generated or analyzed during the current study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Antil, H.; Khatri, R.; Lohner, R.; Verma, D. Fractional deep neural network via constrained optimization. Mach. Learn. Sci. Technol. 2019, 2, 015003. [Google Scholar] [CrossRef]
  2. Graef, J.R.; Ho, S.S.; Kong, L.; Wang, M. A fractional differential equation model for bike share systems. J. Nonlinear Funct. Anal. 2019, 2019, 1830489. [Google Scholar]
  3. Jarad, F.; Abdeljawad, T.; Baleanu, D. Stability of q-fractional non-autonomous systems. Nonlinear Anal. Real World Appl. 2013, 14, 780–784. [Google Scholar] [CrossRef]
  4. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North–Holland Mathematics Studies, 204; Elsevier Science: Amsterdam, The Netherlands, 2006. [Google Scholar]
  5. Goodrich, C.S.; Peterson, A.C. Discrete Fractional Calculus, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2016. [Google Scholar]
  6. Atici, F.M.; Sengul, S. Modeling with fractional difference equations. J. Math. Anal. Appl. 2010, 369, 1–9. [Google Scholar] [CrossRef] [Green Version]
  7. Magin, R.L. Fractional Calculus in Bioengineering; Begell House: Danbury, CT, USA, 2006. [Google Scholar]
  8. Mall, S. Artificial neural network approach for solving fractional order applied problems. In New Paradigms in Computational Modeling and Its Applications; Academic Press: Cambridge, MA, USA, 2021; pp. 163–175. [Google Scholar]
  9. Jackson, F.H. On q-functions and a certain difference operator. Trans. R. Soc. Edinb. 1908, 46, 253–281. [Google Scholar] [CrossRef]
  10. Jackson, F.H. On q-definite integrals. Quart. J. Pure Appl. Math. 1910, 41, 193–203. [Google Scholar]
  11. Agarwal, R.P. Certain fractional q-integrals and q-derivative. Proc. Camb. Philos. Soc. 1969, 66, 365–370. [Google Scholar] [CrossRef]
  12. Annaby, M.; Mansour, Z. q-Fractional Calculus and Equations; Springer: New York, NY, USA, 2012. [Google Scholar]
  13. Abdel-Gawad, H.I.; Aldailami, A.A. On q-dynamic equations modelling and complexity. Appl. Math. Model. 2010, 34, 697–709. [Google Scholar] [CrossRef]
  14. Rajković, P.M.; Marinković, S.D.; Stanković, M.S. Fractional integrals and derivatives in q-calculus. Appl. Anal. Discret. Math. 2007, 1, 311–323. [Google Scholar]
  15. Tariboon, J.; Ntouyas, S.K.; Agarwal, P. New concepts of fractional quantum calculus and applications to impulsive fractional q-difference equations. Adv. Differ. Equ. 2015, 2015, 18. [Google Scholar] [CrossRef] [Green Version]
  16. Allouch, N.; Graef, J.R.; Hamani, S. Boundary Value Problem for Fractional q-Difference Equations with Integral Conditions in Banach Spaces. Fractal Fract. 2022, 6, 237. [Google Scholar] [CrossRef]
  17. Kac, V.; Cheung, P. Quantum Calculus; Springer: New York, NY, USA, 2002. [Google Scholar]
  18. Feng, W.; Zhang, G.; Chai, Y. Existence of positive solutions for second-order differential equations arising from chemical reactor theory. Discrete Contin. Dyn. Syst. 2007, 2007, 373–381. [Google Scholar]
  19. Cabada, A.; Hamdi, Z. Nonlinear fractional differential equations with integral boundary value conditions. J. Appl. Math. Comput. 2014, 228, 251–257. [Google Scholar] [CrossRef]
  20. Webb, J.R.L. A class of positive linear operators and applications to nonlinear boundary value problems. Topol. Methods Nonlinear Anal. 2012, 39, 221–242. [Google Scholar]
  21. Ferreira, R.A.C. Positive solutions for a class of boundary value problems with fractional q-differences. Comput. Math. Appl. 2011, 61, 367–373. [Google Scholar] [CrossRef]
  22. Ferreira, R.A.C. Nontrivial solutions for fractional q-difference boundary value problems. Electron. J. Qual. 2010, 70, 1–10. [Google Scholar]
  23. Li, X.; Han, Z.; Li, X. Boundary value problems of fractional q-difference Schringer equations. J. Appl. Math. Lett. 2015, 46, 100–105. [Google Scholar] [CrossRef]
  24. Cui, Y.Q.; Kang, S.G.; Chen, H.Q. Uniqueness of Solutions for an Integral Boundary Value Problems with Fractional q-differences. J. Appl. Math. Comput. 2018, 8, 524–531. [Google Scholar]
  25. Guo, D. Nonlinear Functional Analysis, 2nd ed.; Science & Technology: Jinan, China, 2000. [Google Scholar]
  26. Nonnenmacher, T.F.; Metzler, R. On the Riemann–Liouville fractional calculus and some recent applications. Fractals 1995, 3, 557–566. [Google Scholar] [CrossRef]
  27. Feng, W.; Zhang, G. New Fixed Point Theorems on Order Intervals and Their Applications. Fixed Point Theory Appl. 2015, 2015, 218. [Google Scholar] [CrossRef] [Green Version]
  28. Jackson, F.H. q-Difference equations. Am. J. Math. 1910, 32, 305–314. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kang, S.; Zhang, Y.; Chen, H.; Feng, W. Positive Solutions for a Class of Integral Boundary Value Problem of Fractional q-Difference Equations. Symmetry 2022, 14, 2465. https://doi.org/10.3390/sym14112465

AMA Style

Kang S, Zhang Y, Chen H, Feng W. Positive Solutions for a Class of Integral Boundary Value Problem of Fractional q-Difference Equations. Symmetry. 2022; 14(11):2465. https://doi.org/10.3390/sym14112465

Chicago/Turabian Style

Kang, Shugui, Yunfang Zhang, Huiqin Chen, and Wenying Feng. 2022. "Positive Solutions for a Class of Integral Boundary Value Problem of Fractional q-Difference Equations" Symmetry 14, no. 11: 2465. https://doi.org/10.3390/sym14112465

APA Style

Kang, S., Zhang, Y., Chen, H., & Feng, W. (2022). Positive Solutions for a Class of Integral Boundary Value Problem of Fractional q-Difference Equations. Symmetry, 14(11), 2465. https://doi.org/10.3390/sym14112465

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