Next Article in Journal
Some Novel Inequalities for LR-(k,h-m)-p Convex Interval Valued Functions by Means of Pseudo Order Relation
Next Article in Special Issue
Existence and Ulam Type Stability for Impulsive Fractional Differential Systems with Pure Delay
Previous Article in Journal
New Bright and Kink Soliton Solutions for Fractional Complex Ginzburg–Landau Equation with Non-Local Nonlinearity Term
Previous Article in Special Issue
Magnetic Field, Variable Thermal Conductivity, Thermal Radiation, and Viscous Dissipation Effect on Heat and Momentum of Fractional Oldroyd-B Bio Nano-Fluid within a Channel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability of a Nonlinear Langevin System of ML-Type Fractional Derivative Affected by Time-Varying Delays and Differential Feedback Control

Department of Mathematics, School of Electronics & Information Engineering, Taizhou University, Taizhou 318000, China
Fractal Fract. 2022, 6(12), 725; https://doi.org/10.3390/fractalfract6120725
Submission received: 25 October 2022 / Revised: 30 November 2022 / Accepted: 5 December 2022 / Published: 8 December 2022

Abstract

:
The Langevin system is an important mathematical model to describe Brownian motion. The research shows that fractional differential equations have more advantages in viscoelasticity. The exploration of fractional Langevin system dynamics is novel and valuable. Compared with the fractional system of Caputo or Riemann–Liouville (RL) derivatives, the system with Mittag–Leffler (ML)-type fractional derivatives can eliminate singularity such that the solution of the system has better analytical properties. Therefore, we concentrate on a nonlinear Langevin system of ML-type fractional derivatives affected by time-varying delays and differential feedback control in the manuscript. We first utilize two fixed-point theorems proposed by Krasnoselskii and Schauder to investigate the existence of a solution. Next, we employ the contraction mapping principle and nonlinear analysis to establish the stability of types such as Ulam–Hyers (UH) and Ulam–Hyers–Rassias (UHR) as well as generalized UH and UHR. Lastly, the theoretical analysis and numerical simulation of some interesting examples are carried out by using our main results and the DDESD toolbox of MATLAB.

1. Introduction

This work mostly takes into account a nonlinear Langevin system of ML-type fractional derivatives affected by time-varying delays and differential feedback control as follows:
M L D 0 + β M L D 0 + α λ u ( t ) = f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) , t ( 0 , b ] , w ( t ) = g ( t , u ( t σ 3 ( t ) ) , w ( t σ 4 ( t ) ) ) , t J = [ 0 , b ] , u ( t ) = ω 1 ( t ) , w ( t ) = ω 2 ( t ) , M L D 0 + α u ( t ) = ω 3 ( t ) , t [ σ , 0 ] ,
where 0 < α , β 1 and b , λ > 0 , M L D 0 + κ is the ML-type fractional derivative with κ -order, for all 0 < κ 1 . f , g C ( J × R 2 , R ) , σ l C ( J , R + ) ( l = 1 , 2 , 3 , 4 ) , σ = max 1 l 4 { max t J σ l ( t ) } , ω m C ( [ σ , 0 ] , R ) ( m = 1 , 2 , 3 ) .
The main motivation comes from the important theory and wide application of the Langevin system. In 1908, the French physicist Paul Langevin put forward the famous Langevin equation to expound the random motion of particles annihilated in the fluid due to the collision between particles and fluid molecules. In [1,2], the authors provided a large number of random examples of the Langevin system as mathematical models. The fractional-order Langevin system makes up for the deficiency of the integer-order Langevin system and has been widely used and studied. The integer-order Langevin system can not meet the accuracy requirements in describing complex viscoelasticity. Therefore, the classical Langevin system has been extended and modified. For example, Kubo [3,4] imitated a complex viscoelastic abnormal diffusion process by applying a general Langevin system. It is worth noting that the derivatives in these generalized equations are of integer order. Because fractional derivatives have advantages in describing the process of memory and viscoelasticity, another modification is to replace the integral derivative with a fractional derivative in the Langevin system. The early research on the fractional Langevin system as a model was conducted by Eab and Lim [5] as well as Sandev and Tomovski [6]. They used the fractional Langevin system to study the one-dimensional diffusion process and the motion of free particles driven by power law noise, respectively. Recently, some new fields have emerged in the application of fractional differential systems (see [7,8,9,10]). Additionally, some important achievements have been made in the research of the fractional Langevin system (see [11,12,13,14,15,16,17]).
Furthermore, for a system with practical application background, its stability is very important. According to practical needs, scientists have put forward many concepts of system stability. UH-stability is one of the most important stabilities, which was posed by Ulam and Hyers [18,19] in the 1940s. In the last ten years, the research on the UH-stability of fractional systems has been highly praised by many scholars. There have many papers dealing with the UH-stability [20,21,22,23,24,25,26] and generalized UH-stability [27,28,29] of fractional system. However, only a few papers [30,31,32,33,34] have discussed that the fractional Langevin system is UH-stable.
It is worth noting that the previous works on the fractional Langevin system basically involve Caputo- or RL-fractional derivatives. However, under certain circumstances, the singularity of Caputo or Riemann–Liouville fractional derivatives is unavoidable. This leads to huge errors and even distortion in the application of the Caputo or Riemann–Liouville fractional system in some physical fields. Hence, some scholars began to modify the previous definition of a fractional derivative to eliminate singularity. For instance, Caputo and Febrizio raised a nonsingular fractional derivative (see [35]), which is called the Caputo–Febrizio fractional derivative. In [36], Atangana and Baleanu proposed another novel nonsingular fractional derivative with an ML-kernel. In both theory [37,38,39,40] and application [41,42,43,44,45,46,47], many scholars have focused on the study of these nonsingular fractional systems since they were put forward. However, only rare published works [48,49] have shown that some ML-type fractional Langevin systems are UH-stable.
Inspired by the aforementioned, the main purpose of this manuscript is to find some sufficient conditions for the existence and stability of the solution of system (1). The main contributions of our work are reflected in the following aspects. (i) Comparison with previous works on fractional Langevin equation such as [30,31,32], our system (1) involves nonsingular fractional derivatives with an ML-kernel. (ii) In system (1), we consider the influence of delayed differential feedback control, which is not found in previous papers. (iii) We obtain some new criteria of existence and UH-type stability of solutions. (iv) By applying the DDESD toolbox of MATLAB, the simulations and numerical solutions are given.
In Section 2, we will review some concepts and lemmas used later. Based on some fixed-point theorems, we prove that system (1) has at least a solution in Section 3. The UH-, UHR-, GUH-, and GUHR-stability of (1) are built in Section 4. As applications, we conduct theoretical analysis and numerical simulation on some examples to verify the correctness and effectiveness of our main results in Section 5. Finally, a brief summary is made in Section 6. In addition, we put lengthy and complex theoretical proofs in the appendices to improve the readability and conciseness of the manuscript.

2. Preliminaries

Definition 1
([50]). Let 0 < α 1 , b > 0 , u : [ 0 , b ] R . Define an integral
M L I 0 + α u ( t ) = 1 α N ( α ) u ( t ) + α N ( α ) Γ ( α ) 0 t ( t x ) α 1 u ( x ) d x ,
where Γ ( α ) = 0 x α 1 e x d x , the normalized constant N ( α ) satisfies N ( 0 ) = N ( 1 ) = 1 . If the above integral exists, then it is called the α-order ML-fractional integral along the left of u.
Definition 2
([36]). Let 0 < α 1 , b > 0 and u C 1 ( 0 , b ) . The α-order ML-fractional derivative along the left of u in sense of Caputo is defined by
M L D 0 + α u ( t ) = N ( α ) ( 1 α ) 0 t E α α 1 α ( t x ) α u ( x ) d x ,
where E α ( z ) = k = 0 z k Γ ( α k + 1 ) is is called the one-parameter ML-function.
Lemma 1
([37]). Given 0 < ν 1 and h C [ 0 , b ] , then the following ML-fractional system
M L D 0 + ν w ( t ) = h ( t ) , t ( 0 , b ) , w ( 0 ) = w 0 ,
is uniquely solved by
w ( t ) = w 0 + 1 ν N ( ν ) [ h ( t ) h ( 0 ) ] + ν N ( ν ) Γ ( ν ) 0 t ( t x ) ν 1 h ( x ) d x .
Remark 1.
Definition 2 and Lemma 1 imply that M L D 0 + α u ( t ) 0 iff u ( t ) is identical to a constant.
The whole paper needs the following fundamental assumptions.
( A 1 )
λ , b > 0 and 0 < α , β 1 are given constants such that Δ 1 λ ( 1 α ) N ( α ) 0 .
( A 2 )
f , g C ( J × R 2 , R ) , σ l C ( J , R + ) ( l = 1 , 2 , 3 , 4 ) , σ = max 1 l 4 { max t J σ l ( t ) } , ω m C ( [ σ , 0 ] , R ) ( m = 1 , 2 , 3 ) .
Consider the below nonlinear integral system
u ( t ) = ω 1 ( 0 ) + 1 Δ { ω 3 ( 0 ) λ ω 1 ( 0 ) N ( α ) Γ ( α ) t α ( 1 α ) ( 1 β ) N ( α ) N ( β ) f ( 0 , ω 1 ( σ 1 ( 0 ) ) , ω ( σ 2 ( 0 ) ) ) + ( 1 α ) ( 1 β ) N ( α ) N ( β ) f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) + λ α N ( α ) Γ ( α ) 0 t ( t s ) α 1 u ( s ) d s + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) 0 t ( t s ) β 1 f ( s , u ( s σ 1 ( s ) ) , w ( s σ 2 ( s ) ) ) d s + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) 0 t ( t s ) α 1 f ( s , u ( s σ 1 ( s ) ) , w ( s σ 2 ( s ) ) ) d s + α β N ( α ) N ( β ) Γ ( α + β ) 0 t ( t s ) α + β 1 f ( s , u ( s σ 1 ( s ) ) , w ( s σ 2 ( s ) ) ) d s } , t J , w ( t ) = ω 2 ( 0 ) + 0 t g ( s , u ( s σ 3 ( s ) ) , w ( s σ 4 ( s ) ) ) d s , t J , u ( t ) = ω 1 ( t ) , w ( t ) = ω 2 ( t ) , t [ σ , 0 ] .
Lemma 2.
If conditions ( A 1 ) and ( A 2 ) are true, then, solving the nonlinear ML-fractional system (1) and solving the nonlinear integral system (2) are equivalent.
Remark 2.
Lemma 2 is the foundation of studying the existence and stability of solutions to system (1). Its proof is shown in Appendix A.

