# Bearing-Only Formation Control for Cascade Multirobots

Authors: Qing Han; Hongan Yang; Hao Lang
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

A new formation control method is proposed, which is used to queue multirobots in a single-direction cascade structure. In the cascade formation, each robot is a follower for the previous robot and a leader for the next robot, and the robots in the middle act as both leader and follower. The follower robot can only observe the bearing information of the leader robot. The observability of the cascade leader-follower formation is studied, which shows that the bearing-only observation meets the observability conditions required for the nonlinear system. Based on the bearing-only observations, the unscented Kalman filter (UKF) is employed for the state estimation of the leader and the follower robots at all levels, which enables the real-time movement control of the follower robots via the input-output feedback control. Simulation results demonstrate that the proposed approach can efficiently control the formation of multirobots as desired.


---

## Body

## 1. Introduction

Multirobots formation control has long been a hot points in both academic research and industrial application, particularly, in the areas such as underwater or outer space exploration, shop-floor transportation, guarding, escorting, and patrolling missions [1, 2]. The advantages of multirobots system over a single robot include greater flexibility, adaptability, and robustness [3–5]. However, multirobots’ formation control is challenging, especially when the observation information is poor and the system is highly nonlinear. A variety of formation control methods have been proposed, such as virtual structure approach [6, 7], behavior-based approach [8], leader-follower approach [9–14], artificial potential approach [15, 16], and graph theory approach [17–19]. Among them, the leader-follower approach has been widely used owing to the simplicity, scalability, and reliability. In this case, the leader follows a predefined trajectory, while the other robots (called follower robots) are keeping the position and direction with a certain distance to the leader [9, 13, 20–22]. A leader-follower formation control is designed in [9] for the follower robot to track the leader robot in the desired separation and bearing angle. In [10], a robust control technique is developed that can effectively handle the unknown parameters and uncertainties in the system. In [11], a receding-horizon leader-follower control framework is presented to solve the formation problem of multirobots with a rapid error convergence rate. To maintain the desired leader-follower relationship, a separation-bearing-orientation scheme (SBOS) for two-robot formation and a separation-separation-orientation scheme (SSOS) for three-robot formation are proposed. In [12], a framework for controlling groups of autonomous mobile robots to achieve predetermined formation based on a leader-follower approach is presented. In [13], the follower position is not rigidly fixed with respect to the leader but varies in proper circle arcs which is centered in the leader reference frame. The formation is implemented using speed measurements from optical encoders and distance measurements from image processing in [14]. In the existing leader-follower approaches, most of them need at least distance-angle information or more information.In real situations, the available observations of the robots might be only bearing observations, which will make a great challenge to the formation control. Until now, only a few special cases regarding bearing-only formation control problems have been solved. In [23], a bearing-only controller which can stabilize a group of mobile robots into a balanced circular formation is presented. In [24–26], bearing-only control laws that guarantee global stability are only applicable to formations of three or four robots. In [27], the bearing rigidity is employed for mobile formations, so that bearings can be used for shape control in mobile formation. In [28], a distributed control law to stabilize bearing-constrained formations and the concept of parallel rigidity are proposed. In [29], the decentralised formation control of multiple robots in the plane when each robot can only measure the local bearings of their neighbours by using bearing-only sensors is studied, and the target formation is locally finite-time stable with collision avoidance guaranteed. Distributed control of multirobot formation with angle constraints using bearing-only observations is investigated in [30]. In [20], observability conditions for position tracking by using bearing-only observations are established, and the localization problem is studied using a new observability condition valid for general nonlinear systems and based on the extended output Jacobian. The bearing-only formation control in [23–27] has special requirements for the number of robots or the form of the formation. A distributed control law in [28–30] is used to stabilize formations. In order to utilize the feedback control law, an off-the-axis point in [20–22] has to be additionally constructed which is a handling point lying on the robot’s axis of orientation, a distance d from the centre.In this paper, we focus no special requirements for the form of the formation and apply situations ofn robots. The contribution is threefold. Firstly, we do not need to construct the off-the-axis points when the input-output feedback control law is used. Secondly, the observability of the cascade leader-follower formation is studied based on the rank of the system observability matrix. This allows us to identify the robot motions that preserve the observability. Finally, based on the bearing-only observations, the unscented Kalman filter (UKF) is employed for the state estimation of the leader and the follower robots at all levels, which enables the real-time and stable movement control of the follower robots via the input-output feedback control.The rest of the paper is organized as follows. In Section2, observability and controllability of the nonlinear system are presented. Section 3 presents the UKF algorithm and the input-output feedback control law. Simulation results are given in Section 4. In Section 5, we offer our conclusions.


## 2. Modeling and Observability Analysis

### 2.1. Modeling

Consider the leader-follower setup in Figure1, where R1 is the leader of R2 and R2 is the leader of R3. The control inputs of R1, R2, and R3 are the linear and angular velocities ν
1
ω
1, ν
2
ω
2, and ν
3
ω
3, respectively. ρ
i is the distance from the centroid of the leader robot to the centroid of the follower robot. φ
i is the view-angle from the y-axis of the follower robot to the centroid of the leader robot. θ
i and θ
j are the orientations of the leader robot and the follower robot with respect to the world frame W, respectively, while α
i is the relative orientation between the leader robot and the follower robot; that is, α
i
≜
θ
i
-
θ
j.Figure 1
Leader-follower coordinate representation.With reference to Figure1, the kinematic model of a cascade leader-follower robots’ formation can be expressed as (1)
S
n
:
s
˙
=
f
s
,
u
=
F
s
U
y
=
h
s
=
h
1
T
s
,
h
2
T
s
T
,where state vector s
≜
s
1
T
,
s
2
T
T, s
1
≜
ρ
1
φ
1
α
1
T, s
2
≜
ρ
2
φ
2
α
2
T, input vector U
≜
ν
1
ω
1
ν
2
ω
2
ν
3
ω
3
T, and output vector y
=
h
s
=
h
1
T
s
,
h
2
T
s
T,  h
1
≜
φ
1
α
1
T
,
h
2
≜
φ
2
α
2
T, and where δ
i
=
φ
i
+
α
i
i
=
1,2,(2)
F
s
=
cos
⁡
δ
1
0
-
cos
⁡
φ
1
0
0
0
-
sin
⁡
δ
1
ρ
1
0
sin
⁡
φ
1
ρ
1
-
1
0
0
0
1
0
-
1
0
0
0
0
cos
⁡
δ
2
0
-
cos
⁡
φ
2
0
0
0
-
sin
⁡
δ
2
ρ
2
0
sin
⁡
φ
2
ρ
2
-
1
0
0
0
1
0
-
1
,the detailed derivation of F
(
s
) is given in Appendix A.The kinematic model of multistage cascade leader-follower robots’ formation can be readily retrieved as an extension of (1):(3)
S
n
:
⁡
s
˙
=
f
s
,
u
=
f
1
T
s
1
,
u
1
,
…
,
f
n
+
1
T
s
n
,
u
n
+
1
T
y
=
h
s
=
h
1
T
s
1
,
…
,
h
n
T
s
n
T
,where s
=
s
1
T
,
s
2
T
,
…
,
s
n
T
T
∈
R
3
n is the state vector, U
=
u
1
T
,
u
2
T
,
…
,
u
n
+
1
T
T
∈
R
2
(
n
+
1
) is the control input vector, h
:
R
2
n
↦
R
2
n, and f
:
R
3
(
n
+
1
)
×
Λ
↦
R
3
(
n
+
1
).


