Next Article in Journal
An Intelligent System for Allocating Times to the Main Activities of Managers
Next Article in Special Issue
On a Nonlinear Mixed Problem for a Parabolic Equation with a Nonlocal Condition
Previous Article in Journal
Macroscopic and Multi-Scale Models for Multi-Class Vehicular Dynamics with Uneven Space Occupancy: A Case Study
Previous Article in Special Issue
Distributed-Order Non-Local Optimal Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Existence of Coupled Fractional Jerk Equations with Multi-Point Boundary Conditions

1
School of Science, Shandong Jiaotong University, Jinan 250357, China
2
School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
3
Department of Mathematics, China University of Mining and Technology, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Axioms 2021, 10(2), 103; https://doi.org/10.3390/axioms10020103
Submission received: 23 April 2021 / Revised: 19 May 2021 / Accepted: 22 May 2021 / Published: 24 May 2021

Abstract

:
By coincidence degree theory due to Mawhin, some sufficient conditions for the existence of solution for a class of coupled jerk equations with multi-point conditions are established. The new existence results have not yet been reported before. Novel coupled fractional jerk equations with resonant boundary value conditions are discussed in detail for the first time. Our work is interesting and complements known results.

1. Introduction

In this paper, we study the following coupled jerk system of the fractional order
D 0 + α x 1 ( t ) = x 2 ( t ) , D 0 + β x 2 ( t ) = x 3 ( t ) , D 0 + γ x 3 ( t ) = f y 1 ( t ) , y 2 ( t ) , y 3 ( t ) , D 0 + λ y 1 ( t ) = y 2 ( t ) , D 0 + μ y 2 ( t ) = y 3 ( t ) , D 0 + ν y 3 ( t ) = g x 1 ( t ) , x 2 ( t ) , x 3 ( t ) ,
with boundary value conditions given as
x 2 ( 0 ) = x 3 ( 0 ) = 0 , x 1 ( 0 ) = i = 1 m 1 a i x 1 ( ξ i ) , y 2 ( 0 ) = y 3 ( 0 ) = 0 , y 1 ( 0 ) = i = 1 m 1 b i y 1 ( η i ) ,
where t ( 0 , 1 ) , 0 ξ 1 ξ 2 ξ m 1 1 , 0 η 1 η 2 η m 1 1 , i = 1 m 1 a i = i = 1 m 1 b i = 1 , i = 1 m 1 a i ξ i 0 , i = 1 m 1 b i η i 0 , 0 < ε < 1 , ε = { α , β , γ , λ , μ , ν } , D 0 + ε denotes the Caputo fractional derivative, and f , g : R 3 R are continuous functions.
In 1978, S. H. Schot [1] presented the definition of a jerk, i.e., the rate of change of acceleration. This involves a third derivative of x and has been found to have numerous applications in many areas of science, such as electrical circuits, laser physics, mechanics, acoustics, and dynamical processes. We refer the reader to [2,3,4,5,6,7,8,9,10,11].
As we all know, a three-dimensional dynamical system x ˙ ( t ) = y , y ˙ ( t ) = z , z ˙ ( t ) = f ( x , y , z ) , can be written in the form x ( t ) = f ( x , x ˙ , x ¨ ) . The third order autonomous differential equation is named as a jerk equation. In recent years, the discussion of jerk equations has attracted much attention because it arises in a variety of different scientific fields, such as the theory of chaos, secure communication, electrical engineering, and economic systems. For recent results, we refer the readers to [12,13,14,15,16,17,18,19,20] and the references therein.
In 2020, Marcelo and Silva [18] used the algebraic criterion to determine a special form for the polynomial jerk function j in order to guarantee the nonchaotic behavior of the following jerk equation:
x ( t ) = j ( x , x ˙ , x ¨ ) .
They provided a simpler proof for the nonchaotic behavior. The algebraic criterion proved in their work is quite general and can be used to study the nonchaotic behavior of other types of ordinary differential equations.
In [19], M. Ismail et al. dealt with an initial value problem of a nonlinear third order jerk equation:
x + f ( x , x ˙ , x ¨ ) = 0 , x ( 0 ) = 0 , x ˙ ( 0 ) = A , x ¨ ( 0 ) = 0 .
The authors succeeded in extending the global error minimization method (GEMM) to obtain analytic approximations. Compared and simulated with the known solutions and the exact numerical ones, their obtained method were proven to be effective and to provide an efficient alternative to the previously known existing methods.
In nature, most nonlinear systems are mutually coupled. Coupled nonlinear systems have rich dynamic behaviors. As an extension of jerk equations, coupled jerk equations also have simple algebraic structures and more complex dynamic characteristics compared with jerk systems.
In [20], Chen et al. considered the following coupled systems of jerk equations:
x = a x ¨ x ˙ + ( | x | 1 ) + ϵ ( x y ) , y = b y ¨ y ˙ + y y 3 + ϵ ( y x ) .
The authors investigated the dynamical evolution of the above coupled system and obtained bifurcation sets in parameter space by the analysis of the equilibrium points and their stabilities.
During the last two decades, great interest has been devoted to the study of fractional differential equations, which can serve as an excellent tool for the mathematical modeling of systems and processes in the fields of physics, chemistry, biology, electromagnetic, mechanics, economics, dynamical processes, etc. For more details regarding fractional differential equations involving initial or boundary conditions, see for instance, [21,22,23,24,25,26,27,28,29,30,31] and the references therein.
At the same time, the theory of fractional jerk equations has been also analytically investigated by some very interesting and novel papers (see [32,33,34]). In 2020, Echenausía-Monroy et al. [32] considered a multi-scroll generator system based on fractional jerk equations:
D 0 + q x x ( t ) = y ( t ) , D 0 + q y y ( t ) = z ( t ) , D 0 + q z z ( t ) = g ( x , y , z ) ,
where g ( x , y , z ) = α [ x + y + z f ( x ) ] , D 0 + q i denotes the Caputo fractional derivatives i = { x , y , z } , and α R is the control parameter. The authors further studied the effects of the above fractional jerk system and provided a physical interpretation based on statistical analysis. Compared to the integer-order system, their results show that the use of fractional-order can decrease the size of the generated attractor and can modify the long-range correlations in the system. The obtained results can not only be used for applications arising in engineering and the sciences, such as mobile surveillance devices, they can also enrich the analysis and understanding of the implications of fractional integration orders.
From the existing results, we can see a fact: although the solutions for fractional jerk equations have been studied by some authors, to the best of our knowledge, coupled jerk differential equations with the fractional order, have not been investigated until now. The objective of this paper is to fill that gap in the relevant literature.
It is easy to check that the following coupled fractional jerk equations:
D 0 + γ ( D 0 + β ( D 0 + α u ) ) ( t ) = f v ( t ) , D 0 + λ v ( t ) , ( D 0 + μ ( D 0 + λ v ) ) ( t ) , D 0 + ν ( D 0 + μ ( D 0 + λ v ) ) ( t ) = g u ( t ) , D 0 + α u ( t ) , ( D 0 + β ( D 0 + α u ) ) ( t ) , u ( 0 ) = i = 1 m 1 a i u ( ξ i ) , D 0 + α u ( 0 ) = D 0 + β ( D 0 + α u ) ( 0 ) = 0 , v ( 0 ) = i = 1 m 1 b i v ( η i ) , D 0 + λ v ( 0 ) = D 0 + μ ( D 0 + λ v ) ( 0 ) = 0 ,
are equivalent to (1) and (2). Due to the conditions (2), BVP (3) happens to be at resonance in the sense that the associated linear homogeneous equations:
D 0 + γ ( D 0 + β ( D 0 + α u ) ) ( t ) = 0 , D 0 + ν ( D 0 + μ ( D 0 + λ v ) ) ( t ) = 0 , u ( 0 ) = i = 1 m 1 a i u ( ξ i ) , D 0 + α u ( 0 ) = D 0 + β ( D 0 + α u ) ( 0 ) = 0 , v ( 0 ) = i = 1 m 1 b i v ( η i ) , D 0 + λ v ( 0 ) = D 0 + μ ( D 0 + λ v ) ( 0 ) = 0 ,
have ( u , v ) = ( c , d ) , c , d R as nontrivial solutions.
In addition, we remark that if all the fractional orders α , β , γ , λ , μ , and ν are equal to 1, then Equation (3) can be rewritten as
u ( t ) = f ( v , v ˙ , v ¨ ) , v ( t ) = g ( u , u ˙ , u ¨ ) , u ( 0 ) = i = 1 m 1 a i u ( ξ i ) , u ˙ ( 0 ) = u ¨ ( 0 ) = 0 , v ( 0 ) = i = 1 m 1 b i v ( η i ) , v ˙ ( 0 ) = v ¨ ( 0 ) = 0 ,
which is a standard nonlinear coupled jerk system and is a complement of [20]. Compared with [20], we can see the nonhomogeneous terms f and g in the above system are more general. Furthermore, if f = g and u = v , then the above problem can be written by u ( t ) = f ( u , u ˙ , u ¨ ) , u ( 0 ) = i = 1 m 1 a i u ( ξ i ) , u ˙ ( 0 ) = u ¨ ( 0 ) = 0 . Compared with previous research works, such as [19] whose conditions include x ( 0 ) = 0 , we find that the boundary value conditions, u ( 0 ) = i = 1 m 1 a i u ( ξ i ) , are nonlocal initial conditions, which were first used by Byszewski [35] and are more appropriate to describe natural phenomena for they contain more additional information.
As far as nonlinear jerk equations are concerned, most of the efforts on this topic are related to the initial value problem. For example, see [15,16] and the references therein. Among the existing results, no one result can be applied to our problem. Therefore, our results enrich the existing literature. This is another reason why we study problem (1).
This paper is organized as follows. In Section 2, we give some notations and lemmas. In Section 3, we establish a theorem of existence of a solution for the problem (1) and (2) by applying the coincidence degree theory due to Mawhin [36]. In Section 4, we give an example to demonstrate our results.

