1. Introduction
In order to realize the common interests of some and even all individuals, the cluster activity is to unify the scattered individuals and social forces. They can adapt to different living conditions by living in groups. The purpose of cluster activities is to control, coordinate and restrict individuals and small groups living in groups. They manage the factors affecting the development of groups. The cluster activities ensure the survival and predation of groups to a certain extent. They can also coordinate with each other to maximize the interests of the group. At the same time, the cluster activities also greatly improve work efficiency.
Through the study and research of biological cluster activities, the swarm motion of bionic robots also develops. Usually, the cluster of robots is called a multi-robot system. The multi-robot system refers to a kind of complex system composed of multiple robots and their corresponding organization rules and information interaction protocols. It can complete specific tasks, such as multiple UAVs, multiple robots and so on. Among them, the organization rules determine the connection relationship between robots. The information interaction protocol is used to determine and update the state of robots. There are a large number of such examples in the real world [
1,
2,
3], as shown in
Figure 1. At the same time, the multi-robot cluster movement also began to develop. The cluster movement can not only improve the energy utilization, but also improve the efficiency of movement. There are also good solutions for the large-scale task requirements.
After more than 20 years of development, the multi-robot collaboration has shown a wide range of applications in the military, national defense, industry, life and other fields [
4]. For example, it has shown great advantages in specific applications such as the formation patrol, environmental exploration, assembly-line production and warehousing and transportation [
5,
6,
7]. Some specific application tasks include multi-robot cooperative map exploration, large object transportation, target capture and so on. They require the robots in the team to form a specific formation structure. According to the environmental constraints and the different task requirements, they can keep the formation unchanged or change the formation shape in the process of task execution. The common formations include the straight line, wedge, column, triangle, double row, etc. [
3], as shown in
Figure 2. Arrows between robots indicate the flow of information. The arrow in front of the robot R
1 indicates the direction of motion of the multi-robot system. The research work of this article will take the mobile robot as the main research object.
Therefore, the research boom of multi-robot formation control has been triggered. It has become one of the key technologies in the field of multi-robot cooperation and coordination. Taking the nonholonomic multi-mobile robot as the research object and the cooperative formation as the application background, the key technologies of multi-mobile robot motion formation control is studied in this article.
Motion control is one of the important research topics in the multi-robot systems. By establishing effective communication control strategies among robots, the multi-robot systems can cooperate to complete a complex task. The multi-robot systems realize the expected formation. At present, the commonly used robot formation control algorithms include the virtual structure method, based behavior method, artificial potential field method, leader-follower method [
8,
9,
10,
11]. The advantage of the leader-follower method is that the formation control structure is simple and easy to implement. In the formation, only the desired path or other behavior of the leader needs to be set, and then the follower can follow the leader with a predetermined position offset to realize the formation control. For this reason, the leader-follower method is widely used in practical engineering. Hence, the leader-follower method can be employed for the formation control.
The concept of leader-follower was first proposed and successfully applied to the formation control of mobile robots [
12,
13]. As the most commonly used formation control method at present, its basic idea is that all formation members are designated as the leader or follower. The leader controls the movement trend of the whole formation by navigating along the predetermined or temporarily set path. According to the distance and azimuth information relative to the leader, the follower follows the leader to realize formation control. The advantage of the leader-follower method is that the formation control structure is simple and easy to implement. In the formation, only the desired path or other behavior of the leader needs to be set, and then the follower can follow the leader with a predetermined position offset to realize the formation control. For the above reasons, the leader-follower method is widely used in practical engineering. The disadvantage of the method is that the formation system depends too much on the leader.
A new multi-robot formation control method is known as the leader-follower method. The stability of the algorithm was proved by the simulation results [
14]. The formation control problem was transformed into the problem of state consistency by the transformation of variables. A controller was designed to make a group of nonholonomic robots gradually converge to the desired formation along the specified trajectory [
15]. When the external input was limited, the navigation and following formation problem of nonholonomic multi-mobile robots was discussed. Under the nonholonomic constraints, the external distributed sensor method was applied. The finite-time convergence problem of multi-robot cooperative control was studied [
16,
17]. The motion consistency problem of complex multi-robots was studied, and the dynamic formation of multi-mobile robots was realized by designing a distributed formation control law [
18,
19]. The time-varying stability problem of planar multi-robots with nonholonomic constraints was studied [
20,
21]. Assuming that there were many external disturbances in the system. The finite-time containment control problem of second-order multi-robot systems with mismatched disturbances was studied [
21,
22]. The adaptive tracking control of omnidirectional mobile robots was analyzed under the conditions of state constraints and input saturation [
23]. Assuming that there were multiple leaders, the cluster containment control problem of the robot system was researched [
24]. The acquisition spacecraft operation of the space robot was analyzed. The collision avoidance compliance compound active disturbance rejection control problem of a robot was studied. Considering the post-capture unstable combined system, collision avoidance and compliant composite error active disturbance rejection control is proposed for the stabilization control. Numerical simulation verifies the effectiveness of the proposed collision avoidance and compliant control method [
25,
26].
Most of the above studies focus on the leader-follower method implementation and external distributed sensor fusion. The multi-robot cooperative formation needs to closely combine the theories and algorithm simulation. This article combines the nonholonomic constrained mobile robot, graph theory, leader-follow formation control and trajectory tracking. The error analysis of the nonholonomic multi-robot system is implemented. Based on information consistency, the formation realization of single formation and multi-combination formation is completed. Through the existing research results, in
Section 2, the communication topology matrixes are designed by using a consistency algorithm based on the graph theory. The formation control law of multi-robot cooperative motion based on the leader-follow model is constructed. In
Section 3, the mathematical model of a wheeled mobile robot under the nonholonomic constraints is established. Then, the leader-follower model is transformed into the local coordinate system error model through the global coordinate transformation. The Lyapunov function is used for the theoretical analysis. In
Section 4, by setting the reasonable communication protocol parameters, the simulation results are competed with another algorithm of leader-follower formation by using MATLAB software. Based on the leader-follower method, the effectiveness and feasibility of the formation control law are verified by compared to other formation control strategies. In the final section, the conclusions are presented.
2. Multi-Robot Cooperative Formation Control Based on Group Formation
2.1. Graph Theory Description
The control cooperative multi-robot dynamic system is a dynamic problem related to each other through a communication diagram, which shows the information flow between each node. The goal of cooperative control is to design control protocols for each node and to ensure the state synchronization behavior of all nodes in a specific sense. In a cooperative system, any control protocol must carry out the distributed control according to its specified graph topology. The control protocol of each node is only allowed to rely on the information about the node and its neighbor nodes in the graph. The communication restrictions imposed by the topology may seriously limit the work. This work is completed by the local distributed control protocol on each node.
In multi-robot cooperative formation, this article focuses on the behavior and interaction of dynamic systems. These dynamic systems are interconnected through the link of a communication network. The communication network is modeled as a graph with directed edges corresponding to the allowed information flow between systems. Each robot in the formation is modeled as a node in the graph. The most basic problem is how the graph topology interacts with the local feedback control protocol of nodes to generate the overall behavior of interconnected nodes. This problem can be solved by the graph theory. The following shows some basic graph theory concepts that are essential in the study of multi-robot cooperative formation.
The topological graph can be represented as G = (V, E), where V = {v1, …, vN} represents the set of N nodes, and E represents the set of E edges. Assuming that the graph is a simple directed graph, that is (vi, vj) ∉ E. At this time, there are no multiple edges between the same pair of nodes. The elements (vi, vj) in the set E represent the edges from node vj to node vi. In the formation, it represents that the information flows from node vj robot to node vi robot. The out degree of node vi robot represents the number of edges with it as the data sending end. The in degree represents the number of edges with it as the data receiving end.
The structure and properties of graphs can be studied by studying the properties of some matrices related to graphs, which is algebraic graph theory. After the weight
aij of the edge is given, the graph can be represented by an adjacency matrix or connected matrix
A = (
aij), where when (
vi,
vj) ∉
E,
aij = 1, otherwise
aij = 0. The weighted penetration of point
i at this time can be expressed as the sum of all values in row
i of matrix
A, as shown in Equation (1). Then
Then, the weighted output can be expressed as Equation (2). Then
Both in degree and out degree are local properties of the graph. In formation topology communication, the out degree can reflect the impact of the node on other robots in the formation, that is, the greater the out degree, the greater the impact of the node robot on the formation consistency.
The diagonal penetration matrix and graph Laplacian matrix are defined as
D =
diag di and
L =
D − A, respectively. Many properties of a graph can be studied as the properties of graph Laplacian matrix. It is of great significance in the study of graph theory of dynamic multi-node systems. In order to explain the relationship between the directed graph and adjacency matrix, between the diagonal matrix and Laplace matrix, a directed graph of a multi-robot system is assumed as shown in
Figure 3. So, its adjoint matrix
A, diagonal matrix
D and Laplacian matrix
L can be shown in Equation (3).
In
Figure 3, the red circle represents the leader robot. The blue circle represents the following robot. The number represents the serial number of the robot. The arrow represents the directed edge. The corresponding matrixes of the directed graph are shown as follows in Equation (3). Then
The characteristic structure of the graph Laplacian matrix L plays a key role in the dynamic system analysis on the graph. The eigenvalue of L is analyzed below, and the approximate form of the graph Laplacian matrix is defined as follows in Equation (4).
In Equation (4), the Jordan matrix
J and the transformation matrix
M are shown as Equation (5). Then
In Equation (5), the eigenvalue
λi and the right eigenvectors
of
M satisfy Equation (6). Then
where
I is the identity matrix. The left eigenvector
of matrix
M satisfies Equation (7).
Then
where
is the transpose vector of the eigenvector
.
Since the sum of all row elements of matrix L is 0, if and only if L has a unique eigenvector c such that Lc = 0 and its eigenvalue is 0. The rank of L is (N − 1). Graph G has a spanning tree, which plays an important role in forming the stability of formation.
2.2. Formation Control Principle
The distributed control structure in the leader-follower formation is adopted in this article. The formation control is realized by using the consistency algorithm given below on the basis of the previous graph theory. The consistency is defined as one or some states of all robots in a multi-robot system that tend to be consistent with the evolution of time. The leader-follower formation is a classical formation form. The main idea of the first-order formation algorithm based on consistency is to control the relative position and speed between the leader and the follower. When these quantities reach a certain stable state, the formation is realized, in which the leader robot plays a vital role. The first-order continuous system model of the formation system is Equation (8) as follows.
Then
where
xi and
ui represent the state quantity and input quantity of node
i respectively, and
xi,
ui∈
Rn,
n represent the dimension of state quantity. Ideally, consider the following forms of control input. Then
where
aij is the formation adjacency matrix element and
Ni is the neighbor set of member
i. By substituting Equation (9) into Equation (8) and introducing the global state vector
x = (
x1, …,
xn)
T ∈
Rn, the global dynamic relationship can be obtained as follows in Equation (10). Then
Based on Equation (10), it can be concluded that the closed-loop dynamic characteristics of the formation depend on the matrix L. In addition, if and only if the topological graph G of the formation has a spanning tree, −L eigenvalues are located in the left half plane of the complex plane. Hence, Equation (10) can ensure the consistency of the system. The final state value is shown as follows.
In discrete time, Equation (8) can be reduced to Equation (12). Then
In this subject, it is assumed that a scene is a scene in which the leader leads other robots to complete a certain distance and move from the starting point to a target point. Ideally, the algorithm of the follower in the formation adopts the following form. Then
Among
ε > 0,
rij represents the relative position between robot
i and robot
j, which is related to the distance between them, while the control algorithm of leader
N is Equation (14). Then
where
rNj represents the relative position between the leader and other robots,
m and
k are constants,
D(
k) represents the distance from the leader to the target point at time
k, and the speed of the leader is affected by the adjacent nodes and the distance from the leader to the target point.
2.3. Implementation of Group Formation Control Algorithm Based on Consistency for Nonholonomic Systems
The basic idea of the consistency algorithm is that each robot uses the robot network to transmit information. It designs an appropriate distributed control algorithm. It couples the robot dynamics with the robot communication topology to form a complex system, so as to realize the continuous or synchronization of the state.
The goal of this section is to design a unified grouping consistent formation (GCF) scheme for a system that can be divided into
k groups and composed of
N incomplete mobile multi-robots to meet the large-scale multi-robot system formation.
V = {
v1,
…,
vN} represents a collection of
N nodes,
Vk represents the
k-th group, and
Vk∈
V, this article adopts the multi-group formation mode to realize a large-scale multi-mobile robot system, as shown in the following figure. R
kL represents the leading robot of the
k-th group, and R
k-s represents the
s-th follower robot of the
k-th group, as shown in
Figure 4. While the leader robot R
kL forms a formation with the inside of the group, it also forms a formation with the leader robots of other groups. The leader robot only has a communication topology diagram with the robots in the group and the leader robots in each group. The group formation greatly increases the size of the formation, and also reduces the amount of information processing of the navigator in the multi-robot system. It is of great significance to a large number of robot systems.
The corresponding communication topology matrix is shown as follows in Equation (15). Then
In order to verify the formation consistency algorithm, five robots are used to form an experimental robot formation in the form of leader-follower to find a stable communication topology matrix
A′ as shown in Equation (16). Then
After using the communication topology matrix to make a simulation experiment in
Section 4, the formation and its communication topology are shown in
Figure 5. The arrow indicates the direction of information flow. In
Figure 5, it can be seen that the outgoing degree of node 5 is the largest. So, node 5 is the key node. That is the leader robot. The outgoing degree of node 1 is 0, which has the least impact on the formation. Both in and out degrees are local properties of the graph. In formation topology communication, their out and in degrees can reflect the influence of nodes on other robots in the formation. The larger the exit degree, the greater the influence of the node robot on the formation consistency. When the penetration is larger, the influence of the node robot on the formation consistency is smaller.
3. Mathematical Model of Formation Operation Control
The basic idea of the leader-follower method is to select a leader of robots and set its trajectory. The following robot tracks the trajectory of the leader robot in real-time. The trajectory of the whole formation is determined by the leader robot. The formation shape is maintained by the relative position relationship between the tracking robot and the leader. In the motion process of a multi-robot system based on the leader-follower model, in order to avoid collision between robots during tracking, a virtual robot generated by the leader robot and completely obtains the state information of the leader robot is introduced to calculate the state error between the follower robot and the virtual robot and make it converge to zero. Based on the above analysis, the formation problem of multi-robot can be transformed into a trajectory tracking problem of the following robot to leader robot.
3.1. Robot System Model
The nonholonomic constrained robots are mainly wheeled mobile structures, including wheeled, wheel leg composite, crawler and so on. In order to facilitate the theoretical analysis of the algorithm, this article takes the mathematical model of the four-wheel mobile robot as an example, as shown in
Figure 6. The front and rear wheelbase of the robot is
L. The distance between the midpoint of the wheel and the midpoint of the robot system is
d, as shown in Equation (17). The motion direction angle of the robot is
θ. It is assumed that the linear velocity and angular velocity of the robot motion are
ν and
ω. The linear speed of the wheel is
u. Then
The influence of sliding can be ignored in this article, so the mobile robot belongs to a pure rolling sliding free robot. For the robot with pure rolling and no sliding, it is necessary to meet the nonholonomic constraints, as shown in Equation (18). Then
By expressing all state quantities of the mobile robot, its mathematical model can be obtained, as shown in Equation (19). Then
where
v and
w are the parameters that can be directly controlled,
x,
y and
θ are the state quantities of the system. The number of parameters that can be controlled is less than the number of state qua tities of the system. Therefore, the system is an underactuated system, also known as a nonholonomic constraint system.
3.2. Mathematical Model of Leader
According to the conditions to be met by the robot model,
Section 3.2 gives the state quantities of the leader robot as
xL,
yL and
θL. Then the ideal state quantities of the follower robot can be expressed as Equation (20).
Then
where the ideal state quantities of the follower robot are as
,
and
.
is the angle between the ideal state of the follower robot and the leader robot in the coordinate system of the leader robot.
is the distance between the ideal state of the follower robot and the leader robot.
The leader formation control algorithm can be introduced virtual leader robots as expectation analysis. In a nonholonomic system, both the virtual leader robot and the navigation robot meet the nonholonomic constraints. The mathematical model of the virtual leader robot can be described as Equation (21). Then
where
and
are the parameter of the virtual leader robot.
xv,
yv and
θv are the state quantities of the virtual leader robot.
,
and
are the first derivatives of
xv,
yv and
θv respectively.
According to Equation (19), the mathematical model of the following robot is shown in Equation (22). Then
where
and
are the parameter of the virtual leader robot.
xf,
yf and
θf are the state quantities of the virtual leader robot.
,
and
are the first derivatives of
xf,
yf and
θf, respectively.
3.3. Leader-Follower System Model
The simple diagram of the leader-follower system model is shown in
Figure 7. According to Equations (20)–(22), the state error between leader formation multi-robot systems in this article is Equation (23).
Then
where
,
and
are the state errors of the system. The ideal state quantities of the follower robot are defined as
,
and
.
xi,
yi and
θi are the state quantities of the
i-th robot.
The virtual robot proposed in the previous section is used as error analysis. Equation (24) can be obtained from the coordinate relative position relationship between the virtual leader robot and the following robot. Then
Sorting out the state error, deriving Equation (24), and substituting Equations (21) and (20), the following Equation (25) can be obtained. Then
The state error after coordinate transformation can be described as Equation (26). Then
The input of the following robot
Rf are
ν and
ω.
ν and
ω can be adjusted to satisfy the following Equation (27). Then
In order to make the tracking error converge to zero, the multi-robot formation control law based on the leader-follower method can be designed as Equation (28). Then
where,
and
are the positive real number, and
+
= 1.
3.4. Convergence Proof
In the field of automatic control, the Lyapunov stability can be used to describe the stability of the dynamic system. According to
Section 3.1,
Section 3.2 and
Section 3.3, the Lyapunov function selected in this article is shown in Equation (29). Then
Obviously, and are greater than or equal to 0 respectively. is also greater than or equal to 0. So, the function positives definite.
According to Equation (29), for function after derivation, the following Equation (30) is obtained.
By introducing
and Equation (26) into Equation (30), Equation (31) can be obtained. Then
The function is known from Equation (31). It negatives definite.
Lyapunov stability main criterion is often used to prove the stability of the system. Let the state equation of the system be:
, if there is a scalar function
with continuous first derivative and it satisfies that
is positive definite and
is negative definite. Then, the system is asymptotically stable in the equilibrium state [
20].
According to Lyapunov stability main criterion, the function positive definite, and is negative definite, so the above control law makes the error system in the equilibrium state asymptotic stability.
5. Conclusions and Future Work
Based on the theory of automatic control and the nonholonomic system, this article studies the formation problem of nonholonomic constrained multi-robot tracking the leader’s movement while maintaining the ideal formation. The formation communication protocol of multi-robot cooperative movement is established. Through the mathematical analysis of nonholonomic constraints, the motion models and motion system error models of the leader robot and follower robot are constructed respectively. Then the convergence of formation control strategy in this article is proved theoretically by the Lyapunov function. The feasibility of formation control of the multi-robot system with the nonholonomic constraints is verified by MATLAB simulation. By comparing with the simulation results of other navigator formation control algorithms, the effectiveness of the GCF algorithm is verified.
In order to illustrate the effectiveness of the group consistency control designed in this article, the comparison of the two with and without group strategy are implemented. Then, two simulation results are obtained. The comparison of the two simulation results shows that the convergence time of formation control using the group strategy is less than that of the algorithm without the group strategy under the same formation. In order to better show the feasibility of the algorithm proposed, the comparison is performed with PFC algorithms. Through the comparison of the two algorithms, the convergence time of the algorithm in this article can be 25% less than that of PFC. Compared with PFC, the number of robots is poor. However, the overall comparison of the formation control algorithm can be acceptable in theoretical analysis and simulation. Hence, the formation control algorithm has effectiveness and feasibility. Because there are many uncertain conditions in the real robot systems, such as external noise and nonlinear interference, the formation control problem of nonholonomic constrained multi-mobile robot systems with external interference would be studied in the future.