### 2.2. Nonlinear Observability of Cascade Multirobots

The functionf can be separated into a summation of independent functions in the special case, and each one is excited by a different component of the control input vector; (1) can be restated as (4)
S
n
:
s
˙
=
f
s
,
u
=
f
ν
1
ν
1
+
f
w
1
ω
1
+
f
ν
2
ν
2
+
f
w
2
ω
2
+
f
ν
3
ν
3
+
f
w
3
ω
3
y
=
h
s
=
h
1
T
s
,
h
2
T
s
T
=
φ
1
α
1
φ
2
α
2
T
.Accordingly, (3) can be also restated as(5)
S
n
:
⁡
s
˙
=
f
s
,
u
=
f
ν
1
ν
1
+
f
w
1
ω
1
+
⋯
+
f
ν
n
+
1
ν
n
+
1
+
f
w
n
+
1
ω
n
+
1
y
=
h
s
=
h
1
T
s
1
,
…
,
h
n
T
s
n
T
.In order to analyse the observability of the entire leader-follower system, the nonlinear observability rank criterion which is developed in [31] is used. The observability matrix of (5) is defined as the matrix with rows,(6)
M
=
∇
L
f
ν
i
,
…
,
f
ν
j
,
f
ω
i
,
…
,
f
ω
j
p
h
k
s
, where i
,
k
=
1
,
…
,
n, j
=
2
,
…
,
n
+
1, p
∈
N, and ∇ represents the gradient operator.Lemma 1.
SystemS
n is locally weakly observable if its observability matrix M which is constructed by the row vector has full rank; for example, in our case r
a
n
k
(
M
)
=
3
n
,
n refers to the number of follower robots.In (4),(7)
f
ν
1
=
cos
⁡
δ
1
-
sin
⁡
δ
1
ρ
1
0
0
0
0
T
,
f
ω
1
=
0
0
1
0
0
0
T
,
f
ν
2
=
-
cos
⁡
φ
1
sin
⁡
φ
1
ρ
1
0
cos
⁡
δ
2
-
sin
⁡
δ
2
ρ
2
0
T
,
f
ω
2
=
0
-
1
-
1
0
0
1
T
,
f
ν
3
=
0
0
0
-
cos
⁡
φ
2
sin
⁡
φ
2
ρ
2
0
T
,
f
ω
3
=
0
0
0
0
-
1
-
1
T
.The necessary Lie derivatives ofh
(
s
) and their gradients are computed; then observability matrix M is obtained.Zeroth-order Lie derivative is(8)
L
0
h
s
=
h
s
=
φ
1
α
1
φ
2
α
2
Tand its gradient is(9)
∇
L
0
h
s
=
∇
L
0
φ
1
α
1
φ
2
α
2
=
∂
φ
1
∂
ρ
1
∂
φ
1
∂
φ
1
∂
φ
1
∂
α
1
∂
φ
1
∂
ρ
2
∂
φ
1
∂
φ
2
∂
φ
1
∂
α
2
∂
α
1
∂
ρ
1
∂
α
1
∂
φ
1
∂
α
1
∂
α
1
∂
α
1
∂
ρ
2
∂
α
1
∂
φ
2
∂
α
1
∂
α
2
∂
φ
2
∂
ρ
1
∂
φ
2
∂
φ
1
∂
φ
2
∂
α
1
∂
φ
2
∂
ρ
2
∂
φ
2
∂
φ
2
∂
φ
2
∂
α
2
∂
α
2
∂
ρ
1
∂
α
2
∂
φ
1
∂
α
2
∂
α
1
∂
α
2
∂
ρ
2
∂
α
2
∂
φ
2
∂
α
2
∂
α
2
=
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
1
.The first-order Lie derivatives of the functionh
(
s
) with respect to f
ν and f
ω are defined as (10)
L
f
ν
1
1
h
1
s
=
∇
L
0
h
1
s
·
f
ν
1
=
-
sin
⁡
δ
1
ρ
1
0
T
,
L
f
ω
1
1
h
1
s
=
∇
L
0
h
1
s
·
f
ω
1
=
0
1
T
,
L
f
ν
2
1
h
1
s
=
∇
L
0
h
1
s
·
f
ν
2
=
sin
⁡
φ
1
ρ
1
0
T
,
L
f
ω
2
1
h
1
s
=
∇
L
0
h
1
s
·
f
ω
2
=
-
1
-
1
T
,
L
f
ν
2
1
h
2
s
=
∇
L
0
h
2
s
·
f
ν
2
=
-
sin
⁡
δ
2
ρ
2
0
T
,
L
f
ω
2
1
h
2
s
=
∇
L
0
h
2
s
·
f
ω
2
=
0
1
T
,
L
f
ν
3
1
h
2
s
=
∇
L
0
h
2
s
·
f
ν
3
=
sin
⁡
φ
2
ρ
2
0
T
,
L
f
ω
3
1
h
2
s
=
∇
L
0
h
2
s
·
f
ω
3
=
-
1
-
1
T
,where denotes the vector inner product, with their gradients given by(11)
∇
L
f
ν
1
1
h
1
s
=
s
i
n
δ
1
ρ
1
2
-
c
o
s
δ
1
ρ
1
-
c
o
s
δ
1
ρ
1
0
0
0
0
0
0
0
0
0
,
∇
L
f
ω
1
1
h
1
s
=
0
0
0
0
0
0
0
0
0
0
0
0
,
∇
L
f
ν
2
1
h
1
s
=
-
s
i
n
φ
1
ρ
1
2
c
o
s
φ
1
ρ
1
0
0
0
0
0
0
0
0
0
0
,
∇
L
f
ω
2
1
h
1
s
=
0
0
0
0
0
0
0
0
0
0
0
0
,
∇
L
f
ν
2
1
h
2
s
=
0
0
0
s
i
n
δ
2
ρ
2
2
-
c
o
s
δ
2
ρ
2
-
c
o
s
δ
2
ρ
2
0
0
0
0
0
0
,
∇
L
f
ω
2
1
h
2
s
=
0
0
0
0
0
0
0
0
0
0
0
0
,
∇
L
f
ν
3
1
h
2
s
=
0
0
0
-
s
i
n
φ
2
ρ
2
2
c
o
s
φ
2
ρ
2
0
0
0
0
0
0
0
,
∇
L
f
ω
3
1
h
2
s
=
0
0
0
0
0
0
0
0
0
0
0
0
.According to the form of (6), row vector matrix M of the cascade leader-follower formation described in Figure 1 can be expressed using gradients of Lie derivatives up to the first order as (12).Theorem 2.
The rank ofM given by (12) is six if(1)
ν
i
>
0
,
ν
j
>
0, where i
=
1
,
…
,
n, j
=
2
,
…
,
n
+
1;
(2)
φ
i
≠
0, i
=
1
,
…
,
n; that is, the follower robot at all levels, which is measuring the bearing, does not move along the line joining the leader robot and the follower robot;
(3)
the leader robot and follower robot at all levels can not do parallel linear motion.Proof.
Meeting the three prerequisite requirements above,M matrix in (12), by means of a finite sequence of elementary row operations, can be transformed and simplified as (13) (the detailed derivation of (13) is given in Appendix B). From the results, we can see that the simplified form of M matrix has six linearly independent rows; therefore, r
a
n
k
(
M
) = 6. Consider(12)
M
=
∇
L
0
h
s
∇
L
f
ν
1
1
h
s
∇
L
f
ω
1
1
h
s
∇
L
f
ν
2
1
h
s
∇
L
f
ω
2
1
h
s
∇
L
f
ν
3
1
h
s
∇
L
f
ω
3
1
h
s
=
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
1
s
i
n
δ
1
ρ
1
2
-
c
o
s
δ
1
ρ
1
-
c
o
s
δ
1
ρ
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-
s
i
n
φ
1
ρ
1
2
c
o
s
φ
1
ρ
1
0
0
0
0
0
0
0
0
0
0
0
0
0
s
i
n
δ
2
ρ
2
2
-
c
o
s
δ
2
ρ
2
-
c
o
s
δ
2
ρ
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-
s
i
n
φ
2
ρ
2
2
c
o
s
φ
2
ρ
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
,
(13)
M
⟹
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
O
14
×
6
.We can see that their gradients of the necessary Lie derivatives ofh
(
s
) are equal to the same order gradients of the time derivatives of h
(
s
) from proposition  2 in [21].Theorem 3.
The rank ofM given by (12) is four if(1)
ν
i
>
0
,
ν
j
>
0, where i
=
1
,
…
,
n, j
=
2
,
…
,
n
+
1;
(2)
φ
i
=
k
1, α
i
=
k
2, where i
=
1
,
…
,
n, and k
1, k
2 are constants; that is, the leader robot and the follower robot at all levels do parallel linear motion.Proof.
Suppose thatφ
i
=
k
1, α
i
=
k
2, and k
1, k
2 are constants; time differentiation of the functions φ
i and α
i is zero, and thus we can get (14)–(16); further we can get simplified form of M matrix as (17),(14)
∇
h
i
0
s
=
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
1
,
(15)
h
i
1
s
=
h
i
2
s
=
⋯
=
h
i
n
-
1
s
=
0
,
(16)
∇
h
i
1
s
=
∇
h
i
2
s
=
⋯
=
∇
h
i
n
-
1
s
=
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
,
(17)
M
⟹
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
1
O
16
×
6
.
From the results of (17), we can see that M matrix has four linearly independent rows; therefore, rank (M) = 4.Clearly, the gradients of the second-order and higher order Lie derivatives are linearly dependent on the row of the observability matrix corresponding to the gradients of first-order and zeroth-order Lie derivatives. Therefore, we have the observability matrix for the cascade formation, using the gradients of the first-order and zeroth-order Lie derivatives, as (12). The rank of the observability matrix does not change with the increasing order of the Lie derivatives; that is, rank (M) = 6 in Theorem 2 and rank (M) = 4 in Theorem 3.When the robots move along the curvilinear trajectory, rank (M) = 6 (see Theorem 2); namely, matrix M has full rank and system S
n is locally weakly observable (as defined in Lemma 1). When the robots move along the rectilinear trajectory, rank (M) = 4 (see Theorem 3); namely, M matrix has no full rank and system S
n is not locally weakly observable.For the formation control problem to be solvable, the system must be observable. If system is locally weakly observable, the system output can convey an information rich enough to allow the observer to provide a correct estimate of the state, thus improving the formation control effectively. For a system not being observable means that the output does not convey an information rich enough to allow the observer to provide a correct estimate of the state, thus affecting the formation control negatively.


### 2.3. Nonlinear Controllability of Cascade Multirobots

The “duality” between controllability and observability (which is, mathematically, just the duality between vector fields and differential forms) is presented in [31]. The design of the input-output feedback controller is a feedback linearizing process. According to [22, 32–34], we know that complete state-controllability and complete observability are equivalent, so the controllability of the robot formation system is simply presented as follows.In order analyze controllability of the cascade robot formation system, the state equation of (5) is restated as (18). According to the concept of nonlinear controllability as defined in [31], system (18) is said to be locally controllable if the vector fields g
i, f and their consecutive Lie brackets f
,
g
i
,
f
,
f
,
g
i
,
…, span the whole state space. Consider(18)
s
˙
=
f
s
+
∑
i
=
1
3
n
g
i
s
u
i
.Suppose ∇
M
(
s
)
=
…
,
g
i
,
…
,
f
,
g
i
,
…
,
f
,
f
,
g
i
,
…. We can easily obtain (19) according to the state equation of (5) and (18); that is, r
a
n
k
(
∇
M
(
s
)
)
=
3
n, and the controllability rank condition is satisfied, so the the cascade robot formation system is said to be controllable. Consider(19)
r
a
n
k
∇
M
s
=
r
a
n
k
…
,
g
i
,
…
,
f
,
g
i
,
…
,
f
,
f
,
g
i
,
…
=
3
n
.


## 2.1. Modeling

Consider the leader-follower setup in Figure1, where R1 is the leader of R2 and R2 is the leader of R3. The control inputs of R1, R2, and R3 are the linear and angular velocities ν
1
ω
1, ν
2
ω
2, and ν
3
ω
3, respectively. ρ
i is the distance from the centroid of the leader robot to the centroid of the follower robot. φ
i is the view-angle from the y-axis of the follower robot to the centroid of the leader robot. θ
i and θ
j are the orientations of the leader robot and the follower robot with respect to the world frame W, respectively, while α
i is the relative orientation between the leader robot and the follower robot; that is, α
i
≜
θ
i
-
θ
j.Figure 1
Leader-follower coordinate representation.With reference to Figure1, the kinematic model of a cascade leader-follower robots’ formation can be expressed as (1)
S
n
:
s
˙
=
f
s
,
u
=
F
s
U
y
=
h
s
=
h
1
T
s
,
h
2
T
s
T
,where state vector s
≜
s
1
T
,
s
2
T
T, s
1
≜
ρ
1
φ
1
α
1
T, s
2
≜
ρ
2
φ
2
α
2
T, input vector U
≜
ν
1
ω
1
ν
2
ω
2
ν
3
ω
3
T, and output vector y
=
h
s
=
h
1
T
s
,
h
2
T
s
T,  h
1
≜
φ
1
α
1
T
,
h
2
≜
φ
2
α
2
T, and where δ
i
=
φ
i
+
α
i
i
=
1,2,(2)
F
s
=
cos
⁡
δ
1
0
-
cos
⁡
φ
1
0
0
0
-
sin
⁡
δ
1
ρ
1
0
sin
⁡
φ
1
ρ
1
-
1
0
0
0
1
0
-
1
0
0
0
0
cos
⁡
δ
2
0
-
cos
⁡
φ
2
0
0
0
-
sin
⁡
δ
2
ρ
2
0
sin
⁡
φ
2
ρ
2
-
1
0
0
0
1
0
-
1
,the detailed derivation of F
(
s
) is given in Appendix A.The kinematic model of multistage cascade leader-follower robots’ formation can be readily retrieved as an extension of (1):(3)
S
n
:
⁡
s
˙
=
f
s
,
u
=
f
1
T
s
1
,
u
1
,
…
,
f
n
+
1
T
s
n
,
u
n
+
1
T
y
=
h
s
=
h
1
T
s
1
,
…
,
h
n
T
s
n
T
,where s
=
s
1
T
,
s
2
T
,
…
,
s
n
T
T
∈
R
3
n is the state vector, U
=
u
1
T
,
u
2
T
,
…
,
u
n
+
1
T
T
∈
R
2
(
n
+
1
) is the control input vector, h
:
R
2
n
↦
R
2
n, and f
:
R
3
(
n
+
1
)
×
Λ
↦
R
3
(
n
+
1
).


## 2.2. Nonlinear Observability of Cascade Multirobots

