Next Article in Journal
Finite-Time Boundedness of Linear Uncertain Switched Positive Time-Varying Delay Systems with Finite-Time Unbounded Subsystems and Exogenous Disturbance
Next Article in Special Issue
Hilfer Fractional Quantum Derivative and Boundary Value Problems
Previous Article in Journal
Inverse Result of Approximation for the Max-Product Neural Network Operators of the Kantorovich Type and Their Saturation Order
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Some Results on the q-Calculus and Fractional q-Differential Equations

Department of Mathematics, Northeastern University, Shenyang 110167, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2022, 10(1), 64; https://doi.org/10.3390/math10010064
Submission received: 1 December 2021 / Revised: 22 December 2021 / Accepted: 23 December 2021 / Published: 25 December 2021
(This article belongs to the Special Issue Q-differential/Difference Equations and Related Applications)

Abstract

:
In this paper, we first discuss some important properties of fractional q-calculus. Then, based on these properties and the q-Laplace transform, we translate a class of fractional q-differential equations into the equivalent q-differential equations with integer order. Thus, we propose a method for solving some linear fractional q-differential equations by means of solving the corresponding integer order equations. Several examples are provided to illustrate our solution method.

1. Introduction

Recently, mathematicians and engineers have paid much attention to the q-calculus and fractional q-differential equations. The q-calculus (also called the quantum calculus) can be dated back to 1908, Jackson’s work [1]. Based on the q-calculus, the q-differential equations were established which can describe some special physical processes appearing in quantum dynamics, discrete dynamical systems and discrete stochastic processes, and so forth. It should be pointed out that the q-differential equations are usually defined on a time scale set T q , where q is the scale index. With the development of the q-calculus theory, some related concepts have also been introduced and studied such as the q-Laplace transform, q-Gamma and q-Beta functions, q-Mittag–Leffler functions, q-Taylor expansion, q-integral transforms theory and so forth. Please refer to articles [2,3,4,5,6,7,8,9,10] for more details on q-calculus and fractional q-differential equations. Up to now, compared with the classical fractional calculus, the study of the fractional q-calculus is still immature.
At present, there have been some studies on the existence and uniqueness of solutions for the fractional q-differential equations. In [11], Abdeljawad et al. proved the uniqueness of an initial value problem involving a nonlinear delay Caputo fractional q-difference system by using a new generalized version of discrete fractional q-Gronwall inequality. In [12], the authors proved that the Caputo q-fractional boundary value problem with the p-Laplacian operator has a unique solution by using the Banach’s contraction mapping principle. In [13], Ren et al. considered the uniqueness of nontrivial solutions by virtue of the contraction mapping principle and also obtained the existence of multiple positive solutions under some appropriate conditions via standard fixed point theorems. In [14], Zhang et al. gave the existence and uniqueness of solution of the Caputo fractional q-differential equations by using the Ascoli–Arzela Theorem and a q-analogue Gronwall inequality. Furthermore, Zhang et al. in [15] discussed the existence of a unique solution in the q-integral space.
Although many research results have been given about the existence and uniqueness of the fractional q-differential equations of various types, there are very few studies about how to solve these problems analytically. The applications of the q-Mittag–Leffler function for presenting the solution of initial value problems of linear Caputo fractional q-differential equations were studied by Abdeljawad in [16,17]. Following the previous work, Abdeljawad et al. [11] provided a particular solution formula for the nonlinear delay Caputo fractional q-difference system by means of the q-Mittag–Leffler function. In [18], the solutions of some Caputo fractional q-difference equations are expressed by means of a new generalized type q-Mittag–Leffler function. In [15], Zhang et al. discussed the unique existence of solutions for the Caputo type nonlinear fractional q-differential equations and obtained the solutions expressed by means of the q-Mittag–Leffler function, which is defined as follows
E α , β q ( λ , z z 0 ) = k = 0 λ k ( z z 0 ) q ( k α ) Γ q ( β + k α ) .
We see that the q-Mittag–Leffler function is an infinite series. Although, in [15,18], the solutions were obtained by means of the q-Mittag–Leffler function, it is very difficult to express the solution in a finite analytic form, because the solution formula concerns the fractional q-integral of q-Mittag–Leffler function.
In this paper, we aim to find a method to obtain an analytic solution in a finite form for the initial value problem of the Caputo type fractional q-differential equation
c D q α u ( t ) = λ u ( t ) + f ( t ) , α > 0 , D q k u ( 0 ) = d k , k = 0 , 1 , · · · , m , m = [ α ] .
We first study some properties of the fractional q-calculus. Then we translate the above fractional q-differential equation into the equivalent q-differential equation with integer order by using the properties derived in this paper and the q-Laplace transform. Thus, we propose a method for solving the fractional q-differential equations by solving the corresponding integer order equations. At last, several examples are provided to illustrate our solution method.
This paper is organized as follows. Section 2 introduces some basic definitions and relevant results on q-calculus. In Section 3, we study some important properties about q-calculus. Section 4 is devoted to a new solution method to find out the analytic solutions of a class of fractional q-differential equations. In Section 5, we provide some examples to illustrate our solution method.

