Next Article in Journal
Fractional-Order Calculus-Based Data Augmentation Methods for Environmental Sound Classification with Deep Learning
Next Article in Special Issue
Systems of Riemann–Liouville Fractional Differential Equations with ρ-Laplacian Operators and Nonlocal Coupled Boundary Conditions
Previous Article in Journal
A Stochastic Bayesian Regularization Approach for the Fractional Food Chain Supply System with Allee Effects
Previous Article in Special Issue
Analytic Resolving Families for Equations with the Dzhrbashyan–Nersesyan Fractional Derivative
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solvability Criterion for Fractional q-Integro-Difference System with Riemann-Stieltjes Integrals Conditions

1
Interdisciplinary Research Institute, Faculty of Science, Beijing University of Technology, Beijing 100124, China
2
College of Sciences, Hebei University of Science and Technology, Shijiazhuang 050018, China
*
Authors to whom correspondence should be addressed.
Fractal Fract. 2022, 6(10), 554; https://doi.org/10.3390/fractalfract6100554
Submission received: 23 August 2022 / Revised: 18 September 2022 / Accepted: 27 September 2022 / Published: 29 September 2022

Abstract

:
Due to the great application potential of fractional q-difference system in physics, mechanics and aerodynamics, it is very necessary to study fractional q-difference system. The main purpose of this paper is to investigate the solvability of nonlinear fractional q-integro-difference system with the nonlocal boundary conditions involving diverse fractional q-derivatives and Riemann-Stieltjes q-integrals. We acquire the existence results of solutions for the systems by applying Schauder fixed point theorem, Krasnoselskii’s fixed point theorem, Schaefer’s fixed point theorem and nonlinear alternative for single-valued maps, and a uniqueness result is obtained through the Banach contraction mapping principle. Finally, we give some examples to illustrate the main results.

1. Introduction

In the early twentieth century, Jackson [1] proposed a new mathematical direction of q-calculus, and it plays an indispensable role in the fields of nuclear, conformal quantum mechanics and dynamics. In the 1960s, Agarwal [2] and Al-Salam [3] put forward a novel concept of fractional q-calculus, its relevant application and development can be seen in the literature [4,5,6]. Compared with classical q-calculus, fractional q-calculus can more accurately describe some phenomena in nature, and many practical problems can be abstracted into fractional q-difference equations or a system of fractional q-difference equations by mathematical modeling. In recent years, abundant theoretical achievements have been made in the research of boundary value problems (BVPs) for fractional q-difference equations, according to the literature [7,8,9,10,11,12,13,14,15,16] and the references therein.
Riemann-Stieltjes integral is a generalization of Riemann integral. As well as we known, the classical Riemann-Stieltjes integral can be widely applied in several areas of analysis, such as probability theory, stochastic processes, physics, econometrics, biometrics and informetrics and so on. BVPs with Riemann-Stieltjes integral boundary condition (BC) have been considered as both multi-point and integral type BCs are treated in a single framework. In recent years, some interesting results about the existence of solutions for nonlinear fractional differential equations with the Riemann-Stieltjes integral BC have been researched, see [17,18] and the references therein.
Nowadays, the system of nonlinear fractional differential equations has important applications in engineering, economy and other fields. This is mainly because the effect of using fractional calculus to solve problems is more practical and efficient than that of classical calculus. Over the years, the BVPs for a system of fractional differential equations have developed rapidly, and numerous mature conclusions have been obtained, which can be referred to the literature [19,20,21,22,23,24,25].
In [24], Tudorache, A. and Luca, R. applied the Guo-Krasnoselskii fixed point theorem to study the existence of solutions for a system of fractional differential equations with p-Laplacian operators
D 0 + α 1 ( φ ϱ 1 ( D 0 + β 1 x ( t ) ) ) + λ f ( t , x ( t ) , y ( t ) ) = 0 , t ( 0 , 1 ) , D 0 + α 2 ( φ ϱ 2 ( D 0 + β 2 y ( t ) ) ) + μ g ( t , x ( t ) , y ( t ) ) = 0 , t ( 0 , 1 ) ,
with the nonlocal BCs
x ( j ) ( 0 ) = 0 , j = 0 , , n 2 ; D 0 + β 1 x ( 0 ) = 0 , D 0 + γ 0 x ( 1 ) = i = 1 p 0 1 D 0 + γ i y ( t ) d H i t , y ( j ) ( 0 ) = 0 , j = 0 , , m 2 ; D 0 + β 2 y ( 0 ) = 0 , D 0 + δ 0 y ( 1 ) = i = 1 q 0 1 D 0 + δ i x ( t ) d K i t .
In [25], Luca, R. considered the existence of solutions of the nonlinear system of fractional differential equations by using a variety of fixed point theorems
D 0 + α x ( t ) + f ( t , x ( t ) , y ( t ) , I 0 + θ 1 x ( t ) , I 0 + σ 1 y ( t ) ) = 0 , t ( 0 , 1 ) , D 0 + β y ( t ) + g ( t , x ( t ) , y ( t ) , I 0 + θ 2 x ( t ) , I 0 + σ 2 y ( t ) ) = 0 , t ( 0 , 1 ) ,
with the nonlocal BCs
x ( 0 ) = x ( 0 ) = = x ( n 2 ) ( 0 ) = 0 , D 0 + γ 0 x ( 1 ) = i = 1 p 0 1 D 0 + γ i y ( t ) d H i ( t ) , y ( 0 ) = y ( 0 ) = = y ( m 2 ) ( 0 ) = 0 , D 0 + δ 0 y ( 1 ) = i = 1 q 0 1 D 0 + δ i x ( t ) d K i ( t ) .
Despite quite a number of contributions dealing with the solvability for the system of classical fractional difference equations. However, as the generalization of the above system, limited work has been done in the nonlinear system of fractional q-difference equations. In particular, there is little research on the existence and uniqueness of solutions for the system of fractional q-difference equations with Riemann-Stieltjes integral BC. To fill this gap, we investigate the system of nonlinear fractional q-difference equations
( D q α u ) ( t ) + P ( t , u ( t ) , v ( t ) , I q ω 1 u ( t ) , I q δ 1 v ( t ) ) = 0 , ( D q β v ) ( t ) + Q ( t , u ( t ) , v ( t ) , I q ω 2 u ( t ) , I q δ 2 v ( t ) ) = 0 ,
with the nonlocal BCs
u ( 0 ) = D q u ( 0 ) = = D q n 2 u ( 0 ) = 0 , D q ζ 0 u ( 1 ) = 0 1 D q ζ v ( t ) d q H ( t ) , v ( 0 ) = D q v ( 0 ) = = D q m 2 v ( 0 ) = 0 , D q ξ 0 v ( 1 ) = 0 1 D q ξ u ( t ) d q K ( t ) ,
where t ( 0 , 1 ) , 0 < q < 1 , α ( n 1 , n ] , β ( z 1 , z ] , n , z N , n 2 and z 2 , ω 1 , ω 2 , δ 1 , δ 2 > 0 , 0 ζ < β 1 , 0 ξ < α 1 , ζ 0 [ 0 , α 1 ) ,   ξ 0 [ 0 , β 1 ) ,   D q i denotes the Riemann-Liouville q-derivative of order i ( i = α , β , ζ 0 , ζ , ξ 0 , ξ ), I q ϖ is the Riemann-Liouville q-integral of order ϖ ( ϖ = ω 1 , ω 2 , δ 1 , δ 2 ), P and Q are nonlinear functions. The BCs include Riemann-Stieltjes integrals, where H ( t ) , K ( t ) are the bounded variation functions. In the case where H ( t ) = K ( t ) = t , the Riemann–Stieltjes integrals in (2) reduce to the classical q-integral.
The present paper is bulit up as follows. The second part offers the necessary definitions, lemmas and theorems needed in the following. The third part obtains the important conclusions by applying various fixed point theorems, including nine theorems or corollaries. In the final part, four examples are provided to verify our main results.

2. Preliminaries

In this section, we present some definitions, lemmas and theorems.
Definition 1
([11]). Let β 0 and f be a function defined on [ 0 , 1 ] . The fractional q-integral of the Riemann-Liouville type is
( I q β f ) ( s ) = 1 Γ q ( β ) 0 s ( s q t ) ( β 1 ) f ( t ) d q t , β > 0 , s [ 0 , 1 ] .
Obviously, ( I q β f ) ( s ) = ( I q f ) ( s ) , when β = 1 .
Definition 2
([11]). The fractional q-derivative of the Riemann-Liouville type of order β 0 is defined by ( D q 0 f ) ( s ) = f ( s ) and
( D q β f ) ( s ) = ( D q l I q l β f ) ( s ) , β > 0 , s [ 0 , 1 ] ,
where l is the smallest integer greater than or equal to β.
Lemma 1
([11]). Let α , β 0 and f be a function defined on [ 0 , 1 ] . Then, the following formulas hold:
1.
( I q β I q α f ) ( x ) = ( I q α + β f ) ( x ) ,
2.
( D q α I q α f ) ( x ) = f ( x ) .
Lemma 2
([11]). Let α > 0 and p be a positive integer. Then, the following equality holds:
( I q α D q p f ) ( x ) = ( D q p I q α f ) ( x ) k = 0 p 1 x α p + k Γ q ( α + k p + 1 ) ( D q k f ) ( 0 ) .
Lemma 3.
If x C [ 0 , 1 ] , then for κ > 0 , we get
| I q κ x ( t ) | x Γ q ( κ ) ,
where x = sup t [ 0 , 1 ] | x ( t ) | .
Proof. 
According to Definition 1, this lemma clearly holds. □
Definition 3
([15]). The function f : I × R 4 R is called an S-Carathéodory function if and only if
 (i)
