Next Article in Journal
Maximum Principles for Fractional Differential Inequalities with Prabhakar Derivative and Their Applications
Next Article in Special Issue
The Effect of a Nonlocal Thermoelastic Model on a Thermoelastic Material under Fractional Time Derivatives
Previous Article in Journal
Systems of Riemann–Liouville Fractional Differential Equations with ρ-Laplacian Operators and Nonlocal Coupled Boundary Conditions
Previous Article in Special Issue
A Stochastic Bayesian Regularization Approach for the Fractional Food Chain Supply System with Allee Effects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Controllability and Hyers–Ulam Stability of Fractional Systems with Pure Delay

1
Department of Mathematical Sciences, College of Sciences, Princess Nourah Bint Abdulrahman University, Riyadh 11671, Saudi Arabia
2
School of Mathematics, Harbin Institute of Technology, Harbin 150001, China
3
Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(10), 611; https://doi.org/10.3390/fractalfract6100611
Submission received: 3 September 2022 / Revised: 22 September 2022 / Accepted: 12 October 2022 / Published: 20 October 2022

Abstract

:
Linear and nonlinear fractional-delay systems are studied. As an application, we derive the controllability and Hyers–Ulam stability results using the representation of solutions of these systems with the help of their delayed Mittag–Leffler matrix functions. We provide some sufficient and necessary conditions for the controllability of linear fractional-delay systems by introducing a fractional delay Gramian matrix. Furthermore, we establish some sufficient conditions of controllability and Hyers–Ulam stability of nonlinear fractional-delay systems by applying Krasnoselskii’s fixed-point theorem. Our results improve, extend, and complement some existing ones. Finally, numerical examples of linear and nonlinear fractional-delay systems are presented to demonstrate the theoretical results.

1. Introduction

The fractional delay differential equations and their applications have gained significant attention owing to their successful modeling in several fields of science and engineering, such as disease, control theory, signal analysis, diffusion processes, biology, forced oscillations, population dynamics, viscoelastic systems, computer engineering, and finance; see, for instance, [1,2,3,4,5,6,7,8]. Recently, the representation of solutions of time-delay systems has been considered. In particular, the pioneering study [9,10] produced several innovative findings on the representations of solutions of time-delay systems, which were used in the control problems and stability analysis; see, for instance, [11,12,13,14,15,16,17,18,19,20,21] and the references therein.
On the one hand, the controllability of systems is one of the most fundamental and significant concepts in modern control theory, which consists of determining the control parameters that steer the solutions of a control system from its initial state to its final state using a set of admissible controls, where initial and final states may vary over an entire space. In recent decades, there has been considerable interest in the controllability analysis of fractional-delay systems of order α 0 , 1 and α 1 , 2 , and several methods for studying the controllability results have been developed, for example, the robust and universal methods [22]; the Laplace transform technique, the Mittag–Leffler function and fixed-point argument [23]; Martelli’s fixed-point theorem, multivalued functions, and cosine and sine families [24]; the Mittag–Leffler matrix functions and the Schauder fixed-point theorem [20,25,26]; the Mittag–Leffler matrix function, the Gramian matrix, and the iterative technique [27]; the solution operator theory, fractional calculations, and fixed point techniques [28]; and the delayed fractional Gram matrix and the explicit solution formula [29]. On the other hand, the Hyers–Ulam stability of fractional delay systems has been studied recently by many authors; see, for example, [19,30,31] and the references therein.
However, to the best of our knowledge, no research has been conducted on the controllability of linear fractional-delay systems of the form
C D 0 + α y x + A y x h = B u ( x ) , x Ω : = 0 , x 1 , y x ψ x , y x ψ x , h x 0 ,
and the controllability and Hyers–Ulam stability of the corresponding nonlinear fractional-delay systems of the form
C D 0 + α y x + A y x h = f x , y x + B u ( x ) , x Ω , y x ψ x , y x ψ x , h x 0 ,
where C D 0 + α is called the Caputo fractional derivative of order α 1 , 2 with the lower index zero, h > 0 is a delay, x 1 > n 1 h , y x R n , ψ C h , 0 , R n , A R n × n and B R n × m are any matrices, f C Ω × R n , R n is a given function, and u x R m shows control vector.
Elshenhab and Wang [11] have presented a novel formulation of solutions to the linear fractional-delay systems
C D 0 + α y x + A y x h = f x , x 0 , y x ψ x , y x ψ x , h x 0 ,
of the following form:
y x = H h , α A x h α ψ 0 + M h , α A x h α ψ 0 A h 0 S h , α A x 2 h α ψ d + 0 x S h , α A x h α f d ,
where H h , α A x α , M h , α A x α , and S h , α A x α are known as the delayed Mittag–Lefflertype matrix functions formulated by
H h , α A x α : = Θ , < x < h , I , h x < 0 , I A x α Γ 1 + α , 0 x < h , I A x α Γ 1 + α + A 2 x h 2 α Γ 1 + 2 α + + ( 1 ) r A r x r 1 h r α Γ 1 + r α , r 1 h x < r h ,
M h , α A x α : = Θ , < x < h , I x + h , h x < 0 , I x + h A x α + 1 Γ 2 + α , 0 x < h , I x + h A x α + 1 Γ 2 + α + A 2 x h 2 α + 1 Γ 2 + 2 α + + ( 1 ) r A r x r 1 h r α + 1 Γ 2 + r α , r 1 h x < r h ,
and
S h , α A x α : = Θ , < x < h , I x + h α 1 Γ α , h x < 0 , I x + h α 1 Γ α A x 2 α 1 Γ 2 α , 0 x < h , I x + h α 1 Γ α A x 2 α 1 Γ 2 α + A 2 x h 3 α 1 Γ 3 α + + ( 1 ) r A r x r 1 h α r + 1 1 Γ α r + 1 , r 1 h x < r h ,
respectively, where the notation Θ and I are the n × n null and identity matrix, respectively, Γ is a gamma function, and r = 0 , 1 , 2 , . . . .
Applying Formula (4), the solution of (2) can be represented as
y x = H h , α A x h α ψ 0 + M h , α A x h α ψ 0 A h 0 S h , α A x 2 h α ψ d + 0 x S h , α A x h α f , y d + 0 x S h , α A x h α B u d ,
Motivated by [11,16], the explicit solutions Formula (8) of (3) combined with the delayed Mittag–Leffler matrix functions are employed as an application to derive controllability results on Ω = 0 , x 1 .
The rest of this paper is structured as follows: in Section 2, we present some preliminaries, some basic notation and definitions, and some useful lemmas. In Section 3, we derive sufficient and necessary conditions for the controllability of (1) by introducing a fractional delay Gramian matrix. In Section 4, we establish sufficient conditions of the controllability of (2) by applying Krasnoselskii’s fixed-point theorem. In Section 5, we discuss the Hyers–Ulam stability of (2) on the finite-time interval 0 , x 1 . Finally, we provide numerical examples of linear and nonlinear fractional-delay systems to demonstrate the theoretical results.