2. Preliminaries

In this section, we present the necessary definitions and lemmas from fractional calculus theory that can be found in the recent literature [30].
Definition 1
([30]). The Riemann-Liouville fractional integral of order α > 0 of a function f : ( 0 , ) R is given by
I 0 + α f ( t ) = 1 Γ ( α ) 0 t ( t s ) α 1 f ( s ) d s ,
provided that the right-hand side is pointwise defined on ( 0 , ) .
Definition 2
([30]). The Caputo fractional derivative of order α > 0 for a function f C n [ 0 , + ) is given by
D 0 + α f ( t ) = 1 Γ ( n α ) 0 t ( t s ) n α 1 f ( n ) ( s ) d s , n 1 < α < n ,
where n = [ α ] + 1 , [ α ] denotes the integer part of α.
Lemma 1
([30]). Let n 1 < α n , D 0 + α u L 1 ( 0 , 1 ) , then
I 0 + α D 0 + α u ( t ) = u ( t ) + c 0 + c 1 t + + c n 1 t n 1 ,
where c i R , i = 0 , 1 , , n 1 .
Lemma 2
([30]). If β > 0 , α + β > 0 , then the equation
I 0 + α I 0 + β f ( x ) = I 0 + α + β f ( x ) ,
is satisfied for a continuous function f.
Now, let us recall notation about the coincidence degree continuation theorem [36].
Definition 3
([36]). Let Y and Z be normed spaces. A linear operator L : d o m ( L ) Y Z is said to be a Fredholm operator of index zero provided that
(i)
Im L is a closed subset of Z;
(ii)
dimker L = codimIm L < + .
It follows from Definition 3, if L is a Fredholm operator of index zero, then there exist continuous projectors P : Y Y , Q : Z Z as continuous projectors such that ker L = Im P , Im L = ker Q and Y = ker L ker P , Z = Im L Im Q . It follows that L | dom L ker P : dom L ker P Im L is invertible. We denote the inverse of this map by K P . If Ω is an open bounded subset of Y, the map N will be called L-compact on Ω ¯ if Q N ( Ω ¯ ) is bounded and K P , Q N = K P ( I Q ) N : Ω ¯ Y is compact.
Theorem 1
([36]). Let L be a Fredholm operator of index zero and N be L-compact on Ω ¯ . Suppose that the following conditions are satisfied:
(1)
L x λ N x for each ( x , λ ) ( dom L ker L ) Ω × ( 0 , 1 ) ;
(2)
N x Im L for each x ker L Ω ; and
(3)
deg ( J Q N | ker L , Ω ker L , 0 ) 0 , where Q : Z Z is a continuous projection as above with Im L = ker Q and J : Im Q ker L is any isomorphism.
Then, the equation L x = N x has at least one solution in dom L Ω ¯ .

