Abstract

This paper focuses on studying containment control problem with switching communication graphs of continuous-time heterogeneous multiagent systems where the control inputs are constrained in a nonconvex set. A nonlinear projection algorithm is proposed to address the problem. We discuss the stability and containment control of the system with switching topologies and nonconvex control input constraints under three different conditions. It is shown that all agents converge to the convex hull of the given leaders ultimately while staying in the nonconvex set under the premise that at least one directed path from leaders to the agents exists in each bounded time interval. Finally, the validity of the results obtained in this paper is verified by simulation.

1. Introduction

With the in-depth research of autonomous control of multiagent systems in recent years, many cooperative control problems of multiagent systems have become the focus of many fields, such as control theory, biology, robotic systems, and spacecraft systems [15]. The study and exploration of multiagent systems provided a unified framework and theoretical basis for various practical problems such as unmanned aerial vehicles, formation aircrafts, multiple robots, and other practical applications. A Lyapunov-based approach was proposed to address the consensus problem in [6]. Lin et al. [7] emphasized on the consensus control problem in consideration of nonconvex control input constraints. The consensus control problem of multiagent systems with time delays was studied in consideration of external interference in [8]. And the work [9] was centered on the containment control problem. Lin et al. [10] studied distributed optimization problems for continuous-time and discrete-time multiagent systems with different constraints.

As an important branch of control theory, there were many articles focusing on the study of containment control [1120]. Ji et al. and Li et al. [11, 13] solved the containment control problem with fixed topology. The output formation-containment problem of heterogeneous system was investigated in [14, 15]. For double-integrator multiagent systems, containment control problem with fixed communication topology and position measurements was addressed in [16, 17]. Two distributed algorithms were proposed for containment control in the case of only using absolute position measurements and relative position measurements, respectively. However, the ‘sign’ function was employed in [17], which may lead to the chattering phenomenon. Cheng et al. [18] avoided employing the ‘sign’ function, and both disturbance and measurement noise were also taken into consideration. Notarstefano et al. and Zhang et al. [19, 20] investigated containment control with switching topologies. Cao et al. [12] studied the solution of containment control in consideration of switching topologies and fixed simultaneously.

In many practical systems, the control inputs were generally constrained in a convex or a nonconvex hull, while most of existing works studied the multiagent problem without considering the constraint of control inputs. The authors of [2124] studied the consensus problem with position constraints on the basis of the property of stochastic matrices. Nevertheless, these approaches could not be applied to the realization of containment control with the nonconvex control inputs’ constraints. The containment control problem with input saturation of the second-order agent system was studied in [25, 26]. However, control input of each agent was supposed to be in a hypercube. In reality, on account of the physical limitation, the control inputs were often constrained in a convex or a nonconvex region. For instance, the maximum driving forces of quadrotors, which formed a nonconvex region, were in the direction of the diagonal axis. A multiagent system model and a projection consensus algorithm were introduced in [23], and it placed emphasis on the effects of control input constraints. The algorithm was executed locally by each agent and its relationship with the alternate projection method was discussed in [23]. Yang et al. and Lin et al. [27, 28] took nonconvex velocity and control inputs’ constraints into consideration, but what they emphasized on was the consensus control problem which meant all the followers had to reach a consensus.

In [29], all agents were assumed to be in the form of the second-order dynamics. In practical applications, the multiagent systems might contain different kinds of agents. Our focus of this paper is to study the containment control for heterogeneous multiagent systems with nonconvex input constraints. Since the agents in this paper have different dynamics, the analysis for the case with all identical agents in [29] cannot be applied to this paper directly. In this study, we expand the results of [29] and mainly focus on the containment control problem of continuous-time heterogeneous multiagent systems, given the nonconvex control input constraints and switching communication graphs. Li et al. [30] studied the containment control of heterogeneous multiagent systems. However, they did not consider the constraints of control inputs and heterogeneous multiagent systems at the same time. To analyze the stability and convergence of the system, a nonlinear projection algorithm is proposed to address the problem. Then, a model transformation is introduced and estimates the distance from the followers to the nonconvex hull by using the Lyapunov function we construct. We prove that the distance decreases and ultimately all followers converge to the nonconvex hull.