The functionf can be separated into a summation of independent functions in the special case, and each one is excited by a different component of the control input vector; (1) can be restated as (4)
S
n
:
s
˙
=
f
s
,
u
=
f
ν
1
ν
1
+
f
w
1
ω
1
+
f
ν
2
ν
2
+
f
w
2
ω
2
+
f
ν
3
ν
3
+
f
w
3
ω
3
y
=
h
s
=
h
1
T
s
,
h
2
T
s
T
=
φ
1
α
1
φ
2
α
2
T
.Accordingly, (3) can be also restated as(5)
S
n
:
⁡
s
˙
=
f
s
,
u
=
f
ν
1
ν
1
+
f
w
1
ω
1
+
⋯
+
f
ν
n
+
1
ν
n
+
1
+
f
w
n
+
1
ω
n
+
1
y
=
h
s
=
h
1
T
s
1
,
…
,
h
n
T
s
n
T
.In order to analyse the observability of the entire leader-follower system, the nonlinear observability rank criterion which is developed in [31] is used. The observability matrix of (5) is defined as the matrix with rows,(6)
M
=
∇
L
f
ν
i
,
…
,
f
ν
j
,
f
ω
i
,
…
,
f
ω
j
p
h
k
s
, where i
,
k
=
1
,
…
,
n, j
=
2
,
…
,
n
+
1, p
∈
N, and ∇ represents the gradient operator.Lemma 1.
SystemS
n is locally weakly observable if its observability matrix M which is constructed by the row vector has full rank; for example, in our case r
a
n
k
(
M
)
=
3
n
,
n refers to the number of follower robots.In (4),(7)
f
ν
1
=
cos
⁡
δ
1
-
sin
⁡
δ
1
ρ
1
0
0
0
0
T
,
f
ω
1
=
0
0
1
0
0
0
T
,
f
ν
2
=
-
cos
⁡
φ
1
sin
⁡
φ
1
ρ
1
0
cos
⁡
δ
2
-
sin
⁡
δ
2
ρ
2
0
T
,
f
ω
2
=
0
-
1
-
1
0
0
1
T
,
f
ν
3
=
0
0
0
-
cos
⁡
φ
2
sin
⁡
φ
2
ρ
2
0
T
,
f
ω
3
=
0
0
0
0
-
1
-
1
T
.The necessary Lie derivatives ofh
(
s
) and their gradients are computed; then observability matrix M is obtained.Zeroth-order Lie derivative is(8)
L
0
h
s
=
h
s
=
φ
1
α
1
φ
2
α
2
Tand its gradient is(9)
∇
L
0
h
s
=
∇
L
0
φ
1
α
1
φ
2
α
2
=
∂
φ
1
∂
ρ
1
∂
φ
1
∂
φ
1
∂
φ
1
∂
α
1
∂
φ
1
∂
ρ
2
∂
φ
1
∂
φ
2
∂
φ
1
∂
α
2
∂
α
1
∂
ρ
1
∂
α
1
∂
φ
1
∂
α
1
∂
α
1
∂
α
1
∂
ρ
2
∂
α
1
∂
φ
2
∂
α
1
∂
α
2
∂
φ
2
∂
ρ
1
∂
φ
2
∂
φ
1
∂
φ
2
∂
α
1
∂
φ
2
∂
ρ
2
∂
φ
2
∂
φ
2
∂
φ
2
∂
α
2
∂
α
2
∂
ρ
1
∂
α
2
∂
φ
1
∂
α
2
∂
α
1
∂
α
2
∂
ρ
2
∂
α
2
∂
φ
2
∂
α
2
∂
α
2
=
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
1
.The first-order Lie derivatives of the functionh
(
s
) with respect to f
ν and f
ω are defined as (10)
L
f
ν
1
1
h
1
s
=
∇
L
0
h
1
s
·
f
ν
1
=
-
sin
⁡
δ
1
ρ
1
0
T
,
L
f
ω
1
1
h
1
s
=
∇
L
0
h
1
s
·
f
ω
1
=
0
1
T
,
L
f
ν
2
1
h
1
s
=
∇
L
0
h
1
s
·
f
ν
2
=
sin
⁡
φ
1
ρ
1
0
T
,
L
f
ω
2
1
h
1
s
=
∇
L
0
h
1
s
·
f
ω
2
=
-
1
-
1
T
,
L
f
ν
2
1
h
2
s
=
∇
L
0
h
2
s
·
f
ν
2
=
-
sin
⁡
δ
2
ρ
2
0
T
,
L
f
ω
2
1
h
2
s
=
∇
L
0
h
2
s
·
f
ω
2
=
0
1
T
,
L
f
ν
3
1
h
2
s
=
∇
L
0
h
2
s
·
f
ν
3
=
sin
⁡
φ
2
ρ
2
0
T
,
L
f
ω
3
1
h
2
s
=
∇
L
0
h
2
s
·
f
ω
3
=
-
1
-
1
T
,where denotes the vector inner product, with their gradients given by(11)
∇
L
f
ν
1
1
h
1
s
=
s
i
n
δ
1
ρ
1
2
-
c
o
s
δ
1
ρ
1
-
c
o
s
δ
1
ρ
1
0
0
0
0
0
0
0
0
0
,
∇
L
f
ω
1
1
h
1
s
=
0
0
0
0
0
0
0
0
0
0
0
0
,
∇
L
f
ν
2
1
h
1
s
=
-
s
i
n
φ
1
ρ
1
2
c
o
s
φ
1
ρ
1
0
0
0
0
0
0
0
0
0
0
,
∇
L
f
ω
2
1
h
1
s
=
0
0
0
0
0
0
0
0
0
0
0
0
,
∇
L
f
ν
2
1
h
2
s
=
0
0
0
s
i
n
δ
2
ρ
2
2
-
c
o
s
δ
2
ρ
2
-
c
o
s
δ
2
ρ
2
0
0
0
0
0
0
,
∇
L
f
ω
2
1
h
2
s
=
0
0
0
0
0
0
0
0
0
0
0
0
,
∇
L
f
ν
3
1
h
2
s
=
0
0
0
-
s
i
n
φ
2
ρ
2
2
c
o
s
φ
2
ρ
2
0
0
0
0
0
0
0
,
∇
L
f
ω
3
1
h
2
s
=
0
0
0
0
0
0
0
0
0
0
0
0
.According to the form of (6), row vector matrix M of the cascade leader-follower formation described in Figure 1 can be expressed using gradients of Lie derivatives up to the first order as (12).Theorem 2.
The rank ofM given by (12) is six if(1)
ν
i
>
0
,
ν
j
>
0, where i
=
1
,
…
,
n, j
=
2
,
…
,
n
+
1;
(2)
φ
i
≠
0, i
=
1
,
…
,
n; that is, the follower robot at all levels, which is measuring the bearing, does not move along the line joining the leader robot and the follower robot;
(3)
the leader robot and follower robot at all levels can not do parallel linear motion.Proof.
Meeting the three prerequisite requirements above,M matrix in (12), by means of a finite sequence of elementary row operations, can be transformed and simplified as (13) (the detailed derivation of (13) is given in Appendix B). From the results, we can see that the simplified form of M matrix has six linearly independent rows; therefore, r
a
n
k
(
M
) = 6. Consider(12)
M
=
∇
L
0
h
s
∇
L
f
ν
1
1
h
s
∇
L
f
ω
1
1
h
s
∇
L
f
ν
2
1
h
s
∇
L
f
ω
2
1
h
s
∇
L
f
ν
3
1
h
s
∇
L
f
ω
3
1
h
s
=
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
1
s
i
n
δ
1
ρ
1
2
-
c
o
s
δ
1
ρ
1
-
c
o
s
δ
1
ρ
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-
s
i
n
φ
1
ρ
1
2
c
o
s
φ
1
ρ
1
0
0
0
0
0
0
0
0
0
0
0
0
0
s
i
n
δ
2
ρ
2
2
-
c
o
s
δ
2
ρ
2
-
c
o
s
δ
2
ρ
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-
s
i
n
φ
2
ρ
2
2
c
o
s
φ
2
ρ
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
,
(13)
M
⟹
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
O
14
×
6
.We can see that their gradients of the necessary Lie derivatives ofh
(
s
) are equal to the same order gradients of the time derivatives of h
(
s
) from proposition  2 in [21].Theorem 3.
The rank ofM given by (12) is four if(1)
ν
i
>
0
,
ν
j
>
0, where i
=
1
,
…
,
n, j
=
2
,
…
,
n
+
1;
(2)
φ
i
=
k
1, α
i
=
k
2, where i
=
1
,
…
,
n, and k
1, k
2 are constants; that is, the leader robot and the follower robot at all levels do parallel linear motion.Proof.
Suppose thatφ
i
=
k
1, α
i
=
k
2, and k
1, k
2 are constants; time differentiation of the functions φ
i and α
i is zero, and thus we can get (14)–(16); further we can get simplified form of M matrix as (17),(14)
∇
h
i
0
s
=
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
1
,
(15)
h
i
1
s
=
h
i
2
s
=
⋯
=
h
i
n
-
1
s
=
0
,
(16)
∇
h
i
1
s
=
∇
h
i
2
s
=
⋯
=
∇
h
i
n
-
1
s
=
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
,
(17)
M
⟹
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
1
O
16
×
6
.
From the results of (17), we can see that M matrix has four linearly independent rows; therefore, rank (M) = 4.Clearly, the gradients of the second-order and higher order Lie derivatives are linearly dependent on the row of the observability matrix corresponding to the gradients of first-order and zeroth-order Lie derivatives. Therefore, we have the observability matrix for the cascade formation, using the gradients of the first-order and zeroth-order Lie derivatives, as (12). The rank of the observability matrix does not change with the increasing order of the Lie derivatives; that is, rank (M) = 6 in Theorem 2 and rank (M) = 4 in Theorem 3.When the robots move along the curvilinear trajectory, rank (M) = 6 (see Theorem 2); namely, matrix M has full rank and system S
n is locally weakly observable (as defined in Lemma 1). When the robots move along the rectilinear trajectory, rank (M) = 4 (see Theorem 3); namely, M matrix has no full rank and system S
n is not locally weakly observable.For the formation control problem to be solvable, the system must be observable. If system is locally weakly observable, the system output can convey an information rich enough to allow the observer to provide a correct estimate of the state, thus improving the formation control effectively. For a system not being observable means that the output does not convey an information rich enough to allow the observer to provide a correct estimate of the state, thus affecting the formation control negatively.


