Abstract

The current paper investigates design and analysis problems of formation tracking for high-order linear time-invariant swarm systems, where communication topologies among agents have leader-following structures and the whole energy supply is limited. Firstly, the communication topology of a swarm system is depicted by a directed graph with a spanning tree, where the communication channels from the leader to followers are directional and the communication channels among followers are bidirectional, and a new formation tracking protocol with an energy integral term and a given upper bound is proposed to achieve formation tracking with the limited energy. Then, sufficient conditions for time-varying formation tracking design and analysis with the limited energy are presented, respectively, which include two/three linear matrix inequality constraints associated with the maximum and minimum nonzero eigenvalues of the Laplacian matrix of a communication topology. Especially, time-invariant formation tracking criteria are further deduced. Finally, two numerical examples are revealed to verify main theoretical conclusions.

1. Introduction

In the last decade, many scientists paid their attentions to design different cooperative control strategies of swarm systems with different application backgrounds, such as explanations for animal flocking phenomena [1, 2], information consensus analysis and design [39], analysis for different source data [10], and formation structure design and maintenance [1013]. Distributed formation is inspired by the biological formation flying without the superintend vertex, where the formation structure maintenance is realized via distributed information transmissions and local interactions. It should be pointed out that distributed formation control has extensive applications in the cooperative operation of satellite clusters, the coordinated assembly of multiple robots, and the collaborative attack of multiple unmanned aerial ships, and Ren [14] summarized the existing formation control approaches including the virtual-structure one, the behavior-based one, and the consensus-based one and pointed out that the consensus-based approach can achieve distributed formation control and can develop the precise mathematical analysis. The consensus-based formation of linear swarm systems as well as their controllability was discussed in [1517].

Based on the time-varying characteristics of formation structures, distributed formation can be divided into time-varying distributed formation and time-invariant distributed formation. For time-varying distributed formation control, the formation geometric structure is time-varying even if the formation structure is formed. For time-invariant distributed formation control, the formation geometric structure is time-invariant once the formation structure is finished. By using the Nyquist stability criterion, time-invariant formation criteria for swarm systems were proposed in [18], where the communication topology was modeled by graph theory and several important features of the communication topology were shown. By a nonlinear formation protocol, sufficient conditions for finite-time time-invariant formation were revealed in [19], where theoretical results were applied into multiple nonholonomic robot systems. Qin et al. [20] studied the impacts of intermittent communications on time-invariant formation control for swarm systems with time delays. In [2123], sufficient and/or necessary conditions of time-varying formation design and analysis were presented, where it was shown that time-varying formation control is more challenging than time-invariant formation control and formation feasible conditions of the time-varying formation are more complex than formation feasible conditions of the time-invariant formation.

According to communication topology structure properties among agents, distributed formation can be divided into two types: the leaderless distributed formation and the leader-following distributed formation. In the leaderless distributed formation, all agents in a swarm system have the same decision weight as shown in [1823]. The leader-following distributed formation is also called formation tracking, where the leader plays a role similar to the superintend vertex and all following vertexes track the leader in a specific formation structure. For first-order swarm systems, Xiao et al. [24] proposed some novel sufficient conditions for formation tracking, where each agent was described as a first-order integrator, which means that the state of one agent is not time-invariant if it does not receive any collaborative information from its neighboring agents. Dong et al. [25] presented several sufficient conditions for formation tracking of second-order swarm systems, where each agent was described as a second-order integrator. In this case, the velocity-similar term converges to a constant and the position-similar term linearly diverges away if the convergence constant of the velocity-similar term is nonzero. In [26], formation tracking criteria for high-order linear time-invariant swarm systems were given, where the formation tracking problem is more challenging since the dynamics of each agent is of high order.