1.1. Notations

Assume that represents the set of -dimensional real matrix. is the transposing matrix of . is the Euclidean norm of . denotes the projection of a vector onto a closed convex set , i.e., .

2. Preliminaries

2.1. Graph Theory

is a directed weighted graph representing the communication graphs among the agents, where is the set of node representing the agents and is the set of directed edges representing the communication between the agents. is a weighted adjacent matrix. Laplacian is a very important matrix in graph theory, which is defined as , where

Then, the Laplacian matrix can be expressed as follows:

For heterogeneous multiagent systems consisted of first-order and second-order agents, the adjacent agents of second-order agent can be denoted as . Similarly, the adjacent agents of first-order agent can be denoted as and represent the second-order and the first-order agents, respectively. Partition the matrix A and D aswhere and represents the adjacent relationship of the second-order and the first-order agents. and represents the adjacent relationship of the first-order and the second-order agents. Then, the Laplacian matrix can be expressed as follows:where , , and and are the Laplacian matrix of second-order agents and first-order agents, respectively.

2.2. Nonconvex Constraints

The control inputs are generally subject to nonconvex constraints in many practical cases. Hence, we have the following assumption.

Assumption 1. (see [7]). is a nonempty bounded closed set, and . is a nonconvex constraints’ operator and is defined as follows:In this study, the operator is used to convert to which has the alignment with the direction of and satisfies and , for all . means that the value cannot be arbitrarily large. We do not have any requirement for to be convex or nonconvex. What we require is that the distance from any point outside to the origin is lower bounded by a positive constant.

3. System Modeling

Consider a continuous-time multiagent system consisted of leaders, first-order followers, and second-order followers. The first-order agents have the following dynamic equation:

The dynamic equations of the second-order agents can be denoted as follows:where represent the second-order and the first-order agents, respectively, and , , and represent the position information, speed information, and control input of the second-order agents separately. and represent the position information and control input of the first-order agents separately. The set of all followers is expressed by . We assume all leaders are stationary and the position states are denoted by .

In this study, if all the followers can converge into the convex set of leaders under the control protocol for any initial value given, that is, , the control protocol can realize the containment control.

4. Model Transformation

To analyze the containment control of continuous-time heterogeneous multiagent system, a control protocol is proposed in this study as follows:where and is speed decay factor of followers. and are the weight of edge at time . And, we assume the weight of edge is always positive and lower bounded by a constant . If follower can get the information from one or more agent directly, then , otherwise, . Assume that is the convex set of the leaders. The Laplacian matrix of the heterogeneous multiagent system can be denoted by the following matrix:

Then, algorithm (8) can be converted into

The following model transformation is introduced for further derivation:when . For the first-order agents, (11) can be denoted as follows:where . For the second-order agents, (11) can be transformed as follows:where . From the above definition, it is easy to know the range of , that is, . We assume that is lower bounded by a positive constant , for all . According to the above definition, system (8) can be converted into the following form.

For the first-order agents,

For the second-order agents,

Denote , where is a positive constant. Then, (15) can be expressed in the following form:

5. Main Result

Prior to the main theorems, we need to give the definition of switching communication graphs. Denote an infinite sequence of nonempty bounded continuous intervals as , and , where is a positive constant. Divide the above interval into a series of subintervals which are represented as with and . We assume there exists a constant that makes . The communication topologies change at and do not change in the subinterval .

Assumption 2. Suppose that there exist at least one path between leaders and followers in the interval for every agent. In other words, each follower can receive information from leaders directly or indirectly in the every interval .

Assumption 3 (see [29]). and for some constant .

Lemma 1 (see [29]). Suppose that is a nonempty closed set. For any vector , and , we haveTo analyze the stability and convergence of the system, we construct a Lyapunov function as follows:where , .