2. Preliminaries

Let N = { 1 , 2 , . . . } be the positive integer set and 0 < q < 1 . Introduce the q-shifted operations
( t s ) q ( 0 ) = 1 , ( t s ) q ( k ) = i = 0 k 1 ( t q i s ) , k N .
Let C be the set of complex numbers. If α C and α N , the q-shifted operation is defined by
( t s ) q ( α ) = t α i = 0 t q i s t q α + i s , 0 s t .
Introduce the notations
[ α ] q = 1 q α 1 q , [ n ] q ! = [ n ] q [ n 1 ] q · · · [ 1 ] q .
For α C \ { n , n N { 0 } } , define the q-analogue Gamma function Γ q ( α ) as
Γ q ( α ) = ( 1 q ) q ( α 1 ) ( 1 q ) 1 α , 0 < q < 1 .
Then, from the definition it is easy to see that
Γ q ( 1 ) = 1 , Γ q ( n + 1 ) = [ n ] q ! , Γ q ( α + 1 ) = [ α ] q Γ q ( α ) .
For α , β C and R e ( α ) > 0 , R e ( β ) > 0 , define the q-Beta function
B q ( α , β ) = 0 1 s α 1 ( 1 q s ) q ( β 1 ) d q s .
The q-Gamma and q-Beta functions have the following relationship [2]
B q ( α , β ) = Γ q ( α ) Γ q ( β ) Γ q ( α + β ) .
For a given q R , a subset A R is called q-geometric if q x A whenever x A . That is, x A , set A includes all geometric sequences { x q n } n = 0 . A typical q-geometric set is T q = { q n : n Z } { 0 } , where 0 < q < 1 , Z is the set of integers.
Definition 1
([1]). Let f ( t ) be a real-valued function on the set T q and 0 < q < 1 . Define the q-derivative of f ( t ) by
D q f ( t ) = d q f ( t ) d q t = f ( t ) f ( q t ) ( 1 q ) t , t T q \ { 0 } ,
D q f ( 0 ) = d q f ( t ) d q t | t = 0 = lim n f ( t q n ) f ( 0 ) t q n , t 0 .
From the definition, we see that the q-derivative is different from the classical derivative. It is a kind of discrete analogue of the classical derivative. On the basis of Definition 1, the higher order q-derivative D q n f ( t ) is defined by D q n f ( t ) = D q ( D q n 1 f ( t ) ) , n 2 .
Furthermore, for two real-valued functions f ( t ) and g ( t ) , we have the following operation rules by a straightforward computation
D q ( a f ( t ) ± b g ( t ) ) = a D q f ( t ) ± b D q g ( t ) , a , b R , D q ( f ( t ) g ( t ) ) = g ( t ) D q f ( t ) + f ( q t ) D q g ( t ) , D q ( f ( t ) g ( t ) ) = g ( t ) D q f ( t ) f ( t ) D q g ( t ) g ( t ) g ( q t ) , g ( t ) 0 , g ( q t ) 0 .
Definition 2
([19]). Let f ( t ) be a real-valued function defined on the set T q and 0 < q < 1 . The q-integral of f ( t ) is defined by
0 t f ( s ) d q s = ( 1 q ) n = 0 t q n f ( t q n ) , t T q ,
and for b > a
a b f ( s ) d q s = 0 b f ( s ) d q s 0 a f ( s ) d q s , a , b T q .
From Definition 2, it holds that
| 0 b f ( s ) d q s | 0 b | f ( s ) | d q s , b > 0 ,
a b f ( s ) d q s = a c f ( s ) d q s + c b f ( s ) d q s , a < c < b .
The operation of q-integration by parts is given in the lemma below.
Lemma 1
([20]). Suppose that f ( t ) and g ( t ) are real-valued functions defined on the set T q , 0 < q < 1 , a , b T q , 0 a < b . Then
a b f ( q t ) D q g ( t ) d q t = ( f g ) ( b ) ( f g ) ( a ) a b g ( t ) D q f ( t ) d q t .
Now, we introduce the concept of fractional q-calculus. Let A be a q-geometric set.
Definition 3
([21]). Suppose that t A and α 1 , 2 , · · · . The α-order fractional q-integral of the Riemann–Liouville type is defined by I q 0 f ( t ) = f ( t ) and
I q α f ( t ) = 1 Γ q ( α ) 0 t ( t q s ) q ( α 1 ) f ( s ) d q s .
Definition 4
([22]). Let n = [ α ] and f ( t ) be a real-valued function defined on ( 0 , ) . The α-order Caputo type fractional q-derivative of function f ( t ) is defined by
c D q α f ( t ) = I q α f ( t ) , α 0 , I q n α D q n f ( t ) , α > 0 ,
where [ α ] represents the smallest integer which is greater than or equal to α.
For α > 0 , we have
c D q α f ( t ) = 1 Γ q ( n α ) 0 t ( t q s ) ( n α 1 ) D q n f ( s ) d q s .
Definition 5
([16]). Suppose that α R and n = [ α ] , f ( t ) is a real-valued function on ( 0 , ) . The α-order Riemann–Liouville type fractional q-derivative of f ( t ) is defined as follows
D q α f ( t ) = I q α f ( t ) , α 0 , D q n I q n α f ( t ) , α > 0 .
For α > 0 , we have
D q α f ( t ) = 1 Γ q ( n α ) D q n 0 t ( t q s ) q ( n α 1 ) f ( s ) d q s .
Under certain conditions, the Riemann–Liouville type fractional q-derivative and Caputo type fractional q-derivative have the following relationship [2,16]
c D q α f ( t ) = D q α ( f ( t ) k = 0 n 1 D q k f ( 0 ) Γ q ( k + 1 ) t k ) , t > 0 , α > 0 , n = [ α ] .