for each ( u , v , x , y ) R 4 , t f ( t , u , v , x , y ) is measurable on I;
 (ii)
for a.e. t I , ( u , v , x , y ) f ( t , u , v , x , y ) is continuous on R 4 ;
 (iii)
for each r > 0 , there exists ψ r ( t ) L 1 ( I , R + ) with t ψ r ( t ) L 1 ( I , R + ) on I such that max { | u | , | v | , | x | , | y | } r implies | f ( t , u , v , x , y ) | ψ r ( t ) , for a.e.I, where L 1 ( I , R + ) = { u X : 0 1 u ( t ) d q t exists } , and normed u L 1 = 0 1 | u ( t ) | d q t for all u L 1 ( I , R + ) .
Theorem 1
([26]). (Schauder fixed point theorem) Let D be a bounded closed convex set in E (D does not necessarily have an interior point), and A : D D is completely continuous, then A must have a fixed point in D.
Theorem 2
([12]). (Krasnoselskii’s fixed point theorem) Let K be a closed convex and nonempty subset of a Banach space X. Let T , S be the operators such that
 (i)
T u + S v K whenever u , v K ;
 (ii)
T is compact and continuous;
 (iii)
S is a contraction mapping.
Then, there exists z K such that z = T z + S z .
Theorem 3
([16]). (Schaefer’s fixed point theorem) Let T be a continuous and compact mapping of a Banach space X into itself, such that the set E = { x | x X : x = λ T x , 0 λ 1 } is bounded. Then T has a fixed point.
Theorem 4
([15]). (Nonlinear alternative for single-valued maps) Let E be a Banach space, let C be a closed and convex subset of E, and let U be an open subset of C and 0 U . Suppose that F : U ¯ C is a continuous, compact (that is, F ( U ¯ ) is a relatively compact subset of C) map. Then either
 (i)
F has a fixed point in U ¯ , or
 (ii)
there is a u U (the boundary of U in C) and λ ( 0 , 1 ) with u = λ F u .
Throughout this paper, we adopt the following assumptions:
( H 1 ) The functions P , Q C ( [ 0 , 1 ] × R 4 , R ) and for x i , y i R , there exist L i ( t ) , l i ( t ) C ( [ 0 , 1 ] , [ 0 , + ) ) , i = 1 , 2 , 3 , 4 , such that
| P ( t , x 1 , x 2 , x 3 , x 4 ) P ( t , y 1 , y 2 , y 3 , y 4 ) | i = 1 4 L i ( t ) | x i y i | , | Q ( t , x 1 , x 2 , x 3 , x 4 ) Q ( t , y 1 , y 2 , y 3 , y 4 ) | i = 1 4 l i ( t ) | x i y i | .
( H 1 ) The functions P , Q C ( [ 0 , 1 ] × R 4 , R ) and for x i , y i R , there exist real constants L i , l i > 0 , i = 1 , 2 , 3 , 4 , such that
| P ( t , x 1 , x 2 , x 3 , x 4 ) P ( t , y 1 , y 2 , y 3 , y 4 ) | i = 1 4 L i | x i y i | , | Q ( t , x 1 , x 2 , x 3 , x 4 ) Q ( t , y 1 , y 2 , y 3 , y 4 ) | i = 1 4 l i | x i y i | .
( H 1 ) The functions P , Q C ( [ 0 , 1 ] × R 4 , R ) and for x i , y i R , there exist real functions ρ i ( t ) , ϱ i ( t ) C ( [ 0 , 1 ] , R + ) , i = 1 , 2 , 3 , 4 , such that
| P ( t , x 1 , x 2 , x 3 , x 4 ) P ( t , y 1 , y 2 , y 3 , y 4 ) | i = 1 4 ρ i ( t ) | x i y i | , | Q ( t , x 1 , x 2 , x 3 , x 4 ) Q ( t , y 1 , y 2 , y 3 , y 4 ) | i = 1 4 ϱ i ( t ) | x i y i | .
( H 2 ) The functions P , Q C ( [ 0 , 1 ] × R 4 , R ) , and for x i R , there exist functions c i ( t ) , d i ( t ) C ( [ 0 , 1 ] , R + ) , and h i , m i ( 0 , 1 ) , i = 1 , 2 , 3 , 4 , such that
| P ( t , x 1 , x 2 , x 3 , x 4 ) | c 0 ( t ) + i = 1 4 c i ( t ) | x i | h i , | Q ( t , x 1 , x 2 , x 3 , x 4 ) | d 0 ( t ) + i = 1 4 d i ( t ) | x i | m i .
( H 3 ) The functions P , Q C ( [ 0 , 1 ] × R 4 , R ) , and for x i R , i = 1 , 2 , 3 , 4 , there exist functions σ 1 ( t ) , σ 2 ( t ) C ( [ 0 , 1 ] , R + ) such that
| P ( t , x 1 , x 2 , x 3 , x 4 ) | σ 1 ( t ) , | Q ( t , x 1 , x 2 , x 3 , x 4 ) | σ 2 ( t ) .
( H 4 ) The functions P , Q : [ 0 , 1 ] × R 4 R and for a.e. t [ 0 , 1 ] , x i R , there exist r 1 ( t ) , r 2 ( t ) , L i ( t ) , l i ( t ) C ( [ 0 , 1 ] , R + ) , i = 1 , 2 , 3 , 4 , such that
| P ( t , x 1 , x 2 , x 3 , x 4 ) | i = 1 4 L i ( t ) | x i | + r 1 ( t ) , | Q ( t , x 1 , x 2 , x 3 , x 4 ) | i = 1 4 l i ( t ) | x i | + r 2 ( t ) .
( H 4 ) The functions P , Q : [ 0 , 1 ] × R 4 R and for a.e. t [ 0 , 1 ] , x i R , there exist non-negative real numbers L i , l i ( i = 1 , 2 , 3 , 4 ) , and r 1 , r 2 , where at least one of r 1 and r 2 is positive, such that
| P ( t , x 1 , x 2 , x 3 , x 4 ) | i = 1 4 L i | x i | + r 1 , | Q ( t , x 1 , x 2 , x 3 , x 4 ) | i = 1 4 l i | x i | + r 2 .
( H 5 ) The functions P , Q : [ 0 , 1 ] × R 4 R and for a.e. t [ 0 , 1 ] , x i R , there exist functions p i ( t ) , q i ( t ) C ( [ 0 , 1 ] , R + ) , where p i ( t ) , q i ( t ) have at least one non-zero function, and there exist nondecreasing functions φ i , η i C ( [ 0 , ) , R + ) , i = 1 , 2 , 3 , 4 , such that
| P ( t , x 1 , x 2 , x 3 , x 4 ) | i = 1 4 p i ( t ) φ i ( | x i | ) + p 0 ( t ) , | Q ( t , x 1 , x 2 , x 3 , x 4 ) | i = 1 4 q i ( t ) η i ( | x i | ) + q 0 ( t ) .
For convenience, we denote
C 1 = 1 + 1 Γ q ( ω 1 ) , C 2 = 1 + 1 Γ q ( δ 1 ) , C 3 = max { C 1 , C 2 } , C 4 = 1 + 1 Γ q ( ω 2 ) , C 5 = 1 + 1 Γ q ( δ 2 ) , C 6 = max { C 4 , C 5 } , C 7 = 1 Γ q ( α ) + Γ q ( β ) | Ω | Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) + Γ q ( β ) | Ω | Γ q ( α ξ ) Γ q ( β ζ ) | 0 1 s β ζ 1 d q H ( s ) | · | 0 1 d q K ( s ) | , C 8 = 1 Γ q ( β ) + Γ q ( α ) | Ω | Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) + Γ q ( α ) | Ω | Γ q ( α ξ ) Γ q ( β ζ ) | 0 1 s α ξ 1 d q K ( s ) | · | 0 1 d q H ( s ) | , C 9 = Γ q ( α ) | Ω | Γ q ( α ζ 0 ) Γ q ( α ξ ) | 0 1 s α ξ 1 d q K ( s ) | + Γ q ( α ) | Ω | Γ q ( α ζ 0 ) Γ q ( α ξ ) | 0 1 d q K ( s ) | , C 10 = Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) | 0 1 d q H ( s ) | , C 11 = C 7 1 Γ q ( α ) , C 12 = C 8 1 Γ q ( β ) . Ω 1 = Γ q ( β ) Γ q ( β ζ ) 0 1 s β ζ 1 d q H ( s ) , Ω 2 = Γ q ( α ) Γ q ( α ξ ) 0 1 s α ξ 1 d q K ( s ) , Ω = Γ q ( α ) Γ q ( β ) Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) Ω 1 Ω 2 .

3. Criterion of Uniqueness and Existence

