#### Abstract

In this paper, we consider the formation tracking problem for multiagent systems with diverse communication time-delays by using iterative learning control (ILC) method based on the frequency domain analysis. A first-order ILC law for multiagent systems with diverse communication time-delays is first proposed and its convergence conditions are given by the general Nyquist stability criterion and Gershgorin’s disk theorem. Then, in order for the system to track accurately, a second-order ILC law is presented. The conditions for system tracking with zero error are established. Numerical simulations show that the proposed ILC laws for multiagent systems with diverse communication time-delays are able to achieve effectively formation tracking. And the convergence speed remains the same as the learning control algorithm without communication delay.

#### 1. Introduction

Multiagent systems can make it more efficient to complete some parallel and complicated tasks than a single agent and more and more scholars have begun to pay attention to multiagent systems cooperation [1–3]. As well known, ILC has been widely utilized for dynamic systems that operate repetitively on the finite time interval, because of its low complexity and high precision. Therefore, the multiagent systems that can operate repetitively using ILC method are a wise choice to achieve flexible task. In 2009, Ahn et al. firstly used ILC to implement circular formation for multiagent systems [4]. Afterward, ILC for multiagent systems was widely used in practical applications. For example, in [5], multiple satellites cooperatively kept the trajectory around the earth; in [6], multiple robots formation moved in a shape; in [7], multiple high-speed trains cooperatively controlled driving speed; in [8, 9], multiple UAVs flew for trajectory tracking. And all the above-mentioned papers were based on ILC method.

However, note that all the above-mentioned results were based on an ideal communication environment, in which data transmitted between multiple agents will not be mistaken. Due to the limited bandwidth of communication channels, noise interference, and signal fading, multiagents inevitably have communication time-delays during signal transmission [10]. Since communication time-delays extensively exit, the performance of systems may be greatly reduced. For example, in [11–14], multiagent systems with communication time-delays have attracted great attention. Since ILC is introduced into multiagent systems, more and more scholars have paid attention to ILC for multiagent systems with communication time-delays. An ILC tracking strategy was proposed for nonlinear multiagent systems with communication time-delays in [15]. Using the 2D analysis, [16] obtained convergence conditions in terms of linear matrix inequality and addressed both exponential convergence and monotonic convergence problems associated with the ILC process corresponding to the proposed consensus learning protocol. In [17], a distributed ILC algorithm was proposed for multiagent systems subject to 2D switching topologies and varying communication time-delays. In [18], an adaptive iterative learning control method, which includes a P-type feedback term and an iterative learning term along the iteration axis, was designed for the consensus problem of homogeneous multiagent systems with state time-delays. In [19], an ILC scheme for multiagent systems with one-step random time-delay was proposed. All the above-mentioned papers have proposed suitable learning law and convergence conditions to guarantee convergence of system errors. However, the bounds of tracking errors were subject to the accuracy of the estimated delay in [15]. The works [16, 18] discussed consensus problem of multiagent systems not formation tracking. And [19] showed that the convergence speed got slower as the increase of the time-delays rate. Therefore, time-delays compensation method for formation tracking of multiagent systems with diverse time-delays is a more challenging issue.

Motivated by the above observations, this paper proposes an ILC time-delays compensation method for multiagent systems. And we can keep convergence speed as multiagent systems without time-delays. A frequency domain ILC model is established for multiagent systems with linear SISO discrete time dynamics. Since ILC is a two-dimensional system with evolution along two independent axes, frequency domain and iteration domain, we discuss the convergence of arbitrary frequency in iteration dimension by exchange variable frequency and parameter iteration. Then we investigate eigenvalues of the characteristic equation to analyze the condition of convergence by taking z-transformation in iteration domain. Because the characteristic equations are in matrix form, it is difficult to analyze their eigenvalues. To this end, we use the generalized Nyquist criteria and Gershgorin disc theorem to analyze the range of eigenvalues and obtain the conditions for system convergence. The results show that the convergence condition of the system has nothing to do with the communication time-delays. The system errors cannot converge to zero but a certain value. In order to compensate for the influence of time-delays, this paper proposes a method using second-order ILC to make the system error converge to zero. The same frequency domain analysis method as first-order ILC was used to analyze the convergence conditions of the second-order ILC for multiagent systems with communication time-delays. Simulations demonstrate that second-order ILC can effectively compensate for the effects of communication delay.

The remainder of the paper is organized as follows. In Section 2, an ILC model is established for multiagent systems with linear SISO discrete time dynamics based on frequency domain. Moreover, we proposed the control objective which is formation tracking for multiagent systems. The conditions of system convergence are obtained in Section 3 when using first-order ILC and second-order ILC, respectively. In Section 4, numerical simulations show the effectiveness of the proposed protocol.