3. Some Properties on q -Calculus

Introduce the following q-integrable and q-continuous differentiable function spaces
q L p ( 0 , b ) = { f ( t ) : f ( t ) q L p ( 0 , b ) = sup 0 < t b , ( t , b ) T q ( 0 t | f ( s ) | p d q s ) 1 p < } ,
C q ( m ) [ 0 , b ] = { f ( t ) : D q k f ( t ) C [ 0 , b ] , k = 0 , 1 , · · · , m . } .
Similarly, we can define C q ( m ) ( 0 , b ] . By the definition of q-derivative, we see that if f ( t ) C [ 0 , b ] , then f ( t ) C q ( k ) ( 0 , b ] , k 0 . Furthermore, if f ( t ) C [ 0 , b ] and f ( 0 ) exists, then f ( t ) C q ( 1 ) [ 0 , b ] .
Define the integer order q-integral operator
I q 0 f ( t ) = f ( t ) , I q f ( t ) = 0 t f ( s ) d q s , I q n f ( t ) = I q n 1 I q f ( t ) , n = 1 , 2 , · · · .
Lemma 2.
Let f ( t ) C [ 0 , b ] , 0 < q < 1 , then we have
D q I q f ( t ) = f ( t ) , 0 < t b ,
I q D q f ( t ) = f ( t ) f ( 0 ) , 0 < t b .
Proof. 
From (5) and (6), we have
D q I q f ( t ) = ( 0 t f ( s ) d q s 0 q t f ( s ) d q s ) / ( 1 q ) t = [ ( 1 q ) k = 0 t q k f ( t q k ) ( 1 q ) k = 0 t q k + 1 f ( t q k + 1 ) ] / ( 1 q ) t = k = 0 q k f ( t q k ) k = 0 q k + 1 f ( t q k + 1 ) = f ( t ) .
Then (17) holds. For (18), from (6) and (5), we obtain
I q D q f ( t ) = 0 t f ( s ) f ( q s ) ( 1 q ) s d q s = ( 1 q ) k = 0 t q k f ( t q k ) f ( t q k + 1 ) ( 1 q ) t q k = k = 0 ( f ( t q k ) f ( t q k + 1 ) ) = f ( t ) lim k f ( t q k ) = f ( t ) f ( 0 + ) .
This gives (18). □
Lemma 3.
Let f ( t ) C [ 0 , b ] , 0 < q < 1 , α 1 , 2 , · · · . Then
lim α 0 + I q α f ( t ) = f ( t ) .
Proof. 
Since D q , s ( t s ) q ( α ) = [ α ] q ( t q s ) q ( α 1 ) (see (2.6) in [20]), we have from (10) that
I q α f ( t ) = 1 Γ q ( α ) 0 t ( t q s ) q ( α 1 ) f ( s ) d q s = 1 Γ q ( α ) [ α ] q 0 t D q , s ( t s ) q ( α ) f ( s ) d q s = 1 Γ q ( α + 1 ) { ( t s ) q ( α ) f ( s ) | 0 t 0 t ( t q s ) q ( α ) D q f ( s ) d q s } = 1 Γ q ( α + 1 ) { t α f ( 0 ) + 0 t ( t q s ) q ( α ) D q f ( s ) d q s } .
So by Lemma 2,
lim α 0 + I q α f ( t ) = f ( 0 ) + 0 t D q f ( s ) d q s = f ( t ) .
Corollary 1.
Let f ( t ) C [ 0 , b ] . Then we have
lim α n c D q α f ( t ) = D q n f ( t ) , 0 < q < 1 , n 1 < α < n , n 1 , t > 0 .
Proof. 
From (12) and Lemma 3, we obtain
lim α n c D q α f ( t ) = lim α n I q n α D q n f ( t ) = D q n f ( t ) .
Lemma 4
([2]). If f ( t ) q L 1 [ 0 , b ] , then the semigroup property holds
I q α I q β f ( t ) = I q β I q α f ( t ) = I q α + β f ( t ) , α 0 , β 0 .
Lemma 5.
Let f ( t ) q L 1 [ 0 , b ] , then
lim α n I q α f ( t ) = I q n f ( t ) , n 1 < α < n .
Proof. 
When n = 1 ,
lim α 1 I q α f ( t ) = lim α 1 1 Γ q ( α ) 0 t ( t q s ) q ( α 1 ) f ( s ) d q s = 0 t f ( s ) d q s = I q f ( t ) .
When n > 1 , let α = α + n 1 and 0 < α < 1 , by Lemma 4,
lim α n I q α f ( t ) = lim α 1 I q n 1 + α f ( t ) = lim α 1 I q α I q n 1 f ( t ) = I q n f ( t ) .
Lemma 6
([2]). Let 0 < α < 1 and f ( t ) C [ 0 , b ] such that D q f ( t ) C ( 0 , b ] . Then
c D q α f ( t ) = D q α f ( t ) t α Γ q ( 1 α ) f ( 0 ) , t > 0 .
Lemma 7.
The following properties hold
(i) 
If | f ( t ) | M , then I q α f ( 0 ) = 0 , α > 0 ,
(ii) 
If | D q f ( t ) | M , then c D q α f ( 0 ) = 0 , 0 < α < 1 .
Proof. 
For property (i), using (8) we can obtain
| I q α f ( t ) | 1 Γ q ( α ) 0 t ( t q s ) q ( α 1 ) | f ( s ) | d q s M Γ q ( α ) 0 t ( t q s ) q ( α 1 ) d q s = M Γ q ( α ) t α 0 1 ( 1 q τ ) q ( α 1 ) d q τ = M Γ q ( α ) t α B q ( α , 1 ) = t α M Γ q ( 1 + α ) 0 , t 0 .
Next, from (i) and c D q α f ( t ) = I q 1 α D q f ( t ) , we see that the property (ii) holds. □
Lemma 8.
For 0 < α < 1 and f ( t ) C [ 0 , b ] , it holds
D q α I q α f ( t ) = f ( t ) , c D q α I q α f ( t ) = f ( t ) .
Proof. 
By Definition 5 and Lemma 2, we see that
D q α I q α f ( t ) = D q I q 1 α I q α f ( t ) = D q I q f ( t ) = f ( t ) .
Then, the first equality holds. From Lemma 6 and noting I q α f ( 0 ) = 0 ,
c D q α I q α f ( t ) = D q α I q α f ( t ) = f ( t ) .
The proof is completed. □
Lemma 9.
For 0 < α < 1 , n 1 < β < n , n 1 < α + β n and n 1 , we have
c D q α c D q β f ( t ) = c D q α + β f ( t ) .
Proof. 
Applying Lemma 4 and Lemma 8, it yields that
c D q α c D q β f ( t ) = c D q α I q n β D q n f ( t ) = c D q α I q α I q n α β D q n f ( t ) = c D q α I q α c D q α + β f ( t ) = c D q α + β f ( t ) .
Remark 1.
In Lemma 9, condition n 1 < α + β n is necessary. For n = 1 , a counterexample is as follows
c D q 0.7 c D q 0.5 t = 1 Γ q ( 4 5 ) t 0.2 , c D q 1.2 t = 0 .

