Next Article in Journal
Multifractal Analysis of the Impact of Fuel Cell Introduction in the Korean Electricity Market
Next Article in Special Issue
Finite-Time Synchronization Criteria for Caputo Fractional-Order Uncertain Memristive Neural Networks with Fuzzy Operators and Transmission Delay Under Communication Feedback
Previous Article in Journal
Fractal Characterization on Three-Dimensional Tortuosity of Fault Tectonic
Previous Article in Special Issue
Mittag-Leffler Synchronization in Finite Time for Uncertain Fractional-Order Multi-Delayed Memristive Neural Networks with Time-Varying Perturbations via Information Feedback
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distributed Consensus Tracking of Incommensurate Heterogeneous Fractional-Order Multi-Agent Systems Based on Vector Lyapunov Function Method

1
School of Integrated Circuits, Wuxi Institute of Technology, Wuxi 214121, China
2
School of Mathematical Sciences, Qufu Normal University, Qufu 273165, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(10), 575; https://doi.org/10.3390/fractalfract8100575
Submission received: 26 July 2024 / Revised: 20 September 2024 / Accepted: 27 September 2024 / Published: 30 September 2024
(This article belongs to the Special Issue Analysis and Modeling of Fractional-Order Dynamical Networks)

Abstract

:
This paper investigates the tracking problem of fractional-order multi-agent systems. Both the order and parameters of the leader are unknown. Firstly, based on the positive system approach, the asymptotically stable criteria for incommensurate linear fractional-order systems are derived. Secondly, the models of incommensurate heterogeneous multi-agent systems are introduced. To cope with incommensurate and heterogeneous situations among followers and the leader, radial basis function neural networks (RBFNNs) and a discontinuous control method are used. Thirdly, the consensus criteria are derived by using the Vector Lyapunov Function method. Finally, a numerical example is presented to illustrate the effectiveness of the proposed theoretical method.

1. Introduction

In recent decades, consensus problems have attracted many people’s attention due to their applications in many fields such as social sciences, physics, control engineering, etc. Consensus of integer-order multi-agent systems (MASs) have been studied in many earlier studies [1,2,3,4,5]. However, there are numerous phenomena that cannot be adequately described by integer-order dynamics, for instance, the synchronized motion of agents in fractional circumstances, like macromolecule fluids and porous media. Consequently, consensus of fractional-order multi-agent systems (FOMASs) were first studied in [6]. In the past decade, FOMASs have attracted much attention, becoming a hot topic in the field of cooperative control [7,8,9,10].
Due to the diversity of each agent, the parameters of every individuality may not be identical, thus called heterogeneous. There are many studies about first-order, second-order or high-order MASs, which have heterogeneous dynamics among their agents [11,12,13,14]. Among these, the complete consensus is difficult to achieve by using static linear feedback control protocol. If the MASs have fractional-order (FO) dynamics, there are also some studies that are similar to those of the integer-order cases [15,16,17,18,19,20]. Note that orders of systems are also parameters in the FOMASs, and the incommensurate order is a distinctive feature of heterogeneous FOMASs. However, there are few studies about FOMASs with different order individualities. In [21], the problem of finite-time consensus tracking for incommensurate nonlinear FOMASs has been investigated, where the communication topology of the MASs is directed and switched. In [22], the problem of practical fixed-time bipartite consensus of a nonlinear incommensurate FOMAS via sliding-mode technique has been discussed, which has a general communication network with a signed directed graph. The above studies promote the finite-time stability theory of incommensurate order FO systems. But the orders of the leader in the above studies are known. In [23], consensus control of incommensurate FOMASs with nonlinear and uncertain dynamics has been studied, in which, the order of the leader is unknown, and the authors have built an analysis frame about the consensus problem of incommensurate systems. However, only bounded consensus has been achieved in the paper.
This manuscript tries to address the tracking problem of heterogeneous FOMASs, in which both the orders and parameters of the leader are unknown. The radial basis function neural networks (RBFNNs) are designed to estimate the errors of the states among different orders, which help to establish the error dynamical equations. Furthermore, the Vector Lyapunov Function method is applied to analyze the stability of the error system. In summary, the contributions of this paper are the following: (1) Based on a positive system approach, the asymptotically stable criteria for incommensurate linear FO systems are derived, and some discussion about the existing literature has been presented. (2) A multi-agent system with incommensurate orders and parameters is used to track the leader with unknown order and parameters, and complete consensus can be achieved under the controllers in this manuscript. (3) An analysis frame of leader-following consensus is established by using the Vector Lyapunov Function method.
The main structure of the remaining part of this paper is organized as follows. In Section 2, some fundamental definitions and lemmas of FO calculus are introduced, and the asymptotic stability of the positive incommensurate FO system is presented. The incommensurate heterogeneous FOMASs and the problem in this paper will be formally stated in Section 3. The main results about sufficient conditions of consensus of the considered FOMASs are presented in Section 4. In Section 5, a simulation example is given to check the theoretical results. Finally, concluding remarks are presented in Section 6.
Notations: Let N be the set of natural integers. Let R be the set of real numbers. R n and R n 1 × n 2 refer to the n-dimensional real vector and n 1 × n 2 real matrices. Let R + n ( R n ) be n-dimensional real vector with non-negative elements. If A = ( a i j ) n 1 × n 2 , B = ( b i j ) n 1 × n 2 , then A B denotes a i j b i j . The superscript “T” denotes matrix transposition. I n denotes the n-dimensional identity matrix. For a vector x R n , x is defined as x = x T x , and x 1 is defined as i = 1 n | x i | . For P R n × n , λ max ( P ) and λ min ( P ) represents the maximum and minimum eigenvalue of P, respectively.