2. Preliminaries

Throughout the paper, we refer to C Ω , R n as the Banach space of a vector-valued continuous function from Ω R n endowed with the norm y C Ω = max x Ω y x for a norm · on R n , and the matrix norm as A = max y = 1 A y , where A : R n R n . We define a space C 1 Ω , R n = y C Ω , R n : y C Ω , R n . Let X, Y be two Banach spaces and L b X , Y be the space of bounded linear operators from X to Y. Now, L p Ω , Y indicates the Banach space of functions f : Ω Y that are Bochner integrable normed by f L p Ω , Y for some 1 < p < . Furthermore, we let ψ C = max s h , 0 ψ s and ψ C = max s h , 0 ψ s .
We mention some basic concepts and lemmas utilized throughout this paper.
Definition 1.
([5]). The Mittag–Leffler function with two parameters is given by
E σ , τ x = r = 0 x r Γ σ r + τ , σ , τ > 0 , x C .
In the case of τ = 1 , then
E σ , 1 x = E σ x = r = 0 x r Γ σ r + 1 , σ > 0 .
Definition 2.
([5]). A function y : h , R n has the Caputo fractional derivative of order α 1 , 2 with a lower index 0 given by
C D 0 + α y x = 1 Γ 2 α 0 x y x α 1 d , x > 0 .
Definition 3.
([32]). The systems (1) or (2) are controllable on Ω = 0 , x 1 if there is a control function u L 2 Ω , R m such that (1) or (2) has a solution y : h , x 1 R n with y 0 = y 0 , y 0 = y 0 satisfies y x 1 = y 1 for all y 0 , y 0 , y 1 R n .
Definition 4.
([33]). The system (2) is Hyers–Ulam stable on [ 0 , x 1 ] if there is, for a given constant ε > 0 , a function φ C Ω , R n satisfying the inequality
C D 0 + α φ x + A φ x h f x , φ x B u ( x ) ε , x 0 , x 1 ,
there exists a solution y C Ω , R n of (2) and a constant M > 0 such that
φ x y x M ε , for all x 0 , x 1 .
Remark 1.
([33]). A function φ C Ω , R n is a solution of the inequality (9) if and only if there is a function π C Ω , R n such that
(i)
π x ε , x Ω .
(ii)
C D 0 + α φ x = A φ x h + f x , φ x + B u ( x ) + π x , x Ω .
Lemma 1.
([17]). The following inequalities hold:
H h , α A x α E α A x α ,
M h , α A x α x + h E α , 2 A x + h α ,
and
S h , α A x α x + h α 1 E α , α A x + h α .
for any x r 1 h , r h , r = 1 , 2 , .
Lemma 2.
Let α > 0 and φ C Ω , R n be a solution of the inequality (9). Then there exists, for a given constant ε > 0 , a solution φ satisfying the inequality
φ x φ x x α ε α E α , α A x α .
where
φ x = H h , α A x h α ψ 0 + M h , α A x h α ψ 0 A h 0 S h , α A x 2 h α ψ d + 0 x S h , α A x h α f , φ d + 0 x S h , α A x h α B u φ d .
Proof. 
From Remark 1, the solution of the equation
C D 0 + α φ x = A φ x h + f x , φ x + B u ( x ) + π x , x Ω ,
can be written as
φ x = H h , α A x h α ψ 0 + M h , α A x h α ψ 0 A h 0 S h , α A x 2 h α ψ d + 0 x S h , α A x h α f , φ d + 0 x S h , α A x h α B u φ d + 0 x S h , α A x h α π d .
From Lemma 1, we obtain
φ x φ x 0 x S h , α A x h α π d ε 0 x x α 1 E α , α A x α d x α ε α E α , α A x α ,
for all x Ω . This ends the proof. □
Lemma 3.
( Krasnoselskii s fixed - point theorem , [34]). Let C be a closed, convex, and non-empty subset of a Banach space X. Suppose that the operators A and B are maps from C into X such that A x + B y C for every pair x, y C . If A is compact and continuous, B is a contraction mapping. Then, there exists z C such that z = A z + B z .

