Abstract

In this work, the consensus problem of fractional-order multiagent systems with the general linear model of fixed topology is studied. Both distributed -type and -type fractional-order iterative learning control (FOILC) algorithms are proposed. Here, a virtual leader is introduced to generate the desired trajectory, fixed communication topology is considered, and only a subset of followers can access the desired trajectory. The convergence conditions are proved using graph theory, fractional calculus, and λ norm theory. The theoretical analysis shows that the output of each agent completely tracks the expected trajectory in a limited time as the iteration number increases for both -type and -type FOILC algorithms. Extensive numerical simulations are given to demonstrate the feasibility and effectiveness.

1. Introduction

In recent years, the research on the coordinating control of multiagent systems (MASs) has become a major task in the field of control because of its application in various areas, such as formation control of unmanned aerial vehicles [1], air traffic control [2], and formation of multiple robots [3]. The consensus is an important and fundamental problem for the coordination control of MASs. One object of the consensus problem is that all agents are expected to reach an agreement in some amounts (such as speed, position, phases, and attitudes) by communicating with their local neighbors [46].

The consensus problem, as an interesting topic for MASs, has achieved many results [510]. However, most studies have focused on the consensus of integer-order MASs (IOMASs). In fact, many systems cannot be well described for some complex problems by integer-order models, such as viscoelastic systems, underwater vehicle, and electrical machines [1113]. Recently, fractional-order calculus is found to be suitable for describing some complex systems, such as macromolecule fluids, porous media, and viscoelastic systems, and have excellent memory and hereditary [14]. Therefore, it is meaningful to study the consensus problem for fractional-order MASs (FOMASs). In particular, integer-order systems are special cases of fractional-order systems, and many researchers have studied the convergence problem and focused on the coordination control of FOMASs. The necessary and sufficient conditions of FOMASs for containment control were presented in [15]. For FOMASs, several consensus problems of FOMASs were studied in [1619]. Du et al. studied the problem of robust consensus for second-order MASs and proposed a class of novel continuous nonsmooth consensus algorithms [20]. Zhu et al. studied the leader-following consensus problem of FOMASs with general linear and nonlinear models with input time delay and designed a control gain matrix and sufficient condition [21]. In [22], group multiple-lag consensus of nonlinear FOMAS with leader-following was investigated, and two kinds of lag consensus were considered.

Notably, the aforementioned studies consider the asymptotic convergence problem of FOMAS, that is, the tracking errors of each agent decreases gradually as time increases. For MAS with repetitive tasks, such as multimechanical arms on industrial production lines, the asymptotic convergence evidently cannot meet the requirements. The complete convergence within a limited time is required [2325]. Among the control methods, iterative learning control (ILC) can achieve full tracking in a limited time for the repetitive system [2628]. Therefore, the ILC algorithm is applied for IOMAS to solve the consensus problem recently, especially for complex MASs [29], MASs with switching topologies and communication time delays [30], and nonlinear MASs [31]. The ILC algorithm can achieve high-precision tracking performance and achieve exact tracking within finite time as the iteration number increases. However, most existing works using ILC are suitable only for IOMAS.

The consensus problem for the fractional MASs with ILC is seldom discussed. In this study, we consider the leader-following consensus problem of FOMAS with general linear leader and design -type and -type fractional-order ILC (FOILC) controllers. The convergence of the proposed controller is analyzed, and the convergence conditions are derived. The convergence conditions can ensure that the tracking errors of all the agents are gradually decreased to a sufficiently small value as the number of iterations increases. Furthermore, the tracking errors can also tend to zero if the initial states of the agents tend to zero.

The rest of the paper is organized as follows. In Section 2, graph theory, fractional calculus, and problem formulation are given. The FOILC and the analysis of convergence for FOMAS are discussed in Section 3. In Section 4, the effectiveness of the proposed algorithm is verified by simulation. Section 5 elaborates the conclusions.

2. Preliminary Knowledge

In this section, some basic definitions and lemmas are introduced, which will be used in the following sections.

2.1. Mathematical Knowledge

