Skip to main content

Existence and multiplicity of positive solutions for nonhomogeneous boundary value problems with fractional q-derivatives

Abstract

In this paper, we study a class of fractional q-difference equations with nonhomogeneous boundary conditions. By applying the classical tools from functional analysis, sufficient conditions for the existence of single and multiple positive solutions to the boundary value problem are obtained in term of the explicit intervals for the nonhomogeneous term. In addition, some examples to illustrate our results are given.

MSC:34A08, 34B18, 39A13.

1 Introduction

Fractional differential equations have attracted considerable interest because of its demonstrated applications in various fields of science and engineering including fluid flow, rheology, diffusive transport akin to diffusion, electrical networks, probability [1, 2]. Many researchers have studied the existence of solutions (or positive solutions) to fractional boundary value problems; for example, see [310] and the references therein.

The early work on q-difference calculus or quantum calculus dates back to Jackson’s papers [11], basic definitions and properties of quantum calculus can be found in the book [12]. For some recent existence results on q-difference equations, we refer to [1315] and the references therein.

The fractional q-difference calculus had its origin in the works by Al-Salam [16] and Agarwal [17]. More recently, there seems to be new interest in the study of this subject and many new developments were made in this theory of fractional q-difference calculus [1822]. Specifically, fractional q-difference equations have attracted the attentions of several researchers. Some recent work on the existence theory of fractional q-difference equations can be found in [20, 2331]. However, the study of boundary value problems for nonlinear fractional q-difference equations is still in the initial stage and many aspects of this topic need to be explored.

By using a fixed-point theorem in a cone, M. El-Shahed and F. Al-Askar [25] were concerned with the existence of positive solutions to nonlinear q-difference equation:

