# Design Hybrid Evolutionary PSO Aiding Miniature Aerial Vehicle Controllers

Authors: Juing-Shian Chiou; Huu-Khoa Tran; Shou-Tao Peng
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

Based on some sort of simplified fuzzy reasoning methods and PID parameters, many fuzzy-PID controller schemes are applied to control the complicated systems, recently. Mathematical optimization calculations are making in encouraging EP with PSO yielding a higher-quality solution. In this paper, a novel HEPSO algorithm as an improved variant of stochastic optimization strategy PSO, which assigns optimization to fuzzy-PID control gains, is established. The benefit integration design of the HEPSO algorithm structure is generating and updating the new parameters for the fuzzy-PID control schemes in the short setting and operation time. The proposed controllers have demonstrated performance to MAV models.


---

## Body

## 1. Introduction

The production systems are difficult in representing their actual behavior as well as in optimizing and planning their performance. For this reason, simulation modeling methods have been used extensively to mimic the behaviors and to provide a flexible platform for optimal parameter design and system improvement. Various optimization methods, based on soft computing, for instance, evolutionary algorithm, swarm intelligence, and ant colony optimization, have been explored and applied to solve the problem in consideration of criterion and attained the sensible results.Introduced by Lawrence J. Fogel in the USA for the first time, 1966 [1], evolutionary programming (EP), with continuing development by his son D. B. Fogel, is parallel to global search techniques that share the concepts of natural genetics and evolution theory. Similarly to genetic programming, which allows evolving numerical parameters, however, the EP program optimization structure is fixed. Its main variation operator is mutation and selection strategy whereas neglect the crossover operation.Inspired by social behavior of bird flocks and fish swarms, particle swarm optimization (PSO) was developed by Kennedy and Eberhart [2] in 1995. PSO is a swarm intelligence technique, founded on stochastic and population-based adaptive optimization. The primary PSO is effortless with simply a few adjusted parameters; however, when comparing with other methods such as neural network, machine learning, and genetic computation, the better performances in computing speed, accuracy, and small memory size are achieved.Different hybrid PSOs have been proposed due to computing optimal objective function, recently. HGAPSO, which is suggested by Juang [3], incorporates GA’s evolutionary operations of elite strategy, where the elites are stored in upper-half of the best-performing individuals in a population. For estimating parameters of nonlinear systems, an effective hybrid optimization strategy by incorporating the jumping property of simulated annealing (SA) into PSO, namely, PSOSA, was presented by Li et al. [4]. A parallel algorithm of PSO and differential evolution (DE) executed in parallel to enhance the population with frequent information sharing was demonstrated by Niu and Li [5]. The combination of PSO and gravitational search algorithm (GSA) was introduced by Mirjalili and Hashim based on gravitation force [6]. A novel hybrid binary particle swarm optimization algorithm HBPSO by combination of the PSO’s concept and GA was presented by Menhas et al. [7]. Cuckoo search (CS) and particle swarm optimization (PSO), inspired by nature metaheuristic algorithm for solving optimization problems, were established by Ghodrati and Lotfi [8]. HPSOM, to integrate the PSO with genetic algorithm mutation method, was introduced by Esmin and Matwin [9].Evolutionary programming (EP) algorithm, which only has a mutation operator and an individual selection strategy, neglects a crossover operator. All of the revealed methods above can increase the convergence speed. This research paper exploits the HEPSO optimization algorithm and implements it to MAV models.Because of the simple control structure, low cost, and ease to design, conventional PID controllers are applied to many industrial fields. Moreover, due to regulated control gains, a new hybrid control approach fuzzy-PID is applied for an autonomous minihelicopter [10]. Since information about the physical process is vague and data available is scarce, fuzzy rule-based model is suitable. The fuzzy logic controller (FLC) has been proved to be one of the most efficient and systematic approaches to deal with such kinds of problems in the fact that its control capability arises from emulating human logic instead of accurate mathematical model [11]. One drawback of designed fuzzy system is that the extracted rules are independent of the membership functions, especially for a complex system problem with a large number of input variables. Therefore no guarantee that the fuzzy system obtained will have sufficiently good performance. In many cases, the system’s performance, however, can be improved by further tuning the membership functions and selecting suitable fuzzification and defuzzification methods [10, 12–14].This paper deals with the design of triangular membership functions of fuzzy-PID controller. The parameter gains are then evaluated on the HEPSO by using fitness functions in association with the system’s performance indices—integral of time multiplied by absolute error (ITAE) [15]. The controller’s designs assign a task to MAV models: rotor speed, servo motor position, and altitude control to show the benefits in enhancing performances.


## 2. Fuzzy-PID Controller Description