3. Existence of Solutions

In this section, by applying the following important fixed-point theorems, we emphasize studying the existence of solutions for system (1).
Lemma 3
(Krasnoselskii’s fixed-point theorem [51]). Let X be a Banach space, and ϕ Y X be closed convex. If the mappings P , Q : Y Y fulfill the following conditions:
( i )
P u + Q v Y , u , v Y .
( ii )
P is contraction and Q is continuous and compact.
Then, there has at least a solution u * Y such that u * = P u * + Q u * .
Lemma 4
(Schauder’s fixed point theorem [52]). Given a Banach space X , let Ω X be closed convex, and define a mapping T : Ω ¯ Ω ¯ . Then, there is a u * Ω ¯ such that T u * = u * , provided that T is completely continuous.
It is easy to see from (2) that the solution of (1) can only meet the continuity. So we define some Banach spaces Y = C ( I , R ) with the norm z C = max t I | z ( t ) | , and X = Y × Y equipped with the norm U = max { u C , w C } , where I = [ σ , T ] and U = ( u , v ) X . We shall study the existence and stability of the solution of (1) in ( X , · ) .
Under the basic assumptions ( A 1 ) and ( A 2 ) , we obtain the following two solvability results of system (1).
Theorem 1.
Assume that ( A 1 ) and ( A 2 ) hold. If conditions ( A 3 ) and ( A 4 ) are satisfied:
( A 3 )
There exist some functions m 1 ( t ) , m 2 ( t ) C ( J , R + ) such that
| f ( t , u , w ) | m 1 ( t ) , | g ( t , u , w ) | m 2 ( t ) , t J , u , w R ;
( A 4 )
0 < κ λ b α | Δ | N ( α ) Γ ( α ) < 1 ;
then system (1) has at least one solution ( u * ( t ) , w * ( t ) ) X .
Theorem 2.
Assume that ( A 1 ) and ( A 2 ) hold. If conditions ( A 5 ) and ( A 6 ) are satisfied:
( A 5 )
There exist some functions M k ( t ) , N k ( t ) C ( J , R + ) ( k = 1 , 2 , 3 ) such that, t J , u , w R ,
| f ( t , u , w ) | M 1 ( t ) + M 2 ( t ) | u | + M 3 ( t ) | w | , | g ( t , u , w ) | N 1 ( t ) + N 2 ( t ) | u | + N 3 ( t ) | w | ;
( A 6 )
0 < Θ 1 , Θ 2 < 1 , where Θ 1 = λ b α | Δ | N ( α ) Γ ( α ) + M 2 b + M 3 b | Δ | N ( α ) N ( β ) 2 ( 1 α ) ( 1 β ) + ( 1 α ) b β Γ ( β ) + ( 1 β ) b α Γ ( α ) + α β b α + β Γ ( α + β + 1 ) , Θ 2 = b N 2 b + N 3 b , · b = max t [ 0 , b ] | · | and · σ = max t [ σ , 0 ] | · | ;
then system (1) has at least one solution ( u * ( t ) , w * ( t ) ) X .
Remark 3.
We apply Lemma 3 to prove Theorem 1 in three steps. Step 1 is to define two appropriate mappings P , Q : X X as
( P U ) ( t ) = ( ( P 1 U ) ( t ) , ( P 2 U ) ( t ) ) , ( Q U ) ( t ) = ( ( Q 1 U ) ( t ) , ( Q 2 U ) ( t ) ) .
Step 2 is to prove that P is contraction. Step 3 is to show that Q is continuous and compact based on the Arzelá–Ascoli theorem. The detailed proof of Theorem 1 can be seen in Appendix B. Similarly, an appropriate mapping T : X X is defined by
( T U ) ( t ) = ( ( T 1 U ) ( t ) , ( T 2 U ) ( t ) ) , t I = [ σ , T ] , U = ( u , w ) X ,
we also employ Schauder’s fixed-point theorem (i.e., Lemma 4) to give the concrete proof of Theorem 2 in the Appendix C.
Remark 4.
Theorem 1 and Theorem 2 do not imply each other. In fact, condition ( A 3 ) is true ⇒ condition ( A 5 ) is also true, otherwise, it may not be true. However, condition ( A 6 ) is true ⇒ condition ( A 4 ) is also true, otherwise, it may not be true.

4. Stability of Ulam–Hyers type Systems

This section mainly discusses the stability of types such as Ulam–Hyers (UH) and Ulam–Hyers–Rassias (UHR) as well as generalized UH and UHR for system (1).
Lemma 5.
[52] Given a Banach space X , let ϕ E X be closed, and define a mapping T : E E . Assume that T is contractive, then, there has a unique u * E such that T u * = u * .
Consider two inequalities as follows:
| M L D 0 + β M L D 0 + α λ x ( t ) f ( t , x ( t σ 1 ( t ) ) , y ( t σ 2 ( t ) ) ) | ϵ , 0 < t b , | y ( t ) g ( t , x ( t σ 3 ( t ) ) , y ( t σ 4 ( t ) ) ) | ϵ , 0 < t b , x ( t ) = ω 1 ( t ) , y ( t ) = ω 2 ( t ) , M L D 0 + α x ( t ) = ω 3 ( t ) , t [ σ , 0 ] ,
and
| M L D 0 + β M L D 0 + α λ x ( t ) f ( t , x ( t σ 1 ( t ) ) , y ( t σ 2 ( t ) ) ) | φ ( t ) ϵ , 0 < t b , | y ( t ) g ( t , x ( t σ 3 ( t ) ) , y ( t σ 4 ( t ) ) ) | φ ( t ) ϵ , 0 < t b , x ( t ) = ω 1 ( t ) , y ( t ) = ω 2 ( t ) , M L D 0 + α x ( t ) = ω 3 ( t ) , t [ σ , 0 ] ,
where ϵ > 0 , 0 < α , β 1 , and φ : [ σ , b ] R + is continuously non-decreasing.
Similar to [18,19], we define the UH-type stability corresponding to (1) as follows.
Definition 3.
Let X = ( x , y ) and U = ( u , w ) in X be any solution of inequality (3) and the unique solution of system (1). For each ϵ > 0 ,
( 1 )
if there exists a constant C 1 > 0 such that
X ( t ) U ( t ) C 1 ϵ ,
then system (1) is called UH-stable;
( 2 )
if there exists a continuous function θ : R R + with θ ( 0 ) = 0 such that
X ( t ) U ( t ) θ ( ϵ ) ,
then system (1) is called generalized UH-stable.
Definition 4.
Let X = ( x , y ) and U = ( u , w ) in X be any solution of inequality (4) and the unique solution of system (1). For each ϵ > 0 ,
( 1 )
if there exists a constant C 2 > 0 such that
X ( t ) U ( t ) C 2 φ ( t ) ϵ , t [ σ , b ] ,
then system (1) is called UHR-stable;
( 2 )
if there exists a constant C 3 > 0 such that
X ( t ) U ( t ) C 3 φ ( t ) , t [ σ , b ] ,
then system (1) is called generalized UHR-stable.
Obviously, UH-stable ⇒ generalized UH-stable, and UHR-stable ⇒ generalized UHR-stable.
Remark 5.
A function X = ( x , y ) X solves the inequality (3) iff there has a continuous vector function ϕ = ( ϕ 1 , ϕ 2 ) in ( 0 , b ] such that
( 1 )
| ϕ 1 ( t ) | ϵ , | ϕ 2 ( t ) | ϵ , 0 < t b .
( 2 )
M L D 0 + β M L D 0 + α λ x ( t ) = f ( t , x ( t σ 1 ( t ) ) , y ( t σ 2 ( t ) ) ) + ϕ 1 ( t ) , 0 < t b .
( 3 )
y ( t ) = g ( t , x ( t σ 3 ( t ) ) , y ( t σ 4 ( t ) ) ) + ϕ 2 ( t ) , 0 < t b .
( 4 )
x ( t ) = ω 1 ( t ) , y ( t ) = ω 2 ( t ) , M L D 0 + α x ( t ) = ω 3 ( t ) , t [ σ , 0 ] .
Remark 6.
A function X = ( x , y ) X solves the inequality (4) iff there has a continuous vector function ψ = ( ψ 1 , ψ 2 ) in ( 0 , b ] such that
( 1 )
| ψ 1 ( t ) | φ ( t ) ϵ , | ψ 2 ( t ) | φ ( t ) ϵ , 0 < t b .
( 2 )
M L D 0 + β M L D 0 + α λ x ( t ) = f ( t , x ( t σ 1 ( t ) ) , y ( t σ 2 ( t ) ) ) + ψ 1 ( t ) , 0 < t b .
( 3 )
y ( t ) = g ( t , x ( t σ 3 ( t ) ) , y ( t σ 4 ( t ) ) ) + ψ 2 ( t ) , 0 < t b .
( 4 )
x ( t ) = ω 1 ( t ) , y ( t ) = ω 2 ( t ) , M L D 0 + α x ( t ) = ω 3 ( t ) , t [ σ , 0 ] .
Theorem 3.
Assume that ( A 1 ) and ( A 2 ) hold. If conditions ( A 7 ) and ( A 8 ) are satisfied:
( A 7 )
There exist a 1 ( t ) , a 2 ( t ) , b 1 ( t ) , b 2 ( t ) C ( J , R + ) such that, t J , u , u ¯ , w , w ¯ R ,
| f ( t , u , w ) | f ( t , u ¯ , w ¯ ) | a 1 ( t ) | u u ¯ | + a 2 ( t ) | w w ¯ | ,
| g ( t , u , w ) g ( t , u ¯ , w ¯ ) | b 1 ( t ) | u u ¯ | + b 2 ( t ) | w w ¯ | ;
( A 8 )
0 < Y 1 , Y 2 < 1 , where Y 1 = λ b α | Δ | N ( α ) Γ ( α ) + a 1 b + a 2 b | Δ | N ( α ) N ( β ) 2 ( 1 α ) ( 1 β ) + ( 1 α ) b β Γ ( β ) + ( 1 β ) b α Γ ( α ) + α β b α + β Γ ( α + β + 1 ) , Y 2 = b b 1 b + b 2 b and · b = max t [ 0 , b ] | · | ;
then there are the following claims:
( a )
System (1) has a unique solution U * ( t ) = ( u * ( t ) , w * ( t ) ) X .
( b )
System (1) is UH-stable and also generalized UH-stable.
( c )
System (1) is UHR-stable and also generalized UHR-stable.
Remark 7.
Similar to the proof of Theorem 2, a mapping T : X X is defined by
( T U ) ( t ) = ( ( T 1 U ) ( t ) , ( T 2 U ) ( t ) ) , t [ σ , b ] , U = ( u , w ) X ,
By the contraction mapping principle, namely, Lemma 5, we prove Theorem 3 in Appendix D. In addition, compared with Theorems 1–2, without the restriction of complete continuity of operator T , Theorem 3 ensures that the unique solution of system (1) is UH-type stable.