For practical swarm systems, the whole energy supply is usually limited and the above literature studies on formation control did not considered the impacts of the limited energy on formation design and analysis problems. Similar to the guaranteed-cost control of isolated systems, the whole guaranteed-cost function was constructed on the basis of the integral principle in [27, 28], where sufficient conditions for the guaranteed-cost formation were proposed via the linear matrix inequality tool and the nonzero eigenvalues of the Laplacian matrix of the communication topology. Yu et al. [29] presented the guaranteed-cost time-varying formation criteria for swarm systems with external disturbances and time-varying delays, but they did not give an approach to determine the upper bound of the performance function. Especially, it should be pointed out that the whole energy supply of a swarm system was assumed to be infinite when distributed formation control was carried out in [2729]. Actually, the whole energy supply is usually limited for practical swarm systems and its impacts on time-varying formation control are critically important. To the best of our knowledge, the design and analysis problems of time-varying formation tracking for high-order linear time-variant swarm systems with the limited energy and fixed topologies are still open and are not comprehensively discussed.

The current paper investigated time-varying formation tracking problems for high-order linear time-variant swarm systems with the limited energy and fixed topologies. Based on the state errors and the formation function errors among neighboring agents including the collaborative information of the leader, a new formation tracking protocol is presented, where an energy integral term is introduced to guarantee that the practical energy consumption is less than or equal to the whole energy supply. Especially, the dynamics of each agent is described by a high-order linear time-invariant model and the communication topology is depicted as a directed graph with a spanning tree. Furthermore, the whole dynamics of the leader-following swarm system is separated as the dynamics of the leader and the relative dynamics among all agents by constructing the specific nonsingular transformation matrix, and design and analysis criteria for time-varying formation tracking with the limited energy are proposed, respectively, where the relationship matrix between the matrix variable and the limited energy is associated with the Laplacian matrix of a star graph with the leader being the central vertex. Moreover, the time-varying formation tracking criteria are extended into the time-invariant ones by simplifying the different formation feasible condition.

The rest part of the current paper is arranged as follows. The communication topology of a leader-following swarm system is depicted and the problem description of formation tracking for a swarm system with the limited energy is revealed in Section 2. Sufficient conditions for formation tracking design and analysis are proposed in Section 3. In Section 4, two numerical examples are illustrated to demonstrate time-varying and time-invariant formation criteria, respectively. Finally, Section 5 summarizes the key features of our conclusions.

1.1. Notations

The symbols and denote the n-dimensional real column vector and the n-dimensional real matrix space, respectively. The symbol represents the N-dimensional column vector with all components 1. The notation denotes the zero vector or zero matrix with compatible dimensions, and the matrix means an identity matrix with dimension . The symbol represents the Kronecker product. The notation shows that the matrix is symmetric and positive definite.

2. Problem Description

2.1. Modeling Communication Topology

For a swarm system with homogenous agents and a leader-following communication structure, the fixed topology can be described by a weighed directed graph , where stands for the vertex set with the element representing agent k and denotes the edge set with the element denoting the communication channel between agent k and agent i. Without loss of generality, agent 0 is set as the leader and agent k are followers. Furthermore, it is assumed that the communication channels from the leader to the followers are directional and the communication channels among followers are bidirectional, which means that the leader does not impacted by the followers and only some followers can obtain the collaborative information of the leader.

The symbol is used to represent the neighboring agent set of agent k, where each neighboring agent is called a neighbor. The symbol is applied to represent the communication weight between agents k and i, where if agents k and i are not connected, if agents k and i are connected. Especially, it is assumed that there does not exist self-loop; that is, . The matrix is applied to denote the Laplacian matrix of the communication topology, where and . The row sum of the Laplacian matrix is equal to zero; that is, . Moreover, it is supposed that the fixed communication topology has a spanning tree. In this case, zero is the simple eigenvalue of the Laplacian matrix and all the other eigenvalues are positive. One can find more basic concepts and conclusions on algebraic graph theory in [30].

2.2. Modeling Formation Control