*Notations. *1 and 0 denote the column vectors of appropriate dimensions whose elements are all ones and all zeros, respectively; denotes an identity matrix with the required dimensions.

#### 2. Problem Formulation

##### 2.1. System

Consider the multiagent systems consisting of agents labeled from 1 through , whose communication topology graph is denoted by . And the dynamic of each agent is a linear SISO discrete time system:where represents the th agent, is the discrete time index, is iteration index, is the th agent state of the th iteration in time , is the control input signal, and is the output. Multiple agents have the same dynamics , , and . Without loss of generality, let system (1) relative degree be 1; that is, . And the dynamics for each agent is stable; that is, . Taking the z-transformation of system (1), we getwhere with is frequency, and is the th agent initial state of the th iteration; let , then (2) can be rewritten as

*Assumption 1. *The initial resetting condition for all agents and all the iteration satisfied desired input. Then, can be abbreviated as .

*Remark 2. **Assumption 1* is common in ILC for multiagent systems such as [16, 19]. If the condition is not met, we can think of it as a robust problem against initial shifts. In our future work, the ILC with initial-state learning of formation tracking control will be one of our objectives to consider.

##### 2.2. Control Objective

For system (1), the objective to realize in this paper is that multiagent systems can move in the desired formation tracking reference trajectory as the increase of iteration. That is, for , , we have , and is desired output of the th agent and denoted aswhere , is desired reference trajectory of multiagent systems. In practice, not all agents can obtain , but a portion of agents can. To this end, we denote the reference-accessibility matrix , which is a diagonal, nonnegative, real matrix. If the th agent has direct information about , then ; otherwise, . is the desired output deviation from desired reference trajectory of the th agent. If multiple agents have the same desired output deviation, then the agents will have the same output. So . represents the desired relative formation between the th agent and the th agent. It is worth pointing out that each agent can obtain its own output and its own desired output deviation . For convenience, we denoteWe denote the output error of the th agent in the th iteration:Taking z-transformation for (6), we getWhen it is satisfied , that is, , multiagent systems realize formation tracking control.

##### 2.3. Preliminaries in Graph Theory

Let = () be directed graph, where is the set of vertices, is the set of edge, and is the weighted adjacency matrix of the graph . If there is an edge from the th agent to the th agent, that is, , then , representing the fact that information is transmitted from th agent to th agent. Otherwise, . Moreover, we assume that . The index set of neighbors of node is . The Laplacian matrix of graph is denoted as , where diag with . A path in the directed graph is a finite sequence , and . If there is a vertex that can be connected to all other vertices through paths, then is said to have a spanning tree, and this vertex is called the root vertex.

*Assumption 3. *The graph of multiagent system (1) is a directed graph. And its graph is connected; that is, each agent has a path to the other agents so that they can exchange information.

Lemma 4. *Vector is a right eigenvector of the Laplacian with eigenvalue ; i.e.,*

Lemma 5 (see [20]). *If an irreducible matrix is weakly generalized diagonally dominant and at least one of the rows is strictly diagonally dominant, is nonsingular.*

*Remark 6. **Assumption 3* implies is a matrix that satisfies the conditions of* Lemma 5*. In , the magnitude of the diagonal entry is , and the sum of the magnitudes of all nondiagonal entries is . Since and at least one of the rows is satisfied , matrix is nonsingular.

##### 2.4. Communication Time-Delays

When multiagent systems cooperate through the network, there are time-delays in transmission due to the limited bandwidth of communication channels, noise interference, and signal fading. Therefore, we introduce communication time-delays into multiagent systems to research more close to the practical systems. A block diagram of multiagent systems with communication time-delays is illustrated in Figure 2 if the communication graph of multiagent systems is shown in Figure 1. When the th agent can obtain from the th agent, we denote as the signal which the th agent received from the th agent:where is communication time-delay when the th agent sends signal to the th agent. .

#### 3. Compensation to ILC for Multiagent System with Communication Time-Delays

In order to make the multiagent systems more close to the practical engineering applications, we analyze the ILC for multiagent systems with communication time-delays first. The result shows the system error of first-order ILC with communication time-delays cannot converge to zero. Second, the method to compensate communication time-delay is given. And the convergence analysis is obtained in Section 3.2.

##### 3.1. First-Order ILC for Multiagent Systems with Communication Time-Delays