5. Applications

In this section, we intend to make a theoretical analysis of several examples by applying the main theorems obtained in this paper. Simultaneously, we conduct some numerical simulations by means of MATLAB.

5.1. Theoretical Analysis

Example 1.
In (1), take b = 1 , α = 0.9 , β = 0.3 , λ = 2 3 , f ( t , u , w ) = sin ( t ) [ sin ( u ) + cos ( w ) ] , g ( t , u , w ) = cos ( t ) arctan ( u + w ) , σ 1 ( t ) = 0.1 , σ 2 ( t ) = 1 + sin ( 2 t ) 6 , σ 3 ( t ) = 2 + cos ( 2 t ) 8 , σ 4 ( t ) = 0.6 , ω 1 ( t ) = 2 cos ( t ) , ω 2 ( t ) = t , ω 3 ( t ) = sin ( t ) , N ( x ) = 1 x + x Γ ( x ) , 0 < x 1 . A simple calculation gives σ = 0.6 , m 1 ( t ) = 2 | sin ( t ) | , m 2 ( t ) = π 2 | cos ( t ) | , N ( 0 ) = N ( 1 ) = 1 and
Δ = 1 λ ( 1 α ) N ( α ) 0.9292 > 0 , κ = λ α | Δ | N ( α ) Γ ( α ) 0.7125 < 1 .
Thus conditions ( A 1 ) ( A 4 ) are all fulfilled. It follows from Theorem 1 that the system of Example 1 has at least a solution U * ( t ) = ( u * ( t ) , w * ( t ) ) C ( [ 0.6 , 1 ] , R ) × C ( [ 0.6 , 1 ] , R ) .
Example 2.
In (1), take b = 1 , α = 0.4 , β = 0.5 , λ = 1 4 , f ( t , u , w ) = log ( 1 + t ) + 1 20 e t u + t 10 w , g ( t , u , w ) = t + 1 10 e t ( u + w ) , σ 1 ( t ) = 0.1 , σ 2 ( t ) = 1 + sin ( 2 t ) 6 , σ 3 ( t ) = 2 + cos ( 2 t ) 8 , σ 3 ( t ) = 0.6 , ω 1 ( t ) = 2 cos ( t ) , ω 2 ( t ) = t , ω 3 ( t ) = sin ( t ) , N ( x ) = 1 x + x Γ ( x ) , 0 < x 1 . By calculation, we get σ = 0.6 , M 1 ( t ) = | log ( 1 + t ) | , M 2 ( t ) = e t 20 , M 3 ( t ) = | t | 10 , N 1 ( t ) = | t | , N 2 ( t ) = N 3 ( t ) = 1 10 e t , N ( 0 ) = N ( 1 ) = 1 , M 2 b = e 20 , M 3 b = 1 10 , N 2 b = N 3 b = 1 10 and
Δ = 1 λ ( 1 α ) N ( α ) 0.8078 > 0 , Θ 2 = b N 2 b + N 3 b = 0.2 < 1 ,
Θ 1 = λ b α | Δ | N ( α ) Γ ( α ) + M 2 b + M 3 b | Δ | N ( α ) N ( β ) [ 2 ( 1 α ) ( 1 β ) + ( 1 α ) b β Γ ( β ) + ( 1 β ) b α Γ ( α ) + α β b α + β Γ ( α + β + 1 ) ] 0.8353 < 1 .
Thus, conditions ( A 1 ) , ( A 2 ) , ( A 5 ) , and ( A 6 ) are all true. One concludes from Theorem 2 that the system of Example 2 has at least a solution U * ( t ) = ( u * ( t ) , w * ( t ) ) C ( [ 0.6 , 1 ] , R ) × C ( [ 0.6 , 1 ] , R ) .
Example 3.
In (1), take b = 1 , α = 0.6 , β = 0.8 , λ = 1 5 , f ( t , u , w ) = cos ( t ) + t 2 20 log ( 1 + u 2 + w 2 ) , g ( t , u , w ) = sin ( t ) + 2 + sin ( 3 t ) 30 arctan ( u + w ) , σ 1 ( t ) = 0.1 , σ 2 ( t ) = 1 + sin ( 2 t ) 6 , σ 3 ( t ) = 2 + cos ( 2 t ) 8 , σ 4 ( t ) = 0.6 , ω 1 ( t ) = 2 cos ( t ) , ω 2 ( t ) = t , ω 3 ( t ) = sin ( t ) , N ( x ) = 1 x + x Γ ( x ) , 0 < x 1 . By a direct calculation, one has N ( 0 ) = N ( 1 ) = 1 , a 1 ( t ) = a 2 ( t ) = t 2 20 , b 1 ( t ) = b 2 ( t ) = 2 + sin ( 3 t ) 30 , a 1 b = a 2 b = 1 20 , b 1 b = b 2 b = 1 10 and
Δ = 1 λ ( 1 α ) N ( α ) 0.9004 > 0 , Y 2 = b b 1 b + b 2 b = 0.2 < 1 ,
Y 1 = λ b α | Δ | N ( α ) Γ ( α ) + a 1 b + a 2 b | Δ | N ( α ) N ( β ) [ 2 ( 1 α ) ( 1 β ) + ( 1 α ) b β Γ ( β ) + ( 1 β ) b α Γ ( α ) + α β b α + β Γ ( α + β + 1 ) ] 0.3455 < 1 .
Thus, we verify that conditions ( A 1 ) , ( A 2 ) , ( A 7 ) , and ( A 8 ) hold. From Theorem 3, we know that the system of Example 3 has a unique solution U * = ( u * ( t ) , w * ( t ) ) C ( [ 0.6 , 1 ] , R ) × C ( [ 0.6 , 1 ] , R ) . Meanwhile, this system is UH-, UHR-, GUH-, and GUHR-stable.

5.2. Numerical Simulation

In this subsection, we first give the numerical simulation algorithm of system (1) as follows:
Step 1: Transform the system into a system of equations. Let v ( t ) = ( M L D 0 + α λ ) u ( t ) , then system (1) becomes the following equations
M L D 0 + α u ( t ) = λ u ( t ) + v ( t ) , t ( 0 , b ] , M L D 0 + β v ( t ) = f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) , t ( 0 , b ] , w ( t ) = g ( t , u ( t σ 3 ( t ) ) , w ( t σ 4 ( t ) ) ) , t J = [ 0 , b ] , u ( t ) = ω 1 ( t ) , w ( t ) = ω 2 ( t ) , v ( t ) = ω 3 ( t ) λ ω 1 ( t ) , t [ σ , 0 ] .
Step 2: Transform Equation (5) into the integral equations. According to the proof of Lemma 2, we have
u ( t ) = ω 1 ( 0 ) + 1 α N ( α ) λ [ u ( t ) ω 1 ( 0 ) ] + [ v ( t ) ( ω 3 ( 0 ) λ ω 1 ( 0 ) ) ] + α N ( α ) Γ ( α ) 0 t ( t s ) α 1 [ λ u ( s ) + v ( s ) ] d s , t ( 0 , b ] , v ( t ) = ω 3 ( 0 ) λ ω 1 ( 0 ) + 1 β N ( β ) [ f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) f ( 0 , u ( σ 1 ( 0 ) ) , w ( σ 2 ( 0 ) ) ) ] , t ( 0 , b ] , w ( t ) = ω 2 ( 0 ) + 0 t g ( s , u ( s σ 3 ( s ) ) , w ( s σ 4 ( s ) ) ) d s , t J = [ 0 , b ] , u ( t ) = ω 1 ( t ) , w ( t ) = ω 2 ( t ) , v ( t ) = ω 3 ( t ) λ ω 1 ( t ) , t [ σ , 0 ] .
Step 3: Transform the integral equations (Equations (6)) into a differential equations of integer order. Find the first derivative at both ends of Equation (6), then Equation (6) becomes a first-order differential equation. Next, we can get a system of delayed ordinary differential equations with some simplifications and arrangements. Finally, the simulation can be carried out with the help of the DDESD toolbox in MATLAB R2019b. In addition, the error estimation and convergence order of the algorithm in the DDESD toolbox are given by Shampine [53].
To simulate and compare with the integer order equation corresponding to system (1), the integer order differential Langevin equation corresponding to system (1) is formulated by
u ( t ) λ u ( t ) = f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) , t ( 0 , b ] , w ( t ) = g ( t , u ( t σ 3 ( t ) ) , w ( t σ 4 ( t ) ) ) , t J = [ 0 , b ] , u ( t ) = ω 1 ( t ) , w ( t ) = ω 2 ( t ) , u ( t ) = ω 3 ( t ) , t [ σ , 0 ] .
Let v ( t ) = ( d d t λ ) u ( t ) , then the system of equations equivalent to Equation (7) is formed as
u ( t ) = λ u ( t ) + v ( t ) , t ( 0 , b ] , v ( t ) = f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) , t ( 0 , b ] , w ( t ) = g ( t , u ( t σ 3 ( t ) ) , w ( t σ 4 ( t ) ) ) , t J = [ 0 , b ] , u ( t ) = ω 1 ( t ) , w ( t ) = ω 2 ( t ) , v ( t ) = ω 3 ( t ) λ ω 1 ( t ) , t [ σ , 0 ] .
Based on the above preparations, we will now discuss and simulate the solutions of Examples 1–3, respectively.
(a) When all parameter values remain invariant, the simulations of solutions of Example 1 and its corresponding integer-order differential equation are shown as Figure 1 and Figure 2, respectively. u ( t ) is the solution of the Langevin equation in the figures. The comparison of numerical solutions of u ( t ) is given in Table 1. It can be seen from the simulation figures and table that the solution of the integer-order equation changes sharply, while the fractional-order equation changes gently.
(b) When all parameter values remain invariant, the simulations of solutions of Example 2 and its corresponding integer-order differential equation are shown as Figure 3 and Figure 4, respectively. u ( t ) is the solution of the Langevin equation in the figures. The comparison of numerical solutions of u ( t ) is given in Table 2. It can be seen from the simulation figures and table that the solution of the integer-order equation changes sharply, while the fractional-order equation changes gently.
(c) When all parameter values remain invariant, the simulations of solutions of Example 3 and its corresponding integer-order differential equation are shown as Figure 5 and Figure 6, respectively. The simulations of Ulam–Hyers stability of Example 3 is shown as Figure 7. u ( t ) is the solution of the Langevin equation in the figures. The comparison of numerical solutions of u ( t ) is given in Table 3 and Table 4. It can be seen from the simulation figures and tables that the solution of the integer-order equation changes sharply, while the fractional-order equation changes gently. When ϵ 0 + , the solution curves of Example 3 and inequality (3) almost coincide. This makes clear that Example 3 is UH-stable.

