Next Article in Journal
Analytical Study on Sodium Alginate Based Hybrid Nanofluid Flow through a Shrinking/Stretching Sheet with Radiation, Heat Source and Inclined Lorentz Force Effects
Next Article in Special Issue
Fractional Order Systems and Their Applications
Previous Article in Journal
The Method of Fundamental Solutions for the 3D Laplace Inverse Geometric Problem on an Annular Domain
Previous Article in Special Issue
Chaos Control for a Fractional-Order Jerk System via Time Delay Feedback Controller and Mixed Controller
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotic Stabilization of Delayed Linear Fractional-Order Systems Subject to State and Control Constraints

College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(2), 67; https://doi.org/10.3390/fractalfract6020067
Submission received: 24 December 2021 / Revised: 19 January 2022 / Accepted: 25 January 2022 / Published: 27 January 2022
(This article belongs to the Special Issue Fractional Order Systems and Their Applications)

Abstract

:
Studies have shown that fractional calculus can describe and characterize a practical system satisfactorily. Therefore, the stabilization of fractional-order systems is of great significance. The asymptotic stabilization problem of delayed linear fractional-order systems (DLFS) subject to state and control constraints is studied in this article. Firstly, the existence conditions for feedback controllers of DLFS subject to both state and control constraints are given. Furthermore, a sufficient condition for invariance of polyhedron set is established by using invariant set theory. A new Lyapunov function is constructed on the basis of the constraints, and some sufficient conditions for the asymptotic stability of DLFS are obtained. Then, the feedback controller and the corresponding solution algorithms are given to ensure the asymptotic stability under state and control input constraints. The proposed solution algorithm transforms the asymptotic stabilization problem into a linear/nonlinear programming (LP/NP) problem which is easy to solve from the perspective of computation. Finally, three numerical examples are offered to illustrate the effectiveness of the proposed method.

1. Introduction

Fractional calculus almost appeared at the same time as classic calculus, but it has not been paid more attention to due to its lack of application background and difficult calculation. Fractional calculus has experienced rapid development during the last few decades both in mathematics and applied sciences. It has been recognized as an excellent tool to describe modern complex dynamics [1,2]. From this perspective, some models governing physical phenomena have been reformulated in light of fractional calculus to better reflect their non-local, frequency- and history-dependent properties. With the rapid development of computer technology, fractional calculus is widely used in many fields, such as image processing [3], fluid mechanics [4], and environmental science [5]. Time delay often occurs in different practical systems, and the delayed fractional-order system can better describe these phenomena [6,7,8].
However, time delay will lead to system performance degradation, poor stability and even failure to work [9,10]. In fact, many scholars have extensively studied the stability of delayed fractional-order systems. For instance, the finite-time stability was discussed in [11,12,13,14]. Using the Laplace transform method, the globally asymptotic stability was studied in [15]. The Mittag-Leffler stability was discussed in [16,17]. The asymptotic stability was investigated using the frequency domain method [18], integral inequality method [19], linear matrix inequality (LMI) method [20], and Lyapunov function method [21,22].
On the other hand, for security reasons or physical constraints, the control input, the state, and/or output variables must be bounded in practice. That is to say, the hard constraints of these variables should be considered. At this time, how to design a feedback controller to ensure the stability of the system under constraints is an interesting topic in system theory and synthesis. For constrained integer-order systems, state feedback controllers were designed to study the asymptotic stability of linear discrete-time systems and linear continuous-time systems under constraints in [23,24]. Based on the LMI method, the optimization problem of state feedback controllers for linear discrete-time systems and continuous-time systems was further analyzed in [25,26]. A more effective method to solve the state feedback controller was developed based on the invariance of polyhedron set in [27,28]. As for unconstrained fractional-order systems, in [29], the global stabilization of fractional-order neural network was investigated by using the positive-system-based method. In [30], the stabilization of fractional-order T-S fuzzy systems was discussed by using LMI method. In [31], based on Lyapunov functions, a state feedback controller was designed to study the stabilization of fractional-order nonlinear systems. However, there are few studies concerning the stabilization problems of fractional-order system with constraints, except [32]. In [32], the stabilization problem of fractional-order linear systems with control input constraints was addressed by using the invariance of polyhedron set, but state constraints are not considered there.
Based on the above discussions, it is not only necessary but also more challenging to study the stabilization problem of DLFS under constraints. In this article, the stabilization problem of DLFS with state and control constraints is studied. Our main contributions include: (1) The sufficient conditions that ensure the state constraint set and/or the control constraint set are positive invariant sets (PIS) are established by using the invariant set theory; (2) A new Lyapunov function is constructed on the basis of the constraints, and the asymptotic stability conditions for DLFS are obtained; and (3) A feedback controller and its solution algorithm are proposed to make DLFS under the state and control input constraints asymptotically stable.
The article is organized as follows. Some preliminaries and the problem formulation are given in Section 2. In Section 3, existence conditions for the PIS and the feedback controller are developed. A feedback controller and its solution algorithm are proposed in Section 4 to ensure the DLFS under state and control constraints asymptotically stable. To illustrative the effectiveness of the proposed method, three examples are presented in Section 5. Section 6 contains some conclusions.
Notations: In this article, R + denotes the set of positive real numbers, R n denotes n dimensional real vector space, R n × n denotes n × n real matrices. z stands for a complex number and Re ( z ) stands for the real part of z. ρ i represents the ith element of vector ρ , and Q i represents the ith row of matrix Q, Q i j represents the ith row and jth column element of matrix Q. For ρ R n , ρ 0 ( ρ > 0 ) means ρ i 0 ( ρ i > 0 ) . For a R n , b R n , a b ( a > b ) means a i b i ( a i > b i ) . A > 0 indicates that each entry of the matrix A is nonnegative.

2. Preliminaries and Problem Formulation

2.1. Preliminaries

Consider the DLFS:
0 C D t α x ( t ) = A x ( t ) + A 0 x ( t τ ) + B u ( t ) , t > 0 , x ( t ) = x 0 , τ t 0 ,
where 0 < α 1 ,   τ = constant is the time delay, A R n × n and A 0 R n × n are the system matrices, B R n × m is the input matrix, x ( t ) R n is the state, and u ( t ) R m is the control input.
Definition 1
([12]). The Caputo fractional derivative of x ( t ) is defined as
0 C D t α x ( t ) = 1 Γ ( 1 α ) 0 t ( t s ) α x ( s ) d s ,
where the order 0 < α 1 and Γ ( z ) = 0 e t t z 1 d t , Re ( z ) R + .
Definition 2
([33]). The fractional integral of x ( t ) is defined as
0 I t α x ( t ) = 1 Γ ( α ) 0 t ( t s ) α 1 x ( s ) d s ,
where the order α R + .
By the above Definition, we obtain the solution of system (1) as
x ( t ) = x 0 + 0 I t α ( 0 C D t α x ( t ) ) = x 0 + 1 Γ ( α ) 0 t ( t s ) α 1 A x ( s ) + A 0 x ( s τ ) + B u ( s ) d s .
Definition 3
([34]). A continuous function β ( x ) : [ 0 , + ) [ 0 , + ) is said to be a class-κ function if the function β ( x ) is strictly increasing and β ( 0 ) = 0 .
Definition 4
([29]). A nonempty set P is called the positive invariant sets (PIS) if and only if
x 0 P implies x ( t ) P , for t > 0 ,
where x ( t ) is the trajectory starting with the initial value x 0 .
Lemma 1
([34]). Assume that there exist class-κ functions β i , i = 1 , 2 , 3 and a continuously differentiable function V ( x ( t ) ) such that:
β 1 ( x ( t ) ) V ( x ( t ) ) β 2 ( x ( t ) ) ,
and
0 C D t α V ( x ( t ) ) β 3 ( x ( t ) ) ,
where the order 0 < α < 1 , then system (1) is guaranteed to be stable. If β 3 ( s ) > 0 for s > 0 , then system (1) is guaranteed to be asymptotically stable.
Lemma 2
([35]). The polyhedron set P ( Q , ρ ) = { x ( t ) R n : Q x ( t ) ρ } , Q R q × n , ρ R q , ρ > 0 , and P ( K , ω ) = { x ( t ) R n : K x ( t ) ω } , K R m × n , ω R m , ω > 0 have the relation
P ( Q , ρ ) P ( K , ω )
if and only if there exists L R m × q , L > 0 such that
L Q = K , L ρ ω .

2.2. Problem Formulation

In this article, the following assumptions are needed.
Assumption 1.
The state variables are constrained by the polyhedron set
P ( Q , ρ ) = { x ( t ) R n : Q x ( t ) ρ } ,
where Q R q × n , q > n , and ρ R q , ρ > 0 , and the polyhedron set P ( Q , ρ ) is closed and nonempty.
Assumption 2.
The control input u ( t ) satisfies the following constraints
w 1 u ( t ) w 2 ,
where w 1 > 0 and w 2 > 0 are real vectors.
The asymptotic stabilization problem of DLFS (1) is to find a state feedback controller u ( t ) = K x ( t ) that makes all trajectories starting with the initial value x 0 asymptotically stable and also satisfy the state constraints (5) and the control constraints (6).
If there exists u ( t ) = K x ( t ) for system (1), then, by linear constraints (6), we obtain
P ( K , w 1 , w 2 ) = { x ( t ) R n : w 1 K x ( t ) w 2 } .
Let A ¯ = A + B K ; clearly, u ( t ) = K x ( t ) is the solution of the asymptotic stabilization problem of system (1) if and only if the closed-loop system
0 C D t α x ( t ) = A ¯ x ( t ) + A 0 x ( t τ ) , t > 0 , x ( t ) = x 0 , τ t 0
is asymptotically stable, and the trajectory x ( t ) starting with x 0 satisfy x ( t ) P ( Q , ρ ) and x ( t ) P ( K , w 1 , w 2 ) for any t 0 .

3. Main Results

3.1. Existence Conditions for the Feedback Controller with the Constraints

