Leader-follower Formation Control and Obstacle Avoidance of Multi-robot Based
on Artificial Potential Field 
ZHANG Ying, LI Xu*  
College of Information Engineering, Shanghai Maritime University, Shanghai 201306 
E-mail: yingzhang@shmtu.edu.cn 
 
Abstract: Leader-Follower formation control of multi-robot was studied in this paper. A formation and obstacle avoidance method with
multi-robot based on the combination of Closed-loop control and Artificial Potential Field was presented. According to the position
information of leader, Closed-loop control was introduced to realize the tracking of the follower to the leader, and the formation control
was achieved. The obstacle avoidance could be achieved by Artificial Potential Field method, and the robots can pass the area of obstacle
smoothly. The simulation result shows that the proposed method can achieve the expected control effect, and it can solve this kind of
problems effectively. 
Key Words: Multi-robot formation, Leader-follower, Artificial Potential Field, Obstacle avoidance 
 
INTRODUCTION 
  1 
In recent years, the control and coordination of multi-robot 
[1] has become an interesting research topic in the field of 
robot. Mobile multi-robot technology has been widely used 
in  the  fields  of  industry,  military,  agriculture,  space  and 
marine development [2]. The robot formation control is an 
important  part  of  multi-robot  coordination.  Faced  with 
complicated tasks and fickle condition, obviously, it’s not 
enough to just depend on the ability of single robot. So it 
attracts people’s attention to develop the coordination and 
cooperation of multi-robot to finish the task that can never 
be done by single robot.  The mobile multi-robot formation 
requires  them  as  one  formation  which  can  arrive  at  the 
target area at the same time, and avoid obstacles in a safe 
manner. This kind of group behavior control is the basis to 
solve the problem of mobile multi-robot coordination, and 
it is of great significance to realize the cooperative missions 
of  multi-robot  in  the  distributed  environment  space. 
Nowadays, the formation control has been applied in many 
fields.  For  example,  in  industry  field,  people  control  the 
mobile  multi-robot  to  carry  large  objects  by  certain 
formation  [3].  In  military,  multiple  autonomous  vehicles 
have been used to patrol or reconnoiter [4]. In the field of 
police, people control the multi-robot to form a cambered 
encirclement or to arrest the invaders and so on [5]. AGV 
(Automated  Guided  Vehicles)  have  been  widely  used  in 
practice.  In  this  paper,  we  do  research  on  multi-robot 
formation on the background of AGV. 
 
  Normally  there  are  three  kinds  of  approaches  for 
formation  control  with  applications  towards  multi-robots: 
leader-follower method, behavior-based method and virtual 
structure  method  [6].  In  the  behavioral  approach,  the 
control  action  for  each  robot  is  derived  by  a  weighted 
                                                           
This  work  is  supported  by  National  Nature  Science  Foundation  of 
China under Grant No. 61273068, Nature Science Foundation of Shanghai 
under  Grant  No.  12ZR1412600,  and  Scientific  Research  Innovation 
Project of Shanghai Education Committee under Grant No. 13YZ084. 
* Correspondence Author: LI Xu 
average of each desired behavior [4]. This method has clear 
formation  feedback,  and  it  also  realizes  the  distributed 
control, but it does not clearly define the group behavior. 
The mathematical analysis cannot be done and the stability 
of  the  formation  cannot  be  guaranteed.  In  the  virtual 
structure approach [7], the entire formation is treated as a 
single rigid body, and the motion of each agent is derived 
from  the  trajectory  of  a  corresponding  point  on  the 
structure. This method can control the motion of the whole 
robot  formation  by  defining  the  rigid  body,  but  it  cannot 
change  the  formation  when  the  environment  has  been 
changed,  which  limits  the  scope  of  its  application.  The 
leader-follower method is also named as master-slave mode 
[8]. In the group formed by multi-mobile robots, a certain 
robot is designated as the leader and the other robots are the 
followers.  It  realizes  the  formation  control  through  the 
followers  keeping  a  certain  angle  and  distance  with  the 
leader,  and  it  can  establish  different  topologies.  The 
leader-follower  method  also  has  a  clear  definition  to  the 
whole formation, so we can make mathematical analysis to 
guarantee the stability of the formation. In this paper, we 
put  forward  a  method  of  robot  formation  and  obstacle 
avoidance  based  on  Closed-loop  control  and  Artificial 
Potential  Field.  According 
location 
information  [9],  we  can  realize  the  followers’  trajectory 
tracking to the leader and formation control by the l ϕ−  
Closed-loop  control,  and  avoid  the  obstacles  effectively 
combined with the Artificial Potential Field method. 
leader’s 
In the second part of this paper, we put forward a robot 
formation  control  and  obstacle  avoidance  method,  a 
tracking control method is in the third part. We change the 
problem of followers’ tracking control into a problem of the 
control in a system with a certain error. In the fourth part, 
we  make  the  simulation  verification  and  finally  make  a 
prospect for the future’s work. 
to 
the 
978-1-4799-7016-2/15/$31.00 c2015 IEEE
4375
 