2. Fractional-Order Calculus with Caputo’s Operator

In this subsection, some preliminary knowledge of FO calculus is introduced first. Then, some introductions about the Mittag–Leffler function are introduced, which will be used in later sections. After which, the asymptotic stability of the positive incommensurate FO system is analyzed.
Definition 1
([24]). For a continuous function r ( t ) : [ t 0 , + ) R , the FO integral of order α ( 0 , + ) is defined by
I t α t 0 r ( t ) = 1 Γ ( α ) t 0 t ( t s ) α 1 r ( s ) d s .
Definition 2
([24]). Let function r ( t ) be differentiable on [ t 0 , + ) ; the Caputo FO derivative of order α ( 0 , 1 ) for r ( t ) is defined as
D t α t 0 C r ( t ) = t 0 I t 1 α r ˙ ( t ) = 1 Γ ( 1 α ) t 0 t r ˙ ( s ) ( t s ) α d s .
Lemma 1
([24]). Let n = [ α ] + 1 for α N or n = α for α N . If y ( t ) C n [ t 0 , + ] , then
I t α t 0 D t α t 0 y ( t ) = y ( t ) j = 0 n 1 y ( j ) ( t 0 ) j ! ( t t 0 ) j .
In particular, if 0 < α < 1 and y ( t ) C 1 [ t 0 , + ] , then
I t α t 0 D t α t 0 y ( t ) = y ( t ) y ( t 0 ) .
Lemma 2
([25]). If a vector function y ( t ) R n is derivable on [ t 0 , + ) , we have the following inequality:
D t α t 0 ( y T ( t ) Q y ( t ) ) 2 y T ( t ) Q t 0 D t α y ( t ) ,
where α ( 0 , 1 ) , t t 0 and Q R n × n is a constant positive definite matrix.
In FO calculus, the Mittag–Leffler function E α , β ( z ) = k = 0 + z k Γ ( α k + β ) is significant, let E α ( · ) = E α , 1 ( · ) . The following differentiation formula is an immediate consequence of the definition of the Mittag–Leffler function E α ( t ) .
Lemma 3.
Let g ( t ) = E α ( t ) ; then, g ˙ ( t ) = 1 α E α , α ( t ) . Consequently, let g ^ ( t ) = g ( ν t α ) = E α ( ν t α ) ; then, g ^ ˙ ( t ) = ν t α 1 E α , α ( ν t α ) ; let g ˜ ( t ) = g ( ν ( s t ) α ) = E α ( ν ( s t ) α ) ; then, g ˜ ˙ ( t ) = ν ( s t ) α 1 E α , α ( ν ( s t ) α ) , where ν ( , + ) , s [ t , + ) , α ( 0 , 1 ) are constants and t 0 in g ( t ) , g ^ ( t ) , t s in g ˜ ( t ) .
Proof. 
According to the definition of E α ( t ) and the property of Γ ( · ) that x Γ ( x ) = Γ ( x + 1 ) , one has
g ˙ ( t ) = k = 1 + k t k 1 Γ ( α k + 1 ) = 1 α k = 0 + t k Γ ( α k + α ) = 1 α E α , α ( t ) .
Let u ( t ) = ν t α ; then, g ^ ˙ ( t ) = d g ( u ) d u d u d t = ν t α 1 E α , α ( ν t α ) . Let u ( t ) = ν ( s t ) α ; then, g ˜ ˙ ( t ) = d g ( u ) d u d u d t = ν ( s t ) α 1 E α , α ( ν ( s t ) α ) . This completes our proof. ☐
The Vector Lyapunov Function method will be the main tool in this paper. The analysis of the asymptotic stability of the positive incommensurate FO system is presented in the following Lemma. Considering the autonomous linear incommensurate FO system, in this paper, we assume the initial time as t 0 = 0 :
D t α 0 y ( t ) = R y ( t ) , y ( 0 ) = y 0 R + n ,
where y ( t ) R n , D α 0 y = [ 0 D t α 1 y 1 ( t ) , 0 D t α 2 y 2 ( t ) , , 0 D t α n y n ( t ) ] T R n , α i ( 0 , 1 ) , R = ( r i j ) n × n R n × n .
Lemma 4
([26]). Assume that the A is a Metzler matrix and also a Hurwitz matrix. Then, there exists λ R n such that A λ R n .
Lemma 5.
Assume that R is a Metzler matrix and also a Hurwitz matrix; then, the system (1) is asymptotically stable for any y 0 R + n .
Proof. 
The system can be represented by components as the following:
D t α i 0 y i ( t ) = r i i y i ( t ) + j i r i j y j ( t ) , i = 1 , 2 , , n .
Then, according to the theory of solution for nonhomogeneous linear FO differential equations (one can see Section 3.2 (pp. 54–56) of [27]), it is easy to obtain
y i ( t ) = y i ( 0 ) E α i ( r i i t α i ) + χ ( t ) ,
where
χ ( t ) = 0 t ( t s ) α i 1 E α i , α i ( r i i ( t s ) α i ) j i r i j y j ( s ) d s .
Integrating by parts, we obtain
χ ( t ) = 1 r i i j i r i j y j ( t ) + 1 r i i j i r i j y j ( 0 ) E α ( r i i t α i ) + 1 r i i 0 t E α i ( r i i ( t s ) α i ) j i r i j y ˙ j ( s ) d s .
According to Lemma 4, one can conclude that there must be an η R + n , which is satisfied by R η R n . By some simple calculations and using Lemma 3, we have
y ˙ i ( t ) = y i ( 0 ) r i i t α i 1 E α i , α i ( r i i t α i ) + χ ˙ ( t ) = 0 t ( t s ) α i 1 E α i , α i ( r i i ( t s ) α i ) j i r i j y ˙ j ( s ) d s + t α i 1 E α i , α i ( r i i t α i ) j = 1 n r i j y j ( 0 ) .
When y ( 0 ) = k η , one can obtain j = 1 n r i j y j ( 0 ) < 0 , i = 1 , 2 , , n . Then,
y ˙ i ( t ) 0 t ( t s ) α i 1 E α i , α i ( r i i ( t s ) α i ) j i r i j y ˙ j ( s ) d s .
According to the system, one has
lim t 0 + D t α i 0 y i ( t ) = lim t 0 + j = 1 n r i j y j ( t ) = j = 1 n r i j y j ( 0 ) < 0 ,
i.e.,
lim t 0 + 0 t y ˙ i ( s ) ( t s ) α i d s < 0 ,
one can conclude that
lim t 0 + y ˙ i ( t ) .
Assume that there is a t i * > 0 such that y ˙ i ( t i * ) > 0 and y ˙ i ( t ) 0 when t [ 0 , t i * ) . Based on (3), we have
0 < y ˙ i ( t i * ) 0 t i * ( t i * s ) α i 1 E α i , α i ( r i i ( t i * s ) α i ) j i r i j y ˙ j ( s ) d s 0 ,
which is contradictory; thus, y ˙ i ( t ) 0 , t [ 0 , + ) . Consequently, y i ( t ) is monotonically non-increasing and y i ( t ) 0 .
Thus, one can conclude that lim t + y i ( t ) is existent and finite; assuming that lim t + y i ( t ) = y ¯ i , one has lim t + D t α 0 y ( t ) = R y ¯ with y ¯ = [ y ¯ 1 , y ¯ 2 , , y ¯ n ] R n . Based on the Final Value Theorem of Laplace Transform, we have lim s 0 s L { 0 D t α y ( t ) } = R y ¯ and lim s 0 s L { y ( t ) } = y ¯ , where L { · } denotes the Laplace Transform of a function. Let Y i ( s ) = L { y i ( t ) } and Y ( s ) = [ Y 1 ( s ) , Y 2 ( s ) , , Y n ( s ) ] R n ; then,
R y ¯ = lim s 0 s L { 0 D t α y ( t ) } = lim s 0 s [ s α 1 Y 1 ( s ) s α 1 1 y 1 ( 0 ) , s α 2 Y 2 ( s ) s α 2 1 y 2 ( 0 ) , , s α n Y n ( s ) s α n 1 y n ( 0 ) ] T = lim s 0 diag { s α 1 , s α 2 , , s α n } ( s Y ( s ) y ( 0 ) ) = lim s 0 diag { s α 1 , s α 2 , , s α n } ( y ¯ y ( 0 ) ) ,
which implies that R y ¯ = 0 ; combined with R in nonsingular, one can obtain lim t + y i ( t ) = 0 with any initial condition y 0 R + n . ☐
Remark 1.
Indeed, the above Lemma is a special case of Lemma 3 in [28]. In [28], to prove y i ( t ) is monotonically non-increasing, the authors assume e ( t ) = y ( t ) y ( t + c ) and then derive e ( t ) 0 by analyzing FO system D t α 0 e ( t ) = R e ( t ) . However, according to the definition of FO derivatives,
D t α 0 e ( t ) = 0 D t α y ( t ) 0 D t α y ( t + c ) .
In general,
D t α 0 y ( t + c ) 0 D t + c α y ( t + c ) = R y ( t + c ) .
Consequently, the proof of Lemma 3 in [28] may not appropriate. Based on the above Lemma, we have proved the monotonically non-increasing nature of y i ( t ) via y ˙ i ( t ) 0 directly.