Theorem 1.
The controller u ( t ) = K x ( t ) , K R m × n is the solution of the asymptotic stabilization problem of system (1) if
(i)
There exists a PIS denoted by M R n for the system (8) satisfies M P ( Q , ρ ) and M P ( K , w 1 , w 2 ) ;
(ii)
There exists a Lyapunov function that makes the system (8) asymptotically stable.
Proof. 
Condition (i) ensures that there exists a PIS
M = { x 0 R n : x 0 M , x ( t ) x ( t ; x 0 ) M }
such that M P ( Q , ρ ) and M P ( K , w 1 , w 2 ) , where x ( t ; x 0 ) is the trajectory starting with x 0 . Then we can obtain Q x ( t ) ρ , w 1 K x ( t ) w 2 . Furthermore, the condition (ii) ensures that system (8) is asymptotically stable. Hence, the controller u ( t ) = K x ( t ) is the solution of the asymptotic stabilization problem of system (1). □
When considering both the state and control constraints, let P ( Q , ρ ) be the PIS of system (8), that is, P ( Q , ρ ) is equal to M, then we obtain:
Corollary 1.
The controller u ( t ) = K x ( t ) , K R m × n is the solution of the asymptotic stabilization problem of system (1) if
(i)
P ( Q , ρ ) is a PIS of the system (8) and P ( Q , ρ ) P ( K , w 1 , w 2 ) ;
(ii)
There exists a Lyapunov function that makes the system (8) asymptotically stable.
Moreover, if considering only the control constraints, let P ( K , w 1 , w 2 ) be the PIS of system (8), that is, P ( K , w 1 , w 2 ) is equal to M , we obtain:
Corollary 2.
The controller u ( t ) = K x ( t ) , K R m × n is the solution of the asymptotic stabilization problem of system (1) if
(i)
P ( K , w 1 , w 2 ) is a PIS of the system (8);
(ii)
There exists a Lyapunov function that makes the system (8) asymptotically stable.

3.2. PIS and Stability Conditions for System (8)

Theorem 2.
If there exists a matrix K R m × n , real matrices F , F 0 R q × q and a scalar ε > 0 such that
Q ( A + B K ) = F Q , Q A 0 = F 0 Q , ( F + F 0 ) ρ ε ρ ,
then the system (8) is asymptotically stable and the polyhedron set P ( Q , ρ ) is a PIS.
Proof. 
Choose the Lyapunov function of the form
V ( x ( t ) ) = max max Q 1 x ( t ) ρ 1 , 0 , , max Q q x ( t ) ρ q , 0 .
If, for 1 i q , it satisfies Q i x ( t ) ρ i < 0 , we obtain V ( x ( t ) ) = 0 . However, due to the fact that Q i x ( t ) cannot always be negative, this case does not exist, so V ( x ( t ) ) > 0 . Hence, there must exist an i, such that 0 < Q i x ( t ) ρ i , that is to say, 0 < V ( x ( t ) ) 1 ; then, according to (9), we obtain
0 C D t α V ( x ( t ) ) = Q i ρ i 0 C D t α x ( t ) = Q i ρ i [ A x ( t ) + A 0 x ( t τ ) + B u ( t ) ] = Q i ρ i [ ( A + B K ) x ( t ) + A 0 x ( t τ ) ] = 1 ρ i [ ( F Q ) i x ( t ) + ( F 0 Q ) i x ( t τ ) ] = 1 ρ i j = 1 n k = 1 q F i k Q k j x j ( t ) + j = 1 n k = 1 q F 0 i k Q k j x j ( t ) = 1 ρ i k = 1 q F i k j = 1 n Q k j x j ( t ) + k = 1 q F 0 i k j = 1 n Q k j x j ( t ) 1 ρ i k = 1 q F i k ρ k + k = 1 q F 0 i k ρ k = 1 ρ i [ ( F + F 0 ) ρ ] i 1 ρ i ( ε ρ i ) ε V ( x ( t ) ) < 0 .
By Lemma 1, we conclude that the system (8) is asymptotically stable.
Assuming that x 0 P ( Q , ρ ) , i.e., Q x 0 ρ and A ¯ = A + B K , according to the Definition of x ( t ) , we obtain
Q x ( t ) = Q x 0 + Q 1 Γ ( α ) 0 t ( t s ) α 1 A ¯ x ( s ) + A 0 x ( s τ ) d s = Q x 0 + 1 Γ ( α ) 0 t ( t s ) α 1 F Q x ( s ) + F 0 Q x ( s τ ) d s ρ + 1 Γ ( α ) 0 t ( t s ) α 1 F Q x ( s ) + F 0 Q x ( s τ ) d s ρ + 1 Γ ( α ) 0 t ( t s ) α 1 ( F ρ + F 0 ρ ) d s ρ ε ρ Γ ( α ) 0 t ( t s ) α 1 d s = ρ ε ρ t α α Γ ( α ) ρ ,
hence the polyhedron set P ( Q , ρ ) is a PIS of system (8). □
Theorem 3.
If there exists a matrix K R m × n , real matrices F , F 0 R q × q , and a scalar ε > 0 such that
K ( A + B K ) = F K , K A 0 = F 0 K , ( F ^ + F 0 ^ ) w ^ ε w ^ ,
where w ^ = w 2 w 1 , F ^ = F + F F F + , and
F i j + = F i j , if i = j , max ( F i j , 0 ) , if i j , F i j = 0 , if i = j , max ( F i j , 0 ) , if i j ,
then the system (8) is asymptotically stable and the polyhedron set P ( K , w 1 , w 2 ) is a PIS.
Proof. 
From (12), we have
F = F + F .
Replace F in (11) with (13), let A ¯ = A + B K ; then, we have
K A ¯ = F K = ( F + F ) K , K A ¯ = F ( K ) = ( F + F ) ( K ) = ( F F + ) K .
Hence,
K K A ¯ = F + F F F + K K .
Let Q = K K , ρ = w ^ = w 2 w 1 , that is, the set P ( K , w 1 , w 2 ) can be reformulated as the form of P ( Q , ρ ) . From (14), there exists a matrix. F ^ = F + F F F + , such that
Q A ¯ = F ^ Q .
In the same way, we have
Q A 0 = F 0 ^ Q .
On the other hand, according to ( F ^ + F 0 ^ ) w ^ ε w ^ , we obtain
( F ^ + F 0 ^ ) ρ ε ρ .
According to (11), we obtain (15)–(17). Then, by Theorem 2, the system (8) is asymptotically stable and the polyhedron set P ( Q , ρ ) is a PIS. Hence, we conclude that the system (8) is asymptotically stable and the polyhedron set P ( K , w 1 , w 2 ) is a PIS. □
Theorem 4.
The polyhedron sets P ( Q , ρ ) and P ( K , w 1 , w 2 ) have the relation
P ( Q , ρ ) P ( K , w 1 , w 2 )
if and only if there exists L R 2 m × 2 q , L > 0 such that
L Q Q = K K , L ρ ρ w 2 w 1 .
Proof. 
According to Lemma 2, P ( Q , ρ ) P ( K , w 2 ) is equivalent to L 1 R m × q , L 1 > 0 such that L 1 Q = K and L 1 ρ w 2 . On the other hand, P ( Q , ρ ) P ( K , w 1 ) is equivalent to L 2 R m × q , L 2 > 0 such that L 2 Q = K and L 2 ρ w 1 .
Hence, P ( Q , ρ ) P ( K , w 1 , w 2 ) if and only if there exists L = L 1 L 2 > 0 , L R 2 m × 2 q such that (18) holds. □
Remark 1.
Theorem 2 proposes a sufficient condition that ensures the state constraint set P ( Q , ρ ) is a PIS and the system (8) is asymptotically stable. On this basis, Theorem 3 proposes a sufficient condition that ensures the control constraint set P ( K , w 1 , w 2 ) is a PIS and the system (8) is asymptotically stable. Theorem 4 gives a sufficient and necessary condition to ensure P ( Q , ρ ) P ( K , w 1 , w 2 ) . When considering both state and control constraints. One can use Theorem 2 and Theorem 4 to find a feedback controller for the asymptotic stabilization problem of system (1) according to Corollary 1. When considering only control constraints. One can use Theorem 3 to find a feedback controller for the asymptotic stabilization problem of system (1) according to Corollary 2.

4. Design Algorithms

In this section, two solution algorithms are designed for the state feedback controllers.
Case 1: Considering both the state and control constraints, by Theorem 2, Theorem 4, and Corollary 1, the solution to matrix inequalities (9) and (18) is the solution to the feedback controller of the asymptotic stabilization problem for system (1). However, for system (1), the rate of convergence to equilibrium is an important index. The largest ε can ensure the fastest convergence rate to the equilibrium. This can be settled by solving the following LP problem with objective function
S ( K , F , F 0 , L , ε ) = ε
and constraints
Q ( A + B K ) = F Q , Q A 0 = F 0 Q , ( F + F 0 ) ρ ε ρ , ε > 0 , L Q Q = K K , L ρ ρ w 2 w 1 .
Case 2: Considering only the control constraints, by Theorem 3 and Corollary 2, the solution to matrix inequality (11) is the solution to the feedback controller of the asymptotic stabilization problem for system (1). This can be found by solving the following NP problem with objective function
S ( K , F , F 0 , ε ) = ε
and constraints
K ( A + B K ) = F K , K A 0 = F 0 K , ( F ^ + F 0 ^ ) w ^ ε w ^ , ε > 0 .
Consider the maximum rate of convergence, from the positive definite function
V ( x ( t ) ) = max max Q 1 x ( t ) ρ 1 , 0 , , max Q q x ( t ) ρ q , 0 ,
we obtain
0 C D t α V ( x ( t ) ) = Q i ρ i 0 C D t α x ( t ) = Q i ρ i [ A x ( t ) + A 0 x ( t τ ) + B u ( t ) ] = Q i ρ i [ ( A + B K ) x ( t ) + A 0 x ( t τ ) ] = 1 ρ i [ ( F Q ) i x ( t ) + ( F 0 Q ) i x ( t τ ) ] = 1 ρ i j = 1 n k = 1 q F i k Q k j x j ( t ) + j = 1 n k = 1 q F 0 i k Q k j x j ( t ) = 1 ρ i k = 1 q F i k j = 1 n Q k j x j ( t ) + k = 1 q F 0 i k j = 1 n Q k j x j ( t ) 1 ρ i k = 1 q F i k ρ k + k = 1 q F 0 i k ρ k = 1 ρ i [ ( F + F 0 ) ρ ] i 1 ρ i ( ε ρ i ) ε V ( x ( t ) ) < 0 .
Therefore, maximizing ε is to maximize the rate of convergence.
Remark 2.
When the system parameters A , A 0 , and B and the constraint parameters Q , ρ , and w ^ are fixed, from (23) we can see that the parameter ε is closely related to the rate of convergence. If 0 < ε < 1 , the largest ε can ensure the fastest convergence rate to the equilibrium.