3. Controllability of Linear Fractional Delay System

In this section, we establish some sufficient and necessary conditions of controllability of (1) by introducing a fractional delay Gramian matrix defined by
W h , α 0 , x 1 = 0 x 1 S h , α A x 1 h α B B T S h , α A T x 1 h α d .
It follows from the definition of the matrix W h , α 0 , x 1 that it is always positive semidefinite for x 1 0 .
Theorem 1.
The linear system (1) is controllable if and only if W h , α 0 , x 1 is positive definite.
Proof. 
Sufficiency. Let W h , α 0 , x 1 be positive definite; then, it will be non-singular and its inverse will be well-defined. As a result, we can derive the associated control input u x , for any finite terminal conditions y 1 , y 1 R n , as
u x = B T S h , α A T x 1 h x α W h , α 1 0 , x 1 β ,
where
β = y 1 H h , α A x 1 h α ψ 0 M h , α A x 1 h α ψ 0 + A h 0 S h , α A x 1 2 h α ψ d .
From (8), the solution y x 1 of (1) can be formulated as
y x 1 = H h , α A x 1 h α ψ 0 + M h , α A x 1 h α ψ 0 A h 0 S h , α A x 1 2 h α ψ d + 0 x 1 S h , α A x 1 h α B u d .
Substituting (11) into (13), we obtain
y x 1 = H h , α A x 1 h α ψ 0 + M h , α A x 1 h α ψ 0 A h 0 S h , α A x 1 2 h α ψ d + 0 x 1 S h , α A x 1 h α B B T S h , α A T x 1 h α d W h , α 1 0 , x 1 β .
From (10), (12), and (14), we obtain
y x 1 = H h , α A x 1 h α ψ 0 + M h , α A x 1 h α ψ 0 A h 0 S h , α A x 1 2 h α ψ d + β = y 1 .
We can see from (3) and (4) that the boundary conditions hold. Thus, (1) is controllable.
Necessity. Assume that (1) is controllable. For the sake of a contradiction, suppose that W h , α 0 , x 1 is not positive definite, and there exists at least a nonzero vector z R n such that z T W h , α 0 , x 1 z = 0 , which implies that
0 = z T W h , α 0 , x 1 z = 0 x 1 z T S h , α A x 1 h α B B T S h , α A T x 1 h α z d = 0 x 1 z T S h , α A x 1 h α B z T S h , α A x 1 h α B T d = 0 x 1 z T S h , α A x 1 h α B z T S h , α A x 1 h α B T d = 0 x 1 z T S h , α A x 1 h α B d .
Hence
z T S h , α A x 1 h α B = 0 , , 0 : = 0 T , for all Ω ,
where 0 denotes the n dimensional zero vector. Consider the initial points y 0 = y 0 = 0 and the final point y 1 = z at x = x 1 . Since (1) is controllable, from Definition 3, there exists a control function u 1 x that steers the response from 0 to y 1 = z at x = x 1 . Then,
y 1 = z = A h 0 S h , α A x 1 2 h α ψ d + 0 x 1 S h , α A x 1 h α B u 1 d .
Multiplying (16) by z T and using (15), we obtain z T z = 0 . This is a contradiction to z 0 . Thus, W h , α 0 , x 1 is positive definite. This ends the proof. □
Remark 2.
We note in the case of α = 2 in (1) that Theorem 1 coincides with the conclusion of Corollary 1 in [16].
Remark 3.
Under condition A, a nonsingular n × n matrix, we note in the case of α = 2 , A = A 2 in (1) that Theorem 1 coincides with the conclusion of Theorem 3.1 in [21] and Corollary 2 in [16].

4. Controllability of Nonlinear Fractional Delay System