3. Preliminaries and Problem Formulation

3.1. Algebraic Graph Theory

In this subsection, we introduce some fundamental definitions pertaining to algebraic graph theory, which will serve as the basis for subsequent sections.
Let G = { V , E , W } be a graph with N nodes, where V = { v 1 , v 2 , , v N } , E V × V denote the set of nodes and edges, respectively. W = { w i j } R N × N is a weighted adjacency matrix, where w i i = 0 , w i j > 0 if ( v j , v i ) E and w i j = 0 ; otherwise, for all i = 1 , 2 , , N . The corresponding Laplacian matrix L is defined as l i j = w i j if j i and l i i = j i , j = 1 N w i j .

3.2. System Model

In this subsection, the model of FOMASs will be introduced, and the problem formulation will be given. Considering the following heterogeneous FOMAS consists of N agents with different orders, the FO dynamical model of the ith agent is described as follows:
D t α i 0 x i ( t ) = A i x i ( t ) + B i f ( x i ( t ) ) + u i ( t ) , i = 1 , 2 , , N ,
where x i ( t ) R n is the state of the ith agent, α i is the order of the ith agent, f ( x i ( t ) ) = [ f 1 ( x i ( t ) ) , f 2 ( x i ( t ) ) , , f n ( x i ( t ) ) ] T R n is a nonlinear vector value function, u i ( t ) is control input or the communication protocol of the ith agent, which would be designed later, and A i and B i are constant matrices. f ( · ) is assumed to satisfy the following Lipschitz condition throughout this paper:
Assumption 1.
There are constants l i j f 0 , i , j = 1 , 2 , , n such that, for any x , y R n ,
f i ( x ) f i ( y ) j = 1 n l i j f x j y j .
Remark 2.
Let L f = ( l i j f ) n × n . Then, for any diagonal matrices Λ f > 0 , according to the Assumption 1, one has
( x y ) T L f T Λ f L f ( x y ) ( f ( x ) f ( y ) ) T Λ f ( f ( x ) f ( y ) ) .
Furthermore, the above assumption is mild, which is fitted for many systems, such as Lorenz’s systems and Hopfield neural network systems, etc.
The leader with unknown order is described as follows:
D t α 0 0 x 0 ( t ) = A 0 x 0 ( t ) + B 0 f ( x 0 ( t ) ) ,
where α 0 is the unknown order, and A 0 and B 0 are also unknown. In this paper, the u i ( t ) is designed as follows:
u i ( t ) = c j = 1 N l i j x j ( t ) c d i ( x i ( t ) x 0 ( t ) ) c d ^ i sign ( x i ( t ) x 0 ( t ) ) .
Let error signal e i ( t ) = x i ( t ) x 0 ( t ) ; one has
u i ( t ) = c j = 1 N l i j e j ( t ) c d i e i ( t ) c d ^ i sign ( e i ( t ) ) .
Let ϕ i ( x 0 ( t ) ) = 0 D t α 0 x 0 ( t ) 0 D t α i x 0 ( t ) ; then, by some simple derivation, one can obtain the following error dynamic equation:
D t α i 0 e i ( t ) = A i e i ( t ) + B i g ( e i ( t ) ) + w i ( x 0 ( t ) ) + ϕ i ( x 0 ( t ) ) c j = 1 N l i j e j ( t ) c d i e i ( t ) c d ^ i sign ( e i ( t ) ) .
where g ( e i ( t ) ) = f ( x i ( t ) ) f ( x 0 ( t ) ) , w i ( x 0 ( t ) ) = ( A i A 0 ) x 0 ( t ) + ( B i B 0 ) f ( x 0 ( t ) ) .
Assumption 2.
The x 0 ( t ) is bounded.
According to the above assumption and Assumption 1, one can conclude that w i ( x 0 ( t ) ) is also bounded, i.e., there are constants ε i > 0 and T > 0 such that w i ( x 0 ( t ) ) ε i , t T , i = 1 , 2 , , N , where · denotes the Euclidean norm of a vector.
Lemma 6
([29]). If φ ( x ) is a continuous function over a compact set Ω z R n , then there is a RBFNN such that
φ ( z ) = ψ T s ( z ) + δ ( z ) , δ     ε ,
where ε > 0 is a prescribed accuracy, and the ideal weight vector ψ satisfies
ψ = arg min ψ * R m sup z Ω z φ ( z ) ψ * T s ( z )
and δ ( z ) represents the approximation error, ψ * = [ ψ 1 * , ψ 2 * , , ψ m * ] T R m indicates the weight vector, in which m > 1 being the number of the RBFNN nodes, and s ( z ) = [ s 1 ( z ) ,   s 2 ( z ) ,   ,   s m ( z ) ] T R m is the basis function vector with
s i ( z ) = exp ( z v i ) T ( z v i ) q i 2 , i = 1 , 2 , , m ,
where v i = [ v i 1 , v i 2 , , v i n ] T stands for the center vector, and q i is the width of the Gaussian function.
Lemma 7
([30]). The FO continuous-time system D t α 0 x ( t ) = A x ( t ) + B u ( t ) is positive if and only if A is a Metzler matrix, where x ( t ) R n , u ( t ) R n are state and input vectors, A R n × n , B R n × m , D t α 0 x = [ 0 D t α 1 x 1 ( t ) , 0 D t α 2 x 2 ( t ) , , 0 D t α n x n ( t ) ] T R n , α i ( 0 , 1 ) .