In this section, we show some existence and uniqueness results for the Systems (1)–(2).
Lemma 4.
Let h , k C ( 0 , 1 ) L 1 ( 0 , 1 ) and Ω 0 , then the system of fractional q-difference equations
D q α u ( t ) + h ( t ) = 0 , t ( 0 , 1 ) , D q β v ( t ) + k ( t ) = 0 , t ( 0 , 1 ) ,
with the coupled BCs (2) has a unique solution ( u ( t ) , v ( t ) ) , namely
u ( t ) = 1 Γ q ( α ) 0 t ( t q s ) ( α 1 ) h ( s ) d q s + t α 1 Ω [ Ω 1 Γ q ( β ξ 0 ) 0 1 ( 1 q s ) ( β ξ 0 1 ) k ( s ) d q s Ω 1 Γ q ( α ξ ) 0 1 0 s ( s q τ ) ( α ξ 1 ) h ( τ ) d q τ d q K ( s ) + Γ q ( β ) Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) 0 1 ( 1 q s ) ( α ζ 0 1 ) h ( s ) d q s Γ q ( β ) Γ q ( β ζ ) Γ q ( β ξ 0 ) 0 1 0 s ( s q τ ) ( β ζ 1 ) k ( τ ) d q τ d q H ( s ) ] , v ( t ) = 1 Γ q ( β ) 0 t ( t q s ) ( β 1 ) k ( s ) d q s + t β 1 Ω [ Ω 2 Γ q ( α ζ 0 ) 0 1 ( 1 q s ) ( α ζ 0 1 ) h ( s ) d q s Ω 2 Γ q ( β ζ ) 0 1 0 s ( s q τ ) ( β ζ 1 ) k ( τ ) d q τ d q H ( s ) + Γ q ( α ) Γ q ( β ξ 0 ) Γ q ( α ζ 0 ) 0 1 ( 1 q s ) ( β ξ 0 1 ) k ( s ) d q s Γ q ( α ) Γ q ( α ξ ) Γ q ( α ζ 0 ) 0 1 0 s ( s q τ ) ( α ξ 1 ) h ( τ ) d q τ d q K ( s ) ] , t [ 0 , 1 ] .
Proof. 
The proof is similar to the Lemma 2.1 in [24]. □
Let U = C [ 0 , 1 ] and V = U × U be the Banach spaces with the norms u = sup t [ 0 , 1 ] | u ( t ) | and ( u , v ) V = u + v , respectively. Nowdays, we introduce the operator T : V V , where T ( x , y ) = ( T 1 ( x , y ) , T 2 ( x , y ) ) for ( x , y ) V , and T 1 , T 2 : V U are defined by
T 1 ( u , v ) ( t ) = 1 Γ q ( α ) 0 t ( t q s ) ( α 1 ) F u v ( s ) d q s + Ω 1 t α 1 Ω Γ q ( β ξ 0 ) 0 1 ( 1 q s ) ( β ξ 0 1 ) · G u v ( s ) d q s Ω 1 t α 1 Ω Γ q ( α ξ ) 0 1 0 s ( s q τ ) ( α ξ 1 ) F u v ( τ ) d q τ d q K ( s ) + Γ q ( β ) t α 1 Ω Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) 0 1 ( 1 q s ) ( α ζ 0 1 ) F u v ( s ) d q s Γ q ( β ) t α 1 Ω Γ q ( β ζ ) Γ q ( β ξ 0 ) 0 1 0 s ( s q τ ) ( β ζ 1 ) G u v ( τ ) d q τ d q H ( s ) ,
and
T 2 ( u , v ) ( t ) = 1 Γ q ( β ) 0 t ( t q s ) ( β 1 ) G u v ( s ) d q s + Ω 2 t β 1 Ω Γ q ( α ζ 0 ) 0 1 ( 1 q s ) ( α ζ 0 1 ) · F u v ( s ) d q s Ω 2 t β 1 Ω Γ q ( β ζ ) 0 1 0 s ( s q τ ) ( β ζ 1 ) G u v ( τ ) d q τ d q H ( s ) + Γ q ( α ) t β 1 Ω Γ q ( β ξ 0 ) Γ q ( α ζ 0 ) 0 1 ( 1 q s ) ( β ξ 0 1 ) G u v ( s ) d q s Γ q ( α ) t β 1 Ω Γ q ( α ξ ) Γ q ( α ζ 0 ) 0 1 0 s ( s q τ ) ( α ξ 1 ) F u v ( τ ) d q τ d q K ( s ) ,
for t [ 0 , 1 ] and ( u , v ) V , where
F u v ( s ) = P ( s , u ( s ) , v ( s ) , I q ω 1 u ( s ) , I q δ 1 v ( s ) ) , G u v ( s ) = Q ( s , u ( s ) , v ( s ) , I q ω 2 u ( s ) , I q δ 2 v ( s ) ) .
According to Lemma 4, it is easy to see that ( u ( t ) , v ( t ) ) is a solution of the Systems (1)–(2) if and only if ( u ( t ) , v ( t ) ) is a fixed point of operator T .
At first, we prove the existence and uniqueness theorem of the Systems (1)–(2) by Banach contraction mapping principle.
Theorem 5.
Suppose that ( H 1 ) holds. If Ω 0 , and
Λ = Λ 1 C 3 ( C 7 + C 9 ) + Λ 2 C 6 ( C 8 + C 10 ) < 1 ,
where Λ 1 = max t [ 0 , 1 ] i = 1 4 L i ( t ) , Λ 2 = max t [ 0 , 1 ] i = 1 4 l i ( t ) . Then the Systems (1)–(2) has a unique solution.
Proof. 
Let r > 0 such that
r = C 0 ( C 7 + C 9 ) + C ˜ 0 ( C 8 + C 10 ) 1 Λ 1 C 3 ( C 7 + C 9 ) Λ 2 C 6 ( C 8 + C 10 ) ,
where C 0 = sup t [ 0 , 1 ] | P ( t , 0 , 0 , 0 , 0 ) | , C ˜ 0 = sup t [ 0 , 1 ] | Q ( t , 0 , 0 , 0 , 0 ) | .
We divide two steps to prove the theorem.
(i) Our first task is to show that T maps bounded sets into bounded sets in V.
Let B r = ( u , v ) V , ( u , v ) V r be a bounded set in V and ( u , v ) B r . Then we show that T ( B r ) B r . By ( H 1 ) and Lemma 3, we get
| F u v ( t ) | | P ( t , u ( t ) , v ( t ) , I q ω 1 u ( t ) , I q δ 1 v ( t ) ) P ( t , 0 , 0 , 0 , 0 ) | + | P ( t , 0 , 0 , 0 , 0 ) | L 1 ( t ) | u ( t ) | + L 2 ( t ) | v ( t ) | + L 3 ( t ) | I q ω 1 u ( t ) | + L 4 ( t ) | I q δ 1 v ( t ) | + C 0 Λ 1 u + v + u Γ q ( ω 1 ) + v Γ q ( δ 1 ) + C 0 = Λ 1 ( C 1 u + C 2 v ) + C 0 Λ 1 C 3 ( u , v ) V + C 0 Λ 1 C 3 r + C 0 ,
similarly,
| G u v ( t ) | Λ 2 C 6 r + C ˜ 0 .
According to the expression of operators T 1 and T 2 , we obtain
| T 1 ( u , v ) ( t ) | 1 Γ q ( α ) 0 t ( t q s ) ( α 1 ) | F u v ( s ) | d q s + Γ q ( β ) t α 1 | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) 0 1 ( 1 q s ) ( β ξ 0 1 ) | G u v ( s ) | d q s · | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) t α 1 | Ω | Γ q ( α ξ ) Γ q ( β ζ ) · | 0 1 0 s ( s q τ ) ( α ξ 1 ) | F u v ( τ ) | d q τ d q K ( s ) | · | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) t α 1 | Ω | Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) 0 1 ( 1 q s ) ( α ζ 0 1 ) | F u v ( s ) | d q s + Γ q ( β ) t α 1 | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) | 0 1 0 s ( s q τ ) ( β ζ 1 ) | G u v ( τ ) | d q τ d q H ( s ) | [ 1 Γ q ( α ) + Γ q ( β ) | Ω | Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) + Γ q ( β ) | Ω | Γ q ( α ξ ) Γ q ( β ζ ) · | 0 1 s β ζ 1 d q H ( s ) | · | 0 1 d q K ( s ) | ] 0 1 | F u v ( s ) | d q s + [ Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) · | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) | 0 1 d q H ( s ) | ] 0 1 | G u v ( s ) | d q s ,
thus, we have
T 1 ( u , v ) C 7 ( Λ 1 C 3 r + C 0 ) + C 10 ( Λ 2 C 6 r + C ˜ 0 ) ,
in like wise,
T 2 ( u , v ) C 9 ( Λ 1 C 3 r + C 0 ) + C 8 ( Λ 2 C 6 r + C ˜ 0 ) .
Using (5) and (6), we obtain that for ( u , v ) B r ,
T ( u , v ) V = T 1 ( u , v ) + T 2 ( u , v ) ( Λ 1 C 3 r + C 0 ) ( C 7 + C 9 ) + ( Λ 2 C 6 r + C ˜ 0 ) ( C 8 + C 10 ) = r ,
that is T ( B r ) B r .
(ii) The next step is to prove that operator T is a contraction.
For ( u i , v i ) B r ( i = 1 , 2 ) , t [ 0 , 1 ] , we get
| T 1 ( u 1 , v 1 ) ( t ) T 1 ( u 2 , v 2 ) ( t ) | 1 Γ q ( α ) 0 t ( t q s ) ( α 1 ) | F u 1 v 1 ( s ) F u 2 v 2 ( s ) | d q s + Γ q ( β ) t α 1 | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) 0 1 ( 1 q s ) ( β ξ 0 1 ) | G u 1 v 1 ( s ) G u 2 v 2 ( s ) | d q s · | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) t α 1 | Ω | Γ q ( α ξ ) Γ q ( β ζ ) | 0 1 0 s ( s q τ ) ( α ξ 1 ) | F u 1 v 1 ( τ ) F u 2 v 2 ( τ ) | d q τ d q K ( s ) | · | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) t α 1 | Ω | Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) 0 1 ( 1 q s ) ( α ζ 0 1 ) · | F u 1 v 1 ( s ) F u 2 v 2 ( s ) | d q s + Γ q ( β ) t α 1 | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) · | 0 1 0 s ( s q τ ) ( β ζ 1 ) | G u 1 v 1 ( τ ) G u 2 v 2 ( τ ) | d q τ d q H ( s ) | .
Since
| F u 1 v 1 ( s ) F u 2 v 2 ( s ) | [ L 1 ( s ) | u 1 ( s ) u 2 ( s ) | + L 2 ( s ) | v 1 ( s ) v 2 ( s ) | + L 3 ( s ) | I q ω 1 u 1 ( s ) I q ω 1 u 2 ( s ) | + L 4 ( s ) | I q δ 1 v 1 ( s ) I q δ 1 v 2 ( s ) | ] Λ 1 ( C 1 u 1 u 2 + C 2 v 1 v 2 ) Λ 1 C 3 ( u 1 , v 1 ) ( u 2 , v 2 ) V ,
and
| G u 1 v 1 ( s ) G u 2 v 2 ( s ) | Λ 2 C 6 ( u 1 , v 1 ) ( u 2 , v 2 ) V .
By (7), we have
| T 1 ( u 1 , v 1 ) ( t ) T 1 ( u 2 , v 2 ) ( t ) | 1 Γ q ( α ) 0 1 | F u 1 v 1 ( s ) F u 2 v 2 ( s ) | d q s + Γ q ( β ) | Ω | Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) 0 1 | F u 1 v 1 ( s ) F u 2 v 2 ( s ) | d q s + Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) 0 1 | G u 1 v 1 ( s ) G u 2 v 2 ( s ) | d q s · | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) | Ω | Γ q ( α ξ ) Γ q ( β ζ ) | 0 1 0 1 | F u 1 v 1 ( τ ) F u 2 v 2 ( τ ) | d q τ d q K ( s ) | · | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) · | 0 1 0 1 | G u 1 v 1 ( τ ) G u 2 v 2 ( τ ) | d q ( τ ) d q H ( s ) | Λ 1 C 3 ( u 1 , v 1 ) ( u 2 , v 2 ) V [ 1 Γ q ( α ) + Γ q ( β ) | Ω | Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) + Γ q ( β ) | Ω | Γ q ( α ξ ) Γ q ( β ζ ) | 0 1 d q H ( s ) | | 0 1 s β ζ 1 d q H ( s ) | ] + Λ 2 C 6 ( u 1 , v 1 ) ( u 2 , v 2 ) V [ Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) | 0 1 d q H ( s ) | ] ,
hence, we deduce
T 1 ( u 1 , v 1 ) T 1 ( u 2 , v 2 ) ( Λ 1 C 3 C 7 + Λ 2 C 6 C 10 ) ( u 1 , v 1 ) ( u 2 , v 2 ) V .
For the same way, we can obtain
T 2 ( u 1 , v 1 ) T 2 ( u 2 , v 2 ) ( Λ 1 C 3 C 9 + Λ 2 C 6 C 8 ) ( u 1 , v 1 ) ( u 2 , v 2 ) V .
From (8) and (9), we have
T ( u 1 , v 1 ) T ( u 2 , v 2 ) V = T 1 ( u 1 , v 1 ) T 1 ( u 2 , v 2 ) + T 2 ( u 1 , v 1 ) T 2 ( u 2 , v 2 ) Λ 1 C 3 ( C 7 + C 9 ) + Λ 2 C 6 ( C 8 + C 10 ) ( u 1 , v 1 ) ( u 2 , v 2 ) V = Λ ( u 1 , v 1 ) ( u 2 , v 2 ) V .
Due to Λ < 1 , it follows that T ( u 1 , v 1 ) T ( u 2 , v 2 ) V < ( u 1 , v 1 ) ( u 2 , v 2 ) V , so operator T is a contraction. Hence, we obtain that the Systems (1)–(2) has a unique solution ( u , v ) B r by using Banach contraction mapping principle. The proof is completed. □
Corollary 1.
Suppose that ( H 1 ) holds. If Ω 0 , and
Λ * = Λ 3 C 3 ( C 7 + C 9 ) + Λ 4 C 6 ( C 8 + C 10 ) < 1 ,
where Λ 3 = i = 1 4 L i , Λ 4 = i = 1 4 l i . Then the Systems (1)–(2) has a unique solution.
Corollary 2.
Suppose that ( H 1 ) holds. If Ω 0 , and
Λ ˜ = Λ 5 C 3 ( C 7 + C 9 ) + Λ 6 C 6 ( C 8 + C 10 ) < 1 ,
where Λ 5 = sup t [ 0 , 1 ] { i = 1 4 ρ i ( t ) } , Λ 6 = sup t [ 0 , 1 ] { i = 1 4 ϱ i ( t ) } . Then the Systems (1)–(2) has a unique solution.
Next, we apply several kinds of fixed point theorems to achieve the existence results of solutions for the Systems (1)–(2).
Theorem 6.
Suppose that ( H 2 ) and Ω 0 hold. Then the System (1)–(2) has at least one solution.
Proof. 
Let B R = { ( u , v ) V , ( u , v ) V R } , and we denote
R 1 = max { [ c 0 + c 1 ( N 1 ) h 1 + c 2 ( N 2 ) h 2 + c 3 N 1 Γ q ( ω 1 ) h 3 + c 4 N 2 Γ q ( δ 1 ) h 4 ] C 7 , [ d 0 + d 1 ( N 1 ) m 1 + d 2 ( N 2 ) m 2 + d 3 N 1 Γ q ( ω 2 ) m 3 + d 4 N 2 Γ q ( δ 2 ) m 4 ] C 10 } , R 2 = max { [ c 0 + c 1 ( N 1 ) h 1 + c 2 ( N 2 ) h 2 + c 3 N 1 Γ q ( ω 1 ) h 3 + c 4 N 2 Γ q ( δ 1 ) h 4 ] C 9 , [ d 0 + d 1 ( N 1 ) m 1 + d 2 ( N 2 ) m 2 + d 3 N 1 Γ q ( ω 2 ) m 3 + d 4 N 2 Γ q ( δ 2 ) m 4 ] C 8 } , R = 2 max R 1 , R 2 .
where there exist N 1 , N 2 R such that | u ( t ) | N 1 , | v ( t ) | N 2 .
Firstly, we show that T maps bounded sets into bounded sets in V. For ( u , v ) B R , we obtain
T 1 ( u , v ) [ c 0 + c 1 ( N 1 ) h 1 + c 2 ( N 2 ) h 2 + c 3 N 1 Γ q ( ω 1 ) h 3 + c 4 N 2 Γ q ( δ 1 ) h 4 ] C 7 + [ d 0 + d 1 ( N 1 ) m 1 + d 2 ( N 2 ) m 2 + d 3 N 1 Γ q ( ω 2 ) m 3 + d 4 N 2 Γ q ( δ 2 ) m 4 ] C 10 2 R 1 ,
similarly, T 2 ( u , v ) 2 R 2 , then
T ( u , v ) V = T 1 ( u , v ) + T 2 ( u , v ) R , ( u , v ) B R ,
as above, we obtain T ( B R ) B R .
Secondly, we prove that T maps bounded sets into equicontinuous sets of V. Let N = max { N 1 , N 2 } , for simplicity of presentation, we denote that
Ψ N = sup t [ 0 , 1 ] | P ( t , u , v , x , y ) | , | u | N , | v | N , | x | N Γ q ( ω 1 ) , | y | N Γ q ( δ 1 ) , Θ N = sup t [ 0 , 1 ] | Q ( t , u , v , x , y ) | , | u | N , | v | N , | x | N Γ q ( ω 2 ) , | y | N Γ q ( δ 2 ) ,
then for ( u , v ) B R and t 1 , t 2 [ 0 , 1 ] with t 1 < t 2 , we have
| T 1 ( u , v ) ( t 2 ) T 1 ( u , v ) ( t 1 ) | Ψ N Γ q ( α + 1 ) ( t 2 α t 1 α ) + Ψ N ( t 2 α 1 t 1 α 1 ) [ Γ q ( β ) | Ω | Γ q ( α ξ ) Γ q ( β ζ ) | 0 1 s β ζ 1 d q H ( s ) | · | 0 1 d q K ( s ) | + Γ q ( β ) | Ω | Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) ] + Θ N ( t 2 α 1 t 1 α 1 ) [ Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) · | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) | 0 1 d q H ( s ) | ] = Ψ N Γ q ( α + 1 ) ( t 2 α t 1 α ) + ( Ψ N C 11 + Θ N C 10 ) ( t 2 α 1 t 1 α 1 ) .
The same can be proved that
| T 2 ( u , v ) ( t 2 ) T 2 ( u , v ) ( t 1 ) | Θ N Γ q ( β + 1 ) ( t 2 β t 1 β ) + ( Ψ N C 9 + Θ N C 12 ) ( t 2 β 1 t 1 β 1 ) .
Hence, we conclude
| T 1 ( u , v ) ( t 2 ) T 1 ( u , v ) ( t 1 ) | 0 , | T 2 ( u , v ) ( t 2 ) T 2 ( u , v ) ( t 1 ) | 0 ,
as t 2 t 1 , ( u , v ) B R . Thus, T ( B R ) is equicontinuous. According to the Arzela-Ascoli theorem, it follows that the set T ( B R ) is relatively compact. Therefore, T is compact on B R . By Theorem 1, we get that the System (1)–(2) has at least one solution. The proof is completed. □
Theorem 7.
Suppose that ( H 1 ) and ( H 3 ) hold. If Ω 0 , and
Λ ¯ = Λ 3 C 3 1 Γ q ( α ) + Λ 4 C 6 1 Γ q ( β ) < 1 .
Then the System (1)–(2) has at least one solution.
Proof. 
Take r 0 > 0 such that
r 0 ( C 7 + C 9 ) σ 1 + ( C 8 + C 10 ) σ 2 .
Let B r 0 = { ( u , v ) V , ( u , v ) V r 0 } , and let the operators be X = ( X 1 , X 2 ) : B r 0 V and Y = ( Y 1 , Y 2 ) : B r 0 V , where X 1 , X 2 , Y 1 , Y 2 : B r 0 U are denoted by
X 1 ( u , v ) ( t ) = 1 Γ q ( α ) 0 t ( t q s ) ( α 1 ) F u v ( s ) d q s , Y 1 ( u , v ) ( t ) = Ω 1 t α 1 Ω Γ q ( β ξ 0 ) 0 1 ( 1 q s ) ( β ξ 0 1 ) G u v ( s ) d q s Ω 1 t α 1 Ω Γ q ( α ξ ) 0 1 0 s ( s q τ ) ( α ξ 1 ) F u v ( τ ) d q τ d q K ( s ) + Γ q ( β ) t α 1 Ω Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) 0 1 ( 1 q s ) ( α ζ 0 1 ) F u v ( s ) d q s Γ q ( β ) t α 1 Ω Γ q ( β ζ ) Γ q ( β ξ 0 ) 0 1 0 s ( s q τ ) ( β ζ 1 ) G u v ( τ ) d q τ d q H ( s ) , X 2 ( u , v ) ( t ) = 1 Γ q ( β ) 0 t ( t q s ) ( β 1 ) G u v ( s ) d q s , Y 2 ( u , v ) ( t ) = Ω 2 t β 1 Ω Γ q ( α ζ 0 ) 0 1 ( 1 q s ) ( α ζ 0 1 ) F u v ( s ) d q s Ω 2 t β 1 Ω Γ q ( β ζ ) 0 1 0 s ( s q τ ) ( β ζ 1 ) G u v ( τ ) d q τ d q H ( s ) + Γ q ( α ) t β 1 Ω Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) 0 1 ( 1 q s ) ( β ξ 0 1 ) G u v ( s ) d q s Γ q ( α ) t β 1 Ω Γ q ( α ξ ) Γ q ( α ζ 0 ) 0 1 0 s ( s q τ ) ( α ξ 1 ) F u v ( τ ) d q τ d q K ( s ) ,
where t [ 0 , 1 ] , ( u , v ) B r 0 . Thus, T 1 = X 1 + Y 1 , T 2 = X 2 + Y 2 and T = X + Y .
By ( H 3 ) , we know that ( u 1 , v 1 ) , ( u 2 , v 2 ) B r 0 ,
X ( u 1 , v 1 ) + Y ( u 2 , v 2 ) V X ( u 1 , v 1 ) V + Y ( u 2 , v 2 ) V = X 1 ( u 1 , v 1 ) + X 2 ( u 1 , v 1 ) + Y 1 ( u 2 , v 2 ) + Y 2 ( u 2 , v 2 ) 1 Γ q ( α ) σ 1 + 1 Γ q ( β ) σ 2 + σ 1 [ Γ q ( β ) | Ω | Γ q ( α ξ ) Γ q ( β ζ ) | 0 1 d q K ( s ) | · | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) | Ω | Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) ] + σ 2 [ Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) · | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) | 0 1 d q H ( s ) | ] + σ 1 · Γ q ( α ) | Ω | Γ q ( α ζ 0 ) Γ q ( α ξ ) | 0 1 s α ξ 1 d q K ( s ) | + Γ q ( α ) | Ω | Γ q ( α ξ ) Γ q ( α ζ 0 ) | 0 1 d q K ( s ) | + σ 2 Γ q ( α ) | Ω | Γ q ( α ξ ) Γ q ( β ζ ) | 0 1 d q H ( s ) | + Γ q ( α ) | Ω | Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) = ( C 7 + C 9 ) σ 1 + ( C 8 + C 10 ) σ 2 r 0 .
For ( u 1 , v 1 ) , ( u 2 , v 2 ) B r 0 , and Λ ¯ < 1 , we have
X ( u 1 , v 1 ) X ( u 2 , v 2 ) V = X 1 ( u 1 , v 1 ) X 1 ( u 2 , v 2 ) + X 2 ( u 1 , v 1 ) X 2 ( u 2 , v 2 ) ( Λ 3 C 3 1 Γ q ( α ) + Λ 4 C 6 1 Γ q ( β ) ) ( u 1 u 2 + v 1 v 2 ) = Λ ¯ ( u 1 , v 1 ) ( u 2 , v 2 ) V < ( u 1 , v 1 ) ( u 2 , v 2 ) V .
Hence, the operator X is a contraction.
Owing to the continuity of P and Q, Y is continuous. Next, we need to verify that Y is a compact operator. Due to ( u , v ) B r 0 ,
Y ( u , v ) V = Y 1 ( u , v ) + Y 2 ( u , v ) ( C 9 + C 11 ) σ 1 + ( C 10 + C 12 ) σ 2 ,
we have derived that the functions from Y are uniformly bounded.
We can show the equicontinuous of the functions from Y ( B r 0 ) . We denote that
Ψ r 0 = sup t [ 0 , 1 ] | P ( t , u , v , x , y ) | , | u | r 0 , | v | r 0 , | x | r 0 Γ q ( ω 1 ) , | y | r 0 Γ q ( δ 1 ) , Θ r 0 = sup t [ 0 , 1 ] | Q ( t , u , v , x , y ) | , | u | r 0 , | v | r 0 , | x | r 0 Γ q ( ω 2 ) , | y | r 0 Γ q ( δ 2 ) ,
for ( u , v ) B r 0 and t 1 , t 2 [ 0 , 1 ] with t 1 < t 2 . An argument similar to the one used in the proof of Theorem 6 shows that
| Y 1 ( u , v ) ( t 2 ) Y 1 ( u , v ) ( t 1 ) | 0 , | Y 2 ( u , v ) ( t 2 ) Y 2 ( u , v ) ( t 1 ) | 0 ,
as t 2 t 1 , ( u , v ) B r 0 . Therefore, Y ( B r 0 ) is equicontinuous. Then, we can see that Y ( B r 0 ) is relatively compact. Hence, Y is compact on B r 0 . Using Theorem 2, we know that the System (1)–(2) has at least one solution. The proof is completed. □
Remark 1.
Evidently, we prove that the operator X is a contraction, the operator Y is compact and continuous in Theorem 7. An alternative method of proof is to show that X is compact and continuous, Y is a contraction, that is Theorem 8.
Theorem 8.
Suppose that ( H 1 ) and ( H 3 ) hold. If Ω 0 , and
Λ ^ = Λ 3 C 3 ( C 9 + C 11 ) + Λ 4 C 6 ( C 8 + C 10 ) < 1 .
Then the Systems (1)–(2) has at least one solution.
Proof. 
On the basis of Remark 1, this theorem can be proved by the same method as employed in Theorem 7. □
Theorem 9.
Suppose that P , Q are S-Carathéodory functions and ( H 4 ) hold. If Ω 0 , and
Ξ = max { C 13 , C 14 } < 1 ,
where C 13 = ( j 1 + j 3 Γ q ( ω 1 ) ) ( C 7 + C 9 ) + ( k 1 + k 3 Γ q ( ω 2 ) ) ( C 8 + C 10 ) , C 14 = ( j 2 + j 4 Γ q ( δ 1 ) ) ( C 7 + C 9 ) + ( k 2 + k 4 Γ 1 ( δ 2 ) ) ( C 8 + C 10 ) , and there exist A 1 , A 2 , j i , k i > 0 such that | r 1 ( t ) | A 1 ,   | r 2 ( t ) | A 2 , | L i ( t ) | j i and | l i ( t ) | k i ( i = 1 , 2 , 3 , 4 ) . Then, the System (1)–(2) has at least one solution.
Proof. 
The main point of Theorem 9 is to prove T is completely continuous. Firstly, for the continuity of functions P and Q, we obtain that the operator T is continuous. Secondly, we show that T is compact.
Let the set Φ V be bounded. Then, there exist integrable functions M 1 ( t ) and M 2 ( t ) L 1 ( [ 0 , 1 ] , R + ) such that for t [ 0 , 1 ] , ( u , v ) Φ , we have
| P ( t , u ( t ) , v ( t ) , I q ω 1 u ( t ) , I q δ 1 v ( t ) ) | M 1 ( t ) , | Q ( t , u ( t ) , v ( t ) , I q ω 2 u ( t ) , I q δ 2 v ( t ) ) | M 2 ( t ) .
According to the Theorem 5, we get
| F u v ( t ) | = | P ( t , u ( t ) , v ( t ) , I q ω 1 u ( t ) , I q δ 1 v ( t ) ) | M 1 L 1 , | G u v ( t ) | = | Q ( t , u ( t ) , v ( t ) , I q ω 2 u ( t ) , I q δ 2 v ( t ) ) | M 2 L 1 ,
where u L 1 = 0 1 | u ( t ) | d q t .
Then
T 1 ( u , v ) M 1 L 1 [ 1 Γ q ( α ) + Γ q ( β ) | Ω | Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) + Γ q ( β ) | Ω | Γ q ( α ξ ) Γ q ( β ζ ) | 0 1 s β ζ 1 d q H ( s ) | · | 0 1 d q K ( s ) | ] + M 2 L 1 [ Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) | 0 1 d q H ( s ) | ] = M 1 L 1 C 7 + M 2 L 1 C 10 ,
in a similar manner, we have
T 2 ( u , v ) M 1 L 1 C 9 + M 2 L 1 C 8 ,
so ( u , v ) Φ ,
T ( u , v ) V = T 1 ( u , v ) + T 2 ( u , v ) M 1 L 1 ( C 7 + C 9 ) + M 2 L 1 ( C 10 + C 8 ) ,
therefore, T ( Φ ) is uniformly bounded.
Another step is to show that T ( Φ ) is equicontinuous. Proceeding as in the proof of Theorem 6, we obtain | T 1 ( u , v ) ( t 2 ) T 1 ( u , v ) ( t 1 ) | 0 and | T 2 ( u , v ) ( t 2 ) T 2 ( u , v ) ( t 1 ) | 0 , as t 2 t 1 , ( u , v ) Φ . Thus, T ( Φ ) is equicontinuous. At the same time, we can also obtain that T is completely continuous.
Finally, we illustrate that S = { ( u , v ) V , ( u , v ) = λ T ( u , v ) , 0 λ 1 } is bounded. Let ( u , v ) S , then t [ 0 , 1 ] , we have u ( t ) = λ T 1 ( u , v ) ( t ) , v ( t ) = λ T 2 ( u , v ) ( t ) . For simplicity, we denote that
F ^ u v ( s ) = r 1 ( s ) + L 1 ( s ) | u ( s ) | + L 2 ( s ) | v ( s ) | + L 3 ( s ) | I q ω 1 u ( s ) | + L 4 ( s ) | I q δ 1 v ( s ) | , G ^ u v ( s ) = r 2 ( s ) + l 1 ( s ) | u ( s ) | + l 2 ( s ) | v ( s ) | + l 3 ( s ) | I q ω 2 u ( s ) | + l 4 ( s ) | I q δ 2 v ( s ) | ,
so,
F ^ u v ( s ) A 1 + j 1 | u ( s ) | + j 2 | v ( s ) | + j 3 | I q ω 1 u ( s ) | + j 4 | I q δ 1 v ( s ) | , G ^ u v ( s ) A 2 + k 1 | u ( s ) | + k 2 | v ( s ) | + k 3 | I q ω 2 u ( s ) | + k 4 | I q δ 2 v ( s ) | ,
then
| u ( t ) | | T 1 ( u , v ) ( t ) | 1 Γ q ( α ) 0 t ( t q s ) ( α 1 ) F ^ u v ( s ) d q s + Γ q ( β ) t α 1 | Ω | Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) 0 1 ( 1 q s ) ( α ζ 0 1 ) · F ^ u v ( s ) d q s + Γ q ( β ) t α 1 | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) | 0 1 s β ζ 1 d q H ( s ) | 0 1 ( 1 q s ) ( β ξ 0 1 ) · G ^ u v ( s ) d q s + Γ q ( β ) t α 1 | Ω | Γ q ( α ξ ) Γ q ( β ζ ) | 0 1 s β ζ 1 d q H ( s ) | | 0 1 [ 0 s ( s q τ ) ( α ξ 1 ) · F ^ u v ( τ ) d q τ ] d q K ( s ) | + Γ q ( β ) t α 1 | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) | 0 1 [ 0 s ( s q τ ) ( β ζ 1 ) · G ^ u v ( τ ) d q τ ] d q H ( s ) | ,
hence,
u 0 1 | F ^ u v ( s ) | d q s [ 1 Γ q ( α ) + Γ q ( β ) | Ω | Γ q ( α ξ ) Γ q ( β ζ ) | 0 1 s β ζ 1 d q H ( s ) | · | 0 1 d q K ( s ) | + Γ a ( β ) | Ω | Γ q ( α ζ 0 ) Γ q ( β ξ 0 ) ] + 0 1 | G ^ u v ( s ) | d q s · [ Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) · | 0 1 s β ζ 1 d q H ( s ) | + Γ q ( β ) | Ω | Γ q ( β ξ 0 ) Γ q ( β ζ ) | 0 1 d q H ( s ) | ] = ( A 1 + j 1 u + j 2 v + j 3 Γ q ( ω 1 ) u + j 4 Γ q ( δ 1 ) v ) C 7 + ( A 2 + k 1 u + k 2 v + k 3 Γ q ( ω 2 ) u + k 4 Γ q ( δ 2 ) v ) C 10 .
Similarly,
v ( A 1 + j 1 u + j 2 v + j 3 Γ q ( ω 1 ) u + j 4 Γ q ( δ 1 ) v ) C 9 + ( A 2 + k 1 u + k 2 v + k 3 Γ q ( ω 2 ) u + k 4 Γ q ( δ 2 ) v ) C 8 ,
by means of (11) and (12), we have
( u , v ) V = u + v A 1 ( C 7 + C 9 ) + A 2 ( C 8 + C 10 ) + C 13 u + C 14 v A 1 ( C 7 + C 9 ) + A 2 ( C 8 + C 10 ) + Ξ ( u , v ) V .
Due to Ξ < 1 , we get
( u , v ) V A 1 ( C 7 + C 9 ) + A 2 ( C 8 + C 10 ) ( 1 Ξ ) 1 , ( u , v ) S ,
thus, S is bounded.
By Theorem 3, it is time to say that the Systems (1)–(2) has at least one solution. Hence, the statements in Theorem 9 are proved. □
Corollary 3.
Suppose that P , Q are S-Carathéodory functions and ( H 4 ) hold. If Ω 0 , and
Ξ ^ = max { C 15 , C 16 } < 1 ,
where C 15 = L 1 ( C 7 + C 9 ) + L 3 Γ q ( ω 1 ) ( C 7 + C 9 ) + l 1 ( C 8 + C 10 ) + l 3 Γ q ( ω 2 ) ( C 8 + C 10 ) , C 16 = L 2 ( C 7 + C 9 ) + L 4 Γ q ( δ 1 ) ( C 7 + C 9 ) + l 2 ( C 8 + C 10 ) + l 4 Γ q ( δ 2 ) ( C 8 + C 10 ) . Then the Systems (1)–(2) has at least one solution.
Theorem 10.
Suppose that P , Q are S-Carathéodory functions and ( H 5 ) hold. If Ω 0 and there exists Π > 0 such that
[ p 0 + p 1 φ 1 ( Π ) + p 2 φ 2 ( Π ) + p 3 φ 3 Π Γ q ( ω 1 ) + p 4 φ 4 Π Γ q ( δ 1 ) ] ( C 7 + C 9 ) + [ q 0 + q 1 η 1 ( Π ) + q 2 η 2 ( Π ) + q 3 η 3 Π Γ q ( ω 2 ) + q 4 η 4 Π Γ q ( δ 2 ) ] ( C 8 + C 10 ) < Π .
Then the Systems (1)–(2) has at least one solution.
Proof. 
Let B Π = { ( u , v ) V , ( u , v ) V Π } . Firstly, we prove that T : B Π B Π . For ( u , v ) B Π and t [ 0 , 1 ] , we have
T 1 ( u , v ) C 7 [ p 0 + p 1 φ 1 ( Π ) + p 2 φ 2 ( Π ) + p 3 φ 3 Π Γ q ( ω 1 ) + p 4 φ 4 Π Γ q ( δ 1 ) ] + C 10 [ q 0 + q 1 η 1 ( Π ) + q 2 η 2 ( Π ) + q 3 η 3 Π Γ q ( ω 2 ) + q 4 η 4 Π Γ q ( δ 2 ) ] ,
and
T 2 ( u , v ) C 9 [ p 0 + p 1 φ 1 ( Π ) + p 2 φ 2 ( Π ) + p 3 φ 3 Π Γ q ( ω 1 ) + p 4 φ 4 Π Γ q ( δ 1 ) ] + C 8 [ q 0 + q 1 η 1 ( Π ) + q 2 η 2 ( Π ) + q 3 η 3 Π Γ q ( ω 2 ) + q 4 η 4 Π Γ q ( δ 2 ) ] .
For ( u , v ) B Π , we have
T ( u , v ) V = T 1 ( u , v ) + T 2 ( u , v ) ( C 7 + C 9 ) [ p 0 + p 1 φ 1 ( Π ) + p 2 φ 2 ( Π ) + p 3 φ 3 Π Γ q ( ω 1 ) + p 4 φ 4 Π Γ q ( δ 1 ) ] + ( C 10 + C 8 ) [ q 0 + q 1 η 1 ( Π ) + q 2 η 2 ( Π ) + q 3 η 3 Π Γ q ( ω 2 ) + q 4 η 4 Π Γ q ( δ 2 ) ] < Π .
Consequently, T ( B Π ) B Π . At the same time, it is easy to see that T is completely continuous, which can be derived in the same way as employed in Theorem 6.
Furthermore, assume that there exists ( u , v ) B Π such that ( u , v ) = λ T ( u , v ) for λ ( 0 , 1 ) , it is simple to get ( u , v ) V T ( u , v ) V < Π , this leads to a contradiction for ( u , v ) B Π . Therefore, by applying Theorem 4, we deduce that T has a fixed point ( u , v ) B Π , which is a solution of the Systems (1)–(2). The proof is completed. □