2  FORMATTING CONTROL AND 
OBSTACLE  AVOIDANCE 
2.1 The Robot Formation Under the l ϕ−  Control 
The  optimized  motion  method  that  used  in  the  robot 
formation control needs a great amount of calculation, but 
much less using the feedback rules, which can also combine 
with the simple advanced motion planning devices. 
l− . 
Two  kinds  of  motion  models  with  a  feedback  control 
are put forward in paper [10]. In the first scenario, through 
controlling the relative distance and orientation between the 
follower and the leader to command the formation l ϕ− 
and in another scenario, a robot maintains its position in the 
formation  by  maintaining  a  specified  distance  from  two 
robotsl
As shown in Figure 1, in the given system composed by 
two robots, the purpose of the l ϕ−  control is to keep the 
relative  distance  l  and  the  relative  angle ϕbetween  the 
following  robot  and  reference  robot.  As  long  as  the  two 
l ϕ ,  we  can  maintain  a  certain 
values  are  fixed  in (
,
d
formation,  and  each  robot’s  position  information  is 
x y θ =
(
,
i
1,
,
)
2
. 
)
i
,
d
i
i
θ
2
2 , 2  
12l
ϕ
12
θ
1
1, 1  
Fig.1 The structure of leader-follower 
The kinematic equations for the system of two mobile 
robots shown in Figure 1 are given by formula (1): 
 
                                    (1) 
cos
sin
θ
i
θ
i
•
 =
x
v
i
i
•
 =
v
y
i
i
•
θ ω
 =
i
i
 
d> , where d is the distance between the castor wheel 
12l
and the center of rear wheels. 
The law of closed-loop control is described as: 
 
•
=
l
)
12
•
ϕ α ϕ ϕ
12
12
−
lα
12
1
=
−
l
(
d
12
(
d
12
2
3 
)
                           
1α 2αin the equation are positive numbers. So we can get 
the  input  control  variables  of  follower  from  the  equation 
(1), (2), (3): 
γω
cos
=
2
d
−
 =
v
d
2
d
12
ρ ω γ
12
1
ϕ ω ρ γ
sin }
12
1
α ϕ ϕ
{
12
l
2 12
tan
4 
sin
l
12
v
1
−
−
+
+
12
(
)
1
1
2
Where, 
ρ
12
=
α
dl
(
1
12
−
cos
ϕ
12
+
v
l
)
1
12
γ
cos
1
 
So we can get the conclusion that as long as the leader’s 
angular velocity, speed, position and sailing directions are 
designated,  the  follower  can  guarantee  to  sail  toward  a 
position  with  a  relative  distance  l  and  angle  ϕto  the 
leader, so as to keep the formation. 
2.2 Obstacle Avoidance
The  robot  may  meet  some  obstacles  as  sailing,  so  we 
mainly consider the formation control in this case. For the 
leader robot, it will take action to avoid obstacles according 
to the position and direction when it meets the obstacles. 
The follower robot  will change the relative position  with 
the leader as the leader’s position and trajectory changing, 
and  finally  pass  the  area  with  obstacles  by  changing  the 
formation. In this part, we use the Artificial Potential Field 
method  to  realize  the  obstacles  avoidance  of  multiple 
robots. 
We  use  the  model  of  Artificial  Potential  Field  in  the 
paper [11]. In this method, robot is considered as a particle 
in  the  potential  field.  The  artificial  potential  field  is 
composed of the  gravity  function and the repulsive force 
function.  We  build  a  gravitational  potential  field  in  the 
target location, and a repulsive force potential field around 
the obstacle. The two potential energy field work together 
and  form  a  composite  artificial  potential  field.  The  robot 
moves  to  the  destination  under  the  influence  of  the  two 
fields. In math, this motion process can be described as a 
path without collision that a particle seeks a path to avoid 
the obstacles in direction of the search potential  function 
falling from the initial position to the target position. The 
path planned by the potential field method is smooth and 
safe. 
The kinematics equation of the follower is: 
•
l
12
•
ϕ
 
12
=
v
2
=
−
cos
1 { sin
l
12
v
1
γ ω γ
1
1
d
ϕ
12
sin
2
−
v
2
sin
−
d
+
ϕ ω ϕ
12
v
cos
1
−
γ ω γ ω ϕ ω
}
1
1
 sin
1
+
 cos
cos
l
12
−
d
d
12
2
1
1
(2) 
+
=
γ θ ϕ θ
1
2
 are  the  linear 
Where, 
and angular velocities at the center of the axle of each robot. 
In  order  to  avoid  collision  between  robots,  we  specify 
iω =
i
iv  (
)1,2
−
, 
12
1
4376
2015 27th Chinese Control and Decision Conference (CCDC)
 
 
gF q
( )
totalF
q
( )
F q
(
rep
obs
)
 
Fig2. Model of the artificial potential field 
)
As shown in figure 2, let q be the position of the robot, 
q gρ
( ,
 be the distance between the robot and the target g, 
gF  at 
the  gravitational  potential  field 
robot q are defined as: 
gU  and  gravitation 
1
2
)
obsj
q g
( ,
           
          
Let 
                                    (5) 
ξρ=
2
g qξρ=
( , )
)
= …    be  the  position  of  the 
m
1,
,
)
                                           (6) 
 