4. Main Results

In this section, the analysis of consensus tracking for the systems (5) and (6) under the protocol (7) is presented.
Theorem 1.
If Assumptions 1 and 2 hold, then the consensus tracking of the heterogeneous FOMAS (5) and (6) could be achieved if there are diagonal matrix P > 0 , constants a i > 0 such that
d ^ i ε i λ max ( P ) + ε ¯ i * c λ min ( P ) ;
i P B i P B i Λ f < a i P 0 0 0 ;
where i = P A i + A i T P + L f T Λ f L f + ( 1 + ζ i c ) P 2 c d i I n , Λ f and L f have been mentioned in Remark 1, i = 1 , 2 , , N .
Proof. 
Consider the Vector Lyapunov Function
V ( e ( t ) ) = [ V 1 ( e 1 ( t ) ) , V 2 ( e 2 ( t ) ) , , V N ( e N ( t ) ) ] T R N ,
where V i ( e i ( t ) ) = 1 2 e i T ( t ) P e i ( t ) , i = 1 , 2 , , N . Denoting the fractional derivative on the time t of V i ( e i ( t ) ) with order α i along with (8) as D t α i 0 V i ( e i ( t ) ) , one has
D t α i 0 V i ( e i ( t ) ) e i T ( t ) P D α i e i ( t ) = e i T ( t ) P { A i e i ( t ) + B i g ( e i ( t ) ) + w i ( x 0 ( t ) ) + ϕ i ( x 0 ( t ) ) c j = 1 N l i j e j ( t ) c d i e i ( t ) c d ^ i sign ( e i ( t ) ) } = 1 2 e i T ( t ) ( P A i + A i T P ) e i ( t ) + e i T ( t ) P B i g ( e i ( t ) ) + e i T ( t ) P w i ( x 0 ( t ) ) c d i e i T ( t ) P e i ( t ) + e i T ( t ) P ϕ i ( x 0 ( t ) ) c d ^ i e i T ( t ) P sign ( e i ( t ) ) c e i T ( t ) P j = 1 N l i j e j ( t )
Due to P > 0 , based on the Cauchy inequality(i.e., ( i = 1 n a k b k ) 2 i = 1 n a k 2 i = 1 n b k 2 ) and the definition of the function sign ( · ) , it is easy to obtain
c d ^ i e i T ( t ) P sign ( e i ( t ) ) = c d ^ i j = 1 n p j e i j ( t ) c d ^ i λ min ( P ) e i ( t ) 1 c d ^ i λ min ( P ) e i ( t ) , e i T ( t ) P w i ( x 0 ( t ) ) ε i λ max ( P ) e i ( t ) .
According to Lemma 6, one can conclude that there must be ψ ^ i R n × m , Δ i ( x 0 ) = [ δ 1 i ( x 0 ) , δ 2 i ( x 0 ) , , δ n i ( x 0 ) ] such that
ϕ i ( x 0 ( t ) ) = ψ ^ i s ( x 0 ) + Δ i ( x 0 ) , i = 1 , 2 , , N ,
where δ j i ( x 0 ) ε ^ i j , j = 1 , 2 , , n . Then, we have
e i T ( t ) P ϕ i ( x 0 ( t ) ) = e i T ( t ) P ψ ^ i s ( x 0 ) + Δ i ( x 0 ) 1 2 e i T ( t ) P e i ( t ) + 1 2 ζ i e i T ( t ) P e i ( t ) + ε ¯ i * e i ( t ) ,
where ζ i = ψ ^ i F 2 , ε ¯ i * = max j ε ^ i j , in which, ψ ^ i F denotes the Frobenius norm of matrix ψ ^ i . Moreover, s T ( x 0 ) s ( x 0 ) 1 has been used. According to the Remark 2 and (9), one can obtain
D t α i 0 V i ( e i t ) 1 2 g T ( e i ( t ) ) Λ f g ( e i ( t ) ) c e i T ( t ) P j = 1 N l i j e j ( t ) 1 2 e i T ( t ) { P A i + A i T P + ( 1 + ζ i ) P 2 c d i I n + L f T Λ f L f } e i ( t )
Noting that l i i 0 and l i j 0 when j i , and l i i = j i l i j , then
c e i T ( t ) P j = 1 N l i j e j ( t ) = c l i i e i T ( t ) P e i ( t ) c j i l i j e i T ( t ) P e j ( t ) c l i i e i T ( t ) P e i ( t ) c 2 j i l i j e i T ( t ) P e i ( t ) + e j T ( t ) P e j ( t ) = c 2 l i i e i T ( t ) P e i ( t ) c 2 j i l i j e j T ( t ) P e j ( t ) .
It is easy to obtain
D t α i 0 V i ( t ) ( a i + c l i i ) V i c j i l i j V j ( t ) .
Let A ¯ R N × N be the matrix with A ¯ i i = c ( a i + c l i i ) , A ¯ i j = c l i j for i , j = 1 , 2 , , N . Let
D t α 0 V ( e ( t ) ) = [ 0 D t α 1 V 1 ( e 1 ( t ) ) , 0 D t α 2 V 2 ( e 2 ( t ) ) , , 0 D t α N V N ( e N ( t ) ) ] T R N .
Then, we have
D t α 0 V ( e ( t ) ) A ¯ V ( e ( t ) ) .
There must be C ^ ( t ) 0 such that D t α 0 V ( e ( t ) ) = A ¯ V ( e ( t ) ) C ^ ( t ) . Consider the following FO auxiliary system
D t α 0 W ( t ) = A ¯ W ( t ) W ( 0 ) = V ( e ( 0 ) ) .
Let e ^ ( t ) = W ( t ) V ( e ( t ) ) . One has
D t α 0 e ^ ( t ) = A ¯ e ^ ( t ) + C ^ ( t ) e ^ ( 0 ) = 0 .
From (10) and j = 1 N l i j = 0 , A ¯ is a Metzler matrix. According to Lemma 7, one has t 0 , e ^ ( t ) 0 , i.e., V ( e ( t ) ) W ( t ) . Based on Lemma 5, it is easy to conclude that the solution of (16) is asymptotically stable, which implies that lim t + W ( t ) = 0 . Consequently, one has lim t + V ( e ( t ) ) = 0 . Because V i ( e i ( t ) ) = 1 2 e i T ( t ) P e i ( t ) and P is a positive definite matrix, lim t + e i ( t ) = 0 , i = 1 , 2 , , N , i.e., lim t + e ( t ) = 0 ; then, the tracking consensus of the FOMAS (5) and (6) can be achieved, which completes our proof. ☐