4. Application Examples

In this section, for the system with the different nonlinearity terms, some examples are appreciated to illustrate our main results.
We consider the following system of fractional q-difference equations:
( D q 3 2 u ) ( t ) + P ( t , u ( t ) , v ( t ) , I q 1 4 u ( t ) , I q 4 3 v ( t ) ) = 0 , t ( 0 , 1 ) , ( D q 5 2 v ) ( t ) + Q ( t , u ( t ) , v ( t ) , I q 9 4 u ( t ) , I q 2 3 v ( t ) ) = 0 , t ( 0 , 1 ) ,
with the nonlocal BCs
u ( 0 ) = 0 , D q 1 5 u ( 1 ) = 0 1 D q 5 4 v ( t ) d q ( t 2 ) , v ( 0 ) = D q v ( 0 ) = 0 , D q 7 5 v ( 1 ) = 0 1 D q 1 6 u ( t ) d q t ,
where α = 3 2 , β = 5 2 , ω 1 = 1 4 , δ 1 = 4 3 , ω 2 = 9 4 , δ 2 = 2 3 , ζ 0 = 1 5 , ξ 0 = 7 5 , ζ = 5 4 , ξ = 1 6 ,   q = 1 2 , H ( t ) = t 2 , K ( t ) = t .
After a simple caculation, we obtain Ω = 2.58954375 0 , C 1 = 1.34100597 ,   C 2 = 2.08201688 , C 3 = C 2 , C 4 = 1.92455621 , C 5 = 1.79251862 , C 6 = C 4 ,   C 7 = 2.29230629 , C 8 = 1.75022309 , C 9 = 0.7590784 , C 10 = 1.42695841 , C 11 = 1.20641206 ,   C 12 = 0.91031044 .
Example 1.
Consider the nonlinear terms of the system
P ( t , x 1 , x 2 , x 3 , x 4 ) = e t + t 36 cos x 1 t 54 sin x 2 + 1 63 + t arctan x 3 x 4 ( t + 9 ) 2 , Q ( t , x 1 , x 2 , x 3 , x 4 ) = 1 5 + t 2 t 48 sin x 1 + t 64 cos x 2 1 36 + t arctan x 3 + x 4 t 2 + 56 ,
where t [ 0 , 1 ] , x i R ( i = 1 , 2 , 3 , 4 ) . For x i , y i R ( i = 1 , 2 , 3 , 4 ) , we obtain
| P ( t , x 1 , x 2 , x 3 , x 4 ) P ( t , y 1 , y 2 , y 3 , y 4 ) | t 36 | x 1 y 1 | + t 54 | x 2 y 2 | + 1 t + 63 | x 3 y 3 | + 1 ( t + 9 ) 2 | x 4 y 4 | Λ 1 i = 1 4 | x i y i | , | Q ( t , x 1 , x 2 , x 3 , x 4 ) Q ( t , y 1 , y 2 , y 3 , y 4 ) | t 48 | x 1 y 1 | + t 64 | x 2 y 2 | + 1 t + 36 | x 3 y 3 | + 1 56 + t 2 | x 4 y 4 | Λ 2 i = 1 4 | x i y i | .
It is obvious that L 1 ( t ) = t 36 , L 2 ( t ) = t 54 , L 3 ( t ) = 1 63 + t , L 4 ( t ) = 1 ( t + 9 ) 2 and l 1 ( t ) = t 48 , l 2 ( t ) = t 64 , l 3 ( t ) = 1 36 + t , l 4 ( t ) = 1 t 2 + 56 . By a simple computation, we obtain Λ 1 = 0.07451499 , Λ 2 = 0.08209325 and Λ = 0.97536897 < 1 , respectively. By Theorem 5, the Systems (12)–(13) has a unique solution.
Example 2.
Consider the nonlinear terms of the system
P ( t , x 1 , x 2 , x 3 , x 4 ) = ( 3 t + 5 ) 2 + t 30 | x 1 | 1 3 + t t 2 + 9 arctan | x 2 | 1 2 + 1 5 + 8 t sin | x 3 | 3 4
1 ( 4 t + 9 ) 2 | x 4 | 1 5 , Q ( t , x 1 , x 2 , x 3 , x 4 ) = e t + 1 46 ( t + 1 ) | x 1 | 3 5 + t 37 | x 2 | 1 6 + 1 t + 29 sin | x 3 | 2 3 8 t arctan | x 4 | 5 6 ,
where t [ 0 , 1 ] , x i R ( i = 1 , 2 , 3 , 4 ) . It is clear that
| P ( t , x 1 , x 2 , x 3 , x 4 ) | ( 3 t + 5 ) 2 + t 30 | x 1 | 1 3 + t t 2 + 9 | x 2 | 1 2 + 1 8 t + 5 | x 3 | 3 4 + 1 ( 4 t + 9 ) 2 | x 4 | 1 5 , | Q ( t , x 1 , x 2 , x 3 , x 4 ) | e t + 1 46 ( t + 1 ) | x 1 | 3 5 + t 37 | x 2 | 1 6 + 1 t + 29 | x 3 | 2 3 + 8 t | x 4 | 5 6 .
Therefore, the assumption ( H 2 ) is satisfied with c 0 ( t ) = ( 3 t + 5 ) 2 , c 1 ( t ) = t 30 ,   c 2 ( t ) = t t 2 + 9 ,   c 3 ( t ) = 1 8 t + 5 , c 4 ( t ) = 1 ( 4 t + 9 ) 2 , d 0 ( t ) = e t , d 1 ( t ) = 1 46 ( t + 1 ) , d 2 ( t ) = t 37 ,   d 3 ( t ) = 1 t + 29 , and d 4 ( t ) = 8 t . By Theorem 6, the Systems (12)–(13) has at least one solution.
Example 3.
Consider the nonlinear terms of the system
P ( t , x 1 , x 2 , x 3 , x 4 ) = e t + t 40 arctan x 1 1 ( t + 6 ) 2 sin x 2 + 1 4 ( t + 9 ) sin x 3 t 32 cos x 4 , Q ( t , x 1 , x 2 , x 3 , x 4 ) = 5 t 6 + t 2 3 t 56 cos x 1 + 1 t + 28 sin x 2 t 72 sin 2 x 3 + t 18 arctan x 4 ,
where t [ 0 , 1 ] , x i R ( i = 1 , 2 , 3 , 4 ) . For t [ 0 , 1 ] , x i , y i R ( i = 1 , 2 , 3 , 4 ) , We obtain
| P ( t , x 1 , x 2 , x 3 , x 4 ) P ( t , y 1 , y 2 , y 3 , y 4 ) | 1 40 | x 1 y 1 | + 1 36 | x 2 y 2 | + 1 36 | x 3 y 3 | + 1 32 | x 4 y 4 | , | Q ( t , x 1 , x 2 , x 3 , x 4 ) Q ( t , y 1 , y 2 , y 3 , y 4 ) | 3 56 | x 1 y 1 | + 1 28 | x 2 y 2 | + 1 36 | x 3 y 3 | + 1 18 | x 4 y 4 | ,
and
| P ( t , x 1 , x 2 , x 3 , x 4 ) | e t + π t 80 + 1 ( t + 6 ) 2 + 1 4 ( t + 9 ) + t 32 , | Q ( t , x 1 , x 2 , x 3 , x 4 ) | 5 t 6 + t 2 + 3 t 56 + 1 t + 28 + t 72 + π t 36 ,
It is obvious that L 1 = 1 40 , L 2 = 1 36 , L 3 = 1 36 , L 4 = 1 32 , l 1 = 3 56 , l 2 = 1 28 ,   l 3 = 1 36 ,   l 4 = 1 18 . By a simple computation, we have Λ 3 = 0.11180556 , Λ 4 = 0.17261905 , and Λ ¯ = 0.53180725 < 1 , respectively. Therefore, the assumptions ( H 1 ) , and ( H 3 ) are satisfied, by Theorem 7, the Systems (12)–(13) has at least one solution.
Example 4.
Consider the nonlinear terms of the system
P ( t , x 1 , x 2 , x 3 , x 4 ) = t 20 + 6 t 35 sin x 1 1 4 t + 9 6 sin 2 x 2 , Q ( t , x 1 , x 2 , x 3 , x 4 ) = 2 t 3 + 5 t 56 sin x 1 3 2 ( 7 5 + t ) sin 2 x 2 ,
where t [ 0 , 1 ] , x i R ( i = 1 , 2 , 3 , 4 ) . It is clear that
| P ( t , x 1 , x 2 , x 3 , x 4 ) | 1 20 + 6 35 | x 1 | + 2 9 6 | x 2 | ,
| Q ( t , x 1 , x 2 , x 3 , x 4 ) | 2 3 + 5 56 | x 1 | + 3 7 5 | x 2 | .
Hence, L 1 = 6 35 , L 2 = 2 9 6 , L 3 = L 4 = 0 , r 1 = 1 20 , l 1 = 5 56 , l 2 = 3 7 5 , l 3 = l 4 = 0 , r 2 = 2 3 . By a simple computation, we obtain C 15 = 0.80677144 , C 16 = 0.88577528 , and Ξ ^ = 0.88577528 < 1 , respectively. By Corollary 3, the Systems (12)–(13) has at least one solution.