4. The Fractional q -Differential Equations

In this section, we establish the solution method for a class of linear fractional q-differential equations.
Let us consider the fractional q-differential equations in the following form
c D q α u ( t ) = λ u ( t ) + f ( t ) , 0 < t b , D q k u ( 0 ) = d k , k = 0 , 1 , · · · , m 1 ,
where { d k } are constants, m = [ α ] , λ is a constant and α = m 1 + r p , r , p N , 1 r < p .
Introduce the q-Mittag–Leffler function
E α , β q ( λ , z z 0 ) = k = 0 λ k ( z z 0 ) q ( k α ) Γ q ( β + k α ) ,
where { λ , z , z 0 , α , β } C and R e ( α ) > 0 , R e ( β ) 0 .
Theorem 1.
Let u ( t ) be the solution of problem (26) and assume that | u ( t ) | M , | λ b α | 1 . Then u ( t ) can be represented as follows
u ( t ) = k = 0 m 1 E α , k + 1 q ( λ , t ) d k t k + 0 t ( t q s ) q ( α 1 ) E α , α q ( λ , t q α s ) f ( s ) d q s .
Proof. 
According to ([15], Theorem 1), the solution of problem (26) satisfies the equivalent integral equation
u ( t ) = k = 0 m 1 d k t k Γ q ( k + 1 ) + 1 Γ q ( α ) 0 t ( t q s ) q ( α 1 ) [ λ u ( s ) + f ( s ) ] d q s .
Let d ( t ) = k = 0 m 1 d k t k Γ q ( k + 1 ) . Then by a circulative iteration and using Lemma 4, we obtain from (28) that
u ( t ) = d ( t ) + λ I q α u ( t ) + I q α f ( t ) = d ( t ) + λ I q α [ d ( t ) + λ I q α u ( t ) + I q α f ( t ) ] + I q α f ( t ) = d ( t ) + λ I q α d ( t ) + λ 2 I q 2 α u ( t ) + λ I q 2 α f ( t ) + I q α f ( t ) = d ( t ) + λ I q α d ( t ) + λ 2 I q 2 α [ d ( t ) + λ I q α u ( t ) + I q α f ( t ) ] + λ I q 2 α f ( t ) + I q α f ( t ) = d ( t ) + λ I q α d ( t ) + λ 2 I q 2 α d ( t ) + λ 3 I q 3 α u ( t ) + λ 2 I q 3 α f ( t ) + λ I q 2 α f ( t ) + I q α f ( t ) · · · · · · · · · · · · · · · · · · · · · · · · = n = 0 λ n I q n α d ( t ) + n = 0 λ n I q ( n + 1 ) α f ( t ) + lim n λ n I q n α u ( t ) .
Since
λ n I q n α d ( t ) = λ n I q n α [ k = 0 m 1 t k Γ q ( k + 1 ) d k ] = k = 0 m 1 λ n t n α + k Γ q ( n α + k + 1 ) d k ,
then
n = 0 λ n I q n α d ( t ) = k = 0 m 1 [ n = 0 λ n t n α + k Γ q ( n α + k + 1 ) ] d k = k = 0 m 1 [ n = 0 λ n t n α Γ q ( n α + k + 1 ) ] t k d k = k = 0 m 1 E α , k + 1 q ( λ , t ) t k d k .
Next, by the identity [21] ( t s ) q ( β + γ ) = ( t s ) q ( β ) ( t q β s ) q ( γ ) , we have
n = 0 λ n I q ( n + 1 ) α f ( t ) = n = 0 λ n Γ q ( n α + α ) 0 t ( t q s ) q ( n α + α 1 ) f ( s ) d q s = 0 t ( t q s ) q ( α 1 ) n = 0 λ n ( t q α s ) q ( n α ) Γ q ( n α + α ) f ( s ) d q s = 0 t ( t q s ) q ( α 1 ) E α , α q ( λ , t q α s ) f ( s ) d q s .
Furthermore, by the argument of Lemma 7
| λ n I q n α u ( t ) | ( λ b α ) n Γ q ( n α + 1 ) M M Γ q ( n α + 1 ) 0 , n .
Then, combining (29)–(32), we complete the proof. □
Although formula (27) gives a solution representation of problem (26), we see that it is very difficult to obtain a finite form solution by formula (27), which concerns the complex fractional q-integral. Below, we will transform problem (26) into an equivalent integer order q-differential equation and then by using this integer order equation to find the solution of problem (26).
We first consider the case of m = 1 and r = 1 in problem (26). Corresponding to problem (26), introduce the auxiliary problem
D q z ( t ) = λ p z ( t ) + w ( t ) , 0 < t b , z ( 0 ) = u ( 0 ) ,
where
w ( t ) = λ p k = 1 p θ k t k p + k = 0 p 1 λ p ( k + 1 ) c D q k p f ( t ) θ p ,
θ k = λ k u ( 0 ) + λ k 1 f ( 0 ) Γ q ( 1 + k p ) , k = 1 , 2 , · · · , p .
Furthermore, we introduce the q-Laplace transform q L s ( · ) which is a linear transform and has the following properties [2]
q L s ( t γ Γ q ( γ + 1 ) ) = ( 1 q ) γ s γ + 1 , γ > 1 , R e ( s ) > 0 ,
q L s ( c D q α f ( t ) ) = s α ( 1 q ) α ( q L s f ) ( s ) f ( 0 ) s α 1 ( 1 q ) α , 0 < α < 1 ,
q L s ( D q f ( t ) ) = s 1 q ( q L s f ) ( s ) f ( 0 ) 1 q .
Theorem 2.
Let functions u ( t ) and z ( t ) satisfy the relationship
u ( t ) = z ( t ) + k = 1 p θ k t k p .
Then, for m = 1 and r = 1 , u ( t ) is a solution of fractional q-differential problem (26) if and only if z ( t ) is a solution of integer order problem (33).
Proof. 
We first assume that u ( t ) is the solution of problem (26) and (39) holds. We need to prove that z ( t ) is the solution of problem (33). Denote the q-Laplace transform of function u ( t ) by u ¯ ( s ) . Applying q-Laplace transform to Equation (26) with m = 1 , r = 1 and using (37), it yields
( S q 1 p λ ) u ¯ ( s ) = 1 1 q S q 1 p 1 u ( 0 ) + f ¯ ( s ) ,
where S q = s 1 q . Multiplying (40) by k = 0 p 1 S q p 1 k p λ k and using the identity
a n b n = ( a b ) k = 0 n 1 a n 1 k b k ,
we obtain
( S q λ p ) u ¯ ( s ) = 1 1 q u ( 0 ) k = 0 p 1 S q 1 p 1 S q p 1 k p λ k + k = 0 p 1 S q p 1 k p λ k f ¯ ( s ) .
Again applying q-Laplace transform to Equation (39) and using (36), we have
u ¯ ( s ) = z ¯ ( s ) + k = 1 p θ k Γ q ( k p + 1 ) ( 1 q ) k p s ( k p + 1 ) = z ¯ ( s ) + k = 1 p θ k Γ q ( k p + 1 ) 1 1 q S q ( k p + 1 ) .
Combining Equation (41) and (42) we obtain
( S q λ p ) z ¯ ( s ) = 1 1 q u ( 0 ) k = 0 p 1 S q k p λ k + k = 0 p 1 S q p 1 k p λ k f ¯ ( s ) k = 1 p θ k Γ q ( k p + 1 ) S q ( k p + 1 ) 1 1 q ( S q λ p ) = λ p 1 q k = 1 p θ k Γ q ( k p + 1 ) S q ( k p + 1 ) 1 1 q k = 1 p θ k Γ q ( k p + 1 ) S q k p + + 1 1 q u ( 0 ) k = 0 p 1 S q k p λ k + k = 0 p 1 S q p 1 k p λ k f ¯ ( s ) = λ p 1 q k = 1 p θ k Γ q ( k p + 1 ) S q ( k p + 1 ) 1 1 q θ p Γ q ( 2 ) S q 1 + u ( 0 ) 1 q + + k = 1 p 1 1 ( 1 q ) S q k / p [ λ k u ( 0 ) θ k Γ q ( k p + 1 ) ] + k = 0 p 1 S q p 1 k p λ k f ¯ ( s ) .
Then, using (35), we have
( S q λ p ) z ¯ ( s ) = λ p 1 q k = 1 p θ k Γ q ( k p + 1 ) S q ( k p + 1 ) 1 1 q θ p S q 1 + + u ( 0 ) 1 q k = 1 p 1 λ k 1 f ( 0 ) 1 q S q k p + k = 0 p 1 λ k S q p 1 k p f ¯ ( s ) .
Since u ( 0 ) = z ( 0 ) (see (39)), S q = s 1 q , we have from (43) that
s 1 q z ¯ ( s ) 1 1 q z ( 0 ) λ p z ¯ ( s ) = λ p k = 1 p θ k Γ q ( k p + 1 ) s k p 1 ( 1 q ) k p + k = 0 p 2 [ λ k s 1 k + 1 p 1 ( 1 q ) p k 1 p f ¯ ( s ) λ k s k + 1 p f ( 0 ) 1 ( 1 q ) p k 1 p ] θ p s 1 + λ p 1 f ¯ ( s ) .
Applying the inverse q-Laplace transform to Equation (44) and using (36)–(38), we arrive at
D q z ( t ) λ p z ( t ) = λ p k = 1 p θ k t k p θ p + k = 0 p 2 c D q p k 1 p f ( t ) λ k + λ p 1 f ( t ) = λ p k = 1 p θ k t k p θ p + k = 0 p 1 c D q p k 1 p f ( t ) λ k .
So z ( t ) satisfies Equation (33). Conversely, assume that z ( t ) is the solution of problem (33) and relation (39) holds. Then, using again the above equalities and a backward deduction, we can prove that u ( t ) is the solution of problem (26). □
Next, we consider the following fractional q-differential system of equations
c D q 1 p U ( t ) = A U ( t ) + F ( t ) , U ( 0 ) = U 0 ,
where U ( t ) = ( u 1 ( t ) , · · · , u n ( t ) ) T , F ( t ) = ( f 1 ( t ) , · · · , f n ( t ) ) T , U 0 = ( u 1 ( 0 ) , · · · , u n ( 0 ) ) T , A R n × n is a matrix. Similar to the argument of Theorem 2, we have the following result.
Theorem 3.
The system of equations (45) has the following solution
U ( t ) = Z ( t ) + k = 1 p Θ k t k p ,
where Z ( t ) is the solution of problem
D q Z ( t ) = A p Z ( t ) + W ( t ) ,
and
W ( t ) = A p k = 1 p Θ k t k p + k = 0 p 1 A p ( k + 1 ) c D q k p F ( t ) Θ p ,
where
Θ k = A k U ( 0 ) + A k 1 F ( 0 ) Γ q ( k + p p ) , k = 1 , 2 , · · · , p .
In order to extend the conclusions of Theorem 2 and Theorem 3 to more general equations, we give the following results, which can be proved by using Lemma 9 and Lemma 7.
Theorem 4.
For m = 1 , 1 < r < p , the initial value problem (26) is equivalent to the following fractional q-differential system
c D q 1 p u 1 ( t ) = u 2 ( t ) , c D q 1 p u 2 ( t ) = u 3 ( t ) , · · · · · · c D q 1 p u r ( t ) = λ u 1 ( t ) + f ( t ) , u 1 ( 0 ) = u ( 0 ) , u k ( 0 ) = 0 , k = 2 , · · · , r .
Theorem 5.
For α = m 1 + r p , 1 r < p , m > 1 , the initial value problem (26) is equivalent to the following fractional q-differential system
c D q 1 p u 1 ( t ) = u 2 ( t ) , c D q 1 p u 2 ( t ) = u 3 ( t ) , · · · · · · c D q 1 p u ( m 1 ) p + r ( t ) = λ u 1 ( t ) + f ( t ) ,
with the initial value conditions
u ( k 1 ) p + 1 ( 0 ) = d k 1 , k = 1 , 2 , · · · , m , u l ( 0 ) = 0 , l ( k 1 ) p + 1 .