thj
 be  the  distance  between  robot  q  and 
repU q  and 
gU q
( )
gF q
( )
(
q
j
obsj
q qρ
obstacle, 
( ,
obsjq
,  the  repulsive  potential  field 
the 
gravitation 
F q
(
rep
1
ξ
(
= 
2
if
0                                   
1
ξ
(
= 
q q
( ,
)  > 
So the resultant force of robot q in the APF is: 
 at the robot q are defined as: 
if
0                                                  
 
ρ
s
       (7) 
obs
1
ρ
q q
( ,
1
ρ ρ
2
s
) 
≤
ρ
 
s
    
1
ρ
q q
( ,
U q
( )
ρ
q q
  ( ,
 
(8)
ρ
q q
( ,
ρ
q q
  ( ,
ρ
q q
  ( ,
1
ρ
s
F q
(
rep
)     
2
ρ
 
s
) > 
( )
ρ
s
if
  
if
obsj
obsj
obsj
−
≤
rep
obsj
obsj
obsj
obsj
−
)
obs
)
)
)
)
)
)
2
F
total
q
( )
=
F q
( )
g
m
+
j
−
1
 
F q
(
rep
obsj
)
                     (9) 
3  TRACKING CONTROL 
The robots are considered as a point mass. The leader’s path 
is  defined  by  the  artificial  potential  functions,  depending 
upon location of static obstacles and the goal position. The 
follower  robot  will  follow  the  leader  by  keeping  the 
ϕ . 
separation distance  ijl  and bearing angle ij
In the system given in Figure 1, each robot meets the 
nonholonomic constraints with pure rolling without sliding. 
The motion model is given as: 
•
q
j
=
j
•
x
•
y
j
•
θ
j
=
θ
cos
j
θ
sin
j
0
θ
j
θ
j
−
d
 sin
d
 cos
1
v
j
ω
 
j
                 (10) 
Where d is the distance from the rear axle to the front of the 
robot. 
According  to  the  tracking  control  system  which  was 
presented in the literature [12], the tracking control points 
could be given as formula (11): 
θ
jr
θ
jr
                                  (11) 
                     
cos
sin
=
v
jr
=
•
x
jr
•
y
jr
jr
•
θ ω
jr
=
v
jr
•
x
jr
•
= 
y
jr
•
θ
jr
         So:  
•
q
jr
                                        (12) 
j
j
j
j
j
2
1
jr
jr
jr
,
,
0
0
1
−
−
−
θ
j
θ
j
cos
sin
0
= −
θ
sin
j
θ
cos
j
0
So the errors of tracking controller can be expressed as: 
e
e
e
j
3
     
x y θ  are the actual position and orientation 
)
Where  (
j
y θ  are  the  position  and 
of  the  robot,  and  (
jr
orientation of a virtual trace point of the robot j [11]. 
x
x
y
y
j
θ θ
j
          (13) 
The basic tracking control problems can be transformed 
to formation control. The navigation robot will be placed to 
the position of trace point.  
Suppose j is the followers, and i is the navigation leader. 
The kinematical equation of the robots can be expressed by 
formula (14).  
x
)
,
,
jr
jr
   
=
                      (14) 
θ
cos
i
θ
sin
i
0
−
d
 sin
d
 cos
1
θ
i
v
θ
i
ω
 
i
i
•
x
i
•
y
 
i
•
θ
i
x
x
i
y
y
i
θ θ
i
The  actual  position  and  orientation  of  the  follower  j 
ϕ θ
+ )
d
i
ij
ϕ θ
+ )
d
i
ij
 cos(
 sin(
                (15) 
cos
sin
θ
j
θ
j
=
=
=
l
d
ij
l
d
ij
−
−
d
d
+
+
jr
jr
jr
with respect to leader i can be defined as: 
j
+
+
d
d
−
−
=
=
=
θ
j
θ
j
cos
sin
ϕ θ
 cos( + )
i
ij
ϕ θ
 sin( + )
i
ij
x
x
i
y
y
i
θ θ
j
Using (14), (16) and the simple trigonometric identities, 
                (16) 
l
ij
l
ij
j
j
j
1
the error system can be rewritten as: 
+
e
+
e
e
θ
sin
j
θ
cos
j
0
= −
cos
sin
0
ϕ θ
l
)
d
ij
i
ij
ϕ θ
l
)
d
i
ij
ij
θ θ
i
j
cos(
sin(
θ
j
θ
j
−
−
−
l
d
ij
l
d
ij
0
0
1
2
3
j
j
cos(
sin(
ϕ θ
)
i
ij
ϕ θ
)
ij
i
+
+
(17)
 
2015 27th Chinese Control and Decision Conference (CCDC)
4377
 
 
 
4  SIMULATION OF FORMATION 
CONTROL 
4.1  Simulation of Trajectory Tracking 
We  will  get  the  position  error  of  the  robot  in  different 
trajectory tracking of the route: curves and circular. As the 
illustrating of the simulation results below, we can find that 
robots  can  well  realize  trajectory  tracking  for  curves  and 
circular. 
30
25
20
y
15
10
5
0
0
r
o
r
r
e
 
x
0.5
0
-0.5
0
2
0
-2
0
2
0
-2
0
r
o
r
r
e
 
y
r
o
r
r
e
 
n
o
i
t
t
a
n
e
i
r
o
5
10
15
x
20
25
30
Fig3. Trajectory tracking for curves 
 
5
5
5
10
10
10
15
time
15
time
15
time
20
25
30
20
25
30
20
25
30
Fig4. Position error of trajectory tracking for curves 
 
After further simplification, the formula (17) becomes 
j
1
as formula (18): 
l
d
ij
l
d
ij
e
e
e
=
=
e
2
3
j
j
j
cos(
sin(
ϕ
d
ij
ϕ
d
ij
+
+
e
e
j
3
)
)
−
−
−
l
ij
l
j
ij
3
θ θ
i
j
cos(
sin(
ϕ
ij
ϕ
ij
+
e
+
e
j
j
)
3
)
3
     (18) 
The  error  system  above  can  be  regarded  as  tracking 
controller.  It  can  keep  the  relative  distances  and  angles 
between the follower and the leader. In order to calculate 
the  system  error  illustrated  as  (18),  it  is  necessary  to 
ϕ ,  and  it  is  apparent 
calculate  the  derivatives  of  ijl  and ij
that  the  desired  separation  distance 
 and  the  desired 
ijϕ  are  constant.  Consider  the  two  robots 
bearing  angle 
formation as shown in Figure1. The x and y components of  
ijl  can be defined as 
l
ijx
l
ijy
−
d
 cos
−
d
 sin
x
jfront
y
−
θ
i
−
θ
i
x
irear
y
irear
                 (19) 
x
i
y
i
=
=
−
−
=
=
x
y
jfront
ijl
d
d
j
j
After  the  derivation  for  (19),  it  becomes  as  formula 
•
(20): 
l
ijx
•
l
ijy
(20) 
=
=
v
i
v
i
cos
sin
θ
i
θ
i
−
v
j
cos
−
v
j
sin
 
j
d
+
sin
θ ω θ
j
j
θ ω θ
j
cos
+
d
j
j
 
             
Where 
=
l
ij
l
2
ijx
+
l
2
ijy
   ,
ϕ
ij
=
− +
θ π
i
,  we  can 
l
arctan ijy
l
ijx
give  derivations  for  the  relative  distance  and  angle.  It  is 
similar to the kinematic equation described as formula (2): 
j
+
•
 =
l
ij
•
ϕ
ij
=
v
j
−
cos
1 { sin
l
ij
v
i
j
d
−
sin
γ ω γ
j
γ ω γ ω ϕ ω
}
i
ϕ ω ϕ
ij
ij
 cos
v
i
−
 sin
j
−
ϕ
ij
cos
cos
sin
i
+
l
ij
−
+
d
d
d
v
ij
j
j
j
i
   
 
 (21)
i
j
j
j
j
j
3
1
2
1
3
3
3
)
)
. 
+
je
v
i
+
γ ϕ=
j
ij
 (22) 