5. Discussion

The system of fractional q-difference equations plays an extremely crucial role in many fields, such as quantum mechanics, dynamical systems, black holes, mathematical physics equations and so on, see [2,3,5,6,27,28,29,30] and the references therein. In this article, we are concerned with the solvability of a system of fractional q-difference equations with Riemann-Stieltjes integrals conditions based on some classical fixed point theorems. We obtain the multiple existence and uniqueness conclusions for the Systems (1)–(2). As a matter of fact, in the limit q 1 , the system studied in this paper reduces to the classical system of fractional differential equations. It follows that the results we have discussed are the generalization of the classical analysis, they can extend classical theory in order to expand the range of the possible applications. In the future, we will devote ourselves to finding new inspirations and outstanding methods to overcome the more complex practical problems associated with the system of fractional q-difference equations. Moreover, we will investigate numerical methods for this kind of system.

Author Contributions

Conceptualization, C.Y. and J.W.; methodology, C.Y. and S.W.; validation, C.Y., J.W. and J.L.; formal analysis, J.W.; resources, S.W.; data curation, C.Y.; writing—original draft preparation, S.W.; writing—review and editing, C.Y., J.W. and J.L.; supervision, C.Y. and J.L.; funding acquisition, C.Y., J.W. and J.L. All authors have read and agreed to the published version of the manuscript.