3. Main Results

In this section, we will prove the existence and uniqueness results for (1) and (2).
We use the Banach space E = C [ 0 , 1 ] with the norm u = max 0 t 1 | u ( t ) | . Define two linear spaces
X = u ( t ) | u ( t ) E , D 0 + α u ( t ) E , ( D 0 + β ( D 0 + α u ) ) ( t ) E , Y = v ( t ) | v ( t ) E , D 0 + λ u ( t ) E , ( D 0 + μ ( D 0 + λ u ) ) ( t ) E .
Clearly, X and Y are Banach spaces with the norm
u X = max u , D 0 + α u ( t ) , ( D 0 + β ( D 0 + α u ) ) ( t ) , v Y = max v , D 0 + λ v ( t ) , ( D 0 + μ ( D 0 + λ v ) ) ( t ) .
We consider the Banach space X × Y endowed with the norm defined by ( u , v ) X × Y = max { u X , v Y } and Z = E × E is a Banach space with the norm defined by ( x , y ) Z = max { x , y } .
We define the linear operator L 1 from dom L 1 X to E by
L 1 u = D 0 + γ ( D 0 + β ( D 0 + α u ) ) ( t ) ,
where dom L 1 = { u X | u ( 0 ) = i = 1 m 1 a i u ( ξ i ) , D 0 + α u ( 0 ) = D 0 + β ( D 0 + α u ) ( 0 ) = 0 } .
We define the linear operator L 2 from dom L 2 Y to E by
L 2 v = D 0 + ν ( D 0 + μ ( D 0 + λ v ) ) ( t ) ,
where dom L 2 = { v Y | v ( 0 ) = i = 1 m 1 b i v ( η i ) , D 0 + λ v ( 0 ) = D 0 + μ ( D 0 + λ v ) ( 0 ) = 0 } .
We define the operator L : dom L ( X × Y ) Z by
L ( u , v ) = ( L 1 u , L 2 v ) ,
where dom L = { ( u , v ) X × Y | u dom L 1 , v dom L 2 } , and we define N : X × Y Z by setting
N ( u , v ) = ( N 1 v , N 2 u ) ,
where N 1 : Y E is defined by
N 1 v ( t ) = f v ( t ) , D 0 + λ v ( t ) , ( D 0 + μ ( D 0 + λ v ) ) ( t )
and N 2 : X E is defined by
N 2 u ( t ) = g u ( t ) , D 0 + α u ( t ) , ( D 0 + β ( D 0 + α u ) ) ( t ) .
Then, the problem (1) and (2) can be written by L ( u , v ) = N ( u , v ) .
Lemma 3.
Let the mapping L be defined by (4), then
Ker L = ( u , v ) X × Y : ( u , v ) = ( c , d ) , c , d R ,
Im L = ( x , y ) E : i = 1 m 1 a i I 0 + α + β + γ x ( ξ i ) = i = 1 m 1 b i I 0 + λ + μ + ν y ( η i ) = 0 .
Proof. 
Clearly, L 1 and L 2 are linear operators. By L 1 u = 0 , i.e., D 0 + γ ( D 0 + β ( D 0 + α u ) ) ( t ) = 0 and Lemma 1, we have
( D 0 + β ( D 0 + α u ) ) ( t ) = c 0 , c 0 R .
Applying the operator of I 0 + β on both sides of the above equation, we have
D 0 + α u ( t ) = I 0 + β c 0 + c 1 c 0 , c 1 R .
By Lemma 2, we obtain
u ( t ) = I 0 + β + α c 0 + I 0 + α c 1 + c , c 0 , c 1 , c R .
In view of D 0 + α u ( 0 ) = D 0 + β ( D 0 + α u ) ( 0 ) = 0 , we obtain c 0 = c 1 = 0 . Then, we have u ( t ) = c . By similar proof, if L 2 v = 0 , then we have v ( t ) = d , d R . Hence, one has that (5) holds. It is clear that Ker L = ( c , d ) R 2 .
Next, we prove that (6) holds. Let x Im L 1 , and thus there exists u dom L 1 such that x ( t ) = D 0 + γ ( D 0 + β ( D 0 + α u ) ) ( t ) . By Lemma 1 and the definition of dom L 1 , we have
u ( t ) = I 0 + α + β + γ x ( t ) + I 0 + β + α c 0 + I 0 + α c 1 + c .
In view of D 0 + α u ( 0 ) = D 0 + β ( D 0 + α u ) ( 0 ) = 0 , we have c 0 = c 1 = 0 . Hence,
u ( t ) = I 0 + α + β + γ x ( t ) + c .
According to u ( 0 ) = i = 1 m 1 a i u ( ξ i ) and i = 1 m 1 a i = 1 , we have
i = 1 m 1 a i I 0 + α + β + γ x ( ξ i ) = 0 .
On the other hand, suppose x satisfies the above equation. Letting u ( t ) = I 0 + α + β + γ x ( t ) , we can prove u ( t ) dom L 1 and L 1 u ( t ) = x .
Similarly, we can obtain
i = 1 m 1 b i I 0 + λ + μ + ν y ( η i ) = 0 .
Thus, (6) holds. The proof is complete. □
Lemma 4.
Let L be defined by (4), then L is a Fredholm operator of index zero. Furthermore, the operator K P : Im L dom L Ker P can be written by
K P ( x , y ) = I 0 + α + β + γ x , I 0 + λ + μ + ν y .
Proof. 
Define the operators P 1 : X X , P 2 : Y Y and P : ( u , v ) ( P 1 u , P 2 v ) , as
P 1 u = u ( 0 ) , P 2 v = v ( 0 ) .
Clearly, Ker L = Im P and P ( u , v ) = P 2 ( u , v ) .
Note that
Ker P = ( u , v ) | u ( 0 ) = 0 , v ( 0 ) = 0 .
Since ( u , v ) = ( u , v ) P ( u , v ) + P ( u , v ) , it is clear that X × Y = Ker P + Ker L . By a simple calculation, we have Ker L Ker P = { ( 0 , 0 ) } . Thus, we obtain
X × Y = Ker L Ker P .
Consider the linear operators Q 1 , Q 2 : E E defined by
Q 1 x ( t ) = Γ ( 1 + α + β + γ ) i = 1 m 1 a i ξ i i = 1 m 1 a i I 0 + α + β + γ x ( ξ i ) , Q 2 y ( t ) = Γ ( 1 + λ + μ + ν ) i = 1 m 1 b i η i i = 1 m 1 b i I 0 + λ + μ + ν y ( η i ) .
Clearly, Q ( x , y ) = Q 1 x ( t ) , Q 2 y ( t ) R 2 . Take x ( t ) E , by a direct computation, we have that
Q 1 ( Q 1 x ( t ) ) = Q 1 x ( t ) · Γ ( 1 + α + β + γ ) i = 1 m 1 a i ξ i i = 1 m 1 a i ( I 0 + α + β + γ 1 ) ( t = ξ i ) = Q 1 x ( t ) .
Similarly, Q 2 2 = Q 2 . This gives that Q 2 ( x , y ) = Q ( x , y ) . It is easy to check from ( x , y ) = ( x , y ) Q ( x , y ) + Q ( x , y ) that Z = Im L + Im Q . Take ( u , v ) Im L Im Q . According to the definitions of Im L and Im Q , we obtain ( u , v ) = ( 0 , 0 ) . This implies that
Z = Im L Im Q .
Now, Ind L = dim ker L codim Im L = 0 . This means that L is a Fredholm mapping of index zero.
Next, we will prove that K P is the inverse of L dom L Ker P .
For ( x , y ) Im L , we have
L K P ( x , y ) = L ( I 0 + α + β + γ x , I 0 + λ + μ + ν y ) = ( D 0 + α + β + γ I 0 + α + β + γ x , D 0 + λ + μ + ν I 0 + λ + μ + ν y ) = ( x , y ) .
For ( u , v ) dom L Ker P , according to the definitions of dom L and Ker P , it is easy to verify that the constants c i , d i , i = 0 , 1 , 2 in the following equations
u ( t ) = I 0 + α + β + γ x ( t ) + I 0 + β + α c 0 + I 0 + α c 1 + c 2 , v ( t ) = I 0 + λ + μ + ν y ( t ) + I 0 + μ + λ d 0 + I 0 + λ d 1 + d 2 ,
are all equal to zero. Thus, we obtain
K p L ( x , y ) = I 0 + α + β + γ D 0 + γ ( D 0 + β ( D 0 + α x ) ) , I 0 + λ + μ + ν D 0 + ν ( D 0 + μ ( D 0 + λ y ) ) = ( x , y ) .
That shows that K P = ( L dom L Ker P ) 1 . The proof is complete. □
To simplify our statement, we write
a = 1 Γ ( γ + 1 ) , b = 1 Γ ( ν + 1 ) .
For every ( u , v ) X × Y ,
P ( u , v ) X × Y = ( P 1 u , P 2 v ) X × Y = max P 1 u X ; P 2 v Y = max { u ( 0 ) X ; v ( 0 ) Y } = max { | u ( 0 ) | ; | v ( 0 ) | } .
For each ( x , y ) Im L , we have
K P ( x , y ) X × Y = ( I 0 + α + β + γ x , I 0 + λ + μ + ν y ) X × Y = max I 0 + α + β + γ x X ; I 0 + λ + μ + ν y Y max { max I 0 + α + β + γ x , D 0 + α I 0 + α + β + γ x , ( D 0 + β ( D 0 + α I 0 + α + β + γ x ) ) ; max I 0 + λ + μ + ν y , D 0 + λ I 0 + λ + μ + ν y , ( D 0 + μ ( D 0 + λ I 0 + λ + μ + ν y ) ) } = max { max I 0 + α + β + γ x , I 0 + β + γ x , I 0 + γ x ; max I 0 + λ + μ + ν y , I 0 + μ + ν y , I 0 + ν y } = max 1 Γ ( γ + 1 ) x ; 1 Γ ( ν + 1 ) y = max a x ; b y .
With the similar arguments to [21], we obtain the following lemma.
Lemma 5.
K P ( I Q ) N : Y Y is completely continuous.
To obtain our main results, we need the following conditions.
( H 1 ) There exist positive continuous functions a i ( t ) , b i ( t ) , i = 1 , 2 , 3 such that for all ( x 1 , x 2 , x 3 ) R 3 , one has
| f ( x 1 , x 2 , x 3 ) | a 1 ( t ) | x 1 | + a 2 ( t ) | x 2 | + a 3 ( t ) | x 3 | , | g ( x 1 , x 2 , x 3 ) | b 1 ( t ) | x 1 | + b 2 ( t ) | x 2 | + b 3 ( t ) | x 3 | .
( H 2 ) There exists a constant A > 0 such that, if for any t [ 0 , 1 ] , | u ( t ) | > A or | v ( t ) | > A , then
u ( t ) · Q 1 ( N 1 v ) > 0 , v ( t ) · Q 2 ( N 2 u ) > 0
or
u ( t ) · Q 1 ( N 1 v ) < 0 , v ( t ) · Q 2 ( N 2 u ) < 0 .
( H 3 ) max p + a i = 1 3 p i , q + b i = 1 3 q i , p + b i = 1 3 q i , q + a i = 1 3 p i < 1 , where p i = a i , q i = b i , i = 1 , 2 , 3 , p = p 1 + p 2 + p 3 Γ α + β + γ + 1 , q = q 1 + q 2 + q 3 Γ λ + μ + ν + 1 .
Lemma 6.
Ω 1 = ( u , v ) dom L Ker L : L ( u , v ) = λ N ( u , v ) , λ [ 0 , 1 ] is bounded.
Proof. 
For ( u , v ) Ω 1 , thus λ 0 . And L ( u , v ) = λ N ( u , v ) Im L = Ker Q , according to the definition of Ker Q , we have
Q 1 ( N 1 v ) = Q 2 ( N 2 u ) = 0 .
By (H 2 ) , there exit t 0 , t 1 [ 0 , 1 ] such that
| u ( t 0 ) | A , | v ( t 1 ) | A .
Again for ( u , v ) Ω 1 , ( u , v ) dom ( L ) Ker ( L ) , then ( I P ) ( u , v ) dom L Ker P and L P ( u , v ) = ( 0 , 0 ) ; thus, from (8), we have
( I P ) ( u , v ) X × Y = K P L ( I P ) ( u , v ) X × Y = K P ( L 1 u , L 2 v ) X × Y max { a L 1 v ; b L 2 u } = max { a N 1 v ; b N 2 u } .
By L u = λ N u and u dom L , we have ( L 1 u , L 2 v ) = ( λ N 1 v , λ N 2 u ) , i.e., L 1 u = λ N 1 v and L 2 v = λ N 2 u . Thus,
u ( t ) = λ I 0 + α + β + γ N 1 v ( t ) + I 0 + β + α c 0 + I 0 + α c 1 + c 2 , v ( t ) = λ I 0 + λ + μ + ν N 2 u ( t ) + I 0 + μ + λ d 0 + I 0 + λ d 1 + d 2 .
By the definition of dom L , we have c i = d i = 0 , i = 0 , 1 . Then,
u ( t ) = λ I 0 + α + β + γ N 1 v ( t ) + c 2 , v ( t ) = λ I 0 + λ + μ + ν N 2 u ( t ) + d 2 .
Furthermore, we have
u ( t 0 ) u ( 0 ) = λ I 0 + α + β + γ N 1 v ( t 0 ) + c 2 λ I 0 + α + β + γ N 1 v ( 0 ) + c 2 = I 0 + α + β + γ N 1 v ( t 0 ) , v ( t 1 ) v ( 0 ) = λ I 0 + λ + μ + ν N 2 u ( t 1 ) + d 2 λ I 0 + λ + μ + ν N 2 u ( 0 ) + d 2 = I 0 + λ + μ + ν N 2 u ( t 1 ) ,
that are,
u ( 0 ) = u ( t 0 ) λ I 0 + α + β + γ N 1 v ( t 0 ) , v ( 0 ) = v ( t 1 ) λ I 0 + λ + μ + ν N 2 u ( t 1 ) .
Together with | u ( t 0 ) | A , we can derive that
| u ( 0 ) | | u ( t 0 ) | + λ I 0 + α + β + γ N 1 v ( t 0 ) A + λ Γ α + β + γ 0 t 0 t 0 s α + β + γ 1 f v ( s ) , D 0 + λ v ( s ) , ( D 0 + μ ( D 0 + λ v ) ) ( s ) d s A + 1 Γ α + β + γ 0 t 0 t 0 s α + β + γ 1 f v ( s ) , D 0 + λ v ( s ) , ( D 0 + μ ( D 0 + λ v ) ) ( s ) d s A + 1 Γ α + β + γ 0 t 0 t 0 s α + β + γ 1 · a 1 ( s ) | v ( s ) | + a 2 ( s ) | D 0 + λ v ( s ) | + a 3 ( s ) | ( D 0 + μ ( D 0 + λ v ) ) ( s ) | d s A + 1 Γ α + β + γ p 1 v ( t ) + p 2 D 0 + λ v ( t ) + p 3 ( D 0 + μ ( D 0 + λ v ) ) ( t ) · 0 t 0 t 0 s α + β + γ 1 d s A + p 1 + p 2 + p 3 Γ α + β + γ 0 t 0 t 0 s α + β + γ 1 d s · v ( t ) Y = A + p v ( t ) Y .
With the similar arguments, we obtain
| v ( 0 ) | A + q 1 + q 2 + q 3 Γ λ + μ + ν + 1 u ( t ) X = A + q u ( t ) X .
From (7) and (9), we have
( u , v ) X × Y = P ( u , v ) + ( I P ) ( u , v ) X × Y P ( u , v ) X × Y + ( I P ) ( u , v ) X × Y max { | u ( 0 ) | + a N 1 v , | u ( 0 ) | + b N 2 u , | v ( 0 ) | + a N 1 v , | v ( 0 ) | + b N 2 u } .
In what follows, the proof can be divided into four cases.
Case 1.
( u , v ) X × Y u ( 0 ) + a N 1 v .
By (10) and (H 1 ), we have
( u , v ) X × Y | u ( 0 ) | + a N 1 v A + p v ( t ) Y + a N 1 v = A + p v ( t ) Y + a f v ( t ) , D 0 + λ v ( t ) , ( D 0 + μ ( D 0 + λ v ) ) ( t ) A + p v ( t ) Y + a a 1 ( t ) | v ( t ) | + a 2 ( t ) | D 0 + λ v ( t ) | + a 3 ( t ) | ( D 0 + μ ( D 0 + λ v ) ) ( t ) | A + p v ( t ) Y + a p 1 v ( t ) + p 2 D 0 + λ v ( t ) + p 3 ( D 0 + μ ( D 0 + λ v ) ) ( t ) = A + p v ( t ) Y + a ( p 1 + p 2 + p 3 ) v ( t ) Y = A + p + a Σ i = 1 3 p i · v ( t ) Y A + p + a Σ i = 1 3 p i · ( u , v ) X × Y .
By (H 4 ) , we have
( u , v ) X × Y A 1 p a Σ i = 1 3 p i : = M .
Therefore, Ω 1 is bounded.
Case 2.
( u , v ) X × Y | u ( 0 ) | + b N 2 u .
From (10) and (H 1 ), we obtain
( u , v ) X × Y | u ( 0 ) | + b N 2 u A + p v ( t ) Y + b N 2 u = A + p v ( t ) Y + b g u ( t ) , D 0 + α u ( t ) , ( D 0 + β ( D 0 + α u ) ) ( t ) A + p v ( t ) Y + b q 1 u ( t ) + q 2 D 0 + λ u ( t ) + q 3 ( D 0 + μ ( D 0 + λ u ) ) ( t ) A + p v ( t ) Y + b ( q 1 + q 2 + q 3 ) u ( t ) Y A + p + b Σ i = 1 3 q i · ( u , v ) X × Y .
By (H 4 ), we can obtain
( u , v ) X × Y 1 1 q b Σ i = 1 3 q i : = M .
Therefore, Ω 1 is bounded.
Case 3.
( u , v ) | v ( 0 ) | + a N 1 v .
According to (11) and (H 4 ) , by a similar proof of Case 2, we can derive
( u , v ) X × Y A 1 q a Σ i = 1 3 p i : = M .
Therefore, Ω 1 is bounded.
Case 4.
( u , v ) Y | v ( 0 ) | + b N 2 u .
By a similar proof of Case 1, we obtain
( u , v ) X × Y A 1 q b Σ i = 1 3 q i : = M .
Therefore, Ω 1 is bounded.
According the above arguments, we prove that Ω 1 is bounded. □
Lemma 7.
Ω 2 = ( u , v ) Ker L : N ( u , v ) Im L is bounded.
Proof. 
Let ( u , v ) Ker L ; thus, we have u = c , v = d , c , d R . In view of N ( u , v ) = ( N 1 v , N 2 u ) Im L = Ker Q , we have
Q 1 ( N 1 v ) = Q 2 ( N 2 u ) = 0 .
By (H 2 ), there exit constants t 0 , t 1 [ 0 , 1 ] such that
| u ( t 0 ) | = | c | A , | v ( t 1 ) | = | d | A .
Hence, Ω 2 is bounded. □
Lemma 8.
Ω 3 = ( u , v ) Ker L : λ ( u , v ) + ( 1 λ ) Q N ( u , v ) = ( 0 , 0 ) , λ [ 0 , 1 ] is bounded.
Proof. 
Let ( u , v ) Ker L ; therefore, we have u = c , v = d , c , d R and
λ c + ( 1 λ ) Q 1 N 1 ( d ) = 0 , λ d + ( 1 λ ) Q 2 N 2 ( c ) = 0 .
If λ = 0 , similar to the proof of Lemma 7, we can obtain | c | A and | d | A . If λ = 1 , we have c = d = 0 . If λ ( 0 , 1 ) , then we can also obtain | c | A and | d | A . Otherwise, if | c | > A or | d | > A , in view of the first part of (H 2 ), we obtain
λ c 2 + c ( 1 λ ) Q 1 N 1 ( d ) = 0 , λ d 2 + d ( 1 λ ) Q 2 N 2 ( c ) = 0 ,
which contradicts (12). Thus, Ω 3 is bounded.
If the second part of (H 2 ) holds, then we can prove the set
Ω 3 = ( u , v ) Ker L : λ ( u , v ) + ( 1 λ ) Q N ( u , v ) = ( 0 , 0 ) , λ [ 0 , 1 ]
is bounded. □
Theorem 2.
Suppose (H 1 )–(H 3 ) hold, then the problem (1) and (2) has at least one solution in Y.
Proof. 
Let Ω be a bounded open subset of Y, such that i = 1 3 Ω ¯ i Ω . It follows from Lemma 5, that N is L-compact on Ω. By Lemmas 6–8, we obtain:
(1)
L ( u , v ) λ N ( u , v ) , for every ( u , v , λ ) [ ( dom L Ker L ) Ω ] × ( 0 , 1 ) .
(2)
N ( u , v ) Im L for every u Ker L Ω .
(3)
Let H ( u , v ) , λ = ± λ I ( u , v ) + ( 1 λ ) J Q N ( u , v ) . Here, we let I and the isomorphism J : Im Q ker L , which are both identical operators. Via the homotopy property of degree, we obtain that
deg J Q N | ker L , Ω ker L , 0 = deg H ( · , 0 ) , Ω ker L , 0 = deg H ( · , 1 ) , Ω ker L , 0 = deg I , Ω ker L , 0 = 1 0 .
Applying Theorem 1, we conclude that L ( u , v ) = N ( u , v ) has at least one solution in dom L Ω ¯ . □