5. Numerical Simulations

Consider a heterogeneous incommensurate FOMAS with four followers: D α i x i ( t ) = A i x i ( t ) + f i ( x i ( t ) ) , where x i R 4 , i = 0 , 1 , 2 , 3 , 4 , and
A 0 = 36 36 0 0 16 28 0 1 0 0 3 0 1 0 0 0 ,   A 1 = 45 45 0 0 1.6 2.8 0 1 0 0 3 0 1 0 0 0
A 2 = 36 36 0 0 16 28 0 1 0 0 2.5 0 1 0 0 1 ,   A 3 = 20 20 0 0 16 25 0 1 0 0 3 0 1 0 0 0
A 4 = 35 35 0 0 15 29 0 1 0 0 0.5 0 1 0 0 0 ,   f i ( x i ( t ) ) = 0 x i 1 x i 3 x i 1 x i 2 g i
in which, g 0 = g 2 = g 4 = 0.5 , g 1 = 1.5 , g 3 = 0.1 ; furthermore, α 0 = 0.99 , α 1 = 0.98 , α 2 = 0.94 , α 3 = 0.95 , α 4 = 0.96 .
Under the above parameters, hyperchaos behaviors will occur of x 0 ( t ) . The phase portraits of the leader are shown in Figure 1 with the initial values [ 10 , 10 , 10 , 10 ] T . Let x i ( 0 ) ( i = 1 , 2 , 3 , 4 ) be random between [ 50 , 50 ] . Then, without control, the states of x i ( t ) ( i = 1 , 2 , 3 , 4 ) are shown in Figure 2. It is obvious that they are chaotic, stable, unstable, etc. Let c = 30 , d i = 25 , d ^ i = 15 , and L = 2 1 0 1 1 2 1 0 0 1 1 0 1 0 0 1 . It is easy to see that, without control, the synchronization error is toward infinity (due to x 3 ( t ) unstable). However, under control, the consensus can be achieved, and the states and errors are shown in Figure 3 and Figure 4, respectively.
Remark 3.
Noting that the order is a significant parameter in the analysis of the dynamics of the fractional-order systems, if the orders are too low, namely, close to zero, the systems usually diverge to infinity. For example, if we set α 0 = 0.3 , then, the states of the leader can be seen in the Figure 5, where the steps of the simulation are selected as 0.0001 . Furthermore, the numerical algorithms of the fractional-order differential equations are in the stage of development; when the order is low, it is difficult to achieve high accuracy. When we set α 0 = 0.6 , α 1 = 0.65 , α 2 = 0.55 , α 3 = 0.65 , α 4 = 0.7 , as one can see in Figure 6. It is easy to see that, without control, the consensus cannot be achieved. However, under control, the consensus can be achieved, as one can see in Figure 7 and Figure 8. Furthermore, the dynamical behaviors are different from the manuscript.

