Research article Special Issues

Positive solutions for a system of fractional q-difference equations with generalized p-Laplacian operators

  • In this paper, we consider the existence of positive solutions for a system of fractional q-difference equations with generalized p-Laplacian operators. By using Guo-Krasnosel'skii fixed point theorem, we obtain some existence results of positive solutions for this system with two parameters under some different combinations of superlinearity and sublinearity of the nonlinear terms. In the end, we give two examples to illustrate our main results.

    Citation: Hongyu Li, Liangyu Wang, Yujun Cui. Positive solutions for a system of fractional q-difference equations with generalized p-Laplacian operators[J]. Electronic Research Archive, 2024, 32(2): 1044-1066. doi: 10.3934/era.2024051

    Related Papers:

    [1] Limin Guo, Weihua Wang, Cheng Li, Jingbo Zhao, Dandan Min . Existence results for a class of nonlinear singular $ p $-Laplacian Hadamard fractional differential equations. Electronic Research Archive, 2024, 32(2): 928-944. doi: 10.3934/era.2024045
    [2] Amila Muthunayake, Cac Phan, Ratnasingham Shivaji . An infinite semipositone problem with a reversed S-shaped bifurcation curve. Electronic Research Archive, 2023, 31(2): 1147-1156. doi: 10.3934/era.2023058
    [3] Xiaoju Zhang, Kai Zheng, Yao Lu, Huanhuan Ma . Global existence and long-time behavior of solutions for fully nonlocal Boussinesq equations. Electronic Research Archive, 2023, 31(9): 5406-5424. doi: 10.3934/era.2023274
    [4] Tao Zhang, Tingzhi Cheng . A priori estimates of solutions to nonlinear fractional Laplacian equation. Electronic Research Archive, 2023, 31(2): 1119-1133. doi: 10.3934/era.2023056
    [5] Yuchen Zhu . Blow-up of solutions for a time fractional biharmonic equation with exponentional nonlinear memory. Electronic Research Archive, 2024, 32(11): 5988-6007. doi: 10.3934/era.2024278
    [6] Messoud Efendiev, Vitali Vougalter . Linear and nonlinear non-Fredholm operators and their applications. Electronic Research Archive, 2022, 30(2): 515-534. doi: 10.3934/era.2022027
    [7] Qingcong Song, Xinan Hao . Positive solutions for fractional iterative functional differential equation with a convection term. Electronic Research Archive, 2023, 31(4): 1863-1875. doi: 10.3934/era.2023096
    [8] Ping Yang, Xingyong Zhang . Existence of nontrivial solutions for a poly-Laplacian system involving concave-convex nonlinearities on locally finite graphs. Electronic Research Archive, 2023, 31(12): 7473-7495. doi: 10.3934/era.2023377
    [9] Seda IGRET ARAZ, Mehmet Akif CETIN, Abdon ATANGANA . Existence, uniqueness and numerical solution of stochastic fractional differential equations with integer and non-integer orders. Electronic Research Archive, 2024, 32(2): 733-761. doi: 10.3934/era.2024035
    [10] Peng Mei, Zhan Zhou, Yuming Chen . Homoclinic solutions of discrete $ p $-Laplacian equations containing both advance and retardation. Electronic Research Archive, 2022, 30(6): 2205-2219. doi: 10.3934/era.2022112
  • In this paper, we consider the existence of positive solutions for a system of fractional q-difference equations with generalized p-Laplacian operators. By using Guo-Krasnosel'skii fixed point theorem, we obtain some existence results of positive solutions for this system with two parameters under some different combinations of superlinearity and sublinearity of the nonlinear terms. In the end, we give two examples to illustrate our main results.



    In this paper, we consider the existence of positive solutions for the following system of fractional q-difference equations with generalized p-Laplacian operators:

    {Dγq(ϕ1(Dαqx))(t)=ηf(t,x(t),y(t)),  0<t<1,Dγq(ϕ2(Dαqy))(t)=ζg(t,x(t),y(t)),  0<t<1,x(0)=Dqx(0)=0, Dqx(1)=β>0, Dαqx(0)=0,y(0)=Dqy(0)=0, Dqy(1)=β>0, Dαqy(0)=0, (1.1)

    where 0<q<1, 2<α<3, 0<γ<1, f, g:[0,1]×[0,+)×[0,+)[0,+) are continuous, η>0 and ζ>0 are two parameters, ϕ1 and ϕ2 are generalized p-Laplacian operators; Dγq and Dαq are the fractional q-derivative of the Riemann-Liouville type, Dq is the q-derivative.

    Due to the extensive application of fractional order equations, many scholars have studied the existence of nontrivial solutions of boundary value problems for fractional order differential equations. In recent years, some authors [1,2,3,4,5] have considered the existence of positive solutions for some Riemann-Liouville type, tempered type, Caputo type and Hadamard type fractional order differential equations. The authors [6,7] have considered the existence of nontrivial solution of Hadamard-type singular fractional differential equations. Some authors [8,9,10,11,12] have considered the existence of nontrivial solutions for some Riemann-Liouville type, tempered type, Caputo type and Hadamard type fractional order differential equations with p-Laplacian operator. Some authors [13,14] have considered the eigenvalue problems of fractional differential equations.

    Meanwhile, after Jackson [15] introduced the q-calculus, Al-Salam [16] and Agarwal [17] developed the fractional q-calculus. Many researchers have studied the existence of nontrivial solutions for fractional q-difference equations these years. The commonly used methods include fixed point theorems, lower-upper solution method, monotone iterative technique, and so on. For example, in [18], Ferreira studied the following boundary value problem of fractional q-difference equation:

    {(Dαqy)(x)=f(x,y(x)), 0<x<1,y(0)=(Dqy)(0)=0, (Dqy)(1)=β0, (1.2)

    where 0<q<1, 2<α3, f:[0,1]×[0,)[0,) is continuous; Dαq is the fractional q-derivative of the Riemann-Liouville type, Dq is the q-derivative. The author obtained the existence of positive solutions about the boundary value problem (1.2) by using Guo-Krasnosel'skii fixed point theorem.

    In [19], Zhai and Ren applied iterative algorithm and lower-upper solution method to study the following fractional q-difference equation:

    {(Dαqu)(t)+f(t,u(t))=0, 0<t<1,u(0)=(Dqu)(0)=0, (Dqu)(1)=0, (1.3)

    where q(0,1), α(2,3). Under some conditions, the authors obtained some existence results of positive or negative solutions for the boundary value problem (1.3).

    In [20], Mao et al. used iterative technique to consider the general fractional q-difference equation of the problem (1.3) as followings:

    {(Dαqu)(t)+f(t,u(t),v(t))=0, 0<t<1,u(0)=(Dqu)(0)=0, (Dqu)(1)=0, (1.4)

    where q(0,1), α(2,3), f may be singular at v=0,t=0,1. The existence of a unique positive solution of the problem (1.4) has been proved.

    In [21], Jiang and Zhong studied the following fractional q-difference equation with p-Laplacian operator:

    {Dβq(ϕp(Dαqx)(t))+f(t,x(t),Dρq(t))=0,x(0)=(Dqx)(0)=0, (Dαqx)(0)=0,x(1)=ζIqx(η), (1.5)

    where α(2,3), β,q,η,ρ(0,1), ϕp(s)=|s|p2s is the p-Laplacian operator (p>1). The authors used Banach's contraction principle to prove the existence and uniqueness of nontrivial solution of the problem (1.5), and also used Guo-Krasnosel'skii fixed point theorem to obtain the existence of positive solutions of the problem (1.5).

    In [22], Li et al. considered the following boundary value problem of nonlinear fractional q-difference equation:

    {Dγq(ϕ(Dαqu)(t))+ηf(u(t))=0,  0<t<1,u(0)=Dqu(0)=0, Dqu(1)=β>0, Dαqu(0)=0, (1.6)

    where 0<q<1, 2<α<3, 0<γ<1, ϕ is the generalized p-Laplacian operator; Dγq and Dαq are the fractional q-derivative of the Riemann-Liouville type, Dq is the q-derivative. The authors used the fixed point theorem to prove the existence of positive solutions of the boundary value problem (1.6).

    In [23], Wang et al. investigated the following boundary value problem of fractional q-difference equation with ϕ-Laplacian:

    {Dβq(ϕ(Dαqu(t)))=λf(u(t)),  0<t<1,u(0)=Dqu(0)=Dqu(1)=0, ϕ(Dαqu(0))=Dq(ϕ(Dαqu(1)))=0, (1.7)

    where 0<q<1, 2<α3, 1<β2, λ>0 is a parameter, and Dβq, Dαq are the standard Riemann-Liouville fractional q-derivatives. The existence and nonexistence of positive solutions of the boundary value problem (1.7) was obtained based on Guo-Krasnosel'skii fixed point theorem on cones.

    Currently, many other authors have studied fractional q-difference equations. Some authors [24,25] have considered the existence of multiple positive solutions for some fractional q-difference equations. The authors [26,27,28,29,30] have considered the existence of nontrivial solutions for fractional q-difference equations with various boundary conditions.

    Meanwhile, many authors have studied the existence of positive solutions of systems of some fractional differential equations with various boundary conditions, see [31,32,33,34,35] and the references therein. For example, in [31], Li et al. investigated the following system of fractional differential equations with p-Laplacian operators:

    {Dα10+(φp1(Dβ10+u(t)))=f(t,v(t)),0<t<1,Dα20+(φp2(Dβ20+v(t)))=g(t,u(t)),0<t<1,u(0)=Dβ10+u(0)=0,  Dγ10+u(1)=m2j=1a1jDγ10+u(ηj),v(0)=Dβ20+v(0)=0  Dγ20+v(1)=m2j=1a2jDγ20+v(ηj), (1.8)

    where αi,γi(0,1],βi(1,2],Dαi0+,Dβi0+ and Dγi0+ are the standard Riemann-Liouville derivatives, i=1,2. The authors derived the conditions for the existence of the maximal and minimal solutions, and obtained the existence of extremal solutions of the system (1.8).

    In [32], He and Song considered the following system of fractional differential equations with p-Laplacian operators and two parameters:

    {Dα10+(φp1(Dβ10+u(t)))=ηf(t,v(t)),0<t<1,Dα20+(φp2(Dβ20+v(t)))=ζg(t,u(t)),0<t<1,u(0)=0,u(1)=a1u(ξ1),Dβ10+u(0)=0,Dβ10+u(1)=b1Dβ10+u(η1),v(0)=0,v(1)=a2v(ξ2),Dβ20+v(0)=0,Dβ20+v(1)=b2Dβ20+v(η2), (1.9)

    where αi,βi(1,2], Dαi0+ and Dβi0+ are the standard Riemann-Liouville derivatives, ξi,ηi(0,1),ai,bi[0,1],i=1,2. η and ζ are positive parameters. By using the Banach contraction mapping principle, the authors gave the existence and uniqueness of the solution for the system (1.9).

    In [33], Hao et al. investigated the following system of fractional boundary value problems with p-Laplacian operators and two parameters:

    {Dα10+(φp1(Dβ10+u(t)))=λf(t,u(t),v(t)),t(0,1),Dα20+(φp2(Dβ20+v(t)))=μg(t,u(t),v(t)),t(0,1),u(0)=u(1)=u(0)=u(1)=0,Dβ10+u(0)=0,Dβ10+u(1)=b1Dβ10+p(η1),v(0)=v(1)=v(0)=v(1)=0,Dβ20+v(0)=0,Dβ20+v(1)=b2Dβ20+v(β2), (1.10)

    where αi(1,2], βi(3,4], Dαi0+ and Dβi0+ are the Riemann-Liouville derivatives, φpi(s)=|s|pi2s,pi>1, f,gC([0,1]×[0,+)×[0,+),[0,+)), λ and μ are positive parameters. By means of Guo-Krasnosel'skii fixed point theorem, the authors obtained various existence results of positive solutions of the system (1.10).

    To the best of our knowledge, limited attention has been devoted to the investigation of systems of fractional q-difference equations. Inspired by the above literatures, in this paper, we consider the existence of positive solutions for the system of fractional q-difference equations (1.1) with generalized p-Laplacian operators and two parameters. Under some sublinear and superlinear conditions, we establish some existence results of positive solutions for the system (1.1) by using Guo-Krasnosel'skii fixed point theorem.

    In this section, we firstly introduce Guo-Krasonsel'skill fixed point theorem. Secondly, we give some knowledge about fractional q-calculus. In the end, we give some lemmas that are used to prove the main results.

    Lemma 2.1. ([36]) Let E be a Banach space, PE be a cone. Assume that Ω1E and Ω2E are bounded open sets with θΩ1Ω2, the operator A:P(ˉΩ2Ω1)P is completely continuous. If the following conditions are satisfied:

    (i)Axx,xPΩ1,Axx,xPΩ2,or(ii)Axx,xPΩ1,Axx,xPΩ2,

    then the operator A has at least one fixed point in P(ˉΩ2Ω1).

    In the following, we give some definitions and lemmas about fractional q-calculus. For the detailed knowledge about fractional q-derivative and fractional q-integral, we can refer to [15,16,17,18].

    Define

    [a]q=1qa1q,aR, q(0,1).

    The q-analogue of the power function (ab)(α) with αR is

    (ab)(α)=aαn=0abqnabqα+n, nN.

    The q-gamma function is defined by

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

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

    The q-derivative of a function f is defined by

    (Dqf)(x)=f(x)f(qx)(1q)x,(Dqf)(0)=limx0(Dqf)(x),

    and q-derivative of higher order by

    (D0qf)(x)=f(x),(Dnqf)(x)=Dq(Dn1qf)(x),nN.

    The q-integral of a function f is given by

    (Iqf)(x)=x0f(t)dqt=x(1q)n=0f(xqn)qn,x[0,b],

    where f is defined in the interval [0,b].

    Definition 2.1. ([17,18]) The fractional q-integral of the Riemann-Liouville type is defined by (Iq0f)(x)=f(x), and

    (Iαqf)(x)=1Γq(α)x0(xqt)(α1)f(t)dqt,α>0,x[0,1].

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

    (Dαqf)(x)=(DmqImαqf)(x),α>0,

    where m=[α].

    The paper [37] introduced the definition of a generalized p-Laplacian operator ϕ, which included two important cases ϕ(u)=u and ϕ(u)=|u|p2u(p1). In this paper, we assume that ϕ1 and ϕ2 are generalized p-Laplacian operators, namely, ϕ1 and ϕ2 satisfy the following condition:

    (H0) ϕi:RR(i=1,2) is an odd and increasing homeomorphism, and there exist increasing homeomorphisms ψ1,ψ2,ψ3,ψ4:(0,)(0,) such that

    ψ1(x)ϕ1(y)ϕ1(xy)ψ2(x)ϕ1(y),  x,y>0,
    ψ3(x)ϕ2(y)ϕ2(xy)ψ4(x)ϕ2(y),  x,y>0.

    In the following, we give the other important lemmas. We list the following fractional q-difference equation with homogeneous boundary conditions:

    {Dγq(ϕ(Dαqv))(t)+ηf(v(t)+φ(t))=0, 0<t<1,v(0)=Dqv(0)=0,Dqv(1)=0,Dαqv(0)=0, (2.1)

    where φ(t)=βtα1[α1]q, and φ(t) is the unique solution of the following fractional q-difference equation with nonhomogeneous boundary condition

    {Dαqφ(t)=0, 0<t<1,φ(0)=Dqφ(0)=0,Dqφ(1)=β>0, Dαqφ(0)=0, (2.2)

    where 0<q<1,2<α<3,0<γ<1, ϕ is a generalized p-Laplacian operator.

    Lemma 2.2. ([22]) Let v(t) be a solution of the boundary value problem (2.1). Then u(t)=v(t)+φ(t) is the solution of the boundary value problem (1.6).

    Lemma 2.3. ([22]) Let 2<α<3,0<γ<1,yC[0,1] be a given function. Then the following boundary value problem of fractional q-difference equation

    {Dγq(ϕ(Dαqx))(t)+ηy(t)=0,x(0)=Dqx(0)=0,Dqx(1)=0,Dαqx(0)=0

    has a unique solution

    x(t)=10G(t,qs)ϕ1(ηΓq(γ)s0(sqτ)(γ1)y(τ)dqτ)dqs,

    where

    G(t,qs)=1Γq(α){(1qs)(α2)tα1(tqs)(α1),0qst1,(1qs)(α2)tα1,                   0tqs1. (2.3)

    Lemma 2.4. ([18]) The function G(t,qs) defined by (2.3) has the following properties:

    (1) G(t,qs)0,  G(t,qs)G(1,qs),  0t,s1, (2.4)
    (2) G(t,qs)tα1G(1,qs),  0t,s1. (2.5)

    Lemma 2.5. ([37]) Let (H0) hold. Then we have

    ψ21(x)yϕ11(xϕ1(y))ψ11(x)y,  x,y>0,
    ψ41(x)yϕ21(xϕ2(y))ψ31(x)y,  x,y>0.

    Let E=C[0,1]×C[0,1] with the norm (x,y)E=x+y, where x=maxt[0,1]|x(t)| and y=maxt[0,1]|y(t)|. It is obvious that E is a Banach space.

    Set P={(x,y)E:x(t)0,y(t)0,mint[θ,1](x(t)+y(t))θα1(x,y)E}, where θ is a real constant and 0<θ<1.

    By [22], we define the following operators Aη and Aζ:

    Aη(x,y)(t)=10G(t,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs, t[0,1],
    Aζ(x,y)(t)=10G(t,qs)ϕ21(ζΓq(γ)s0(sqτ)(γ1)g(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs,t[0,1],

    where G(t,qs) is defined by (2.3). Let A(x,y)=(Aη(x,y),Aζ(x,y)),(x,y)E. Then by the literature [22], we easily know that the fixed points of the operator A are solutions of the system of fractional q-difference equations (1.1).

    Lemma 2.6. A:PP is completely continuous.

    Proof. For (x,y)P, we easily have Aη(x,y)(t)0,Aζ(x,y)(t)0,t[0,1].

    By (2.5), for t[θ,1], we have

    Aη(x,y)(t)=10G(t,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs10tα1G(1,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqsθα110G(1,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs=θα1Aη(x,y). (2.6)

    Similar to the proof of (2.6), when t[θ,1], we easily have

    Aζ(x,y)(t)=10G(t,qs)ϕ12(ζΓq(γ)s0(sqτ)(γ1)g(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqstα110G(1,qs)ϕ12(ζΓq(γ)s0(sqτ)(γ1)g(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqsθα1Aζ(x,y). (2.7)

    By (2.6) and (2.7), we get

    mint[θ,1](Aη(x,y)(t)+Aζ(x,y)(t))θα1(Aη(x,y)+Aζ(x,y))=θα1A(x,y)E. (2.8)

    From (2.8), we have A(P)P.

    In the following, we prove A:PP is completely continuous. Firstly, we prove A is bounded.

    Let DP be bounded. Namely, there exists K>0 such that (x,y)E<K,(x,y)D. By the continuity of f and g, we know that there exists M>0 such that

    maxt[0,1],(x,y)D|f(t,x(t)+φ(t),y(t)+φ(t))|<M, (2.9)
    maxt[0,1],(x,y)D|g(t,x(t)+φ(t),y(t)+φ(t))|<M. (2.10)

    By (2.9), (2.10) and Lemma 2.4, for (x,y)D, we get

    |Aη(x,y)(t)|10G(1,qs)ϕ11(ηMΓq(γ)s0(sqτ)(γ1)dqτ)dqs, (2.11)
    |Aζ(x,y)(t)|10G(1,qs)ϕ12(ζMΓq(γ)s0(sqτ)(γ1)dqτ)dqs. (2.12)

    By (2.11) and (2.12), we easily know that A(D) is bounded. Secondly, we prove A is equicontinuous on D. Namely, for each (x,y)D,ε>0,δ>0 such that |t2t1|<δ, we have

    |Aη(x,y)(t2)Aη(x,y)(t1)|<ε,|Aζ(x,y)(t2)Aζ(x,y)(t1)|<ε.

    In fact, assume that 0<t1<t2<1, then we have

    |Aη(x,y)(t2)Aη(x,y)(t1)|=|10(G(t2,qs)G(t1,qs))ϕ11(ηΓq(γ)s0(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs|10|G(t2,qs)G(t1,qs)|ϕ11(ηMΓq(γ)s0(sqτ)(γ1)dqτ)dqst10|(1qs)(α2)(t2α1t1α1)|1Γq(α)ϕ11(ηMΓq(γ)s0(sqτ)(γ1)dqτ)dqs+t2t1|(1qs)(α2)(t2α1t1α1)(t2qs)(α1)|1Γq(α)ϕ11(ηMΓq(γ)s0(sqτ)(γ1)dqτ)dqs+1t2|(1qs)(α2)(t2α1t1α1)|1Γq(α)ϕ11(ηMΓq(γ)s0(sqτ)(γ1)dqτ)dqs. (2.13)

    By (2.13), we have

    |Aη(x,y)(t2)Aη(x,y)(t1)|0(t1t2).

    Similarly, we also have

    |Aζ(x,y)(t2)Aζ(x,y)(t1)|0(t1t2).

    Hence, by Arzela-Ascoli theorem and the continuity of f and g, we have A:PP is completely continuous.

    In the following, we give the denotations that we need in this section.

    Let

    f0=limsupx+y0+maxt[0,1]f(t,x+φ(t),y+φ(t))ϕ1(x+y), g0=limsupx+y0+maxt[0,1]g(t,x+φ(t),y+φ(t))ϕ2(x+y),
    f=liminfx+ymint[θ,1]f(t,x+φ(t),y+φ(t))ϕ1(x+y), g=liminfx+ymint[θ,1]g(t,x+φ(t),y+φ(t))ϕ2(x+y).

    For f0,g0,f,g(0,), we denote that

    D1=ψ2(θ22α2M3)f,  D2=ψ1(12M1)f0,
    D3=ψ4(θ22α2M4)g,  D4=ψ3(12M2)g0.

    Let

    ˉf0=liminfx+y0+mint[θ,1]f(t,x+φ(t),y+φ(t))ϕ1(x+y), ˉg0=liminfx+y0+mint[θ,1]g(t,x+φ(t),y+φ(t))ϕ2(x+y),
    ˉf=limsupx+ymaxt[0,1]f(t,x+φ(t),y+φ(t))ϕ1(x+y), ˉg=limsupx+ymaxt[0,1]g(t,x+φ(t),y+φ(t))ϕ2(x+y).

    For ˉf0,ˉg0,ˉf,ˉg(0,), we give the following denotations:

    Z1=ψ2(θ22α2M3)ˉf0,  Z2=ψ1(12M1)ˉf,
    Z3=ψ4(θ22α2M4)ˉg0,  Z4=ψ3(12M2)ˉg,

    where

    M1=10G(1,qs)ψ11(sγΓq(γ+1))dqs,  M2=10G(1,qs)ψ31(sγΓq(γ+1))dqs,
    M3=1θG(1,qs)ψ21(1Γq(γ)sθ(sqτ)(γ1)dqτ)dqs,
    M4=1θG(1,qs)ψ41(1Γq(γ)sθ(sqτ)(γ1)dqτ)dqs.

    Theorem 3.1. (1) Assume that f0,g0,f,g(0,),D1<D2,D3<D4, then for each η(D1,D2) and ζ(D3,D4), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (2) Assume that f0=0,g0,f,g(0,),D3<D4, then for each η(D1,) and ζ(D3,D4), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (3) Assume that f0,f,g(0,),g0=0,D1<D2, then for each η(D1,D2) and ζ(D3,), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (4) Assume that f0=g0=0,f,g(0,), then for each η(D1,) and ζ(D3,), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (5) Assume that f0,g0(0,),f= or f0,g0(0,),g=, then for each η(0,D2) and ζ(0,D4), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (6) Assume that f0=0,g0(0,),g= or f0=0,g0(0,),f=, then for each η(0,) and ζ(0,D4), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (7) Assume that f0(0,),g0=0,g= or f0(0,),g0=0,f=, then for each η(0,D2) and ζ(0,), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (8) Assume that f0=g0=0,g= or f0=g0=0,f=, then for each η(0,) and ζ(0,), the system of fractional q-difference equations (1.1) has at least one positive solution.

    Proof. In the following, we only prove the Cases (1) and (6).

    Case (1): Since η(D1,D2) and ζ(D3,D4), we easily know that there exists ε>0 such that

    0<ψ2(θ22α2M3)fεηψ1(12M1)f0+ε, (3.1)
    0<ψ4(θ22α2M4)gεζψ3(12M2)g0+ε. (3.2)

    For the above ε>0 in (3.1) and (3.2), there exists r1>0 such that

    f(t,x+φ,y+φ)(f0+ε)ϕ1(x+y), t[0,1], 0x+yr1, (3.3)
    g(t,x+φ,y+φ)(g0+ε)ϕ2(x+y), t[0,1], 0x+yr1. (3.4)

    Let W1={(x,y)E:(x,y)E<r1}. For (x,y)PW1, we have

    0x(t)+y(t)x+y=(x,y)E=r1, t[0,1].

    By Lemmas 2.4, 2.5 and (3.1)–(3.4), we have

    Aη(x,y)(t)=10G(t,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs10G(t,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)(f0+ε)ϕ1(x(τ)+y(τ))dqτ)dqs10G(1,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)(f0+ε)ϕ1(r1)dqτ)dqs10G(1,qs)ψ11(ηΓq(γ)s0(sqτ)(γ1)(f0+ε)dqτ)dqsr1ψ11(η(f0+ε))10G(1,qs)ψ11(sγΓq(γ+1))dqsr1r12=(x,y)E2. (3.5)

    and

    Aζ(x,y)(t)=10G(t,qs)ϕ12(ζΓq(γ)s0(sqτ)(γ1)g(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs10G(t,qs)ϕ12(ζΓq(γ)s0(sqτ)(γ1)(g0+ε)ϕ2(x(τ)+y(τ))dqτ)dqs10G(1,qs)ϕ12(ζ(g0+ε)Γq(γ)s0(sqτ)(γ1)ϕ2(r1)dqτ)dqs10G(1,qs)ψ13(ζ(g0+ε)Γq(γ)s0(sqτ)(γ1)dqτ)dqsr1ψ13(ζ(g0+ε))10G(1,qs)ψ13(sγΓq(γ+1))dqsr1r12=(x,y)E2. (3.6)

    By (3.5) and (3.6), we have

    A(x,y)E=Aη(x,y)+Aζ(x,y)(x,y)E,(x,y)PW1. (3.7)

    For ε>0 in (3.1) and (3.2), from the definitions of f and g, there exists ˉr2>0 such that

    f(t,x+φ,y+φ)(fε)ϕ1(x+y),t[θ,1],x+yˉr2, (3.8)
    g(t,x+φ,y+φ)(gε)ϕ2(x+y),t[θ,1],x+yˉr2. (3.9)

    Take r2=max{2r1,θ1αˉr2}. Let W2={(x,y)E:(x,y)E<r2}. For (x,y)PW2, we have mint[θ,1](x(t)+y(t))θα1(x,y)E=θα1r2ˉr2.

    By (3.8), (3.9) and Lemmas 2.4 and 2.5, we have

    Aη(x,y)Aη(x,y)(θ)=10G(θ,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs1θG(θ,qs)ϕ11(ηΓq(γ)sθ(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs1θG(θ,qs)ϕ11(ηΓq(γ)sθ(sqτ)(γ1)(fε)ϕ1(x(τ)+y(τ))dqτ)dqs1θθα1G(1,qs)ϕ11(ηΓq(γ)sθ(sqτ)(γ1)(fε)ϕ1(θα1(x,y)E)dqτ)dqsθ2α21θG(1,qs)ψ12(ηΓq(γ)(fε)sθ(sqτ)(γ1)dqτ)(x,y)Edqs=θ2α2ψ12(η(fε))1θG(1,qs)ψ12(1Γq(γ)sθ(sqτ)(γ1)dqτ)dqsr2(x,y)E2, (3.10)

    and

    Aζ(x,y)Aζ(x,y)(θ)=10G(θ,qs)ϕ12(ζΓq(γ)s0(sqτ)(γ1)g(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs1θG(θ,qs)ϕ12(ζΓq(γ)sθ(sqτ)(γ1)g(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqsθα11θG(1,qs)ϕ12(ζΓq(γ)sθ(sqτ)(γ1)(gε)ϕ2(x(τ)+y(τ))dqτ)dqsθα11θG(1,qs)ϕ12(ζ(gε)Γq(γ)sθ(sqτ)(γ1)ϕ2(θα1(x,y)E)dqτ)dqsθ2α21θG(1,qs)ψ14(ζ(gε)Γq(γ)sθ(sqτ)(γ1)dqτ)dqs(x,y)E=θ2α2ψ14(ζ(gε))1θG(1,qs)ψ14(1Γq(γ)sθ(sqτ)(γ1)dqτ)dqs(x,y)E=θ2α2ψ14(ζ(gε))M4r212r2=(x,y)E2. (3.11)

    From (3.10) and (3.11), we have

    A(x,y)E=Aη(x,y)+Aζ(x,y)(x,y)E,(x,y)PW2. (3.12)

    By (3.7), (3.12) and Lemma 2.1, we know that A has at least one fixed point (x,y)P(¯W2W1). So the system of fractional q-difference equations (1.1) has at least one positive solution. The proof of the case (1) is completed.

    Case (6): Since η(0,) and ζ(0,D4), we easily know that there exists ε>0 such that

    0<η<ψ1(12M1)1ε,  ψ4(θ22αM4)ε<ζ<ψ3(12M2)1g0+ε. (3.13)

    Since f0=0 and g0(0,), for the above ε>0 in (3.13), we know that there exists r3>0 such that

    f(t,x+φ,y+φ)<εϕ1(x+y),t[0,1],0x+yr3, (3.14)
    g(t,x+φ,y+φ)<(g0+ε)ϕ2(x+y),t[0,1],0x+yr3. (3.15)

    Let W3={(x,y)E:(x,y)E<r3}. By (3.13), (3.14) and Lemma 2.5, for any (x,y)PW3,t[0,1], we have

    Aη(x,y)(t)=10G(t,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs10G(t,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)εϕ1(x(τ)+y(τ))dqτ)dqs10G(1,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)εϕ1(r3)dqτ)dqs10G(1,qs)ψ11(ηεΓq(γ)s0(sqτ)(γ1)dqτ)dqsr3=ψ11(ηε)10G(1,qs)ψ11(sγΓq(γ+1))dqsr3<r32=(x,y)E2. (3.16)

    By (3.16), we have

    Aη(x,y)(x,y)E2,(x,y)PW3. (3.17)

    By (3.13), (3.15) and Lemma 2.5, similar to the proof of (3.16), we easily obtain

    Aζ(x,y)(x,y)E2,(x,y)PW3. (3.18)

    By (3.17) and (3.18), we have

    A(x,y)E=Aη(x,y)+Aζ(x,y)(x,y)E,(x,y)PW3. (3.19)

    Since g=, for ε>0 in (3.13), we know that there exists ˉr4>0 such that

    g(t,x+φ,y+φ)1εϕ2(x,y),t[θ,1],x,y0,x+yˉr4. (3.20)

    Take r4=max{3r3,ˉr4θ1α}. Let W4={(x,y)E:(x,y)E<r4}. For any (x,y)PW4, we can easily know that

    mint[θ,1](x(t)+y(t))θα1(x,y)E=θα1r4ˉr4. (3.21)

    Hence, by (3.20), (3.21) and Lemma 2.5, for any (x,y)PW4, we have

    Aζ(x,y)(θ)=10G(θ,qs)ϕ12(ζΓq(γ)s0(sqτ)(γ1)g(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs1θG(θ,qs)ϕ12(ζΓq(γ)sθ(sqτ)(γ1)g(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqsθα11θG(1,qs)ϕ12(ζΓq(γ)1εsθ(sqτ)(γ1)ϕ2(x(τ)+y(τ))dqτ)dqsθα11θG(1,qs)ϕ12(ζΓq(γ)1εsθ(sqτ)(γ1)ϕ2(θα1(x,y)E)dqτ)dqsθα11θG(1,qs)ψ14(ζεΓq(γ)sθ(sqτ)(γ1)dqτθα1(x,y)E)dqs=θ2α2ψ14(ζε)1θG(1,qs)ψ14(1Γq(γ)sθ(sqτ)(γ1)dqτ)dqsr4=θ2α2ψ14(ζε)M4r4r4=(x,y)E. (3.22)

    By (3.22), we have

    A(x,y)EAζ(x,y)(x,y)E,(x,y)PW4. (3.23)

    Hence, by (3.19), (3.23) and Lemma 2.1, we can obtain that A has at least one fixed point (x,y)P(¯W4W3). So the system of fractional q-difference equations (1.1) has at least one positive solution.

    Theorem 3.2. (1) Assume that ˉf0,ˉg0,ˉf,ˉg(0,), and Z1<Z2,Z3<Z4, then for each η(Z1,Z2) and ζ(Z3,Z4), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (2) Assume that ˉf0,ˉg0,ˉf(0,),ˉg=0, and Z1<Z2, then for each η(Z1,Z2) and ζ(Z3,), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (3) Assume that ˉf0,ˉg0,ˉg(0,),ˉf=0, and Z3<Z4, then for each η(Z1,) and ζ(Z3,Z4), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (4) Assume that ˉf0,ˉg0(0,),ˉf=ˉg=0, then for each η(Z1,) and ζ(Z3,), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (5) Assume that ˉf,ˉg(0,),ˉf0= or ˉf,ˉg(0,),ˉg0=, then for each η(0,Z2) and ζ(0,Z4), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (6) Assume that ˉf0=,ˉg=0,ˉf(0,) or ˉf(0,),ˉg=0,ˉg0=, then for each η(0,Z2) and ζ(0,), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (7) Assume that ˉf0=,ˉg(0,),ˉf=0 or ˉg(0,),ˉg0=,ˉf=0, then for each η(0,) and ζ(0,Z4), the system of fractional q-difference equations (1.1) has at least one positive solution.

    (8) Assume that ˉf0=,ˉf=ˉg=0 or ˉf=ˉg=0,ˉg0=, then for each η(0,) and ζ(0,), the system of fractional q-difference equations (1.1) has at least one positive solution.

    Proof. We will only prove the Cases (1) and (6). Since the other proofs are similar, so we omit.

    We firstly prove the Case (1). Since η(Z1,Z2) and ζ(Z3,Z4), there exists ε>0 such that

    0<ψ2(θ22α2M2)ˉf0εηψ1(12M1)ˉf+ε,0<ψ4(θ22α2M4)ˉg0εζψ3(12M2)ˉg+ε (3.24)

    From the definitions of ˉf0 and ˉg0, we easily know that there exists R1>0 such that

    f(t,x+φ,y+φ)(ˉf0ε)ϕ1(x+y),t[θ,1],x,y0,x+yR1, (3.25)
    g(t,x+φ,y+φ)(ˉg0ε)ϕg(x+y),t[θ,1],x,y0,x+yR1. (3.26)

    Let W1={(x,y)E:(x,y)E<R1}. By (3.24), (3.25) and Lemma 2.5, for any (x,y)PW1, we can get

    Aη(x,y)(θ)=10G(θ,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs1θG(θ,qs)ϕ11(ηΓq(γ)sθ(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs1θG(θ,qs)ϕ11(η(ˉf0ε)Γq(γ)sθ(sqτ)(γ1)ϕ1(x(τ)+y(τ))dqτ)dqsθα11θG(1,qs)ϕ11(η(ˉf0ε)Γq(γ)sθ(sqτ)(γ1)ϕ1(θα1(x,y)E)dqτ)dqsθα11θG(1,qs)ϕ11(η(ˉf0ε)Γq(γ)sθ(sqτ)(γ1)dqτϕ1(θα1(x,y)E)dqτ)dqsθα11θG(1,qs)ψ12(η(ˉf0ε)Γq(γ)sθ(sqτ)(γ1)dqτ)θα1(x,y)Edqs=θ2α2ψ12(η(ˉf0ε))1θG(1,qs)ψ12(1Γq(γ)sθ(sqτ)(γ1)dqτ)dqsR1=θ2α2ψ12(η(ˉf0ε))M3R1(x,y)E2. (3.27)

    By (3.24), (3.26) and Lemma 2.5, for any (x,y)PW1, we have

    Aζ(x,y)(θ)=10G(θ,qs)ϕ12(ζΓq(γ)s0(sqτ)(γ1)g(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs1θG(θ,qs)ϕ12(ζΓq(γ)sθ(sqτ)(γ1)g(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs1θG(θ,qs)ϕ12(ζΓq(γ)sθ(sqτ)(γ1)(ˉg0ε)ϕ2(x(τ)+y(τ))dqτ)dqsθα11θG(1,qs)ϕ12(ζΓq(γ)sθ(sqτ)(γ1)(ˉg0ε)ϕ2(θα1(x,y)E)dqτ)dqsθα11θG(1,qs)ψ14(ζ(ˉg0ε)Γq(γ)sθ(sqτ)(γ1)dqτ)θα1(x,y)Edqs=θ2α2ψ14(ζ(ˉg0ε))1θG(1,qs)ψ14(1Γq(γ)sθ(sqτ)(γ1)dqτ)dqs(x,y)E(x,y)E2. (3.28)

    By (3.27) and (3.28), we have

    A(x,y)E=Aη(x,y)+Aζ(x,y)(x,y)E,(x,y)PW1. (3.29)

    Let F(t,u)=max0x+yuf(t,x+φ,y+φ), G(t,u)=max0x+yug(t,x+φ,y+φ). Then we have

    f(t,x+φ,y+φ)F(t,u),t[0,1], x,y0, x+yu,
    g(t,x+φ,y+φ)G(t,u),t[0,1], x,y0, x+yu.

    Similar to the proof of [33], we know that

    limsupu+maxt[0,1]F(t,u)ϕ1(u)ˉf,  limsupu+maxt[0,1]G(t,u)ϕ2(u)ˉg.

    Clearly, we know that there exists ¯R2>0 such that

    F(t,u)ϕ1(u)limsupu+maxt[0,1]F(t,u)ϕ1(u)+εˉf+ε,u¯R2,t[0,1],
    G(t,u)ϕ2(u)limsupu+maxt[0,1]G(t,u)ϕ2(u)+εˉg+ε,u¯R2,t[0,1].

    Hence, we have

    F(t,u)(ˉf+ε)ϕ1(u),  G(t,u)(ˉg+ε)ϕ2(u)  t[0,1],u¯R2. (3.30)

    Let R2= max {2R1,¯R2}, and W2={(x,y)E:(x,y)E<R2}, for any (x,y)PW2, we get

    f(t,x+φ,y+φ)F(t,(x,y)E),t[0,1], (3.31)
    g(t,x+φ,y+φ)G(t,(x,y)E),t[0,1]. (3.32)

    By (3.30)–(3.32), for any (x,y)PW2, we have

    Aη(x,y)(t)=10G(t,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs10G(1,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)F(τ,(x,y)E)dqτ)dqs10G(1,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)(ˉf+ε)ϕ1((x,y)E)dqτ)dqs10G(1,qs)ψ11(η(ˉf+ε)Γq(γ)s0(sqτ)(γ1)dqτ)(x,y)Edqs=ψ11(η(ˉf+ε))10G(1,qs)ψ11(1Γq(γ)s0(sqτ)(γ1)dqτ)dqsR2=ψ11(η(ˉf+ε))10G(1,qs)ψ11(sγΓq(γ+1))dqsR2(x,y)E2, (3.33)

    and

    Aζ(x,y)(t)=10G(t,qs)ϕ12(ζΓq(γ)s0(sqτ)(γ1)g(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs10G(1,qs)ϕ12(ζΓq(γ)s0(sqτ)(γ1)G(τ,(x,y)E)dqτ)dqs10G(1,qs)ϕ12(ζΓq(γ)s0(sqτ)(γ1)(ˉg+ε)ϕ2((x,y)E)dqτ)dqs10G(1,qs)ψ13(ζ(ˉg+ε)Γq(γ)s0(sqτ)(γ1)dqτ)(x,y)Edqs=ψ13(ζ(ˉg+ε))10G(1,qs)ψ13(sγΓq(γ+1))dqsR2(x,y)E2. (3.34)

    By (3.33) and (3.34), we have

    A(x,y)E=Aη(x,y)+Aζ(x,y)(x,y)E,(x,y)PW2. (3.35)

    By (3.29), (3.35) and Lemma 2.1, we know that A has at least one fixed point (x,y)P(¯W2W1), so the system of fractional q-difference equations (1.1) has at least one positive solution. The proof of the case (1) is completed.

    In the following, we prove the Case (6). Since ˉf0=,ˉf(0,),ˉg=0, we can easily know that there exist ε>0 and R3>0 such that

    ψ2(θ22αM3)ε<η<ψ1(12M1)¯f+ε, (3.36)
    0<ζ<ψ3(12M2)1ε, (3.37)

    and

    f(t,x+φ,y+φ)1εϕ1(x+y),t[θ,1],x,y>0,0x+yR3. (3.38)

    Let W3={(x,y)E:(x,y)E<R3}. For t[θ,1], (x,y)PW3, we easily know that

    mint[θ,1](x(t)+y(t))θα1(x,y)E.

    By (3.36) and (3.38), we have

    Aη(x,y)(θ)=10G(θ,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs1θG(θ,qs)ϕ11(ηΓq(γ)sθ(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqsθα11θG(1,qs)ϕ11(ηΓq(γ)εsθ(sqτ)(γ1)ϕ1(x(τ)+y(τ))dqτ)dqsθα11θG(1,qs)ϕ11(ηΓq(γ)εsθ(sqτ)(γ1)ϕ1(θα1(x,y)E)dqτ)dqsθ2α2ψ12(ηε)1θG(1,qs)ψ12(1Γq(γ)sθ(sqτ)(γ1)dqτ)dqs=θ2α2ψ12(ηε)M3(x,y)E(x,y)E. (3.39)

    So by (3.39), we have

    A(x,y)EAη(x,y)(x,y)E,(x,y)PW3. (3.40)

    Similar to the proof of [33], we obtain

    limsupu+maxt[0,1]F(t,u)ϕ1(u)ˉf, limsupu+maxt[0,1]G(t,u)ϕ2(u)=0.

    So we know that for above ε>0 in (3.36) and (3.37), there exists ¯R4>0 such that

    F(t,u)ϕ1(u)limsupu+maxt[0,1]F(t,u)ϕ1(u)+εˉf+ε, t[0,1],u¯R4,
    G(t,u)ϕ2(u)limsupu+maxt[0,1]G(t,u)ϕ2(u)+εε, t[0,1],u¯R4,

    so we have

    F(t,u)¯(f+ε)ϕ1(u),t[0,1],u¯R4,
    G(t,u)εϕ2(u),t[0,1],u¯R4.

    Let R4= max {2R3,¯R4} and W4={(x,y)E:(x,y)E<R4}. We easily have

    f(t,x+φ,y+φ)F(t,(x,y)E), t[0,1],
    g(t,x+φ,y+φ)G(t,(x,y)E), t[0,1].

    Hence, for any t[0,1] and (x,y)PW4, we get

    Aη(x,y)(t)=10G(t,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)f(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs10G(1,qs)ϕ11(ηΓq(γ)s0(sqτ)(γ1)F(τ,(x,y)E)dqτ)dqs10G(1,qs)ϕ11(η(ˉf+ε)Γq(γ)s0(sqτ)(γ1)ϕ1((x,y)E)dqτ)dqs10G(1,qs)ψ11(η(ˉf+ε)Γq(γ)s0(sqτ)(γ1)dqτ)dqs(x,y)E=ψ11(η(ˉf+ε))10G(1,qs)ψ11(sγΓq(γ+1))dqs(x,y)E(x,y)E2, (3.41)

    and

    Aζ(x,y)(t)=10G(t,qs)ϕ12(ζΓq(γ)s0(sqτ)(γ1)g(τ,x(τ)+φ(τ),y(τ)+φ(τ))dqτ)dqs10G(1,qs)ϕ12(ζΓq(γ)s0(sqτ)(γ1)G(τ,(x,y)E)dqτ)dqs10G(1,qs)ϕ12(ζεΓq(γ)s0(sqτ)(γ1)ϕ2((x,y)E)dqτ)dqs10G(1,qs)ψ13(ζεΓq(γ)s0(sqτ)(γ1)dqτ)dqs(x,y)E=ψ13(ζε)10G(1,qs)ψ13(sγΓq(γ+1))dqs(x,y)E(x,y)E2. (3.42)

    So by (3.41) and (3.42), we have

    A(x,y)E=Aη(x,y)+Aζ(x,y)(x,y)E,(x,y)PW4. (3.43)

    By (3.40), (3.43) and Lemma 2.1, we know that A has at least one fixed point (x,y)P(¯W4W3). Hence the system of fractional q-difference equations (1.1) has at least one positive solution.

    Example 4.1. We consider the following system of fractional q-difference equations:

    {D12q(ϕ1(D52qx))(t)=ηf(t,x(t),y(t)),  0<t<1,D12q(ϕ2(D52qy))(t)=ζg(t,x(t),y(t)),  0<t<1,x(0)=Dqx(0)=0, Dqx(1)=1, D52qx(0)=0,y(0)=Dqy(0)=0, Dqy(1)=1, D52qy(0)=0, (4.1)

    where q=12, ϕ1(u)=u,ϕ2(u)=|u|1u. Take f(t,x,y)=t(x+y2φ(t))2, g(t,x,y)=t(x+y), where φ(t)=4+27t32. By a simple calculation we get

    f0=limsupx+y0+maxt[0,1]f(t,x+φ(t),y+φ(t))ϕ1(x+y)=limsupx+y0+maxt[0,1]t(x+y)=0,
    g0=limsupx+y0+maxt[0,1]g(t,x+φ(t),y+φ(t))ϕ2(x+y)=limsupx+y0+maxt[0,1]t(x+y+2φ(t))=27(4+2),
    g=liminfx+ymint[θ,1]g(t,x+φ(t),y+φ(t))ϕ2(x+y)=liminfx+ymint[θ,1]t(x+y+2φ(t))=;

    and ψ1(x)=ψ2(x)=x, ψ3(x)=ψ4(x)=1, D40.6465.

    Then, for each η(0,) and ζ(0,0.6465), by Theorem 3.1 Case (6) we obtain that the system (4.1) has at least one positive solution.

    Example 4.2. We consider the following system of fractional q-difference equations:

    {D12q(ϕ1(D52qx))(t)=ηf(t,x(t),y(t)),  0<t<1,D12q(ϕ2(D52qy))(t)=ζg(t,x(t),y(t)),  0<t<1,x(0)=Dqx(0)=0, Dqx(1)=1, D52qx(0)=0,y(0)=Dqy(0)=0, Dqy(1)=1, D52qy(0)=0, (4.2)

    where q=12, ϕ1(u)=u,ϕ2(u)=|u|1u. Take f(t,x,y)=t(x+y2φ(t))arctan(x+y2φ(t)), g(t,x,y)=tx+y, where φ(t)=4+27t32. By a simple calculation we get

    ˉf0=liminfx+y0+mint[θ,1]f(t,x+φ(t),y+φ(t))ϕ1(x+y)=liminfx+y0+mint[θ,1]tarctan(x+y)=,
    ˉg=limsupx+ymaxt[0,1]g(t,x+φ(t),y+φ(t))ϕ2(x+y)=limsupx+ymaxt[0,1]tx+y+2φ(t)=0,
    ˉf=limsupx+ymaxt[0,1]f(t,x+φ(t),y+φ(t))ϕ1(x+y)=limsupx+ymaxt[0,1]tarctan(x+y)=2π;

    and ψ1(x)=ψ2(x)=x, ψ3(x)=ψ4(x)=1, Γ12(52)1.1906, Γ12(32)0.9209, M10.2991, Z22.6259.

    Then, for each η(0,2.6259) and ζ(0,), by Theorem 3.2(6) we obtain that the system (4.2) has at least one positive solution.

    The system of fractional q-difference equations plays an important role in the study of many fields, such as quantum mechanics, mathematical physics equations and so on, for example, see [16,17,24,35] and the references therein. In [35], by using some classical fixed point theorems, the authors studied the existence of nontrivial solutions of a system of fractional q-difference equations with Riemann-Stieltjes integrals conditions. In this paper, we investigate the existence of positive solutions for a system of fractional q-difference equations with generalized p-Laplacian operators and two parameters. The system in this paper is different from that of [35]. We give some assumptions which are combinations of superlinearity and sublinearity of the nonlinear terms f and g. Under those assumptions, by using Guo-Krasnosel'skii fixed point theorem, we obtain some existence results of positive solutions in terms of different values of the parameters η and ζ. In fact, since the system studied in this paper contains generalized p-Laplacian operators, the obtained results in this paper can enrich the relevant knowledge of theories for the system of fractional q-difference equations and expand the range of the possible applications. However, this study still has certain limitations, as we only investigated the existence of positive solutions. In the future, some further work can continue to be considered such as the uniqueness and multiplicity of positive solutions and iterative sequences of positive solutions, the case where the nonlinear terms may be changing sign or the generalized p-Laplacian operator becomes a p(t) -Laplacian operator, etc.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The project is supported by the National Natural Science Foundation of China (11801322; 12371173) and Shandong Natural Science Foundation(ZR2021MA064).

    The authors would like to thank reviewers for their valuable comments, which help to enrich the content of this paper.

    The authors declare there is no conflicts of interest.



    [1] Y. Zhao, S. Sun, Z. Han, Q. Li, The existence of multiple positive solutions for boundary value problems of nonlinear fractional differential equations, Commun. Nonlinear Sci. Numer. Simul., 16 (2011), 2086–2097. https://doi.org/10.1016/j.cnsns.2010.08.017 doi: 10.1016/j.cnsns.2010.08.017
    [2] X. Zhang, J. Jiang, L. Liu, Y. Wu, Extremal solutions for a class of tempered fractional turbulent flow equations in a porous medium, Math. Probl. Eng., 2020 (2020), 2492193. https://doi.org/10.1155/2020/2492193 doi: 10.1155/2020/2492193
    [3] Z. Bai, H. L¨u, Positive solutions for boundary value problem of nonlinear fractional differential equation, J. Math. Anal. Appl., 311 (2005), 495–505. https://doi.org/10.1016/j.jmaa.2005.02.052 doi: 10.1016/j.jmaa.2005.02.052
    [4] X. Zhang, L. Yu, J. Jiang, Y. Wu, Y. Cui, Positive solutions for a weakly singular Hadamard-type fractional differential equation with changing-sign nonlinearity, J. Funct. Space, 2020 (2020), 5623589. https://doi.org/10.1155/2020/5623589 doi: 10.1155/2020/5623589
    [5] X. Zhang, P. Chen, H. Tian, Y. Wu, The iterative properties for positive solutions of a tempered fractional equation, Fractal Fract., 7 (2023), 761. https://doi.org/10.3390/fractalfract7100761 doi: 10.3390/fractalfract7100761
    [6] X. Zhang, L. Yu, J. Jiang, Y. Wu, Y. Cui, Solutions for a singular Hadamard-type fractional differential equation by the spectral construct analysis, J. Funct. Space, 2020 (2020), 8392397. https://doi.org/10.1155/2020/8392397 doi: 10.1155/2020/8392397
    [7] X. Zhang, P. Xu, Y. Wu, B. Wiwatanapataphee, The uniqueness and iterative properties of solutions for a general Hadamard-type singular fractional turbulent flow model, Nonlinear Anal.-Model. Control, 27 (2022), 428–444. https://doi.org/10.15388/namc.2022.27.25473 doi: 10.15388/namc.2022.27.25473
    [8] X. Zhang, P. Chen, H. Tian, Y. Wu, Upper and lower solution method for a singular tempered fractional equation with a p-Laplacian operator, Fractal Fract., 7 (2023), 522. https://doi.org/10.3390/fractalfract7070522 doi: 10.3390/fractalfract7070522
    [9] Y. Li, G. Li, Positive solutions of p-Laplacian fractional differential equations with integral boundary value conditions, J. Nonlinear Sci. Appl., 9 (2016), 717–726.
    [10] T. Chen, W. Liu, Z. Hu, A boundary value problem for fractional differential equation with p-Laplacian operator at resonance, Nonlinear Anal. Theory Methods Appl., 75 (2012), 3210–3217. https://doi.org/10.1016/j.na.2011.12.020 doi: 10.1016/j.na.2011.12.020
    [11] L. Zhang, W. Zhang, X. Liu, M. Jia, Positive solutions of fractional p-laplacian equations with integral boundary value and two parameters, J. Inequal. Appl., 2020 (2020), 2. https://doi.org/10.1186/s13660-019-2273-6 doi: 10.1186/s13660-019-2273-6
    [12] A. Ahmadkhanlu, On the existence and multiplicity of positive solutions for a p-Laplacian fractional boundary value problem with an integral boundary condition, Filomat, 37 (2023), 235–250. https://doi.org/10.2298/FIL2301235A doi: 10.2298/FIL2301235A
    [13] X. Zhang, D. Kong, H. Tian, Y. Wu, B. Wiwatanapataphee, An upper-lower solution method for the eigenvalue problem of Hadamard-type singular fractional differential equation, Nonlinear Anal.-Model. Control, 27 (2022), 789–802. https://doi.org/10.15388/namc.2022.27.27491 doi: 10.15388/namc.2022.27.27491
    [14] Z. Han, H. Lu, C. Zhang, Positive solutions for eigenvalue problems of fractional differential equation with generalized p-Laplacian, Appl. Math. Comput., 257 (2015), 526–536. https://doi.org/10.1016/j.amc.2015.01.013 doi: 10.1016/j.amc.2015.01.013
    [15] F. H. Jackson, On q-definite integrals, Quart. J. Pure Appl. Math., 41 (1910), 193–203.
    [16] W. A. Al-Salam, Some fractional q-integrals and q-derivatives, Proc. Edinb. Math. Soc., 15 (1966), 135–140. https://doi.org/10.1017/S0013091500011469 doi: 10.1017/S0013091500011469
    [17] R. P. Agarwal, Certain fractional q-integrals and q-derivatives, Proc. Camb. Philos. Soc., 66 (1969), 365–370. https://doi.org/10.1017/S0305004100045060 doi: 10.1017/S0305004100045060
    [18] R. A. C. Ferreira, Positive solutions for a class of boundary value problems with fractional q-differences, Comput. Math. Appl., 61 (2011), 367–373. https://doi.org/10.1016/j.camwa.2010.11.012 doi: 10.1016/j.camwa.2010.11.012
    [19] C. Zhai, J. Ren, Positive and negative solutions of a boundary value problem for a fractional q-difference equation, Adv. Differ. Equations, 2017 (2017). https://doi.org/10.1186/s13662-017-1138-x doi: 10.1186/s13662-017-1138-x
    [20] J. Mao, Z. Zhao, C. Wang, The unique iterative positive solution of fractional boundary value problem with q-difference, Appl. Math. Lett., 100 (2020), 106002. https://doi.org/10.1016/j.aml.2019.106002 doi: 10.1016/j.aml.2019.106002
    [21] M. Jiang, S. Zhong, Existence of solutions for nonlinear fractional q-difference equations with Riemann-Liouville type q-derivatives, J. Appl. Math. Comput., 47 (2015), 429–459. https://doi.org/10.1007/s12190-014-0784-3 doi: 10.1007/s12190-014-0784-3
    [22] X. Li, Z. Han, S. Sun, L. Sun, Eigenvalue problems of fractional q-difference equations with generalized p-Laplacian, Appl. Math. Lett., 57 (2016), 46–53. https://doi.org/10.1016/j.aml.2016.01.003 doi: 10.1016/j.aml.2016.01.003
    [23] J. Wang, C. Yu, B. Zhang, S. Wang, Positive solutions for eigenvalue problems of fractional q-difference equation with ϕ-Laplacian, Adv. Differ. Equations, 2021 (2021), 499. https://doi.org/10.1186/s13662-021-03652-x doi: 10.1186/s13662-021-03652-x
    [24] C. Yu, S. Li, J. Li, J. Wang, Triple-positive solutions for a nonlinear singular fractional q-difference equation at resonance, Fractal Fract., 6 (2022), 689. https://doi.org/10.3390/fractalfract6110689 doi: 10.3390/fractalfract6110689
    [25] G. Wang, Twin iterative positive solutions of fractional q-difference Schr¨odinger equations, Appl. Math. Lett., 76 (2018), 103–109. https://doi.org/10.1016/j.aml.2017.08.008 doi: 10.1016/j.aml.2017.08.008
    [26] X. Li, Z. Han, S. Sun, P. Zhao, Existence of solutions for fractional q-difference equation with mixed nonlinear boundary conditions, Adv. Differ. Equations, 2014 (2014), 326. https://doi.org/10.1186/1687-1847-2014-326 doi: 10.1186/1687-1847-2014-326
    [27] S. Liang, J. Zhang, Existence and uniqueness of positive solutions for three-point boundary value problem with fractional q-differences, J. Appl. Math. Comput., 40 (2012), 277–288. https://doi.org/10.1007/s12190-012-0551-2 doi: 10.1007/s12190-012-0551-2
    [28] X. Li, Z. Han, S. Sun, Existence of positive solutions of nonlinear fractional q-difference equation with parameter, Adv. Differ. Equations, 2013 (2013), 260. https://doi.org/10.1186/1687-1847-2013-260 doi: 10.1186/1687-1847-2013-260
    [29] J. Ma, J. Yang, Existence of solutions for multi-point boundary value problem of fractional q-difference equation, Electron. J. Qual. Theory Differ. Equations, 2011 (2011), 1–10.
    [30] B. Ahmad, S. Etemad, M. Ettefagh, S. Rezapour, On the existence of solutions for fractional q-difference inclusions with q-antiperiodic boundary conditions, Math. Bull. Soc. Math. Sci. Rom., 59 (2016), 119–134. https://www.jstor.org/stable/26407454
    [31] S. Li, X. Zhang, Y. Wu, L. Caccetta, Extremal solutions for p-Laplacian differential systems via iterative computation, Appl. Math. Lett., 26 (2013), 1151–1158. https://doi.org/10.1016/j.aml.2013.06.014 doi: 10.1016/j.aml.2013.06.014
    [32] J. He, X. Song, The uniqueness of solution for a class of fractional order nonlinear systems with p-Laplacian operator, Abstr. Appl. Anal., 2014 (2014), 921209. https://doi.org/10.1155/2014/921209 doi: 10.1155/2014/921209
    [33] X. Hao, H. Wang, L. Liu, Y. Cui, Positive solutions for a system of nonlinear fractional nonlocal boundary value problems with parameters and p-Laplacian operator, Bound. Value Probl., 2017 (2017), 182. https://doi.org/10.1186/s13661-017-0915-5 doi: 10.1186/s13661-017-0915-5
    [34] X. Zhang, L. Liu, Y. Wu, The uniqueness of positive solution for a singular fractional differential system involving derivatives, Commun. Nonlinear Sci., 18 (2013), 1400–1409. https://doi.org/10.1016/j.cnsns.2012.08.033 doi: 10.1016/j.cnsns.2012.08.033
    [35] C. Yu, S. Wang, J. Wang, J. Li, Solvability criterion for fractional q-integro-difference system with Riemann-Stieltjes integrals conditions, Fractal Fract., 6 (2022), 554. https://doi.org/10.3390/fractalfract6100554 doi: 10.3390/fractalfract6100554
    [36] D. Guo, V. Lakshmikantham, Nonlinear Problems in Abstract Cones, Academic Press, San Diego, 5 (2014).
    [37] H. Wang, On the number of positive solutions of nonlinear systems, J. Math. Anal. Appl., 281 (2003), 287–306. https://doi.org/10.1016/S0022-247X(03)00100-8 doi: 10.1016/S0022-247X(03)00100-8
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1167) PDF downloads(63) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog