Collaboration Control of Fractional-Order Multiagent Systems with Sampling Delay

Because of the complexity of the practical environments,many distributedmultiagent systems cannot be illustratedwith the integerorder dynamics and can only be described with the fractional-order dynamics. In this paper, collaboration control problems of continuous-time networked fractional-order multiagent systems via sampled control and sampling delay are investigated. Firstly, the sampled-data control ofmultiagent systemswith fractional-order derivative operator is analyzed in a directedweighted network ignoring sampling delay. Then, the collaborative control of fractional-order multiagent systems with sampled data and sampling delay is studied in a directed and symmetrical network. Many sufficient conditions for reaching consensus with sampled data and sampling delay are obtained. Some numerical simulations are presented to illustrate the utility of our theoretical results.


Introduction
In recent years, consensus problems in distributed networked multiagent systems have attracted increasing attention of more and more researches including control theory, mathematics, biology, physics, computer science, and robotics.The applications of multiagent systems are extensive, ranging from multiple space-craft alignment, heading direction in flocking behavior, distributed computation, and rendezvous of multiple vehicles.Based on certain quantities of interest, collaboration control problems of agent systems have been studied by many researchers and many important results have been achieved in much literature [1][2][3][4][5][6][7][8].
With the development of digital sensors and controllers, in many cases that the system itself is a continuous process, the synthesis of control law can only use the data sampled at the discrete sampling instants.Compare to continuoustime systems with continuous-time controller, continuoustime systems via sampled control have many advantages, such as flexibility, robustness, and low cost.Therefore, sampled control for continuous-time system is more coincident with applications in our real life.Robots, vehicles, airplanes, satellites, and almost all of the modern artificial products are controlled by digital controller where continuous signals are transferred into discrete ones.For consensus problems of continuous-time multiagent systems via sampled control, some interesting results about consensus problem for multiagent system have been reported [9][10][11][12][13][14][15].However, in the many real applications, we always want to find how large the sampling period should be chosen to guarantee that the system runs well.This requires us to look for an upper bound of sampling period.Moreover, sampling delay of the system cannot be ignored and sometimes may play a key role in the stability analysis of the networks.Therefore, we will also study the case when sampling delay exists.
The important results of the above literature focus the consensus problems of multiagent systems with integer-order dynamical equation.In the complex environment, many phenomena cannot be explained by the framework of integerorder dynamics, for example, the synchronized motion of agents in complex environments such as macromolecule fluids and porous media [16][17][18].Under these circumstances, many dynamic characteristics of natural phenomena can only be described in the dynamics of fractional-order (noninteger order) behavior, for example: flocking movement and food searching by means of the individual secretions Mathematical Problems in Engineering and microbial, submarine underwater robots in the bottom of the sea with a large number of microorganisms and viscous substances, unmanned aerial vehicles running in the complex space environment [19,20].Cao et al. [21,22] studied distribution coordination of multiagent systems with fractional-order dynamics firstly and gave the relationship between the number of individuals and the fractional order in the stable multiagent systems.However, to the best of our knowledge, there are few researches done on the coordination control of fractional-order multiagent systems via sampled data.
In this paper, we investigate the consensus of fractionalorder multiagent systems (FOMAS) with sampled-data control.Because few methods are presented to analyze the fractional-order systems with sampling delay, the problems of the fractional-order systems with sampling delay and sampled data will become more difficult.The main innovation of this paper lies in the study on the distributed coordination of FOMAS with sampled data and sampling delay.The rest of the paper is organized as follows.In Section 2, we recall some basic definitions about fractional calculus.In Section 3, some preliminaries about graph theory, fractional-order coordination model of multiagent systems are shown out.A distributed coordination algorithm for FOMAS with sampled data control is studied in Section 4. Section 5 presents the consensus of FOMAS with sampled data and sampling delay.In Section 6, numerical examples are simulated to verify the theoretical analysis.Conclusions are finally drawn in Section 7.

Fractional Calculus
Fractional derivatives provide an excellent instrument to describe the memory and hereditary features of various materials and processes.Fractional calculus also appears in the control of dynamical systems, when the controlled system and the controller are described by a fractionalorder differential equation.This is the main advantage of fractional derivatives in comparison with classical integerorder models, in which such effects are in fact neglected.The advantages of fractional-order derivatives become evident in modeling mechanical and electrical characteristics of real materials, as well as in many fields to describe the rheological properties of rocks.
Fractional operator plays an important role in modern science, which is used as a generalization of integration and differentiation with noninteger order fundamental operator     , where a and t are the limits of the operation and  ∈ .The continuous integrodifferential operator is defined as Three definitions most frequently used for the general fractional operators are the Grünwald-Letnikov (GL) definition, the Riemann-Liouville (RL), and the Caputo definition [16][17][18].The GL definition is given by where [⋅] means the integer part, and (   ) is fractional binomial coefficients.The RL definition is given as for ( − 1 <  < ) and where Γ(⋅) is the Gamma function.
The Caputo definition can be written as The initial conditions for the fractional order differential equations with the Caputo derivatives are in the same form as that for the integer-order differential equations.In this paper, a simple notation  () is used to replace     ().

Problem Statement
Assume that multiagent system consists of  autonomous agents then connected relations between the agents constitute a network topology G. Assume G = {,,} represents a directed weighted graph, in which  = {V 1 , V 2 , . . ., V  } represents a collection of  nodes, and its set of edges is  ⊆  × .The node indexes belong to a finite index set  = {1, 2, . . ., }, with adjacency matrix  = [  ] ∈  × with weighted adjacency elements   ≥ 0. An edge of the weighted diagraph G is denoted by   = (V  , V  ) ∈ .We assume that the adjacency element   > 0 when   ∈ ; otherwise,   = 0.The set of neighbors of a node  is denoted by   = { ∈  :   > 0}.
Let G be a weighted digraph without self-loops, that is, let   = 0, and matrix  = diag{ 1 ,  2 , . . .,   } be the diagonal matrix with the diagonal elements   = ∑  =1   representing the sum of the elements in the th row of matrix .The Laplacian matrix of the weighted digraph G is defined as  = −.For two nodes  and , there is subscript set { 1 ,  2 , . . .  } satisfying   1 > 0,   1  2 > 0, . . .,     > 0, and then there is a directed linked path between node  and node  which is used for the information transmission, also we can say that node  can receive the information from node .If node  can find a path to reach any node of the graph, then node  is globally reachable from every other node in the digraph.In this paper, the directed graph and directed symmetrical graph for fractional-order multiagent systems will be considered.
Given that the dynamics of multiagent systems indicated with fractional derivative in the complex environments, the fractional dynamical equations are defined as where   () ∈  and   () ∈  represent the th agent's state and control input, respectively, and  ()  represents the  ( > 0) order fractional derivative.Assume that the following control protocols are used in multiagent systems: where   represents the (, ) elements of adjacency matrix ,  > 0 is control gain, and   represents the neighbors collection of the th agent.Suppose that for any initial value of the system, the states of autonomous agents meet lim  → ∞ (  () −   ()) = 0, for ,  ∈ , and then multiagent systems asymptotically reach consensus.
Theorem 2. Suppose that multiagent systems are composed of  independent agents, whose connection network topology is directed, and there is a global reachable node.Then fractionalorder multiagent system (7) with sampled data can asymptotically reach consensus, if  < 1 and where   is the eigenvalue of the Laplacian matrix .
Proof.Since the spectral radius  of the matrix Ψ is satisfying  ≤ ||Ψ||, we will require the spectral radius  < 1.By the definition (Ψ) = max{|  |,  = 1, . . ., } where   is the characteristic value of the matrix Ψ, we should calculate the characteristic value   of the matrix Ψ with |  | < 1.
For any matrix , there exists a unitary matrix  satisfying where conjugate matrix   =  −1 , Λ is an upper triangular matrix, and its diagonal elements  1 , . . .,   are the eigenvalues of matrix .From Ψ  =   − ℎ  Λ, we can obtain the characteristic values of matrix Ψ being Suppose that the connection network topology of FOMAS is directed, and there is a global reachable node; then Rank() =  − 1 and  = 0 is a single eigenvalue of Laplacian matrix .
For other characteristic values of Laplacian matrix  with   ̸ = 0( ̸ = 1), the corresponding characteristic values of matrix Ψ are being   =  − ℎ    ,  = 2, . . ., .Let  = Re(  )+ Im(  ) (where  is complex number unit), and then the corresponding characteristic values, for  = 2, . . ., , are From |  | < 1, we can obtain Then, the condition of the fractional-order multiagent system (7) is obtained.The proof is finished.Corollary 3. Suppose that multiagent systems are composed of  independent agents, whose connection network topology is directed and symmetrical, and there is a global reachable node.Then FOMAS (7) can asymptotically reach consensus via sampled data, if  < 1 and where   is the maximum eigenvalue of the Laplacian matrix .
Corollary 4. Suppose multiagent systems are composed of  independent agents, whose connection network topology is directed and symmetrical, and there is a global reachable node.
Then multiagent system (7) with  = 1 can asymptotically reach consensus via sampled data, if where   is the maximum eigenvalue of the Laplacian matrix .

