Abstract

This paper discusses limited-budget time-varying formation design and analysis problems for a high-order linear swarm system with a fixed communication topology. Firstly, the communication topology among agents is modeled as an undirected and connected graph, and a new formation control protocol with an energy integral term is proposed to realize formation control and to guarantee the practical energy assumption is less than the limited energy budget. Then, by the matrix inequality tool, sufficient conditions for limited-budget formation design and analysis are proposed, respectively, which are scalable and checkable since they are independent of the number of agents of a swarm system and can be transformed into linear matrix inequality constraints. Moreover, an explicit expression of the formation center function is given, which contains the formation function part and the cooperative state part and is not associated with the derivatives of the formation functions. Finally, a numerical simulation is shown to demonstrate the effectiveness of theoretical results.

1. Introduction

Recently, many researchers from system and control fields paid their attentions to the distributed coordination of swarm systems, which have potential applications in many aspects, such as flocking algorithm and theory [1, 2], synchronization analysis and design [311], formation control [1214], and distributed parallel computation [1519]. Formation control of a distributed swarm system is inspired by the biological formation behaviors without the superintend node, where it is needed that some specific geometric structure is achieved and maintained by local interactions among all animals. Formation control of distributed swarm systems has extensive application in the coordinated maneuver of multiple androids and coordinated attack of multiple unmanned surface ships. In [20], formation control methods were divided into three classes: the virtual structure method, the behavior-based method, and the synchronization-based method, and it was pointed out that the synchronization-based method is distributed and can overcome some limitations of the virtual structure method and the behavior-based method.

According to the time-varying property of geometric structures, formation is divided into time-invariant geometric structures and time-varying geometric structures. For time-invariant geometric structures, the time derivative of the formation function of a swarm system as a whole is zero; that is, the formation structure is time-invariant once a geometric structure is formed by all agents. For time-varying geometric structures, the time derivative of the formation function of a swarm system as a whole is nonzero, which means that the formation structure is time-varying even if a geometric structure is formed by all agents. By the Nyquist stability criterion, time-invariant formation criteria were proposed in [21]. Du et al. [22] presented time-invariant formation control criteria by using nonlinear formation protocol and applied theoretical results into multiple androids systems. Qin et al. [23] investigated the influences of time delays on the time-invariant formation of swarm systems and gave several sufficient conditions for time-invariant formation control. In [2426], time-varying formation analysis and design for swarm systems were addressed and some important conclusions for time-varying formation were presented. Overall, time-varying formation control is more complex than time-invariant formation control, and some specific formation feasible conditions are needed.

In [2729], formation control design and analysis problems for swarm systems with different topology structures were investigated, where the influences of the limited energy supply were not dealt with, which is critically important in practical applications. In [30, 31], the concept of the guaranteed-cost control was introduced to analyze the impacts of the energy constraint and gave several guaranteed-cost formation criteria. Yu et al. [32] investigated the guaranteed-cost time-varying formation control problem, where the impacts of unknown external disturbances and time delays were considered and the guaranteed-cost time-varying formation criteria were proposed, but the performance cost cannot be determined by their methods. In [3032], the energy supply of a swarm system as a whole cannot be given previously. However, the whole energy supply is usually limited in practical swarm systems and it is critically important to discuss the influences of the whole energy supply on time-varying formation design and analysis problems. To the best of our knowledge, limited-budget formation control for high-order linear swarm systems with fixed communication topologies is still open and is not comprehensively addressed and analyzed.

This paper deals with time-varying formation design and analysis for high-order linear swarm systems with the limited energy budget. A formation control protocol is proposed on the basis of the state errors and the formation function errors among neighboring agents, where an energy integral term is introduced to ensure that the whole energy budget is larger than or equal to the practical energy assumption. Furthermore, the dynamics of a swarm system with a fixed communication topology is divided into two parts by the state space decomposition, which can be used to describe the macroscopic motion as a whole and the microcosmic motions between any two agents. By constructing the relationship between the matrix variable and the energy budget, limited-budget formation design and analysis criteria are proposed, respectively, where the scalability property and the checkable property are discussed in detail. Moreover, the formation center function is determined, which involves two independent parts associated with the average value of the time-varying formation functions and the average value of the initial conditions of cooperative states, respectively.

Compared with the recent works on the time-varying formation control, the innovation of this paper is twofold. Firstly, this paper considers the limited-budget time-varying formation control, where the energy budget of the whole swarm system is given and is limited previously. The energy budget was not considered and the practical energy consumption cannot be limited by the energy budget in the time-varying formation as shown in [2429]. Secondly, the effect of the energy budget to the matrix variable of the matrix inequality conditions is considered and determined, which can introduce the energy constraint into the design of the gain matrix of the time-varying formation control protocol. In contrast, the formation control methods in [2429] cannot give the specific interaction mechanism for the energy budget.

