# Trafficability Analysis at Traffic Crossing and Parameters Optimization Based on Particle Swarm Optimization Method

Authors: Bin He; Qiang Lu
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

In city traffic, it is important to improve transportation efficiency and the spacing of platoon should be shortened when crossing the street. The best method to deal with this problem is automatic control of vehicles. In this paper, a mathematical model is established for the platoon’s longitudinal movement. A systematic analysis of longitudinal control law is presented for the platoon of vehicles. However, the parameter calibration for the platoon model is relatively difficult because the platoon model is complex and the parameters are coupled with each other. In this paper, the particle swarm optimization method is introduced to effectively optimize the parameters of platoon. The proposed method effectively finds the optimal parameters based on simulations and makes the spacing of platoon shorter.


---

## Body

## 1. Introduction

With the rapid development of economy and growth of population, the number of vehicles becomes more and more in China and the increase in vehicle number is unprecedented. As a result, the traffic congestion has been becoming a serious problem. Therefore, it attracts many researchers’ attention to improve the platoon’s control which can shorten the spacing of platoon and increase traffic flow.The best method to deal with this problem is automatic control of vehicles. One of the transportation automatic control systems is the Automated Highway System (AHS) which includes the longitudinal control, the lateral control, and the comprehensive control. An AHS is a proposed intelligent transportation system technology designed to provide driverless cars with specific rights-of-way. It is most often touted as a means of traffic congestion relief, as it would drastically reduce following distances and headway, thus allowing more cars to occupy a given stretch of road [1]. It belongs to longitudinal control to cross the traffic intersection at a fast speed. The aim of longitudinal control is to ensure a safe distance and maintain a relatively stable speed between vehicles.The idea of longitudinal vehicle control has developed very quickly and become very attractive with the increasing issues of traffic congestion and road safety. The researchers showed many technical considerations in the design of longitudinal control systems, such as external forces, process and measurement noise, and sampling and quantization of measurements [2]. Many models which include linear model [3, 4] and nonlinear model [5, 6] had been established and the researchers had presented systematic analysis of longitudinal control for the platoon. The scholars applied many controllers, such as classical proportional integral controller [7, 8] and the fuzzy controller [9] in longitudinal control of platoon. Wang et al. [10] analyzed the local stability of platoon and studied the relationship between the safe headway and the global stability of platoon. However, there are many parameters to be set for the platoon. When the global optimization method is applied to find the optimum parameters of platoon, the method evaluates the fitness of application and numerous fitness evaluations are needed. As such, the convergence is also an important factor in the selection of a method for preventing a platoon from numerous iterations of the method.The particle swarm optimization (PSO) method is a population based stochastic optimization method proposed by Kennedy and Eberhart in 1995 and is inspired by social behavior such as flocks of birds or schools of fish [11]. The main advantages of PSO method are simple to understand, easy to implement, and quick in convergence [12]. The PSO method has been widely studied. Pan et al. [13] analyzed the kinetic characteristic of three models of PSO method. Ren and Wang [14] proposed an accelerated convergence particle swarm optimization algorithm based on analyzing the convergence of basal PSO method. Coelho [15] presented a novel quantum-behaved PSO method using chaotic mutation operator. Chen and Zhao [16] proposed a particle swarm optimization method that uses an adaptive variable population size and periodic partial increasing or declining individuals in the form of ladder function. The PSO method has been successfully applied in the CPG parameter optimization [17], continuous nonlinear function optimization [11], reactive power and voltage control [18], and parameter tuning of controller for a power system [19].The remainder of this paper is organized as follows. Section2 builds a dynamic platoon model and presents a systematic analysis of a longitudinal control law for the platoon. Section 3 shows comprehensive and detailed analysis on parameters optimization based on the PSO method. Section 4 shows the simulation. The conclusion and some ideas about further research are discussed in Section 5.


## 2. Dynamic Model of Platoon