4. Example

Let us consider the following coupled system of fractional jerk equations at resonance
D 0 + 0.1 x 1 ( t ) = x 2 ( t ) , D 0 + 0.3 x 2 ( t ) = x 3 ( t ) , D 0 + 0.5 x 3 ( t ) = f y 1 , y 2 , y 3 , D 0 + 0.7 y 1 ( t ) = y 2 ( t ) , D 0 + 0.4 y 2 ( t ) = y 3 ( t ) , D 0 + 0.2 y 3 ( t ) = g x 1 , x 2 , x 3 , x 1 ( 0 ) = i = 1 3 a i x 1 ( ξ i ) , x 2 ( 0 ) = x 3 ( 0 ) = 0 , y 1 ( 0 ) = i = 1 3 b i y 1 ( η i ) , y 2 ( 0 ) = y 3 ( 0 ) = 0 ,
where
f ( y 1 , y 2 , y 3 ) = | y 1 | 10 + sin y 2 20 + cos y 3 20 ,
g ( x 1 , x 2 , x 3 ) = sin 2 x 1 16 + sin x 2 32 + | x 3 | 32 ,
and a 1 = a 2 = 1 4 , a 3 = 1 2 , b 1 = b 2 = b 3 = 1 3 , ξ 1 = 0.3 , ξ 2 = 0.5 , ξ 3 = 0.7 , η 1 = 0.2 , η 2 = 0.4 , and η 1 = 0.6 . Corresponding to BVP (1) and (2), we have that α = 0.1 , β = 0.5 , γ = 0.3 , λ = 0.7 , μ = 0.4 , and ν = 0.2 .
It is easy to see that
p 1 = 1 10 , = p 2 = 1 20 , p 3 = 1 20 , a = 1 Γ ( γ + 1 ) = 1 Γ ( 1.3 ) 1.11 q 1 = 1 16 , q 2 = 1 32 , q 3 = 1 16 , b = 1 Γ ( ν + 1 ) = 1 Γ ( 1.2 ) 1.09 .
By direct calculation, we can obtain
p = p 1 + p 2 + p 3 Γ α + β + γ + 1 = 0.2 Γ ( 1.9 ) 0.21 , p + a i = 1 3 p i 0.65 , q + b i = 1 3 q i 0.28 , q = q 1 + q 2 + q 3 Γ λ + μ + ν + 1 = 0.125 Γ ( 2.3 ) 0.11 , p + b i = 1 3 q i 0.38 , q + a i = 1 3 p i 0.55 .
We take A = 18 , then we can obtain that (H 1 )–(H 3 ) hold. Hence, all the conditions of Theorem 2 are satisfied, and consequently BVP (13) has at least one solution.