6. Conclusions

It is well known that the Langevin equation is a powerful tool for describing the random motion of particles in fluids. In a particularly complex viscous liquid, the integer-order Langevin equation describing the motion of particles is no longer accurate. Some scholars began to use the fractional Langevin equation as a model to study this problem and achieved good results. Most previous research on the theory and application of fractional differential systems mainly involve the Caputo- or RL-derivative. However, under certain circumstances, the Caputo- or RL-derivative is singular. To overcome this singularity, the ML-type fractional derivative is introduced. In this manuscript, we investigate the existence, uniqueness, and UH-type stability of solutions for the nonlinear ML-type fractional Langevin equation (Equation (1)). Theoretical analysis and numerical simulation of some examples verify the correctness and effectiveness of our main results. Compared with the numerical simulation of the integer-order differential equation corresponding to these examples, we find that the fractional-order differential equation is more detailed and accurate than the integer-order differential equation in describing the random motion of free particles. In practical applications, the fractional Langevin system is often affected by impulsive and random effects. Sometimes it is necessary to consider the asymptotic stability of the fractional Langevin system in the sense of Lyapunov. These are not considered in this work and need to be studied in the future.

Funding

The APC was funded by research start-up funds for high-level talents of Taizhou University.

Data Availability Statement

Not applicable.

Acknowledgments

The author would like to express his heartfelt gratitude to the editors and reviewers for their constructive comments.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A

 Proof of Lemma 2.
Let v ( t ) = M L D 0 + α λ u ( t ) , it follows from (1) that
M L D 0 + β v ( t ) = f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) , t ( 0 , b ] , M L D 0 + α u ( t ) = λ u ( t ) + v ( t ) , t ( 0 , b ] .
Adequacy.Let u , w C 1 [ 0 , b ] , if ( u ( t ) , w ( t ) ) solves system (1), then, from Lemma 1 and (A1), one has, for t J = [ 0 , b ] ,
v ( t ) = v ( 0 ) + 1 β N ( β ) f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) f ( 0 , u ( σ 1 ( 0 ) ) , w ( σ 2 ( 0 ) ) ) + β N ( β ) Γ ( β ) 0 t ( t τ ) β 1 f ( τ , u ( τ σ 1 ( τ ) ) , w ( τ σ 2 ( τ ) ) ) d τ ,
and
u ( t ) = u ( 0 ) + 1 α N ( α ) λ [ u ( t ) u ( 0 ) ] + [ v ( t ) v ( 0 ) ] + α N ( α ) Γ ( α ) 0 t ( t s ) α 1 [ λ u ( s ) + v ( s ) ] d s .
Together with (A2) and (A3), u ( 0 ) = ω 1 ( 0 ) and v ( 0 ) = M L D 0 + α u ( 0 ) λ u ( 0 ) = ω 3 ( 0 ) λ ω 1 ( 0 ) , one yields
u ( t ) = ω 1 ( 0 ) + 1 α N ( α ) [ λ [ u ( t ) ω 1 ( 0 ) ] + 1 β N ( β ) [ f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) f ( 0 , ω 1 ( σ 1 ( 0 ) ) , ω 2 ( σ 2 ( 0 ) ) ) ] + β N ( β ) Γ ( β ) 0 t ( t τ ) β 1 f ( τ , u ( τ σ 1 ( τ ) ) , w ( τ σ 2 ( τ ) ) ) d τ ] + α N ( α ) Γ ( α ) 0 t ( t s ) α 1 [ ( ω 3 ( 0 ) λ ω 1 ( 0 ) ) + λ u ( s ) + 1 β N ( β ) f ( s , u ( s σ 1 ( s ) ) , w ( t σ 2 ( s ) ) ) + β N ( β ) Γ ( β ) 0 s ( s τ ) β 1 × f ( τ , u ( τ σ 1 ( τ ) ) , w ( τ σ 2 ( τ ) ) ) d τ ] d s = 1 λ ( 1 α ) N ( α ) ω 1 ( 0 ) + λ ( 1 α ) N ( α ) u ( t ) + ( 1 α ) ( 1 β ) N ( α ) N ( β ) [ f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) f ( 0 , ω 1 ( σ 1 ( 0 ) ) , ω 2 ( σ 2 ( 0 ) ) ) ] + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) 0 t ( t τ ) β 1 f ( τ , u ( τ σ 1 ( τ ) ) , w ( τ σ 2 ( τ ) ) ) d τ + ω 3 ( 0 ) λ ω 1 ( 0 ) N ( α ) Γ ( α ) t α + λ α N ( α ) Γ ( α ) 0 t ( t s ) α 1 u ( s ) d s + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) 0 t ( t s ) α 1 f ( s , u ( s σ 1 ( s ) ) , w ( s σ 2 ( s ) ) ) d s + α β N ( α ) N ( β ) Γ ( α ) Γ ( β ) × 0 t ( t s ) α 1 0 s ( s τ ) β 1 f ( τ , u ( τ σ 1 ( τ ) ) , w ( τ σ 2 ( τ ) ) ) d τ d s .
To calculate the last quadratic integral of (A4), we exchange the order of integrals to get
0 t ( t s ) α 1 0 s ( s τ ) β 1 f ( τ , u ( τ σ 1 ( τ ) ) , w ( τ σ 2 ( τ ) ) ) d τ d s = 0 t f ( τ , u ( τ σ 1 ( τ ) ) , w ( τ σ 2 ( τ ) ) ) τ t ( t s ) α 1 ( s τ ) β 1 d s d τ = Γ ( α ) Γ ( β ) Γ ( α + β ) 0 t ( t τ ) α + β 1 f ( τ , u ( τ ) ) d τ .
Combining (A4) and (A5), one has
u ( t ) = ω 1 ( 0 ) + 1 Δ { ω 3 ( 0 ) λ ω 1 ( 0 ) N ( α ) Γ ( α ) t α ( 1 α ) ( 1 β ) N ( α ) N ( β ) f ( 0 , ω 1 ( σ 1 ( 0 ) ) , ω ( σ 2 ( 0 ) ) ) + ( 1 α ) ( 1 β ) N ( α ) N ( β ) f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) + λ α N ( α ) Γ ( α ) 0 t ( t s ) α 1 u ( s ) d s + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) 0 t ( t τ ) β 1 f ( τ , u ( τ σ 1 ( τ ) ) , w ( τ σ 2 ( τ ) ) ) d τ + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) 0 t ( t s ) α 1 f ( s , u ( s σ 1 ( s ) ) , w ( s σ 2 ( s ) ) ) d s + α β N ( α ) N ( β ) Γ ( α + β ) 0 t ( t τ ) α + β 1 f ( τ , u ( τ σ 1 ( τ ) ) , w ( τ σ 2 ( τ ) ) ) d τ } .
Obviously, Equation (A6) is the first equation of (2).
Now, integrating both sides of the second equation of (1) from 0 to t, we have
w ( t ) = w ( 0 ) + 0 t g ( s , u ( s σ 3 ( s ) ) , w ( s σ 4 ( s ) ) ) d s = ω 2 ( 0 ) + 0 t g ( s , u ( s σ 3 ( s ) ) , w ( s σ 4 ( s ) ) ) d s , t J .
When t [ σ , 0 ] , it is evident that u ( t ) = ω 1 ( t ) and w ( t ) = ω 2 ( t ) hold. Thus, one has completely derived the nonlinear integral system (2). In other words, ( u ( t ) , w ( t ) ) solves the nonlinear integral system (2).
Necessity. For t [ σ , 0 ] , let M L D 0 + α u ( t ) = ω 3 ( t ) , then u ( t ) = ω 1 ( t ) , w ( t ) = ω 2 ( t ) and M L D 0 + α u ( t ) = ω 3 ( t ) meet with system (2) ⇒ u ( t ) = ω 1 ( t ) , w ( t ) = σ 2 ( t ) and M L D 0 + α u ( t ) = ω 3 ( t ) meet with system (1). For t J , u , w C 1 [ 0 , b ] , if ( u ( t ) , w ( t ) ) solves the integral system (2), then, one deduces (A1) and the second equation of (1) by finding the one-order derivative in (A7), β -order ML-derivative in (A2), and α -order ML-derivative of (A3). Thus, for u , w C 1 [ 0 , b ] , ( u ( t ) , w ( t ) ) also solves system (1). This completes the proof. □