In this section, we estabilish sufficient conditions of controllability of (2) using Krasnoselskii’s fixed point theorem.
We impose the following assumptions:
(G1)
The function f : Ω × R n R n is continuous, and there exists a constant L f L q Ω , R + and q > 1 such that
f x , y 1 f x , y 2 L f x y 1 y 2 , for all x Ω , y 1 , y 2 R n ,
let sup x Ω f x , 0 = M f < .
(G2)
The linear operator Y : L 2 Ω , R m R n defined by
Y = 0 x 1 S h , α A x 1 h α B u d .
Suppose that Y 1 exists and takes values in L 2 Ω , R m / ker Y , and there exists a constant M 1 > 0 such that Y 1 M 1 .
To establish our result, we now employ Krasnoselskii’s fixed point theorem.
Theorem 2.
Let ( G 1 ) and ( G 2 ) hold. Then, the nonlinear system (2) is controllable if
M 2 1 + M 1 x 1 α α E α , α A x 1 α B < 1 ,
where
M 2 = x 1 α 1 q α p p + 1 1 p E α , α A x 1 α L f L q Ω , R + and 1 p + 1 q = 1 , p , q > 1 .
Proof. 
Before we start to prove this theorem, we shall use the following assumptions and estimates: We consider the set
B ϵ = y C h , x 1 , R n : y C h , x 1 = sup x h , x 1 y x ϵ .
Let x 0 , x 1 . From ( G 1 ) and Hölder inequality, we obtain
0 x x α 1 E α , α A x α L f d 0 x x α 1 E α , α A x α p d 1 p 0 x L f q d 1 q E α , α A x α 0 x x α 1 p d 1 p 0 x L f q d 1 q = x α 1 q α p p + 1 1 p E α , α A x α L f L q Ω , R + .
Furthermore, consider the following control function u y :
u y x = Y 1 y 1 H h , α A x 1 h α ψ 0 M h , α A x 1 h α ψ 0 + A h 0 S h , α A x 1 2 h α ψ d 0 x 1 S h , α A x 1 h α f , y d x ,
for x Ω . From (18), (19), ( G 1 ) , ( G 2 ) , and Lemma 1, we obtain
u y x Y 1 y 1 + H h , α A x 1 h α ψ 0 + M h , α A x 1 h α ψ 0 + A h 0 S h , α A x 1 2 h α ψ d + 0 x 1 S h , α A x 1 h α f , y d ) M 1 y 1 + M 1 E α A x 1 h α ψ C + M 1 x 1 E α , 2 A x 1 α ψ C + M 1 A ψ C h 0 x 1 h α 1 E α , α A x 1 h α d + M 1 0 x 1 x 1 α 1 E α , α A x 1 α L f y d + M 1 0 x 1 x 1 α 1 E α , α A x 1 α f , 0 d M 1 y 1 + M 1 E α A x 1 h α ψ C + M 1 x 1 E α , 2 A x 1 α ψ C + M 1 A ψ C x 1 α α E α , α A x 1 α + M 1 x α 1 q α p p + 1 1 p E α , α A x α L f L q Ω , R + y C Ω + M 1 M f x 1 α α E α , α A x 1 α M 1 y 1 + M 1 M 2 ϵ + M 1 θ x 1 ,
where
θ x = E α A x h α ψ C + x E α , 2 A x α ψ C + x α A ψ C + M f α E α , α A x α .
Furthermore,
u y x u z x M 1 0 x 1 S h , α A x 1 h α f , y f , z d M 1 0 x 1 S h , α A x 1 h α L f y z d M 1 M 2 y z C Ω .
We also define the operators L 1 , L 2 on B ϵ as follows:
L 1 y x = H h , α A x h α ψ 0 + M h , α A x h α ψ 0 A h 0 S h , α A x 2 h α ψ d + 0 x S h , α A x h α B u y d ,
L 2 y x = 0 x S h , α A x h α f , y d .
Now, we see that B ϵ is a closed, bounded, and convex set of C h , x 1 , R n . Therefore, our proof is divided into three main steps.
Step 1. We prove L 1 y + L 2 z B ϵ for all y, z B ϵ .
For each x Ω and y, z B ϵ , using (20), we obtain
L 1 y + L 2 z C h , x 1 = sup x h , x 1 L 1 y + L 2 z x sup x h , x 1 H h , α A x h α ψ 0 + M h , α A x h α ψ 0 + A h 0 S h , α A x 2 h α ψ d + 0 x S h , α A x h α B u y d + 0 x S h , α A x h α f , z d E α A x h α ψ C + x E α , 2 A x α ψ C + x α A ψ C α E α , α A x α + M f x α α E α , α A x α + x α α E α , α A x α B M 1 y 1 + M 1 M 2 ϵ + M 1 θ x 1 d + x α 1 q α p p + 1 1 p E α , α A x α L f L q Ω , R + z C Ω θ x 1 + M 2 ϵ + M 1 x α α E α , α A x 1 α B y 1 + M 1 M 2 ϵ x α α E α , α A x 1 α B + M 1 θ x 1 x α α E α , α A x 1 α B θ x 1 1 + M 1 x 1 α α E α , α A x 1 α B + M 1 x 1 α α E α , α A x 1 α B y 1 + M 2 1 + M 1 x 1 α α E α , α A x 1 α B ϵ .
Thus, for some ϵ sufficiency large, and from (17), we have L 1 y + L 2 z B ϵ .
Step 2. We prove L 1 : B ϵ C h , x 1 , R n is a contraction.
For each x Ω and y, z B ϵ , using (21), we obtain
L 1 y x L 1 z x 0 x S h , α A x h α B u y u z d B M 1 M 2 y z C Ω 0 x S h , α A x h α d x 1 α B M 1 M 2 α E α , α A x 1 α y z C Ω μ y z C Ω ,
where μ : = x 1 α B M 1 M 2 α E α , α A x 1 α . From (17), note μ < 1 ; we conclude that L 1 is a contraction mapping.
Step 3. We prove L 2 : B ϵ C h , x 1 , R n is a continuous compact operator.
Firstly, we show that L 2 is continuous. Let y n be a sequence such that y n y as n in B ϵ . Thus, for each x Ω , using (23) and Lebesgue’s dominated convergence theorem, we obtain
L 2 y n x L 2 y x 0 x S h , α A x h α f , y n f , y d 0 x x α 1 E α , α A x α L f y n y d 0 , as n .
Hence L 2 : B ϵ C h , x 1 , R n is a continuous.
Next, we prove that L 2 is uniformly bounded on B ϵ . For each x Ω , y B ϵ , we have
L 2 y = sup x Ω L 2 y x sup x Ω 0 x S h , α A x h α f , y d x α 1 q α p p + 1 1 p E α , α A x α L f L q Ω , R + y C Ω + M f x α α E α , α A x α M 2 ϵ + M f x 1 α α E α , α A x 1 α ,
which implies that L 2 is uniformly bounded on B ϵ .
It remains to be shown that L 2 is equicontinuous. For each x 2 , x 3 Ω , 0 < x 2 < x 3 x 1 and y B ϵ , using (23), we obtain
L 2 y x 3 L 2 y x 2 0 x 3 S h , α A x 3 h α f , y d 0 x 2 S h , α A x 2 h α f , y d = Ψ 1 + Ψ 2 ,
where
Ψ 1 = x 2 x 3 S h , α A x 3 h α f , y d ,
and
Ψ 2 = 0 x 2 S h , α A x 3 h α S h , α A x 2 h α f , y d .
Thus
L 2 y x 3 L 2 y x 2 Ψ 1 + Ψ 2 .
Now, we can check Ψ i 0 as x 2 x 3 , i = 1 , 2. For Ψ 1 , we obtain
Ψ 1 x 2 x 3 x 3 α 1 E α , α A x 3 α L f y d + x 2 x 3 x 3 α 1 E α , α A x 3 α f , 0 d x 3 x 2 α 1 q α p p + 1 1 p E α , α A x 3 α L f L q Ω , R + y C Ω + M f x 3 x 2 α α E α , α A x 3 α 0 , as x 2 x 3 .
For Ψ 2 , we obtain
Ψ 2 ϵ 0 x 2 S h , α A x 3 h α S h , α A x 2 h α L f d + M f 0 x 2 S h , α A x 3 h α S h , α A x 2 h α d .
From (7), we know that the delayed Mittag–Leffler type matrix function S h , α A x α is uniformly continuous for x Ω . Thus,
S h , α A x 3 h α S h , α A x 2 h α 0 , as x 2 x 3 .
Therefore, we have Ψ i 0 as x 2 x 3 , i = 1 , 2, which implies that, using (24),
L 2 y x 3 L 2 y x 2 0 , as x 2 x 3 ,
for all y B ϵ . Thus, the Arzelà–Ascoli theorem tells us that L 2 is compact on B ϵ .
Therefore, according to Krasnoselskii’s fixed-point theorem (Lemma 3), L 1 + L 2 has a fixed point y on B ϵ . In addition, y is also a solution of (2) and L 1 y + L 2 y x 1 = y 1 . This means that u y steers the system (2) from y 0 to y 1 in finite time x 1 , which implies that (2) is controllable on Ω . This completes the proof. □
Remark 4.
We note in the case of α = 2 in (2) that Theorem 2 coincides with the conclusion of Corollary 3 in [16].
Remark 5.
Under condition A, there is a nonsingular n × n matrix; we note in the case of α = 2 and A = A 2 in (2) that Theorem 2 coincides with the conclusion of Theorem 4.1 in [21] and Corollary 4 in [16].