## 2.3. Nonlinear Controllability of Cascade Multirobots

The “duality” between controllability and observability (which is, mathematically, just the duality between vector fields and differential forms) is presented in [31]. The design of the input-output feedback controller is a feedback linearizing process. According to [22, 32–34], we know that complete state-controllability and complete observability are equivalent, so the controllability of the robot formation system is simply presented as follows.In order analyze controllability of the cascade robot formation system, the state equation of (5) is restated as (18). According to the concept of nonlinear controllability as defined in [31], system (18) is said to be locally controllable if the vector fields g
i, f and their consecutive Lie brackets f
,
g
i
,
f
,
f
,
g
i
,
…, span the whole state space. Consider(18)
s
˙
=
f
s
+
∑
i
=
1
3
n
g
i
s
u
i
.Suppose ∇
M
(
s
)
=
…
,
g
i
,
…
,
f
,
g
i
,
…
,
f
,
f
,
g
i
,
…. We can easily obtain (19) according to the state equation of (5) and (18); that is, r
a
n
k
(
∇
M
(
s
)
)
=
3
n, and the controllability rank condition is satisfied, so the the cascade robot formation system is said to be controllable. Consider(19)
r
a
n
k
∇
M
s
=
r
a
n
k
…
,
g
i
,
…
,
f
,
g
i
,
…
,
f
,
f
,
g
i
,
…
=
3
n
.


## 3. Input-Output Feedback Control

In order to achieve the feedback control, state estimation is necessary. We choose the UKF algorithm to estimate the state of the robots in this paper. Combining the UKF algorithm with the input-output state feedback control, the formation control can make full use of bearing-only observations of the cascade robot formation to realize motion control, so the desired formation control can be achieved.


### 3.1. UKF Algorithm

In order to control the formation, an estimates
^ of the true state s is required. The UKF is designed to estimate angle information, that is, φ
i
α
i
T of the state s. The UKF uses the input vector U and the output vector y
=
h
(
s
)
=
h
1
T
s
,
…
,
h
n
T
(
s
)
T. Combining the Euler forward method with the UKF algorithm, we can obtain the more accurate information. We assume additive noise in both the process equation (20) and the measurement equation (21):(20)
s
˙
=
F
s
U
+
O
,
(21)
y
=
G
s
+
N
,where G is the output transition matrix and O and N are white Gaussian noises with zero mean and covariance matrices P
L and P
N, respectively. We assume that s
(
0
), O, and N are uncorrelated. We apply the Euler forward method with sampling time T
c to discretize (20) and we can obtain (22)
s
k
+
1
=
Γ
s
k
,
u
k
+
T
c
O
,where Γ
(
s
(
k
)
,
u
(
k
)
)
=
T
c
F
(
s
)
U
+
s
(
k
) and k
∈
N.In [35, 36], the UKF is based on the Unscented Transformation, which includes a prediction and correction steps. In order to improve the robot localization, the process equation (20) and the measurement equation (21) have to be sampled, which are detailed in Appendix C.


### 3.2. Input-Output State Feedback Control

The control law is designed for R1and R2, which is also similarly designed to other robots of cascade formation. Changing the form of the state equation in (1), the first level kinematic system model of cascade formation is equivalent to (23) and (24) is obtained by taking the derivative of α
1
=
θ
1
-
θ
2, (23)
s
˙
r
=
A
s
U
1
+
B
s
U
2
,
(24)
α
˙
1
=
ω
1
-
ω
2
,where s
r
≜
ρ
1
φ
1
T is the reduced state-space vector; more details can be found in Appendix D. Matrix E
s
2
×
4 is the upper-left submatrix of F
(
s
), and matrices A
s
2
×
2 and B
(
s
)
2
×
2 are the two upper-left and right submatrices of E
s
2
×
4, respectively.Using standard techniques of I/O linearization in [21, 37, 38], we propose an input-output state feedback control for the robot formation control. According to (23), let us consider the following control input of R2:(25)
U
2
≜
ν
ω
2
T
=
B
s
-
1
C
-
A
s
U
1
,where A
(
s
)
=
cos
⁡
δ
0
-
sin
⁡
δ
1
/
ρ
1
0, B
(
s
)
=
-
cos
⁡
φ
1
0
sin
⁡
φ
1
/
ρ
1
-
1, and C
=
-
K
(
s
r
-
s
r
i
d
e
)
,
K
=
d
i
a
g
k
1
k
2, with k
1, k
2
>
0. The superscript “ide” refers to the desired state, and C is auxiliary control parameter. Equation (25) acts in (23) as a feedback linearizing control, combining with (24), so that the closed-loop dynamics become(26)
s
˙
r
=
C
=
-
K
s
r
-
s
r
i
d
e
=
k
1
ρ
1
i
d
e
-
ρ
1
k
2
φ
1
i
d
e
-
φ
1
,
α
˙
1
=
ω
1
-
ω
2
.


## 3.1. UKF Algorithm

In order to control the formation, an estimates
^ of the true state s is required. The UKF is designed to estimate angle information, that is, φ
i
α
i
T of the state s. The UKF uses the input vector U and the output vector y
=
h
(
s
)
=
h
1
T
s
,
…
,
h
n
T
(
s
)
T. Combining the Euler forward method with the UKF algorithm, we can obtain the more accurate information. We assume additive noise in both the process equation (20) and the measurement equation (21):(20)
s
˙
=
F
s
U
+
O
,
(21)
y
=
G
s
+
N
,where G is the output transition matrix and O and N are white Gaussian noises with zero mean and covariance matrices P
L and P
N, respectively. We assume that s
(
0
), O, and N are uncorrelated. We apply the Euler forward method with sampling time T
c to discretize (20) and we can obtain (22)
s
k
+
1
=
Γ
s
k
,
u
k
+
T
c
O
,where Γ
(
s
(
k
)
,
u
(
k
)
)
=
T
c
F
(
s
)
U
+
s
(
k
) and k
∈
N.In [35, 36], the UKF is based on the Unscented Transformation, which includes a prediction and correction steps. In order to improve the robot localization, the process equation (20) and the measurement equation (21) have to be sampled, which are detailed in Appendix C.


