Next Article in Journal
Detection of Internal Wire Broken in Mining Wire Ropes Based on WOA–VMD and PSO–LSSVM Algorithms
Next Article in Special Issue
Exotic Particle Dynamics Using Novel Hermitian Spin Matrices
Previous Article in Journal
A Time-Fractional Differential Inequality of Sobolev Type on an Annulus
Previous Article in Special Issue
Existence of Multiple Weak Solutions to a Discrete Fractional Boundary Value Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exponential Stability and Relative Controllability of Nonsingular Conformable Delay Systems

School of Mathematics and Statistics, Guizhou University of Finance and Economics, Guiyang 550025, China
Axioms 2023, 12(10), 994; https://doi.org/10.3390/axioms12100994
Submission received: 20 September 2023 / Revised: 10 October 2023 / Accepted: 17 October 2023 / Published: 20 October 2023

Abstract

:
In this paper, we investigate a delayed matrix exponential and utilize it to derive a representation of solutions to a linear nonsingular delay problem with permutable matrices. To begin with, we present a novel definition of α -exponential stability for these systems. Subsequently, we put forward several adequate conditions to ensure the α -exponential stability of solutions for such delay systems. Moreover, by constructing a Grammian matrix that accounts for delays, we provide a criterion to determine the relative controllability of a linear problem. Additionally, we extend our analysis to nonlinear problems. Lastly, we offer several examples to verify the effectiveness of our theoretical findings.

1. Introduction

Delay systems are common in various engineering applications, ranging from chemical processes to mechanical systems. The presence of time delays in these systems can cause instability and oscillatory behavior, which can be destructive to the performance and reliability of the system. Consequently, the investigation into the stability and controllability of time-delay systems has garnered substantial research attention over the years. In [1], the concept of delayed matrix exponential was introduced by the authors. They utilized this concept to derive a representation of solutions for a linear delay problem, specifically when permutable matrices are taken into account. In a subsequent work by [2], by using the concepts presented in [1], the authors extended those ideas to a discrete matrix delayed exponential function and investigated the representation of solutions to linear discrete delay systems. Motivated by the works of [1,2], numerous scholars have focused on investigating a wide range of delay systems involving permutable matrices. For a comprehensive understanding of delay systems, we refer interested readers to the following papers: [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]. These references delve into various aspects related to the aforementioned systems.
Exponential stability is a very important property of a control system. It is related to the behavior of the system under continuous parameter variations. As for exponential stability, it is crucial for linear control systems because it ensures that the output signals remain approximately constant in the presence of noise and parameter variations. Here, we present a new definition of exponential stability for the investigated systems [4,5,7,11,13,16,18,21,22,23,24].
Control theory has attracted a lot of research interest for many years, with the aim of designing efficient and robust control strategies for various engineering applications. Relative controllability is an important concept in control theory, which measures the ability of a control input to guide the system from one state to another. Specifically, we examine the relative controllability of a conformable delay system by employing the concept of null controllability. Null controllability is a fundamental property of linear systems that has been widely studied in the literature and provides a powerful tool for analyzing the controllability of nonlinear systems [3,8,13,14,15,17,25,26,27,28,29].
In this paper, we build upon the concepts introduced in [1,3,5] to investigate the α -exponential stability and relative controllability of investigated systems. Initially, we focus on the α -exponential stability analysis of systems
B D α 0 z ( ζ ) = F z ( ζ ) + P z ( ζ κ ) , ζ 0 , κ 0 , z ( ζ ) = ψ ( ζ ) , κ ζ 0 ,
and
B D α 0 z ( ζ ) = F z ( ζ ) + P z ( ζ κ ) + m ( ζ , z ( ζ ) ) , ζ 0 , κ 0 , z ( ζ ) = ψ ( ζ ) , κ ζ 0 ,
and
B D α 0 z ( ζ ) = F z ( ζ ) + P z ( ζ κ ) + l ( ζ , z ( ζ κ ) ) , ζ 0 , κ 0 , z ( ζ ) = ψ ( ζ ) , κ ζ 0 ,
here D α 0 ( 0 < α < 1 ) represents the conformable derivative with lower index zero (see Definition 1 [30]), B , F , P R n × n are constant permutable matrices and B is a nonsingular matrix and ψ C 1 ( [ κ , 0 ] , R n ) . Together with, m , l C ( [ 0 , + ] × R n , R n ) . Also we are going to discuss the relatively controllability of
B D α 0 z ( ζ ) = F z ( ζ ) + P z ( ζ κ ) + y ( ζ , z ( ζ ) ) + J v ( ζ ) , ζ I , κ 0 , z ( ζ ) = ψ ( ζ ) , κ ζ 0 ,
where I : = [ 0 , ζ 1 ] , ζ 1 > 0 , J R n , y C ( I × R n , R n ) and the control function v ( · ) takes values from L 2 ( I , R n ) .
We investigate sufficient conditions to ensure the α -exponential stability of solutions to Equations (1)–(3) using Gronwall inequality techniques. Additionally, we explore the application of a delay Grammian matrix in establishing both sufficient and necessary conditions for relative controllability in linear delayed systems. The discussion also extends to nonlinear problems, employing Krasnoselskii’s fixed point theorem.
The remainder of this paper is structured as follows. In Section 2, we provide relevant notations, concepts, and lemmas, along with a representation of solutions to investigated systems. Section 3 is dedicated to the study of α -exponential stability of solutions for Equations (1)–(3). The relative controllability of system (4) is investigated in Section 4. Finally, we present numerical examples in the Section 5 to illustrate our main findings.

2. Preliminaries