The set of real numbers is denoted by , and the set of complex numbers is denoted by . The set of natural numbers is denoted by , and is the iteration number. For a given vector , denotes any vector norm, where . In particular, , , and . For any matrix , is the induced matrix norm. is the spectral radius. Moreover, denotes the Kronecker product, and is the identity matrix.

Remark 1. denotes a set consisting of all functions with mth derivatives that are continuous on the finite time interval .

Definition 1. Continuous vector function , where , is given. The norm of the function is defined aswhere is a positive real number.

Definition 2 (see [32]). The -order fractional integral of the function over the interval is defined aswhere and is the Gamma function defined as .
For the function , the -order Caputo derivative over the interval is defined aswhere and represents the integer part of .

Lemma 1 (see [33, 34]). The function is assumed continuous. Then, the initial value problem is described asThe function in equation (4) is equivalent to the following nonlinear Volterra integral equation:

2.2. Graph Theory

In this study, N FOMAS and a virtual leader are considered, and indicates the agent, . The communications among the agents can be described by a directed graph , where represents the node set, is the set of edges, and denotes the weighted adjacency matrix, that is, , if and otherwise [35]. is the set of neighborhood of the lth node. The Laplacian matrix of the graph is defined aswhere , .

When joining a virtual leader, the relationship between the followers and the leader is expressed by a diagonal matrix , that is, , where if the agent can receive information from the virtual leader. The trajectory of the leader is only accessible to a small portion of the followers due to communication or sensor limitations. The communication among followers is described by the graph . If the leader is labeled by vertex 0, then the complete information flow among all the agents can be characterized by a new graph , where is the new edge set. If the graph is connected or a spanning tree, it can be described as follows.

Definition 3 (see [36]). Among all agents, any agent can receive indirectly or directly the information of the virtual leader in a path. In other words, we do not consider the isolated agent, which means the agent cannot accept any information of the leader, and self-loop agents, which means the agent accepts their own information only.

Lemma 2. If is a connected graph, then the matrix associated with the graph is a positive definite matrix, where is the Laplacian matrix and .
The proof of Lemma 2 is described in [36].

3. Problem Description

A FOMAS consists of agents indexed by and a leader. The dynamic of each follower can be described as follows:where and i are the time and number of iterations; is the state vector of the jth agent; and are the input and output vectors, respectively; are constant matrices with compatible dimensions; and is the -order Caputo derivative of .

The dynamic equations described by equation (7) are rewritten in the following compact form:where , , , , and is the Kronecker product.

The expected trajectory is defined on a finite time interval and is generated by the leader. The dynamics of the leader agent are as follows:where is the continuous and unique desired control input.

We define that is the continuous and unique desired control input for all followers, that is, .

This study aims to find the appropriate control input of each agent with ILC to ensure that the trajectories of each agent are exactly the same as in a finite time , that is, finding appropriate satisfies

In other words, the goal of the function in equation (10) is that each agent in the FOMAS converges to the desired reference trajectory in a finite time with FOILC as the number of iterations increases. The major task is to design a set of distributed FOILC rules for guaranteeing that each individual agent in the network can track the trajectory of the leader under the sparse communication graph .

denotes the distributed information measured or received by the jth agent at the ith iteration. Specifically,where is the element of the weighted adjacency matrix , that is, and is the element of the diagonal matrix and describes the relationship between the followers and the leader, that is, . We define the tracking error between the j-th and the leader fractional-order agents as follows:

Consequently, equation (15) can be rewritten in a compact form:where , , L is the Laplacian matrix of graph G, and .

3.1. Open-Loop -Type Iterative Learning Control

To solve the consensus tracking problem of FOMASs, the following distributed -type ILC for the jth agent is designed as

Based on (12), the distributed measurement in equation (11) can be rewritten as

By using equation (13), the updating law in equation (14) can be rewritten in terms of the tracking errors for all the followers as

For convenience, we define and as the jth eigenvalue of .

To simplify the controller design and convergence analysis, the following assumptions are imposed.

Assumption 1. When the system described by (7) is repeatedly run over the interval , the initial state of the each subsystem in equation (7) is equal to the expected initial state. In particular, is satisfied for all and i.