## 3.2. Input-Output State Feedback Control

The control law is designed for R1and R2, which is also similarly designed to other robots of cascade formation. Changing the form of the state equation in (1), the first level kinematic system model of cascade formation is equivalent to (23) and (24) is obtained by taking the derivative of α
1
=
θ
1
-
θ
2, (23)
s
˙
r
=
A
s
U
1
+
B
s
U
2
,
(24)
α
˙
1
=
ω
1
-
ω
2
,where s
r
≜
ρ
1
φ
1
T is the reduced state-space vector; more details can be found in Appendix D. Matrix E
s
2
×
4 is the upper-left submatrix of F
(
s
), and matrices A
s
2
×
2 and B
(
s
)
2
×
2 are the two upper-left and right submatrices of E
s
2
×
4, respectively.Using standard techniques of I/O linearization in [21, 37, 38], we propose an input-output state feedback control for the robot formation control. According to (23), let us consider the following control input of R2:(25)
U
2
≜
ν
ω
2
T
=
B
s
-
1
C
-
A
s
U
1
,where A
(
s
)
=
cos
⁡
δ
0
-
sin
⁡
δ
1
/
ρ
1
0, B
(
s
)
=
-
cos
⁡
φ
1
0
sin
⁡
φ
1
/
ρ
1
-
1, and C
=
-
K
(
s
r
-
s
r
i
d
e
)
,
K
=
d
i
a
g
k
1
k
2, with k
1, k
2
>
0. The superscript “ide” refers to the desired state, and C is auxiliary control parameter. Equation (25) acts in (23) as a feedback linearizing control, combining with (24), so that the closed-loop dynamics become(26)
s
˙
r
=
C
=
-
K
s
r
-
s
r
i
d
e
=
k
1
ρ
1
i
d
e
-
ρ
1
k
2
φ
1
i
d
e
-
φ
1
,
α
˙
1
=
ω
1
-
ω
2
.


## 4. Simulation Validation

The effectiveness of the cascade robot formation control approach proposed in this paper is verified by three Pioneer-3at robots with laser radar and panoramic camera, and the angle information is easily obtained by these sensors. In order to validate the proposed formation control approach in a real scenario, we set up 3D simulation platform using the powerful Webots 7, by which the real environment can be simulated and simulation datum can also be imported into Matlab and analysed there. The cascade formation of three robots involves one leader following controller at all levels (R2 following R1 and R3following R2). Simulation scenarios are shown as in Figure 2.Figure 2
The simulation scene of mobile robots.The cascade robot formation undergoes the piecewise rectilinear-curvilinear combination trajectory 1, trajectory 2 looking like 8, and spiral trajectory 3 that are particularly suited for checking the observability conditions discussed in Section2.2. Simulation results are shown as in Figures 3, 4, and 5, respectively.Figure 3
Simulation results of trajectory 1.
(a)
Trajectory 1 of the robots
(b)
Time history of observation angle error
(c)
Time history of direction angle error
(d)
Time history of velocities of R2 and R3Figure 4
Simulation results of trajectory 2.
(a)
Trajectory 2 of the robots
(b)
Time history of observation angle error
(c)
Time history of direction angle error
(d)
Time history of velocities of R2and R3Figure 5
Simulation results of trajectory 3.
(a)
Trajectory 3 of the robots
(b)
Time history of observation angle error
(c)
Time history of direction angle error
(d)
Time history of distance error
(e)
Time history of velocities of R2 and R3In the simulation, the follower robot at all levels of cascade formation can adjust its movement according to movement of the leader robot. UKF is employed to estimate the state of robots. The input-output state feedback control is used to stabilize the cascade robot formation.


### 4.1. Initial Conditions

The following velocity inputs have been assigned to the leader robot when the cascade robot formation follows trajectory 1. Consider(27)
v
1
t
=
4.5
m
/
s
,
ω
1
t
=
0
r
a
d
/
s
,
t
∈
0,2
,
6,8
-
π
4
r
a
d
/
s
,
t
∈
2,6
π
4
r
a
d
/
s
,
t
∈
8,14
.The following velocity inputs have been assigned to the leader robot when the cascade robot formation follows trajectory 2. Consider(28)
v
1
t
=
4.5
m
/
s
,
ω
1
t
=
π
3
r
a
d
/
s
,
t
∈
0,5.5
-
π
3
r
a
d
/
s
,
t
∈
5.5,11
.The linear and angular velocities of leader robot are the functions of timet when the cascade robot formation follows trajectory 3.The initial configuration vectors of leader robot and follower robots are(29)
x
1
0
y
1
0
θ
1
0
T
=
4.5
2.0
π
2
T
,
x
2
0
y
2
0
θ
2
0
T
=
4.2
1.7
π
2
T
,
x
3
0
y
3
0
θ
3
0
T
=
3.9
1.4
π
2
T
.Parameters of UKF and the input-output state feedback law are as follows.The formation considered here consists of three robots. We sets
(
0
)
=
0.4474
2.4223
0
0.4360
2.4469
0
T and s
r
i
d
e
=
0.45
5
π
/
4
0.45
5
π
/
4, where distances are in meters and angles in radians. The gains of the controller are k
1
=
0.15
,
k
2
=
0.25. The other parameters of the UKF are T
c
=
0.01
s, P
L
=
d
i
a
g
h
h, P
N
=
d
i
a
g
h
h, and P
=
1.13
1.13, where h
=
3.0
×
1
0
-
2
r
a
d
2. White Gaussian noise has been injected into the observations.In the initial moment, the robots need to adjust the original postures to the ideal postures, and the original position error is relatively large, which has no effect on the robot formation after the robots enter the normal movement, so we do not consider position error at the initial moment when all kinds of simulations are analysed.


### 4.2. Simulation Analysis

