Abstract
Group consensus seeking is investigated for a class of discrete-time heterogeneous multiagent systems composed of first-order and second-order agents with both communication and input time delays. Considering two types of system topologies, novel protocols based on the competitive and cooperative relationships among the agents are presented, respectively. By matrix theory and frequency domain analysis method, the sufficient conditions solving consensus problem are obtained. The results show that the achievement of group consensus is bound up with the input time delays, coupling weights between the agents and the system’s control parameters, but it is irrelevant to the communication delays. Finally, numerical simulations are presented to illustrate the correctness of the theoretical results.
1. Introduction
Multiagent systems (MASs) have recently attracted widespread concern. It not only has extensive application in distributed sensor networks, formation control, and artificial intelligence, but also contains great potential value for homeland security and military affairs. Consensus, a typical and crucial issue of the coordination of MASs, refers to that all the agents asymptotically achieve certain global agreement only via the information interactions between them. As an extended problem of consensus, group consensus means that different consensus states can be reached by the different subgroups of a complex system while there is no consensus between those subgroups. It has been widely used to handle multitasking problems and the decomposition of large-scale tasks. So far, there has been a lot of research on consensus and group consensus of MASs, such as [1–5].
Note that most of the existing work is established mainly on homogeneous agents, given that all agents in one complex system possess identical dynamics. However, the difference of the agent dynamics is ubiquitous in practical applications, which is the result of external impacts or interaction restrictions of system. Furthermore, we usually have the demand of using multiple agents with different dynamics to reduce control costs. Thus, the consensus problem of heterogeneous MASs is increasingly being studied, which is also a more complex and challenging task, such as [6], where certain conditions were presented to guarantee the consensus of heterogeneous MASs under fixed and switching topology. In [7], Sun considers the consensus tracking control problem for general linear multiagent systems with unknown dynamics. Jiang [8] investigated the consensus of the heterogeneous networks with the influence of time delays. In [9], for discrete-time heterogeneous systems, Sun discussed both the sufficient conditions to ensure consensus and the consensus equilibrium point based on the graph theory and the matrix theory. Liu [10] proposed some sufficient conditions for heterogeneous multiagent systems to achieve consensus with sampled bounded communication delays. In [11], Bernardo et al. discussed the consensus of time-varying heterogeneous communication delays and applied this strategy for platooning of vehicles.
At the same time, due to the application demands of multitasking, the group consensus problem of MASs is now occupying research hotspots [12–14]. Currently, methods for achieving group consensus include common cooperation-based strategies [15–19] and fewer competition-based strategies [20]. Since the analysis of the group consensus problem of heterogeneous systems is a more complex work than that of homogeneous system, there are few related articles. This incites us to study the group consensus problem for heterogeneous multiagent systems. The main contributions of this paper are as follows. First, a heterogeneous system is established which is composed of agents with first-order and second-order dynamics. And the sufficient algebraic conditions are correspondingly deduced to ensure the consensus of the system. The dynamics of the first-order agents in our system model does not include the virtual velocity, which exists in [15, 17, 19] but reduces the flexibility of the system. Second, utilizing the characteristics of two types of common topologies, effective control algorithms to achieve group consensus are proposed based on competitive and cooperative relationships, respectively. Finally, note that there exists communication delay when the channel is congested and the input delay when the sensor is aging or the computing power is insufficient, while most of articles like [10, 15, 18] did not consider any or just consider one of them. The two kinds of delays are reasonably considered in the agents’ dynamics of our system. Utilizing frequency domain analysis and Gerschgorin circle theorem, we derived the upper bound of input time delay to guarantee the consensus of system. Meanwhile, discrete-time based systems are common in the real world, and current related research is focused on continuous systems [10–19]. Hence, it is also of practical values to investigate the consensus problem of the discrete-time systems. The above innovations make our system model closer to reality and having a wider scope of application.
The remaining of this paper is organized as follows. In Section 2, we first list some preliminaries about graph theory, lemmas, and problem statements. Section 3 contains the main analysis of the group consensus of discrete-time heterogeneous MASs with multiple time delays. In Section 4, several simulations are demonstrated to validate the theoretical results. Finally, some conclusions are drawn based on our work.
Note. Throughout this paper, we use , , denoting the one-dimensional Euclidean space, the n-dimensional real vector space, and the dimensional real matrices, respectively; represents the identity matrix with -dimension; indicates complex numbers sets; the modulus and the real parts of are represented as and ; and indicate the determinant and the th eigenvalue of matrix .
2. Preliminaries and Problems Statements
To assist our work, some definitions, lemmas, and statements are first introduced.
2.1. Graph Theory and Interconnection Topology
Based on graph theory, in a MAS with agents, the agents and the information exchange between them can be described by a weighted directed graph (digraph) , where denotes the node set, denotes the edge set, and denotes the adjacency matrix. denotes that there exists a directed edge from node to node in the digraph so that node can receive information from node . if ; in addition, we assume for all .
Meanwhile, the neighbor set and the in-degree of th node are represented as and , respectively. Denote in-degree matrix as , then Laplacian matrix .
A directed spanning tree consists of nodes with one and only one directed edge for receiving information from other nodes more than a node called the root. A graph which has a directed spanning tree means that there exists such a tree that is composed of several edges and all the nodes of this graph.
Definition 1. There is a graph which can be called bipartite graph when the following two conditions are satisfied:(i)The vertex can be completely divided into two disjoint subsets and .(ii)The two vertexes , associated with each edge of the graph belong to different subset , .
Lemma 2 (see [21]). Consider a complex network with bipartite digraph topology which has a directed spanning tree, it then has when and .
Lemma 3 (see [20]). Suppose that and is a Laplacian matrix. Then, the following four conditions are equivalent:(i) all have positive real parts except a simple zero eigenvalue;(ii) indicates that ;(iii)the consensus is asymptotically achieved for a system ;(iv)the directed graph represented by has at least one directed spanning tree.
2.2. Discrete-Time Heterogeneous Multiagent System
Consider a discrete-time heterogeneous MASs including agents. For the purposes of discussion, we assume the second-order agents are and the first-order agents are represented by the remaining . Hence, the dynamics of the system are as follows:where , , , using , , to represent control input, position state, and velocity state of th agent, respectively.
For heterogeneous MASs, each agent may have second-order and first-order neighbors, which can be described as and separately. Accordingly, all the neighbors of are denoted as . Based on the different orders of the dynamics, we partition the adjacency matrix of the system as , where and denote adjacency matrix consisting only of second-order and first-order nodes, is composed of the adjacency weights among the nodes that are from the second-order to first-order, and indicate the opposite meaning. Meanwhile, we rewritten the Laplacian matrix as , where and indicate the Laplacian matrix considering second-order and first-order nodes with their connections, , and . Thus, we can rewrite matrix as .
Definition 4. The discrete-time heterogeneous MAS (1) and (2) can asymptotically reach consensus if and only if the following two conditions hold for any initial states:(i), for ;(ii), for .
Lemma 5 (see [22]). The inequality , holds for all nonnegative integers .
3. Main Results
Next, we will investigate the group consensus problem for the discrete-time heterogeneous MASs with communication and input time delays under two typical topologies.
3.1. Couple-Group Consensus for Discrete-Time Heterogeneous MASs with Multiple Time Delays under Bipartite Digraph
Most of the existing works, such as [15–19], are based on the agents’ cooperative relationship. From a different point of view, a novel group consensus protocol was proposed in [9] referring the characteristic of the bipartite digraph, which is based on the competitive relationship between the agents. Inspired by their work, two novel delayed protocols based on the competitive relationship are designed, which are listed as follows:andwhere denotes the communication time delays from the agent to agent and denotes the input time delays. are the control parameters of the systems. And we assume if there exists an interaction among nodes and . The closed form of the heterogeneous systems (1) and (2) with the protocols (3) and (4) is shown asand
Theorem 6. Couple-group consensus for heterogeneous MASs (5) and (6) under bipartite digraph topology with a spanning tree can be achieved asymptotically if the following conditions are satisfied:(i)There exists a real number satisfying the equation so that the inequality holds;(ii), , , where and .
Proof. Let (5) and (6) do -transformation. Assuming , it yields thatwhere and represent and after -transformation.
Rewriting (7) to vector form, letBased on the equivalent representation mentioned above, we can rewrite matrix . Then we get Define , then (9) can be rewritten aswhereThe characteristic equation of system (10) is given byAccording to Lyapunov stability theory we know that if the roots of (12) are either at or in the unit circle of the complex plane, the group consensus of the system can be realized. Next, we discuss these two cases separately.
When , . Based on Lemma 2, we know that zero is a simple eigenvalue of . Then the roots of (12) is at .
When , define , where . It follows that where . Based on general Nyquist criteria, one can know that the roots of (12) can be located in the unit circle of the complex plane if and only if the point is not enclosed by the Nyquist curve of .
According to Gerschgorin circle theorem, an area containing all the eigenvalues of the matrix (the union of a set of circles ) is described as . Next, we separately discuss the two cases of the second-order () and the first-order () nodes.
When , the circles are identified as Define , and denotes the center of the disc . Suppose the first cross of the curve on the real axis locates at . According to (15), the following equation is established Since the point , cannot be enclosed in , we obtainBased on Euler’s formula and from (17), we have After some manipulation, we getIt is easy to see that (19) holds for if and only if the following two inequalities are satisfied:As , we can obtainThe inequation (22) can be satisfied if the following two inequations hold,It is easy to know that and from (23) we have . According to (24) and based on Lemma 5, we have .
In addition, since , we get When , the following inequality can also be obtained by Gerschgorin circle theorem: Define . If the point is not enclosed in , it follows thatAfter some manipulations, we haveFrom (28), we getCombining the analysis aforementioned, Theorem 6 is proved.
In Theorem 6, agents with identical input time delay and different communication time delay are investigated. Next, we will extend our work to the case of different input and communication delays. Based on (5) and (6), considering the following systems listed as (30) and (31):andwhere denotes the communication time delay from the th agent to the th agent, denotes the input time delay of the th agent. are the control parameters of the systems.
Corollary 7. Couple-group consensus for the heterogeneous MASs (30) and (31) under bipartite digraph topology with a spanning tree can be achieved if the following conditions are satisfied:(i)there exists a real number satisfying the equation so that the inequality holds;(ii), ;(iii) and , where and , respectively.
Based on Lemmas 2 and 5, as the proof is similar to Theorem 6, we omit it here.
Remark 8. The results in Theorem 6 and Corollary 7 indicate that the control parameters of the system , , and the coupling weight of the agent interaction are critical factors for achieving group consensus. Furthermore, the reduction in either control parameters or coupling weights can make the system more tolerant of input time delays. However, communication time delays do not play a key role in achieving the group consensus.
Remark 9. Corollary 7 shows that the maximum input time delay that each node with respective dynamics can tolerate is different; it is decided by the control parameter of the system and the coupling weights impacted by its neighbors, which has the same dynamics. Similarly, the realization of the group consensus is not affected by the communication delays.
Remark 10. In order to achieve group consensus, many control protocols have been proposed generally based on the cooperative relationship among agents. According to the feature of the bipartite graph, that is, the information exchange of agents only occuring between different subgroups, we design protocols (3) and (4) based on the competitive relationship. This relaxes the condition that system topology needs to satisfy the in-degree balance when using cooperative protocols. And this also gives us a new perspective to investigate the group consensus problems.
3.2. Group Consensus for Discrete-Time Heterogeneous MASs with Multiple Time Delays under Digraph with In-Degree Balance
In order to meet the demand that information exchange of agents exists within the same subgroup, next, we will discuss the case of the heterogeneous MASs with multiple time delays under directed graph with in-degree balance.
Assumption 11. and , where and denote two subgroups.
Remark 12. As in [3, 5, 18], Assumption 11 means that the interaction between the two subgroups is balanced, which is called the in-degree balance. This is necessary when group consensus is discussed in most situations.
Inspired by the control protocol based on the cooperative relationship, the system dynamic is designed and the closed form is listed as follows:andwhere and are the neighbor sets of the agent , which belong to the same and different subgroups of , respectively.
Theorem 13. For the heterogeneous MASs (32) and (33) under directed graph topology with a spanning tree, the group consensus can be achieved if the following conditions are satisfied:(i)Assumption 11 holds;(ii)there exists a real number satisfying the equation so that the inequality holds;(iii), , , where and .
Proof. Let the heterogeneous systems’ dynamics (32) and (33) do -transformation. Assuming , we haveDefine , then (34) can be rewritten aswhereThe characteristic equation of system (35) isAccording to Lyapunov stability theory, the consensus of the system can be realized if the roots of (37) are either at or in the unit circle of the complex plane. Next, we discuss these two cases separately.
When , . Based on Lemma 3, we know that zero is a simple eigenvalue of . Then the roots of (37) are at .
When , define , where , . Based on general Nyquist criteria, one can know that the roots of (37) can be located in the unit circle of the complex plane if and only if the point is not enclosed by the Nyquist curve of . Hence, according to the Gerschgorin circle theorem, we have When , it yields thatNoting that this inequation is similar to those in Theorem 6, then after some derivation processes, we can obtainWhen , the following inequality can also be obtained by the Gerschgorin circle theorem, Then we can get Combining the above analysis, this finishes the proof of Theorem 13.
Based on (32) and (33), consider the following dynamics as (43) and (44):and
Corollary 14. For the heterogeneous MASs (43) and (44) under directed graph topology with a spanning tree, the group consensus can be achieved if the following conditions are satisfied: (i)Assumption 11 holds;(ii)there exists a real number satisfying the equation so that the inequality holds;(iii), ;(iv) and , where and , respectively.
We omit the proof of Corollary 14 as it is similar to Theorem 13.
Remark 15. The dynamics of the agents in our work do not contain the virtual velocity estimation, which exists in many related works [15, 17, 19] for the sake of analyzing heterogeneous systems. In those works, with the view of imitating the dynamics of the second-order agents, the velocity estimation is added to the first-order agents, but this requires extra computation and wastes system’s resources.
Remark 16. In this paper, the topology of the systems we considered is either a bipartite digraph or a digraph satisfying in-degree balance, both of them has a spanning tree. They seem to be more specific. But in fact, group consensus for heterogeneous systems usually can hardly be achieved unless adding some stronger conditions. As in [15–19], the systems’ topology is undirected or containing a spanning tree; meanwhile, in those papers in order to achieve the group consensus, additional assumptions are specified, like in-degree balance and the geometric multiplicity of zero eigenvalues of the system’s Laplacian matrix which are at least two. It will be a challenge to discuss the problem of the group consensus for heterogeneous MASs under a more general condition.
Remark 17. Note that all the results in this paper are theoretically obtained in the case , , . Utilizing the Kronecker product, they can be generalized to , , .
4. Simulation
In this subsection, the effectiveness of the obtained results is demonstrated by several simulations.
Example 18. Figure 1 shows a heterogeneous system with the bipartite digraph topology which has a spanning tree. As known, the whole system is divided into two subgroups: and . Hereinto, agents 1 and 2 and 3, 4, and 5 belong to these two subgroups. Without loss of generality, we further assume that agents 1 and 5 are second-order and that the dynamics of the remaining agents are first-order. In our scenario, the dynamics of agents in the same subgroup do not need to be the same, which is different from the case in [16–19]. Obviously, those are special cases of ours.