For multiagent systems with communication time-delays, a first-order ILC law is proposed aswhere and is a gain to be determined. This ILC law is multiple agents joint learning law. The control input of the th agent needs previous iterative control input, self-error, and adjacent agent information. In (10), is communication time-delay mentioned above.

Taking the z-transformation of (10), we getBy (7), (11) and* Assumption 1*, we have thatLet ; then (12) can be rewritten asFurther, we can show thatwhere is

Theorem 7. *Consider multiagent system (1); let Assumptions 1 and 3 hold and the learning law (10) be applied. Then the system error of multiagent systems can converge, given that any one of the following conditions is satisfied for all :where*

*Proof. *As we know, ILC is employed for dealing with the repeated tracking control. For this kind of machines, let them learn as people by using previous errors correct this control input. Thus, by repeating the learning process, the tracking performance can be achieved with high accuracy. Since ILC is a two-dimensional system with evolution along two independent axes, time domain and iteration domain, the scholars Fang et al. proposed an ILC 2D model [21, 22]. This paper considers a two-dimensional system of frequency domain and iteration domain after taking z-transformation. Considering as two-dimensional function , we discuss the convergence of arbitrary frequency in iteration dimension by exchange variable frequency and parameter iteration. Thus, we can exchange variable and parameter . And can be rewritten as ; that is, . Then, (14) can be rewritten asTaking the z-transformation of (18), we getFurther, we get . The characteristic equation isWhen , . From* Remark 6*, we can get that is nonsingular; that is, . Since , is not the root of characteristic equation.

When , both sides of (20) are divided by :Next, we need to prove that all roots of (21) have modulus less than unity. Let . Based on the general Nyquist stability criterion [23], the roots of (21) have modulus less than unity, if the eigenloci ofdo not enclose the point for , , and . By Gershgorin’s disk theorem, we have for all , , and , whereSo the eigenloci do not enclose the point for , , and when the point with is not in the disc for all , , , and . That is, for all , , and , when . We denote asNote that for all , , and . So, as long as for all , , , and . Further, we can show thatUsing the conditions of* Theorem 7*, we can prove for all , , , and in Appendix. Then eigenloci of for all , , and do not enclose the point . Thus the roots of (21) have modulus less than unity. That is, the system achieves a formation tracking asymptotically.

In addition, let the final value of be . Multiply on both sides of and take the limit:By final value theorem, we getSince , is invertible. And is satisfied; thus . This completes the proof.

When there is not communication time-delay in multiagent systems, we can get system error converging to zero combining equation (27) and* Lemma 4*. Otherwise, system error cannot converge to zero and converges to a certain value about desired reference trajectory, graph, and time-delays as shown in (27).

##### 3.2. Second-Order ILC for Multiagent Systems with Time-Delays

In this subsection, we analyze convergence condition of second-order ILC for multiagent systems. Using previous two errors to correct the control input can compensate the influence of communication time-delays. Therefore, we propose the second-order ILC law:where is first-order learning gain, is second-order learning gain, , and .

Theorem 8. *Consider multiagent system (1); let Assumptions 1 and 3 hold and the learning law (28) be applied. Then the formation tracking objective is achieved, given that , , and any one of the following conditions are satisfied for all :where*

*Proof. *Taking the z-transformation of (28), we getBy (7), (31), and* Assumption 1*, we have thatLet ; (32) can be rewritten asExchanging variable frequency and parameter iteration , can be rewritten as in (33):Taking z-transformation of (34), we getFrom* Theorem 8*, . Further, we have thatThe characteristic equation of system isLet . Since and , is the only root of (37).

When , both sides of (37) are divided by :Next, we need to prove that all roots of (38) have modulus less than unity. Similar to the proof of* Theorem 7*, the conditions of* Theorem 8* are easily verified.

Taking the limit of (36) yieldsBy final value theorem, we getSince as proved above, . Based on Sylvester inequality, we get . Therefore, ; that is .* Theorem 8* is thus proved.

The result of* Theorem 7* shows that the error of multiagent system cannot converge to zero considering diverse communication time-delays. In* Theorem 8*, using second-order ILC and the fact that the second-order learning gain is opposite of the first-order learning gain, the error of multiagent system can converge to zero considering diverse communication time-delays. In addition, the convergence speed can be kept as convergence speed of multiagent systems without time-delays.

#### 4. Simulation

In this section, numerical simulations are presented to illustrate the first-order and second-order ILC for multiagent systems. Consider a multiagent system consisting of 4 agents in the directed graph as shown in Figure 3. It can be seen that the graph is connected graph. In addition, only the 1st agent can obtain the desired reference trajectory.