5. Numerical Examples

Example 1.
Consider the delayed fractional-order electrical circuit [36] shown in Figure 1 with the given resistance R = 1 Ω , inductance L = 0.4167 mH , capacitance C = 1.67 mF , delay element and source voltage u.
Let v ( t ) represent the voltage across the capacitor and i ( t ) represent the current passing through the inductor, using Kirchhoff’s laws, we can write the circuit equations
d α v ( t ) d t α d α i ( t ) d t α = A v ( t ) i ( t ) + A 0 v ( t τ ) i ( t τ ) + B u ( t )
where A = 0 1 C 1 L R L , A 0 = 0 1 R C L 0 , B = 0 1 L .
Let x 1 ( t ) = v ( t ) , x 2 ( t ) = i ( t ) . Then, the above circuit equations can be written in the form of the DLFS (1), with α = 0.7 , τ = 0.1 , A = 0 0.6 2.4 2.4 , A 0 = 0 1 4 0 , and B = 0 2.4 .
Suppose the control input P ( K , w ̲ , w ¯ ) satisfies the constraint
100 K x ( t ) 50 .
Without the control input, system (1) is unstable; for the initial condition ( 4 , 1 ) T , the time response of the state is shown in Figure 2.
The asymptotic stabilization problem of DLFS (1) is to find a controller u ( t ) = K x ( t ) that makes all trajectories starting with the initial condition x 0 asymptotically stable and meanwhile satisfy the control constraints (24).
According to Case 2, solving the NP problem with the objective function
S ( K , F , F 0 , ε ) = ε
and under the constraints
K 0 0.6 2.4 2.4 + 0 2.4 K = F K K 0 1 4 0 = F 0 K F + F F F + + F 0 + F 0 F 0 F 0 + 50 100 ε 50 100 ε > 0 .
By calculation, we obtain ε = 0.8 ,   K = 2 1 ,   F = 1.2 , and F 0 = 2 .
It indicates that there exists a u ( t ) = 2 1 x ( t ) such that all trajectories originating from the initial condition are asymptotically stable to the origin, and meanwhile the corresponding trajectory satisfies the control constraints (24).
For the initial condition ( 4 , 1 ) T , the time response of the state is shown in Figure 3.
Additionally, u ( t ) satisfies the control constraints (24), which can be seen from Figure 4.
So, u ( t ) = 2 1 x ( t ) is a solution to the asymptotic stabilization problem of DLFS (1).
Example 2.
Consider another delayed fractional-order electrical circuits, which can be written in the form of the DLFS (1), with α = 0.7 , τ = 0.5 , A = 3.25 10 1 7.1 , A 0 = 1 2.2 0 2.1 , and B = 1 0.8 .
The state constraints P ( Q , ρ ) are
2 x 2 ( t ) 0.8 , x 1 ( t ) + 2 x 2 ( t ) 2.4 , 0.5 x 1 ( t ) + 4 x 2 ( t ) 1.8 ,
where Q = 0 2 1 2 0.5 4 R 3 × 2 and ρ = 0.8 2.4 1.8 R 3 , ρ > 0 . It can be easily seen that the polyhedron set P ( Q , ρ ) is closed and nonempty.
Suppose the control input P ( K , w ̲ , w ¯ ) satisfies the constraint
6 K x ( t ) 9 .
Without the control input, system (1) is unstable, the time response of the state with the initial condition ( 6.8 , 0.4 ) T is shown in Figure 5.
The asymptotic stabilization problem of DLFS (1) is to find a controller u ( t ) = K x ( t ) that makes all trajectories starting from the initial value x 0 asymptotically stable and meanwhile satisfy the state constraints (25) and the control constraints (26).
According to Case 1, solving the LP problem with the objective function
S ( K , F , F 0 , L , ε ) = ε
and under the constraints
0 2 1 2 0.5 4 3.25 10 1 7.1 + 1 0.8 K = F 0 2 1 2 0.5 4 , 0 2 1 2 0.5 4 1 2.2 0 2.1 = F 0 0 2 1 2 0.5 4 , F + F 0 0.8 2.4 1.8 ε 0.8 2.4 1.8 , L 0 1 0.5 0 1 0.5 2 2 4 2 2 4 T = K K , L 0.8 2.4 1.8 0.8 2.4 1.8 T 9 6 , ε > 0 .
By computation, we obtain
ε = 0.9 ,
K = 1.25 10 ,
F = 0.9 0 0 1.1 2 0 0 0.7 3 ˙ 0.5 3 ˙ ,
F 0 = 2.1 0 0 0 1 0 0 0.55 2.1 ,
and
L = 0 0 0 0 0 2.5 0 0 0 3.75 1.25 0 .
It suggests that there exists u ( t ) = 1.25 x 1 ( t ) 10 x 2 ( t ) such that all trajectories starting from the initial condition are asymptotically stable to the origin, and meanwhile the corresponding state trajectory satisfies the state constraints (25) and the control constraints (26).
For the initial condition ( 6.8 , 0.4 ) T , the time response of the state is shown in Figure 6.
For the initial condition ( 6.8 , 0.4 ) T , the corresponding phase trajectory which satisfies the state constraints (25) and the control constraints (26), is shown in Figure 7.
Furthermore, u ( t ) satisfies the control constraints (26), which can be seen in Figure 8.
Hence, u ( t ) = 1.25 x 1 ( t ) 10 x 2 ( t ) is a solution to the asymptotic stabilization problem of DLFS (1).
Clearly, from Example 2, it can be seen that the control input is not saturated because the state constraint set is a PIS. Next, we give an example without state constraints to verify the saturation of the input control.
Example 3.
Consider the DLFS (1), with α = 0.8 , τ = 0.8 , A = 1 3 1 5 , A 0 = 1 0 0 2 , and B = 1 1 1 0.5 .
Suppose the control input P ( K , w 1 , w 2 ) satisfies the constraint
1 2 K x ( t ) 2 4 .
Without the control input, system (1) is unstable; for the initial condition ( 1.84 , 0.64 ) T , the time response of the state is shown in Figure 9.
The asymptotic stabilization problem of DLFS (1) is to find a controller u ( t ) = K x ( t ) that makes all trajectories starting with the initial condition x 0 asymptotically stable and meanwhile satisfy the control constraints (27).
According to Case 2, solving the NP problem with the objective function
S ( K , F , F 0 , ε ) = ε
and under the constraints
K 1 3 1 5 + 1 1 1 0.5 K = F K K 1 0 0 2 = F 0 K F + F F F + + F 0 + F 0 F 0 F 0 + 2 4 1 2 ε 2 4 1 2 ε > 0 .
By calculation, we obtain
ε = 0.8908 ,
K = 2.051 7.444 1.255 2.629 ,
F = 4.11 0.43 0.33 1.96
and
F 0 = 1.634 1.0362 0.2239 1.3660 .
It indicates that there exists a u ( t ) = 2.051 7.444 1.255 2.629 x ( t ) such that all trajectories originating from the initial condition are asymptotically stable to the origin, and meanwhile the corresponding trajectory satisfies the control constraints (27).
For the initial condition ( 1.84 , 0.64 ) T , the time response of the state is shown in Figure 10.
For the initial condition ( 1.84 , 0.64 ) T , the corresponding phase trajectory is shown in Figure 11. It can be seen that the control constraints (27) is satisfied.
In addition, u ( t ) satisfies the control constraints (27), which can be seen from Figure 12.
So, u ( t ) = 2.051 7.444 1.255 2.629 x ( t ) is a solution to the asymptotic stabilization problem of DLFS (1).