For simplicity, we assume each edge’s weight as , and it is easy to obtain the in-degree of each agent: , , , , and . Based on Theorem 6, we set , , and , then the upper bound of input time delays is obtained, . In the simulation, we choose and all the conditions in Theorem 6 are satisfied now. Since the realization of the group consensus is independent from the communication time delays, for convenience, we choose their values to be the same as . The trajectories of the agents of the systems (5) and (6) are shown in Figure 2. The results indicate that the systems’ couple-group consensus is reached asymptotically.

(a)

(b)
From Corollary 7 and according to the control parameters given above, the upper bound of input delays can be calculated for each node: , , , , and . In this case, we choose , , , . It is clear that the conditions in Corollary 7 are all satisfied here. Figure 3 demonstrates that the couple-group consensus is realized. Furthermore, if we set one of the input delays, for instance, to violate the upper bound we proposed, from the trajectory in Figure 4, we know that the system is divergent.

(a)

(b)

(a)

(b)
Next, we will testify Theorem 13 and Corollary 14.
Example 19. Figure 5 shows an in-degree balance graph and , , , , and . Set , , , and . Obviously, the conditions of Theorem 13 are satisfied. The trajectories of the agents are shown in Figure 6. Consider the systems (34) and (35), and choose different input time delays as , , , , and . Here, Corollary 14 is satisfied, and the trajectories are shown in Figure 7. From Figures 6 and 7 we know that the couple-group consensus is achieved asymptotically. When set exceeds the upper bound, the system diverges (Figure 8).