5. Examples

In this section, we give some examples to illustrate the validity of the proposed method for solving the fractional q-differential equations.
Example 1.
Consider the fractional q-differential equation
c D q 1 2 u ( t ) = u ( t ) + f ( t ) , u ( 0 ) = 1 ,
where f ( t ) = 1 t 1 Γ q ( 3 2 ) t 1 2 . According to Theorem 2, we only need to solve the integer order equation
c D q z ( t ) = z ( t ) + w ( t ) , z ( 0 ) = 1 ,
where
w ( t ) = k = 1 2 θ k t k 2 + k = 0 1 c D q k 2 f ( t ) θ 2 = t 2 .
It is easy to see that the solution z ( t ) = 1 t . So the solution of problem (46) is
u ( t ) = z ( t ) + k = 1 2 θ k t k p = 1 + t + 2 Γ q ( 3 2 ) t 1 2 .
On the other hand, by formula (27), the solution of problem (46) is
u ( t ) = E 1 2 , 1 q ( 1 , t ) + 0 t ( t q s ) ( 1 2 ) E 1 2 , 1 2 q ( 1 , t q 1 2 s ) ( 1 s 1 Γ q ( 3 2 ) s 1 2 ) d q s .
Compare these two solutions given by (47) and (48), respectively. Solution (47) is a finite analytic function, but solution (48) concerns the fractional q-integral of an infinite function series.
Example 2.
Consider a higher order problem
c D q 3 2 u ( t ) = D q u ( t ) + f ( t ) , u ( 0 ) = 1 , D q u ( 0 ) = 1 ,
where f ( t ) = 1 t 1 Γ q ( 3 2 ) t 1 2 . Let u 1 ( t ) = D q u ( t ) . Since c D q 3 2 u ( t ) = c D q 1 2 ( D q u ( t ) ) = c D q 1 2 u 1 ( t ) , then from (49), we see that u 1 ( t ) satisfies
c D q 1 2 u 1 ( t ) = u 1 ( t ) + f ( t ) , u 1 ( 0 ) = 1 .
According to the result of Example 1, we obtain
u 1 ( t ) = 1 + t + 2 Γ q ( 3 2 ) t 1 2 .
Then, by Lemma 2 we have
u ( t ) u ( 0 ) = I q D q u ( t ) = I q u 1 ( t ) .
That is,
u ( t ) = u ( 0 ) + 0 t 1 + s + 2 Γ q ( 3 2 ) s 1 2 d q s = 1 + t + 1 Γ q ( 3 ) t 2 + 2 Γ q ( 5 2 ) t 3 2 .
Note that, for problem (49), the solution formula (27) is not applicable.
Example 3.
Consider the following problem
c D q 2 3 u ( t ) = 2 u ( t ) + f ( t ) , u ( 0 ) = 0 .
where f ( t ) = 1 2 t 2 3 Γ q ( 5 3 ) . Let u 1 ( t ) = u ( t ) , u 2 ( t ) = D q 1 3 u 1 ( t ) . Then, problem (50) is equivalent to the following system of equations
c D q 1 3 u 1 ( t ) = u 2 ( t ) , u 1 ( 0 ) = 0 , c D q 1 3 u 2 ( t ) = 2 u 1 ( t ) + f ( t ) , u 2 ( 0 ) = 0 ,
or
c D q 1 3 U ( t ) = A U ( t ) + F ( t ) , U ( 0 ) = 0 ,
where
U ( t ) = ( u 1 , u 2 ) T , F ( t ) = ( 0 , f ) T , A = 0 1 2 0 .
According to Theorem 3, system (51) is equivalent to the following problem
D q Z ( t ) = A 3 Z ( t ) + W ( t ) , Z ( 0 ) = 0 ,
where W ( t ) = ( 4 t , 2 ) T is given by Theorem 3. It is easy to see that the solution of problem (52) is Z ( t ) = ( 0 , 2 t ) T . So the solution of problem (51) is
U ( t ) = Z ( t ) + k = 1 3 Θ k t k 3 = ( t 2 3 Γ q ( 5 3 ) , t 1 3 Γ q ( 4 3 ) ) T .
Thus, we obtain the solution of problem (50)
u ( t ) = u 1 ( t ) = t 2 3 Γ q ( 5 3 ) .
On the other hand, by formula (27), the solution of problem (50) is
u ( t ) = 0 t ( t q s ) ( 1 3 ) E 2 3 , 2 3 q ( 2 , t q 2 3 s ) ( 1 2 s 2 3 Γ q ( 5 3 ) ) d q s .
Obviously, the expression of solution (53) is more concise than that of solution (54).