The dynamics of each agent is modeled as the following high-order linear time-invariant system:where , is the collaborative state of the leader, and and represent the collaborative state and the control input, respectively. Furthermore, a vector-valued function is applied to design a specific geometric structure for the followers to achieve and maintain, where the element is the formation function of agent , which has the piecewise continuous differentiable property. The formation structure is time-varying if the vector-valued function is time-varying, and the formation structure is fixed if the derivative of the vector-valued function is zero.

In the sequel, a new formation tracking protocol with an energy integral term and a fixed communication topology is presented aswhere , , and denotes the practical energy consumption of swarm system (1) as a whole. Actually, the control input contains two components and , where the component denotes the impacts of the leader on the control input of agent and the component represents the impacts of neighboring followers on the control input of agent .

Let be the energy supply of swarm system (1) as a whole; that is, the energy supply is limited. Now, we, respectively, give the analysis and design definitions of the limited-energy formation tracking of swarm system (1) with tracking protocol (2) as follows.

Definition 1. For any given and the control gain , swarm system (1) with tracking protocol (2) is said to be limited-energy formation tracking achievable if and for any bounded disagreement initial conditions .

Definition 2. For any given , swarm system (1) is said to be limited-energy formation tracking achievable by tracking protocol (2) if there exists a gain matrix such that and for any bounded disagreement initial conditions .
The key objective of the current paper is to give design and analysis criteria for limited-energy formation tracking of swarm system (1) with tracking protocol (2), where both time-varying formation and time-invariant formation are considered.

Remark 1. Tracking protocol (2) has three new properties. The first one is that tracking protocol (2) has an energy integral term associated with all following agents to guarantee that the practical energy consumption is less than or equal to the limited energy. In this case, it is the key challenge to introduce the limited energy into the dynamics of the whole system and to determine the constrained relationship between the energy supply and the control gain. The second one is that the impacts of the leader and neighboring followers are given, respectively, which can be used to decompose the collaborative state of the leader from the whole dynamics of swarm system (1) with tracking protocol (2). The third one is that the formation functions of neighboring agents are involved in tracking protocol (2), which can determine arbitrary piecewise continuous differentiable structures among followers to track the collaborative state of the leader, but it should be pointed out that the formation feasibility is dependent on the dynamics of each agent.

3. Main Results

By the linear matrix inequality, this section, respectively, presents sufficient conditions for limited-energy formation tracking design and analysis of swarm system (1) with tracking protocol (2) and time-varying geometric structures. Then, limited-energy time-varying formation tracking criteria are extended into time-invariant formation tracking cases with limited energy.

Let with and , then one can obtain by tracking protocol (2) thatwhere and

Thus, the dynamics of follower can be rewritten as

Let and , and denotes the Laplacian matrix of the communication topology among followers; then, one can find by (3) and (5) thatwhere

Let denote nonzero eigenvalues of the Laplacian matrix , which are also the eigenvalues of the matrix by the structure of the Laplacian matrix . By the linear matrix inequality technique, the following theorem presents a sufficient condition of limited-energy formation tracking design; that is, a design method of the gain matrix is proposed to make swarm system (1) with tracking protocol (2) and a fixed communication topology achieve limited-energy formation tracking.

Theorem 1. For any given , if and there exists such thatthen swarm system (1) is limited-energy time-varying formation tracking achievable by tracking protocol (2) with .