5. Hyers–Ulam Stability of Nonlinear Fractional Delay System

In this section, we discuss the Hyers–Ulam stability of (2) on the finite time interval 0 , x 1 .
Theorem 3.
Let ( G 1 ) , ( G 2 ) , and (17) be satisfied. Then, the system (2) is Hyers–Ulam stable.
Proof. 
With the help of Theorem 2, let ϖ C Ω , R n be a solution of the inequality (9) and y be the unique solution of (2), that is,
y x = H h , α A x h α ψ 0 + M h , α A x h α ψ 0 A h 0 S h , α A x 2 h α ψ d + 0 x S h , α A x h α f , y d + 0 x S h , α A x h α B u y d .
From Lemma 2, and by a similar way in the proof of Theorem 2 and by virtue of (21), we obtain
ϖ x y x ϖ x ϖ x + ϖ x y x x α ε α E α , α A x α + 0 x S h , α A x h α B u ϖ u y d + 0 x S h , α A x h α f , ϖ f , y d x 1 α ε α E α , α A x 1 α + x 1 α B M 1 M 2 α E α , α A x 1 α ϖ y C Ω + M 2 ϖ y C Ω = x 1 α ε α E α , α A x 1 α + M 2 1 + x 1 α B M 1 α E α , α A x 1 α ϖ y C Ω .
So,
ϖ y C Ω x 1 α ε α 1 ρ E α , α A x 1 α ,
where
ρ : = M 2 1 + x 1 α B M 1 α E α , α A x 1 α .
Thus,
ϖ x y x N ε , N = x 1 α α 1 ρ E α , α A x 1 α .
This completes the proof. □
Remark 6.
Let α = 2 in (2). Then, Theorem 3 coincides with the conclusion of Theorem 3 in [16].
Remark 7.
We note that Theorems 1–3 improve, extend, and complement some existing results in [16,19,21,35].

6. Examples