− +
v
j
ω
e
j
e
cos
j
3
v
e
sin
i
= −
−
+
ω ω ϕ
e
e
l
sin(
d
d
j
i ij
j
ij
j
2
3
+
−
+
ω ω
ϕ
d
l
e
cos(
d
d
i ij
ij
−
ω ω
j
Where 
Now, by performing derivation for (18), and integrating 
formula (21), applying simple trigonometric identities, the 
dynamic errors could be got by formula (22): 
•
e
j
•
e
j
•
e
We  use  the  l ϕ−  method  to  realize  the  formation 
control. Only the velocity of the robot is taken as the control 
input. We assume that the robot's velocity is not much high, 
so when we calculate the control input, we can ignore the 
dynamic  effects.  The  main  objective  of  this  paper  is  to 
study  the  effectiveness  of  multi-robot  formation  keeping 
and  obstacle  avoidance  control  in  the  leader-follower 
framework. So we ignore the influence of the factors such 
as  nonlinear  and  disturbances  in  the  process  of  research. 
Therefore, the appropriate velocity as the required control 
input can be calculated by using kinematic controller. 
4378
2015 27th Chinese Control and Decision Conference (CCDC)
 
1.5
1
0.5
y
0
-0.5
-1
-1.5
-1.5
r
o
r
r
e
 
x
0.5
0
-0.5
0
r
o
r
r
e
 
y
1
0
-1
0
r
o
r
r
e
 
n
o
i
t
a
t
n
e
i
r
o
0.1
0
-0.1
0
-1
-0.5
0
x
0.5
1
1.5
Fig5. Trajectory tracking for circular 
2
4
6
8
2
4
6
8
2
4
6
8
10
time
10
time
10
time
12
14
16
18
20
12
14
16
18
20
12
14
16
18
20
Fig6. Position error of trajectory tracking for circular 
 
4.2  Simulation  of  Robots  Formation  and  Obstacle 
Avoidance   
For the obstacle avoidance simulation, the leader-follower 
formation is combined with artificial potential field, and the 
leading  robot  is  regarded  as  the  only  robot  in  artificial 
potential  field.  The  other  robots  follow  the  leader,  and 
march with formation. The experiment region is a rectangle 
area of 30*40. The point coordinates of obstacles are [22, 
15], [25, 13], [21, 10], [20, 13] respectively, and their radius 
is 1. The coordinate of the target point is [35, 25], its radius 
is 1. The next position of the leader is determined by the 
artificial potential field method. In the region of obstacle, 
followers will follow the trajectory of the leader with the 
column to pass this area. When leaving the obstacle region, 
the followers will recover their original formation to march 
again. Figure 6 and Figure 7 are the simulation results with 
two followers and tour followers respectively. 
30
25
20
y
15
10
5
0
 
5
Obstacle Avoidance
the leader robot
the follower robot
 
Goal
F3
F1
F2
leader
10
15
20
25
30
35
40
x
 
Fig7. Obstacle avoidance with three robots 
30
25
20
y
15
10
5
0
 
5
Obstacle Avoidance
the leader robot
the follower robot
 
Goal
F3
F1
F2
F4
leader
10
15
20
25
30
35
40
x
Fig8. Obstacle avoidance with five robots 
 
         
             
 
As  shown  above,  after 
forming 
formation  of 
multi-robots, when the leader robot detects the obstacles in 
the  process  of  navigation  with  formation,  the  repulsive 
potential field force will be formed between the leader and 
the  obstacles,  and  considering  the  gravitational  potential 
field  force  between  the  leader  and  the  target  point,  the 
leader will avoid the obstacles to navigate to the target point 
under 
the  action  of  resultant  potential  field  force. 
According  to  the  trajectory  of  the  leader  robot,  the 
formation  is  changed  correspondingly  to  avoid  obstacles. 
Once through the obstacle region, the robots will regain the 
formation and continue to voyage to the target point. 
5 CONCLUSIONS 
Based  on  the  kinematics  model  of  robot,  the  paper 
establishes  the  new  system  dynamic  model  of  robot 
formation.  The 
feedback 
linearization of inputs and outputs was designed. It can be 
used to transform the tracking control problem into control 
problem  of  specific  error  system.  Combined  with  the 
obstacle avoidance algorithm of artificial potential field, we 
can successfully realize the function of formation tracking 
and obstacle avoidance for multi-robots. 
REFERENCES 
tracking  controller  with 
2015 27th Chinese Control and Decision Conference (CCDC)
4379
 
 
[8]  X. Wu, Z. Feng, J Zhu.  A novel method for formation control of 
multiple  autonomous  underwater  vehicles.  Ship  Science  and 
Technology. Vol.30, No.2,  128-134, 2008 
[9]  R.  Cui,  D.  Xu,  W.  Yan.  Leader-Follower  Formation  Control  of 
Autonomous Underwater Vehicles Using Position Information Only. 
Acta Armamentarii, Vol.29, No.8, 980-984, 2008 
[10] J. P. Desai, J. P. Ostrowski, and V. Kumar, "Controlling formations 
of multiple mobile robots," in Proceedings 1998 IEEE International 
Conference  on  Robotics  and  Automation,  Leuven,  Belgium,  pp. 
2864-2869, 1998. 
[11] O.  Khatib,  "Real-Time  Obstacle  Avoidance  for  Manipulators  and 
International 
Mobile  Robots", 
Conference on Robotics and Automation,  pp. 500-505, 1985 
in  Proceedings  1985 
[12] R.  Fierro,  and  F.L.  Lewis,  "Control  of  a  Nonholonomic  Mobile 
Robot:  Backstepping  Kinematics  Into  Dynamics",  Proc.  IEEE 
Conference on Decision Control., Kobe, Japan,  pp. 1722-1727,1996 
IEEE 
 
 
[1]  K. Yu, G. Xu, Z. Xiao, et al. Current Status and Development Trend 
of  Multi-AUV  Cooperation  System.  Ship and  Ocean  Engineering. 
Vol.38, No.5, 134-137, 2009. 
[2]  T.  Arai,  E.  Pagello  and  L.  E.  Parker,  "Editorial:  Advances  in 
Multi-Robot  Systems",  IEEE  Transactions  on  Robotics  and 
Automation, Vol. 18, No.5, pp. 655-661,2002 
[3]  Yamaguchi  A,Fukuehi  M,  Ota  J,et  al.  Motion  planning  for 
cooperative  transportation  of  a  large  object  by  multiple  mobile 
robots in a 3d environment. Lenven, Belgium: Proceedings of IEEE 
International Conference on Robotics & Automation, 2000. 
[4]  Baleh  T,  Arkin  R  C.  Behavior-based  formation  control  for 
multi-robot teams. IEEE Transactions on Robotics & Automation, 
Vol.14, No.6, 926-939, 1998. 
[5]  Sugihara K, Suzuki I. Distributed motion coordination of multiple 
the  IEEE 
mobile  robots.  Philadelphia,  PA:  Proceedings  of 
International Symposium on Intelligent Control, 1990. 
[6]  R.M.  Murray,  Recent 
in  cooperative  control  of 
multi-vehicle  systems,  J.  Dyn.  Syst.  Meas.  Control-Trans.  ASME 
Vol.129, No.5, 571-583, 2007. 
research 
[7]  Lewis M A, Tan K H. High precision formation control of mobile 
robots  using  virtual  structures.  Autonomous  Robots,  Vol.4,  No.4,   
387-403, 1997. 
 
4380
2015 27th Chinese Control and Decision Conference (CCDC)