The main arrangement of this paper is given as follows. The communication topology of a swarm system is modeled and the problem description of limited-budget formation for a swarm system is shown in Section 2. The main results of limited-budget formation are presented in Section 3. A numerical simulation is given to demonstrate the limited-budget formation criterion in Section 4. Finally, Section 5 summarizes the main conclusions of this paper.

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

2. Problem Description

2.1. Communication Topology Model

For a swarm system with identical agents, the fixed communication topology is depicted a weighed graph , where is the node set with representing agent k and is the edge set with denoting the interaction channel from agent k to agent i. The notation is used to denote the neighbor set of agent k. The symbol denotes the communication weight from agent k to agent i, where and when agent k and agent i are not connected and when agent k and agent i are connected. The Laplacian matrix of the communication topology is denoted as with and . If the fixed communication topology is connected, then the Laplacian matrix is symmetric and positive semidefinite; that is, zero is its simple eigenvalue and all its nonzero eigenvalues are positive. More basic concepts and conclusions on graph theory can be found in [33].

2.2. Swarm System Model

The dynamics of each agent is depicted in the following high-order linear system:where , , and , and and , respectively, denote the cooperative state and the control input. Moreover, a vector-valued function is used to depict a specific formation structure for swarm system (1) to maintain, where the piecewise continuous differentiable components are formation functions of agent . The formation is time-varying when is time-varying and the formation is fixed when is time-invariant.

In the following, a limited-budget formation control protocol with a fixed communication topology is proposed as follows:where , , and , and denotes the practical energy consumption of swarm system (1) as a whole. Let be the maximum energy supply of swarm system (1); that is, the energy budget is limited. The definition of the limited-budget formation control of a swarm system is proposed as follows.

Definition 1. For any given , swarm system (1) is said to be limited-budget formation achievable by control protocol (2) if there exists a gain matrix such that and for any bounded disagreement initial conditions , where is said to be the formation center function.
It can be found from Definition 1 that the formation control is equivalent to the consensus control if , which means that the consensus control can be regarded as a special case of the formation control. Besides, the vector denotes the desired formation shape of the swarm systems and the state of all agents should track with the formation center function.
The main objective of this paper contains two aspects. The first one is to design the gain matrix such that swarm system (1) with control protocol (2) achieves limited-budget formation. The second one is to determine an explicit expression of the formation center function.

Remark 1. It should be pointed out that the vector-valued function can be used to design the arbitrary geometric structure with the piecewise continuous differentiable property, but it may be unfeasible for some agent dynamics of a swarm system; that is, the formation achievable property is associated with the mechanism structure of practical dynamic agents in a swarm system. Moreover, the energy consumption of a practical swarm system is limited and control protocol (2) imposes an integral term to guarantee that the actual energy consumption is less than the maximum energy supply. This makes it challenging to construct the relationship between the gain matrix and the limited-budget , that is, how to design the gain matrix such that swarm system (1) with control protocol (2) achieves formation under the condition that the maximum energy supply is larger than or equal to the actual energy consumption.

Remark 2. Notice that the practical energy consumption is constructed according to the weight matrix and is an integral quadratic associated with . In the design of the gain matrix, the matrix is often considered as a diagonal matrix, whose th element in the diagonal line denotes the weight to the th channel of the control protocol. By adjusting the value of the matrix , the contribution of the control protocol to the practical energy consumption will be changed, and is utilized to determine the matrix variable of the matrix inequality conditions. Moreover, compared with the formation control protocols of [2429], control protocol (2) introduces the practical energy consumption to save the control energy under the limited-budget constraint; that is, the control gain should be designed to satisfy that with the weight matrix .

3. Main Results

Based on the matrix inequality tool, this section gives sufficient conditions for limited-budget formation design and analysis, respectively. Then, an explicit expression of the formation center function is proposed, which contains two parts: the formation function one and the cooperative state one.

For , let , then it can be obtained by (1) and (2) that

Let then one can transform swarm system (3) into

It is assumed that the communication topology of swarm system (1) is undirected and connected, so the Laplacian matrix is symmetric and positive semidefinite, which owns a simple zero eigenvalue and positive eigenvalues. Hence, one can find an orthonormal matrix such that , where denote the positive eigenvalues of the Laplacian matrix . One can set that then one can transform swarm system (4) intowhere represent N-dimensional column vectors with the ith part being 1 and zero elsewhere.

The following theorem shows a sufficient condition of limited-budget formation design for swarm system (1) with a fixed communication topology on the basis of the matrix inequality tool, which proposes a design approach of the gain matrix such that swarm system (1) with control protocol (2) achieves limited-budget formation.

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