Appendix B

 Proof of Theorem 1.
Based on Lemma 2, for all U = ( u , w ) X , we define two mappings P , Q : X X as follows:
( P U ) ( t ) = ( ( P 1 U ) ( t ) , ( P 2 U ) ( t ) ) , ( Q U ) ( t ) = ( ( Q 1 U ) ( t ) , ( Q 2 U ) ( t ) ) ,
where
( P 1 U ) ( t ) = ω 1 ( 0 ) + 1 Δ ω 3 ( 0 ) λ ω 1 ( 0 ) N ( α ) Γ ( α ) t α + λ α N ( α ) Γ ( α ) 0 t ( t s ) α 1 u ( s ) d s , t J , ω 1 ( t ) , t [ σ , 0 ] ,
( P 2 U ) ( t ) = ω 2 ( 0 ) , t J , ω 2 ( t ) , t [ σ , 0 ] ,
( Q 1 U ) ( t ) = 1 Δ { ( 1 α ) ( 1 β ) N ( α ) N ( β ) f ( 0 , ω 1 ( σ 1 ( 0 ) ) , ω ( σ 2 ( 0 ) ) ) + ( 1 α ) ( 1 β ) N ( α ) N ( β ) × f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) 0 t ( t s ) β 1 × f ( s , u ( s σ 1 ( s ) ) , w ( s σ 2 ( s ) ) ) d s + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) 0 t ( t s ) α 1 × f ( s , u ( s σ 1 ( s ) ) , w ( s σ 2 ( s ) ) ) d s + α β N ( α ) N ( β ) Γ ( α + β ) × 0 t ( t s ) α + β 1 f ( s , u ( s σ 1 ( s ) ) , w ( s σ 2 ( s ) ) ) d s } , t J , 0 , t [ σ , 0 ] ,
and
( Q 2 U ) ( t ) = 0 t g ( s , u ( s σ 3 ( s ) ) , w ( s σ 4 ( s ) ) ) d s , t J , 0 , t [ σ , 0 ] .
It is easy to see from (A8)–(A12) that P U + Q V X , U , V X , that is, condition ( i ) in Lemma 3 holds. Next, we need to verify that condition ( ii ) in Lemma 3 also holds. In fact, t I , U = ( u 1 , w 1 ) , V = ( u 2 , w 2 ) X , when t [ 0 , b ] , we have
| ( P 1 U ) ( t ) ( P 1 V ) ( t ) | = | λ α Δ N ( α ) Γ ( α ) 0 t ( t s ) α 1 [ u 1 ( s ) u 2 ( s ) ] d s |
λ α | Δ | N ( α ) Γ ( α ) 0 t ( t s ) α 1 | u 1 ( s ) u 2 ( s ) | d s λ α | Δ | N ( α ) Γ ( α ) 0 t ( t s ) α 1 d s · u 1 u 2 C = λ t α | Δ | N ( α ) Γ ( α ) u 1 u 2 C λ b α | Δ | N ( α ) Γ ( α ) U V = κ U V ,
and
| ( P 2 U ) ( t ) ( P 3 V ) ( t ) | = | ω 2 ( 0 ) ω 2 ( 0 ) | = 0 .
When t [ σ , 0 ] , we yield
| ( P 1 U ) ( t ) ( P 1 V ) ( t ) | = | ω 1 ( t ) ω 1 ( t ) | = 0 .
and
| ( P 2 U ) ( t ) ( P 3 V ) ( t ) | = | ω 2 ( t ) ω 2 ( t ) | = 0 .
It follows from (A13)–(A16) that
P U P V λ b α | Δ | N ( α ) Γ ( α ) U V = κ U V .
In light of ( A 4 ) and (A17), one concludes that P : X X is contractive.
Next, we shall apply the Arzelá–Ascoli theorem to prove that Q : X X is completely continuous. For all t I = [ σ , b ] , U = ( u , v ) X , when t J = [ 0 , b ] , we derive from ( A 3 ) that
| ( Q 1 U ) ( t ) | 1 | Δ | { ( 1 α ) ( 1 β ) N ( α ) N ( β ) | f ( 0 , ω 1 ( σ 1 ( 0 ) ) , ω ( σ 2 ( 0 ) ) ) | + ( 1 α ) ( 1 β ) N ( α ) N ( β ) × | f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) | + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) 0 t ( t s ) β 1 × | f ( s , u ( s σ 1 ( s ) ) , w ( s σ 2 ( s ) ) ) | d s + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) 0 t ( t s ) α 1 × | f ( s , u ( s σ 1 ( s ) ) , w ( s σ 2 ( s ) ) ) | d s + α β N ( α ) N ( β ) Γ ( α + β ) × 0 t ( t s ) α + β 1 | f ( s , u ( s σ 1 ( s ) ) , w ( s σ 2 ( s ) ) ) | d s } 1 | Δ | { ( 1 α ) ( 1 β ) N ( α ) N ( β ) m 1 ( t ) + ( 1 α ) ( 1 β ) N ( α ) N ( β ) m 1 ( t ) + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) × 0 t ( t s ) β 1 m 1 ( s ) d s + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) 0 t ( t s ) α 1 m 1 ( s ) d s + α β N ( α ) N ( β ) Γ ( α + β ) 0 t ( t s ) α + β 1 m 1 ( s ) d s } max t J { m 1 ( t ) } | Δ | { 2 ( 1 α ) ( 1 β ) N ( α ) N ( β ) + ( 1 α ) t β N ( α ) N ( β ) Γ ( β ) + ( 1 β ) t α N ( α ) N ( β ) Γ ( α ) + α β t α + β N ( α ) N ( β ) Γ ( α + β + 1 ) } max t J { m 1 ( t ) } | Δ | { 2 ( 1 α ) ( 1 β ) N ( α ) N ( β ) + ( 1 α ) b β N ( α ) N ( β ) Γ ( β ) + ( 1 β ) b α N ( α ) N ( β ) Γ ( α ) + α β b α + β N ( α ) N ( β ) Γ ( α + β + 1 ) } ,
and
| ( Q 2 U ) ( t ) | 0 t | g ( s , u ( s σ 3 ( s ) ) , w ( s σ 4 ( s ) ) ) | d s b max t J { m 2 ( t ) } .
When t [ σ , 0 ] , we obtain
| ( Q 1 U ) ( t ) | = | ( Q 2 U ) ( t ) | = 0 .
From (A18)–(A20), one knows that Q : X X is uniformly bounded.
Furthermore, for all U = ( u , w ) X , t 1 , t 2 I = [ σ , b ] with t 1 < t 2 , we will argue the equicontinuity of operator Q in three cases. For the convenience of later writing, we denote
F u w ( t ) = f ( t , u ( t σ 1 ( t ) ) , w ( t σ 2 ( t ) ) ) , G u w ( t ) = g ( t , u ( t σ 3 ( t ) ) , w ( t σ 4 ( t ) ) ) .
Case 1: When 0 t 1 < t 2 b , it follows from (A21), ( A 1 ) , ( A 2 ) , and ( A 3 ) that
| ( Q 1 U ) ( t 2 ) ( Q 1 U ) ( t 1 ) | = 1 | Δ | | ( 1 α ) ( 1 β ) N ( α ) N ( β ) [ F u w ( t 2 ) F u w ( t 1 ) ] + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) 0 t 2 ( t 2 s ) β 1 F u w ( s ) d s 0 t 1 ( t 1 s ) β 1 F u w ( s ) d s + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) 0 t 2 ( t 2 s ) α 1 F u w ( s ) d s 0 t 1 ( t 1 s ) α 1 F u w ( s ) d s + α β N ( α ) N ( β ) Γ ( α + β ) 0 t 2 ( t 2 s ) α + β 1 F u w ( s ) d s 0 t 1 ( t 1 s ) α + β 1 F u w ( s ) d s | = 1 | Δ | | ( 1 α ) ( 1 β ) N ( α ) N ( β ) [ F u w ( t 2 ) F u w ( t 1 ) ] + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) × t 1 t 2 ( t 2 s ) β 1 F u w ( s ) d s + 0 t 1 ( t 2 s ) β 1 ( t 1 s ) β 1 f ( s , u ( s ) ) d s + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) [ t 1 t 2 ( t 2 s ) α 1 F u w ( s ) d s + 0 t 1 ( t 2 s ) α 1 ( t 1 s ) α 1 × F u w ( s ) d s ] + α β N ( α ) N ( β ) Γ ( α + β ) [ t 1 t 2 ( t 2 s ) α + β 1 F u w ( s ) d s + 0 t 1 ( t 2 s ) α + β 1 ( t 1 s ) α + β 1 F u w ( s ) d s ] | 1 | Δ | { ( 1 α ) ( 1 β ) N ( α ) N ( β ) | F u w ( t 2 ) F u w ( t 1 ) | + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) × t 1 t 2 ( t 2 s ) β 1 | F u w ( s ) | d s + 0 t 1 | ( t 2 s ) β 1 ( t 1 s ) β 1 | · | F u w ( s ) | d s + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) [ t 1 t 2 ( t 2 s ) α 1 | F u w ( s ) | d s + 0 t 1 | ( t 2 s ) α 1 ( t 1 s ) α 1 | × | F u w ( s ) | d s ] + α β N ( α ) N ( β ) Γ ( α + β ) [ t 1 t 2 ( t 2 s ) α + β 1 | F u w ( s ) | d s + 0 t 1 | ( t 2 s ) α + β 1 ( t 1 s ) α + β 1 | | F u w ( s ) | d s ] } 1 | Δ | { ( 1 α ) ( 1 β ) N ( α ) N ( β ) | F u w ( t 2 ) F u w ( t 1 ) | + max t J { m 1 ( t ) } ( 1 α ) β N ( α ) N ( β ) Γ ( β ) [ t 1 t 2 ( t 2 s ) β 1 d s + 0 b | ( t 2 s ) β 1 ( t 1 s ) β 1 | d s ] + max t J { m 1 ( t ) } α ( 1 β ) N ( α ) N ( β ) Γ ( α ) [ t 1 t 2 ( t 2 s ) α 1 d s + 0 b | ( t 2 s ) α 1 ( t 1 s ) α 1 | d s ] + max t J { m 1 ( t ) } α β N ( α ) N ( β ) Γ ( α + β ) [ t 1 t 2 ( t 2 s ) α + β 1 d s + 0 b | ( t 2 s ) α + β 1 ( t 1 s ) α + β 1 | d s ] } = 1 | Δ | { ( 1 α ) ( 1 β ) N ( α ) N ( β ) | F u w ( t 2 ) F u w ( t 1 ) | + max t J { m 1 ( t ) } ( 1 α ) β N ( α ) N ( β ) Γ ( β ) [ 1 β ( t 2 t 1 ) β + 0 b | ( t 2 s ) β 1 ( t 1 s ) β 1 | d s ] + max t J { m 1 ( t ) } α ( 1 β ) N ( α ) N ( β ) Γ ( α ) [ 1 α ( t 2 t 1 ) α + 0 b | ( t 2 s ) α 1 ( t 1 s ) α 1 | d s ] + max t J { m 1 ( t ) } α β N ( α ) N ( β ) Γ ( α + β ) [ 1 α + β ( t 2 t 1 ) α + β 1 + 0 b | ( t 2 s ) α + β 1 ( t 1 s ) α + β 1 | d s ] } 0 , as t 2 t 1 ,
and
| ( Q 2 U ) ( t 2 ) ( Q 2 U ) ( t 1 ) | = | 0 t 2 G u w ( s ) d s 0 t 1 G u w ( s ) d s | = | t 1 t 2 G u w ( s ) d s | t 1 t 2 | G u w ( s ) | d s max t J { m 2 ( t ) } ( t 2 t 1 ) 0 , as t 2 t 1 .
Case 2: When σ t 1 < 0 < t 2 b , then t 2 t 1 means that t 1 0 and t 2 0 + . From (A11) and (A12), we obtain
( Q U ) ( t 2 ) ( Q U ) ( t 1 ) ( Q U ) ( 0 + ) ( Q U ) ( 0 ) = ( 0 , 0 ) , as t 2 t 1 .
Case 3:When σ t 1 < t 2 0 , then
( Q U ) ( t 2 ) ( Q U ) ( t 1 ) ( 0 , 0 ) ( 0 , 0 ) = ( 0 , 0 ) ( 0 , 0 ) , as t 2 t 1 .
From (A22)–(A25), we conclude that, ϵ > 0 (small enough), ζ = ζ ( ϵ ) > 0 , for all t 1 , t 2 I , u X , when | t 2 t 1 | < ζ , there has ( Q u ) ( t 2 ) ( Q u ) ( t 1 ) < ϵ , namely, Q : X X is equicontinuous. Thereby, according to Lemmas 2 and 3, one asserts that there has at least a fixed point U * ( t ) = ( u * ( t ) , w * ( t ) ) X such that U * ( t ) = ( P U * ) ( t ) + ( Q U * ) ( t ) , which meets with system (1). This completes the proof. □