Theorem 1. If and , under assumption 1 and 2, is a nonincreasing function relative to time and , namely, the limit of exits. And the following inequalities hold:

Proof. For the second-order agents, we can deducewhere is time increment and . is a positive constant. Then, it is evident that , , and . According to Lemma 1, (22) is converted toAnalogously, we can do the same conversion for :To facilitate later calculations, we have simplified the upper formula to the following form:whereFor all , we haveWhen , we obtainAnd . Then, (24) can be converted into the following forms:Combined with the following conditions, , , and ; it follows thatSimilarly, for the first-order agents, we can deducewhere is time increment and . According to (31) and Lemma 1, we haveTo facilitate later calculations, we have simplified the upper formula to the following form:where , , and . We can deduce from (33)Similar to the second-order agents, according to the condition and , we have . Then,It is evident that every is a convex hull that is composed of . From (23), (30), and (35), we can obtain thatThus, can be deduced which indicate that is a nonincreasing function relative to time . It is obvious that . Thus, it is easy to know that the limit of exists.
For the second-order agents, according to the definition of derivative and (23), we haveThat is,Similarly, according to (30), we haveThat is,whereFor the first-order agents, similarly, according to (35), we haveThat is,where .

Theorem 2. Under assumption 1 and 2, all the followers can move into the convex hull spanned by static leaders with algorithm (3), i.e., .

Proof. We will discuss the convergence of the heterogeneous multiagent system with switching graphs and nonconvex control input constraints at intervals in three cases.Case 1: suppose that there exists a agent can receive the information from leaders at time . At this point, . First of all, we assume is a second-order agent. From (40), we haveAccording to the definition of calculus, for all , we haveFor all , thus, we havewhere It is obvious that . According to (38) and (46), we obtainFor all , we havewhere . Obviously, .
If is a first-order agent, we haveThen,where . Obviously, .Case 2: If there exists a follower such that , for the second-order agents, from (23), we obtainFor all , we havewhere . Obviously, . Then, we can obtain thatAndwhereHence,whereFor the first-second agents, similarly, we haveThen,where .Case 3: suppose that there exists a follower such that , , where . According to the above assumption, for the second-order agents, we haveHence,For all , we havewhere , . Furthermore,where , . And similarly, for the first-order agents, there exists that makesUnder assumption 2, we know that there exists at least one follower which can receive the information directly from leaders. From the analysis of the above cases, in the time interval , we have the following conclusions for all agents:where the constant . And there exists another follower which can receive the information from leaders or in the interval . Hence,where the constant . Then, we can draw the following conclusions:where . Hence, there exists a constant that makes . Then, we can know the limit of exists and . Therefore, we can draw the conclusion that the control protocol can realize the containment control.

6. Simulation

Consider a continuous-time heterogeneous multiagent system with 6 followers and 4 leaders. According to the topology assumptions, we design four different communication topologies which are shown in Figure 1. And they switch every 1.5s. The control inputs are limited in the following set: . Suppose the weight of every edge is 0.6. Initial position values of followers are , , , , , and . Initial position values of leaders are , , , and . Initial velocity of agents are , , , , and . We suppose each parameter is taken as . Figure 2 shows the trajectory of followers. As you can see in Figures 2 and 3, all the followers finally converge to the set composed of multiple stationary leaders while the control input of each followers remains in their corresponding constraint set. The simulation result is consistent with theorem.

7. Conclusion

This paper focuses on distributed containment control with nonconvex control input constraints of heterogeneous multiagent systems. For this system, a protocol with nonconvex control input is proposed. Based on the Lyapunov function and the definition of the derivative, we can theoretically infer the distance between each follower and the set consisted of leaders in each time interval. Then, we can obtain that each follower is driven to converge to a convex set composed of multiple stationary leaders. Future work could be directed to the containment control of the system which consisted of more different order multiagent.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Key R&D Program of China (no. 2021YFB3900300).