(a)

(b)

(a)

(b)

(a)

(b)
Remark 20. In this paper, we designed a competition-based control protocol under bipartite graph and a cooperation-based protocol under in-degree balance graph to realize group consensus. Competitive protocol will lead to a couple of agents which belong to two different subgroups apart from each other. In this case, agents in the same subgroup have no interactions, and the final states of the system are only two and opposite, which is called the couple-group consensus. Cooperative protocol will lead the agents in one subgroup to be close to each other, while the different subgroups need to satisfy the in-degree balance, which means all agents from one subgroup impose a weight of zero on any node within another subgroup. In this case, more than two system states can be achieved. Both approaches have their own limitation of topology. Next, we plan to combine the competition and cooperation ideas to study the group consensus problem under more general system topology.
5. Conclusion
Under two forms of typical topologies, the couple-group consensus for the discrete-time heterogeneous MASs with input and communication time delays is investigated. For the bipartite graph, group consensus algorithms based on the competitive relationship between agents are designed; for the graph satisfying in-degree balance, we propose the algorithms based on cooperative relationship accordingly. Utilizing matrix theory and general Nyquist criteria, we theoretically propose some sufficient algebraic criteria and the upper bound of the input time delays for ensuring the consensus. From them, we find that the achievement of the couple-group consensus for the heterogeneous systems depends on the input time delays of agents, the control parameters of the systems, and the coupling weights between the agents, but it is independent of the communication time delays. The simulations verify the correctness of our theoretical results. Our future work will extend to more complicated group consensus problems of the heterogeneous multiagent systems, for example, considering sampled information or switching topology which commonly exist in reality.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported in part by the National Natural Science Foundation of China under Grant no. 61876200, in part by the Postgraduate Scientific Research and Innovation Project of Chongqing under Grant no. CYS18246, in part by the Natural Science Foundation Project of Chongqing Science and Technology Commission under Grant no. cstc2018jcyjAX0112, and in part by the Key Theme Special Project of Chongqing Science and Technology Commission under Grant nos. cstc2017zdcy-zdyfx0091 and cstc2017rgzn-zdyfx0022.