Appendix C

Proof of Throrem 2.
According to Lemma 2, a self-mapping T on X is defined as
( T U ) ( t ) = ( ( T 1 U ) ( t ) , ( T 2 U ) ( t ) ) , t I = [ σ , b ] , U = ( u , w ) X ,
( T 1 U ) ( t ) = ω 1 ( 0 ) + 1 Δ { ω 3 ( 0 ) λ ω 1 ( 0 ) N ( α ) Γ ( α ) t α ( 1 α ) ( 1 β ) N ( α ) N ( β ) F u w ( 0 ) + ( 1 α ) ( 1 β ) N ( α ) N ( β ) F u w ( t ) + λ α N ( α ) Γ ( α ) 0 t ( t s ) α 1 u ( s ) d s + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) 0 t ( t s ) β 1 F u w ( s ) d s + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) 0 t ( t s ) α 1 F u w ( s ) d s + α β N ( α ) N ( β ) Γ ( α + β ) 0 t ( t s ) α + β 1 F u w ( s ) d s } , t J , ω 1 ( t ) , t [ σ , 0 ] ,
( T 2 U ) ( t ) = ω 2 ( 0 ) + 0 t G u w ( s ) d s , t J , ω 2 ( t ) , t [ σ , 0 ] ,
where F u w ( t ) and G u w ( t ) are denoted as (A21).
Take Ω = { U = ( u , w ) X : U < R } , where R = max { ω 1 σ , ω 2 σ , Σ 1 1 Θ 1 , Σ 2 1 Θ 2 } , Σ 1 = ω 1 σ + ω 3 σ + λ ω 1 σ | Δ | N ( α ) Γ ( α ) b α + M 1 b | Δ | N ( α ) N ( β ) 2 ( 1 α ) ( 1 β ) + ( 1 α ) b β Γ ( β ) + ( 1 β ) b α Γ ( α ) + α β b α + β Γ ( α + β + 1 ) and Σ 2 = ω 2 σ + b N 1 b . Now let us prove that T ( Ω ¯ ) Ω ¯ . In fact, for all U = ( u , w ) Ω ¯ , namely, U R , when t [ σ , 0 ] , we have
| ( T 1 U ) ( t ) | = | ω 1 ( t ) | ω 1 σ R , | ( T 2 U ) ( t ) | = | ω 2 ( t ) | ω 2 σ R .
When t [ 0 , b ] , it follows from ( A 5 ) that
| F u w ( t ) | M 1 ( t ) + M 2 ( t ) | u ( t σ 1 ( t ) ) | + M 3 ( t ) | w ( t σ 2 ( t ) ) | M 1 b + R M 2 b + R M 3 b ,
and
| G u w ( t ) | N 1 ( t ) + N 2 ( t ) | u ( t σ 3 ( t ) ) | + N 3 ( t ) | w ( t σ 4 ( t ) ) | N 1 b + R N 2 b + R N 3 b .
By (A30), (A31), ( A 1 ) , ( A 2 ) , and ( A 6 ) , we derive
| ( T 1 U ) ( t ) | | ω 1 ( 0 ) | + 1 | Δ | { | ω 3 ( 0 ) | + λ | ω 1 ( 0 ) | N ( α ) Γ ( α ) b α + ( 1 α ) ( 1 β ) N ( α ) N ( β ) | F u w ( 0 ) | + ( 1 α ) ( 1 β ) N ( α ) N ( β ) | F u w ( t ) | + λ α N ( α ) Γ ( α ) 0 t ( t s ) α 1 | u ( s ) | d s + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) 0 t ( t s ) β 1 | F u w ( s ) | d s + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) 0 t ( t s ) α 1 | F u w ( s ) | d s + α β N ( α ) N ( β ) Γ ( α + β ) 0 t ( t s ) α + β 1 | F u w ( s ) | d s } ω 1 σ + 1 | Δ | { ω 3 σ + λ ω 1 σ N ( α ) Γ ( α ) b α + ( 1 α ) ( 1 β ) N ( α ) N ( β ) M 1 b + R M 2 b + R M 3 b + ( 1 α ) ( 1 β ) N ( α ) N ( β ) M 1 b + R M 2 b + R M 3 b + λ α N ( α ) Γ ( α ) 0 t ( t s ) α 1 R d s + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) 0 t ( t s ) β 1 M 1 b + R M 2 b + R M 3 b d s + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) 0 t ( t s ) α 1 M 1 b + R M 2 b + R M 3 b d s + α β N ( α ) N ( β ) Γ ( α + β ) 0 t ( t s ) α + β 1 M 1 b + R M 2 b + R M 3 b d s } ω 1 σ + 1 | Δ | { ω 3 σ + λ ω 1 σ N ( α ) Γ ( α ) b α + 2 ( 1 α ) ( 1 β ) N ( α ) N ( β ) M 1 b + R M 2 b + R M 3 b + λ b α R N ( α ) Γ ( α ) + ( 1 α ) b β N ( α ) N ( β ) Γ ( β ) M 1 b + R M 2 b + R M 3 b + ( 1 β ) b α N ( α ) N ( β ) Γ ( α ) M 1 b + R M 2 b + R M 3 b + α β b α + β N ( α ) N ( β ) Γ ( α + β + 1 ) M 1 b + R M 2 b + R M 3 b } = ω 3 σ + λ ω 1 σ | Δ | N ( α ) Γ ( α ) b α + M 1 b | Δ | N ( α ) N ( β ) [ 2 ( 1 α ) ( 1 β ) + ( 1 α ) b β Γ ( β ) + ( 1 β ) b α Γ ( α ) + α β b α + β Γ ( α + β + 1 ) ] + { λ b α | Δ | N ( α ) Γ ( α ) + M 2 b + M 3 b | Δ | N ( α ) N ( β ) [ 2 ( 1 α ) ( 1 β ) + ( 1 α ) b β Γ ( β ) + ( 1 β ) b α Γ ( α ) + α β b α + β Γ ( α + β + 1 ) ] } R = Σ 1 + Θ 1 R R ,
and
| ( T 2 U ) ( t ) | = | ω 2 ( 0 ) + 0 t G u w ( s ) d s | | ω 2 ( 0 ) | + 0 t | G u w ( s ) | d s ω 2 σ + b N 1 b + R N 2 b + R N 3 b = Σ 2 + Θ 2 R R .
From (A29), (A32) and (A33), one has ( T U ) ( t ) R , which implies that ( T U ) ( t ) Ω ¯ , and T is uniformly bounded in Ω ¯ . Similar to (A22)–(A25) in the proof of Theorem 1, it is easy to prove that T : Ω ¯ Ω ¯ is equicontinuous, so we omit it. Therefore, the complete continuity of T : Ω ¯ Ω ¯ is established by the Arzelá–Ascoli theorem. Thus, it follows from Lemmas 2 and 4 that the mapping T has a fixed point U * = ( u * , w * ) Ω ¯ , which meets with system (1). This completes the proof. □

Appendix D