6. Conclusions

In this paper, an incommensurate heterogeneous fractional-order multi-agent system has been studied. As a main theory tool, the Vector Lyapunov Function method for incommensurate fractional-order systems has been discussed. Then, a leader (target) with unknown parameters and orders has been tracked under a discontinuous distributed control protocol. Different from some existing studies, the complete consensus can be achieved with the proposed control strategy instead of a bounded consensus. In general, this manuscript has constructed the theoretical analytical framework for the leader-following consensus of incommensurate fractional-order multi-agent systems. As we all know, there exist many constraints in the communications among agents, such as time delay, stochastic disturbances, bandwidth limitation, etc. In the future, we will extend the framework of this manuscript to fit a more practical networked environment.

Author Contributions

C.H.: writing—original draft preparation, funding acquisition. F.W.: writing—review and editing, visualization, funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the project funded by the High end Training Program for Professional Leaders of Vocational College Teachers in Jiangsu Province (Sugao Peihan [2022] No.11), Jiangsu Province Demonstration Virtual Simulation Training Base Cultivation Project (Sujiao Zhihan [2023] No.30), the Natural Science Foundation of Shandong Province (No: ZR2022QF075), the Project Supported by the Foundation of Key Laboratory of Advanced Process Control for Light Industry (Jiangnan University) Ministry of Education, P.R.China (No: APCLI2403).

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Olfati-Saber, R.; Murray, R. Consensus problems in networks of agents with switching topology and time-delays. IEEE Trans. Autom. Control 2004, 49, 1520–1533. [Google Scholar] [CrossRef]
  2. Cao, Y.; Yu, W.; Ren, W.; Chen, G. An overview of recent progress in the study of distributed multi-agent coordination. IEEE Trans. Ind. Inform. 2012, 9, 427–438. [Google Scholar] [CrossRef]
  3. Oh, K.; Park, M.; Ahn, H. A survey of multi-agent formation control. Automatica 2015, 53, 424–440. [Google Scholar] [CrossRef]
  4. Qin, J.; Ma, Q.; Shi, Y.; Wang, L. Recent advances in consensus of multi-agent systems: A brief survey. IEEE Trans. Ind. Electron. 2016, 64, 4972–4983. [Google Scholar] [CrossRef]
  5. Yu, W.; Zheng, W.; Chen, G.; Ren, W.; Cao, J. Second-order consensus in multi-agent dynamical systems with sampled position data. Automatica 2011, 47, 1496–1503. [Google Scholar] [CrossRef]
  6. Cao, Y.; Li, Y.; Ren, W.; Chen, Y. Distributed coordination of networked fractional-order systems. IEEE Trans. Syst. Man Cybern. Part Cybern. 2009, 40, 362–370. [Google Scholar]
  7. Yu, Z.; Jiang, H.; Hu, C.; Yu, J. Necessary and sufficient conditions for consensus of fractional-order multiagent systems via sampled-data control. IEEE Trans. Cybern. 2017, 47, 1892–1901. [Google Scholar] [CrossRef]
  8. Su, H.; Ye, Y.; Chen, X.; He, H. Necessary and sufficient conditions for consensus in fractional-order multiagent systems via sampled data over directed graph. IEEE Trans. Syst. Man Cybern. Syst. 2021, 51, 2501–2511. [Google Scholar] [CrossRef]
  9. Gong, P.; Lan, W.; Han, Q. Robust adaptive fault-tolerant consensus control for uncertain nonlinear fractional-order multi-agent systems with directed topologies. Automatica 2020, 117, 109011. [Google Scholar] [CrossRef]
  10. Gong, P.; Wang, K. Output feedback consensus control for fractional-order nonlinear multi-agent systems with directed topologies. J. Frankl. Inst. 2020, 357, 1473–1493. [Google Scholar] [CrossRef]
  11. Du, H.; Wen, G.; Wu, D.; Cheng, Y.; Lu¨, J. Distributed fixed-time consensus for nonlinear heterogeneous multi-agent systems. Automatica 2020, 113, 108797. [Google Scholar] [CrossRef]
  12. Li, X.; Yu, Z.; Li, Z.; Wu, N. Group consensus via pinning control for a class of heterogeneous multi-agent systems with input constraints. Inf. Sci. 2021, 542, 247–262. [Google Scholar] [CrossRef]
  13. Zhao, Y.; Xian, C.; Wen, G.; Huang, P.; Ren, W. Design of distributed event-triggered average tracking algorithms for homogeneous and heterogeneous multi-agent systems. IEEE Trans. Autom. Control 2021, 67, 1269–1284. [Google Scholar] [CrossRef]
  14. Han, T.; Zheng, W. Bipartite output consensus for heterogeneous multi-agent systems via output regulation approach. IEEE Trans. Circuits Syst. II Express Briefs 2020, 68, 281–285. [Google Scholar] [CrossRef]
  15. Gong, P. Distributed tracking of heterogeneous nonlinear fractional-order multi-agent systems with an unknown leader. J. Frankl. Inst. 2017, 354, 2226–2244. [Google Scholar] [CrossRef]
  16. Wang, F.; Yang, Y. Quasi-synchronization for fractional-order delayed dynamical networks with heterogeneous nodes. Appl. Math. Comput. 2018, 339, 1–14. [Google Scholar] [CrossRef]
  17. Gao, Z.; Zhang, H.; Zhang, J.; Sun, S. Semi-global leader-following output consensus for heterogeneous fractional-order multi-agent systems with input saturation via observer-based protocol. Neurocomputing 2020, 402, 298–306. [Google Scholar] [CrossRef]
  18. Wang, F.; Zheng, Z.; Yang, Y. Quasi-synchronization of heterogenous fractional-order dynamical networks with time-varying delay via distributed impulsive control. Chaos, Solitons Fractals 2021, 142, 110465. [Google Scholar] [CrossRef]
  19. Fan, H.; Rao, Y.; Shi, K.; Wen, H. Time-varying function matrix projection synchronization of Caputo fractional-order uncertain memristive neural networks with multiple delays via mixed open loop feedback control and impulsive control. Fractal Fract. 2024, 8, 301. [Google Scholar] [CrossRef]
  20. Fan, H.; Chen, X.; Shi, K.; Liang, Y.; Wang, Y.; Wen, H. Mittag-Leffler synchronization in finite time for uncertain fractional-order multi-delayed memristive neural networks with time-varying perturbations via information feedback. Fractal Fract. 2024, 8, 422. [Google Scholar] [CrossRef]
  21. Gong, P.; Han, Q.; Lan, W. Finite-time consensus tracking for incommensurate fractional-order nonlinear multiagent systems with directed switching topologies. IEEE Trans. Cybern. 2020, 52, 65–76. [Google Scholar] [CrossRef] [PubMed]
  22. Gong, P.; Han, Q. Practical fixed-time bipartite consensus of nonlinear incommensurate fractional-order multiagent systems in directed signed networks. SIAM J. Control Optim. 2020, 58, 3322–3341. [Google Scholar] [CrossRef]
  23. Shahvali, M.; Naghibi-Sistani, M.; Modares, H. Distributed consensus control for a network of incommensurate fractional-order systems. IEEE Control Syst. Lett. 2019, 3, 481–486. [Google Scholar] [CrossRef]
  24. Kilbas, A.; Srivastava, H.; Trujillo, J. Theory and Applications of Fractional Differential Equations; Elsevier: New York, NY, USA, 2006. [Google Scholar]
  25. Xu, Q.; Zhuang, S.; Liu, S.; Xiao, J. Decentralized adaptive coupling synchronization of fractional-order complex-variable dynamical networks. Neurocomputing 2016, 186, 119–126. [Google Scholar] [CrossRef]
  26. Rami, M.; Tadeo, F. Controller synthesis for positive linear systems with bounded controls. IEEE Trans. Circuits Syst. II Express Briefs 2007, 54, 151–155. [Google Scholar] [CrossRef]
  27. Balachandran, K. An Introduction to Fractional Differential Equations; Springer: Singapore, 2023. [Google Scholar]
  28. Shen, J.; Lam, J. Stability and performance analysis for positive fractional-order systems with time-varying delays. IEEE Trans. Autom. Control 2016, 61, 2676–2681. [Google Scholar] [CrossRef]
  29. Lu, S.; Wang, X. Adaptive neural network output feedback control of incommensurate fractional-order PMSMs with input saturation via command filtering and state observer. Neural Comput. Appl. 2021, 33, 5631–5644. [Google Scholar] [CrossRef]
  30. Kaczorek, T. Selected Problems of Fractional Systems Theory; Springer: Berlin, Germany, 2011. [Google Scholar]