Let R n represent the n-dimensional Euclidean space with the vector norm · , and R n × n denote the n × n matrix space with real-valued elements. For z R n and F R n × n , we define z = max 1 i n | z i | and F = max 1 i n j = 1 n | f i j | , here, z i z and f i j F . Let z lies in Banach space C ( I , R n ) , z C = sup ζ I z ( ζ ) . Furthermore, we define the norm ψ C = sup ζ [ κ , 0 ] ψ ( ζ ) . Let Y 1 , Y 2 be two Banach spaces, L b ( Y 1 , Y 2 ) stands for the space of all bounded linear operators from Y 1 to Y 2 . Additionally, L g ( I , Y 2 ) stands for the Banach space of functions x : I Y 2 that are Bochner integrable, with the norm y L g ( I , Y 2 ) (where 1 < g < ).
Consider the matrices Θ and I , where Θ represents the zero matrix and I represents the identity matrix. Similar to Lemma 2.1 (You and Wang 2019 [21]), we introduce the following lemma without proof.
Lemma 1.
Suppose Ω lie in a constant n × n matrix. If Ω β e β κ , β R + , then
e κ , α Ω ( ζ κ ) e β α ζ α , ζ R ,
where
e κ , α Ω ζ = Θ , < ζ < κ , κ > 0 I , κ ζ 0 , I + Ω ζ α α + Ω 2 1 2 ! ( ζ κ ) α α 2 + Ω 3 1 3 ! ( ζ 2 κ ) α α 3 + + Ω k 1 k ! ( ζ ( k 1 ) κ ) α α k , ( k 1 ) κ ζ < k κ , k = 1 , 2 , ,
(this is referred to as the delayed matrix exponential (see [1])).
Remark 1.
Clearly, (5) is reflated to α-time power function, which is different from 1-time power function in Lemma 2.1 (You and Wang 2019 [21]).
From (6), the uniform continuity of e κ , α Ω · on a compact interval can be observed easily. By employing the transformation x = B z in Equation (1), we obtain
D α 0 x ( ζ ) = F B 1 x ( ζ ) + P B 1 x ( ζ κ ) , ζ 0 , κ 0 , x ( ζ ) = B ψ ( ζ ) , κ ζ 0 ,
Using Theorem 3.4 [31], for (7), we have
x ( ζ ) = P ( ζ ) e F B 1 κ α α B ψ ( κ ) + κ 0 P ( ζ κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s ,
where P ( ζ ) = e F B 1 ζ α α e κ , α P 1 B 1 ζ and P 1 = e F B 1 ( ζ κ ) α ζ α α P satisfying initial condition ([31]) P ( ζ ) = e A ( ζ ) α α , κ ζ < 0 , where
e A ζ α ζ : = e A ζ α α , ζ 0 , e A ( ζ ) α α , κ ζ 0 .
Now substituting x = B z into (8), by solving it, we can observe that any solution of Equation (1) takes the form
z ( ζ ) = P ( ζ ) e F B 1 κ α α ψ ( κ ) + B 1 κ 0 P ( ζ κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s ,
where F B = B F , F P = P F , P B = B P are employed.
Likewise, by utilizing the transformation x = B z in
B D α 0 z ( ζ ) = F z ( ζ ) + P z ( ζ κ ) + ϖ ( ζ ) , ζ 0 , κ 0 , z ( ζ ) = ψ ( ζ ) , κ ζ 0 ,
one has
D α 0 x ( ζ ) = F B 1 x ( ζ ) + P B 1 x ( ζ κ ) + ϖ ( ζ ) , ζ 0 , κ 0 , x ( ζ ) = B ψ ( ζ ) , κ ζ 0 ,
(here x : [ 0 , + ) R n is continuous).
Using Theorem 3.5 [31], for (12), we have
x ( ζ ) = P ( ζ ) e F B 1 κ α α B ψ ( κ ) + κ 0 P ( ζ κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s + 0 ζ P ( ζ κ s ) e F B 1 κ α α ϖ ( s ) s α 1 d s .
Now substituting x = B z into (13), upon solving it, it becomes evident that the form of any solution to Equation (11) is given by
z ( ζ ) = P ( ζ ) e G B 1 κ α α ψ ( κ ) + B 1 κ 0 P ( ζ κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s + B 1 0 ζ P ( ζ κ s ) e F B 1 κ α α ϖ ( s ) s α 1 d s .
Likewise, the solution to Equation (2) is given by
z ( ζ ) = P ( ζ ) e F B 1 κ α α ψ ( κ ) + B 1 κ 0 P ( ζ κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s + B 1 0 ζ P ( ζ κ s ) e F B 1 κ α α m ( s , z ( s ) ) s α 1 d s ,
and the form of any solution to Equation (3) is given by
z ( ζ ) = P ( ζ ) e F B 1 κ α α ψ ( κ ) + B 1 κ 0 P ( ζ κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s + B 1 0 ζ P ( ζ κ s ) e F B 1 κ α α l ( s , z ( s κ ) ) s α 1 d s .
Based on the representation (15) since D α 0 ψ exists (i.e., ψ D α 0 ( [ κ , 0 ] , R n ) ), here is the presented definition.
Definition 1.
The α-exponential stability is attributed to the trivial solution of (1) (or (2), (3)), provided that positive constants c 1 , c 2 , δ exist, depending on B , F , P and ψ 1 , where ψ 1 : = max [ κ , 0 ] ψ + max [ κ , 0 ] D α 0 ψ , catering to z ( ζ ) c 1 e c 2 ζ α , ζ 0 , for ψ 1 < δ , here, z represents any solution to (1) (or (2), (3)).
Definition 2
(see [3]). Equation (4) to be relatively controllable, suppose there exists an arbitrary initial vector function ψ C 1 ( [ κ , 0 ] , R n ) , and a final state of the vector z 1 R n at ζ 1 . If there is a control v L 2 ( I , R n ) such that Equation (4) has a solution z C ( [ κ , ζ 1 ] , R n ) satisfying the boundary condition z ( ζ 1 ) = z 1 .

3. α -Exponential Stability

This section focuses on examining the α -exponential stability of solutions in conformable systems. We take into account the following hypotheses.
( H 1 ) Set σ ( F B 1 ) = { λ 1 , λ 2 , · · · , λ n } be the eigenvalues of F B 1 with
R e λ 1 R e λ 2 · · · R e λ n k < 0 , k > 0 ,
i.e., for 0 < α < 1 , there exist W , w > 0 such that
e F B 1 ζ α α W e w ζ α α f o r a l l ζ 0 .
( H 2 ) Suppose P 1 B 1 β e β κ for β R + .
( H 3 ) For ζ 0 and z R n , there is a positive constant L such that the inequality m ( ζ , z ) L z holds.
( H 4 ) For ζ 0 and z R n , there is a positive constant Q such that the inequality l ( ζ , z ) Q z holds.
( H 5 ) Suppose β w < 0 .
( H 6 ) Suppose B 1 W 2 L e w κ α α + β w < 0 .
( H 7 ) Suppose B 1 W 2 Q e w κ α α + β w < 0 .
Lemma 2.
Under the hypotheses ( H 1 ) and ( H 2 ) , the inequality
P ( ζ ) e F B 1 κ α α W 2 e β κ α α e ( β w ) ζ α α , β , ζ R + .
holds.
Proof. 
From ( H 1 ) and ( H 2 ) via Lemma 1, we have
P ( ζ ) e F B 1 κ α α e F B 1 ζ α α e κ , α P 1 B 1 ζ e F B 1 κ α α W e w ζ α α e β ( ζ + κ ) α α W e w κ α α W 2 e ( β w ) ζ α + β | ( ζ + κ ) α ζ α | α W 2 e β κ α α e ( β w ) ζ α α .
Thus, the proof is concluded. □
At present, we are prepared to provide adequate conditions for the α -exponential stability of solutions in investigated systems.
Theorem 1.
Under the hypotheses ( H 1 ) , ( H 2 ) and ( H 5 ) , the solution of Equation (1) exhibits α-exponential stability.
Proof. 
It is worth noting that the solution of (1) takes the form (10). Now from (10) we get
z ( ζ ) W 2 e β κ α α e ( β w ) ζ α α ψ ( κ ) + B 1 κ 0 W 2 e w κ α α e ( β w ) ( ζ s ) α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s W 2 e β κ α α e ( β w ) ζ α α ψ ( κ ) + B 1 κ 0 W 2 e w κ α α e ( β w ) [ ( ζ s ) α ζ α ] + ( β w ) ζ α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s W 2 e β κ α α e ( β w ) t α α ψ ( κ ) + B 1 κ 0 W 2 e w κ α α e ( β w ) | ( ζ s ) α ζ α | + ( β w ) ζ α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s
W 2 e β κ α α e ( β w ) ζ α α ψ ( κ ) + B 1 κ 0 W 2 e w κ α α e ( β w ) ( s ) α + ( β w ) ζ α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s U ( ψ , D α 0 ψ ) e ( β w ) ζ α α ,
where
U ( ψ , D α 0 ψ ) = W 2 e β κ α α ψ ( κ ) + B 1 κ 0 e w κ α α e ( β w ) ( s ) α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s > 0 .
Now if we choose D α 0 ψ small enough, then 0 M ( ψ , D α 0 ψ ) < δ by ( H 5 ) . □
Theorem 2.
Under the hypotheses ( H 1 ) , ( H 2 ) , ( H 3 ) and ( H 6 ) , the solution of Equation (2) exhibits α-exponential stability.
Proof. 
It is important to note that the solution of Equation (2) takes the form (14). Let t 0 . From (14), using ( H 1 ) ( H 3 ) via (16), we get
z ( ζ ) P ( ζ ) e F B 1 κ α α ψ ( κ ) + B 1 κ 0 P ( ζ κ s ) e F B 1 κ α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 0 ζ P ( ζ κ s ) e F B 1 κ α α m ( s , z ( s ) ) s α 1 d s W 2 e β κ α α e ( β w ) t α α ψ ( κ ) + B 1 κ 0 W 2 e w κ α α e ( β w ) ζ s ) α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 0 ζ W 2 e w κ α α e ( β w ) ( ζ s ) α α L | | z ( s ) | | s α 1 d s W 2 e β κ α α e ( β w ) ζ α α ψ ( κ ) + B 1 κ 0 W 2 e w κ α α e ( β w ) ( s ) α + ( β w ) ζ α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 0 ζ W 2 e w κ α α e ( β w ) ( s ) α ] + ( β w ) ζ α α L | | z ( s ) | | s α 1 d s W 2 e β κ α α e ( β w ) ζ α α ψ ( κ ) + B 1 κ 0 W 2 e w κ α α e ( β w ) ζ α α e ( β w ) ( s ) α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 0 ζ W 2 e w κ α α e ( β w ) ζ α α e ( β w ) ( s ) α α L | | z ( s ) | | s α 1 d s W 2 e β κ α α e ( β w ) ζ α α ψ ( κ ) + B 1 W 2 e w κ α α e ( β w ) ζ α α κ 0 e ( β w ) ( s ) α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 L W 2 e w κ α α e ( β w ) ζ α α 0 ζ e ( β w ) ( s ) α α | | z ( s ) | | s α 1 d s ,
so,
e ( w β ) ζ α α z ( ζ ) U ( ψ , D α 0 ψ ) + B 1 L W 2 e w κ α α 0 ζ e ( β w ) ( s ) α α z ( s ) s α 1 d s .
By employing the classical Gronwall inequality and Equation (9), we obtain
e ( w β ) ζ α α z ( ζ ) U ( ψ , D α 0 ψ ) e B 1 W 2 L e w κ α α ζ α α ,
this results in
z ( ζ ) U ( ψ , D α 0 ψ ) e ( B 1 W 2 L e w κ α α + β w ) ζ α α .
Based on ( H 6 ) , the desired outcome is valid. □
Theorem 3.
Under the hypotheses ( H 1 ) , ( H 2 ) , ( H 4 ) and ( H 7 ) , the solution of Equation (3) exhibits α-exponential stability.
Proof. 
It is known that the solution of Equation (3) can be expressed in the form of (15). Let ζ 0 , from (15), utilizing ( H 1 ) , ( H 2 ) and ( H 4 ) via (16), we get
z ( ζ ) W 2 e β κ α α e ( β w ) ζ α α ψ ( κ ) + B 1 W 2 e w κ α α κ 0 e ( β w ) [ ( ζ s ) α ζ α ] + ( β w ) ζ α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 W 2 Q e w κ α α 0 ζ e ( β w ) [ ( ζ s ) α ζ α ] + ( β w ) ζ α α | | z ( s κ ) | | s α 1 d s W 2 e β κ α α e ( β w ) ζ α α ψ ( κ ) + B 1 W 2 e w κ α α κ 0 e ( β w ) ( s ) α + ( β w ) ζ α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 W 2 Q e w κ α α 0 ζ e ( β w ) ( s ) α + ( β w ) ζ α α | | z ( s κ ) | | s α 1 d s W 2 e β κ α α e ( β w ) ζ α α ψ ( κ ) + B 1 W 2 e w κ α α κ 0 e ( β w ) ζ α α e ( β w ) ( s ) α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 W 2 Q e w κ α α 0 ζ e ( β w ) ζ α α e ( β w ) ( s ) α α | | z ( s κ ) | | s α 1 d s W 2 e β κ α α e ( β w ) ζ α α ψ ( κ ) + B 1 W 2 e w κ α α e ( β w ) ζ α α κ 0 e ( β w ) ( s ) α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 W 2 Q e w κ α α e ( β w ) ζ α α 0 ζ e ( β w ) ( s ) α α | | z ( s κ ) | | s α 1 d s .
This suggests that
e ( w β ) ζ α α z ( ζ ) U ( ψ , D α 0 ψ ) + B 1 W 2 Q e w κ α α 0 ζ e ( β w ) ( s ) α α z ( s κ ) s α 1 d s .
Take note e ( w β ) ζ α α z ( ζ ) = e ( w β ) ζ α α ψ ( ζ ) , ζ [ κ , 0 ] . Applied Lemma 2.4 of [21] to (17), we can obtain
e ( w β ) ζ α α z ( ζ ) U ( ψ , D α 0 ψ ) + B 1 W 2 Q e w κ α α 0 κ e ( β w ) ( s ) α α ψ ( s κ ) s α 1 d s e B 1 W 2 Q e w κ α α ( ζ α α κ α α ) ,
which yields (from (9))
z ( ζ ) U * ( ψ , D α 0 ψ ) e ( B 1 W 2 Q e w κ α α + β w ) ζ α α .
where U * ( ψ , D α 0 ψ ) = [ U ( ψ , D α 0 ψ ) + B 1 W 2 Q e w κ α α 0 κ e ( β w ) ( s ) α α ψ ( s κ ) s α 1 d s ] e B 1 W 2 Q e w κ α α κ α α > 0 . It can be observed from ( H 7 ) that Equation (3) exhibits α -exponential stability. Therefore, the proof is concluded. □

4. Relative Controllability

We will examine the controllability relative to (4) in this section. It should be noted that any solution to (4) can be displayed in the pattern of
z ( ζ ) = P ( ζ ) e F B 1 κ α α ψ ( κ ) + B 1 κ 0 P ( ζ κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s + B 1 0 ζ P ( ζ κ s ) e F B 1 κ α α [ y ( s , z ( s ) ) + J v ( s ) ] s α 1 d s .

4.1. Linear Problem

Assuming y 0 , we can turn (4) into the system
B D α 0 z ( ζ ) = F z ( ζ ) + P z ( ζ κ ) + J v ( ζ ) , ζ I , κ 0 , z ( ζ ) = ψ ( ζ ) , κ ζ 0 .
The introduction of a delay Grammian matrix is taken the form:
M c [ 0 , ζ 1 ] = 0 ζ 1 B 1 P ( ζ 1 κ s ) e F B 1 κ α α J J e F B 1 κ α α P ( ζ 1 κ s ) B 1 s α 1 d s .
where · represents the matrix’s transpose.
Theorem 4.
The linear problem (19) exhibits relative controllability if and only if the matrix M c [ 0 , ζ 1 ] is nonsingular.
Proof. 
Sufficiency. As the matrix M c [ 0 , ζ 1 ] is nonsingular, the existence of its inverse is guaranteed.
A control function can be chosen in the following manner:
v ( ζ ) = J e F B 1 κ α α P ( ζ 1 κ ζ ) B 1 M c 1 [ 0 , ζ 1 ] η ,
where
η = z 1 P ( ζ 1 ) e F B 1 κ α α ψ ( κ ) B 1 κ 0 P ( ζ 1 κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s .
Then
z ( ζ 1 ) = P ( ζ 1 ) e F B 1 κ α α ψ ( κ ) + B 1 κ 0 P ( ζ 1 κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s + B 1 0 ζ 1 P ( ζ 1 κ s ) e F B 1 κ α α J J e F B 1 κ α α P ( ζ 1 κ s ) B 1 M c 1 [ 0 , ζ 1 ] η s α 1 d s = z 1 .
It is evident that the initial condition is satisfied as well. By applying Definition 2 and Equation (21), we can conclude that (19) exhibits relative controllability.
Necessity. We will present a proof by contradiction. Let us assume that M c [ 0 , ζ 1 ] is a singular matrix, implying the existence of at least a non-zero state z ˜ R n meeting
z ˜ M c [ 0 , ζ 1 ] z ˜ = 0 .
Moreover, one can achieve
0 = z ˜ M c [ 0 , ζ 1 ] z ˜ = 0 ζ 1 z ˜ B 1 P ( ζ 1 κ s ) e F B 1 κ α α J J e F B 1 κ α α P ( ζ 1 κ s ) B 1 z ˜ s α 1 d s = ( e F B 1 κ α α ) 2 0 ζ 1 z ˜ B 1 P ( ζ 1 κ s ) J 2 s α 1 d s ,
which implies that
z ˜ B 1 P ( ζ 1 κ s ) J = ( 0 , , 0 n ) : = 0 , s I .
Given that (19) exhibits relative controllability, according to Definition 2, there is one control v 1 ( ζ ) that can steer the initial state to 0 at ζ 1 , i.e.,
z ( ζ 1 ) = P ( ζ 1 ) e F B 1 κ α α ψ ( κ ) + B 1 κ 0 P ( ζ 1 κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s + B 1 0 ζ 1 P ( ζ 1 κ s ) e F B 1 κ α α J v 1 ( s ) s α 1 d s = 0 .
Likewise, there is one control v 2 ( ζ ) capable of directing the initial state to the state z ˜ at ζ 1 , which can be represented as
z ( ζ 1 ) = P ( ζ 1 ) e F B 1 κ α α ψ ( κ ) + B 1 κ 0 P ( ζ 1 κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s + B 1 0 ζ 1 P ( ζ 1 κ s ) e F B 1 κ α α J v 2 ( s ) s α 1 d s = z ˜ .
From (23) and (24), one has
z ˜ = B 1 0 ζ 1 P ( ζ 1 κ s ) e F B 1 κ α α J [ v 2 ( s ) v 1 ( s ) ] s α 1 d s .
Multiply both sides of (25) by z ˜ and via (22) we have
z ˜ z ˜ = 0 ζ 1 z ˜ B 1 P ( ζ 1 κ s ) e F B 1 κ α α J [ v 2 ( s ) v 1 ( s ) ] s α 1 d s = e F B 1 κ α α 0 ζ 1 z ˜ B 1 P ( ζ 1 κ s ) J [ v 2 ( s ) v 1 ( s ) ] s α 1 d s = 0 .
This leads to the conclusion that z ˜ =  0, which contradicts the non-zero nature of z ˜ . Hence, M c [ 0 , ζ 1 ] is nonsingular. Thus, the proof is concluded. □

4.2. Nonlinear Problem

Let us examine the following requirements:
[ A W ] : The operator M : L 2 ( I , R n ) R n defined as
M v = B 1 0 ζ 1 P ( ζ 1 κ s ) e F B 1 κ α α J v ( s ) s α 1 d s .
The operator M has an inverse, denoted as M 1 , which operates on the space L 2 ( I , R n ) / k e r M .
Set U 1 = M 1 L b ( R n , L 2 ( I , R n ) / k e r M ) . According to the Remark in Wang et al. (2017) [15], it can be inferred that
U 1 = M c 1 [ 0 , ζ 1 ] ,
here, M c 1 [ 0 , ζ 1 ] is determined to (20). According to Theorem 4, M c 1 [ 0 , ζ 1 ] will be well-defined if (19) exhibits relative controllability.
[ A F ] : The function y : I × R n R n is continuous and there is one constant h > 1 and L y ( · ) L h ( I , R + ) meeting
y ( ζ , z ) y ( ζ , v ) L y ( ζ ) z v , z , v R n .
Theorem 5.
Assuming that ( H 1 ) , ( H 2 ) , [ A W ] and [ A F ] hold true, it can be concluded that system (4) exhibits relative controllability under the condition that
b 1 + W 2 U 1 B 1 J e w κ α α e ( β w ) ζ 1 α α ζ 1 α α < 1 ,
where b = B 1 W 2 e w κ α α e ( β w ) ( ζ 1 + κ ) α α 1 g α g + 1 ζ 1 g α g + 1 1 g L y L h ( I , R + ) , and 1 g + 1 h = 1 , g , h > 1 .
Proof. 
By exploiting hypothesis [ A W ] for any z ( · ) C ( I , R n ) , defining one control function v z ( t ) as follows:
v z ( ζ ) = M 1 [ z 1 P ( ζ 1 ) e F B 1 κ α α ψ ( κ ) B 1 κ 0 P ( ζ 1 κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s B 1 0 ζ 1 P ( ζ 1 κ s ) e F B 1 κ α α y ( s , z ( s ) ) s α 1 d s ] ( ζ ) , ζ I .
We demonstrate that by employing this control, the operator T : C ( I , R n ) C ( I , R n ) defined as follows:
( T z ) ( ζ ) = P ( ζ ) e F B 1 κ α α ψ ( κ ) + B 1 κ 0 P ( ζ κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s + B 1 0 ζ P ( ζ κ s ) e F B 1 κ α α y ( s , z ( s ) ) s α 1 d s + B 1 0 ζ P ( ζ κ s ) e F B 1 κ α α J v z ( s ) s α 1 d s ,
The operator T has one fixed point z, which corresponds to a solution of (4).
We must verify ( T z ) ( ζ 1 ) = z 1 and ( T z ) ( 0 ) = z 0 , indicating that v z guides the system (4) from z 0 to z 1 within a finite time ζ 1 . This confirmation implies that system (4) exhibits relative controllability over I.
For every positive number r, consider the set S r = { z C ( I , R n ) : z C r } . Let N = sup ζ I y ( ζ , 0 ) . The proof will be divided into a couple of procedures.
  • Step 1. It is asserted that there is one positive constant r meeting T ( S r ) S r .
Note that
0 ζ s α 1 L y ( s ) d s 1 g α g + 1 ζ 1 g α g + 1 1 g L y L h ( I , R + ) ,
and
0 ζ y ( s , 0 ) s α 1 d s ζ 1 α α N .
Taking into account (28), and using ( H 1 ) , ( H 2 ) , [ A W ] and [ A F ] , we can get
v z M 1 L ζ 1 ( R n , L 2 ( I , R n ) / k e r M ) ( z 1 + P ( ζ 1 ) e F B 1 κ α α ψ ( κ ) + B 1 κ 0 P ( ζ 1 κ s ) e F B 1 κ α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 0 ζ 1 P ( ζ 1 κ s ) e F B 1 κ α α y ( s , z ( s ) ) s α 1 d s ) U 1 ( z 1 + W 2 e β κ α α e ( β w ) ζ 1 α α ψ ( κ ) + B 1 W 2 e w κ α α e ( β w ) ( ζ 1 + κ ) α α κ 0 B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 W 2 e w κ α α e ( β w ) ( ζ 1 + κ ) α α 0 ζ 1 L y ( s ) z ( s ) s α 1 d s + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α 0 ζ 1 y ( s , 0 ) s α 1 d s ) U 1 ( z 1 + W 2 e β κ α α e ( β w ) ζ 1 α α ψ ( κ ) + B 1 W 2 e w κ α α e ( β w ) ( ζ 1 + κ ) α α κ 0 B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 W 2 e w κ α α e ( β w ) ( ζ 1 + κ ) α α 1 g α g + 1 ζ 1 g α g + 1 1 g L y L h ( I , R + ) z C + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α N ζ 1 α α ) U 1 z 1 + U 1 a + U 1 b z C ,
where
a = W 2 e β κ α α e ( β w ) ζ 1 α α ψ ( κ ) + B 1 W 2 e w κ α α e ( β w ) ( ζ 1 + κ ) α α κ 0 B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α N ζ 1 α α .
Based on the assumptions [ A W ] and [ A F ] , we can derive
( T z ) ( ζ ) W 2 e β κ α α e ( β w ) ζ 1 α α ψ ( κ ) + B 1 W 2 e w κ α α e ( β w ) ( ζ 1 + κ ) α α κ 0 B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α 1 g α g + 1 ζ 1 g α g + 1 1 g L y L h ( I , R + ) z C + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α N ζ 1 α α + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α J [ U 1 z 1 + U 1 a + U 1 b z C ] ζ 1 α α
W 2 e β κ α α e ( β w ) ζ 1 α α ψ ( κ ) + B 1 W 2 e w κ α α e ( β w ) ( ζ 1 + κ ) α α κ 0 B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α 1 g α g + 1 ζ 1 g α g + 1 1 g L y L h ( I , R + ) z C + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α N ζ 1 α α + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α J U 1 z 1 ζ 1 α α + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α J U 1 a ζ 1 α α + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α J U 1 b z C ζ 1 α α a [ 1 + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α J U 1 ζ 1 α α ] + E 1 W 2 e w κ α α e ( β w ) ζ 1 α α J U 1 ζ 1 α α z 1 + b [ 1 + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α J U 1 ζ 1 α α ] r a 1 + U 1 J N a + U 1 J N a z 1 + b 1 + U 1 J N a r = r ,
where
r = a 1 + U 1 J N a + U 1 J N z 1 1 b 1 + U 1 J N a .
Thus, we deduce that T ( S r ) is a subset of S r for this particular value of r.
Throughout the remaining proof, the aforementioned value of r will be the one under consideration. Operators T 1 and T 2 are defined as follows:
( T 1 z ) ( ζ ) = P ( ζ ) e F B 1 κ α α ψ ( κ ) + B 1 κ 0 P ( ζ κ s ) e F B 1 κ α α [ B D α 0 ψ ( s ) F ψ ( s ) ] s α 1 d s + B 1 0 ζ P ( ζ κ s ) e F B 1 κ α α J v z ( s ) s α 1 d s , ( T 2 z ) ( ζ ) = B 1 0 ζ P ( ζ κ s ) e F B 1 κ α α y ( s , z ( s ) ) s α 1 d s ,
for ζ I .
  • Step 2. The operator T 1 exhibits the properties of a contraction map**.
Let z , ω S r . In view of [ A W ] and [ A F ] , for every ζ I , we obtain
u z ( ζ ) u ω ( ζ ) U 1 B 1 0 ζ 1 P ( ζ 1 κ s ) e F B 1 κ α α y ( s , z ( s ) ) y ( s , w ( s ) ) s α 1 d s U 1 B 1 0 ζ 1 W 2 e w κ α α e ( β w ) ( ζ s ) α α L y ( s ) z ( s ) ω ( s ) s α 1 d s U 1 b z ω C .
Thus,
( T 1 z ) ( ζ ) ( T 1 ω ) ( ζ ) B 1 0 ζ P ( ζ κ s ) e F B 1 κ α α J u z ( s ) u ω ( s ) s α 1 d s B 1 0 ζ W 2 e w κ α α e ( β w ) ( ζ s ) α α s α 1 d s J U 1 b z ω C B 1 W 2 e w κ α α e ( β w ) ζ 1 α α J U 1 b ζ 1 α α z ω C = b W 2 U 1 B 1 J e w κ α α e ( β w ) ζ 1 α α ζ 1 α α z ω C ,
so we obtain
T 1 z T 1 ω C Λ z ω C ,
where Λ = b W 2 U 1 B 1 J e w κ α α e ( β w ) ζ 1 α α ζ 1 α α .
Since (27) implies Λ < 1 , it follows that the operator T is a contraction map**.
  • Step 3. The operator T 2 : S r C ( I , R n ) is both compact and continuous.
Consider z n S r with z n z in S r as n . By utilizing [ A F ] , we obtain y ( · , z n ) y ( · , z ) in C ( I , R n ) as n and thus
( T 2 z n ) ( ζ ) ( T 2 z ) ( ζ ) B 1 0 ζ W 2 e w κ α α e ( β w ) ( ζ s ) α α y ( s , z n ( s ) ) y ( s , z ( s ) ) s α 1 d s B 1 W 2 e w κ α α e ( β w ) ζ 1 α α ζ 1 α α y ( s , z n ( s ) ) y ( s , z ( s ) ) C 0 w h i l e n .
This indicates that T 2 is continuous on S r .
To demonstrate the compactness of T 2 on S r , we need to establish that T 2 ( S r ) is both equicontinuous and bounded. For any z S r , ζ 1 ζ + ν ζ > 0 , we proceed with the following proof:
( T 2 z ) ( ζ + ν ) ( T 2 z ) ( ζ ) = B 1 ζ ζ + ν P ( ζ + ν κ s ) e F B 1 κ α α y ( s , z ( s ) ) s α 1 d s + B 1 0 ζ e F B 1 κ α α e F B 1 ( ζ + ν κ s ) α α ( e κ , α P 1 B 1 ( ζ + ν κ s ) e κ , α P 1 B 1 ( ζ κ s ) ) y ( s , z ( s ) ) s α 1 d s + B 1 0 ζ ( e F B 1 ( ζ + ν κ s ) α α e F B 1 ( ζ κ s ) α α ) e κ , α P 1 B 1 ( ζ κ s ) e F B 1 κ α α y ( s , z ( s ) ) s α 1 d s .
Let
Q 1 = B 1 ζ ζ + ν P ( ζ + ν κ s ) e F B 1 κ α α y ( s , z ( s ) ) s α 1 d s , Q 2 = B 1 0 ζ e F B 1 κ α α e F B 1 ( ζ + ν κ s ) α α ( e κ , α P 1 B 1 ( ζ + ν κ s ) e κ , α P 1 B 1 ( ζ κ s ) ) y ( s , z ( s ) ) s α 1 d s , Q 3 = B 1 0 ζ ( e F B 1 ( ζ + ν κ s ) α α e F B 1 ( ζ κ s ) α α ) e κ , α P 1 B 1 ( ζ κ s ) e F B 1 κ α α y ( s , z ( s ) ) s α 1 d s .
Based on the aforementioned analysis, we can conclude that
( T 2 z ) ( ζ + ν ) ( T 2 z ) ( ζ ) Q 1 + Q 2 + Q 3 .
Now, it remains to verify that Q i 0 as ν 0 for i = 1 , 2 , 3 . Note
Q 1 B 1 ζ ζ + ν P ( ζ + ν κ s ) e F B 1 κ α α y ( s , z ( s ) ) s α 1 d s B 1 ζ ζ + ν W 2 e w κ α α e ( β w ) ( ζ + ν s ) α α L y ( s ) z ( s ) s α 1 d s + B 1 ζ ζ + ν W 2 e w κ α α e ( β w ) ( ζ + ν s ) α α y ( s , 0 ) ) s α 1 d s B 1 W 2 e w κ α α e ( β w ) ν α α 1 g α g + 1 ζ + ν ) g α g + 1 ζ g α g + 1 1 g L y L h ( I , R + ) r + B 1 W 2 e w κ α α e ( β w ) ν α α N 1 α [ ( ζ + ν ) α ζ α ] 0 a s ν 0 ,
and
Q 2 B 1 e F B 1 κ α α e F B 1 ( ζ + ν κ ) α α 0 ζ e κ , α P 1 B 1 ( ζ + ν κ s ) e κ , α P 1 B 1 ( ζ κ s ) · ( L y ( s ) z ( s ) + y ( s , 0 ) ) s α 1 d s B 1 e F B 1 κ α α e F B 1 ( ζ + ν κ ) α α L y L h ( I , R + ) ζ 1 α α r × 0 ζ e κ , α P 1 B 1 ( ζ + ν κ s ) e κ , α P 1 B 1 ( ζ κ s ) g d s 1 g + B 1 e F B 1 κ α α e F B 1 ( ζ + ν κ ) α α N ζ 1 α α 0 ζ e κ , α P 1 B 1 ( ζ + ν κ s ) e κ , α P 1 B 1 ( ζ κ s ) d s 0 a s ν 0 ,
where the uniformly continuity of e κ , α P 1 B 1 · is used. Next
Q 3 B 1 W e β ζ 1 α α 0 ζ e F B 1 ( ζ + ν κ s ) α α e F B 1 ( ζ κ s ) α α ( L y ( s ) z C + N ) s α 1 d s .
Note that lim ν 0 e F B 1 ( ζ + ν κ s ) α α e F B 1 ( ζ κ s ) α α = 0 . So
Q 3 B 1 W e β ζ 1 α α 0 ζ e F B 1 ( ζ + ν κ s ) α α e F B 1 ( ζ κ s ) α α ( L y ( s ) z C + N ) s α 1 d s 0 w h i l e ν 0 ,
from utilizing the dominated convergence theorem of Lebesgue.
Consequently, we receive
( T 2 z ) ( ζ + ν ) ( T 2 z ) ( ζ ) 0 w h i l e ν 0 ,
for all z S r , thus establishing the equicontinuity of T 2 ( S r ) .
Continuing with the aforementioned calculations, we can further deduce that
( T 2 z ) ( ζ ) B 1 W 2 e w κ α α e ( β w ) ζ 1 α α 0 ζ s α 1 L y ( s ) z ( s ) d s + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α 0 ζ y ( s , 0 ) s α 1 d s B 1 W 2 e w κ α α e ( β w ) ζ 1 α α 1 g α g + 1 ζ 1 g α g + 1 1 g L y L h ( I , R + ) r + B 1 W 2 e w κ α α e ( β w ) ζ 1 α α ζ 1 α α N .
Hence T 2 ( S r ) is bounded. From the Arzela-Ascoli theorem, T 2 ( S r ) is relatively compact in C ( I , R n ) . Hence, T 2 : S r C ( I , R n ) is an operator that is both compact and continuous.
By invoking Krasnoselskii’s fixed point theorem, we can ensure the existence of one fixed point z for the operator T on S r . Obviously, z is a solution to the system (4) that satisfies z ( ζ 1 ) = z 1 . Additionally, the boundary condition z ( ζ ) = ψ ( ζ ) , for κ ζ 0 holds true based on Equation (18). Thus, the proof is concluded. □

5. Examples

This section will take a few numerical examples to verify our theoretical results. For simulation purposes, we utilize the infinite-norm.
Example 1.
Let us examine the subsequent linear conformable delay system, which is nonsingular:
B D 0.5 0 z ( ζ ) = F z ( ζ ) + P z ( ζ 0.2 ) , ζ 0 , z ( ζ ) = ψ ( ζ ) = ( 0.2 , 0.1 ) , 0.2 ζ 0 .
where κ = 0.2 and we set
F = 1.2 0.4 0 2 , P = 0.6 0.1 0 0.8 , B = 1.4 0 0 1.4 .
Note
F P = 0.72 0.44 0 1.6 = P F , F B = 1.68 0.56 0 2.8 = B F ,
P B = 0.84 0.14 0 1.12 = B P , B 1 = 0.7143 0 0 0.7143 ,
P 1 B 1 = 9.6 7.44 0 2.16 β e 0.2 β , c h o o s i n g β = 0.9354 .
Clearly, for α = 1 2 , e F B 1 ζ α α = e 2.8 ζ W e 2 w ζ with W = 1 and w = 1 . Next ψ 1 = 0.2 < δ : = 0.25 and
U ( ψ , D α 0 ψ ) = W 2 ( e β κ α α ψ ( κ ) + B 1 κ 0 e w κ α α e ( β w ) ( s ) α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s ) = 0.1616 > 0 .
and β w = 0.0646 < 0 .
At present, those requirements stipulated in Theorem 1 are fulfilled. Hence,
z ( ζ ) U ( ψ , D α 0 ψ ) e ( β w ) 2 ζ = 0.1616 e 0.1292 ζ 0 w h i l e ζ ,
Therefore, the solution to (29) exhibits α-exponential stability.
Example 2.
Let us contemplate the given linear conformable system:
B D 0.5 0 z ( ζ ) = F z ( ζ ) + P z ( ζ 0.2 ) + l ( ζ , z ( ζ 0.2 ) ) , ζ 0 , z ( ζ ) = ψ ( ζ ) = ( 0.1 , 0.2 ) , 0.2 ζ 0 .
where κ = 0.2 and suppose
F = 2.3 0.5 0 2.4 , P = 1.6 0.5 0 1.5 , B = 1.6 0 0 1.6 .
Note
l ( ζ , z ( ζ 0.2 ) ) = 0.2 z 1 ( ζ 0.2 ) s i n ζ 0.15 z 2 ( ζ 0.2 ) s i n ζ .
Note
F P = 3.68 0.4 0 3.6 = P F , F B = 3.68 0.8 0 3.84 = B F ,
P B = 2.56 0.8 0 2.4 = B P , B 1 = 0.625 0 0 0.625 ,
P 1 B 1 = 1.5527 1.0946 0 1.8264 β e 0.4 β , c h o o s i n g β = 0.8211 .
Clearly, for α = 1 2 , e F B 1 ζ α α = e 3 ζ W e 2 w ζ with W = 1 and w = 1.3 . l ( ζ , z ( ζ 0.2 ) ) Q z ( ζ 0.2 ) with Q = 0.2 .
Subsequently ψ 1 = 0.2 < δ : = 0.25 and
U ( ψ , D α 0 ψ ) = W 2 ( e β κ α α ψ ( κ ) + B 1 κ 0 e w κ α α e ( β w ) ( s ) α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s ) = 0.1587 > 0 .
U * ( ψ , D α 0 ψ ) = [ U ( ψ , D α 0 ψ ) + B 1 0 κ W 2 Q e ( w β ) s ψ ( s κ ) s α 1 d s ] e B 1 W 2 Q ζ 1 α α e ( w β ) κ κ = 0.1690 > 0 and B 1 W 2 Q e w κ α α + β w = 0.1445 < 0 .
It is obvious those requirements stipulated in Theorem 3 are fulfilled. Hence,
z ( ζ ) U * ( ψ , D α 0 ψ ) e ( B 1 W 2 Q e w κ α α + β w ) ζ α α = 0.1690 e 0.2890 ζ 0 w h i l e ζ ,
Hence, the solution of (30) exhibits α-exponential stability.
Example 3.
Let us contemplate the given conformable system:
B D 0.5 0 z ( ζ ) = F z ( ζ ) + P z ( ζ 0.2 ) + m ( ζ , z ( ζ ) ) , ζ 0 , z ( ζ ) = ψ ( ζ ) = ( 0.2 , 0.15 , 0.1 ) , 0.2 ζ 0 .
where κ = 0.2 and setting
F = 3.6 0 4 0 6.8 0.8 0 0 2.2 , P = 2.8 0 1 0 2.45 0 0 0 2.45 ,
B = 1.8 0 0 0 1.8 0 0 0 1.8 , m ( ζ , z ( ζ ) ) = 0.2 z 1 ( ζ ) sin ζ 0.1 z 2 ( ζ ) sin ζ 0.2 z 3 ( ζ ) sin ζ .
Note
F P = 10.08 0 13.4 0 16.66 1.96 0 0 5.39 = P F ,
F B = 6.48 0 7.2 0 12.24 1.44 0 0 3.96 = B F ,
P B = 5.04 0 1.8 0 4.41 0 0 0 4.41 = B P , B 1 = 0.5556 0 0 0 0.5556 0 0 0 0.5556 ,
P 1 B 1 = 3.6248 0 3.4163 0 3.7828 0.2463 0 0 2.5999 β e 0.8 β , c h o o s i n g β = 0.3996 .
Clearly, for α = 1 2 , e F B 1 ζ α α = e 7.556 ζ W e 2 w ζ where W = 1 and w = 3 . m ( ζ , z ( ζ ) ) L z ( ζ ) , here L = 0.2 . Together with ψ 1 = 0.2 < δ : = 0.25 ,
U ( ψ , D α 0 ψ ) = W 2 ( e β κ α α ψ ( κ ) + B 1 κ 0 e w κ α α e ( β w ) ( s ) α α B D α 0 ψ ( s ) F ψ ( s ) s α 1 d s ) = 0.1641 > 0 ,
with B 1 W 2 L e w κ α α + β w = 0.9744 < 0 .
It is evident those requirements stipulated in Theorem 2 are fulfilled. Hence,
z ( ζ ) U ( ψ , D α 0 ψ ) e ( B 1 W 2 L e w κ α α + β w ) ζ α α = 0.1641 e 1.9488 ζ 0 w h i l e ζ ,
Thus, the solution of (31) exhibits α-exponential stability.
Example 4.
Set t 1 = 0.9 . Let us contemplate the given conformable differential controlled system:
B D 0.6 0 z ( ζ ) = F z ( ζ ) + P z ( ζ 0.2 ) + y ( ζ , z ( ζ ) ) + J v ( ζ ) , ζ I : = [ 0 , 0.9 ] , z ( ζ ) = ψ ( ζ ) = ( 0.3 , 0.2 ) , 0.2 ζ 0 .
where κ = 0.2 and let
F = 2 0.2 0 1.8 , P = 0.8 0.2 0 0.6 , B = 1 0.2 0 1.2 ,
J = 1 0 0 1 , y ( ζ , z ( ζ ) ) = 0.2 ζ z 1 ( ζ ) 0.1 ζ z 2 ( ζ ) .
Note
F P = 1.6 0.52 0 1.08 = P F , F B = 2 0.16 0 2.16 = B F ,
P B = 0.8 0.08 0 0.72 = B P , B 1 = 1 0.1667 0 0.8333 ,
P 1 B 1 = 1.8221 0.3644 0 1.4577 β e 0.4 β , c h o o s i n g β = 0.7993 .
Clearly, for α = 0.6 , e F B 1 ζ α α = e 3.3 ζ 0.6 W e 2 w ζ 0.6 where W = 1 and w = 1.2 . Together with J = 1 and B 1 = 1 .
Subsequently, we will utilize Equation (26) to estimate U 1 . To achieve this, we first need to obtain M c [ 0 , ζ 1 ] , followed by the calculation of its inverse, M c [ 0 , ζ 1 ] 1 . The delay Grammian matrix takes the form of
M c [ 0 , ζ 1 ] = 0 ζ 1 B 1 P ( ζ 1 κ s ) e F B 1 κ α α J J e F B 1 κ α α P ( ζ 1 κ s ) B 1 s α 1 d s = 0 ζ 1 B 1 e F B 1 ( ζ 1 κ s ) α α e κ , α P 1 B 1 ( ζ 1 κ s ) e F B 1 κ α α J × J e F B 1 κ α α e κ , α ( P 1 B 1 ) ( ζ 1 κ s ) e ( F B 1 ) ( ζ 1 κ s ) α α B 1 s α 1 d s = 0 ζ 1 B 1 e F ( ζ 1 κ s ) α α e κ , α P ( ζ 1 κ s ) e F κ α α J J e F κ α α e κ , α P ( ζ 1 κ s ) e F ( ζ 1 κ s ) α α B 1 s α 1 d s = 0 0.9 B 1 e F ( 0.7 s ) α α e ζ , α P ( 0.7 s ) e F κ α α e F κ α α e κ , α P ( 0.7 s ) e F ( 0.7 s ) α α s α 1 d s = M 1 + M 2 + M 3 .
where F = F B 1 and P = P 1 B 1 ,
M 1 = 0 0.3 B 1 e F ( 0.7 s ) α α I + P ( 0.7 s ) α α + P 2 1 2 ! ( 0.5 s ) α α 2 e F κ α α e F κ α α × I + P ( 0.7 s ) α α + ( P ) 2 1 2 ! ( 0.5 s ) α α 2 e F ( 0.7 s ) α α B 1 s α 1 d s ,
and
M 2 = 0.3 0.6 B 1 e F ( 0.7 s ) α α I + P ( 0.7 s ) α α e F κ α α e F κ α α × I + P ( 0.7 s ) α α e F ( 0.7 s ) α α B 1 s α 1 d s ,
and
M 3 = 0.6 0.9 B 1 e F ( 0.7 s ) α α I 2 e F κ α α e F κ α α e F ( 0.7 s ) α α B 1 s α 1 d s .
Note
M 1 = 19.4212 11.0558 12.2618 7.6964 , M 2 = 6.6098 3.9537 4.2875 2.9241 ,
M 3 = 1.0284 0.7608 0.7767 0.4579 .
Therefore, we obtain
M c [ 0 , ζ 1 ] = 27.0594 15.7703 17.3260 11.0784 , M c 1 [ 0 , ζ 1 ] = 0.4174 0.5942 0.6529 1.0196 .
and
U 1 = M c 1 [ 0 , ζ 1 ] = 0.4254 = 0.6522 .
Moreover, u , v R n ,
y ( ζ , u ) y ( ζ , v ) = max { 0.2 ζ u 1 v 1 , 0.1 ζ u 2 v 2 } 0.2 ζ max { u 1 v 1 , u 2 v 2 } = 0.2 ζ u v .
Now, setting L y ( ζ ) = 0.2 ζ L g ( I , R + ) with g = h = 2 .
Observing that L y L 2 ( I , R + ) = ( 0 0.9 ( 0.2 s ) 2 d s ) 1 2 = 0.0441 ,
b = B 1 W 2 e w κ α α e ( β w ) ( ζ 1 + κ ) α α 1 g α g + 1 ζ 1 g α g + 1 1 g L y L h ( I , R + ) = 0.1030 ,
and
b 1 + W 2 U 1 B 1 J e w κ α α e ( β w ) ζ 1 α α ζ 1 α α = 0.2814 < 1 .
As a result, all the requirements stipulated in Theorem 5 are satisfied, indicating that (32) exhibits relative controllability on the interval [ 0 , 0.9 ] .

6. Conclusions

This article is a generalization of literature [20]. Based on literature [20], this paper introduces a class of nonsingular conformable delay systems with non-singular term B. By giving an estimate of the delayed matrix exponential with non-singular term B, which makes our analysis of systems (1)–(3) more refined. In the first part, we give the sufficient conditions for α -exponential stability of systems (1)–(3) through introducing the definition of α -exponential stability. In the second part, by constructing a Grammian matrix with non-singular terms B, the relative controllability of the linear and nonlinear problems discussed is provided. Finally, we validate our theoretical results with several examples.
Further work will be able to discuss some issues such as periodic solutions and their stability of related systems.

Funding

This research was funded by Guizhou University of Finance and Economics grant number 2019XYB11 and National Natural Science Foundation of China (12261018).

Data Availability Statement

The authors confirm that the data supporting the findings of this study are available within the article.

Acknowledgments

The authors are grateful to the referees for their careful reading of the manuscript and valuable comments. The authors thank the help from the editor too.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Khusainov, D.Y.; Shuklin, G.V. Linear autonomous time-delay system with permutation matrices solving. Stud. Univ. Žilina 2003, 17, 101–108. [Google Scholar]
  2. Diblík, J.; Khusainov, D.Y. Representation of solutions of discrete delayed system x(k + 1)=Ax(k) + Bx(km) + f(k) with commutative matrices. J. Math. Anal. Appl. 2006, 318, 63–76. [Google Scholar] [CrossRef]
  3. Khusainov, D.Y.; Shuklin, G.V. Relative controllability in systems with pure delay. Int. J. Appl. Math. 2005, 2, 210–221. [Google Scholar] [CrossRef]
  4. Medveď, M.; Pospišil, M.; Škripková, L. Stability and the nonexistence of blowing-up solutions of nonlinear delay systems with linear parts defined by permutable matrices. Nonlinear Anal. Theory Methods Appl. 2011, 74, 3903–3911. [Google Scholar] [CrossRef]
  5. Medveď, M.; Pospíšil, M. Sufficient conditions for the asymptotic stability of nonlinear multidelay differential equations with linear parts defined by pairwise permutable matrices. Nonlinear Anal. Theory Methods Appl. 2012, 75, 3348–3363. [Google Scholar] [CrossRef]
  6. Diblík, J.; Fečkan, M.; Pospišil, M. Representation of a solution of the Cauchy problem for an oscillating system with two delays and permutable matrices. Ukrainian Math. J. 2013, 65, 58–69. [Google Scholar] [CrossRef]
  7. Diblík, J.; Khusainov, D.Y.; Užičková, M. Controllability of linear discrete systems with constant coefficients and pure delay. SIAM J. Control Optim. 2008, 47, 1140–1149. [Google Scholar] [CrossRef]
  8. Diblík, J.; Fečkan, M.; Pospišil, M. On the new control functions for linear discrete delay systems. SIAM J. Control Optim. 2014, 52, 1745–1760. [Google Scholar] [CrossRef]
  9. Diblík, J.; Morávkovxax, B. Discrete matrix delayed exponential for two delays and its property. Adv. Differ. Equ. 2013, 2013, 139. [Google Scholar] [CrossRef]
  10. Diblík, J.; Morávkovxax, B. Representation of the solutions of linear discrete systems with constant coefficients and two delays. Abstr. Appl. Anal. 2014, 2014, 320476. [Google Scholar] [CrossRef]
  11. Diblík, J.; Khusainov, D.Y.; Baštinec, J.; Sirenko, A.S. Exponential stability of linear discrete systems with constant coefficients and single delay. Appl. Math. Lett. 2016, 51, 68–73. [Google Scholar] [CrossRef]
  12. You, Z.; Fečkan, M.; Wang, J.; O’Regan, D. Relative controllability of impulsive multi-delay differential systems. Nonlinear Anal. Model. Control 2022, 27, 70–90. [Google Scholar] [CrossRef]
  13. Pospíšil, M. Representation and stability of solutions of systems of functional differential equations with multiple delays. Electron. J. Qual. Theory Differ. Equ. 2012, 54, 1–30. [Google Scholar] [CrossRef]
  14. Pospíšil, M. Representation of solutions of delayed difference equations with linear parts given by pairwise permutable matrices via Z-transform. Appl. Math. Comput. 2017, 294, 180–194. [Google Scholar] [CrossRef]
  15. Wang, J.; Luo, Z.; Fečkan, M. Relative controllability of semilinear delay differential systems with linear parts defined by permutable matrices. Eur. J. Control 2017, 38, 39–46. [Google Scholar] [CrossRef]
  16. Luo, Z.; Wei, W.; Wang, J. Finite time stability of semilinear delay differential equations. Nonlinear Dyn. 2017, 89, 713–722. [Google Scholar] [CrossRef]
  17. Liang, C.; Wang, J.; O’Regan, D. Controllability of nonlinear delay oscillating systems. Electron. J. Qual. Theory Differ. Equ. 2017, 2017, 1–18. [Google Scholar] [CrossRef]
  18. You, Z.; Wang, J. On the exponential stability of nonlinear delay systems with impulses. IMA J. Math. Control Inf. 2018, 35, 773–803. [Google Scholar] [CrossRef]
  19. Wang, J.; Fečkan, M.; Li, M. Stability and Controls Analysis for Delay Systems; Elsevier: Amsterdam, The Netherlands, 2023. [Google Scholar]
  20. Zhou, A.; Wang, J. Relative controllability of conformable delay differential systems with linear parts defined by permutable matrices. Filomat 2023, 37, 2659–2673. [Google Scholar] [CrossRef]
  21. You, Z.; Wang, J.; O’Regan, D. Exponential stability and relative controllability of nonsingular delay systems. Bull. Brazilian Math. Soc. New Ser. 2019, 50, 457–479. [Google Scholar] [CrossRef]
  22. Li, M.; Wang, J. Exploring delayed Mittag-Leffler type matrix functions to study finite time stability of fractional delay differential equations. Appl. Math. Comput. 2018, 324, 254–265. [Google Scholar] [CrossRef]
  23. Cao, X.; Wang, J. Finite-time stability of a class of oscillating systems with two delays. Math. Meth. Appl. Sci. 2018, 41, 4943–4954. [Google Scholar] [CrossRef]
  24. You, Z.; Wang, J. Stability of impulsive delay differential equations. J. Appl. Math. Comput. 2018, 56, 253–268. [Google Scholar] [CrossRef]
  25. Wang, J.; Sathiyaraj, T.; O’Regan, D. Relative controllability of a stochastic system using fractional delayed sine and cosine matrices. Nonlinear Anal. Model. Control 2021, 26, 1031–1051. [Google Scholar] [CrossRef]
  26. Sathiyaraj, T.; Wang, J.; O’Regan, D. Controllability of stochastic nonlinear oscillating delay systems driven by the Rosenblatt distribution. Proc. R. Soc. Edinb. A Math. 2021, 151, 217–239. [Google Scholar] [CrossRef]
  27. **, X.; Wang, J.; Shen, D. Convergence analysis for iterative learning control of impulsive linear discrete delay systems. J. Differ. Equ. Appl. 2021, 27, 739–762. [Google Scholar] [CrossRef]
  28. Wang, X.; Wang, J.; Fečkan, M. Controllability of conformable differential systems. Nonlinear Anal. Model. Control 2020, 25, 658–674. [Google Scholar] [CrossRef]
  29. Malik, M.; Sajid, M.; Kumar, V. Controllability of singular dynamic systems on time scales. Asian J. Control 2022, 24, 2771–2777. [Google Scholar] [CrossRef]
  30. Khalil, R.; Horani, M.A.; Yousef, A.; Sababheh, M. A new definition of fractional derivative. J. Comput. Appl. Math. 2014, 264, 65–70. [Google Scholar] [CrossRef]
  31. **ao, G.; Wang, J. Representation of solutions of linear conformable delay differential equations. Appl. Math. Lett. 2021, 117, 107088. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhou, A. Exponential Stability and Relative Controllability of Nonsingular Conformable Delay Systems. Axioms 2023, 12, 994. https://doi.org/10.3390/axioms12100994

AMA Style

Zhou A. Exponential Stability and Relative Controllability of Nonsingular Conformable Delay Systems. Axioms. 2023; 12(10):994. https://doi.org/10.3390/axioms12100994

Chicago/Turabian Style

Zhou, Airen. 2023. "Exponential Stability and Relative Controllability of Nonsingular Conformable Delay Systems" Axioms 12, no. 10: 994. https://doi.org/10.3390/axioms12100994

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

Article Metrics

Back to TopTop