Remark 2. Assumption 1 is the standard condition for ILC to ensure perfect tracking performance. If the assumption is removed, then the tracking performance will be scarified and the extra system information or additional control mechanisms are required. For instance, the initial state learning control is discussed in [3537]. Notably, perfect tracking can never be achieved without perfect initial condition.

Assumption 2. The communication graph contains a spanning tree with the leader and the followers being the root.

Remark 3. Assumption 3 is the necessary condition of the consensus tracking problem for the FOMAS. If the graph does not contain a spanning tree, then an isolated agent exists and cannot follow the trajectory of the leader because it does not know the control objective. However, the graph does not necessarily contain a spanning tree because the isolated agent can communicate with the selected leader.

Theorem 1. We consider the FOMAS in equation (7) under Assumptions 1 and 2, the communication graph , and -type updating rule in equation (14). If the learning gains and are chosen, thenwhere is the identity matrix with the subscript denoting its dimension, is the Laplacian matrix of , and . The control input and output converge to and , respectively, as the iteration number tends to infinity.

Proof. The convergence analysis is presented as follows. We assume thatwhere and is a vector in which all entries are 1. From (1), we can obtainLemma 1 yieldsAccordingly, we can obtainBy taking the norm on both sides of (11) and noting Definition 1, we haveBy contrast, we haveSubstituting equation (23) into equation (22) yieldsThus, we can find sufficiently large to makeTherefore, equation (24) is reformulated asEquations (7) and (18) yieldBy substituting (19) and (27) into (28), we can obtainBy taking the norm on both sides of equation (29), we haveBy substituting equation (26) into equation (30), we can obtainwhereGiven that , sufficiently large is obtained, thereby satisfying . Thus, from equation (31), we haveWhen the number of iterations increases, that is, , we haveTaking the norm on both sides yieldsAccordingly, we can obtainThe tracking errors of the agents tend to zero with . Thus, converges to . The proof is completed.

3.2. Open-Loop -Type Iterative Learning Control

For the FOMASs described in (7), if we let in (14), then the distributed open-loop -type FOILC becomes a distributed open-loop -type FOILC as follows:

Similar to the -type FOILC, the updating law (37) can be rewritten in the vector compact form

Then, a corollary can be obtained as follows:

Theorem 2. Consider the FOMAS (7) under the directed graph , suppose Assumption 1 and 2 are satisfied. Let the distributed open-loop -type ILC scheme (38) be applied for the system with learning gain satisfyingwhere is an expression about the system parameters and the control parameter as follows: , and then . Namely, for the achievable desired trajectory , the outputs of the FOMASs converge to the desired trajectory when , i.e., .

Proof. The proof process of Theorem 2 is similar to the proof process of Theorem 1.

Remark 4. In the proof, the constant is only a tool for analysis and can be very large, that is, . The reason is that is not used in the implementation and does not affect the actual control performance.

Remark 5. The convergence of FOILC for the FOMASs is uninfluenced by the initial states. We define . Suitable and can be designed such that for all . Thus, an appropriate matrix norm can be designed such that .

Remark 6. According to Lemma 2, the matrix is a positive definite matrix when the directed graph is a connected graph. Then, the matrix is the Hurwitz stable matrix. Accordingly, the gain matrix can be found to satisfy the condition in equation (17).

Remark 7. Compared with -type FOILC, -type FOILC can reduce the deviation of the limit trajectory from the expected trajectory. So, the increased P-type component reduces the adverse effect of the -type component on the disturbance sensitivity and improves the control performance.

4. Simulations

A network consisting of four follower fractional-order agents and a virtual leader fractional-order agent is used to verify the effectiveness of the proposed method. The dynamic function of the agent is a MIMO system. In equation (7), the matrix , , and and are given as follows:and the desired reference trajectory is

The virtual leader and the four followers are labeled as 0, 1, 2, 3, and 4. The information topology-directed graph is shown in Figure 1.

From Figure 1, the weighted adjacency matrix and the in-degree can be obtained as follows on the basis of graph theory:

The Laplacian matrix and the relationship between the followers and the leader are shown as follows:

4.1. The Results of the Open-Loop -Type FOILC

For the open-loop -type FOILC described in (16), and are designed to satisfy the convergence conditions, i.e., and .