Figure3(a) shows the trajectories of the leader and the follower robots. It shows that R1, R2, and R3 can maintain the desired formation when three robots undergo trajectory 1. Figure 3(b) shows the observation angle error is almost zero when the robot moves along trajectory 1. The maximum error is −0.0003 even when the trajectory changes suddenly. Figure 3(c) shows the direction angle error is also very small when the robot moves along trajectory 1, and the maximum error is −0.008 when the trajectory changes suddenly. Figure 3(d) shows that the velocities of the follower robots change slightly when the follower robots move along rectilinear trajectory (or curvilinear trajectory). When the trajectories of robot formation change, the velocities of robot formation change significantly. The velocities of the follower robots are almost the same as the initial velocity of the leader robot when the robot formation moves along rectilinear trajectory.Figure4(a) shows the trajectories of the three robots following trajectory 2. According to observable theory described in Section 2.2, the cascade leader-follower formation system is observable when three robots undergo trajectory 2, so it is evident that R1, R2, and R3 maintain the desired formation. Figure 4(b) shows that observation angle error is almost zero when the robot formation moves along trajectory 2. The maximum error is −0.0005 even when the trajectory changes suddenly. Figure 4(c) shows that direction angle error is also almost zero when the robot formation moves along trajectory 2, and the maximum error is −0.02107 even when the direction changes suddenly. Figure 4(d) shows that the velocities of the follower robots have a big change when direction of the robot formation motion has great change. The velocity of the follower robot outside of the leader robot is greater than the initial velocity of the leader robot, and the velocity of the follower robot inside the leader robot is less than the initial velocity of the leader robot.Figure5(a) shows the trajectories of the three robots following trajectory 3. According to observable theory described in Section 2.2, the cascade leader-follower formation system is observable when three robots undergo trajectory 3. We can see that R1, R2, and R3 maintain the desired formation. Figure 5(b) shows observation angle error is almost zero when the robot formation moves along trajectory 3. Figure 5(c) shows that direction angle error almost converged to zero after t
=
2
s when the robot formation moves along trajectory 3. Figure 5(d) shows that distance error is also almost zero when the robot formation moves along trajectory 3. Figure 5(e) shows that the velocities of the follower robots have a trend of increase when the robot formation moves along trajectory 3.


## 4.1. Initial Conditions

The following velocity inputs have been assigned to the leader robot when the cascade robot formation follows trajectory 1. Consider(27)
v
1
t
=
4.5
m
/
s
,
ω
1
t
=
0
r
a
d
/
s
,
t
∈
0,2
,
6,8
-
π
4
r
a
d
/
s
,
t
∈
2,6
π
4
r
a
d
/
s
,
t
∈
8,14
.The following velocity inputs have been assigned to the leader robot when the cascade robot formation follows trajectory 2. Consider(28)
v
1
t
=
4.5
m
/
s
,
ω
1
t
=
π
3
r
a
d
/
s
,
t
∈
0,5.5
-
π
3
r
a
d
/
s
,
t
∈
5.5,11
.The linear and angular velocities of leader robot are the functions of timet when the cascade robot formation follows trajectory 3.The initial configuration vectors of leader robot and follower robots are(29)
x
1
0
y
1
0
θ
1
0
T
=
4.5
2.0
π
2
T
,
x
2
0
y
2
0
θ
2
0
T
=
4.2
1.7
π
2
T
,
x
3
0
y
3
0
θ
3
0
T
=
3.9
1.4
π
2
T
.Parameters of UKF and the input-output state feedback law are as follows.The formation considered here consists of three robots. We sets
(
0
)
=
0.4474
2.4223
0
0.4360
2.4469
0
T and s
r
i
d
e
=
0.45
5
π
/
4
0.45
5
π
/
4, where distances are in meters and angles in radians. The gains of the controller are k
1
=
0.15
,
k
2
=
0.25. The other parameters of the UKF are T
c
=
0.01
s, P
L
=
d
i
a
g
h
h, P
N
=
d
i
a
g
h
h, and P
=
1.13
1.13, where h
=
3.0
×
1
0
-
2
r
a
d
2. White Gaussian noise has been injected into the observations.In the initial moment, the robots need to adjust the original postures to the ideal postures, and the original position error is relatively large, which has no effect on the robot formation after the robots enter the normal movement, so we do not consider position error at the initial moment when all kinds of simulations are analysed.


## 4.2. Simulation Analysis

Figure3(a) shows the trajectories of the leader and the follower robots. It shows that R1, R2, and R3 can maintain the desired formation when three robots undergo trajectory 1. Figure 3(b) shows the observation angle error is almost zero when the robot moves along trajectory 1. The maximum error is −0.0003 even when the trajectory changes suddenly. Figure 3(c) shows the direction angle error is also very small when the robot moves along trajectory 1, and the maximum error is −0.008 when the trajectory changes suddenly. Figure 3(d) shows that the velocities of the follower robots change slightly when the follower robots move along rectilinear trajectory (or curvilinear trajectory). When the trajectories of robot formation change, the velocities of robot formation change significantly. The velocities of the follower robots are almost the same as the initial velocity of the leader robot when the robot formation moves along rectilinear trajectory.Figure4(a) shows the trajectories of the three robots following trajectory 2. According to observable theory described in Section 2.2, the cascade leader-follower formation system is observable when three robots undergo trajectory 2, so it is evident that R1, R2, and R3 maintain the desired formation. Figure 4(b) shows that observation angle error is almost zero when the robot formation moves along trajectory 2. The maximum error is −0.0005 even when the trajectory changes suddenly. Figure 4(c) shows that direction angle error is also almost zero when the robot formation moves along trajectory 2, and the maximum error is −0.02107 even when the direction changes suddenly. Figure 4(d) shows that the velocities of the follower robots have a big change when direction of the robot formation motion has great change. The velocity of the follower robot outside of the leader robot is greater than the initial velocity of the leader robot, and the velocity of the follower robot inside the leader robot is less than the initial velocity of the leader robot.Figure5(a) shows the trajectories of the three robots following trajectory 3. According to observable theory described in Section 2.2, the cascade leader-follower formation system is observable when three robots undergo trajectory 3. We can see that R1, R2, and R3 maintain the desired formation. Figure 5(b) shows observation angle error is almost zero when the robot formation moves along trajectory 3. Figure 5(c) shows that direction angle error almost converged to zero after t
=
2
s when the robot formation moves along trajectory 3. Figure 5(d) shows that distance error is also almost zero when the robot formation moves along trajectory 3. Figure 5(e) shows that the velocities of the follower robots have a trend of increase when the robot formation moves along trajectory 3.


## 5. Conclusion

A new formation control method is proposed for multirobots which are queued in a single-direction cascade structure. The observability of the cascade leader-follower formation is studied, which shows that the bearing-only observation meets the observability conditions required for the nonlinear system. Based on the cascade bearing-only observations, the UKF is employed for the state estimation of the leader and the follower robots at all levels, which enables the movement control of the follower robots via the input-output feedback control. Simulation results show that the multirobots can rapidly form a cascade formation, run along complex trajectory with very small error, and maintain the desired formation. Further research will include formation control of multirobots in the environment of obstacles.



---

*Source: 5320954-2016-07-17.xml*