In this section, we present applications of the results derived.
Example 1.
Consider the following linear delay fractional controlled system:
C D 0 + 1.5 y x + A y x 0.5 = B u ( x ) , for x Ω : = 0 , 1 , y x ψ x , y x ψ x for 0.5 x 0 ,
where
A = 1 2 0 1 , B = 1 2 , ψ x = 2 x x , ψ x = 2 1 .
Wenote that B R 2 × 1 and u x R show the control vector. Constructing the corresponding fractional delay Gramian matrix of (25) via (10), we obtain
W 0.5 , 1.5 0 , 1 = 0 1 S 0.5 , 1.5 A 0.5 1.5 B B T S 0.5 , 1.5 A T 0.5 1.5 d = : O 1 + O 2 ,
where
O 1 = 0 0.5 S 0.5 , 1.5 A 0.5 1.5 B B T S 0.5 , 1.5 A T 0.5 1.5 d ,
for 0.5 0 , 0.5 ,
O 2 = 0.5 1 S 0.5 , 1.5 A 0.5 1.5 B B T S 0.5 , 1.5 A T 0.5 1.5 d ,
for 0.5 0.5 , 0 , and
H 0.5 , 1.5 A x 1.5 : = Θ , < x < 0.5 , I , 0.5 x < 0 , I A x 1.5 Γ 2.5 0 x < 0.5 , I A x 1.5 Γ 2.5 + A 2 x 0.5 3 Γ 4 , 0.5 x < 1 ,
M 0.5 , 1.5 A x 1.5 : = Θ , < x < 0.5 , I x + 0.5 , 0.5 x < 0 , I x + 0.5 A x 2.5 Γ 3.5 , 0 x < 0.5 , I x + 0.5 A x 2.5 Γ 3.5 + A 2 x 0.5 4 Γ 5 , 0.5 x < 1 ,
and
S 0.5 , 1.5 A x 1.5 : = Θ , < x < 0.5 , I x + 0.5 0.5 Γ 1.5 , 0.5 x < 0 , I x + 0.5 0.5 Γ 1.5 A x 2 Γ 3 , 0 x < 0.5 , I x + 0.5 0.5 Γ 1.5 A x 2 Γ 3 + A 2 x 0.5 3.5 Γ 4.5 , 0.5 x < 1 .
Next, we can calculate that
O 1 = 50961 400000 503601 1000000 503601 1000000 57031 500000 , O 2 = 3183 20000 3183 10000 3183 10000 3183 5000 .
Then, we obtain
W 0.5 , 1.5 0 , 1 = O 1 + O 2 = 114621 400000 821901 1000000 821901 1000000 261269 500000 ,
and
W 0.5 , 1.5 1 0 , 1 = 130634500000 131446670889 68491750000 43815556963 68491750000 43815556963 23879375000 43815556963 .
Therefore, we see that W 0.5 , 1.5 0 , 1 is positive definite. Furthermore, for any finite terminal conditions y 1 , y 1 R 2 such that y x 1 = y 1 = y 11 , y 12 T , y x 1 = y 1 = y 11 , y 12 T , as a result we can establish the corresponding control as follows:
u x = B T S 0.5 , 1.5 A T 0.5 x 1.5 W 0.5 , 1.5 1 0 , 1 β ,
where
β = y 1 M 0.5 , 1.5 A 0.5 1.5 ψ 0 + A 0.5 0 S 0.5 , 1.5 A 1.5 ψ d = y 11 1127 500 y 12 5337 5000 .
Hence, the system (25) is controllable on 0 , 1 by Theorem 1.
Example 2.
Consider the following nonlinear delay fractional controlled system:
C D 0 + 1.8 y x + A y x 0.6 = f x , y x + B u ( x ) , for x Ω 1 : = 0 , 1.2 , y x ψ x , y x ψ x for 0.6 x 0 ,
where
A = 1 0 0 2 , B = I 2 × 2 , ψ x = 3 x + 1 x 2 , ψ x = 3 2 x ,
f x , y x = 0.5 x 0.6 cos y 1 x 0.5 x 0.6 cos y 2 x .
Now, we set u ( x ) = y ˜ , where y ˜ R 2 . From the definition of Υ in ( G 2 ) , we obtain
Y = 0 1.2 S 0.6 , 1.8 A 0.6 1.8 B d y ˜ = 0 0.6 S 0.6 , 1.8 A 0.6 1.8 d y ˜ + 0.6 1.2 S 0.6 , 1.8 A 0.6 1.8 d y ˜ = 578469 1000000 0 0 141647 250000 y ˜ + 23783 100000 0 0 23783 100000 y ˜ = 816299 1000000 0 0 402209 500000 y ˜ .
Define the inverse Y 1 : R 2 L 2 Ω 1 , R 2 by
Y 1 y ˜ x : = 1000000 816299 0 0 500000 402209 y ˜ .
Then, we obtain
Y 1 y ˜ x 1000000 816299 0 0 500000 402209 y ˜ = 1.243 y ˜ ,
and thus we obtain Y 1 1.243 = : M 1 . Hence, the assumption ( G 2 ) is satisfied by Y.
Next, keep in mind that cos λ cos δ λ δ , for all λ, δ R , we obtain
f x , y f x , z = 0.5 x 0.6 cos y 1 x cos z 1 x 2 + cos y 2 x cos z 2 x 2 0.5 x 0.6 y 1 x z 1 x 2 + y 2 x z 2 x 2 = 0.5 x 0.6 y z ,
for all x Ω 1 , and y x , z x R 2 . We set L f x = 0.5 x 0.6 such that L f L q Ω 1 , R + in ( G 1 ) . By choosing p = q = 2 , we have
L f L 2 Ω 1 , R + = 0 1.2 0.5 0.6 2 d 1 2 = 0.18974 .
Then, we obtain
M 2 = 1.2 1.3 2.6 1 2 E 1.8 , 1.8 2 1.2 1.8 L f L q Ω , R + = 0.301 .
Finally, we calculate that
M 2 1 + M 1 1.2 1.8 1.8 E 1.8 , 1.8 2 1.2 1.8 B = 0.8815 < 1 ,
which implies that all the conditions of Theorems 2 and 3 are satisfied. Therefore, the system (26) is controllable and Hyers–Ulam stable.