The Laplacian matrix of graph and the reference-accessibility matrix areFrom and , we get . Let desired reference trajectory be . The desired output deviation of each agent is chosen asSet the initial state of 4 agents as 875, 575, 175, and -125. It satisfies* Assumption 1* that the initial resetting conditions for all agents and all the iterations meet desired input. When dynamics of system is Case 1, , we set . It can be seen that above learning gains satisfy condition 1 of* Theorem 7*. When dynamics of system is Case 2, , we set . It can be seen that above learning gains satisfy condition 2 of* Theorem 7*. In Case 1 or Case 2, when and iteration is 200, the trajectory of the multiagent system without communication time-delays is shown as Figure 4. From this figure, it can be seen that the desired formation is well achieved to track desired reference trajectory.

In order to measure the formation accuracy quantitatively, the disagreement among all agents on their output errors is defined as . When , multiagent system achieves formation tacking. Based on this, Figure 5 shows formation performance of the system at the first 200 iterations in Cases 1 and 2. From Figure 5, when iteration is around 120, the multiagent systems error in Case 1 can converge to zero and when iteration is around 150, the multiagent systems error in Case 2 can converge to zero.

When there are communication time-delays in the multiagent system, we set 3 time-delays graphs of the multiagent system as shown in Figure 6.

**(a)**

**(b)**

**(c)**

From Figure 7, the system errors cannot converge to zero but a certain value as the increase of the iteration when there are time-delays. In addition, the bigger the communication time-delays are, the bigger the bounds of tracking errors are. Multiagent system cannot achieve formation tracking and the trajectory of multiagent system is as shown in Figure 8 when time-delays graph 1 holds.

**(a)**

**(b)**

By* Theorem 8*, we use second-order ILC and set when the time-delays graphs of multiagent systems are as shown in Figure 6. It can be seen that Case 1, , satisfies conditions 1 and 3 of* Theorem 8* and Case 2, , satisfies condition 2 of* Theorem 8*. The system errors are as shown in Figure 9. From Figure 9, when iteration is around 120 in Case 1, all the multiagent systems errors with communication time-delays and without communication time-delays can converge to zero. When iteration is around 150 in Case 2, all the multiagent systems errors with communication time-delays and without communication time-delays can converge to zero. Therefore, second-order ILC formula can keep the same learning error convergence speed, no matter what time-delays are. In Case 1 or Case 2, multiagent systems can achieve formation tracking and the trajectory of multiagent systems is shown as Figure 10.

**(a)**

**(b)**

#### 5. Conclusion

When using first-order ILC to implement the formation tracking for multiagent system with communication time-delays, the system error cannot converge to zero but a certain value. Using second-order ILC and making second-order learning gain be opposite of first-order learning gain can compensate influence of communication time-delays. Based on the frequency domain ILC model for multiagent systems, we analyze the convergence condition by considering the convergence of arbitrary frequency in iteration dimension. We use the generalized Nyquist criteria and Gershgorin disc theorem to analyze the range of eigenvalues of the system characteristic equation and obtain the conditions for system convergence. The second-order ILC proposed in this paper for compensating for the influence of time-delays not only can make the system error converge to zero, but also can ensure that the convergence speed is the same as the convergence speed for system without time-delays. The simulation results also verify the effectiveness of the conclusion.

#### Appendix

is a quadratic function. Since and with* Assumption 1*, and . That is, and has minimum value for . The minimum value of for isTaking and for , we get . This example can prove that is not valid for all , and . Thus, for all as long aswhereWe denote is function of two variables. We need to obtain the extremum of . Taking the partial derivative with respect to , we getTaking the partial derivative with respect to , we getThen we take second derivative of .Next we denote : is function of two variables. Similarly, we do the same operation of as . Taking the partial derivative with respect to , we getTaking the partial derivative with respect to , we getThen we take second derivative of .Fromandwe get () and () which are the stationary points of both and for all , , and .

When , second derivative of isThere is no such and satisfied and . Then (A.2) is not satisfied.

When , the second derivative of isand the second derivative of isIf condition 1 and condition 2 of* Theorem 7* are satisfied, thenis satisfied and , are satisfied. That is, only and are minimum values of and , respectively. Overall, if any one of the two conditions in Theorem holds, then (A.2) is satisfied. That is, for all , , , and is proved.

#### Data Availability

No data were used to support this study.

#### Conflicts of Interest

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

#### Acknowledgments

This work is supported by research grants from the National Natural Science Foundation of China (Nos. 61673253, 61271213) and the Key Support Projects of Shanghai Science and Technology Committee (16010500100).