By checking the convergence condition and defining , we obtain that

Therefore, the convergence condition is satisfied.

The initial states of all the agents, the initial errors, and the initial inputs are 0, that is, . The input signals at the zeroth iteration are set to zero, that is, , for all agents. Figures 2 and 3 show the tracking performance by the open-loop -type FOILC at 5th, 20th, and 70th iterations. At the 5th iteration, the deviations are very large between the trajectories of the follower and the desired one. As the iteration number increases, the outputs y1(t) and y2(t) of the followers nearly completely track the reference trajectory over the time [0, 1]. By 70 iterations, the maximum tracking errors of y1(t) of the four agents are 0.000326, 0.000843, 0.000916, and 0.000276, and the maximum tracking errors of y2(t) of the four agents are 0.0016, 0.0041, 0.0045, and 0.0014, respectively. And moreover, compared with Agents 2 and 4, both Agents 1 and 3 can get information from the leader directly, so they can converge to the expected trajectory more quickly. We define the maximum errors at the ith iteration as and . Compared with Figure 2(d) and Figure 3(d), the outputs y2(t) converge to the expected trajectory more quickly than y1(t). That is because the trajectory y1d(t) is more complicated than y2d(t) (Figure 4).

4.2. The Results of the Open-Loop -Type FOILC

In order to compare with the performance of the open-loop -type FOILC, the model and the structure of the FOMASs with open-loop -type FOILC are the same as example 1. The control parameters of the open-loop -type FOILC described in (38) are designed as which are the same as of the open-loop -type FOILC. Thus, if is selected, the convergence condition (39) can be calculated as follows:

Thus, the convergence condition in Theorem 1 is satisfied.

The simulation results are plotted in (Figure 4 and 5). Similar to results of the open-loop -type FOILC, as iteration number increases, the outputs y1(t) and y2(t) of the followers nearly completely track the reference trajectory over the time [0, 1]. But, the maximum tracking errors of y1(t) of the four agents are 0.013, 0.034, 0.037, and 0.011, and the maximum tracking errors of y2(t) of the four agents are 0.1, 0.26, 0.28, and 0.085 by 70 iterations, respectively. Obviously, compared with the open-loop -type FOILC, the convergence speed is slower than the open-loop -type FOILC.

4.3. The Results of the Open-Loop -Type and -Type IOILC

In order to reflect the superiority of the fractional controller, we designed an integer-order iterative learning control (IOILC) to apply the FOMASs. The -type and -type IOILC are as follows:

The parameters of the FOMASs are the same as Figure 1 and equations (40) and (41). The learning gains of -type are , , and the learning gain of -type is , respectively, which are the same as the parameters of the FOILC proposed in this paper. The results are shown in Figures 6 and 7. It can be seen that the errors of y1(t) and y2(t) decrease first as the iteration number increases when employing the -type and -type IOILC. But, when the iteration number is greater than 60, the error gradually increases. Compared with IOILC, we can see that the FOILC has better performance.

5. Conclusion

In this study, a consensus tracking problem of fractional-order multiagent systems with general linear models is formulated. Distributed -type and -type fractional-order ILC control algorithms are proposed to solve the consensus tracking problem of the FOMAS. The convergence is proved, and the convergence condition is presented based on graph theory, fractional calculus, and norm theory. It is shown theoretically, together with numerical simulation, that the proposed method can make the agents reach agreement exactly on a common output trajectory over a finite time interval. Through comparison, we find that the -type FOILC can obtain much better performance in terms of convergence speed. Further research may include the closed-loop PDα-type iterative learning control schemes for fractional-order nonlinear multiagent systems, as well as for experimental verifications.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

All authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant nos. 61501155, 61871164, 51777170, and 61871163), the Zhejiang Provincial Natural Science Foundation of China (no. LQ19E070003) the Zhejiang Provincial Key Lab of Equipment Electronics, the Key Research and Development Plan of Shaanxi Province (2018ZDCXL-GY-05-07-01), the Aeronautical Science Foundation of China (20162853026), and the Natural Science Basic Research Plan in Shannxi Province of China (Program no. 2019JM-462).