6. Conclusions

The asymptotic stabilization problem of DLFS subject to state and control constraints is studied. Based on the invariant set theory and fractional-order Lyapunov stability theory, the feedback controller and the corresponding solution algorithms that ensure the asymptotic stability of the DLFS under constraints are given. Numerical examples show that the proposed method is effective. Stability of delayed fractional-order neural network systems with constraints are a very interesting topic. In the future, we will study the stabilization of delayed fractional-order neural network systems subject to constraints by using the invariant set method.

Author Contributions

Conceptualization, X.S. and Z.W.; methodology, X.S. and Z.W.; software, X.S.; validation, X.S., Z.W., Z.S. and Z.Z.; investigation, Z.W.; writing—original draft preparation, X.S.; writing—review and editing, Z.W., Z.S. and Z.Z.; supervision, Z.W.; funding acquisition, Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China under Grants 61973199, 62003794, 62173214, the Shandong Provincial Natural Science Foundation ZR2020QF050, ZR2021MF003 and the Taishan Scholar Project of Shandong Province of China.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable to this article as no datasets were generated or analysed during the current study.

Acknowledgments

We would like to express our great appreciation to the editors and reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DLFSdelayed linear fractional-order systems
PISpositive invariant sets
LPlinear programming
NPnonlinear programming
LMIlinear matrix inequality