Various models for vehicles dynamics have been used in the study of longitudinal control of platoon. For a platoon travelling at a constant speed in a fixed direction, we adopt the following third-order model [2]. The state equation of the vehicle is given by
(1)
x
i
˙
=
v
i
,
v
i
˙
=
a
i
,
a
˙
=
1
τ
i
(
c
i
-
a
i
)
,
where x
i, v
i, and a
i denote the absolute position, velocity, and acceleration of the vehicle, respectively. The parameter τ
i is the time constant of the vehicle propulsion system and the parameter c
i is the control input. The assumed configuration for a platoon of five vehicles is shown in Figure 1.Figure 1
Platoon of five vehicles.The parameterδ
d is assigned position of vehicle along the road and the parameter δ
i is the deviation of the ith vehicle position from its assigned position. And it is described by
(2)
δ
i
=
x
i
-
1
-
x
i
-
δ
d
.Given the direction of platoon from right to left, the platoon variables are the velocityv
0 of the lead vehicle, the velocity v
i of the ith vehicle, the acceleration a
0 of the lead vehicle and the acceleration a
i of the ith vehicle.From the input/output point of view of theith linearized model [3, 4], we can obtain
(3)
x
i
⃛
=
c
i
.The objective of longitudinal control is to maintain the spacing error below a predetermined level or, if possible, at zero. As to the first vehicle, the controller [6] can be written as
(4)
c
1
=
c
x
1
δ
1
+
c
v
1
δ
1
˙
+
c
a
1
δ
1
¨
+
c
v
L
1
(
v
0
-
v
L
)
+
c
a
L
1
(
a
0
-
a
L
)
,
where c
x
1, c
v
1, c
a
1, c
v
L
1, and c
a
L
1 are design constants [6]. The parameter v
L denotes the steady-state velocity of the lead vehicle and the parameter a
L is the steady-state acceleration of the lead vehicle.Whenw
0
=
v
0
-
v
L and w
0
˙
=
a
0
-
a
L, we can obtain
(5)
c
1
=
c
x
1
δ
1
+
c
v
1
δ
1
˙
+
c
a
1
δ
1
¨
+
c
v
L
1
w
0
+
c
a
L
1
w
0
˙
.Differentiating both sides of (2) three times with respect to the time variable, we obtain
(6)
δ
⃛
1
=
x
0
⃛
-
x
1
⃛
.Taking Laplace transforms, we obtain(7)
h
δ
1
ω
0
(
s
)
=
s
2
-
c
a
L
1
s
-
c
v
L
1
s
3
+
c
a
1
s
2
+
c
v
1
s
+
c
x
1
.For theith vehicle, we can also obtain
(8)
c
i
=
c
x
δ
i
+
c
v
δ
i
˙
+
c
a
δ
i
¨
+
c
v
L
(
v
0
-
v
i
)
+
c
a
L
(
a
0
-
a
i
)
,
where c
x, c
v, c
a, c
v
L, and c
a
L are design constants [6]. Differentiating both sides of (2) three times with respect to the time variable, we obtain
(9)
δ
i
⃛
=
x
i
-
1
⃛
-
x
i
⃛
.Taking Laplace transforms, we obtain(10)
h
δ
i
δ
i
-
1
(
s
)
=
c
a
s
2
+
c
v
s
+
c
x
s
3
+
(
c
a
+
c
a
L
)
s
2
+
(
c
v
+
c
v
L
)
s
+
c
x
.For the second vehicle, the transfer function is(11)
h
δ
2
δ
1
(
s
)
=
(
c
a
1
-
c
a
L
)
s
2
+
(
c
v
1
-
c
v
L
)
s
+
c
x
1
s
3
+
(
c
a
+
c
a
L
)
s
2
+
(
c
v
+
c
v
L
)
s
+
c
x
.Therefore, the main design objective for the longitudinal control law is shown in Figure2, where X
(
s
)
=
(
c
a
L
1
s
+
c
v
L
1
)
/
(
s
3
+
(
c
a
+
c
a
L
)
s
2
+
(
c
v
+
c
v
L
)
s
+
c
x
).Figure 2
Block diagram for the platoon.We use the block diagram in Figure2 to analyze the platoon. Some considerations suggest the main design objective for the longitudinal control law.(
1
) In order to keep the platoon stable, all poles of the transfer function are required to be in the left half plane. According to Routh-Hurwitz stability criterion, we obtain c
a
1
>
0, c
x
1
>
0, c
a
1
c
v
1
>
c
x
1, and (
c
a
+
c
a
L
)
(
c
v
+
c
v
L
)
>
c
x.(
2
) Since the perturbations in the lead vehicle’s velocity from its steady-state value should not get magnified from one vehicle to the next as one goes down the platoon, we require that |
δ
i
(
j
w
)
/
δ
i
-
1
(
j
w
)
|
<
1 for all w
>
0.Therefore, these parametersc
x
1, c
v
1, c
a
1, c
v
L
1, c
a
L
1, c
x, c
v, c
a, c
v
L, and c
a
L are needed to be selected correctly.