PID controllers’ main functions are improving the dynamic response and reducing or eliminating the steady-state errors. In this paper, an improvement to such widely used fuzzy-PID controllers by incorporating a fast learning PID controller with parameters yielding a high-quality solution is proposed.The tracking errore
(
t
) and the differential tracking error d
e
(
t
), as shown in Figure 1, are the inputs of the fuzzy inference system [16]. The formal triangular membership functions are used for the two inputs e
(
t
) and d
e
(
t
) and the output C
I.Figure 1
Inputse and d
e and output C
I of fuzzification interface.The dynamic behavior of fuzzy logic controllers (FLC) is described by a set of linguistic rules, based on expert knowledge. Mandani’s MIN-MAX inference engine type and center of area method (COA) defuzzification are employed in this paper. The sliding range of ZE limited the error range and the optimal system condition regulated the PID controller.Since its combination yields the basic implementation parameters of the fuzzy control algorithm, seven partitions are decomposed: negative big (NB), negative medium (NM), negative small (NS), zero (ZE), positive small (PS), positive medium (PM), and positive big (PB) partitions. The fuzzy control rules are then designed as shown in Table1.Table 1
The fuzzy rule table.
C
I
(
t
)
e
(
t
)
NB
NM
NS
ZE
PS
PM
PB
d
e
(
t
)
NB
ZE
NS
NM
NM
NB
NB
NB
NM
PS
ZE
NS
NM
NM
NB
NB
NS
PM
PS
ZE
NS
NM
NM
NB
ZE
PM
PM
PS
ZE
NS
NM
NM
PS
PB
PM
PM
PS
ZE
NS
NM
PM
PB
PB
PM
PM
PS
ZE
NS
PB
PB
PB
PB
PM
PM
PS
ZEMost physical empirical systems have inherently intractable characteristics such as high order and nonlinearity. Thus, the satisfied performance requirement obtained by the fuzzy-PID controllers has been tackled in many studies [10, 12, 13, 17, 18]. In a fuzzy-PID structure, the endeavor controller can be written as
(1)
u
=
C
I
+
u
PID
=
C
I
+
k
P
×
e
+
k
I
∫
e
+
k
D
×
d
e
d
t
.


## 3. Hybrid Evolutionary PSO

The primary inspiration of fundamental PSO algorithm that a great number of particles move around in a multidimensional space was the social behavior of animals, such as the schooling of fish, the flocking of birds, and the swarm theory [2]. The fundamental PSO algorithm can be verified by [2–9, 11, 19, 20]
(2)
v
i
D
t
+
1
=
ω
×
v
i
D
t
+
c
1
×
ran
d
1
·
×
p
i
D
t
-
x
i
D
t
+
c
2
×
ran
d
2
(
·
)
×
p
g
D
t
-
x
i
D
t
(3)
x
i
D
t
+
1
=
x
i
D
t
+
v
i
D
t
+
1
,
where v
i
D is the velocity of particle i and x
i
D is the position of particle i with objective value fitness and p
i
D is the best historical position of particle i itself, while p
g
D is the global best position. The inertia weight ω is a user-defined that manages the previous values of particle velocities in relation to the current value, along with c
1 and c
2 as learning rates. The items ran
d
1
· and ran
d
2
· are uniformly distributed random numbers between 0 and 1. The c
1
×
ran
d
1
·
×
p
i
D
t
-
x
i
D
t term of 2 refers to the cognitive component, which reflects the distance at which the best solution P
i
(
t
) of a particle is located. The distance of particle is verified on its own.Although evolutionary programming is similar to genetic programming, while its numerical parameters are allowed to evolve, the program structure to be optimized is fixed. The different point between EP and GA is that EP neglects crossover operator while it keeps mutation operator and selection strategy. The mutation is the main varying operator; the population members are viewed as a part of specific species rather than the same species members; therefore, each parent generates a new offspring [1].Mutate everyx
i
(
t
) and generate X
j
(
t
)
=
[
x
j
1
(
t
)
,
x
j
2
(
t
)
,
…
,
x
j
D
(
t
)
] by
(4)
x
j
D
t
+
1
=
x
j
D
t
+
N
0,1
,
where x
j
D is the jth individual error data and N
(
0,1
) is a normal distribution, when a random variable x is distributed normally with mean of zero and one. The mutation rate is noted as mr.The new individual is selected by using the roulette wheel technique for every particlej. Select x
j
D from the set of all X
j; then, update v
i and x
i by 2 and 3.In this paper, by combination of PSO and EP, the proposed algorithm creates an elite HEPSO. These results usually achieve the better one than either the PSO or the existing algorithms alone.The optimization issue is to determine a group of[
k
P
,
k
I
,
k
D
,
k
e
,
k
d
e
]
∈
f
(
E
); therefore, the performance index of f
(
E
) is the minima. After determining the parameters of f
(
E
) for all particles, each new particle was said to represent a group of solutions.The progression can be assigned as illustrated in Figure2.Figure 2
The MAV controllers were aided by HEPSO algorithm.The fitness function integral of time multiplied by absolute error (ITAE), in this research, integrated four items such as the rise time, settling time, peak time, and maximum overshoot:(5)
f
E
=
α
1
×
Rise
Time
+
α
2
×
Settling
Time
+
α
3
×
Peak
Time
+
α
4
×
r
-
Overshoot
,
where α
i are the weighting factors. The rising time, settling time, peak time, and maximum overshoot are estimated via the output performance, and then its values are recorded. Afterward, the particle groups which contain large error can be eliminated. Thus, the convergence speed of the system also accelerated. The fitness function is settled in range: f
E
∈
[
0,500
].The error behaviors of transient response and steady-state error occur at the same time, which enables the designer to offer the systematic optimization ideas and to differentiate the control system quality. Due to evaluation of the performance quality of the control system, some specific performance indicators, which combine transient response and steady-state error in the control system, are usually designed to determine error criterion. The ITAE performance index [15] is mathematically given by
(6)
ITAE
=
∫
0
∞
t
·
e
t
·
d
t
.


## 4. Miniature Aerial Vehicle Model

The MAV model is developed based on the Newton-Euler formulation. The MAV model as shown in Figure3, which is derived from [16, 21–23], is depicted using right hand generalized earth coordinate systems of axes and body frame. The two front rotors (rotors 1 and 2) point towards positive x-axis, toward right rotor 2 direct y-axis and positive directed z-axis is downwards. The right handed rotation of x-, y-, and z-axis sign arrow of three positive angle variables of roll (
φ
), pitch (
θ
), and yaw (
ψ
), respectively. The coordinate y
-
z-axis gauges the tilt angle μ.Figure 3
The single tilting MAV model.The tilting MAV has six degrees of freedom (DOF). It has four inputs: three rotorcraft speeds and one tilt angle. The angle position is driven by one servo motor. Since two front rotors are maneuvered in different speed, the roll(
φ
) control is created such as when speed of rotor 1 is down and rotor 2 is up and makes the MAV toward the left and vice versa. The third tail rotor generates the pitch (
θ
) control when it changes velocity. The tilting angle μ which created yaw (
ψ
) control is driven by the small servo motor.


## 5. Simulation Results

The HEPSO algorithm parameters given by the number of population generations for speed control, servo angle position control, and altitude control are set on 30, 20, and 50 generations, respectively. The fitness weighting factor[
α
1
,
α
2
,
α
3
,
α
4
]
=
[
30
,
25
,
10
,
5
] and PID in range are k
PID
∈
[
0
,
100
]. The mutation rate is mr
=
0.1. The sampling time in this simulation is 0.01.We compared two different optimization methods which have been applied to fuzzy-PID controller gains: PSO versus HEPSO. The simulation results show that the convergence of the parameters in hybrid algorithm is yielding a higher-quality solution.The transfer functions applied to controllers are as follows.(i)
Speed control (for three DC rotorcrafts):11.4286
/
(
s
+
57.14
)
(
s
+
2.004
) (Figure 4).
(ii)
Servo motor position control (for the tilting angle DC servo motor with attached on the third axis):106.7
/
s
(
s
+
281.8
)
(
s
+
2.582
) (Figure 5).
(iii)
Altitude control (forZ motion of MAV): 1
/
(
0.84
s
2
+
0.4286
+
19.376
) (Figure 6).The results of two controller strategies in comparison are examined: the first is PSO optimized fuzzy-PID controller and the proposed controller is HEPSO optimized algorithm aiding fuzzy-PID controller. The created generation of minimum HEPSO algorithm fitness function is faster with more than 5 generations; thus, it can save the operation time. Both of them have shown the effectiveness of the controllers design scheme.Figure 4
DC rotor speed control. (a) The fitness functions of the proposed controller design. (b) The convergence of ITAE error. (c) The controllers of DC rotor speed control response.
(a)
(b)
(c)Figure 5
Servo motor position control. (a) The fitness functions of the proposed controller design. (b) The convergence of ITAE error. (c) The controllers of servo motor position control response.
(a)
(b)
(c)Figure 6
The altitude control. (a) The fitness functions of the proposed controller design. (b) The convergence of ITAE error. (c) The controllers of altitude control response.
(a)
(b)
(c)


## 6. Conclusions

In this paper, two optimization methods have been applied to MAV system models. The proposed HEPSO optimization algorithm aiding fuzzy-PID controllers have been implemented to altitude control and the MAV devices: DC motor speed control and DC servo motor position control. As obviously shown in the results, in which two controller configurations were tested, the most powerful proposed controllers display the best performance in terms of saving the settling time, reliability, and stability, as well. Moreover, when comparing two different optimization methods, PSO versus HEPSO, the convergence of the parameters in hybrid algorithm is yielding a higher-quality solution. The generation evolution of minimum HEPSO algorithm fitness function is faster; thus, it can save the setting and operation time. The effectiveness of the controllers design scheme for MAV is achieved.



---

*Source: 684094-2015-10-12.xml*