{ D q α C u ( t ) + a ( t ) f ( u ( t ) ) = 0 , 0 < t < 1 , 2 < α 3 , u ( 0 ) = D q 2 u ( 0 ) = 0 , a D q u ( 1 ) + b D q 2 u ( 1 ) = 0 ,

where a,b0 and D q α C is the fractional q-derivatives of the Caputo type.

In [27], Graef and Kong investigated the boundary value problem with fractional q-derivatives

{ ( D q α u ) ( t ) + f ( t , u ( t ) ) = 0 , 0 < t < 1 , n 1 < α n , n N , ( D q i u ) ( 0 ) = 0 , i = 0 , , n 2 , b D q u ( 1 ) = j = 1 m a j D q u ( t j ) + λ ,

where λ0 is a parameter, and the uniqueness, existence and nonexistence of positive solutions are considered in terms of different ranges of λ.

By applying the Banach contraction principle, Krasnoselskii’s fixed-point theorem, and the Leray-Schauder nonlinear alternative, Ahmad, Ntouyas and Purnaras [29] studied the existence of solution for the following nonlinear fractional q-difference equation with nonlocal boundary conditions:

{ ( C D q α u ) ( t ) = f ( t , u ( t ) ) , 0 t 1 , 1 < α 2 , a 1 u ( 0 ) b 1 D q u ( 0 ) = c 1 u ( η 1 ) , a 2 u ( 1 ) + b 2 D q u ( 1 ) = c 2 u ( η 2 ) ,

where D q α C is the fractional q-derivative of the Caputo type, and a i , b i , c i , η i R.

Recently, in [32], the authors investigate the following singular semipositone integral boundary value problem for fractional q-derivatives equation:

{ ( D q α u ) ( t ) + f ( t , u ( t ) ) = 0 , t ( 0 , 1 ) , 2 < α 3 , u ( 0 ) = ( D q u ) ( 0 ) = 0 , u ( 1 ) = μ 0 1 u ( s ) d q s ,

where 0<μ< [ α ] q , D q α is the q-derivative of Riemann-Liouville type of order α, f:[0,1]×(0,+)(,+) is continuous and semipositone, and may be singular at u=0.

Since finding positive solutions of boundary value problems is interest in various fields of sciences, fractional q-calculus equations has tremendous potential for applications. In this paper, we will deal with the following nonhomogeneous boundary value problem with fractional q-derivatives:

{ ( D q α u ) ( t ) + f ( t , u ( t ) ) = 0 , t ( 0 , 1 ) , u ( 0 ) = ( D q u ) ( 0 ) = 0 , γ ( D q u ) ( 1 ) + β ( D q 2 u ) ( 1 ) = λ ,
(1.1)

where q(0,1), 2<α3, γ0, β>0, and λ is a parameter, D q α is the q-derivative of Riemann-Liouville type of order α, f:[0,1]×RR is continuous. In the present work, we gave the corresponding Green’s function of the boundary value problem (1.1) and its properties. By using the generalized Banach contraction principle and Krasnoselskii’s fixed-point theorem, the uniqueness, existence, and multiplicity of positive solution to the BVP (1.1) are obtained in term of the explicit intervals for the nonhomogeneous term. Our results are different from those of [25, 27].

2 Preliminaries on q-calculus and lemmas

For the convenience of the reader, below we cite some definitions and fundamental results on q-calculus as well as the fractional q-calculus. The presentation here can be found in, for example, [12, 18, 20, 22].

Let q(0,1) and define

[ a ] q = 1 q a 1 q ,aR.

The q-analogue of the power function ( a b ) n with n N 0 :={0,1,2,} is

( a b ) ( 0 ) =1, ( a b ) ( n ) = k = 0 n 1 ( a b q k ) ,nN,a,bR.

More generally, if αR, then

( a b ) ( α ) = a α k = 0 a b q k a b q α + k ,a0.
(2.1)

Clearly, if b=0, then a ( α ) = a α . The q-gamma function is defined by

Γ q (x)= ( 1 q ) ( x 1 ) ( 1 q ) x 1 ,xR{0,1,2,},

and satisfies Γ q (x+1)= [ x ] q Γ q (x).

The q-derivative of a function f is defined by

( D q f)(x)= f ( q x ) f ( x ) ( q 1 ) x ,( D q f)(0)= lim x 0 ( D q f)(x),

and q-derivatives of higher order by

( D q 0 f ) (x)=f(x), ( D q n f ) (x)= D q ( D q n 1 f ) (x),nN.

The q-integral of a function f defined in the interval [0,b] is given by

( I q f)(x)= 0 x f(s) d q s=x(1q) k = 0 f ( x q k ) q k ,x[0,b].

If a[0,b] and f is defined in the interval [0,b], then its integral from a to b is defined by

a b f(s) d q s= 0 b f(s) d q s 0 a f(s) d q s.

Similar to that for derivatives, an operator I q n is given by

( I q 0 f ) (x)=f(x), ( I q n f ) (x)= I q ( I q n 1 f ) (x),nN.

The fundamental theorem of calculus applies to these operators I q and D q , i.e.,

( D q I q f)(x)=f(x),

and if f is continuous at x=0, then

( I q D q f)(x)=f(x)f(0).
(2.2)

The following formulas will be used later, namely, the integration by parts formula:

0 x f(s)( D q g)(s) d q s= [ f ( s ) g ( s ) ] s = 0 s = x 0 x ( D q f)(s)g(qs) d q s,

and

(2.3)
(2.4)
(2.5)
(2.6)

where D q t denotes the derivative with respect to the variable t.

Definition 2.1 Let α0 and f be a function defined on [0,1]. The fractional q-integral of Riemann-Liouville type is ( I q 0 f)(x)=f(x) and

( I q α f ) (x)= 1 Γ q ( α ) 0 x ( x q s ) ( α 1 ) f(s) d q s,α>0,x[0,1].

Definition 2.2 The fractional q-derivative of the Riemann-Liouville type of order α0 is defined by ( D q 0 f)(x)=f(x) and

( D q α f ) (x)= ( D q [ α ] I q [ α ] α f ) (x),α>0,

where [α] is the smallest integer greater than or equal to α.

Lemma 2.3 ([20])

Assume that α0 and abt, then ( t a ) ( α ) ( t b ) ( α ) .

Lemma 2.4 Let α,β0 and f be a function defined on [0,1]. Then the following formulas hold:

  1. (1)

    ( I q β I q α f)(x)=( I q α + β f)(x),

  2. (2)

    ( D q α I q α f)(x)=f(x).

Lemma 2.5 ([20])

Let α>0 and n be a positive integer. Then the following equality holds:

( I q α D q n f ) (x)= ( D q n I q α f ) (x) k = 0 n 1 x α n + k Γ q ( α + k n + 1 ) ( D q k f ) (0).

Lemma 2.6 ([22])

Let α R + , λ(1,+), the following is valid:

I q α ( ( t a ) ( λ ) ) = Γ q ( λ + 1 ) Γ q ( α + λ + 1 ) ( t a ) ( α + λ ) ,0<a<t<b.

Particularly, for λ=0, a=0, using q-integration by parts, we have

( I q α 1 ) ( t ) = 1 Γ q ( α ) 0 t ( t q s ) ( α 1 ) d q s = 1 Γ q ( α ) 0 t D q s ( ( t s ) ( α ) ) [ α ] q d q s = 1 Γ q ( α + 1 ) 0 t s D q ( ( t s ) ( α ) ) d q s = 1 Γ q ( α + 1 ) t ( α ) .

Obviously, we have 0 t ( t q s ) ( α 1 ) d q s= 1 [ α ] q t ( α ) , and

0 t ( 1 q s ) ( α 1 ) d q s = 0 t D q s ( ( 1 s ) ( α ) ) [ α ] q d q s = 1 [ α ] q 0 t s D q ( ( 1 s ) ( α ) ) d q s = 1 [ α ] q [ 1 ( 1 t ) ( α ) ] .

In order to define the solution for the problem (1.1), we need the following lemmas.

Lemma 2.7 For given yC[0,1], the unique solution of the boundary value problem

( D q α u ) (t)+y(t)=0,t(0,1),2<α3,
(2.7)

subject to the boundary conditions

u(0)=( D q u)(0)=0,γ( D q u)(1)+β ( D q 2 u ) (1)=λ,
(2.8)

is given by

u(t)= 0 1 G(t,qs)y(s) d q s+ λ t α 1 ( γ + [ α 2 ] q β ) [ α 1 ] q ,
(2.9)

where

G(t,s)={ γ t α 1 ( 1 s ) ( α 2 ) + [ α 2 ] q β t α 1 ( 1 s ) ( α 3 ) ( γ + [ α 2 ] q β ) Γ q ( α ) ( t s ) ( α 1 ) Γ q ( α ) , 0 s t 1 , γ t α 1 ( 1 s ) ( α 2 ) + [ α 2 ] q β t α 1 ( 1 s ) ( α 3 ) ( γ + [ α 2 ] q β ) Γ q ( α ) , 0 t s 1 .
(2.10)

Proof Since 2<α3, we put n=3. In view of Definition 2.1 and Lemma 2.4, we see that

( D q α u ) (t)=y(t) ( I q α D q 3 I q 3 α u ) (t)= ( I q α y ) (t).

Then it follows from Lemma 2.5 that the solution u(t) of (2.7) and (2.8) is given by

u(t)= c 1 t α 1 + c 2 t α 2 + c 3 t α 3 0 t ( t q s ) ( α 1 ) Γ q ( α ) y(s) d q s,
(2.11)

for some constants c 1 , c 2 , c 3 R. From u(0)=0, we have c 3 =0.

Differentiating both sides of (2.11) and with the help of (2.4) and (2.6), we obtain,

( D q u)(t)= [ α 1 ] q c 1 t α 2 + [ α 2 ] q c 2 t α 3 0 t [ α 1 ] q ( t q s ) ( α 2 ) Γ q ( α ) y(s) d q s,

and

( D q 2 u ) ( t ) = [ α 1 ] q [ α 2 ] q c 1 t α 3 + [ α 2 ] q [ α 3 ] q c 2 t α 4 0 t [ α 1 ] q [ α 2 ] q ( t q s ) ( α 3 ) Γ q ( α ) y ( s ) d q s .

Then by the boundary condition ( D q u)(0)=0, we get c 2 =0. Using the boundary condition γ( D q u)(1)+β( D q 2 u)(1)=λ, we get

c 1 = 1 γ + [ α 2 ] q β × ( γ 0 1 ( 1 q s ) ( α 2 ) Γ q ( α ) y ( s ) d q s + β 0 1 [ α 2 ] q ( 1 q s ) ( α 3 ) Γ q ( α ) y ( s ) d q s + λ [ α 1 ] q ) .

Hence, we have

u ( t ) = t α 1 γ + [ α 2 ] q β ( γ 0 1 ( 1 q s ) ( α 2 ) Γ q ( α ) y ( s ) d q s + β 0 1 [ α 2 ] q ( 1 q s ) ( α 3 ) Γ q ( α ) y ( s ) d q s + λ [ α 1 ] q ) 0 t ( t q s ) ( α 1 ) Γ q ( α ) y ( s ) d q s = 0 1 G ( t , q s ) y ( s ) d q s + λ t α 1 ( γ + [ α 2 ] q β ) [ α 1 ] q .

This completes the proof of the lemma. □

Lemma 2.8 The function G(t,s) defined by (2.10) satisfies the following conditions:

  1. (i)

    G(t,qs)0, and G(t,qs)G(1,qs) for all 0t,s1.

  2. (ii)

    G(t,qs) t α 1 G(1,qs) for all 0t,s1.

Proof We start by defining the following two functions:

g 1 ( t , s ) = γ t α 1 ( 1 s ) ( α 2 ) + [ α 2 ] q β t α 1 ( 1 s ) ( α 3 ) ( γ + [ α 2 ] q β ) Γ q ( α ) , 0 t s 1 , g 2 ( t , s ) = g 1 ( t , s ) ( t s ) ( α 1 ) Γ q ( α ) , 0 s t 1 .

Obviously, g 1 (t,qs)0. Now g 2 (0,qs)=0, and for t0

g 2 ( t , q s ) = γ t α 1 ( 1 q s ) ( α 2 ) + [ α 2 ] q β t α 1 ( 1 q s ) ( α 3 ) ( γ + [ α 2 ] q β ) Γ q ( α ) ( t q s ) ( α 1 ) Γ q ( α ) = t α 1 ( γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ( γ + [ α 2 ] q β ) ( 1 q s t ) ( α 1 ) ) ( γ + [ α 2 ] q β ) Γ q ( α ) t α 1 ( γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ( γ + [ α 2 ] q β ) ( 1 q s ) ( α 1 ) ) ( γ + [ α 2 ] q β ) Γ q ( α ) = t α 1 ( γ ( ( 1 q s ) ( α 2 ) ( 1 q s ) ( α 1 ) ) + [ α 2 ] q β ( ( 1 q s ) ( α 3 ) ( 1 q s ) ( α 1 ) ) ) ( γ + [ α 2 ] q β ) Γ q ( α ) 0 .

Therefore, G(t,qs)0.

Moreover, for s(0,1], it follows from (2.4) and Lemma 2.3 that

D q t g 2 ( t , q s ) = [ α 1 ] q ( γ + [ α 2 ] q β ) Γ q ( α ) ( γ t α 2 ( 1 q s ) ( α 2 ) + [ α 2 ] q β t α 2 ( 1 q s ) ( α 3 ) ( γ + [ α 2 ] q β ) ( t q s ) ( α 2 ) ) [ α 1 ] q t α 2 ( γ + [ α 2 ] q β ) Γ q ( α ) ( γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ( γ + [ α 2 ] q β ) ( 1 q s ) ( α 2 ) ) = [ α 1 ] q [ α 2 ] q β t α 2 ( γ + [ α 2 ] q β ) Γ q ( α ) ( ( 1 q s ) ( α 3 ) ( 1 q s ) ( α 2 ) ) 0 ,

which implies that g 2 (t,qs) is an increasing function with respect to t. It is clear that g 1 (t,qs) is increasing in t. Therefore, G(t,qs) is an increasing function of t for all s(0,1], and so G(t,qs)G(1,qs).

When 0tqs1, then

G(t,qs)= γ t α 1 ( 1 q s ) ( α 2 ) + [ α 2 ] q β t α 1 ( 1 q s ) ( α 3 ) ( γ + [ α 2 ] q β ) Γ q ( α ) G(qs,qs)G(1,qs).

Finally, we prove part (ii). When 0qst1, we have

G ( t , q s ) G ( 1 , q s ) = γ t α 1 ( 1 q s ) ( α 2 ) + [ α 2 ] q β t α 1 ( 1 q s ) ( α 3 ) ( γ + [ α 2 ] q β ) ( t q s ) ( α 1 ) γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ( γ + [ α 2 ] q β ) ( 1 q s ) ( α 1 ) γ t α 1 ( 1 q s ) ( α 2 ) + [ α 2 ] q β t α 1 ( 1 q s ) ( α 3 ) ( γ + [ α 2 ] q β ) t α 1 ( 1 q s ) ( α 1 ) γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ( γ + [ α 2 ] q β ) ( 1 q s ) ( α 1 ) = t α 1 .

If 0tqs1, then we have

G ( t , q s ) G ( 1 , q s ) = γ t α 1 ( 1 q s ) ( α 2 ) + [ α 2 ] q β t α 1 ( 1 q s ) ( α 3 ) γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ( γ + [ α 2 ] q β ) ( 1 q s ) ( α 1 ) t α 1 [ γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ( γ + [ α 2 ] q β ) ( 1 q s ) ( α 1 ) ] γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ( γ + [ α 2 ] q β ) ( 1 q s ) ( α 1 ) = t α 1 ,

which implies that part (ii) holds. This completes the proof of the lemma. □

Remark 2.9 If we let 0<τ<1, then

min t [ τ , 1 ] G(t,qs) τ α 1 G(1,qs),for s[0,1].

According to [20], we may take τ= q n , nN.

3 The main results

Let X=C([0,1]) be a Banach space endowed with the norm u X = max 0 t 1 |u(t)|. Define the cone PX by P={uX:u(t)0,0t1}.

Define the operator T:PX as follows:

(Tu)(t)= 0 1 G(t,qs)f ( s , u ( s ) ) d q s+ λ t α 1 ( γ + [ α 2 ] q β ) [ α 1 ] q .
(3.1)

Theorem 3.1 Assume that f:[0,1]×[0,+)[0,+) is continuous and there exists a nonnegative function hC[0,1] such that

| f ( t , u ) f ( t , v ) | h(t)|uv|,t[0,1],u,v[0,+).
(3.2)

Then the BVP (1.1) has a unique positive solution for any λ(0,+), provided

0 1 s α 1 ( γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ) h(s) d q s< ( γ + [ α 2 ] q β ) Γ q ( α ) 2 .
(3.3)

If, in addition, f(t,0)0 on [0,1], then the conclusion is true for λ=0.

Proof We will show that under the assumptions (3.2) and (3.3), T m is a contraction operator for m sufficiently large.

By the definition of G(t,qs), for u,vP, we have

| ( T u ) ( t ) ( T v ) ( t ) | 0 1 G ( t , q s ) | f ( s , u ( s ) ) f ( s , v ( s ) ) | d q s 0 1 t α 1 ( γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ) ( γ + [ α 2 ] q β ) Γ q ( α ) h ( s ) d q s u v X = t α 1 u v X ( γ + [ α 2 ] q β ) Γ q ( α ) 0 1 ( γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ) h ( s ) d q s = Λ 1 t α 1 ( γ + [ α 2 ] q β ) Γ q ( α ) u v X ,

where Λ 1 = 0 1 (γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) )h(s) d q s.

Consequently,

| ( T 2 u ) ( t ) ( T 2 v ) ( t ) | 0 1 G ( t , q s ) | f ( s , ( T u ) ( s ) ) f ( s , ( T v ) ( s ) ) | d q s Λ 1 u v X ( γ + [ α 2 ] q β ) Γ q ( α ) 0 1 G ( t , q s ) s α 1 h ( s ) d q s Λ 1 t α 1 u v X ( γ + [ α 2 ] q β ) 2 [ Γ q ( α ) ] 2 0 1 s α 1 ( γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ) h ( s ) d q s = Λ 1 Λ 2 t α 1 ( γ + [ α 2 ] q β ) 2 [ Γ q ( α ) ] 2 u v X ,

where Λ 2 = 0 1 s α 1 (γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) )h(s) d q s.

By introduction, we get

| ( T m u ) ( t ) ( T m v ) ( t ) | Λ 1 Λ 2 m 1 t α 1 ( γ + [ α 2 ] q β ) m [ Γ q ( α ) ] m u v X .

From the condition (3.3), we have

Λ 1 Λ 2 m 1 ( γ + [ α 2 ] q β ) m [ Γ q ( α ) ] m = Λ 1 Λ 2 [ Λ 2 ( γ + [ α 2 ] q β ) Γ q ( α ) ] m Λ 1 Λ 2 ( 1 2 ) m < 1 4 ,

for m sufficiently large. So, we get

( T m u ) ( t ) ( T m v ) ( t ) X < 1 4 u v X .

Hence, it follows from the generalized Banach contraction principle that the BVP (1.1) has a unique positive solution for any λ(0,+). If λ=0, then the condition f(t,0)0 on [0,1] and Lemma 2.8 imply that u(t)>0 in (0,1]. This completes the proof of the theorem. □

Remark 3.2 When h(t) h 0 is a constant, the condition (3.2) reduces to a Lipschitz condition.

Our next existence results is based on Krasnoselskii’s fixed-point theorem [33].

Lemma 3.3 (Krasnoselskii’s)

Let E be a Banach space, and let PE be a cone. Assume Ω 1 , Ω 2 are open subsets of E with θ Ω 1 , Ω ¯ 1 Ω 2 and let T:P( Ω ¯ 2 Ω 1 )P be a completely continuous operator such that, either

  1. (1)

    Tuu, uP Ω 1 and Tuu, uP Ω 2 , or

  2. (2)

    Tuu, uP Ω 1 and Tuu, uP Ω 2 .

Then T has at least one fixed point in P( Ω ¯ 2 Ω 1 ).

Define a cone KX by

K= { u X : u ( t ) 0 , u ( t ) t α 1 u , t [ 0 , 1 ] } .

Obviously, K is a cone of nonnegative functions in X.

Lemma 3.4 The operator T:KK is completely continuous.

Proof Firstly, we prove that T(K)K. By (2.9) and Lemma 2.8, we have

T u = max 0 t 1 { 0 1 G ( t , q s ) f ( s , u ( s ) ) d q s + λ t α 1 ( γ + [ α 2 ] q β ) [ α 1 ] q } 0 1 G ( 1 , q s ) f ( s , u ( s ) ) d q s + λ ( γ + [ α 2 ] q β ) [ α 1 ] q .

On the other hand,

(Tu)(t) t α 1 ( 0 1 G ( 1 , q s ) f ( s , u ( s ) ) d q s + λ ( γ + [ α 2 ] q β ) [ α 1 ] q ) t α 1 Tu.

Hence, we have T(K)K.

Next, we show that T is uniformly bounded. For fixed r>0, consider a bounded subset K r of K defined by K r ={uK:ur,r>0}, and let M= max 0 u r |f(t,u)|+1. Then for u K r , we get

| ( T u ) ( t ) | 0 1 G ( t , q s ) | f ( s , u ( s ) ) | d q s + | λ | t α 1 ( γ + [ α 2 ] q β ) [ α 1 ] q M 0 1 G ( 1 , q s ) d q s + | λ | ( γ + [ α 2 ] q β ) [ α 1 ] q < + ,

which implies that T( K r ) is bounded.

Finally, we show that T is equicontinuous. For all ε>0, setting

δ=min { ε ω ( α 1 ) , 1 2 ( ε ω ) 1 α 1 } ,

where

ω= M ( γ + [ α 2 ] q β ) + | λ | Γ q ( α ) ( γ + [ α 2 ] q β ) Γ q ( α ) [ α 1 ] q .

For any u K r , we can prove that if t 1 , t 2 [0,1] and 0< t 2 t 1 <δ, then

| ( T u ) ( t 2 ) ( T u ) ( t 2 ) | <ε.

In fact, we have

| ( T u ) ( t 2 ) ( T u ) ( t 1 ) | 0 1 | G ( t 2 , q s ) G ( t 1 , q s ) | f ( s , u ( s ) ) d q s + | λ | ( t 2 α 1 t 1 α 1 ) ( γ + [ α 2 ] q β ) [ α 1 ] q M 0 1 | G ( t 2 , q s ) G ( t 1 , q s ) | d q s + | λ | ( t 2 α 1 t 1 α 1 ) ( γ + [ α 2 ] q β ) [ α 1 ] q M ( 0 t 1 | G ( t 2 , q s ) G ( t 1 , q s ) | d q s + t 1 t 2 | G ( t 2 , q s ) G ( t 1 , q s ) | d q s + t 2 1 | G ( t 2 , q s ) G ( t 1 , q s ) | d q s ) + | λ | ( t 2 α 1 t 1 α 1 ) ( γ + [ α 2 ] q β ) [ α 1 ] q M γ + [ α 2 ] q β ( 0 t 1 γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) Γ q ( α ) ( t 2 α 1 t 1 α 1 ) d q s + t 1 t 2 γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) Γ q ( α ) ( t 2 α 1 t 1 α 1 ) d q s + t 2 1 γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) Γ q ( α ) ( t 2 α 1 t 1 α 1 ) d q s ) + | λ | ( t 2 α 1 t 1 α 1 ) ( γ + [ α 2 ] q β ) [ α 1 ] q = M ( t 2 α 1 t 1 α 1 ) ( γ + [ α 2 ] q β ) Γ q ( α ) 0 1 ( γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ) d q s + | λ | ( γ + [ α 2 ] q β ) [ α 1 ] q ( t 2 α 1 t 1 α 1 ) = M ( γ + [ α 2 ] q β ) + | λ | Γ q ( α ) ( γ + [ α 2 ] q β ) Γ q ( α ) [ α 1 ] q ( t 2 α 1 t 1 α 1 ) = ω ( t 2 α 1 t 1 α 1 ) .

If δ t 1 < t 2 <1, then

| ( T u ) ( t 2 ) ( T u ) ( t 1 ) | ω ( t 2 α 1 t 1 α 1 ) <ω(α1)( t 2 t 1 )<ω(α1)δε.

If 0 t 1 <δ, t 2 <2δ, then

| ( T u ) ( t 2 ) ( T u ) ( t 1 ) | ω ( t 2 α 1 t 1 α 1 ) <ω t 2 α 1 <ω ( 2 δ ) α 1 ε.

By means of Arzela-Ascoli theorem, T:KK is completely continuous.

For the sake of convenience, we introduce the following weight functions:

ϕ ( r ) = max { f ( t , u ( t ) ) : ( t , u ) [ 0 , 1 ] × [ 0 , r ] } , φ ( r ) = min { f ( t , u ( t ) ) : ( t , u ) [ τ , 1 ] × [ τ α 1 r , r ] } ,

and set

l= ( 0 1 G ( 1 , q s ) d q s ) 1 ,L= ( τ α 1 τ 1 G ( 1 , q s ) d q s ) 1 .

 □

Theorem 3.5 Suppose that there exists two positive numbers ξ 1 < ξ 2 such that one of the following conditions is satisfied

( H 1 ) φ( ξ 1 ) ξ 1 L, ϕ( ξ 2 ) l 2 ξ 2 ;

( H 2 ) ϕ( ξ 1 ) l 2 ξ 1 , φ( ξ 2 ) ξ 2 L.

Then the BVP (1.1) has at least one positive solution u K, such that ξ 1 u ξ 2 for λ(0, ( γ + [ α 2 ] q β ) [ α 1 ] q 2 ξ 1 ]. If, in addition, f(t,0)0 on [0,1], then the conclusion is true for λ=0.

Proof Because the proofs are similar, we prove only the case ( H 1 ). Denote Ω ξ 1 ={uX:u< ξ 1 }. Then for any uK Ω ξ 1 , we get u= ξ 1 , 0u(t) ξ 1 , 0t1, and τ α 1 ξ 1 min τ t 1 t α 1 uu(t)u= ξ 1 , τt1. By assumption ( H 1 ), we have

f(t,u)φ( ξ 1 ) ξ 1 L,τt1.

In view of (2.9) and Lemma 2.8, we have

T u min τ t 1 { 0 1 G ( t , q s ) f ( s , u ( s ) ) d q s + λ t α 1 ( γ + [ α 2 ] q β ) [ α 1 ] q } min τ t 1 0 1 G ( t , q s ) f ( s , u ( s ) ) d q s τ α 1 τ 1 G ( 1 , q s ) f ( s , u ( s ) ) d q s τ α 1 τ 1 G ( 1 , q s ) d q s ξ 1 L = ξ 1 = u .

On the other hand, define Ω ξ 2 ={uX:u< ξ 2 }. For any t[0,1] and uK Ω ξ 2 , we have u= ξ 2 and 0u(t) ξ 2 , 0t1. Thus,

f(t,u)ϕ( ξ 2 ) l 2 ξ 2 ,for 0t1,0u ξ 2 .

It follows

T u = max 0 t 1 { 0 1 G ( t , q s ) f ( s , u ( s ) ) d q s + λ t α 1 ( γ + [ α 2 ] q β ) [ α 1 ] q } 0 1 G ( 1 , q s ) d q s l 2 ξ 2 + λ ( γ + [ α 2 ] q β ) [ α 1 ] q l 1 l 2 ξ 2 + ξ 1 2 < ξ 2 2 + ξ 2 2 = u .

By Lemma 3.3, the operator T has at least one fixed point u K( Ω ¯ ξ 2 Ω ξ 1 ), and ξ 1 u ξ 2 . Since u (t) t α 1 u ξ 1 t α 1 >0, 0<t<1, then, the solution u is positive for λ>0. As in the proof of Theorem 3.1, u (t) is a positive solution for λ=0. This completes the proof of the theorem. □

Theorem 3.6 Suppose that there exists three positive numbers ξ 1 < ξ 2 < ξ 3 such that one of the following conditions is satisfied

( H 3 ) φ( ξ 1 ) ξ 1 L, ϕ( ξ 2 )< l 2 ξ 2 , φ( ξ 3 ) ξ 3 L;

( H 4 ) ϕ( ξ 1 ) l 2 ξ 1 , φ( ξ 2 )> ξ 2 L, ϕ( ξ 3 ) l 2 ξ 3 .

Then the BVP (1.1) has at least two positive solutions u 1 , u 2 K such that ξ 1 u 1 < ξ 2 u 2 ξ 3 for λ(0, ( γ + [ α 2 ] q β ) [ α 1 ] q 2 ξ 1 ]. If, in addition, f(t,0)0 on [0,1], then the conclusion is true for λ=0.

Proof We prove only the case ( H 4 ). Since φ:[0,+)[0,+) is continuous and φ( ξ 2 )> ξ 2 L, there exist two positive numbers η 1 , η 2 such that ξ 1 < η 1 < ξ 2 < η 2 < ξ 3 and φ( η 1 ) η 1 L, φ( η 2 ) η 2 L. Thus, it follows from the assumption ( H 4 ) that

ϕ( ξ 1 ) l 2 ξ 1 ,φ( η 1 ) η 1 L,andφ( η 2 ) η 2 L,ϕ( ξ 3 ) l 2 ξ 3 .

From Theorem 3.5, the operator T has two fixed point u 1 K( Ω ¯ η 1 Ω ξ 1 ), u 2 K( Ω ¯ ξ 3 Ω η 2 ) with ξ 1 u 1 < ξ 2 < u 2 ξ 3 . Therefore, the BVP (1.1) has at least two positive solutions for λ(0, ( γ + [ α 2 ] q β ) [ α 1 ] q 2 ξ 1 ]. As in the proof of Theorem 3.1, u 1 , u 2 are two positive solutions for λ=0. This completes the proof of the theorem. □

Denote the integer part of m by [m]. Generally, we have the following theorem.

Theorem 3.7 Suppose that there exists m+1 positive numbers ξ 1 < ξ 2 << ξ m + 1 such that one of the following conditions is satisfied:

( H 5 ) φ( ξ 2 j 1 )> ξ 2 j 1 L, ϕ( ξ 2 j )< l 2 ξ 2 j , j=1,2,,[ m + 2 2 ];

( H 6 ) ϕ( ξ 2 j 1 )< l 2 ξ 2 j 1 , φ( ξ 2 j )> ξ 2 j L, j=1,2,,[ m + 2 2 ].

Then the BVP (1.1) has at least m positive solutions u i K, i=1,2,,m, such that ξ i < u i < ξ i + 1 for λ(0, ( γ + [ α 2 ] q β ) [ α 1 ] q 2 ξ 1 ]. If, in addition, f(t,0)0 on [0,1], then the conclusion is true for λ=0.

4 Examples

Example 4.1 The fractional q-difference boundary value problem

{ D 0.5 2.5 u ( t ) + 2 e t 5 ( 1 + e t ) ( tan 1 u + t 2 + t sin 2 t + 1 ) = 0 , 0 < t < 1 , u ( 0 ) = ( D 0.5 u ) ( 0 ) = 0 , 0.25 ( D 0.5 u ) ( 1 ) + 0.75 ( D 0.5 2 u ) ( 1 ) = λ
(4.1)

has a unique positive solution for any λ(0,+).

Proof In this case, α=2.5, q=0.5, γ=0.25, β=0.75, λ>0. Let

f(t,u)= 2 e t 5 ( 1 + e t ) ( tan 1 u + t 2 + t sin 2 t + 1 ) ,(t,u)[0,1]×(0,+),

and h(t)= 2 e t 5 ( 1 + e t ) . It is easy to prove that

| f ( t , u ) f ( t , v ) | h ( t ) | tan 1 u tan 1 v | h ( t ) | u v | , for  ( t , u ) , ( t , v ) [ 0 , 1 ] × [ 0 , + ) .

A simple computation showed

( γ + [ α 2 ] q β ) Γ q (α)= ( 0.25 + [ 0.5 ] 0.5 × 0.75 ) Γ 0.5 (2.5)0.9168,

and

Λ 2 = 0 1 s α 1 ( γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ) 2 e s 5 ( 1 + e s ) d q s 2 5 0 1 s α 1 ( γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ) d q s 2 5 0 1 ( γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ) d q s 0.3774 ,

which implies that

0 1 s α 1 ( γ ( 1 q s ) ( α 2 ) + [ α 2 ] q β ( 1 q s ) ( α 3 ) ) h(s) d q s< ( γ + [ α 2 ] q β ) Γ q ( α ) 2 .

Obviously, for any m2, we have

Λ 1 Λ 2 m 1 ( γ + [ α 2 ] q β ) m [ Γ q ( α ) ] m 0.3774 0.9168 × 2 m 1 <0.2058< 1 4 .

Thus, Theorem 3.1 implies that the boundary value problem (4.1) has a unique positive solution for any λ(0,+). □

Example 4.2 Consider the following fractional boundary value problem:

{ D 0.5 2.5 u ( t ) + u 2 ( 1 2 cos ( π t π 2 ) + 1 4 ) = 0 , 0 < t < 1 , u ( 0 ) = ( D 0.5 u ) ( 0 ) = 0 , 2 ( D 0.5 u ) ( 1 ) + 7 ( D 0.5 2 u ) ( 1 ) = λ ,
(4.2)

where α=2.5, q=0.5, γ=2, β=7. Choosing n=1, then τ= q n =0.5.

By calculation, we get (γ+ [ α 2 ] q β) Γ q (α)8.1132. By Lemma 2.6, Lemma 2.8 and with the aid of a computer, we obtain that

l= ( 0 1 G ( 1 , 0.5 s ) d 0.5 s ) 1 1.6756,

and

L= ( ( 0.5 ) α 1 0.5 1 G ( 1 , 0.5 s ) d 0.5 s ) 1 0.0657.

Let f(t,u)= u 2 ( 1 2 cos(πt π 2 )+ 1 4 ). Take ξ 1 = 3 4 , ξ 2 =3, then ( γ + [ α 2 ] q β ) [ α 1 ] q 2 ξ 1 2.9578, and f(t,u) satisfies

  1. (i)

    ϕ( 3 4 )=max{ u 2 ( 1 2 cos(πt π 2 )+ 1 4 ):(t,u)[0,1]×[0, 3 4 ]}= 27 64 < l 2 ξ 1 0.6283;

  2. (ii)

    φ(3)=min{ u 2 ( 1 2 cos(πt π 2 )+ 1 4 ):(t,u)[0.5,1]×[1.0608,3]}0.2652>L ξ 2 0.1971.

So, by Theorem 3.5, the problem (4.2) has one positive solution u such that 1 u X 3 for λ(0,2.9578].

References

  1. Podlubny I: Fractional Differential Equations. Academic Press, San Diego; 1999.

    Google Scholar 

  2. Kibas AA, Srivastava HM, Trujillo JJ: Theory and Applications of Fractional Differential Equations. Elsevier, Amsterdam; 2006.

    Google Scholar 

  3. Al-Refai M, Ali Hajji M: Monotone iterative sequences for nonlinear boundary value problems of fractional order. Nonlinear Anal. 2011, 74: 3531-3539. 10.1016/j.na.2011.03.006

    Article  MathSciNet  Google Scholar 

  4. Ahmad B, Sivasundaram S: On four-point nonlocal boundary value problems of nonlinear integro-differential equations of fractional order. Appl. Math. Comput. 2010, 217: 480-487. 10.1016/j.amc.2010.05.080

    Article  MathSciNet  Google Scholar 

  5. Stanék S: The existence of positive solutions of singular fractional boundary value problems. Comput. Math. Appl. 2011, 62: 1379-1388. 10.1016/j.camwa.2011.04.048

    Article  MathSciNet  Google Scholar 

  6. Zhao Y, Chen H, Huang L: Existence of positive solutions for nonlinear fractional functional differential equation. Comput. Math. Appl. 2012, 64: 3456-3467. 10.1016/j.camwa.2012.01.081

    Article  MathSciNet  Google Scholar 

  7. Agarwal RP, O’Regan D, Stanek S: Positive solutions for Dirichlet problems of singular nonlinear fractional differential equations. J. Math. Anal. Appl. 2010, 371: 57-68. 10.1016/j.jmaa.2010.04.034

    Article  MathSciNet  Google Scholar 

  8. Goodrich CS: On discrete sequential fractional boundary value problems. J. Math. Anal. Appl. 2012, 385: 111-124. 10.1016/j.jmaa.2011.06.022

    Article  MathSciNet  Google Scholar 

  9. Zhang S: Positive solutions to singular boundary value problem for nonlinear fractional differential equation. Comput. Math. Appl. 2010, 59: 1300-1309. 10.1016/j.camwa.2009.06.034

    Article  MathSciNet  Google Scholar 

  10. Zhang X, Liu L, Wu Y: Multiple positive solutions of a singular fractional differential equation with negatively perturbed term. Math. Comput. Model. 2012, 55: 1263-1274. 10.1016/j.mcm.2011.10.006

    Article  MathSciNet  Google Scholar 

  11. Jackson FH: q -Difference equations. Am. J. Math. 1970, 32: 305-314.

    Article  Google Scholar 

  12. Kac V, Cheung P: Quantum Calculus. Springer, New York; 2002.

    Book  Google Scholar 

  13. Ernst, T: The history of q-calculus and a new method. UUDM Report 2000:16, Department of Mathematics, Uppsala University, 2000, ISSN:1101-3591

  14. Ahmad B: Boundary-value problems for nonlinear third-order q -difference equations. Electron. J. Differ. Equ. 2011., 2011: Article ID 94

    Google Scholar 

  15. Goodrich CS: Existence of a positive solution to a system of discrete fractional boundary value problems. Appl. Math. Comput. 2011, 217: 4740-4753. 10.1016/j.amc.2010.11.029

    Article  MathSciNet  Google Scholar 

  16. Al-Salam WA: Some fractional q -integrals and q -derivatives. Proc. Edinb. Math. Soc. 1966/1967, 15(2):135-140. 10.1017/S0013091500011469

    Article  MathSciNet  Google Scholar 

  17. Agarwal RP: Certain fractional q -integrals and q -derivatives. Proc. Camb. Philos. Soc. 1969, 66: 365-370. 10.1017/S0305004100045060

    Article  Google Scholar 

  18. Annaby MH, Mansour IS: q-Fractional Calculus and Equations. Springer, Berlin; 2012.

    Book  Google Scholar 

  19. Atici FM, Eloe PW: Fractional q -calculus on a time scale. J. Nonlinear Math. Phys. 2007, 14: 333-344.

    Article  MathSciNet  Google Scholar 

  20. Ferreira RAC: Nontrivial solutions for fractional q -difference boundary value problems. Electron. J. Qual. Theory Differ. Equ. 2010, 70: 1-10.

    Article  Google Scholar 

  21. Goodrich CS: Continuity of solutions to discrete fractional initial value problem. Comput. Math. Appl. 2010, 59: 3489-3499. 10.1016/j.camwa.2010.03.040

    Article  MathSciNet  Google Scholar 

  22. Rajković PM, Marinković SD, Stanković MS: Fractional integrals and derivatives in q -calculus. Appl. Anal. Discrete Math. 2007, 1: 311-323. 10.2298/AADM0701311R

    Article  MathSciNet  Google Scholar 

  23. Atici FM, Eloe PW: Two-point boundary value problems for finite fractional difference equations. J. Differ. Equ. Appl. 2011, 17: 445-456. 10.1080/10236190903029241

    Article  MathSciNet  Google Scholar 

  24. Ferreira RAC: Positive solutions for a class of boundary value problems with fractional q -differences. Comput. Math. Appl. 2011, 61: 367-373. 10.1016/j.camwa.2010.11.012

    Article  MathSciNet  Google Scholar 

  25. El-Shahed M, Al-Askar F: Positive solutions for boundary value problem of nonlinear fractional q -difference equation. ISRN Math. Anal. 2011., 2011: Article ID 385459

    Google Scholar 

  26. Liang S, Zhang J: Existence and uniqueness of positive solutions for three-point boundary value problem with fractional q -differences. J. Appl. Math. Comput. 2012, 40: 277-288. 10.1007/s12190-012-0551-2

    Article  MathSciNet  Google Scholar 

  27. Graef JR, Kong L: Positive solutions for a class of higher order boundary value problems with fractional q -derivatives. Appl. Math. Comput. 2012, 218: 9682-9689. 10.1016/j.amc.2012.03.006

    Article  MathSciNet  Google Scholar 

  28. Ma J, Yang J: Existence of solutions for multi-point boundary value problem of fractional q -difference equation. Electron. J. Qual. Theory Differ. Equ. 2011, 92: 1-10.

    Article  MathSciNet  Google Scholar 

  29. Ahmad B, Ntouyas S, Purnaras I: Existence results for nonlocal boundary value problems of nonlinear fractional q -difference equations. Adv. Differ. Equ. 2012., 2012: Article ID 140

    Google Scholar 

  30. El-Shahed M, Gaber M: Two-dimensional q -differential transformation and its application. Appl. Math. Comput. 2011, 217: 9165-9172. 10.1016/j.amc.2011.03.152

    Article  MathSciNet  Google Scholar 

  31. Zhao Y, Chen H, Zhang Q: Existence results for fractional q -difference equations with nonlocal q -integral boundary conditions. Adv. Differ. Equ. 2013., 2013: Article ID 48

    Google Scholar 

  32. Zhao Y, Ye G, Chen H: Multiple positive solutions of a singular semipositone integral boundary value problem for fractional q -derivatives equation. Abstr. Appl. Anal. 2013., 2013: Article ID 643571

    Google Scholar 

  33. Guo DJ, Lakshmikantham V: Nonlinear Problems in Abstract Cones. Academic Press, San Diego; 1988.

    Google Scholar 

Download references

Acknowledgements

Dedicated to Professor Hari M Srivastava.

The authors are highly grateful for the referees’ careful reading and comments on this paper. The research is supported by the National Natural Science Foundation of China (Grant No. 11271372, 11201138); it is also supported by the Hunan Provincial Natural Science Foundation of China (Grant No. 13JJ3106, 12JJ2004), and the Scientific Research Fund of Hunan Provincial Education Department (Grant No. 12B034).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yulin Zhao.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The authors declare that the study was realized in collaboration with the same responsibility. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Zhao, Y., Chen, H. & Zhang, Q. Existence and multiplicity of positive solutions for nonhomogeneous boundary value problems with fractional q-derivatives. Bound Value Probl 2013, 103 (2013). https://doi.org/10.1186/1687-2770-2013-103

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-103

Keywords