Mathematical Problems in Engineering
Remark 5.The system of Corollary 4 forfractional-order  = 1 becomes the first-order multiagent system.The consensus condition obtained in Corollary 4 is same as that in [9].

Consensus of FOMAS with Sampled Data and Sampling Delay
In the practical application, the sampled-data transferring will result in the communication delays.The sampling delays will affect the control features of the system and sometimes may play a key role in the stability analysis of the network.
In this section, we will study the consensus of multiagent systems with sampled data and sampling delay.Suppose that land the sampling period is ℎ, the sampling delay is .The sampled control protocols are used in multiagent systems as follows: Based on ( 7) and ( 17), the dynamics of FOMAS with sampled data and sampling delay can be rewritten as where and   is a unit matrix with -dimensions.If the norm of matrix Φ is satisfying ‖Φ‖ < 1, the discrete-time FOMAS (18) will asymptotically reach consensus.
Theorem 7. Suppose that multiagent systems are composed of  independent agents, whose connection network topology is directed and symmetrical, and there is a global reachable node.
Then FOMAS (18) with sampled data and sampling delay can asymptotically reach consensus, if  < 1 and where   is the maximum eigenvalue of the Laplacian matrix .

Simulations
Suppose that the system is composed of four fractional-order dynamical agents (Figure 1).The connection weights between individuals are  21 =  12 = 0.7,  42 =  24 = 0.8,  31 =  13 = 0.9, and  14 =  41 = 1.Through the network topology of the system, we can get the adjacency matrix  = ( 0 0.7 0.9 1 0.7 0 0 0.8 0.9 0 0 0 1 0.8 0 0 ) . (30) Suppose that the order of the fractional multiagent dynamics is  = 0.8 and the system control gain is  = 1; then we can obtain the relationship between the sampling period and the upper bound of sampling delays (Figure 2) from the conditions in Theorem 7. In order to make the system meet the condition of reaching consensus, we can set the sampling period according to the sampling delay of the system or decide the stable fields of sampling delays by means of the sampling period.Suppose the order of the fractional multiagent dynamics is  = 0.8 and the control gain  = 1; then we can obtain that the upper bound of sampling period ℎ max = 1.0993 s corresponding delay  = 0.2632 s from Figure 2.   Simulation 1. Assume the sampling delay of multiagent system is  = 0.25 s, we can obtain the upper bound of sampling period is 1.0272 s.In computer simulation, selecting the sampling period ℎ = 0.90 s, the consensus can be asymptotically reached (Figure 3) through fractional-order coordination algorithm.
Simulation 2. Assume the sampled delay of multiagent system is  = 0.49 s, the upper bound of sampling period is 0.5054 in Figure 2. In computer simulation, selecting the sampling period ℎ = 0.50 s, the consensus can be asymptotically reached much more slowly with increasing of time delay (Figure 4).
Simulation 3. Assume the sampling delay of multiagent system is  = 0.50 s, the upper bound of sampling period is 0.4928 in Figure 2. In computer simulation, selecting the sampling period ℎ = 0.50 s, the movement trajectories of the multiagent systems will be asymptotically diverged and the consensus cannot be reached (Figure 5) through fractionalorder coordination algorithm.
Simulation 4. Let time delay continue increasing, suppose the sampling delay is  = 0.80 s, the upper bound of sampling period will be 0.2739 in Figure 2. In computer simulation, selecting the sampling period ℎ = 0.25 s, the movement trajectories of the multiagent systems will be asymptotically converged and the consensus can be reached (Figure 6).Although the sampling period is less than the sampling delay, the consensus can be still achieved under the condition of Theorem 7.

Conclusions
This paper studies distributed coordination of fractionalorder multiagent system with sampled control and sampling delay.By applying the stability theory of discretetime domain, sampled-data control of FOMAS with directed network topology is investigated, and the upper bound of the sampling period is obtained.Based on the Hermite-Biehler Theorem, the collaborative control of fractional-order multiagent systems with sampled data and sampling delay is studied.The relations between sampling delay and sampled period are obtained to ensure the consensus of FOMAS.
Research of the robust stability of FOMAS will be carried out in the following work.

Figure 1 :
Figure 1: Network topology of the multiagent systems.

Figure 2 :
Figure 2: Relationship between the communication delay and the upper bound of sampling period.

Figure 3 :Figure 4 :
Figure 3: Movement trajectories of the multiagent systems with delay 0.25 s and sampling period 0.90 s.

Figure 5 :Figure 6 :
Figure 5: Movement trajectories of the multiagent systems with delay 0.50 s and sampling period 0.50 s.