Figure 1. (ad) Hyperchaos behaviors of leader x 0 ( t ) .
Figure 1. (ad) Hyperchaos behaviors of leader x 0 ( t ) .
Fractalfract 08 00575 g001
Figure 2. (ad) States of x i ( t ) ( i = 1 , 2 , 3 , 4 ) without control.
Figure 2. (ad) States of x i ( t ) ( i = 1 , 2 , 3 , 4 ) without control.
Fractalfract 08 00575 g002
Figure 3. (ad) Consensus behaviors of x i ( t ) ( i = 1 , 2 , 3 , 4 ) with control.
Figure 3. (ad) Consensus behaviors of x i ( t ) ( i = 1 , 2 , 3 , 4 ) with control.
Fractalfract 08 00575 g003
Figure 4. (ad) Consensus errors e i ( t ) ( i = 1 , 2 , 3 , 4 ) with control.
Figure 4. (ad) Consensus errors e i ( t ) ( i = 1 , 2 , 3 , 4 ) with control.
Fractalfract 08 00575 g004
Figure 5. States of leader x 0 ( t ) when α 0 = 0.3 .
Figure 5. States of leader x 0 ( t ) when α 0 = 0.3 .
Fractalfract 08 00575 g005
Figure 6. (ad) States of x i ( t ) ( i = 1 , 2 , 3 , 4 ) without control, α 0 = 0.6 , α 1 = 0.65 , α 2 = 0.55 , α 3 = 0.65 , α 4 = 0.7 .
Figure 6. (ad) States of x i ( t ) ( i = 1 , 2 , 3 , 4 ) without control, α 0 = 0.6 , α 1 = 0.65 , α 2 = 0.55 , α 3 = 0.65 , α 4 = 0.7 .
Fractalfract 08 00575 g006
Figure 7. (ad) Consensus behaviors of x i ( t ) ( i = 1 , 2 , 3 , 4 ) with control, α 0 = 0.6 , α 1 = 0.65 , α 2 = 0.55 , α 3 = 0.65 , α 4 = 0.7 .
Figure 7. (ad) Consensus behaviors of x i ( t ) ( i = 1 , 2 , 3 , 4 ) with control, α 0 = 0.6 , α 1 = 0.65 , α 2 = 0.55 , α 3 = 0.65 , α 4 = 0.7 .
Fractalfract 08 00575 g007
Figure 8. (ad) Consensus errors e i ( t ) ( i = 1 , 2 , 3 , 4 ) with control, α 0 = 0.6 , α 1 = 0.65 , α 2 = 0.55 , α 3 = 0.65 , α 4 = 0.7 .
Figure 8. (ad) Consensus errors e i ( t ) ( i = 1 , 2 , 3 , 4 ) with control, α 0 = 0.6 , α 1 = 0.65 , α 2 = 0.55 , α 3 = 0.65 , α 4 = 0.7 .
Fractalfract 08 00575 g008
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

Huang, C.; Wang, F. Distributed Consensus Tracking of Incommensurate Heterogeneous Fractional-Order Multi-Agent Systems Based on Vector Lyapunov Function Method. Fractal Fract. 2024, 8, 575. https://doi.org/10.3390/fractalfract8100575

AMA Style

Huang C, Wang F. Distributed Consensus Tracking of Incommensurate Heterogeneous Fractional-Order Multi-Agent Systems Based on Vector Lyapunov Function Method. Fractal and Fractional. 2024; 8(10):575. https://doi.org/10.3390/fractalfract8100575

Chicago/Turabian Style

Huang, Conggui, and Fei Wang. 2024. "Distributed Consensus Tracking of Incommensurate Heterogeneous Fractional-Order Multi-Agent Systems Based on Vector Lyapunov Function Method" Fractal and Fractional 8, no. 10: 575. https://doi.org/10.3390/fractalfract8100575

APA Style

Huang, C., & Wang, F. (2024). Distributed Consensus Tracking of Incommensurate Heterogeneous Fractional-Order Multi-Agent Systems Based on Vector Lyapunov Function Method. Fractal and Fractional, 8(10), 575. https://doi.org/10.3390/fractalfract8100575

Article Metrics

Back to TopTop