Proof of Theorem 3.
Firstly, we utilize Lemma 5 to prove claim ( a ) . Define a mapping T : X X as (A26)–(A28), then, for all U 1 = ( u 1 , w 1 ) , U 2 = ( u 2 , w 2 ) X , when t J = [ σ , 0 ] , one has
| ( T 1 U 1 ) ( t ) ( T 1 U 2 ) ( t ) | = | ω 1 ( t ) ω 1 ( t ) | = 0 ,
| ( T 2 U 1 ) ( t ) ( T 2 U 2 ) ( t ) | = | ω 2 ( t ) ω 2 ( t ) | = 0 .
When t J = [ 0 , b ] , similar to (A30) and (A31), one derives from ( A 1 ) , ( A 2 ) , and ( A 7 ) that
| F u 1 w 1 ( t ) F u 2 w 2 ( t ) | a 1 b + a 2 b U 1 U 2 ,
| G u 1 w 1 ( t ) G u 2 w 2 ( t ) | b 1 b + b 2 b U 1 U 2 .
From (A36) and (A37), one gets
| ( T 1 U 1 ) ( t ) ( T 1 U 2 ) ( t ) | 1 | Δ | { ( 1 α ) ( 1 β ) N ( α ) N ( β ) | F u 1 w 1 ( 0 ) F u 2 w 2 ( 0 ) | + ( 1 α ) ( 1 β ) N ( α ) N ( β ) | F u 1 w 1 ( t ) F u 2 w 2 ( t ) | + λ α N ( α ) Γ ( α ) 0 t ( t s ) α 1 | u 1 ( s ) u 2 ( s ) | d s + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) 0 t ( t s ) β 1 | F u 1 w 1 ( s ) F u 2 w 2 ( s ) | d s + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) 0 t ( t s ) α 1 | F u 1 w 1 ( s ) F u 2 w 2 ( s ) | d s + α β N ( α ) N ( β ) Γ ( α + β ) 0 t ( t s ) α + β 1 | F u 1 w 1 ( s ) F u 2 w 2 ( s ) | d s } 1 | Δ | { ( 1 α ) ( 1 β ) N ( α ) N ( β ) a 1 b + a 2 b U 1 U 2 + ( 1 α ) ( 1 β ) N ( α ) N ( β ) a 1 b + a 2 b U 1 U 2 + λ α N ( α ) Γ ( α ) 0 t ( t s ) α 1 d s · U 1 U 2 + ( 1 α ) β N ( α ) N ( β ) Γ ( β ) 0 t ( t s ) β 1 d s · a 1 b + a 2 b U 1 U 2 + α ( 1 β ) N ( α ) N ( β ) Γ ( α ) 0 t ( t s ) α 1 d s · a 1 b + a 2 b U 1 U 2 + α β N ( α ) N ( β ) Γ ( α + β ) 0 t ( t s ) α + β 1 d s · a 1 b + a 2 b U 1 U 2 } { λ b α | Δ | N ( α ) Γ ( α ) + a 1 b + a 2 b | Δ | N ( α ) N ( β ) [ 2 ( 1 α ) ( 1 β ) + ( 1 α ) b β Γ ( β ) + ( 1 β ) b α Γ ( α ) + α β b α + β Γ ( α + β + 1 ) ] } U 1 U 2 = Y U 1 U 2 ,
and
| ( T 2 U 1 ) ( t ) ( T 2 U 2 ) ( t ) | 0 t | G u 1 w 1 ( s ) G u 2 w 2 ( s ) | d s b b 1 b + b 2 b U 1 U 2 .
Taking k = max { Y 1 , Y 2 } , we derive from (A34), (A35), (A38) and (A39) that
T U 1 T U 2 k U 1 U 2 .
By ( A 8 ) and (A40), one knows that T : X X is contractive. Thus it follows from Lemmas 2 and 5 that the operator T has a unique fixed point U * ( t ) = ( u * ( t ) , w * ( t ) ) X , which is a unique solution of system (1).
Next, we shall prove claim ( b ) . Define the operator T : X X as (A26)–(A28). From claim ( a ) , system (1) exists a unique solution U * ( t ) = ( u * ( t ) , w * ( t ) X satisfying T U * = U * . Let θ ( t ) ( 0 , 0 ) and X ( t ) = ( x ( t ) , y ( t ) ) X be a solution of inequality (3). Based upon Remark 5 and Lemma 2, we easily verify that X ( t ) = ( x ( t ) , y ( t ) ) satisfies the operator equation as below:
X ( t ) = ( T X ) ( t ) + ( T ϕ ) ( t ) ( T θ ) ( t ) .
As an analogue to (A40), one derives from (A41) and ( A 8 ) that
X U * = ( T X ) ( t ) ( T U * ) ( t ) + ( T ϕ ) ( t ) ( T θ ) ( t ) ( T X ) ( t ) ( T U * ) ( t ) + ( T ϕ ) ( t ) ( T θ ) ( t ) k X U * + k ϕ θ k X U * + k ϵ ,
which implies that
X U * k 1 k ϵ .
Thus, (A43) shows that system (1) is UH-stable and also generalized UH-stable. The proof of claim ( c ) resembles claim ( b ) , so we omit it. The proof is completed. □

References

  1. Beck, C.; Roepstorff, G. From dynamical systems to the Langevin equation. Physica A 1987, 145, 1–14. [Google Scholar] [CrossRef]
  2. Coffey, C.; Kalmykov, Y.; Waldron, J. The Langevin Equation; World Scientific: Singapore, 2004. [Google Scholar]
  3. Kubo, R. The fluctuation-dissipation theorem. Rep. Prog. Phys. 1966, 29, 255–284. [Google Scholar] [CrossRef] [Green Version]
  4. Kubo, R.; Toda, M.; Hashitsume, N. Statistical Physics II; Springer: Berlin/Heidelberg, Germany, 1991. [Google Scholar]
  5. Eab, C.; Lim, S. Fractional generalized langevin equation approach to single-file diffusion. Physica A 2010, 389, 2510–2521. [Google Scholar] [CrossRef] [Green Version]
  6. Sandev, T.; Tomovski, Z. Langevin equation for a free particle driven by power law type of noises. Phys. Lett. A 2014, 378, 1–9. [Google Scholar] [CrossRef]
  7. Jajarmi, A.; Baleanu, D.; Zarghami, V.; Pirouz, H. A new and general fractional Lagrangian approach: A capacitor microphone case study. Results Phys. 2021, 31, 104950. [Google Scholar] [CrossRef]
  8. Erturk, V.; Godwe, E.; Baleanu, D.; Kumar, P. Novel fractional-order Lagrangian to describe motion of beam on nanowire. Acta Phys. Pol. A 2021, 140, 265–272. [Google Scholar] [CrossRef]
  9. Qureshi, S.; Yusuf, A.; Aziz, S. Fractional numerical dynamics for the logistic population growth model under Conformable Caputo: A case study with real observations. Phys. Scr. 2021, 96, 114002. [Google Scholar] [CrossRef]
  10. Tanriverdi, T.; Baskonus, H.; Mahmud, A.; Muhamad, K.A. Explicit solution of fractional order atmosphere-soil-land plant carbon cycle system. Ecol. Complex. 2021, 48, 100966. [Google Scholar] [CrossRef]
  11. Abbas, M.; Ragusa, M. Solvability of Langevin equations with two Hadamard fractional derivatives via Mittag–Leffler functions. Appl. Anal. 2022, 101, 3231–3245. [Google Scholar] [CrossRef]
  12. Saeed, A.; Almalahi, M.; Abdo, M. Explicit iteration and unique solution for φ-Hilfer type fractional Langevin equations. AIMS Math. 2022, 7, 3456–3476. [Google Scholar] [CrossRef]
  13. Omaba, M.; Nwaeze, E. On a nonlinear fractional Langevin equation of two fractional orders with a multiplicative noise. Fractal Fract. 2022, 6, 290. [Google Scholar] [CrossRef]
  14. Almalahi, M.; Panchal, S.; Jarad, F. Multipoint BVP for the Langevin equation under φ-Hilfer fractional operator. J. Funct. Space 2022, 2022, 2798514. [Google Scholar] [CrossRef]
  15. Kumar, V.; Stamov, G.; Stamova, I. Controllability results for a class of piecewise nonlinear impulsive fractional dynamic systems. Appl. Math. Comput. 2022, 439, 127625. [Google Scholar] [CrossRef]
  16. Kou, Z.; Kosari, S. On a generalization of fractional Langevin equation with boundary conditions. AIMS Math. 2022, 7, 1333–1345. [Google Scholar] [CrossRef]
  17. Batiha, I.; Ouannas, A.; Albadarneh, R.; Al-Nana, A.; Momani, S. Existence and uniqueness of solutions for generalized Sturm–Liouville and Langevin equations via Caputo–Hadamard fractional-order operator. Eng. Computation 2022, 39, 2581–2603. [Google Scholar] [CrossRef]
  18. Ulam, S. A Collection of Mathematical Problems-Interscience Tracts in Pure and Applied Mathmatics; Interscience: New York, NY, USA, 1906. [Google Scholar]
  19. Hyers, D. On the stability of the linear functional equation. Proc. Nat. Acad. Sci. USA 1941, 27, 2222–2240. [Google Scholar] [CrossRef] [Green Version]
  20. Wang, J.; Li, X. A uniform method to Ulam–Hyers stability for some linear fractional equations. Mediterr. J. Math. 2016, 13, 625–635. [Google Scholar] [CrossRef]
  21. Rezaei, H.; Jung, S.; Rassias, T. Laplace transform and Hyers–Ulam stability of linear differential equations. J. Math. Anal. Appl. 2013, 403, 244–251. [Google Scholar] [CrossRef]
  22. Wang, C.; Xu, T. Hyers–Ulam stability of fractional linear differential equations involving Caputo fractional derivatives. Appl. Math. Czech. 2015, 60, 383–393. [Google Scholar] [CrossRef] [Green Version]
  23. Haq, F.; Shah, K.; Rahman, G. Hyers–Ulam stability to a class of fractional differential equations with boundary conditions. Int. J. Appl. Comput. Math. 2017, 3, 1135–1147. [Google Scholar] [CrossRef]
  24. Wang, J.; Zhou, Y. Nonlinear impulsive problems for fractional differential equations and Ulam stability. Comput. Math. Appl. 2012, 64, 3389–3405. [Google Scholar] [CrossRef] [Green Version]
  25. Wang, J.; Zhou, Y. Ulam’s type stability of impulsive ordinary differential equations. J. Math. Anal. Appl. 2012, 395, 258–264. [Google Scholar] [CrossRef] [Green Version]
  26. Huang, H.; Zhao, K.; Liu, X. On solvability of BVP for a coupled Hadamard fractional systems involving fractional derivative impulses. AIMS Math. 2022, 7, 19221–19236. [Google Scholar] [CrossRef]
  27. Wang, J.; Li, X. Eα-Ulam type stability of fractional order ordinary differential equations. J. Appl. Math. Comput. 2014, 45, 449–459. [Google Scholar] [CrossRef]
  28. Ibrahim, R. Generalized Ulam–Hyers stability for fractional differential equations. Int. J. Math. 2014, 23, 9. [Google Scholar] [CrossRef] [Green Version]
  29. Yu, X. Existence and β-Ulam–Hyers stability for a class of fractional differential equations with non-instantaneous impulses. Adv. Differ. Equ. 2015, 2015, 104. [Google Scholar] [CrossRef] [Green Version]
  30. Feckan, M.; Wang, J.; Zhou, Y. Presentation of solutions of impulsive fractional Langevin equations and existence results. Eur. Phys. J. Spec. Top. 2013, 222, 1857–1874. [Google Scholar]
  31. Wang, J.; Li, X. Ulam–Hyers stability of fractional Langevin equations. Appl. Math. Comput. 2015, 258, 72–83. [Google Scholar] [CrossRef]
  32. Gao, Z.; Yu, X. Stability of nonlocal fractional Langevin differential equations involving fractional integrals. J. Appl. Math. Comput. 2017, 53, 599–611. [Google Scholar] [CrossRef]
  33. Develi, F. Existence and Ulam–Hyers stability results for nonlinear fractional Langevin equation with modified argument. Math. Method Appl. Sci. 2022, 45, 3417–3425. [Google Scholar] [CrossRef]
  34. Baitiche, Z.; Derbazi, C.; Matar, M. Ulam stability for nonlinear-Langevin fractional differential equations involving two fractional orders in the psi-Caputo sense. Appl. Anal. 2021, 101, 4866–4881. [Google Scholar] [CrossRef]
  35. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 2015, 1, 73–85. [Google Scholar]
  36. Atangana, A.; Baleanu, D. New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model. Therm. Sci. 2016, 20, 763–769. [Google Scholar] [CrossRef] [Green Version]
  37. Sadeghi, S.; Jafari, H.; Nemati, S. Operational matrix for Atangana–Baleanu derivative based on Genocchi polynomials for solving FDEs. Chaos Soliton Fract. 2020, 135, 109736. [Google Scholar] [CrossRef]
  38. Ganji, R.; Jafari, H.; Baleanu, D. A new approach for solving multi variable orders differential equations with Mittag–Leffler kernel. Chaos Soliton Fract. 2020, 130, 109405. [Google Scholar] [CrossRef]
  39. Tajadodi, H. A numerical approach of fractional advection-diffusion equation with Atangana-Baleanu derivative. Chaos Soliton Fract. 2020, 130, 109527. [Google Scholar] [CrossRef]
  40. Khan, H.; Khan, A.; Jarad, F.; Shah, A. Existence and data dependence theorems for solutions of an ABC-fractional order impulsive system. Chaos Soliton Fract. 2020, 131, 109477. [Google Scholar] [CrossRef]
  41. Khan, H.; Li, Y.; Khan, A. Existence of solution for a fractional-order Lotka–Volterra reaction-diffusion model with Mittag–Leffler kernel. Math. Meth. Appl. Sci. 2019, 42, 3377–3387. [Google Scholar] [CrossRef]
  42. Acay, B.; Bas, E.; Abdeljawad, T. Fractional economic models based on market equilibrium in the frame of different type kernels. Chaos Soliton Fract. 2020, 130, 109438. [Google Scholar] [CrossRef]
  43. Khan, A.; Gómez-Aguilar, J.; Khan, T.; Khan, H. Stability analysis and numerical solutions of fractional order HIV/AIDS model. Chaos Soliton Fract. 2019, 122, 119–128. [Google Scholar] [CrossRef]
  44. Khan, H.; Gómez-Aguilar, J.; Alkhazzan, A.; Khan, A. A fractional order HIV-TB coinfection model with nonsingular Mittag–Leffler Law. Math. Meth. Appl. Sci. 2020, 43, 3786–3806. [Google Scholar] [CrossRef]
  45. Koca, I. Modelling the spread of Ebola virus with Atangana–Baleanu fractional operators. Eur. Phys. J. Plus 2018, 133, 100. [Google Scholar] [CrossRef]
  46. Morales-Delgado, V.; Gómez-Aguilar, J.; Taneco-Hernández, M.; Escobar-Jiménez, R.F.; Olivares-Peregrino, V.H. Mathematical modeling of the smoking dynamics using fractional differential equations with local and nonlocal kernel. J. Nonlinear Sci. Appl. 2018, 11, 1004–1014. [Google Scholar] [CrossRef] [Green Version]
  47. Sene, N.; Gautam, S. Generalized Mittag–Leffler input stability of the fractional differential equations. Symmetry 2019, 11, 608. [Google Scholar] [CrossRef] [Green Version]
  48. Zhao, K. Stability of a nonlinear ML-nonsingular kernel fractional Langevin system with distributed lags and integral control. Axioms 2022, 11, 350. [Google Scholar] [CrossRef]
  49. Zhao, K. Existence, stability and simulation of a class of nonlinear fractional Langevin equations involving nonsingular Mittag–Leffler kernel. Fractal Fract. 2022, 6, 469. [Google Scholar] [CrossRef]
  50. Jarad, F.; Abdeljawad, T.; Hammouch, Z. On a class of ordinary differential equations in the frame of Atangana–Baleanu fractional derivative. Chaos Soliton Fract. 2018, 117, 16–20. [Google Scholar] [CrossRef]
  51. Granas, A.; Dugundji, J. Fixed Point Theory; Springer: New York, NY, USA, 2003. [Google Scholar]
  52. Guo, D.; Lakshmikantham, V. Nonlinear Problems in Abstract Cone; Academic Press: Orlando, FL, USA, 1988. [Google Scholar]
  53. Shampine, L. Solving ODEs and DDEs with residual control. Appl. Numer. Math. 2005, 52, 113–127. [Google Scholar] [CrossRef]
Figure 1. Simulation of solutions in Example 1.
Figure 1. Simulation of solutions in Example 1.
Fractalfract 06 00725 g001
Figure 2. Simulation of solutions of the integer-order equation corresponding to Example 1.
Figure 2. Simulation of solutions of the integer-order equation corresponding to Example 1.
Fractalfract 06 00725 g002
Figure 3. Simulation of solutions in Example 2.
Figure 3. Simulation of solutions in Example 2.
Fractalfract 06 00725 g003
Figure 4. Simulation of solutions of the integer-order equation corresponding to Example 2.
Figure 4. Simulation of solutions of the integer-order equation corresponding to Example 2.
Fractalfract 06 00725 g004
Figure 5. Simulation of solutions in Example 3.
Figure 5. Simulation of solutions in Example 3.
Fractalfract 06 00725 g005
Figure 6. Simulation of solutions of the integer-order equation corresponding to Example 3.
Figure 6. Simulation of solutions of the integer-order equation corresponding to Example 3.
Fractalfract 06 00725 g006
Figure 7. Simulation of UH-stability in Example 3 with ϵ = 0, 0.005.
Figure 7. Simulation of UH-stability in Example 3 with ϵ = 0, 0.005.
Fractalfract 06 00725 g007
Table 1. Comparison of numerical solutions of u ( t ) in Example 1.
Table 1. Comparison of numerical solutions of u ( t ) in Example 1.
t00.20.40.60.81
Fractional order2.00001.99711.98581.95821.90331.8187
Integer order2.00001.81011.61931.44441.30581.2184
Table 2. Comparison of numerical solutions of u ( t ) in Example 2.
Table 2. Comparison of numerical solutions of u ( t ) in Example 2.
t00.20.40.60.81
Fractional order2.00001.47080.86500.44820.21470.1024
Integer order2.00001.81001.61901.43551.26861.1276
Table 3. Comparison of numerical solutions of u ( t ) in Example 3.
Table 3. Comparison of numerical solutions of u ( t ) in Example 3.
t00.20.40.60.81
Fractional order2.00001.66141.25760.93850.70520.5470
Integer order2.00001.82981.69241.58711.51151.4627
Table 4. Comparison of numerical solutions for the UH-stability of u ( t ) in Example 3.
Table 4. Comparison of numerical solutions for the UH-stability of u ( t ) in Example 3.
t00.20.40.60.81
ε = 0 2.00001.66141.25760.93850.70520.5470
ε = 0.005 2.00001.66111.25690.93760.70420.5459
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, K. Stability of a Nonlinear Langevin System of ML-Type Fractional Derivative Affected by Time-Varying Delays and Differential Feedback Control. Fractal Fract. 2022, 6, 725. https://doi.org/10.3390/fractalfract6120725

AMA Style

Zhao K. Stability of a Nonlinear Langevin System of ML-Type Fractional Derivative Affected by Time-Varying Delays and Differential Feedback Control. Fractal and Fractional. 2022; 6(12):725. https://doi.org/10.3390/fractalfract6120725

Chicago/Turabian Style

Zhao, Kaihong. 2022. "Stability of a Nonlinear Langevin System of ML-Type Fractional Derivative Affected by Time-Varying Delays and Differential Feedback Control" Fractal and Fractional 6, no. 12: 725. https://doi.org/10.3390/fractalfract6120725

APA Style

Zhao, K. (2022). Stability of a Nonlinear Langevin System of ML-Type Fractional Derivative Affected by Time-Varying Delays and Differential Feedback Control. Fractal and Fractional, 6(12), 725. https://doi.org/10.3390/fractalfract6120725

Article Metrics

Back to TopTop