Proof. of Theorem 1. DefineOne can find thatSince , it can be shown thatthat is,In this case, the transformation matrix is orthonormal and , so it can be found from (8)–(11) that the component can be used to describe the formation center function as the time tends to infinity and swarm system (1) with control protocol (2) achieves formation if and only if all the components tend to zero as the time tends to infinity.
In the following, an approach is shown to design the gain matrix such that all the components tend to zero as the time tends to infinity. Let then one can construct a Lyapunov function candidate as follows:By taking the time derivative of along the trajectories of subsystems (6), one can deduce thatLet , then it can be found thatDue to , it can be shown that . Hence, it can be found by (16) thatOne can set thatSinceit can be derived by (18) thatIfthen one can show by (18) thatIn this case, since is invertible, one can find thatThus, if , then by (15), (17), and (23). Hence, swarm system (1) with a fixed communication topology is formation achievable by control protocol (2) with .
Finally, the influences of the limited budget are analyzed. By (2), one can find thatBy and , it can be deduced by (27) thatLet , then one hasByit can be derived from (27)–(29) thatSince is the maximum eigenvalue of the Laplacian matrix and , one can derive by (31) thatIfthen all the components tend to zero as the time tends to infinity, and tend to zero as the parameter the time tends to infinity; that is, and . Thus, it can be found by (27) and (30) thatDue toone can derive thatSinceit can be shown thatthat is,Since , from (34)–(36), it can be found thatwhere . It is assumed that where are disagreement, so there always exists some . Hence, it can be obtained thatIn this case, there exists a positive scalar such thatThe matrix has a simple zero eigenvalue and nonzero eigenvalues are positive, so it can be derived by (37) and (39) that can guarantee that , which means thatThus, the Proof of Theorem 1 is completed.
In Theorem 1, a limited-budget formation design criterion is proposed, which gives a design approach of the gain matrix of control protocol (2) to make swarm system (1) achieve limited-budget formation. If the gain matrix is given, then the following theorem shows a limited-budget formation analysis criterion, which can be obtained directly according to the Proof of Theorem 1.

Theorem 2. For any given , and , swarm system (1) with control protocol (2) achieves limited-budget formation if and there exists such that

There exist two difficulties to check the matrix inequalities in Theorem 2. The first one is that it is not scalable since the number of constraint conditions increases as the number of agents increases. The second one is that the term makes matrix inequalities not linear, which are difficult to be checked. By the convex feature of linear matrix inequalities and the Schur lemma in [34], the following theorem can be obtained on the basis of Theorem 2, which is scalable and linear and can be checked by the feasp solver in the Matlab’s LMI toolbox given in [35] and will bring no conservatism from Theorem 2 to Theorem 3.

Theorem 3. For any given , and , swarm system (1) with control protocol (2) achieves limited-budget formation if and there exists such that

Remark 3. Intuitively speaking, by designing piecewise continuous differentiable , arbitrarily required formation structures for swarm system (1) with control protocol (2) to maintain can be obtained. However, the achievability of a geometric structure is dependent on the structure property of a swarm system; that is, the achievability is closely related to the system matrix . If the formation functions are not time-varying, that is, are constant, then the condition is needed for swarm system (1) with control protocol (2) to achieve time-variant formation structures as shown in [15]. If the formation functions are time-varying, then Theorem 1 requires that the condition for swarm system (1) with control protocol (2) to achieve time-varying formation structures. The literature studies [25, 26] showed some similar feasible conditions for swarm systems to achieve time-varying formation structures, where feasible conditions are different from the one in Theorem 1 since the limited budget was not taken into consideration.

Remark 4. Because the whole energy budget must be limited in practical engineering swarm systems, it is critically important to analyze the impacts of the limited budget for swarm system (1) with control protocol (2) to achieve specific formation; that is, swarm system (1) with control protocol (2) can obtain and maintain specific geometric structures under the condition that the whole energy consumption is restricted. Note that the condition introduces the maximum energy supply into the gain matrix design, where the initial condition and the interaction relationship matrix are used. The initial condition is often available in practical applications and the interaction relationship matrix shows the impact of the topology. By the structure property of the transformation matrix , Theorem 1 proposes an approach to construct the interaction relationship between the limited budget and the matrix variable, which introduces the limited budget into the formation criterion. It should be pointed out that the interaction relationship matrix can be regarded as the Laplacian matrix of a complete graph with all edge weights being . In this case, the interaction relationship matrix has a simple zero eigenvalue and identical nonzero eigenvalues. Moreover, the sufficient conditions in Theorem 1 for swarm system (1) with control protocol (2) to achieve limited-budget formation are scalable since they are independent of the number of agents and the two eigenvalues and of the Laplacian matrix can be estimated by the approaches in [36, 37].
By the analysis of Theorem 1, the component can be applied to depict the formation center function as the time tends to infinity, which is associated with subsystem (5). The following theorem shows a method to give an explicit expression of the formation center function.