5. Conclusions

In this paper, we investigated a resonant boundary value problem of the system of jerk differential equations with the fractional order. The interesting point is that two fractional jerk equations are coupled. By coincidence degree theory due to Mawhin, the existence result is proved. Our result obtained in this paper is new and complements the existing literature on the topic. As far as our work is concerned, we mainly concentrated on the existence of solutions. To the best of our knowledge, some results were not considered for fractional coupled jerk equations with resonant conditions, such as unique solutions, positive solutions, and numerical solutions. In future research, we will study the corresponding problem, and we hope to be able to make some progress. The corresponding physical interpretation of fractional coupled jerk equations, in order to compare the presented results on integer-order jerk system, is proposed as future work.

Author Contributions

All authors contributed equally in writing this article. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China under Grant number 61803230, 12071302; the Fundamental Research Funds for the Central Universities (2020YJSLX01); Shandong Provincial Natural Science Foundation (ZR2020MG001); a Project of Shandong Province Higher Educational Science and Technology Program under Grant number (No. J17KA073, J18KA330) and University outstanding youth innovation team development plan of Shandong Province under Grant number 2019KJN023. This research was also supported by the Research Fund for the Doctoral Program of Shandong Jiaotong University, Shandong Jiaotong University Climbing Research Innovation Team Program (No. SDJTUC1803).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schot, S.H. Jerk: The time rate of change of acceleration. Am. J. Phys. 1978, 46, 1090–1094. [Google Scholar] [CrossRef]
  2. Rothbart, H.A.; Wahl, A.M. Mechanical Design and Systems Handbook. J. Appl. Mech. 1965, 32, 478. [Google Scholar] [CrossRef] [Green Version]
  3. El-Nabulsi, R. Nonlocal-in-time kinetic energy in nonconservative fractional systems, disordered dynamics, jerk and snap and oscillatory motions in the rotating fluid tube. Int. J. Nonlinear Mech. 2017, 93, 65–81. [Google Scholar] [CrossRef]
  4. Tedesco, L. Ellipsoidal expansion of the Universe, cosmic shear, acceleration and jerk parameter. Eur. Phys. J. Plus 2018, 133, 1–9. [Google Scholar] [CrossRef] [Green Version]
  5. El-Nabulsi, R. Jerk in Planetary Systems and Rotational Dynamics, Nonlocal Motion Relative to Earth and Nonlocal Fluid Dynamics in Rotating Earth Frame. Earth Moon Planets 2018, 122, 15–41. [Google Scholar] [CrossRef]
  6. El-Nabulsi, R. Time-nonlocal kinetic equations, jerk and hyperjerk in plasmas and solar physics. Adv. Space Res. 2018, 61, 2914–2931. [Google Scholar] [CrossRef]
  7. El-Nabulsi, R. Free variable mass nonlocal systems, jerks, and snaps, and their implications in rotating fluids in rockets. Acta Mech. 2021, 232, 89–109. [Google Scholar] [CrossRef]
  8. Eager, D.; Pendrill, A.; Reistad, N. Beyond velocity and acceleration: Jerk, snap and higher derivatives. Eur. J. Phys. 2016, 37, 1–11. [Google Scholar] [CrossRef]
  9. Gómez-Aguilar, J.F.; Rosales-García, J.; Escobar-Jiménez, R.F.; López-López, M.G.; Alvarado-Martínez, V.M.; Olivares-Peregrino, V.H. On the Possibility of the Jerk Derivative in Electrical Circuits. Adv. Math. Phys. 2016, 2016, 1–8. [Google Scholar] [CrossRef] [Green Version]
  10. Faires, V.M. Design of Machine Elements, 4th ed.; Macmillan: New York, NY, USA, 1965. [Google Scholar]
  11. Linz, S.J. Nonlinear dynamical models and jerky motion. Am. J. Phys. 1997, 65, 523–526. [Google Scholar] [CrossRef]
  12. Heidel, J.; Zhang, F. Nonchaotic behavior in three-dimensional quadratic systems II. The conservative case. Nonlinearity 1999, 12, 617–633. [Google Scholar] [CrossRef] [Green Version]
  13. Yang, X.S. On non-chaotic behavior of a class of jerk systems. Far East J. Dyn. Syst. 2002, 4, 27–38. [Google Scholar]
  14. Ma, X.; Wei, L.; Guo, Z. He’s homotopy perturbation method to periodic solutions of nonlinear Jerk equations. J. Sound Vib. 2008, 314, 217–227. [Google Scholar] [CrossRef]
  15. Rahman, M.S.; Hasan, A.M. Modified harmonic balance method for the solution of nonlinear jerk equations. Results Phys. 2018, 8, 893–897. [Google Scholar] [CrossRef]
  16. Gottlieb, H.P.W. Simple nonlinear jerk functions with periodic solutions. Am. J. Phys. 1998, 66, 903–906. [Google Scholar] [CrossRef]
  17. Gottlieb, H.P.W. Harmonic balance approach to limit cycles for nonlinear jerk equations. J. Sound Vib. 2006, 297, 243–250. [Google Scholar] [CrossRef] [Green Version]
  18. Messias, M.; Silva, R.P. Determination of Nonchaotic Behavior for Some Classes of Polynomial Jerk Equations. Int. J. Bifurc. Chaos 2020, 30, 1–12. [Google Scholar] [CrossRef]
  19. Ismail, G.; Abu-zinadah, H.H. Analytic Approximations to Non-linear Third Order Jerk Equations via Modified Global Error Minimization Method. J. King Saud Univ. Sci. 2020, 33, 1–5. [Google Scholar] [CrossRef]
  20. Chen, Z.Y.; Bi, Q.S. Bifurcations and chaos of coupled Jerk systems. Acta Phys. Sin. 2010, 59, 7669–7678. [Google Scholar]
  21. Zhang, Y.L.; Bai, Z. Existence of solutions for nonlinear fractional three-point boundary value problems at resonance. J. Appl. Math. Comput. 2011, 36, 417–440. [Google Scholar] [CrossRef]
  22. Guan, T.; Wang, G. Maximum Principle for the Space-Time Fractional Conformable Differential System Involving the Fractional Laplace Operator. J. Math. 2020, 2020, 1–8. [Google Scholar] [CrossRef]
  23. Zhang, W.; Liu, W.B. Variational approach to fractional Dirichlet problem with instantaneous and non-instantaneous impulses. Appl. Math. Lett. 2020, 99, 1–7. [Google Scholar] [CrossRef]
  24. Cui, Y. Uniqueness of solution for boundary value problems for fractional differential equations. Appl. Math. Lett. 2016, 51, 48–54. [Google Scholar] [CrossRef]
  25. Wang, Y.; Wang, H. Triple positive solutions for fractional differential equation boundary value problems at resonance. Appl. Math. Lett. 2020, 106, 1–7. [Google Scholar] [CrossRef]
  26. Su, X. Boundary value problem for a coupled system of nonlinear fractional differential equations. Appl. Math. Lett. 2009, 22, 64–69. [Google Scholar] [CrossRef] [Green Version]
  27. Sun, B.; Jiang, W. Existence of solutions for functional boundary value problems at resonance on the half-line. Bound. Value Probl. 2020, 2020, 1–15. [Google Scholar] [CrossRef]
  28. Kosmatov, N. A boundary value problem of fractional order at resonance. Electron. J. Differ. Equ. 2010, 2010, 1–10. [Google Scholar]
  29. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  30. Kilbas, A.A.; Srivastava, H.H.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  31. Shah, A.; Khan, R.A.; Khan, A.; Khan, H.; Gómez-Aguilar, J.F. Investigation of a system of nonlinear fractional order hybrid differential equations under usual boundary conditions for existence of solution. Math. Meth. Appl. Sci. 2020, 44, 1–11. [Google Scholar] [CrossRef]
  32. Echenausía-Monroy, J.L.; Gilardi-Velázquez, H.E.; Jaimes-Reátegui, R.; Aboites, V.; Huerta-Cuéllar, G. A physical interpretation of fractional-order-derivatives in a jerk system: Electronic approach. Commun. Nonlinear Sci. Numer. Simul. 2020, 90, 105413. [Google Scholar] [CrossRef]
  33. Liu, C.; Chang, J.R. The periods and periodic solutions of nonlinear jerk equations solved by an iterative algorithm based on a shape function method. Appl. Math. Lett. 2020, 102, 1–9. [Google Scholar] [CrossRef]
  34. Prakash, P.; Singh, J.P.; Roy, B.K. Fractional-order memristor-based chaotic jerk system with no equilibrium point and its fractional-order backstepping control. IFAC-PapersOnLine 2018, 51, 1–6. [Google Scholar] [CrossRef]
  35. Byszewski, L. Theorems about existence and uniqueness of solutions of a semi-linear evolution nonlocal Cauchy problem. J. Math. Anal. Appl. 1991, 162, 494–505. [Google Scholar] [CrossRef] [Green Version]
  36. Mawhin, J. Topological degree and boundary value problems for nonlinear differential equations in topological methods for ordinary differential equations. Lect. Notes Math. 1993, 1537, 74–142. [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

Hu, L.; Han, Y.; Zhang, S. On the Existence of Coupled Fractional Jerk Equations with Multi-Point Boundary Conditions. Axioms 2021, 10, 103. https://doi.org/10.3390/axioms10020103

AMA Style

Hu L, Han Y, Zhang S. On the Existence of Coupled Fractional Jerk Equations with Multi-Point Boundary Conditions. Axioms. 2021; 10(2):103. https://doi.org/10.3390/axioms10020103

Chicago/Turabian Style

Hu, Lei, Yaozhen Han, and Shuqin Zhang. 2021. "On the Existence of Coupled Fractional Jerk Equations with Multi-Point Boundary Conditions" Axioms 10, no. 2: 103. https://doi.org/10.3390/axioms10020103

APA Style

Hu, L., Han, Y., & Zhang, S. (2021). On the Existence of Coupled Fractional Jerk Equations with Multi-Point Boundary Conditions. Axioms, 10(2), 103. https://doi.org/10.3390/axioms10020103

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