7. Conclusions

In this work, we established some sufficient and necessary conditions for the controllability of linear fractional-delay systems by using a fractional delay Gramian matrix and the representation of solutions of these systems with the help of their delayed Mittag–Leffler matrix functions. Furthermore, we established some sufficient conditions for the controllability and Hyers–Ulam stability of nonlinear fractional-delay systems by applying Krasnoselskii’s fixed-point theorem and the representation of the solutions of these systems. Finally, the effectiveness of the obtained results was illustrated by numerical examples.
Our future work includes extending and complementing the results of this paper to derive the controllability and Hyers–Ulam stability results of fractional stochastic delay systems with compact analytic semigroups or using the delayed Mittag–Leffler matrix functions with various behaviors such as impulses and delays in multi-states.

Author Contributions

Conceptualization, B.A., X.W. and A.M.E.; data curation, B.A., X.W. and A.M.E.; formal analysis, B.A., X.W. and A.M.E.; software, A.M.E.; supervision, X.W.; validation, B.A., X.W. and A.M.E.; visualization, B.A., X.W. and A.M.E.; writing—original draft, A.M.E.; writing—review and editing, B.A., X.W. and A.M.E.; investigation, X.W. and A.M.E.; methodology, B.A., X.W. and A.M.E.; funding acquisition, B.A. All authors have read and agreed to the published version of the manuscript.