Theorem 4. If swarm system (1) with control protocol (2) achieves limited-budget formation, then the formation center function satisfies that where and

Proof. of Theorem 4. Due toit can be shown thatFurthermore, the dynamic response of subsystem (5) introduced by initial conditions is thatBy (47) and (48), one can see thatMoreover, the dynamic response of subsystem (5) introduced by the formation function is thatIn particular, the dynamic response of subsystem (5) introduced by the derivative of the formation function is thatThe component can be applied to depict the formation center function as the time tends to infinity, so from (50) and (51), one can show thatwhereThus, the Proof of Theorem 4 is completed.

Remark 5. Limited-budget formation motions of a swarm system with a fixed communication topology contain two parts. The first part is the macroscopic motion as a whole and the second part is the microcosmic motions between any two agents. Limited-budget formation is obtained if and only if the microcosmic motions between any two agents are asymptotically stable. If limited-budget formation is achieved, then all agents in a swarm system move along the trajectories of the formation center function, which can be used to depict the macroscopic motion as a whole. From Theorem 4, the formation center function involves two parts, which are associated with the average value of the time-varying formation functions and the average value of the initial conditions of cooperative states, which are called the formation function part and the cooperative state part, but it is not associated with the derivatives of the formation functions; that is, the variance of the formation structure does not impact the macroscopic motion as a whole. In particular, according to the explicit expression of the formation center function in Theorem 4, the limited budget of energy assumptions does not also influence the macroscopic motion as a whole.

4. Numerical Simulations

This section provides the numerical simulation to demonstrate the correctness of our main conclusions about limited-budget formation control for swarm systems with fixed communication topologies.

Consider a swarm system with five agents, where the dynamics of each agent is modeled by (1) with

The communication topology of this swarm system is depicted in Figure 1, where 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 conditions of the cooperative states are presented as follows:

The formation functions are chosen as follows:

One can find that the formation functions satisfy the feasibility conditions in Theorem 1. Let and . By the feasp solver in the Matlab’s LMI toolbox, it can be found from Theorem 1 that

In this case, the gain matrix is .

In Figure 2, the trajectories of for this swarm system are presented, where the curves depicted by blue star markers represent the trajectories of the formation center function. Figure 3 shows the state snapshots of five agents and the formation center at  s,  s,  s, and  s, where five agents are depicted by black circles, green plusses, blue asterisks, pink x-marks, and red hexagrams, orderly. By Figure 3, one can find that all the agents achieve a time-varying pentagon around the formation center at different times.

The formation center of this swam system counterclockwise rotates about 1.62 cycles as shown in Figure 4, where the initial state is denoted by a red pentacle and the final state is represented by a black pentacle.

Figure 5 shows that the practical energy consumption converges to a finite value less than . From Figures 25, one can see that this swam system achieves the desired pentagonal formation in the form of the circular movement and the limited budget is guaranteed.

5. Conclusions

For high-order linear swarm systems with fixed communication topologies, a new formation control protocol with an energy integral term was proposed by using the state errors and the formation function errors among neighboring agents, which can guarantee that the practical energy assumption is less than the whole energy budget. Furthermore, by an orthonormal transformation associated with the Laplacian matrix of the fixed communication topology, limited-budget formation control problems were converted into asymptotic stability ones, and sufficient conditions for swam systems to achieve limited-budget formation were proposed, where the relationship between the matrix variable and the limited budget was constructed on the basis of the specific structure feature of the orthonormal transformation matrix. In particular, those criteria are scalable and checkable since their constraints are independent of the number of agents and can be converted into linear matrix inequalities. Moreover, an explicit expression of the formation center function was presented, which contains two independent parts associated with the formation functions and the cooperative states, but it is independent of the variances of the formation functions.

The further works will focus on extending the current results to the swarm systems with heterogeneous dynamics, and the fixed topology in this paper will be changed to the switching connected topologies and the jointly switching topologies.

Data Availability

The data used to support the findings of the study are available in Section 4 of this paper.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Hongtao Dang and Le Wang conceptualized the study and prepared the original draft; Hongtao Dang investigated the data and was responsible for methodology and funding acquisition and involved in project administration; Yan Zhang validated the data; Jianye Yang performed formal analysis; Yan Zhang and Jianye Yang reviewed and edited the manuscript; Le Wang supervised the study. 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), also funded by the National Natural Science Foundation of China under Grants 61867005.