## 3. Particle Swarm Optimization Method

The particle swarm optimization method uses the concept called particle and swarm [17, 20]. The particles correspond to an animal, bird, and insect in a herd, flock, and swarm, respectively. Let us consider a swarm including m particles which are seeking the optimum value of the objective function in an n-dimensional search space, each particle having a vector of position Y
i
=
(
y
i
1
,
y
i
2
,
…
,
y
i
n
), i
=
1,2
,
…
,
m which is associated with a solution, a vector of velocity U
i
=
(
u
i
1
,
u
i
2
,
…
,
u
i
n
), i
=
1
,
2
,
…
,
m which determines the movement value of a particle in each dimension to improve its current position, and a vector of particle best position P
best
=
(
p
best
1
,
p
best
2
,
…
,
p
best
n
) which is associated with most fitted positions of a particle from the first step of the algorithm. It is notable that the fitness of a position can easily be calculated considering the objective function of the optimization problem. A vector of the global best particle G
best
=
(
G
best
1
,
G
best
2
,
…
,
G
best
n
) is reserved for knowledge sharing among all particles of a swarm [21]. Using the aforementioned notations, each argument of the velocity and position vector for each particle in the swarm is updated through the multiple iterations of the algorithm using the following model:
(12)
U
i
(
t
+
1
)
=
b
U
i
(
t
)
+
φ
1
(
P
best
-
Y
i
(
t
)
)
+
φ
2
(
G
best
-
Y
i
(
t
)
)
,
i
=
1,2
,
…
,
m
Y
i
(
t
+
1
)
=
Y
i
(
t
)
+
U
i
(
t
+
1
)
,
i
=
1,2
,
…
,
m
,
where U
i
(
t
) is a velocity and Y
i
(
t
) is a position of the particle at t iteration. P
best is a previous best position and G
best is a global best position of each particle obtained so far. φ
1 and φ
2 are determined as φ
1
=
rand
(
0
,
k
1
) and φ
2
=
rand
(
0
,
k
2
). k
1 is a cognition learning factor and k
2 is a social learning factor. b is the inertia weight which determines the particle speed prior to the current velocity and thus functions as a balancing algorithm between global search and local search capacities.


## 4. Simulations

### 4.1. Simulation Scheme

In this paper, the PSO method is used to optimize the parametersc
v
L
1 and c
a
L
1. In order to reduce the particle dimension and avoid complicated analysis, we choose c
a
1
=
c
a
+
c
a
L, c
v
1
=
c
v
+
c
v
L, and c
x
1
=
c
x. Then (7), (10), and (11) have the same poles. The values of parameters are set as c
a
1
=
15, c
v
1
=
74, c
x
1
=
120, c
a
=
5, c
v
=
49, c
x
=
120, c
a
L
=
10, and c
v
L
=
25 [6]. The spacing between these vehicles is selected as fitness function. To examine the behavior of the platoon with optimized parameters, the simulation for platoon consisting of 5 vehicles is run in MATLAB. In the simulation, all the vehicles are assumed to be initially travelling at the velocity of 0 m/s. When t
=
0 s, the lead vehicle’s velocity increases from initial value and the acceleration is 2 m/s2. Finally it reaches the value of 20 m/s. Taking inverse Laplace transforms for (7), the fitness function is obtained:
(13)
F
(
c
a
L
1
,
c
v
L
1
,
t
)
=
(
2
c
a
L
1
-
c
v
L
1
/
2
+
8
e
4
t
+
3
c
a
L
1
-
c
v
L
1
/
2
+
18
e
6
t
-
5
c
a
L
1
-
c
v
L
1
+
25
e
5
t
)
×
20
,
where t
=
10 s.The size of the swarm is 50 and the number of iterations is 300. Parametersc
v
L
1 and c
a
L
1 are varied in the interval [
-
3.1,0
] in step of 0.01. The parameters k
1 and k
2 are set as k
1
=
0.5 and k
2
=
2.5. The inertia weight [22] is described by
(14)
b
=
b
max
⁡
-
b
max
⁡
-
b
min
⁡
M
j
,
where b
max
⁡ and b
min
⁡ are the maximal and minimum values of the inertia weight, and they are set as b
max
⁡
=
1.2 and b
min
⁡
=
0.4. M is the maximal iteration time and j is current iteration time.The optimization procedure of the PSO method is given by the following steps.Step 1.
Initialize a population of particles with random positions and velocities.Step 2.
Evaluate each particle’s fitness value.Step 3.
Compare each particle’s fitness with the particle’sP
best. If the current value is better than P
best, then set the P
best value equal to the current value and the P
best location equal to the current location.Step 4.
Compare the fitness with the population’s overall previous bestG
best. If the current value is better than G
best, then reset the G
best value to the current particle’s array index and value.Step 5.
Update each particle’s velocity and position according to (12).Step 6.
Return to Step2 until the maximum number of iterations is reached.


### 4.2. Simulations

Through the particle swarm optimization method, we obtain the values of parameters which arec
a
L
1
=
-
3.0315 and c
v
L
1
=
-
0.0492. The curve of G
best is obtained, as shown in Figure 3.Figure 3
Curve ofG
best.Then transfer functions are described by(15)
h
δ
1
ω
0
(
s
)
=
s
2
+
3.0315
s
+
0.0492
s
3
+
15
s
2
+
74
s
+
120
,
h
δ
2
δ
1
(
s
)
=
h
δ
3
δ
2
(
s
)
=
h
δ
4
δ
3
(
s
)
=
5
s
2
+
49
s
+
120
s
3
+
15
s
2
+
74
s
+
120
.The Bode Diagrams of two transfer functions are shown in Figure4. We can see that the frequence response meets the second condition of the platoon.Bode diagrams.
(a)
(b)According to the longitudinal control law, the diagram of platoon is shown in Figure5.Diagrams of the platoon.
(a)
Spacing of the platoon
(b)
Position of vehicles
(c)
Velocity of vehicles
(d)
Acceleration of vehiclesThe simulation results show that the deviations of vehicles from their preassigned positions do not exceed 0.06 m. The accelerations of vehicles in the platoon are within the range of acceptable comfort limits.


## 4.1. Simulation Scheme

In this paper, the PSO method is used to optimize the parametersc
v
L
1 and c
a
L
1. In order to reduce the particle dimension and avoid complicated analysis, we choose c
a
1
=
c
a
+
c
a
L, c
v
1
=
c
v
+
c
v
L, and c
x
1
=
c
x. Then (7), (10), and (11) have the same poles. The values of parameters are set as c
a
1
=
15, c
v
1
=
74, c
x
1
=
120, c
a
=
5, c
v
=
49, c
x
=
120, c
a
L
=
10, and c
v
L
=
25 [6]. The spacing between these vehicles is selected as fitness function. To examine the behavior of the platoon with optimized parameters, the simulation for platoon consisting of 5 vehicles is run in MATLAB. In the simulation, all the vehicles are assumed to be initially travelling at the velocity of 0 m/s. When t
=
0 s, the lead vehicle’s velocity increases from initial value and the acceleration is 2 m/s2. Finally it reaches the value of 20 m/s. Taking inverse Laplace transforms for (7), the fitness function is obtained:
(13)
F
(
c
a
L
1
,
c
v
L
1
,
t
)
=
(
2
c
a
L
1
-
c
v
L
1
/
2
+
8
e
4
t
+
3
c
a
L
1
-
c
v
L
1
/
2
+
18
e
6
t
-
5
c
a
L
1
-
c
v
L
1
+
25
e
5
t
)
×
20
,
where t
=
10 s.The size of the swarm is 50 and the number of iterations is 300. Parametersc
v
L
1 and c
a
L
1 are varied in the interval [
-
3.1,0
] in step of 0.01. The parameters k
1 and k
2 are set as k
1
=
0.5 and k
2
=
2.5. The inertia weight [22] is described by
(14)
b
=
b
max
⁡
-
b
max
⁡
-
b
min
⁡
M
j
,
where b
max
⁡ and b
min
⁡ are the maximal and minimum values of the inertia weight, and they are set as b
max
⁡
=
1.2 and b
min
⁡
=
0.4. M is the maximal iteration time and j is current iteration time.The optimization procedure of the PSO method is given by the following steps.Step 1.
Initialize a population of particles with random positions and velocities.Step 2.
Evaluate each particle’s fitness value.Step 3.
Compare each particle’s fitness with the particle’sP
best. If the current value is better than P
best, then set the P
best value equal to the current value and the P
best location equal to the current location.Step 4.
Compare the fitness with the population’s overall previous bestG
best. If the current value is better than G
best, then reset the G
best value to the current particle’s array index and value.Step 5.
Update each particle’s velocity and position according to (12).Step 6.
Return to Step2 until the maximum number of iterations is reached.


## 4.2. Simulations

Through the particle swarm optimization method, we obtain the values of parameters which arec
a
L
1
=
-
3.0315 and c
v
L
1
=
-
0.0492. The curve of G
best is obtained, as shown in Figure 3.Figure 3
Curve ofG
best.Then transfer functions are described by(15)
h
δ
1
ω
0
(
s
)
=
s
2
+
3.0315
s
+
0.0492
s
3
+
15
s
2
+
74
s
+
120
,
h
δ
2
δ
1
(
s
)
=
h
δ
3
δ
2
(
s
)
=
h
δ
4
δ
3
(
s
)
=
5
s
2
+
49
s
+
120
s
3
+
15
s
2
+
74
s
+
120
.The Bode Diagrams of two transfer functions are shown in Figure4. We can see that the frequence response meets the second condition of the platoon.Bode diagrams.
(a)
(b)According to the longitudinal control law, the diagram of platoon is shown in Figure5.Diagrams of the platoon.
(a)
Spacing of the platoon
(b)
Position of vehicles
(c)
Velocity of vehicles
(d)
Acceleration of vehiclesThe simulation results show that the deviations of vehicles from their preassigned positions do not exceed 0.06 m. The accelerations of vehicles in the platoon are within the range of acceptable comfort limits.


## 5. Conclusion

In this paper, a mathematical model is built for the platoon’s longitudinal movement and a longitudinal control law is analyzed in detail. It is well known that the parameter calibration of the platoon is a difficult problem. However, the PSO method effectively finds the parameters by storing the previous knowledge of particles and estimating the best positions and achieves the computation time of 2.7 s. The Automated Highway System is seen as a better way to deal with the problem of the traffic congestion. In general, the road conditions are complex. A number of studies should be done to analyze the effects of disturbances and modeling errors which may need further investigation.



---

*Source: 974398-2014-03-03.xml*