References

  1. Chen, L.; Chen, Y.; Lopes, A.M.; Kong, H.; Wu, R. State of charge estimation of lithium-ion batteries based on fuzzy fractional-order unscented kalman filter. Fractal Fract. 2021, 5, 91. [Google Scholar] [CrossRef]
  2. Tian, Y.; Xia, Q.; Chai, Y.; Chen, L.; Lopes, A.M.; Chen, Y. Guaranteed cost leaderless consensus protocol design for fractional-order uncertain multi-agent systems with state and input delays. Fractal Fract. 2021, 5, 141. [Google Scholar] [CrossRef]
  3. Gonzalez-Acuna, R.G.; Davila, A.; Gutierrez-Vega, J.C. Optical flow of non-integer order in particle image velocimetry techniques. Signal Process. 2019, 155, 317–322. [Google Scholar] [CrossRef]
  4. Lashkarian, E.; Hejazi, S.R.; Habibi, N.; Motamednezhad, A. Symmetry properties, conservation laws, reduction and numerical approximations of time-fractional cylindrical-burgers equation. Commun. Nonlinear Sci. Numer. Simul. 2019, 67, 176–191. [Google Scholar] [CrossRef]
  5. Yin, C.; Dadras, S.; Cheng, Y.H.; Huang, X.G.; Cao, J.W.; Malek, H. Multidimensional fractional-order newton-based extremum seeking for online light-energy saving technique of lighting system. IEEE Trans. Ind. Electron. 2019, 67, 8576–8586. [Google Scholar] [CrossRef]
  6. Mohsenipour, R.; Jegarkandi, M.F. Robust stability analysis of fractional-order interval systems with multiple time delays. Int. J. Robust Nonlinear Control 2019, 29, 1823–1839. [Google Scholar] [CrossRef]
  7. Sakthivel, R.; Mohanapriya, S.; Ahn, C.K.; Karimi, H.R. Output tracking control for fractional-order positive switched systems with input time delay. IEEE Trans. Circuits Syst. II-Express Briefs 2019, 66, 1013–1017. [Google Scholar] [CrossRef]
  8. Jia, J.; Huang, X.; Li, Y.; Cao, J. Global stabilization of fractional-order memristor-based neural networks with time delay. IEEE Trans. Neural Netw. Learn. Syst. 2020, 31, 997–1009. [Google Scholar] [CrossRef]
  9. Wang, Z.; Wang, X.; Xia, J.; Shen, H.; Meng, B. Adaptive sliding mode output tracking control based-FODOB for a class of uncertain fractional-order nonlinear time-delayed systems. Sci. China-Technol. Sci. 2020, 63, 1854–1862. [Google Scholar] [CrossRef]
  10. Liu, H.; Pan, Y.; Cao, J.; Zhou, Y.; Wang, H. Positivity and stability analysis for fractional-order delayed systems: A T-S fuzzy model approach. IEEE Trans. Fuzzy Syst. 2020, 29, 927–939. [Google Scholar] [CrossRef]
  11. Du, F.; Lu, J. Finite-time stability of neutral fractional order time delay systems with Lipschitz nonlinearities. Appl. Math. Comput. 2020, 375, 125079. [Google Scholar] [CrossRef]
  12. Naifar, O.; Nagy, A.M.; Makhlouf, A.B.; Kharrat, M. Finite-time stability of linear fractional-order time-delay systems. Int. J. Robust Nonlinear Control 2018, 29, 180–187. [Google Scholar] [CrossRef] [Green Version]
  13. Liu, P.; Zeng, Z.; Wang, J. Asymptotic and finite-time cluster synchronization of coupled fractional-order neural networks with time delay. IEEE Trans. Neural Netw. Learn. Syst. 2020, 31, 4956–4967. [Google Scholar] [CrossRef] [PubMed]
  14. Thanh, N.T.; Phat, V.N.; Niamsup, P. New finite-time stability analysis of singular fractional differential equations with time-varying delay. Fract. Calc. Appl. Anal. 2020, 23, 504–519. [Google Scholar] [CrossRef]
  15. Deng, W.; Li, C.; Lu, J. Stability analysis of linear fractional differential system with multiple time delays. Nonlinear Dyn. 2007, 48, 409–416. [Google Scholar] [CrossRef]
  16. You, X.; Song, Q.; Zhao, Z. Global Mittag-Leffler stability and synchronization of discrete-time fractional-order complex-valued neural networks with time delays. Neural Netw. 2020, 122, 382–394. [Google Scholar] [CrossRef] [PubMed]
  17. Hu, T.; He, Z.; Zhang, X.; Zhong, S. Finite-time stability for fractional-order complex-valued neural networks with time delay. Appl. Math. Comput. 2020, 365, 124715. [Google Scholar] [CrossRef]
  18. Bonnet, C.; Partington, J.R. Stabilization of some fractional delay systems of neutral type. Automatica 2007, 43, 2047–2053. [Google Scholar] [CrossRef]
  19. He, B.; Zhou, H.; Kou, C.; Chen, Y.Q. New integral inequalities and asymptotic stability of fractional-order systems with unbounded time delay. Nonlinear Dyn. 2018, 94, 1523–1534. [Google Scholar] [CrossRef]
  20. Chen, L.; Wu, R.; Cheng, Y.; Chen, Y.Q. Delay-dependent and order-dependent stability and stabilization of fractional-order linear systems with time-varying delay. IEEE Trans. Circuits Syst. II-Express Briefs 2020, 67, 1064–1068. [Google Scholar] [CrossRef]
  21. Tavazoei, M.S.; Badri, V. Stability analysis of fractional order time-delay systems: Constructing new lyapunov functions from those of integer order counterparts. IET Control Theory Appl. 2019, 13, 2476–2481. [Google Scholar]
  22. Nie, X.; Liu, P.; Liang, J.; Cao, J. Exact coexistence and locally asymptotic stability of multiple equilibria for fractional-order delayed Hopfield neural networks with Gaussian activation function. Neural Netw. 2021, 142, 690–700. [Google Scholar] [CrossRef] [PubMed]
  23. Hennet, J.C.; Beziat, J.P. A class of invariant regulators for the discrete-time linear constrained regulation problem. Automatica 1991, 27, 549–554. [Google Scholar] [CrossRef]
  24. Vassilaki, M.; Bitsoris, G. Constrained regulation of linear continuous-time dynamical systems. Syst. Control Lett. 1989, 13, 247–252. [Google Scholar] [CrossRef]
  25. Bitsoris, G.; Olaru, S. Further results on the linear constrained regulation problem. In Proceedings of the 21st Mediterranean Conference on Control and Automation, Crete, Greece, 25–28 June 2013; pp. 824–830. [Google Scholar]
  26. Bitsoris, G.; Olaru, S.; Vassilaki, M. On the linear constrained regulation problem for continuous-time systems. In Proceedings of the 19th World Congress International Federation of Automatic Control, Cape Town, South Africa, 24–29 August 2014; pp. 4004–4009. [Google Scholar]
  27. Castelan, E.B. On invariant polyhedra of continuous-time linear systems. IEEE Trans. Autom. Control 1991, 38, 1680–1685. [Google Scholar] [CrossRef]
  28. Blanchini, F. Set invariance in control. Automatica 1999, 35, 1747–1767. [Google Scholar] [CrossRef]
  29. Jia, J.; Wang, F.; Zeng, Z. Global stabilization of fractional-order memristor-based neural networks with incommensurate orders and multiple time-varying delays: A positive-system-based approach. Nonlinear Dyn. 2021, 104, 2303–2329. [Google Scholar] [CrossRef]
  30. Zhang, X.; Zhao, Z. Normalization and stabilization for rectangular singular fractional order T-S fuzzy systems. Fuzzy Sets Syst. 2019, 381, 140–153. [Google Scholar] [CrossRef]
  31. Zhang, Z.; Zhang, J. Asymptotic stabilization of general nonlinear fractional-order systems with multiple time delays. Nonlinear Dyn. 2020, 102, 605–619. [Google Scholar]
  32. Lim, Y.H.; Ahn, H.S. On the positive invariance of polyhedral sets in fractional-order linear systems. Automatica 2013, 49, 3690–3694. [Google Scholar] [CrossRef]
  33. Yan, L.; Chen, Y.Q.; Podlubny, I. Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized Mittag-Leffler stability. Comput. Math. Appl. 2010, 59, 1810–1821. [Google Scholar]
  34. Wen, Y.; Zhou, X.F.; Zhang, Z.; Liu, S. Lyapunov method for nonlinear fractional differential systems with delay. Nonlinear Dyn. 2015, 82, 1015–1025. [Google Scholar] [CrossRef]
  35. Athanasopoulos, N.; Bitsoris, G.; Vassilaki, M. Stabilization of bilinear continuous-time systems. In Proceedings of the 18th Mediterranean Conference on Control and Automation, MED’10, Marrakech, Morocco, 23–25 June 2010; pp. 442–447. [Google Scholar]
  36. Kaczorek, T.; Rogowski, K. Fractional Linear Systems and Electrical Circuits; Springer: New York, NY, USA, 2014; p. 65. [Google Scholar]