Proof of Theorem 1. A nonsingular matrix is introduced as follows:Whose invertible matrix isIt can be shown that , so one can deduce thatwhere . Because the communication channels among followers are bidirectional, the Laplacian matrix is symmetric. Since it is assumed that the whole communication topology has a spanning tree, the matrix is nonzero. In this case, the matrix is symmetric and positive definite. Thus, one can find an orthonormal matrix such that , where are nonzero eigenvalues of the Laplacian matrix of the whole communication topology. Let and , then swarm system (1) with tracking protocol (2) can be transformed by (6) intowhere and the column vector is N-dimensional and its k-th component is 1 and zero elsewhere. Because and is orthonormal, it can be found that if , then one can obtain thatwhich means that swarm system (1) with tracking protocol (2) achieves formation tracking.
Furthermore, we give an approach to design such that . Let then the following Lyapunov function candidate is adopted:By taking the time derivative of , it can be derived thatLet , then one hasSince the nonzero eigenvalues of the Laplacian matrix satisfy that , one can deduce that . Therefore, it can be derived by (16) thatDue toone can obtain thatBy (15), (17), and (19), ifthen one can deduce thatThat is, swarm system (1) with tracking protocol (2) achieves formation tracking. In the other work, swarm system (1) is formation tracking achievable by tracking protocol (2) with .
In the sequel, the impacts of the limited energy are dealt with. One can show thatDue to , one can deduce thatSince , it can be found by (22) and (23) thatLet , then one can see thatAccording to , from (22) to (25), one can find thatBecause is the maximum nonzero eigenvalue of the Laplacian matrix of the whole communication topology, it can be obtained by (26) thatThus, one can find that ifthen one has and . In this case, one can deduce by (27) thatSince , one hasDue toIt can be found thatThus, one can find by (30) thatBecause it is supposed that are disagreement, there must exist is nonzero. In this case, one can show thatHence, there exists such thatSince the matrix has a simple zero eigenvalue and positive eigenvalues, one can find by (33) and (35) that can ensure that ; that is, it is required thatLet , then the conclusion of Theorem 1 can be deduced
Theorem 1 gives a sufficient condition for limited-energy formation design, which proposes a design method of the gain matrix of tracking protocol (2). For the case that the gain matrix is given, by the Schur lemma in [31] and the convex property of the linear matrix inequality, the following theorem can be obtained directly on the basis of the Proof of Theorem 1, which presents a sufficient condition for limited-energy formation tracking analysis.

Theorem 2. For any given and the gain matrix , if and there exists such thatthen swarm system (1) with tracking protocol (2) achieves limited-energy time-varying formation tracking.

Moreover, if the formation structure among followers is time-invariant, then one can obtain the following two conclusions by the above analysis, which present limited-energy formation tracking design and analysis criteria, respectively.

Theorem 3. For any given , if , and there exists such thatthen swarm system (1) is limited-energy time-invariant formation tracking achievable by tracking protocol (2) with .

Theorem 4. For any given , the gain matrix , if , and there exists such thatthen swarm system (1) with tracking protocol (2) achieves limited-energy time-invariant formation tracking.

Remark 2. By choosing different formation functions with the piecewise continuous differentiable property, the different formation structures can be achieved for followers in swarm systems with leader-following structures to maintain. From Theorems 1 to 4, it can be found that some formation structures may be unfeasible, and the feasible property depends on the dynamics of each follower; that is, the feasible condition must be satisfied for time-varying formation cases and the feasible conditions and should hold for time-invariant formation cases. Furthermore, the relationship matrix between the limited energy and the matrix variable essentially is the Laplacian matrix of a star graph with all edge weights one, where the leader is the central vertex and followers are not connected with each other. This is coincident with the leader-following communication structure of the whole swarm system. Especially, sufficient conditions in Theorems 1 to 4 include the minimum and maximum nonzero eigenvalues of the Laplacian matrix, whose precise values are difficult to be solved, but they can be estimated by the proposed methods in [32, 33], respectively. It should be pointed out that sufficient conditions in Theorems 1 to 4 can be checked by the FEASP solver of the Matlab’s LMI toolbox in [34].

4. Numerical Examples

This section presents two numerical examples to demonstrate the effectiveness of main conclusions about limited-energy formation tracking for swarm systems with fixed communication topologies.