Funding

The research project is supported by National Natural Science Foundation of China (12272011, 11772007), Beijing Natural Science Foundation (Z180005, 1172002), Natural Science Foundation of Hebei Province (A2015208114) and the Foundation of Hebei Education Department (QN2017063).

Data Availability Statement

Not applicable.

Acknowledgments

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

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jackson, F.H. On q-functions and a certain difference operator. Trans. Roy. Soc. Edin. 1909, 46, 253–281. [Google Scholar] [CrossRef]
  2. Agarwal, R.P. Certain fractional q-integrals and q-derivatives. Proc. Camb. Philol. Soc. 1969, 66, 365–370. [Google Scholar] [CrossRef]
  3. Al-Salam, W.A. Some fractional q-integrals and q-derivatives. Proc. Edinburgh. Math. Soc. 1966, 15, 135–140. [Google Scholar] [CrossRef]
  4. Koornwinder, T.H.; Swarttow, R.F. On q-analogues of the Fourier and Hankel transforms. Proc. Trans. Am. Math. Soc. 1992, 333, 445–461. [Google Scholar] [CrossRef]
  5. Atici, F.M.; Eloe, P.W. Fractional q-calculus on a time scale. J. Nonlinear Math. Phys. 2007, 14, 341–352. [Google Scholar] [CrossRef]
  6. Rajkovi<i>c</i>´, P.M.; Marinkovi<i>c</i>´, S.D.; Stankovi<i>c</i>´, M.S. Fractional integrals and derivatives in q-calculus. Appl. Anal. Discret. Math. 2007, 1, 311–323. [Google Scholar]
  7. Ahmad, B. Boundary-value problems for nonlinear third-order q-difference equations. Proc. Electron. J. Differ. Equ. 2011, 2011, 107384. [Google Scholar] [CrossRef]
  8. Ahmad, B.; Ntouyas, S.K. Boundary value problems for q-difference inclusions. Abstr. Appl. Anal. 2011, 2011, 292860. [Google Scholar] [CrossRef]
  9. Ahmad, B.; Alsaedi, A.; Ntouyas, S.K. A study of second-order q-difference equations with boundary conditions. Adv. Differ. Equ. 2012, 2012, 35. [Google Scholar] [CrossRef]
  10. Ferreira, R.A.C. Nontrivial solutions for fractional q-difference boundary value problems. Electron. J. Qual. Theory Differ. Equ. 2010, 70, 1–10. [Google Scholar] [CrossRef]
  11. 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]
  12. Ma, J.; Yang, J. Existence of solutions for multi-point boundary value problem of fractional q-difference equation. Electron. J. Qual. Theory 2011, 92, 1–10. [Google Scholar] [CrossRef]
  13. Yu, C.; Wang, J. Positive solutions of nonlocal boundary value problem for high-order nonlinear fractional q-difference equations. Abstr. Appl. Anal. 2013, 2013, 928147. [Google Scholar] [CrossRef]
  14. 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. [Google Scholar] [CrossRef]
  15. Yu, C.; Wang, J. Existence of solutions for nonlinear second-order q-difference equations with first-order q-derivatives. Adv. Differ. Equ. 2013, 2013, 124. [Google Scholar] [CrossRef]
  16. Yu, C.; Li, J.; Wang, J. Existence and uniqueness criteria for nonlinear quantum difference equations with p-Laplacian. AIMS Math. 2022, 7, 10439–10453. [Google Scholar] [CrossRef]
  17. Wang, Y. Positive solutions for fractional differential equation involving the Riemann-Stieltjes integral conditions with two parameters. J. Nonlinear Sci. Appl. 2016, 9, 5733–5740. [Google Scholar] [CrossRef]
  18. Ahmad, B.; Alsaedi, A.; Alruwaily, Y. On Riemann-Stieltjes integral boundary value problems of Caputo-Riemann-Liouville type fractional integro-differential. Filomat 2020, 34, 2723–2738. [Google Scholar] [CrossRef]
  19. Hu, Z.; Liu, W.; Liu, J. Existence of solutions for a coupled system of fractional p-Laplacian equations at resonance. Adv. Differ. Equ. 2013, 2013, 312. [Google Scholar] [CrossRef]
  20. Zhao, X.; Kang, S.; Gao, Y. Existence and uniqueness of solution to a coupled system of fractional difference equations with nonlocal conditions. J. Biomath. 2013, 28, 302–306. [Google Scholar]
  21. Hu, L.; Zhang, S. Existence and uniqueness of solutions for a higher-order coupled fractional differential equations at resonance. Adv. Differ. Equ. 2015, 2015, 202. [Google Scholar] [CrossRef]
  22. Lyu, P.; Vong, S. A linearized and second-order unconditionally convergent scheme for coupled time fractional Klein-Gordon-Schrodinger equation. Numer. Meth. Part Differ. Equ. 2018, 34, 2153–2179. [Google Scholar] [CrossRef]
  23. Nouara, A.; Amara, A.; Kaslik, E.; Etemad, S.; Rezapour, S.; Martinez, F.; Kaabar, M.K. A study on multiterm hybrid multi-order fractional boundary value problem coupled with its stability analysis of Ulam-Hyers type. Adv. Differ. Equ. 2021, 2021, 343. [Google Scholar] [CrossRef]
  24. Tudorache, A.; Luca, R. Positive solutions for a system of Riemann-Liouville fractional boundary value problems with p-Laplacian operators. Adv. Differ. Equ. 2020, 2020, 292. [Google Scholar] [CrossRef]
  25. Luca, R. On a system of Riemann-Liouville fractional differential equations with coupled nonlocal boundary conditions. Adv. Differ. Equ. 2021, 2021, 134. [Google Scholar] [CrossRef]
  26. Guo, D. Nonlinear Functional Analysis; Shandong Science and Technology Press: Jinan, China, 2001; pp. 1–559. [Google Scholar]
  27. Page, D.N. Information in black hole radiation. Phys. Rev. Lett. 1993, 71, 3743–3746. [Google Scholar] [CrossRef]
  28. Youm, D. Q-deformed conformal quantum mechanics. Phys. Rev. D 2000, 62, 276–284. [Google Scholar] [CrossRef]
  29. Lavagno, A.; Swamy, P.N. Q-deformed structures and nonextensive statistics: A comparative study. Physica A 2002, 305, 310–315. [Google Scholar] [CrossRef]
  30. Annaby, M.H.; Mansour, Z.S. Q-Fractional Calculus and Equations. In Lecture Notes in Mathematics 2056; Springer: Berlin, Germany, 2012; pp. 1–318. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yu, C.; Wang, S.; Wang, J.; Li, J. Solvability Criterion for Fractional q-Integro-Difference System with Riemann-Stieltjes Integrals Conditions. Fractal Fract. 2022, 6, 554. https://doi.org/10.3390/fractalfract6100554

AMA Style

Yu C, Wang S, Wang J, Li J. Solvability Criterion for Fractional q-Integro-Difference System with Riemann-Stieltjes Integrals Conditions. Fractal and Fractional. 2022; 6(10):554. https://doi.org/10.3390/fractalfract6100554

Chicago/Turabian Style

Yu, Changlong, Si Wang, Jufang Wang, and Jing Li. 2022. "Solvability Criterion for Fractional q-Integro-Difference System with Riemann-Stieltjes Integrals Conditions" Fractal and Fractional 6, no. 10: 554. https://doi.org/10.3390/fractalfract6100554

APA Style

Yu, C., Wang, S., Wang, J., & Li, J. (2022). Solvability Criterion for Fractional q-Integro-Difference System with Riemann-Stieltjes Integrals Conditions. Fractal and Fractional, 6(10), 554. https://doi.org/10.3390/fractalfract6100554

Article Metrics

Back to TopTop