Funding

Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R216), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors sincerely appreciate the editors and anonymous referees for their carefully reading and helpful comments for improving this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Coimbra, C.F.M. Mechanics with variable-order differential operators. Ann. Phys. 2003, 12, 692–703. [Google Scholar] [CrossRef]
  2. Shi, K.; Cai, X.; She, K.; Zhong, S.; Soh, Y.; Kwon, O. Quantized memory proportional–integral control of active power sharing and frequency regulation in island microgrid under abnormal cyber attacks. Appl. Energy 2022, 322, 119540. [Google Scholar] [CrossRef]
  3. Diethelm, K. The Analysis of Fractional Differential Equations; Springer: Berlin, Germany, 2010. [Google Scholar]
  4. Heymans, N.; Podlubny, I. Physical interpretation of initial conditions for fractional differential equations with Riemann-Liouville fractional derivatives. Rheol. Acta 2006, 45, 765–771. [Google Scholar] [CrossRef] [Green Version]
  5. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  6. Obembe, A.D.; Hossain, M.E.; Abu-Khamsin, S.A. Variable-order derivative time fractional diffusion model for heterogeneous porous media. J. Pet. Sci. Eng. 2017, 152, 391–405. [Google Scholar] [CrossRef]
  7. Sweilam, N.H.; Al-Mekhlafi, S.M. Numerical study for multi-strain tuberculosis (TB) model of variable-order fractional derivatives. J. Adv. Res. 2016, 7, 271–283. [Google Scholar] [CrossRef] [Green Version]
  8. Tarasov, V. Handbook of Fractional Calculus with Applications; Applications in Physics, Part A; de Gruyter: Berlin, Germany, 2019. [Google Scholar]
  9. Khusainov, D.Y.; Shuklin, G.V. Linear autonomous time-delay system with permutation matrices solving. Stud. Univ. Zilina. Math. Ser. 2003, 17, 101–108. [Google Scholar]
  10. Khusainov, D.Y.; Diblík, J.; Růžičková, M.; Lukáčová, J. Representation of a solution of the Cauchy problem for an oscillating system with pure delay. Nonlinear Oscil. 2008, 11, 276–285. [Google Scholar] [CrossRef]
  11. Elshenhab, A.M.; Wang, X.T. Representation of solutions for linear fractional systems with pure delay and multiple delays. Math. Meth. Appl. Sci. 2021, 44, 12835–12850. [Google Scholar] [CrossRef]
  12. Elshenhab, A.M.; Wang, X.T. Representation of solutions of linear differential systems with pure delay and multiple delays with linear parts given by non-permutable matrices. Appl. Math. Comput. 2021, 410, 126443. [Google Scholar] [CrossRef]
  13. Elshenhab, A.M.; Wang, X.T. Representation of solutions of delayed linear discrete systems with permutable or nonpermutable matrices and second-order differences. RACSAM Rev. R. Acad. Cienc. Exactas Fís. Nat. Ser. A Mat. 2022, 116, 58. [Google Scholar] [CrossRef]
  14. Elshenhab, A.M.; Wang, X.T.; Mofarreh, F.; Bazighifan, O. Exact solutions and finite time stability of linear conformable fractional systems with pure delay. Comput. Model. Eng. Sci. 2023, 134, 927–940. [Google Scholar] [CrossRef]
  15. Elshenhab, A.M.; Wang, X.T.; Bazighifan, O.; Awrejcewicz, J. Finite-time stability analysis of linear differential systems with pure delay. Mathematics 2022, 10, 1359. [Google Scholar] [CrossRef]
  16. Elshenhab, A.M.; Wang, X.T. Controllability and Hyers–Ulam stability of differential systems with pure delay. Mathematics 2022, 10, 1248. [Google Scholar] [CrossRef]
  17. Elshenhab, A.M.; Wang, X.T.; Cesarano, C.; Almarri, B.; Moaaz, O. Finite-Time Stability Analysis of fractional-delay systems. Mathematics 2022, 10, 1883. [Google Scholar] [CrossRef]
  18. Huseynov, I.T.; Mahmudov, N.I. Delayed analogue of three-parameter Mittag–Leffler functions and their applications to Caputo-type fractional time delay differential equations. Math. Meth. Appl. Sci. 2020, 1–25. [Google Scholar] [CrossRef]
  19. Liu, L.; Dong, Q.; Li, G. Exact solutions and Hyers–Ulam stability for fractional oscillation equations with pure delay. Appl. Math. Lett. 2021, 112, 106666. [Google Scholar] [CrossRef]
  20. Nawaz, M.; Jiang, W.; Sheng, J. The controllability of nonlinear fractional differential system with pure delay. Adv. Differ. Equ. 2020, 2020, 183. [Google Scholar] [CrossRef]
  21. Liang, C.; Wang, J.; O’Regan, D. Controllability of nonlinear delay oscillating systems. Electron. J. Qual. Theory Differ. Equ. 2017, 47, 1–18. [Google Scholar] [CrossRef]
  22. Klamka, J. Controllability of Dynamical Systems; Kluwer Academic: Dordrecht, The Netherlands, 1993. [Google Scholar]
  23. Nirmala, R.J.; Balachandran, K.; Rodríguez-Germa, L.; Trujillo, J.J. Controllability of nonlinear fractional delay dynamical systems. Rep. Math. Phys. 2016, 77, 87–104. [Google Scholar] [CrossRef]
  24. Dineshkumar, C.; Udhayakumar, R.; Vijayakumar, V.; Shukla, A.; Nisar, K.S. A note on approximate controllability for nonlocal fractional evolution stochastic integrodifferential inclusions of order r∈1,2 with delay. Chaos Solitons Fractals 2021, 153, 111565. [Google Scholar] [CrossRef]
  25. Balachandran, K.; Park, J.Y.; Trujillo, J.J. Controllability of nonlinear fractional dynamical systems. Nonlinear Anal. 2012, 75, 1919–1926. [Google Scholar] [CrossRef]
  26. Balachandran, K.; Govindaraj, V.; Rodriguez-Germa, L.; Trujillo, J.J. Controllability results for nonlinear fractional-order dynamical systems. J. Optim. Theory Appl. 2012, 156, 33–44. [Google Scholar] [CrossRef]
  27. Balachandran, K.; Govindaraj, V.; Rivero, M.; Trujillo, J.J. Controllability of fractional damped dynamical systems. Appl. Math. Comput. 2015, 257, 66–73. [Google Scholar] [CrossRef]
  28. Sakthivel, R.; Ganesh, R.; Ren, Y.; Anthoni, S.M. Approximate controllability of nonlinear fractional dynamical systems. Commun. Nonlinear Sci. Numer. Simul. 2013, 18, 3498–3508. [Google Scholar] [CrossRef]
  29. Li, M.; Debbouche, A.; Wang, J. Relative controllability in fractional differential equations with pure delay. Math. Methods Appl. Sci. 2018, 41, 8906–8914. [Google Scholar] [CrossRef]
  30. Agarwal, R.; Hristova, S.; O’Regan, D. Existence and Ulam type stability for nonlinear Riemann–Liouville fractional differential equations with constant delay. Electron. J. Qual. Theory Differ. Equ. 2020, 67, 1–18. [Google Scholar] [CrossRef]
  31. Liu, K.; Wang, J.; Zhou, Y.; O’Regan, D. Hyers–Ulam stability and existence of solutions for fractional differential equations with Mittag–Leffler kernel. Chaos Solitons Fractals 2020, 132, 109534. [Google Scholar] [CrossRef]
  32. Sharma, J.P.; George, R.K. Controllability of matrix second order systems: A trigonometric matrix approach. Electron. J. Diff. Equ. 2007, 80, 1–14. [Google Scholar]
  33. Wang, J.; Lv, L.; Zhou, Y. Ulam stability and data dependence for fractional differential equations with Caputo derivative. Electron. J. Qual. Theory Differ. Equ. 2011, 63, 1–10. [Google Scholar]
  34. Smart, D.R. Fixed-Point Theorems; University Press: Cambridge, UK, 1980. [Google Scholar]
  35. Li, M.; Wang, J.R. Finite time stability of fractional delay differential equations. Appl. Math. Lett. 2017, 64, 170–176. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Almarri, B.; Wang, X.; Elshenhab, A.M. Controllability and Hyers–Ulam Stability of Fractional Systems with Pure Delay. Fractal Fract. 2022, 6, 611. https://doi.org/10.3390/fractalfract6100611

AMA Style

Almarri B, Wang X, Elshenhab AM. Controllability and Hyers–Ulam Stability of Fractional Systems with Pure Delay. Fractal and Fractional. 2022; 6(10):611. https://doi.org/10.3390/fractalfract6100611

Chicago/Turabian Style

Almarri, Barakah, Xingtao Wang, and Ahmed M. Elshenhab. 2022. "Controllability and Hyers–Ulam Stability of Fractional Systems with Pure Delay" Fractal and Fractional 6, no. 10: 611. https://doi.org/10.3390/fractalfract6100611

APA Style

Almarri, B., Wang, X., & Elshenhab, A. M. (2022). Controllability and Hyers–Ulam Stability of Fractional Systems with Pure Delay. Fractal and Fractional, 6(10), 611. https://doi.org/10.3390/fractalfract6100611

Article Metrics

Back to TopTop