Example 1 (time-varying formation tracking). Consider a swarm system with six agents and a leader-following communication structure, where the dynamics of each agent is described by (1) withThe communication topology with a leader-following structure is revealed by Figure 1, where agent 0 is the leader, agents 1 to 5 are followers, and the edge weight is set to be 0-1; that is, the weight of the connected edge is one and the weight of the unconnected edge is zero. The initial states of this swarm system are given as follows:The formation structure among followers is determined byIt is clear that the formation functions satisfy that in Theorem 1. Let and . By the FEASP solver, one can obtain by Theorem 1 thatIn this case, one can obtain thatFigure 2 shows the trajectories of , where the curves formed by circle markers are the trajectories of the collaborative state of the leader.
In Figure 3, the state snapshots of the leader and five followers are shown at , , , and , where the leader is depicted by pentagrams and five followers are described by plusses, circles, x-marks, asterisks, and hexagrams, respectively, where it can be seen that five followers achieve a time-varying pentagon to track the collaborative state of the leader. Figure 4 reveals that the leader counterclockwise rotates about 3.18 cycles, where the initial state is represented by a blue pentagram and the final state is depicted by a black pentagram. In Figure 5, it is shown that the practical energy consumption converges to a finite value less than . From Figures 2 to 5, it can be found that this swam system achieves limited-energy time-varying formation tracking.

Example 2 (time-invariant formation tracking). In this example, the leader-following swarm system consists of one leader and eight followers, whose communication topology is shown in Figure 6. In Figure 6, agent 0 is identified as the leader and agents 1–8 as followers, and without loss of generality, the edge weight among the nine agents is also set as 0-1.
The matrix pair of the dynamics equation of each follower is given aswhere is the system matrix of the leader. The initial values of the states for the leader are given as , and eight followers are initialized asThe purpose of this example is to make eight followers track the leader with a time-invariant formation, which is described by a regular octagon geometric structure. The corresponding time-invariant formation function is described byOne can find that the formation function satisfies the condition that and in Theorem 2. Let and ; then, by using the FEASP solver, it can be found by Theorem 2 thatIn this case, it can be deduced thatIn Figure 7, the trajectories of are depicted, where the curves formed by circle markers are the trajectories of the collaborative state of the leader. Figure 8 describes the state snapshots of the leader and eight followers at , , , and , and the leader and eight followers are depicted by pentagram, asterisks, plusses, circles, x-marks, triangles, hexagrams, squares, and diamonds, respectively, where eight followers achieve a regular time-invariant octagon to track the leader. Figure 9 shows that the leader moves along a straight line, where the initial and final states are depicted by a blue pentagram and a black pentagram, respectively. Figure 10 shows the trajectories of the energy consumption and the limited energy. From Figures 7 to 10, one can find that this swam system achieves limited-energy time-invariant formation tracking.

5. Conclusions

For high-order linear time-invariant swarm systems with the limited energy and fixed topologies, a new formation tracking protocol with an energy integral term was proposed to realize formation tracking under the condition that the whole energy consumption is less than or equal to the limited energy supply. By constructing some nonsingular transformation matrix, the whole dynamics of a swarm system with a leader-following topology structure was divided into two independent parts: the dynamics of the leader and the relative dynamics among all agents, and it was shown that the formation tracking is achieved if the relative dynamics is asymptotically stable. Furthermore, the limited energy was introduced into the formation tracking criteria by the matrix variable and disagreement initial states and sufficient conditions for swam systems with the limited energy to achieve time-varying formation tracking and time-invariant formation tracking were proposed, respectively. Especially, these criteria are independent of the number of agents and only include linear matrix inequality constraints, so they are scalable and checkable; that is, the computation complexity does not increase as the number of agents of a swarm system is added. The further works will focus on extending main conclusions for swarm systems with the homogenous dynamics and the fixed topology to swarm systems with the heterogeneous dynamics and switching topologies.

Data Availability

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

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Conceptualization was carried out by Jianye Yang and Xinzhong Han. Methodology was developed by Jianye Yang. Cheng Wang was responsible for formal analysis and validation. Investigation and project administration were performed by Hongtao Dang. Writing-original draft preparation and writing-review and editing were carried out by Cheng Wang and Jianye Yang. Xinzhong Han supervised the study. Cheng Wang and Xinzhong Han were involved in funding acquisition. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

This research was funded by the Key Research and Development Program of Shaanxi (no. 2019GY-025).