Figure 1. The delayed fractional-order electrical circuit in [36].
Figure 1. The delayed fractional-order electrical circuit in [36].
Fractalfract 06 00067 g001
Figure 2. The time response of the system state without feedback control.
Figure 2. The time response of the system state without feedback control.
Fractalfract 06 00067 g002
Figure 3. The time response of the system state with feedback control.
Figure 3. The time response of the system state with feedback control.
Fractalfract 06 00067 g003
Figure 4. The control input u ( t ) .
Figure 4. The control input u ( t ) .
Fractalfract 06 00067 g004
Figure 5. The time response of the system state without feedback control.
Figure 5. The time response of the system state without feedback control.
Fractalfract 06 00067 g005
Figure 6. The time response of the system state with feedback control.
Figure 6. The time response of the system state with feedback control.
Fractalfract 06 00067 g006
Figure 7. The phase trajectory with the initial condition ( 6.8 , 0.4 ) T .
Figure 7. The phase trajectory with the initial condition ( 6.8 , 0.4 ) T .
Fractalfract 06 00067 g007
Figure 8. The control input u ( t ) .
Figure 8. The control input u ( t ) .
Fractalfract 06 00067 g008
Figure 9. The time response of the system state without feedback control.
Figure 9. The time response of the system state without feedback control.
Fractalfract 06 00067 g009
Figure 10. The time response of the system state with feedback control.
Figure 10. The time response of the system state with feedback control.
Fractalfract 06 00067 g010
Figure 11. The phase trajectory with the initial condition ( 1.84 , 0.64 ) T .
Figure 11. The phase trajectory with the initial condition ( 1.84 , 0.64 ) T .
Fractalfract 06 00067 g011
Figure 12. The control input u ( t ) .
Figure 12. The control input u ( t ) .
Fractalfract 06 00067 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Si, X.; Wang, Z.; Song, Z.; Zhang, Z. Asymptotic Stabilization of Delayed Linear Fractional-Order Systems Subject to State and Control Constraints. Fractal Fract. 2022, 6, 67. https://doi.org/10.3390/fractalfract6020067

AMA Style

Si X, Wang Z, Song Z, Zhang Z. Asymptotic Stabilization of Delayed Linear Fractional-Order Systems Subject to State and Control Constraints. Fractal and Fractional. 2022; 6(2):67. https://doi.org/10.3390/fractalfract6020067

Chicago/Turabian Style

Si, Xindong, Zhen Wang, Zhibao Song, and Ziye Zhang. 2022. "Asymptotic Stabilization of Delayed Linear Fractional-Order Systems Subject to State and Control Constraints" Fractal and Fractional 6, no. 2: 67. https://doi.org/10.3390/fractalfract6020067

APA Style

Si, X., Wang, Z., Song, Z., & Zhang, Z. (2022). Asymptotic Stabilization of Delayed Linear Fractional-Order Systems Subject to State and Control Constraints. Fractal and Fractional, 6(2), 67. https://doi.org/10.3390/fractalfract6020067

Article Metrics

Back to TopTop