6. Conclusions

We consider how to solve a class of linear fractional q-differential equations. By transforming this class of equations into the equivalent integer order q-differential equations, we establish a solution method that can find out the analytical solutions under certain conditions. Our work further develops the solution theory of fractional q-differential equations.

Author Contributions

Conceptualization, Y.S. and T.Z.; Formal analysis, Y.S. and T.Z.; Methodology, Y.S. and T.Z.; Validation, T.Z.; Writing–original draft, Y.S.; Writing–review & editing, Y.S. and T.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the editor and the referees for their positive comments and useful suggestions which have improved this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jackson, F.H. On q-functions and a certain difference operator. Trans. R. Soc. Edinb. 1908, 46, 64–72. [Google Scholar] [CrossRef]
  2. Annaby, M.H.; Mansour, Z.S. q-Fractional Calculus and Equations; Springer: New York, NY, USA, 2012. [Google Scholar]
  3. Agarwal, R.P. Certain fractional q-integrals and q-derivatives. Proc. Camb. Phil. Soc. 1969, 66, 365–370. [Google Scholar] [CrossRef]
  4. Aral, A.; Gupta, V.; Agarwal, R.P. Applications of q-Calculus in Operator Theory; Springer: New York, NY, USA, 2013. [Google Scholar]
  5. Kac, V.; Cheung, P. Quantum Calculus; Springer: New York, NY, USA, 2002. [Google Scholar]
  6. Abdi, W.H. On q-Laplace transforms. Proc. Natl. Acad. Sci. USA India (Sect. A) 1960, 29, 89–408. [Google Scholar]
  7. Abdi, W.H. Application of q-Laplace transform to the solution of certain q-integral equations. Rend. Circ. Mat. Palermo. 1962, 11, 245–257. [Google Scholar] [CrossRef]
  8. Annaby, M.H.; Mansour, Z.S. q-Taylor and interpolation series for Jackson q-difference operators. J. Math. Anal. Appl. 2008, 334, 472–483. [Google Scholar] [CrossRef] [Green Version]
  9. Askey, R. The q-gamma and q-beta functions. Appl. Anal. 1979, 8, 125–141. [Google Scholar] [CrossRef]
  10. Andrews, G.E.; Askey, R.; Roy, R. Special Functions; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
  11. Abdeljawad, T.; Alzabut, J.; Baleanu, D. A generalized q-fractional Gronwall inequality and its applications to nonlinear delay q-fractional difference systems. J. Inequal. Appl. 2016, 2016, 240. [Google Scholar] [CrossRef] [Green Version]
  12. Aktuğlu, H.; Özarslan, M.A. On the solvability of Caputo q-fractional boundary value probelm involving p-Laplacian operator. Abstr. Appl. Anal. 2013, 2013, 658617. [Google Scholar] [CrossRef]
  13. Ren, J.; Zhai, C.B. Nonlacal q-fractional boundary value problem with Stieltjes integral conditions. Nonlinear Anal. Model. 2019, 24, 582–602. [Google Scholar]
  14. Zhang, T.; Guo, Q.X. The solution theory of the nonlinear q-fractional differential equations. Appl. Math. Lett. 2020, 104, 106282. [Google Scholar] [CrossRef]
  15. Zhang, T.; Wang, Y.Z. The unique existence of solution in the q-integrable space for the nonlinear q-fractional differential equations. Fractals 2021, 29, 2150050. [Google Scholar] [CrossRef]
  16. Abdeljawad, T.; Benli, D. Caputo q-fractional initial value problems and a q-analogue Mittag-Leffler function. Commun. Nonlinear Sci. 2011, 16, 4682–4688. [Google Scholar] [CrossRef] [Green Version]
  17. Abdeljawad, T.; Benli, B.; Baleanu, D. A generalized q-Mittag-Leffler function by q-Caputo fractional linear equations. In Abstract and Applied Analysis; Hindawi Publishing Corporation: New York, NY, USA, 2012; p. 546062. [Google Scholar]
  18. Abdeljawad, T.; Benli, B. A quantum generalized Mittag-Leffler function via Caputo q-fractional equations. arXiv 2011, arXiv:1102.1585. [Google Scholar]
  19. Jackson, F.H. On q-definite integral. Quart. J. Pure Appl. Math. 1910, 41, 193–203. [Google Scholar]
  20. Zhang, T.; Tong, C. A difference method for solving the nonlinear q-fractional differential equations on time scales. Fractals 2020, 28, 2050121. [Google Scholar] [CrossRef]
  21. Atici, F.M.; Eloe, P.W. Fractional q-calculus on a time scales. J. Nonlinear Math. Phys. 2007, 14, 341–352. [Google Scholar] [CrossRef] [Green Version]
  22. Rajkovic, P.M.; Marinkovic, S.D.; Stankovic, M.S. On q-analogues of Caputo derivative and Mittag-Leffler function. Fract. Calc. Appl. Anal. 2007, 10, 359–373. [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

Sheng, Y.; Zhang, T. Some Results on the q-Calculus and Fractional q-Differential Equations. Mathematics 2022, 10, 64. https://doi.org/10.3390/math10010064

AMA Style

Sheng Y, Zhang T. Some Results on the q-Calculus and Fractional q-Differential Equations. Mathematics. 2022; 10(1):64. https://doi.org/10.3390/math10010064

Chicago/Turabian Style

Sheng, Ying, and Tie Zhang. 2022. "Some Results on the q-Calculus and Fractional q-Differential Equations" Mathematics 10, no. 1: 64. https://doi.org/10.3390/math10010064

APA Style

Sheng, Y., & Zhang, T. (2022). Some Results on the q-Calculus and Fractional q-Differential Equations. Mathematics, 10(1), 64. https://doi.org/10.3390/math10010064

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop