# Cloud Particles Evolution Algorithm

Authors: Wei Li; Lei Wang; Qiaoyong Jiang; Xinhong Hei; Bin Wang
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

Many evolutionary algorithms have been paid attention to by the researchers and have been applied to solve optimization problems. This paper presents a new optimization method called cloud particles evolution algorithm (CPEA) to solve optimization problems based on cloud formation process and phase transformation of natural substance. The cloud is assumed to have three states in the proposed algorithm. Gaseous state represents the global exploration. Liquid state represents the intermediate process from the global exploration to the local exploitation. Solid state represents the local exploitation. The cloud is composed of descript and independent particles in this algorithm. The cloud particles use phase transformation of three states to realize the global exploration and the local exploitation in the optimization process. Moreover, the cloud particles not only realize the survival of the fittest through competition mechanism but also ensure the diversity of the cloud particles by reciprocity mechanism. The effectiveness of the algorithm is validated upon different benchmark problems. The proposed algorithm is compared with a number of other well-known optimization algorithms, and the experimental results show that cloud particles evolution algorithm has a higher efficiency than some other algorithms.


---

## Body

## 1. Introduction

Many real-world problems which are classified as global optimization problems are very complex and are quite difficult to solve. Various optimization algorithms are developed to solve these problems. Most of these algorithms learn and imitate a variety of intelligent behavior in nature. These algorithms exploit a set of potential solutions and detect the optimal solution through cooperation and competition among the individuals of the population [1]. For example, Genetic algorithm [2] proposed by Professor Holland is a computational model of biological evolution which imitates natural selection and genetic mechanism. Ant colony algorithm [3] is a heuristic bionic optimization algorithm which imitates the foraging behavior of the ant colony. In the ant colony algorithm, an individual is called an artificial ant. The artificial ants search probabilistically in the solution space to create candidate solutions. The evaluating and the updating of these candidate solutions are based on the pheromone associated with each one of them. The candidate solution with the maximum amount of pheromone is considered to be the optimal solution of the problem. Particle swarm optimization (PSO) [4] is a bionic intelligent computing method proposed by Kennedy and Eberhart which imitates the flying and the foraging behavior of birds. In PSO, the optimal solution is called a leader. The leader guides each particle to search. Besides, the particles cooperate with each other to realize exploration. These behaviors are directed by a combination of swarm’s previous best (gBest) and their own previous best (pBest) [5]. The artificial bee colony algorithm (ABC) [6, 7] is a nonderivative optimization method based on self-organizational models of bee and swarm intelligence. In ABC, bees are divided into three groups which are called “employed” bees, “onlooker” bees, and “scouts” bees. “Onlookers” bees look for food sources. Food source is the search goal of bees. “Employed” bees work on food source. In addition, “scouts” bees search for the new food source near the hive. Artificial immune system (AIS) [8] is an adaptive system inspired by the immunology which imitates human immunological functions, principles, and models to solve complex problems. In AIS, the external organisms are called antigen and the internal organisms are called antibodies. Immune clonal operation provides multistrategy conditions for recombination and mutation [9]. The clonal selection selects the excellent antibodies from the subpopulation generated by clonal proliferation.This paper which is inspired by the natural phenomena of cloud formation presents a novel optimization algorithm to solve optimization problems called cloud particles evolution algorithm (CPEA). The proposed algorithm uses phase transformation mechanism to imitate the formation and diversification of the cloud to solve the optimization problems. It is well known that the species evolution is influenced by individual fitness, living environment, competition, and cooperation between the individuals. Therefore, the algorithm introduces competition and reciprocal mechanism in the evolution process so as to strengthen the viability of the cloud particles. The performance of the cloud particles evolution algorithm (CPEA) is tested on the classical optimization problems. The results are compared with other algorithms in terms of best function value and the number of function evaluations.The remainder of this paper is organized as follows. In Section2, the proposed CPEA and the concepts are introduced in detail. In Section 3, the performance of the proposed algorithm is validated by different optimization problems. The conclusion and future research are drawn in Section 4.


## 2. Cloud Particles Evolution Algorithm

### 2.1. Cloud Formation

The idea of the proposed CPEA is inspired from cloud formation. The state of the cloud particles’ changing process [10–12] is shown in Figure 1.Figure 1
Cloud particles transformation process.After the wet air rises, it will go through a set of changes because of some certain reasons temperature and pressure. Then, the water vapor is formed. Water vapor plays an important role in this process. When the water vapor is saturated, it will adsorb on the cloud condensation nuclei (CNN) and the initial cloud particles are formed. Cloud particles are condensed and sublimated continuously by absorbing water vapor around the cloud. In the process of the condensation, the cloud particles are getting closer and colliding. Then the larger cloud particles are formed. Ice phase is generated when the temperature is below 0°C at this time. The emergence of ice crystals damages the stable status of cloud phase structure. Ice crystals grow up quickly to be the snow crystals because of the condensation of the water vapor. On the one hand, water vapor transfers to ice crystals and ice crystals begin to grow up as supercooled water droplets evaporate. On the other hand, the supercooled water droplets will be frozen after the colliding to make the snow particles grow bigger. If the supercooled water vapor which has been collided is too many, the snow crystals will be transferred to be solid substances. Then, the solid substances such as snow and hail will be dropped. If the solid substances descent in the area which has a higher temperature over 0°C, they will melt into raindrop. The water in our planet will be irradiated by the sunshine, and it will become the water vapor again in the air. The cloud is formed and changed like this.


### 2.2. Phase and Phase Transformation

There are various kinds of substance in nature. Substances commonly exist in three states: solid, liquid, and gaseous. Figure2 is a simple diagram to show how water transforms among the three states. Evaporation, condensation, and melting are the three ways for water transforming. In order to describe the different forms of the substance, the phase can be used as a symbol for states of substances. In thermodynamics, phase is defined as some parts which have the same quality [13]. A phase transformation is the transition which is from one phase to another in thermodynamics.Figure 2
Three phases of water.Phase transformation is relevant to Gibbs free energy. In thermodynamics, the Gibbs free energy [14] is a thermodynamic potential that measures the “useful” or process-initiating work obtainable from a thermodynamic system at a constant temperature and pressure.


### 2.3. Population Evolution Mechanism

Darwin pointed out that mutation and natural selection are the two important factors for the evolution, which explained many natural phenomena successfully [15]. Nowak from Harvard University ranked the cooperation as the third important factor of the species evolution. Evolution is a unity of the cooperation and the opposition [16]. Every gene, cell, and organism of the individual should strengthen its own evolution process at the expense of beating its rivals. Therefore, individuals often conflict with each other because of limited resources [17]. However, the complex life organizational structure in different levels is produced in the process of evolution, which indicates that the cooperation exists in the evolution. Therefore, Nowak summarizes five mechanisms for the cooperation evolution such as kin selection, direct reciprocity, indirect reciprocity, network reciprocity, and group selection (shown in Figure 3).Figure 3
Five mechanisms of cooperation evolution [16].Haldane believes that the altruism exists in the evolution [18]. According to the study of the social insect, Hamilton [19] used mathematic model for quantitative description and defined the Hamilton rules. The rule states natural selection favors cooperation if the donor and the recipient have kin relationships. Direct reciprocity requires several repeated encounters between two individuals which can offer help with each other. The direct reciprocity theory framework is repeated Prisoner’s Dilemma. The classical game strategy is tit for tat, TFT [20]. The interactions among individuals in nature are often asymmetrical and fleeting. Therefore, indirect reciprocity is a more prevalent reciprocity form. Indirect reciprocity refers to establishing a good reputation by helping others, which will be rewarded by others later. In Indirect reciprocity, one individual acts as donor while the other acts as recipient. Indirect reciprocity will favor the cooperation evolution, if q; the probability of one individual’s reputation exceeds the cost-to-benefit ratio of the altruistic act, namely, q
>
c
/
b [21]. Network reciprocity is a new form of reciprocity. Cooperators can form network clusters in order to help each other. The cooperation mechanism of this spatial reciprocity is called network reciprocity. Group selection is a minimalist stochastic model [22]. The population is divided into several groups. Cooperators help each other in their own groups while defectors do not help. The pure cooperator groups grow faster than the pure defector groups while the defector groups reproduce faster than the cooperator groups in any mixed group. Therefore, defectors are often selected in the lower level (within groups) while cooperators are often selected in the higher (between groups) level.


### 2.4. Description of the Proposed Algorithm

The main goal of the cloud particles evolution system is to realize the evolution of the cloud particles through phase transformation. In cloud particles evolution model, cloud particles have three phases: cloud gaseous phase, cloud liquid phase, and cloud solid phase. Cloud gaseous phase indicates gaseous cloud particles which are distributed in the entire search space. Cloud liquid phase refers to the rain in nature. Cloud solid phase represents hailstones in nature.The proposed algorithm which is similar to other metaheuristic algorithms begins with an initial population called the cloud particles. At the beginning, let the initial state be cloud gaseous phase which is composed of many cloud particles (shown in Figure4). The best individual (best cloud particle) is chosen as a seed or nucleus. Then, the cloud particles are going to have a condensational growth or collisions and coalescence growth (shown in Figure 5). If the population evolves successfully, the cloud particles are going to change from cloud gaseous state to liquid state by condensation. Cloud changes from gaseous phase to cloud liquid phase. This process is called phase transformation. When cloud particles are in liquid state, cloud particles transform from global exploration to local exploitation gradually. In this process, according to the idea of the indirect reciprocity mentioned before, if q
>
c
/
b the algorithm will favor the cooperation evolution or it will change from the liquid phase to solid phase by phase transformation. When cloud particles are in solid state (shown in Figure 6), the cloud particles condense into ice crystals at the beginning, and then ice crystals curdle into hailstones. The hailstone is the best solution.Figure 4
Cloud gaseous phase.Figure 5
Cloud liquid phase.Figure 6
Cloud solid phase.It can be seen from the biological evolutionary behavior that the individuals compete to survive. However, competitive ability is always asymmetry because the different individual has different fitness. Usually, it is easier for disadvantaged individuals to compete with the individuals in the same situation while they seldom compete with the advantaged individuals [23]. Therefore, competition and cooperation which are two different survival strategies for living organisms are formed in order to adapt to the environment. In the CPEA model, according to the natural phenomena that the population develops cooperatively when they evolve, these two strategies can be shown as follows. The populations carry out the competitive evolution in cloud gaseous phase. In cloud liquid phase, the populations employ competitive mechanism for improving local exploitation ability and employ reciprocal mechanism to ensure the diversity of the population. Cloud solid phase indicates the algorithm has found the optimal solution area; therefore the algorithm realizes fast convergence by competitive evolution. The optimization process of CPEA model is shown in Figure 7.Figure 7
Optimization process of CPEA model.In CPEA, the initialization, phase transformation driving force, condensation operation, reciprocity operation, and solidification operation are described as follows.


#### 2.4.1. Initialization

When the metaheuristic methods are used to solve the optimization problem, one population is formed at first. Then, the variables which are involved in the optimization problem can be represented as an array. In this algorithm, one single solution is called one “cloud particle.” The cloud particles have three parameters [24]. They are Expectation (E
x), Entropy (E
n), and Hyperentropy (H
e). E
x is the expectation values of the distribution for all cloud particles in the domain. E
n is the range of domain which can be accepted by linguistic values (qualitative concept). In another word, E
n is ambiguity. H
e is the dispersion degree of entropy (E
n). That is to say, H
e is the entropy of entropy. In addition, H
e can be also defined as the uncertainty measure for entropy. In m dimensional optimization problems, a cloud particle is an array of 1
×
m. This array is defined as follows:
(1)
Cloud
particle
=
x
1
,
x
2
,
…
,
x
m
.To start the optimization algorithm, a matrix of cloud particles which size isN
×
m is generated by cloud generator (i.e., population of cloud particles). Cloud particles generation is described as follows: Φ
(
μ
,
δ
) is the normal random variable which has an expectation μ and a variance δ; N is the number of cloud particle; X is the cloud particle generated by cloud generator. Consider
(2)
S
=
{
s
i
∣
Φ
(
E
n
,
H
e
)
,
i
=
1
⋯
N
}
,
X
=
{
x
i
∣
Φ
(
E
x
,
s
i
)
,
s
i
∈
S
,
i
=
1
⋯
N
}
.The matrixX generated by cloud generator is given as (rows and column are the number of population and the number of design variables, resp.)(3)
Population
of
cloud
particles
=
Cloudparticl
e
1
Cloudparticl
e
2
⋮
Cloudparticl
e
N
=
X
1
1
X
2
1
X
3
1
⋯
X
m
1
X
1
2
X
2
2
X
3
2
⋯
X
m
2
⋮
⋮
⋮
⋮
⋮
X
1
N
X
2
N
X
3
N
⋯
X
m
N
.Each of the decision variable values (x
1
,
x
2
,
…
,
x
m) can represent real values for continuous or discrete problems. The fitness of a cloud particle is calculated by the evaluation of objective function (F) given as
(4)
F
i
=
f
x
1
1
,
x
2
1
,
…
,
x
i
1
i
=
1,2
,
3
,
…
,
m
.In the function above,N and m are the number of cloud particles and the number of decision variables. The cloud particles which number is N are generated by the cloud generator. The cloud particle which has the optimal value is selected as the nucleus. The total population consists of N
p subpopulations. Each subpopulation of cloud particles selects a cloud particle as the nucleus. Therefore, there are N
p nuclei.


#### 2.4.2. Phase Transformation Driving Force

The phase transformation driving force (PT) can reflect the evolution extent of the population. In the proposed algorithm, phase transformation driving force is the factor of deciding phase transform. During the evolution, the energy possessed by a cloud particle is called its evolution potential. Phase transformation driving force is the evolution potential energy difference between the cloud particles of new generation and the cloud particles of old generation. The phase transformation driving force is defined as follows:(5)
EP
⁡
i
=
b
e
s
t
i
n
d
i
v
i
d
u
a
l
i
i
=
1,2
,
3
,
…
,
g
e
n
(6)
PT
=
EP
⁡
genindex
-
1
-
EP
⁡
-
2
+
E
P
genindex
-
EP
⁡
-
2
2
EP
⁡
-
=
EP
⁡
genindex
-
1
+
EP
⁡
genindex
2
.The cloud evolution system realizes phase transformation ifPT
≤
δ. EP represents the evolution potential of cloud particles.gen represents the maximum generation of the algorithm.bestindividual represents the best solution of the current generation. genindex is the number of the current evolution.


#### 2.4.3. Condensation Operation

The population realizes the global exploration in cloud gaseous phase by condensation operation. The population takes the best individual as the nucleus in the process of the condensation operation. In the cloud gaseous phase, the condensation growth space of cloud particle is calculated by the following equations. cd is the condensation factor.generation is the number of current evolution. E
n
i is the entropy of the ith subpopulation. H
e
i is hyperentropy of the ith subpopulation. Consider
(7)
E
n
i
=
n
u
c
l
e
u
s
i
g
e
n
e
r
a
t
i
o
n
,
(8)
H
e
i
=
E
n
i
cd
.


#### 2.4.4. Reciprocity Operation

If the evolution potential of cloud particles between the two generations is less than a given threshold in cloud liquid phase, the algorithm will determine whether to employ the reciprocal evolution. If the reciprocal evolution condition is met, the reciprocal evolution is employed. Otherwise the algorithm carries out phase transformation. In reciprocity operation, the probability of one individual’s reputation is calculated byu. u
∈
[
0
,
1
] is a uniformly random number. The cost-to-benefit ratio of the altruistic act is (
1
-
0.618
)
≈
0
.
38 and 0.618 is the golden ratio. N
p represents the number of the subpopulations. n
u
c
l
e
u
s
i is the nucleus of the ith subpopulation. Consider
(9)
n
u
c
l
e
u
s
N
p
=
∑
i
=
1
N
p
-
1
n
u
c
l
e
u
s
i
N
p
-
1
.Generally, the optima value will act as the nucleus. However, the nucleus of the last subpopulation is the mean of the optima values coming from the rest of the subpopulations in reciprocity operation. The purpose of the reciprocity mechanism is to ensure the diversity of the cloud particles. For example, circle represents the location of the global optimum and pentagram represents the location of the local optima in Figure8. Figure 8(a) displays the position of nucleus before reciprocity operator. Figure 8(b) shows the individuals generated by cloud generator based on Figure 8(a). Figure 8(c) shows the position of nucleus with reciprocity operator based on Figure 8(a). Figure 8(d) shows the individuals generated by cloud generator based on Figure 8(c). By comparing Figure 8(b) with Figure 8(d), it can be seen that the distribution of the individuals in Figure 8(d) is more widely than those in Figure 8(b).Reciprocity operator.
(a)
(b)
(c)
(d)


#### 2.4.5. Solidification Operation

The algorithm realizes convergence operation in cloud solid phase by solidification operation. sf is the solidification factor:(10)
E
n
i
j
=
n
u
c
l
e
u
s
i
j
sf
,
if
n
u
c
l
e
u
s
i
j
>
10
-
2
,
n
u
c
l
e
u
s
i
j
,
otherwise
,
(11)
H
e
i
j
=
E
n
i
j
100
,
where i
=
1,2
,
…
,
N
p and j
=
1,2
,
…
,
D. N
p represents the number of the subpopulations. D represents the number of dimensions.The pseudocode of CPEA is illustrated in Algorithms1, 2, 3, and 4. N
p represents the number of the subpopulation.Childnum represents the number of individuals of every subpopulation. The search space is [
-
r
,
r
].Algorithm 1:Evolutionary algorithm of the cloud gaseous state.
(1)    if State == 0 then
(2)        Calculate the evolution potential of cloud particles according to (5)
(3)        Calculate the phase transformation driving force according to (6)
(4)        if PT ≤δ then
(5)% Phase Transformation from gaseous to liquid
(6)N
p
=
N
p
/
2
(7)Childnum = Childnum × 2
(8)            Select the nucleus of each sub-population
(9)δ
=
δ
/
10
(10)          State = 1
(11)      else
(12)% Condensation Operation
(13)          fori
=
1: N
P
(14)              UpdateE
n
i for each sub-population according to (7)
(15)              UpdateH
e
i for each sub-population according to (8)
(16)          endfor
(17)      endif
(18)  endif
State = 0: Cloud Gaseous Phase,  State = 1: Cloud Liquid Phase,  State = 2: Cloud Solid PhaseAlgorithm 2:Evolutionary algorithm of the cloud liquid state.
(1)  if State == 1 then
(2)      Calculate the evolution potential of cloud particles according to (5)
(3)      Calculate the phase transformation driving force according to (6)
(4)      fori
=
1: N
p
(5)         forj
=
1: D
(6)            if|
n
u
c
l
e
u
s
i
j
|
>
1
(7)E
n
i
j
=
|
n
u
c
l
e
u
s
i
j
|
/
20
(8)            elseif|
n
u
c
l
e
u
s
i
j
|
>
0.1
(9)E
n
i
j
=
|
n
u
c
l
e
u
s
i
j
|
/
3
(10)            else
(11)E
n
i
j
=
|
n
u
c
l
e
u
s
i
j
|
(12)            endif
(13)         endfor
(14)      endfor
(15)      fori
=
1: N
p
(16)Update H
e
i for each sub-population according to (8)
(17)      endfor
(18)      if PT ≤δ then
(19)if rand
>
c
/
b then
(20)Reciprocity Operation according to (9)% Reciprocity Operation
(21)else
(22)
% Phase Transformation from liquid to solid
(23)
N
p
=
N
p
/
2
(24)Childnum = Childnum × 2
(25)Select the nucleus of each sub-population
(26)
δ
=
δ
/
100
(27)State = 2
(28)endif
(29)       endif
(30)  endifAlgorithm 3:Evolutionary algorithm of the cloud solid state.
(1)  if State == 2 then
(2)     Calculate the evolution potential of cloud particles according to (5)
(3)     Calculate the phase transformation driving force according to (6)
(4)% Solidification Operation
(5)     if PT ≤δ then
(6)          sf = sf × 2
(7)δ
=
δ
/
10
(8)     endif
(9)     fori
=
1: N
P
(10)        forj
=
1: D
(11)            UpdateE
n
i
j according to (10)
(12)            UpdateH
e
i
j according to (11)
(13)        endfor
(14)     endfor
(15)  endifAlgorithm 4:Cloud particles evolution algorithm: CPEA.
(1)  InitializeD (number of dimensions), N
p
=
5; Childnum = 30; E
n
=
2/3r; H
e
=
E
n
/1000;
(2)  fori
=
1: N
p
(3)     forj
=
1: D
(4)        ifj is even
(5)E
x
i
j
=
r
/
2
(6)        else
(7)E
x
i
j
=
-
r
/
2
(8)        end if
(9)     end for
(10)  end for
(11)  Use cloud generator to produce cloud particles
(12)  Evaluate Subpopulation
(13)  Selectnucleus
(14)  while requirements are not satisfied do
(15)Use cloud generator to produce cloud particles
(16)Evaluate Subpopulation
(17)Selectnucleus
(18)Evolutionary algorithm of the cloud gaseous state
(19)Evolutionary algorithm of the cloud liquid state
(20)Evolutionary algorithm of the cloud solid state
(21)  end while


## 2.1. Cloud Formation

The idea of the proposed CPEA is inspired from cloud formation. The state of the cloud particles’ changing process [10–12] is shown in Figure 1.Figure 1
Cloud particles transformation process.After the wet air rises, it will go through a set of changes because of some certain reasons temperature and pressure. Then, the water vapor is formed. Water vapor plays an important role in this process. When the water vapor is saturated, it will adsorb on the cloud condensation nuclei (CNN) and the initial cloud particles are formed. Cloud particles are condensed and sublimated continuously by absorbing water vapor around the cloud. In the process of the condensation, the cloud particles are getting closer and colliding. Then the larger cloud particles are formed. Ice phase is generated when the temperature is below 0°C at this time. The emergence of ice crystals damages the stable status of cloud phase structure. Ice crystals grow up quickly to be the snow crystals because of the condensation of the water vapor. On the one hand, water vapor transfers to ice crystals and ice crystals begin to grow up as supercooled water droplets evaporate. On the other hand, the supercooled water droplets will be frozen after the colliding to make the snow particles grow bigger. If the supercooled water vapor which has been collided is too many, the snow crystals will be transferred to be solid substances. Then, the solid substances such as snow and hail will be dropped. If the solid substances descent in the area which has a higher temperature over 0°C, they will melt into raindrop. The water in our planet will be irradiated by the sunshine, and it will become the water vapor again in the air. The cloud is formed and changed like this.


## 2.2. Phase and Phase Transformation

There are various kinds of substance in nature. Substances commonly exist in three states: solid, liquid, and gaseous. Figure2 is a simple diagram to show how water transforms among the three states. Evaporation, condensation, and melting are the three ways for water transforming. In order to describe the different forms of the substance, the phase can be used as a symbol for states of substances. In thermodynamics, phase is defined as some parts which have the same quality [13]. A phase transformation is the transition which is from one phase to another in thermodynamics.Figure 2
Three phases of water.Phase transformation is relevant to Gibbs free energy. In thermodynamics, the Gibbs free energy [14] is a thermodynamic potential that measures the “useful” or process-initiating work obtainable from a thermodynamic system at a constant temperature and pressure.


## 2.3. Population Evolution Mechanism

Darwin pointed out that mutation and natural selection are the two important factors for the evolution, which explained many natural phenomena successfully [15]. Nowak from Harvard University ranked the cooperation as the third important factor of the species evolution. Evolution is a unity of the cooperation and the opposition [16]. Every gene, cell, and organism of the individual should strengthen its own evolution process at the expense of beating its rivals. Therefore, individuals often conflict with each other because of limited resources [17]. However, the complex life organizational structure in different levels is produced in the process of evolution, which indicates that the cooperation exists in the evolution. Therefore, Nowak summarizes five mechanisms for the cooperation evolution such as kin selection, direct reciprocity, indirect reciprocity, network reciprocity, and group selection (shown in Figure 3).Figure 3
Five mechanisms of cooperation evolution [16].Haldane believes that the altruism exists in the evolution [18]. According to the study of the social insect, Hamilton [19] used mathematic model for quantitative description and defined the Hamilton rules. The rule states natural selection favors cooperation if the donor and the recipient have kin relationships. Direct reciprocity requires several repeated encounters between two individuals which can offer help with each other. The direct reciprocity theory framework is repeated Prisoner’s Dilemma. The classical game strategy is tit for tat, TFT [20]. The interactions among individuals in nature are often asymmetrical and fleeting. Therefore, indirect reciprocity is a more prevalent reciprocity form. Indirect reciprocity refers to establishing a good reputation by helping others, which will be rewarded by others later. In Indirect reciprocity, one individual acts as donor while the other acts as recipient. Indirect reciprocity will favor the cooperation evolution, if q; the probability of one individual’s reputation exceeds the cost-to-benefit ratio of the altruistic act, namely, q
>
c
/
b [21]. Network reciprocity is a new form of reciprocity. Cooperators can form network clusters in order to help each other. The cooperation mechanism of this spatial reciprocity is called network reciprocity. Group selection is a minimalist stochastic model [22]. The population is divided into several groups. Cooperators help each other in their own groups while defectors do not help. The pure cooperator groups grow faster than the pure defector groups while the defector groups reproduce faster than the cooperator groups in any mixed group. Therefore, defectors are often selected in the lower level (within groups) while cooperators are often selected in the higher (between groups) level.


## 2.4. Description of the Proposed Algorithm

The main goal of the cloud particles evolution system is to realize the evolution of the cloud particles through phase transformation. In cloud particles evolution model, cloud particles have three phases: cloud gaseous phase, cloud liquid phase, and cloud solid phase. Cloud gaseous phase indicates gaseous cloud particles which are distributed in the entire search space. Cloud liquid phase refers to the rain in nature. Cloud solid phase represents hailstones in nature.The proposed algorithm which is similar to other metaheuristic algorithms begins with an initial population called the cloud particles. At the beginning, let the initial state be cloud gaseous phase which is composed of many cloud particles (shown in Figure4). The best individual (best cloud particle) is chosen as a seed or nucleus. Then, the cloud particles are going to have a condensational growth or collisions and coalescence growth (shown in Figure 5). If the population evolves successfully, the cloud particles are going to change from cloud gaseous state to liquid state by condensation. Cloud changes from gaseous phase to cloud liquid phase. This process is called phase transformation. When cloud particles are in liquid state, cloud particles transform from global exploration to local exploitation gradually. In this process, according to the idea of the indirect reciprocity mentioned before, if q
>
c
/
b the algorithm will favor the cooperation evolution or it will change from the liquid phase to solid phase by phase transformation. When cloud particles are in solid state (shown in Figure 6), the cloud particles condense into ice crystals at the beginning, and then ice crystals curdle into hailstones. The hailstone is the best solution.Figure 4
Cloud gaseous phase.Figure 5
Cloud liquid phase.Figure 6
Cloud solid phase.It can be seen from the biological evolutionary behavior that the individuals compete to survive. However, competitive ability is always asymmetry because the different individual has different fitness. Usually, it is easier for disadvantaged individuals to compete with the individuals in the same situation while they seldom compete with the advantaged individuals [23]. Therefore, competition and cooperation which are two different survival strategies for living organisms are formed in order to adapt to the environment. In the CPEA model, according to the natural phenomena that the population develops cooperatively when they evolve, these two strategies can be shown as follows. The populations carry out the competitive evolution in cloud gaseous phase. In cloud liquid phase, the populations employ competitive mechanism for improving local exploitation ability and employ reciprocal mechanism to ensure the diversity of the population. Cloud solid phase indicates the algorithm has found the optimal solution area; therefore the algorithm realizes fast convergence by competitive evolution. The optimization process of CPEA model is shown in Figure 7.Figure 7
Optimization process of CPEA model.In CPEA, the initialization, phase transformation driving force, condensation operation, reciprocity operation, and solidification operation are described as follows.


### 2.4.1. Initialization

When the metaheuristic methods are used to solve the optimization problem, one population is formed at first. Then, the variables which are involved in the optimization problem can be represented as an array. In this algorithm, one single solution is called one “cloud particle.” The cloud particles have three parameters [24]. They are Expectation (E
x), Entropy (E
n), and Hyperentropy (H
e). E
x is the expectation values of the distribution for all cloud particles in the domain. E
n is the range of domain which can be accepted by linguistic values (qualitative concept). In another word, E
n is ambiguity. H
e is the dispersion degree of entropy (E
n). That is to say, H
e is the entropy of entropy. In addition, H
e can be also defined as the uncertainty measure for entropy. In m dimensional optimization problems, a cloud particle is an array of 1
×
m. This array is defined as follows:
(1)
Cloud
particle
=
x
1
,
x
2
,
…
,
x
m
.To start the optimization algorithm, a matrix of cloud particles which size isN
×
m is generated by cloud generator (i.e., population of cloud particles). Cloud particles generation is described as follows: Φ
(
μ
,
δ
) is the normal random variable which has an expectation μ and a variance δ; N is the number of cloud particle; X is the cloud particle generated by cloud generator. Consider
(2)
S
=
{
s
i
∣
Φ
(
E
n
,
H
e
)
,
i
=
1
⋯
N
}
,
X
=
{
x
i
∣
Φ
(
E
x
,
s
i
)
,
s
i
∈
S
,
i
=
1
⋯
N
}
.The matrixX generated by cloud generator is given as (rows and column are the number of population and the number of design variables, resp.)(3)
Population
of
cloud
particles
=
Cloudparticl
e
1
Cloudparticl
e
2
⋮
Cloudparticl
e
N
=
X
1
1
X
2
1
X
3
1
⋯
X
m
1
X
1
2
X
2
2
X
3
2
⋯
X
m
2
⋮
⋮
⋮
⋮
⋮
X
1
N
X
2
N
X
3
N
⋯
X
m
N
.Each of the decision variable values (x
1
,
x
2
,
…
,
x
m) can represent real values for continuous or discrete problems. The fitness of a cloud particle is calculated by the evaluation of objective function (F) given as
(4)
F
i
=
f
x
1
1
,
x
2
1
,
…
,
x
i
1
i
=
1,2
,
3
,
…
,
m
.In the function above,N and m are the number of cloud particles and the number of decision variables. The cloud particles which number is N are generated by the cloud generator. The cloud particle which has the optimal value is selected as the nucleus. The total population consists of N
p subpopulations. Each subpopulation of cloud particles selects a cloud particle as the nucleus. Therefore, there are N
p nuclei.


### 2.4.2. Phase Transformation Driving Force

The phase transformation driving force (PT) can reflect the evolution extent of the population. In the proposed algorithm, phase transformation driving force is the factor of deciding phase transform. During the evolution, the energy possessed by a cloud particle is called its evolution potential. Phase transformation driving force is the evolution potential energy difference between the cloud particles of new generation and the cloud particles of old generation. The phase transformation driving force is defined as follows:(5)
EP
⁡
i
=
b
e
s
t
i
n
d
i
v
i
d
u
a
l
i
i
=
1,2
,
3
,
…
,
g
e
n
(6)
PT
=
EP
⁡
genindex
-
1
-
EP
⁡
-
2
+
E
P
genindex
-
EP
⁡
-
2
2
EP
⁡
-
=
EP
⁡
genindex
-
1
+
EP
⁡
genindex
2
.The cloud evolution system realizes phase transformation ifPT
≤
δ. EP represents the evolution potential of cloud particles.gen represents the maximum generation of the algorithm.bestindividual represents the best solution of the current generation. genindex is the number of the current evolution.


### 2.4.3. Condensation Operation

The population realizes the global exploration in cloud gaseous phase by condensation operation. The population takes the best individual as the nucleus in the process of the condensation operation. In the cloud gaseous phase, the condensation growth space of cloud particle is calculated by the following equations. cd is the condensation factor.generation is the number of current evolution. E
n
i is the entropy of the ith subpopulation. H
e
i is hyperentropy of the ith subpopulation. Consider
(7)
E
n
i
=
n
u
c
l
e
u
s
i
g
e
n
e
r
a
t
i
o
n
,
(8)
H
e
i
=
E
n
i
cd
.


### 2.4.4. Reciprocity Operation

If the evolution potential of cloud particles between the two generations is less than a given threshold in cloud liquid phase, the algorithm will determine whether to employ the reciprocal evolution. If the reciprocal evolution condition is met, the reciprocal evolution is employed. Otherwise the algorithm carries out phase transformation. In reciprocity operation, the probability of one individual’s reputation is calculated byu. u
∈
[
0
,
1
] is a uniformly random number. The cost-to-benefit ratio of the altruistic act is (
1
-
0.618
)
≈
0
.
38 and 0.618 is the golden ratio. N
p represents the number of the subpopulations. n
u
c
l
e
u
s
i is the nucleus of the ith subpopulation. Consider
(9)
n
u
c
l
e
u
s
N
p
=
∑
i
=
1
N
p
-
1
n
u
c
l
e
u
s
i
N
p
-
1
.Generally, the optima value will act as the nucleus. However, the nucleus of the last subpopulation is the mean of the optima values coming from the rest of the subpopulations in reciprocity operation. The purpose of the reciprocity mechanism is to ensure the diversity of the cloud particles. For example, circle represents the location of the global optimum and pentagram represents the location of the local optima in Figure8. Figure 8(a) displays the position of nucleus before reciprocity operator. Figure 8(b) shows the individuals generated by cloud generator based on Figure 8(a). Figure 8(c) shows the position of nucleus with reciprocity operator based on Figure 8(a). Figure 8(d) shows the individuals generated by cloud generator based on Figure 8(c). By comparing Figure 8(b) with Figure 8(d), it can be seen that the distribution of the individuals in Figure 8(d) is more widely than those in Figure 8(b).Reciprocity operator.
(a)
(b)
(c)
(d)


### 2.4.5. Solidification Operation

The algorithm realizes convergence operation in cloud solid phase by solidification operation. sf is the solidification factor:(10)
E
n
i
j
=
n
u
c
l
e
u
s
i
j
sf
,
if
n
u
c
l
e
u
s
i
j
>
10
-
2
,
n
u
c
l
e
u
s
i
j
,
otherwise
,
(11)
H
e
i
j
=
E
n
i
j
100
,
where i
=
1,2
,
…
,
N
p and j
=
1,2
,
…
,
D. N
p represents the number of the subpopulations. D represents the number of dimensions.The pseudocode of CPEA is illustrated in Algorithms1, 2, 3, and 4. N
p represents the number of the subpopulation.Childnum represents the number of individuals of every subpopulation. The search space is [
-
r
,
r
].Algorithm 1:Evolutionary algorithm of the cloud gaseous state.
(1)    if State == 0 then
(2)        Calculate the evolution potential of cloud particles according to (5)
(3)        Calculate the phase transformation driving force according to (6)
(4)        if PT ≤δ then
(5)% Phase Transformation from gaseous to liquid
(6)N
p
=
N
p
/
2
(7)Childnum = Childnum × 2
(8)            Select the nucleus of each sub-population
(9)δ
=
δ
/
10
(10)          State = 1
(11)      else
(12)% Condensation Operation
(13)          fori
=
1: N
P
(14)              UpdateE
n
i for each sub-population according to (7)
(15)              UpdateH
e
i for each sub-population according to (8)
(16)          endfor
(17)      endif
(18)  endif
State = 0: Cloud Gaseous Phase,  State = 1: Cloud Liquid Phase,  State = 2: Cloud Solid PhaseAlgorithm 2:Evolutionary algorithm of the cloud liquid state.
(1)  if State == 1 then
(2)      Calculate the evolution potential of cloud particles according to (5)
(3)      Calculate the phase transformation driving force according to (6)
(4)      fori
=
1: N
p
(5)         forj
=
1: D
(6)            if|
n
u
c
l
e
u
s
i
j
|
>
1
(7)E
n
i
j
=
|
n
u
c
l
e
u
s
i
j
|
/
20
(8)            elseif|
n
u
c
l
e
u
s
i
j
|
>
0.1
(9)E
n
i
j
=
|
n
u
c
l
e
u
s
i
j
|
/
3
(10)            else
(11)E
n
i
j
=
|
n
u
c
l
e
u
s
i
j
|
(12)            endif
(13)         endfor
(14)      endfor
(15)      fori
=
1: N
p
(16)Update H
e
i for each sub-population according to (8)
(17)      endfor
(18)      if PT ≤δ then
(19)if rand
>
c
/
b then
(20)Reciprocity Operation according to (9)% Reciprocity Operation
(21)else
(22)
% Phase Transformation from liquid to solid
(23)
N
p
=
N
p
/
2
(24)Childnum = Childnum × 2
(25)Select the nucleus of each sub-population
(26)
δ
=
δ
/
100
(27)State = 2
(28)endif
(29)       endif
(30)  endifAlgorithm 3:Evolutionary algorithm of the cloud solid state.
(1)  if State == 2 then
(2)     Calculate the evolution potential of cloud particles according to (5)
(3)     Calculate the phase transformation driving force according to (6)
(4)% Solidification Operation
(5)     if PT ≤δ then
(6)          sf = sf × 2
(7)δ
=
δ
/
10
(8)     endif
(9)     fori
=
1: N
P
(10)        forj
=
1: D
(11)            UpdateE
n
i
j according to (10)
(12)            UpdateH
e
i
j according to (11)
(13)        endfor
(14)     endfor
(15)  endifAlgorithm 4:Cloud particles evolution algorithm: CPEA.
(1)  InitializeD (number of dimensions), N
p
=
5; Childnum = 30; E
n
=
2/3r; H
e
=
E
n
/1000;
(2)  fori
=
1: N
p
(3)     forj
=
1: D
(4)        ifj is even
(5)E
x
i
j
=
r
/
2
(6)        else
(7)E
x
i
j
=
-
r
/
2
(8)        end if
(9)     end for
(10)  end for
(11)  Use cloud generator to produce cloud particles
(12)  Evaluate Subpopulation
(13)  Selectnucleus
(14)  while requirements are not satisfied do
(15)Use cloud generator to produce cloud particles
(16)Evaluate Subpopulation
(17)Selectnucleus
(18)Evolutionary algorithm of the cloud gaseous state
(19)Evolutionary algorithm of the cloud liquid state
(20)Evolutionary algorithm of the cloud solid state
(21)  end while


## 2.4.1. Initialization

When the metaheuristic methods are used to solve the optimization problem, one population is formed at first. Then, the variables which are involved in the optimization problem can be represented as an array. In this algorithm, one single solution is called one “cloud particle.” The cloud particles have three parameters [24]. They are Expectation (E
x), Entropy (E
n), and Hyperentropy (H
e). E
x is the expectation values of the distribution for all cloud particles in the domain. E
n is the range of domain which can be accepted by linguistic values (qualitative concept). In another word, E
n is ambiguity. H
e is the dispersion degree of entropy (E
n). That is to say, H
e is the entropy of entropy. In addition, H
e can be also defined as the uncertainty measure for entropy. In m dimensional optimization problems, a cloud particle is an array of 1
×
m. This array is defined as follows:
(1)
Cloud
particle
=
x
1
,
x
2
,
…
,
x
m
.To start the optimization algorithm, a matrix of cloud particles which size isN
×
m is generated by cloud generator (i.e., population of cloud particles). Cloud particles generation is described as follows: Φ
(
μ
,
δ
) is the normal random variable which has an expectation μ and a variance δ; N is the number of cloud particle; X is the cloud particle generated by cloud generator. Consider
(2)
S
=
{
s
i
∣
Φ
(
E
n
,
H
e
)
,
i
=
1
⋯
N
}
,
X
=
{
x
i
∣
Φ
(
E
x
,
s
i
)
,
s
i
∈
S
,
i
=
1
⋯
N
}
.The matrixX generated by cloud generator is given as (rows and column are the number of population and the number of design variables, resp.)(3)
Population
of
cloud
particles
=
Cloudparticl
e
1
Cloudparticl
e
2
⋮
Cloudparticl
e
N
=
X
1
1
X
2
1
X
3
1
⋯
X
m
1
X
1
2
X
2
2
X
3
2
⋯
X
m
2
⋮
⋮
⋮
⋮
⋮
X
1
N
X
2
N
X
3
N
⋯
X
m
N
.Each of the decision variable values (x
1
,
x
2
,
…
,
x
m) can represent real values for continuous or discrete problems. The fitness of a cloud particle is calculated by the evaluation of objective function (F) given as
(4)
F
i
=
f
x
1
1
,
x
2
1
,
…
,
x
i
1
i
=
1,2
,
3
,
…
,
m
.In the function above,N and m are the number of cloud particles and the number of decision variables. The cloud particles which number is N are generated by the cloud generator. The cloud particle which has the optimal value is selected as the nucleus. The total population consists of N
p subpopulations. Each subpopulation of cloud particles selects a cloud particle as the nucleus. Therefore, there are N
p nuclei.


## 2.4.2. Phase Transformation Driving Force

The phase transformation driving force (PT) can reflect the evolution extent of the population. In the proposed algorithm, phase transformation driving force is the factor of deciding phase transform. During the evolution, the energy possessed by a cloud particle is called its evolution potential. Phase transformation driving force is the evolution potential energy difference between the cloud particles of new generation and the cloud particles of old generation. The phase transformation driving force is defined as follows:(5)
EP
⁡
i
=
b
e
s
t
i
n
d
i
v
i
d
u
a
l
i
i
=
1,2
,
3
,
…
,
g
e
n
(6)
PT
=
EP
⁡
genindex
-
1
-
EP
⁡
-
2
+
E
P
genindex
-
EP
⁡
-
2
2
EP
⁡
-
=
EP
⁡
genindex
-
1
+
EP
⁡
genindex
2
.The cloud evolution system realizes phase transformation ifPT
≤
δ. EP represents the evolution potential of cloud particles.gen represents the maximum generation of the algorithm.bestindividual represents the best solution of the current generation. genindex is the number of the current evolution.


## 2.4.3. Condensation Operation

The population realizes the global exploration in cloud gaseous phase by condensation operation. The population takes the best individual as the nucleus in the process of the condensation operation. In the cloud gaseous phase, the condensation growth space of cloud particle is calculated by the following equations. cd is the condensation factor.generation is the number of current evolution. E
n
i is the entropy of the ith subpopulation. H
e
i is hyperentropy of the ith subpopulation. Consider
(7)
E
n
i
=
n
u
c
l
e
u
s
i
g
e
n
e
r
a
t
i
o
n
,
(8)
H
e
i
=
E
n
i
cd
.


## 2.4.4. Reciprocity Operation

If the evolution potential of cloud particles between the two generations is less than a given threshold in cloud liquid phase, the algorithm will determine whether to employ the reciprocal evolution. If the reciprocal evolution condition is met, the reciprocal evolution is employed. Otherwise the algorithm carries out phase transformation. In reciprocity operation, the probability of one individual’s reputation is calculated byu. u
∈
[
0
,
1
] is a uniformly random number. The cost-to-benefit ratio of the altruistic act is (
1
-
0.618
)
≈
0
.
38 and 0.618 is the golden ratio. N
p represents the number of the subpopulations. n
u
c
l
e
u
s
i is the nucleus of the ith subpopulation. Consider
(9)
n
u
c
l
e
u
s
N
p
=
∑
i
=
1
N
p
-
1
n
u
c
l
e
u
s
i
N
p
-
1
.Generally, the optima value will act as the nucleus. However, the nucleus of the last subpopulation is the mean of the optima values coming from the rest of the subpopulations in reciprocity operation. The purpose of the reciprocity mechanism is to ensure the diversity of the cloud particles. For example, circle represents the location of the global optimum and pentagram represents the location of the local optima in Figure8. Figure 8(a) displays the position of nucleus before reciprocity operator. Figure 8(b) shows the individuals generated by cloud generator based on Figure 8(a). Figure 8(c) shows the position of nucleus with reciprocity operator based on Figure 8(a). Figure 8(d) shows the individuals generated by cloud generator based on Figure 8(c). By comparing Figure 8(b) with Figure 8(d), it can be seen that the distribution of the individuals in Figure 8(d) is more widely than those in Figure 8(b).Reciprocity operator.
(a)
(b)
(c)
(d)


## 2.4.5. Solidification Operation

The algorithm realizes convergence operation in cloud solid phase by solidification operation. sf is the solidification factor:(10)
E
n
i
j
=
n
u
c
l
e
u
s
i
j
sf
,
if
n
u
c
l
e
u
s
i
j
>
10
-
2
,
n
u
c
l
e
u
s
i
j
,
otherwise
,
(11)
H
e
i
j
=
E
n
i
j
100
,
where i
=
1,2
,
…
,
N
p and j
=
1,2
,
…
,
D. N
p represents the number of the subpopulations. D represents the number of dimensions.The pseudocode of CPEA is illustrated in Algorithms1, 2, 3, and 4. N
p represents the number of the subpopulation.Childnum represents the number of individuals of every subpopulation. The search space is [
-
r
,
r
].Algorithm 1:Evolutionary algorithm of the cloud gaseous state.
(1)    if State == 0 then
(2)        Calculate the evolution potential of cloud particles according to (5)
(3)        Calculate the phase transformation driving force according to (6)
(4)        if PT ≤δ then
(5)% Phase Transformation from gaseous to liquid
(6)N
p
=
N
p
/
2
(7)Childnum = Childnum × 2
(8)            Select the nucleus of each sub-population
(9)δ
=
δ
/
10
(10)          State = 1
(11)      else
(12)% Condensation Operation
(13)          fori
=
1: N
P
(14)              UpdateE
n
i for each sub-population according to (7)
(15)              UpdateH
e
i for each sub-population according to (8)
(16)          endfor
(17)      endif
(18)  endif
State = 0: Cloud Gaseous Phase,  State = 1: Cloud Liquid Phase,  State = 2: Cloud Solid PhaseAlgorithm 2:Evolutionary algorithm of the cloud liquid state.
(1)  if State == 1 then
(2)      Calculate the evolution potential of cloud particles according to (5)
(3)      Calculate the phase transformation driving force according to (6)
(4)      fori
=
1: N
p
(5)         forj
=
1: D
(6)            if|
n
u
c
l
e
u
s
i
j
|
>
1
(7)E
n
i
j
=
|
n
u
c
l
e
u
s
i
j
|
/
20
(8)            elseif|
n
u
c
l
e
u
s
i
j
|
>
0.1
(9)E
n
i
j
=
|
n
u
c
l
e
u
s
i
j
|
/
3
(10)            else
(11)E
n
i
j
=
|
n
u
c
l
e
u
s
i
j
|
(12)            endif
(13)         endfor
(14)      endfor
(15)      fori
=
1: N
p
(16)Update H
e
i for each sub-population according to (8)
(17)      endfor
(18)      if PT ≤δ then
(19)if rand
>
c
/
b then
(20)Reciprocity Operation according to (9)% Reciprocity Operation
(21)else
(22)
% Phase Transformation from liquid to solid
(23)
N
p
=
N
p
/
2
(24)Childnum = Childnum × 2
(25)Select the nucleus of each sub-population
(26)
δ
=
δ
/
100
(27)State = 2
(28)endif
(29)       endif
(30)  endifAlgorithm 3:Evolutionary algorithm of the cloud solid state.
(1)  if State == 2 then
(2)     Calculate the evolution potential of cloud particles according to (5)
(3)     Calculate the phase transformation driving force according to (6)
(4)% Solidification Operation
(5)     if PT ≤δ then
(6)          sf = sf × 2
(7)δ
=
δ
/
10
(8)     endif
(9)     fori
=
1: N
P
(10)        forj
=
1: D
(11)            UpdateE
n
i
j according to (10)
(12)            UpdateH
e
i
j according to (11)
(13)        endfor
(14)     endfor
(15)  endifAlgorithm 4:Cloud particles evolution algorithm: CPEA.
(1)  InitializeD (number of dimensions), N
p
=
5; Childnum = 30; E
n
=
2/3r; H
e
=
E
n
/1000;
(2)  fori
=
1: N
p
(3)     forj
=
1: D
(4)        ifj is even
(5)E
x
i
j
=
r
/
2
(6)        else
(7)E
x
i
j
=
-
r
/
2
(8)        end if
(9)     end for
(10)  end for
(11)  Use cloud generator to produce cloud particles
(12)  Evaluate Subpopulation
(13)  Selectnucleus
(14)  while requirements are not satisfied do
(15)Use cloud generator to produce cloud particles
(16)Evaluate Subpopulation
(17)Selectnucleus
(18)Evolutionary algorithm of the cloud gaseous state
(19)Evolutionary algorithm of the cloud liquid state
(20)Evolutionary algorithm of the cloud solid state
(21)  end while


## 3. Experiments and Discussions

In this section, the performance of the proposed CPEA is tested with a comprehensive set of benchmark functions which include 40 different global optimization problems and CEC2013 benchmark problems. The definition of the benchmark functions and their global optimum(s) are listed in the Appendix. The benchmark functions are chosen from Ali et al. [25] and Rahnamayan et al. [26]. They are unimodal functions, multimodal functions, and rotate functions. All the algorithms are implemented on the same machine with a T2450 2.4 GHz CPU, 1.5 GB memory, and windows XP operating system with Matlab R2009b.


### 3.1. Parameter Study

In this section, the setting basis for initial values of some parameters involved in the cloud particles evolution algorithm model is described. Firstly, we hope that all the cloud particles which are generated by the cloud generator can fall into the search space as much as possible. Secondly, the diversity of the population can be ensured effectively if the cloud particles are distributed in the search space uniformly. Therefore, some experiments are designed to analyze the feasibility of the selected scheme in order to illustrate the rationality of the selected parameters.


#### 3.1.1. Definition ofE
n

Let the search rangeS be [
-
r
,
r
] and r is the search radius. Set r
=
2, E
x
1
=
-
2
/
3, E
x
2
=
2
/
3, H
e
=
0.001
E
n. Five hundred cloud particles are generated in every experiment. The experiment is independently simulated 30 times. The cloud particle is called a valid cloud particle if it falls into the search space. Figure 9 is a graph which shows the number of the valid cloud particles generated with different E
n, and Figure 10 is a graph which shows the distribution of valid cloud particles generated with different E
n.Figure 9
Valid cloud particles of differentE
n.Figure 10
Distribution of cloud particles in differentE
n.It can be seen from Figure9 that the number of the valid cloud particles reduces as the E
n increases. The average number of valid cloud particles is basically the same when E
n is (
1
/
3
)
r, (
1
/
2
)
r, and (
2
/
3
)
r, respectively. However, the number of valid cloud particles is relatively fewer when E
n is (
3
/
4
)
r, (
4
/
5
)
r, and r. It can be seen from Figure 10 that all the cloud particles fall into the research space when E
n is (
1
/
3
)
r, (
1
/
2
)
r, and (
2
/
3
)
r, respectively. However, the cloud particles concentrate in a certain area of the search space when E
n is (
1
/
3
)
r or (
1
/
2
)
r. By contrast, the cloud particles are widely distributed throughout all the search space [
-
2
,
2
] when E
n is (
2
/
3
)
r. Therefore, set E
n
=
(
2
/
3
)
r. The dashed rectangle in Figure 10 represents the decision variables space.


#### 3.1.2. Definition ofH
e

Ackley Function (multimode) and Sphere Function (unimode) are selected to calculate the number of fitness evaluations (NFES) of the algorithm when the solution reachesf_stop with different H
e. Let the search space be [
-
32,32
] (Ackley) and [
-
100,100
] (Sphere), f
_
stop
=
10
-
3. H
e
1
=
E
n, H
e
2
=
0.1
E
n, H
e
3
=
0.01
E
n, H
e
4
=
0.001
E
n, H
e
5
=
0.0001
E
n, H
e
6
=
0.00001
E
n, and H
e
7
=
0.000001
E
n. The reported values are the average of the results for 50 independent runs. As shown in Figure 11, for Ackley Function, the number of fitness evaluations is least when H
e is equal to 0.001
E
n. For Sphere Function, the number of fitness evaluations is less when H
e is equal to 0.001
E
n. Therefore, H
e is equal to 0.001
E
n in the cloud particles evolution algorithm.NFES of CPEA in differentH
e.
(a)
(b)


#### 3.1.3. Definition ofE
x

Ackley Function (multimode) and Sphere Function (unimode) are selected to calculate the number of fitness evaluations (NFES) of the algorithm when the solution reachesf_stop with E
x being equal to r, (
1
/
2
)
r, (
1
/
3
)
r, (
1
/
4
)
r, and (
1
/
5
)
r, respectively. Let the search space be [
-
r
,
r
], f
_
stop
=
10
-
3. As shown in Figure 12, for Ackley Function, the number of fitness evaluations is least when E
x is equal to (
1
/
2
)
r. For Sphere Function, the number of fitness evaluations is less when E
x is equal to (
1
/
2
)
r or (
1
/
3
)
r. Therefore, E
x is equal to (
1
/
2
)
r in the cloud particles evolution algorithm.NFES of CPEA in differentE
x.
(a)
(b)


#### 3.1.4. Definition ofc
/
b

Many things in nature have the mathematical proportion relationship between their parts. If the substance is divided into two parts, the ratio of the whole and the larger part is equal to the ratio of the larger part to the smaller part; this ratio is defined as golden ratio, expressed algebraically(12)
a
+
b
a
=
a
b
≝
φ
,
in which φ represents the golden ration. Its value is
(13)
φ
=
1
+
5
2
=
1.6180339887
,
1
φ
=
φ
-
1
=
0.618
⋯
.The golden section has a wide application not only in the fields of painting, sculpture, music, architecture, and management but also in nature. Studies have shown that plant leaves, branches, or petals are distributed in accordance with the golden ratio in order to survive in the wind, frost, rain, and snow. It is the evolutionary result or the best solution which adapts its own growth after hundreds of millions of years of long-term evolutionary process.In this model, competition and reciprocity of the population depend on the fitness and are a dynamic relationship which constantly adjust and change. With competition and reciprocity mechanism, the convergence speed and the population diversity are improved. Premature convergence is avoided, and large search space and complex computational problems are solved. Indirect reciprocity is selected in this paper. Here, setc
/
b
=
(
1
-
0.618
)
≈
0
.
38.


### 3.2. Experiments Settings

For fair comparison, we set the parameters to be fixed. For CPEA,population
size
=
5, c
h
i
l
d
n
u
m
=
30, E
n
=
(
2
/
3
)
r, (the search space is [
-
r
,
r
]), H
e
=
E
n
/
1000, E
x
=
(
1
/
2
)
r, cd
=
1000, δ
=
10
-
2, and sf
=
1000. For CMA-ES, the parameter values are chosen from [27]. The parameters of JADE come from [28]. For PSO, the parameter values are chosen from [29]. The parameters of DE are set to be F
=
0
.
3 and CR
=
0.9 used in [30]. For ABC, the number of colony size NP
=
20 and the number of food sources FoodNumber
=
NP
/
2 and limit
=
100. A food source will not be exploited anymore and is assumed to be abandoned when limit is exceeded for the source. For CEBA, population
size
=
10, c
h
i
l
d
n
u
m
=
50, E
n
=
1, and H
e
=
E
n
/
6. The parameters of PSO-cf-Local and FDR-PSO come from [31, 32].


### 3.3. Comparison Strategies

Five comparison strategies come from the literature [33] and literature [26] to evaluate the algorithms. The comparison strategies are described as follows.Successful Run. A run during which the algorithm achieves the fixed accuracy level within the max number of function evaluations for the particular dimension is shown as follows:
(14)
SR
=
successful
runs
total
runs
.Error. The error of a solution x is defined as f
(
x
)
-
f
(
x
*
), where x
* is the global minimum of the function. The minimum error is written when the max number of function evaluations is reached in 30 runs. Both the average and standard deviation of the error values are calculated.Convergence Graphs. The convergence graphs show the median performance of the total runs with termination by either the max number of function evaluations or the termination error value.Number of Function Evaluations (NFES). The number of function evaluations is recorded when the value-to-reach is reached. The average and standard deviation of the number of function evaluations are calculated.Acceleration Rate (AR). The acceleration rate is defined as follows, based on the NFES for the two algorithms:
(15)
AR
=
NFE
S
a
1
NFE
S
a
2
,
in which AR
>
1 means Algorithm 2 (a
2) is faster than Algorithm 1 (a
1).The average acceleration rate (A
R
ave) and the average success rate (S
R
ave) over n test functions are calculated as follows:
(16)
A
R
ave
=
1
n
∑
i
=
1
n
A
R
i
,
SR
ave
=
1
n
∑
i
=
1
n
SR
i
.


### 3.4. Comparison of CMA-ES, DE, RES, LOS, ABC, CEBA, and CPEA

This section compares CMA-ES, DE, RES, LOS, ABC, CEBA [34], and CPEA in terms of mean number of function evaluations for successful runs. Results of CMA-ES, DE, RES, and LOS are taken from literature [35]. Each run is stopped and regarded as successful if the function value is smaller than f_stop.As shown in Table1, CPEA can reach f_stop faster than other algorithms with respect to f
07 and f
17. However, CPEA converges slower than CMA-ES with respect to f
18. When the dimension is 30, CPEA can reach f_stop relatively fast with respect to f
33. Nonetheless, CPEA converges slower while the dimension increases.Table 1
Average number of function evaluations to reachf_stop of CPEA versus CMA-ES, DE, RES, LOS, ABC, and CEBA on f
07, f
17, f
18, and f
33 function.
F
f_stop
init
D
CMA-ES
DE
RES
LOS
ABC
CEBA
CPEA
f
07
1
e
-
3
[
-600, 600
]
D
20
3111
8691
—
—
12463
94412
2950
30
4455
11410
2
e
5
—
29357
163545
4270
100
12796
31796
—
—
34777
458102
7260
f
17
0.9
[
-5.12, 5.12]
D
20
68586
12971
—
9.2
e
4
11183
30436
2850
DE:[
-600, 600]
D
30
147416
20150
1.0
e
5
2.3
e
5
26990
44138
6187
100
1010989
73620
—
—
449404
60854
10360
f
18
1
e
-
3
[
-30, 30]
D
20
2667
—
—
6.0
e
4
10183
64400
3330
30
3701
12481
1.1
e
5
9.3
e
4
17810
91749
4100
100
11900
36801
—
—
59423
205470
12248
f
33
0.9
[
-5.12, 5.12]
D
30
152000
1.25
e
6
—
—
—
50578
13100
100
240899
5522
—
—
—
167324
14300
The bold entities indicate the best results obtained by the algorithm.D represents dimension.


### 3.5. Comparison of CMA-ES, DE, PSO, ABC, CEBA, and CPEA

This section compares the error of a solution of CMA-ES, DE, PSO, ABC [36], CEBA, and CPEA under a given NFES. The best solution of every single function in different algorithms is highlighted with boldface.As shown in Table2, CPEA obtains better results in function f
0
1, f
0
2, f
0
4, f
0
5, f
0
7, f
0
9, f
10, f
11, and f
12 while it obtains relatively bad results in other functions. The reason for this phenomenon is that the changing rule of E
n is set according to the priori information and background knowledge instead of adaptive mechanism in CPEA. Therefore, the results of some functions are relatively worse. Besides, E
n which can be optimized in order to get better results is an important controlling parameter in the algorithm.Table 2
Comparison on the ERROR values of CPEA versus CMA-ES, DE, PSO, ABC, and CEBA on functionf
01~f
13.
F
NFES
D
CMA-ESmean/SD
DEmean/SD
PSOmean/SD
ABCmean/SD
CEBAmean/SD
CPEAmean/SD
f
01
30 k
30
5.87e − 29/1.56e − 29
3.20e − 07/2.51e − 07
5.44e − 12/1.56e − 11
9.71e − 16/2.26e − 16
1.74e − 02/1.55e − 02
3.1
9
e
-
33/
6.2
1
e
-
33
30 k
50
1.02e − 28/1.87e − 29
2.95e − 03/2.88e − 03
3.93e − 03/6.62e − 02
2.35e − 08/2.47e − 08
4.07e − 02/2.11e − 02
1.09
e
-
48/
1.41
e
-
48
f
02
10 k
30
1.50e − 05/3.66e − 05
2.30e + 00/1.78e + 00
6.24e − 06/8.74e − 05
0.00
e
+
00/
0.00
e
+
00
6.16e + 00/3.22e + 00
0.00
e
+
00/
0.00
e
+
00
50 k
50
0.00
e
+
00/
0.00
e
+
00
1.65e + 00/2.05e + 00
3.87e + 002.65e + 00
0.00
e
+
00/
0.00
e
+
00
6.08e + 00/5.34e + 00
0.00
e
+
00/
0.00
e
+
00
f
03
10 k
2
2.65e − 30/8.08e − 31
0.00
e
+
00/
0.00
e
+
00
5.01e − 14/8.52e − 14
4.83e − 04/1.40e − 03
5.83e − 04/8.32e − 04
1.52e − 13/2.44e − 13
f
04
100 k
30
8.2
8
e
-
27/
1.6
9
e
-
27
2.47e − 03/2.50e − 03
1.69e − 02/2.57e − 02
1.36e + 00/1.93e + 00
1.59e + 00/7.93e + 00
3.31e − 14/1.42e − 13
100 k
50
2.29e − 06/2.22e − 06
7.37e − 01/1.93e − 01
5.01e + 00/3.88e + 00
1.19e + 02/4.54e + 01
3.13e + 02/3.94e + 02
3.98
e
-
07/
8.79
e
-
07
f
05
30 k
30
8.19e − 09/5.08e − 09
1.27e − 12/3.94e − 12
8.66e − 34/2.88e − 33
8.00e − 13/2.24e − 12
2.54e − 09/4.13e − 09
2.11
e
-
37/
4.20
e
-
36
50 k
50
2.68e − 08/1.67e − 09
1.93e − 13/2.01e − 13
4.37
e
-
54/
7.58
e
-
54
9.54e − 12/1.55e − 11
1.34e − 15/2.01e − 15
1.07e − 23/1.85e − 23
f
06
30 k
2
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
3.60e − 09/1.05e − 08
1.1e − 03/8.2e − 04
7.71e − 09/4.52e − 09
f
07
30 k
30
1.88e − 12/1.40e − 12
3.08e − 06/5.41e − 06
1.64e − 09/2.85e − 08
1.54e − 09/2.78e − 08
2.07e − 06/5.57e − 05
5.32
e
-
16/
7.59
e
-
16
50 k
50
0.00
e
+
00/
0.00
e
+
00
3.28e − 05/5.69e − 05
1.79e − 02/2.49e − 02
1.27e − 09/2.21e − 09
9.49e − 03/4.23e − 03
0.00
e
+
00/
0.00
e
+
00
f
08
30 k
3
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
2.80e − 03/2.13e − 03
3.85e − 06/6.21e − 06
f
09
30 k
2
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
f
10
5 k
2
6.99e − 47/2.21e − 46
1.28e − 26/3.9e − 26
6.59e − 07/7.71e − 07
2.32e − 04/2.60e − 04
8.24e − 08/1.37e − 07
8.29
e
-
59/
3.42
e
-
58
f
11
30 k
30
1.23e − 04/1.82e − 04
3.74e − 01/2.91e − 01
4.78e − 01/1.94e − 01
2.39e + 00/3.27e + 00
1.42e + 00/8.13e + 00
5.70
e
-
40/
6.61
e
-
39
50 k
50
5.58e − 01/1.55e − 01
3.36e + 00/1.55e + 00
3.15e + 01/1.11e + 01
5.71e + 01/6.98e + 00
2.05e + 01/9.40e + 00
2.33
e
-
05/
3.58
e
-
05
f
12
30 k
30
3.19e − 14/1.45e − 14
1.20e − 08/2.24e − 08
2.15e − 05/5.18e − 04
3.56e − 10/1.22e − 10
1.93e − 07/6.45e − 08
2.33
e
-
23/
2.50
e
-
22
50 k
50
1.41e − 14/2.43e − 14
6.77e − 05/1.11e − 05
6.08e − 03/3.35e − 03
6.37e − 10/6.83e − 10
9.76e − 03/3.08e − 03
1.90
e
-
72/
3.29
e
-
72
f
13
30 k
30
3.94e − 11/1.15e − 10
1.34
e
-
31/3.21e − 31
1.34
e
-
31/
0.00
e
+
00
6.57e − 04/3.77e − 03
2.09e − 02/7.34e − 02
1.20e − 12/2.82e − 12
50 k
50
3.08e − 03/3.80e − 03
1.34
e
-
31/
0.00
e
+
00
1.34
e
-
31/
0.00
e
+
00
1.41e − 02/9.29e − 03
4.22e − 02/2.53e − 02
2.36e − 12/1.04e − 12
“SD” stands for standard deviation.


### 3.6. Comparison of JADE, DE, PSO, ABC, CEBA, and CPEA

This section compares JADE, DE, PSO, ABC, CEBA, and CPEA in terms of SR and NFES for most functions atD
=
30 except f
03 (D
=
2). The termination criterion is to find a value smaller than 10
-
8. The best result of different algorithms for each function is highlighted in boldface.As shown in Table3, CPEA can reach the specified precision successfully except for function f
14 which successful run is 96. In addition, CPEA can reach the specified precision faster than other algorithms in function f
0
1, f
0
2, f
0
4, f
11, f
14, and f
17. For other functions, CPEA spends relatively long time to reach the precision. The reason for this phenomenon is that, for some functions, the changing rule of E
n given in the algorithm makes CPEA to spend more time to get rid of the local optima and converge to the specified precision.Table 3
Comparison on the SR and the NFES of CPEA versus JADE, DE, PSO, ABC, and CEBA for each function.
F
f
01  mean/SD
f
02  mean/SD
f
03  mean/SD
f
04  mean/SD
f
07  mean/SD
f
11  mean/SD
f
14  mean/SD
f
16  mean/SD
f
17  mean/SD
JADE
SR
1
1
1
1
1
0.98
0.94
1
1
NFES
2.8e + 4/6.0e + 2
1.0e + 4/3.5e + 2
5.7e + 3/6.9e + 2
6.9e + 4/3.0e + 3
3
.
3
e
+
3/
2.6
e
+
2
6.6e + 4/4.2e + 3
4.0e + 4/2.1e + 3
1
.
3
e
+
3/
3.8
e
+
2
1.3e + 5/1.8e + 3
DE
SR
1
1
1
1
0.96
1
0.94
1
0
NFES
3.3e + 4/6.6e + 2
1.1e + 4/2.3e + 3
1
.
6
e
+
3/
1.1
e
+
2
1.9e + 5/1.2e + 4
3.4e + 4/1.6e + 3
1.4e + 5/1.0e + 4
2.6e + 4/1.6e + 4
2.1e + 4/8.1e + 3
—
PSO
SR
1
1
1
0
0.86
0
0
1
0
NFES
2.9e + 4/7.5e + 2
1.6e + 4/1.4e + 3
7.7e + 3/9.6e + 2
—
2.8e + 4/9.8e + 2
—
—
1.8e + 4/7.2e + 2
—
ABC
SR
1
1
1
0
0.96
0
0
1
1
NFES
1.7e + 4/1.0e + 3
6.3e + 3/1.5e + 3
1.6e + 4/2.6e + 3
—
1.9e + 4/2.3e + 3
—
—
6.5e + 3/3.9e + 2
3.9e + 4/1.5e + 4
CEBA
SR
1
1
1
0
1
1
0.96
1
1
NFES
1.3e + 5/5.8e + 3
5.9e + 4/4.6e + 3
4.2e + 4/6.1e + 3
—
2.2e + 5/1.2e + 4
1.3e + 5/5.8e + 3
2.9e + 5/1.4e + 4
5.4e + 4/3.8e + 3
1.1e + 5/7.3e + 3
CPEA
SR
1
1
1
1
1
1
0.96
1
1
NFES
4.3
e
+
3/
5
.
1
e
+
2
4.5
e
+
3/
2
.
7
e
+
2
5.7e + 3/3.0e + 3
6.1
e
+
4/
2.6
e
+
3
7.7e + 3/3.4e + 2
2.3
e
+
4/
3.4
e
+
3
3.4
e
+
4/
1.6
e
+
3
2.4e + 3/4.3e + 2
1
.
0
e
+
4/
1
.
3
e
+
3As shown in Table4, the solution obtained by CPEA is closer to the best solution compared with other algorithms except for functions f
13 and f
20. The reason is that the phase transformation mechanism introduced in the algorithm can improve the convergence of the population greatly. Besides, the reciprocity mechanism of the algorithm can improve the adaptation and the diversity of the population in case the algorithm falls into the local optimum. Therefore, the algorithm can obtain better results in these functions.Table 4
Comparison on the ERROR values of CPEA versus JADE, DE, PSO, and ABC atD
=
50 except f
21 (D
=
2) and f
22 (D
=
2) for each function.
F
NFES
JADEmean/SD
DEmean/SD
PSOmean/SD
ABCmean/SD
CPEAmean/SD
f
13
40 k
1.34
e
-
31
/
0.00
e
+
00
1.34
e
-
31
/
0.00
e
+
00
1.35e − 31/0.00
e
+
00
3.52e − 03/6.33e − 02
3.64e − 12/4.68e − 12
f
15
40 k
9.18e − 11/7.04e − 11
2.94e − 05/3.66e − 05
1.81e − 09/4.18e − 08
4.11e − 14/2.62e − 14
4.5
9
e
-
53
/
1.36
e
-
52
f
16
40 k
1.15e − 20/1.42e − 20
6.62e − 08/1.24e − 07
6.91e − 19/7.55e − 19
3.41e − 15/5.23e − 15
1.77
e
-
76
/
0.00
e
+
00
f
17
200 k
6.71e − 07/1.32e − 07
2.58e − 01/6.15e − 01
6.73e + 00/2.06e + 00
6.07e − 09/1.85e − 08
0.00
e
+
00
/
0.00
e
+
00
f
18
50 k
1.16e − 07/4.27e − 08
3.89e − 04/2.10e − 04
1.70e − 09/3.13e − 08
2.93e − 09/2.34e − 08
8.88
e
-
16
/
0.00
e
+
00
f
19
40 k
1.14e − 02/1.68e − 02
9.90e − 03/6.66e − 03
9.00e − 06/1.62e − 05
2.22e − 05/1.95e − 04
1.64
e
-
64
/
5.13
e
-
64
f
20
40 k
0.00
e
+
00
/
0.00
e
+
00
0.00
e
+
00
/
0.00
e
+
00
1.40e − 10/2.81e − 10
4.38e − 03/3.33e − 03
1.21e − 10/3.11e − 09
f
21
30 k
0.00
e
+
00
/
0.00
e
+
00
0.00
e
+
00
/
0.00
e
+
00
0.00
e
+
00
/
0.00
e
+
00
1.06e − 17/3.34e − 17
0.00
e
+
00
/
0.00
e
+
00
f
22
30 k
1.28e − 45/1.86e − 45
4.99e − 164/0.00
e
+
00
2.63e − 51/6.82e − 51
6.62e − 18/5.84e − 18
0.00
e
+
00
/
0.00
e
+
00
f
23
30 k
0.00
e
+
00
/
0.00
e
+
00
0.00
e
+
00
/
0.00
e
+
00
0.00
e
+
00
/
0.00
e
+
00
0.00
e
+
00/4.24e − 04
0.00
e
+
00
/
0.00
e
+
00As shown in Table5, CPEA can find the optimal solution for the test functions. However, the stability of CPEA is worse than JADE and DE (f
25, f
26, f
27, f
28 and f
31). This result shows that the setting of H
e needs to be further improved.Table 5
Comparison on the optimum solution of CPEA versus JADE, DE, PSO, and ABC for each function.
F
NFES
JADE mean/SD
DE mean/SD
PSO mean/SD
ABC mean/SD
CPEA mean/SD
f
08
40 k
-
3.86278
/
0.00
e
+
00
-
3.86278
/
0.00
e
+
00
-
3.86278
/
6.20
e
-
17
-
3.86278
/
0.00
e
+
00
-
3.86278
/
0.00
e
+
00
f
24
40 k
-
3.3
1044
/
3.60
e
-
02
-
3.24608
/
5.70
e
-
02
-
3.25800
/
5.90
e
-
02
-
3.33539
/
0.00
e
+
00
-
3.3353
9
/
3.03
e
-
06
f
25
40 k
-
10.1532
/
4.00
e
-
14
-
10.1532
/
4.20
e
-
16
-
5.79196
/
3.30
e
+
00
-
10
.
153
2
/
1.25
e
-
13
-
10
.
153
2
/
1.58
e
-
09
f
26
40 k
-
10.4029
/
9.40
e
-
16
-
10.4029
/
2.50
e
-
16
-
7.14128
/
3.50
e
+
00
-
10
.
4029
/
1.77
e
-
15
-
10
.
402
9
/
1.94
e
-
10
f
27
40 k
-
10.5364
/
8.10
e
-
12
-
10.5364
/
7.10
e
-
16
-
7.02213
/
3.60
e
+
00
-
10
.
536
4
/
4.03
e
-
05
-
10
.
536
4
/
2.02
e
-
09
f
28
40 k
3.00000
/
1.10
e
-
15
3.00000
/
6.4
e
-
16
3.00000
/
1.30
e
-
15
3.0000
0
/
6.17
e
-
05
3.00000
/
2
.
01
e
-
15
f
29
40 k
2.55
e
-
20
/
4.79
e
-
20
0.00
e
+
00
/
0.00
e
+
00
2.10
e
-
03
/
4.28
e
-
03
1.14
e
-
18
/
1.62
e
-
18
0.00
e
+
00
/
0.00
e
+
00
f
30
40 k
-
2.3458
/
0.00
e
+
00
-
2.3458
/
0.00
e
+
00
-
2.3458
/
0.00
e
+
00
-
2.3458
/
0.00
e
+
00
-
2.3458
/
0.00
e
+
00
f
31
40 k
4.78
e
-
31
/
1.07
e
-
30
3.59
e
-
55
/
1.13
e
-
54
1.22
e
-
36
/
3.85
e
-
36
1.40
e
-
11
/
1.38
e
-
11
8.69
e
-
25
/
1.14
e
-
24
f
32
40 k
1
.
49
e
-
02
/
4
.
23
e
-
03
3.46
e
-
02
/
5.14
e
-
03
6.58
e
-
02
/
1.22
e
-
02
2.95
e
-
01
/
1.13
e
-
01
5
.
70
e
-
03
/
3
.
58
e
-
03
Data of functionsf
08, f
24, f
25, f
26, f
27, and f
28 in JADE, DE, and PSO are from literature [28].


### 3.7. Comparison of CMAES, JADE, DE, PSO, ABC, CEBA, and CPEA

The convergence curves of different algorithms on some selected functions are plotted in Figures13 and 14. The solutions obtained by CPEA are better than other algorithms with respect to functions f
0
1, f
11, f
16, f
32, and f
33. However, the solution obtained by CPEA is worse than CMA-ES with respect to function f
0
4. The solution obtained by CPEA is worse than JADE for high dimensional problems such as function f
0
5 and function f
0
7. The reason is that the setting of initial population of CPEA is insufficient for high dimensional problems.Figure 13
Convergence graph for test functionsf
0
1, f
0
4, f
16, and f
32 (D
=
30).Figure 14
Convergence graph for test functionsf
0
5, f
11, f
33 (D
=
50), and f
0
7 (D
=
100).


### 3.8. Comparison of DE, ODE, RDE, PSO, ABC, CEBA, and CPEA

The algorithms DE, ODE, RDE, PSO, ABC, CBEA, and CPEA are compared in terms of convergence speed. The termination criterion is to find error values smaller than10
-
8. The results of solving the benchmark functions are given in Tables 6 and 7. Data of all functions in DE, ODE, and RDE are from literature [26]. The best result of NEFS and the SR for each function are highlighted in boldface. The average success rates and the average acceleration rate on test functions are shown in the last row of the tables.Table 6
Comparison of DE, ODE, RDE, PSO, ABC, CEBA, and CPEA. The best result for each case is highlighted in boldface.
F
D
DE
ODE
ARODE
RDE
ARRDE
PSO
ARPSO
ABC
ARABC
CEBA
ARCEBA
CPEA
ARCPEA
NFES
SR
NFES
SR
NFES
SR
NFES
SR
NFES
SR
NFES
SR
NFES
SR
f
0
1
30
87748
1
47716
1
1.83
115096
1
0.76
20290
1
4.32
18370
1
4.77
133381
1
0.66
5266
1
16.66
f
0
2
30
41588
1
23124
1
1.80
54316
1
0.77
11541
1
3.60
6247
1
6.66
58416
1
0.71
2966
1
14.02
f
0
3
2
4324
1
4776
1
0.90
4904
1
0.88
6110
1
0.71
18050
1
0.24
38278
1
0.11
4350
1
0.99
f
0
4
20
177880
1
168680
1
1.05
231152
1
0.77
—
0
—
—
0
—
—
0
—
11068
1
16.07
f
0
5
30
25140
1
8328
1
3.01
29096
1
0.86
5910
1
4.25
6134
1
4.10
18037
1
1.39
1700
1
14.79
f
0
6
2
8016
1
6608
1
1.21
8248
1
0.97
7185
1
1.12
25407
1
0.32
66319
1
0.12
5033
1
1.59
f
0
7
30
113428
1
69342
0.96
1.64
149744
1
0.76
20700
0.84
5.48
20930
0.96
5.42
209278
1
0.54
9941
1
11.41
f
0
8
3
3376
1
3580
1
0.94
3652
1
0.92
4090
1
0.83
1492
1
2.26
206859
1
0.02
31874
1
0.11
f
0
9
2
5352
1
4468
1
1.20
6292
1
0.85
2830
1
1.89
1077
1
4.97
9200
1
0.58
3856
1
1.39
f
10
2
3608
1
3288
1
1.09
3924
1
0.92
5485
1
0.66
10891
1
0.33
6633
1
0.54
1920
1
1.88
f
11
30
385192
1
369104
1
1.04
498200
1
0.77
—
0
—
—
0
—
142564
1
2.70
12662
1
30.42
f
12
30
187300
1
155636
1
1.20
244396
1
0.76
—
0
—
25787
1
7.23
157139
1
1.19
10183
1
18.39
f
13
30
101460
1
70408
1
1.44
138308
1
0.73
6512
1
15.58
—
0
—
—
0
—
4300
1
23.60
f
14
30
570290
0.28
72250
0.88
7.89
285108
1
2.00
—
0
—
—
0
—
292946
1
1.95
31112
0.98
18.33
f
15
30
96488
1
53304
1
1.81
126780
1
0.76
18757
1
5.14
15150
1
6.37
103475
1
0.93
6466
1
14.92
f
16
2
1016
1
996
1
1.02
1084
1
0.94
428
1
2.37
301
1
3.38
3292
1
0.31
610
1
1.67
f
17
10
328844
1
70389
0.76
4.67
501875
0.96
0.66
—
0
—
7935
1
41.44
48221
1
6.82
3725
1
88.28
f
18
30
169152
1
98296
1
1.72
222784
1
0.76
27250
0.92
6.21
32398
1
5.22
184041
1
0.92
13400
1
12.62
f
19
30
411164
1
337532
1
1.22
927230
0.24
0.44
—
0
—
156871
1
2.62
148014
1
2.78
6233
1
65.97
Ave
0.96
0.98
0.98
0.96
0.86
0.67
4.01
0.79
6.36
0.89
1.31
0.99
18.58Table 7
Comparison of DE, ODE, RDE, PSO, ABC, CEBA, and CPEA. The best result for each case is highlighted in boldface.
F
D
DE
ODE
ARODE
RDE
ARRDE
PSO
ARPSO
ABC
ARABC
CEBA
ARCEBA
CPEA
ARCPEA
NFES
SR
NFES
SR
NFES
SR
NFES
SR
NFES
SR
NFES
SR
NFES
SR
f
20
5
2163
0.88
2024
1
1.07
2904
0.96
0.74
11960
1
0.18
—
0
—
58252
1
0.04
4216
0.92
0.51
f
21
2
7976
1
5092
1
1.56
8332
1
0.96
8372
0.98
0.95
15944
1
0.50
11129
1
0.72
1025
1
7.78
f
22
2
3780
1
3396
1
1.11
3820
1
0.98
5755
1
0.66
2050
1
1.84
6635
1
0.57
855
1
4.42
f
23
10
19528
1
15704
1
1.24
23156
1
0.84
9290
1
2.10
5845
1
3.34
17907
1
1.09
1130
1
17.28
f
24
6
—
0
—
0
—
—
0
—
7966
0.94
—
3018
1
—
—
0
—
42716
0.98
—
f
25
4
—
0
—
0
—
—
0
—
—
0
—
4280
1
—
—
0
—
5471
0.94
—
f
26
4
—
0
—
0
—
—
0
—
—
0
—
8383
1
—
—
0
—
5391
0.96
—
f
27
4
4576
1
4500
1
1.02
5948
1
0.77
—
0
—
4850
1
0.94
—
0
—
5512
0.92
0.83
f
28
2
4780
1
4684
1
1.02
4924
1
0.97
6770
1
0.71
6228
1
0.77
14248
1
0.34
4495
1
1.06
f
29
2
8412
1
5772
1
1.46
10860
1
0.77
—
0
—
1914
1
4.40
37765
1
0.22
4950
1
1.70
f
30
2
3256
1
3120
1
1.04
3308
1
0.98
3275
1
0.99
1428
1
2.28
114261
1
0.03
3025
1
1.08
f
31
4
1108
1
1232
1
0.90
1388
1
0.80
990
1
1.12
9944
1
0.11
131217
1
0.01
2662
1
0.42
f
34
2
5232
1
5956
0.92
0.88
5768
1
0.91
3495
1
1.50
1852
1
2.83
9926
1
0.53
4370
1
1.20
f
35
5
38532
1
16340
1
2.36
46080
1
0.84
11895
1
3.24
25303
1
1.52
35292
1
1.09
2095
1
18.39
f
36
2
7888
1
4272
1
1.85
9148
1
0.86
10910
1
0.72
8359
1
0.94
13529
1
0.58
1340
1
5.89
f
37
2
5284
1
4728
1
1.12
5488
1
0.97
7685
1
0.69
1038
1
5.09
9056
1
0.58
1022
1
5.17
f
38
2
5280
1
4804
1
1.10
5540
1
0.95
7680
1
0.69
1598
1
3.30
8905
1
0.59
990
1
5.33
f
39
10
37824
1
24206
1
1.56
46800
1
0.80
—
0
—
—
0
—
153161
1
0.25
38420
0.94
0.98
f
40
2
2976
1
2872
1
1.03
3200
1
0.93
2790
1
1.07
666
1
4.47
9032
1
0.33
3085
1
0.96
Ave
0.84
0.84
1.27
0.84
0.88
0.73
1.12
0.89
2.31
0.77
0.46
0.98
4.56It can be seen from Tables6 and 7 that CPEA outperforms DE, ODE, RDE, PSO, ABC, and CBEA on 38 of the functions. CPEA performs better than other algorithms on 23 functions. In Table 6, the average AR is 0.98, 0.86, 4.01, 6.36, 1.31, and 18.58 for ODE, RDE, PSO, ABC, CBEA, and CPEA, respectively. In Table 7, the average AR is 1.27, 0.88, 1.12, 2.31, 0.46, and 4.56 for ODE, RDE, PSO, ABC, CBEA, and CPEA, respectively. In Table 6, the average SR is 0.96, 0.98, 0.96, 0.67, 0.79, 0.89, and 0.99 for DE, ODE, RDE, PSO, ABC, CBEA, and CPEA, respectively. In Table 7, the average SR is 0.84, 0.84, 0.84, 0.73, 0.89, 0.77, and 0.98 for DE, ODE, RDE, PSO, ABC, CBEA, and CPEA, respectively. These results clearly demonstrate that CPEA has better performance.


### 3.9. Test on the CEC2013 Benchmark Problems

To further verify the performance of CPEA, a set of recently proposed CEC2013 rotated benchmark problems are used. In the experiment,t-test [37] has been carried out to show the differences between CPEA and the other algorithm. “General merit over contender” displays the difference between the number of better results and the number of worse results, which is used to give an overall comparison between the two algorithms. Table 8 gives the information about the average error, standard deviation, and t-value of 30 runs of 7 algorithms over 300,000 FEs on 20 test functions with D
=
10. The best results among the algorithms are shown in bold.Table 8
Results for CEC2013 benchmark functions over 30 independent runs.
Function
Result
DE
PSO
ABC
CEBA
PSO-FDR
PSO-cf-Local
CPEA
F
0
1
F
mean SD
1.64e − 02(+)/3.98e − 02
1.51e − 14(≈)/5.76e − 14
1.51e − 14(≈)/5.76e − 14
4.99e + 02(+)/1.46e + 02
0
.
00
e
+
00
(
≈
)/
0
.
00
e
+
00
0
.
00
e
+
00
(
≈
)/
0
.
00
e
+
00
0.00e+00/
0
.
00
e
+
00
F
0
2
F
mean SD
1.20e + 05(+)/1.50e + 05
1.73e + 05(+)/2.02e + 05
3.31e + 06(+)/1.13e + 06
3.95e + 06(+)/1.60e + 06
4.37e + 04(+)/3.76e + 04
2.11e + 04(≈)/1.61e + 04
1.82
e
+
0
4/
1.12
e
+
0
4
F
0
3
F
mean SD
1.75e + 06(+)/6.09e + 05
4.63e + 05(+)/8.77e + 05
1.35e + 07(+)/1.30e + 07
6.79e + 08(+)/3.27e + 08
1.38e + 06(+)/1.89e + 06
2.92e + 04(+)/7.10e + 04
4.91
e
+
00/
2.52
e
+
0
1
F
0
4
F
mean SD
1.30e + 03(+)/1.43e + 03
1.97e + 02(−)/
1.11
e
+
0
2
1.16e + 04(+)/3.19e + 03
6.25e + 03(+)/2.25e + 03
1.54
e
+
0
2
(
-
)/1.60e + 02
8.18e + 02(≈)/5.26e + 02
7.26e + 02/3.26e + 02
F
0
5
F
mean SD
1.38e − 01(+)/3.62e − 01
5.30e − 14(−)/5.76e − 14
1.62e − 13(−)/5.72e − 14
1.15e + 02(+)/9.45e + 01
7.57e − 15(−)/2.88e − 14
0
.
00
e
+
00
(
-
)/
0
.
00
e
+
00
7.72e − 05/2.23e − 05
F
0
6
F
mean SD
6.30e + 00(≈)/4.47e + 00
6.21e + 00(≈)/4.36e + 00
1.51
e
+
00
(
-
)/
2.90
e
+
00
5.72e + 01(+)/1.94e + 01
5.28e + 00(≈)/4.92e + 00
9.26e + 00(+)/2.51e + 00
6.21e + 00/4.80e + 00
F
0
7
F
mean SD
7.25e − 01(+)/1.93e + 00
1.63e + 00(+)/1.89e + 00
4.08e + 01(+)/1.18e + 01
4.01e + 01(+)/5.31e + 00
1.89e + 00(+)/1.94e + 00
6.87e − 01(+)/1.23e + 00
1.69
e
-
03/
3.05
e
-
03
F
0
8
F
mean SD
2.02
e
+
01
(
≈
)/6.33e − 02
2.02
e
+
01
(
≈
)/6.88e − 02
2.04e + 01(+)/8.13e − 02
2.02
e
+
0
1
(
≈
)/5.86e − 02
2.02
e
+
0
1
(
≈
)/6.46e − 02
2.02
e
+
0
1
(
≈
)/7.73e − 02
2.02
e
+
0
1/
4.62
e
-
02
F
0
9
F
mean SD
1.03
e
+
00
(
-
)/
8.41
e
-
01
2.97e + 00(≈)/1.26e + 00
5.30e + 00(+)/8.55e − 01
6.83e + 00(+)/5.71e − 01
2.44e + 00(≈)/1.37e + 00
2.26e + 00(≈)/8.62e − 01
2.32e + 00/1.29e + 00
F
10
F
mean SD
5.64e − 02(≈)/4.11e − 02
3.88e − 01(+)/2.90e − 01
1.42e + 00(+)/3.27e − 01
6.24e + 01(+)/1.40e + 01
3.39e − 01(+)/2.25e − 01
1.12e − 01(+)/5.53e − 02
5
.
3
5
e
-
02/
2
.
81
e
-
02
F
1
1
F
mean SD
7.29e − 01(−)/8.63e − 01
5.63e − 01(−)/7.24e − 01
0
.
00
e
+
00
(
-
)/
0
.
00
e
+
00
4.56e + 01(+)/5.89e + 00
6.63e − 01(−)/6.57e − 01
1.22e + 00(−)/1.29e + 00
7.95e + 00/3.30e + 00
F
12
F
mean SD
5.31e + 00(≈)/2.27e + 00
1.38e + 01(+)/6.58e + 00
2.80e + 01(+)/8.22e + 00
4.59e + 01(+)/5.13e + 00
1.38e + 01(+)/4.80e + 00
6.36e + 00(+)/2.19e + 00
5.00
e
+
00/
1.21
e
+
00
F
13
F
mean SD
7.74
e
+
00
(
-
)/4.59e + 00
1.89e + 01(≈)/6.14e + 00
3.39e + 01(+)/9.86e + 00
4.68e + 01(+)/8.05e + 00
1.53e + 01(≈)/6.07e + 00
9.04e + 00(−)/
4.00
e
+
00
1.70e + 01/8.72e + 00
F
14
F
mean SD
2.83e + 00(−)/3.36e + 00
4.82e + 01(−)/5.11e + 01
1.46
e
-
01
(
-
)/
6.09
e
-
02
9.74e + 02(+)/1.12e + 02
5.65e + 01(−)/8.34e + 01
2.00e + 02(−)/1.00e + 02
3.17e + 02/1.53e + 02
F
15
F
mean SD
3.04e + 02(≈)/2.03e + 02
6.88e + 02(+)/2.91e + 02
7.35e + 02(+)/1.94e + 02
7.21e + 02(+)/1.65e + 02
6.38e + 02(+)/1.92e + 02
3.90e + 02(+)/1.24e + 02
3.02
e
+
0
2/
1.15
e
+
0
2
F
16
F
mean SD
9.56e − 01(+)1.42e − 01
7.13e − 01(+)/1.86e − 01
8.93e − 01(+)/1.91e − 01
1.02e + 00(+)/1.14e − 01
7.60e − 01(+)/2.38e − 01
5.88e − 01(+)/2.59e − 01
4.47
e
-
02/
4.10
e
-
02
F
17
F
mean SD
1.45e + 01(≈)/4.11e + 00
8.74e + 00(−)/4.71e + 00
8.55
e
+
00
(
-
)/2.61e + 00
4.64e + 01(+)/7.21e + 00
1.32e + 01(≈)/5.44e + 00
1.37e + 01(≈)/
1.15
e
+
00
1.43e + 01/1.62e + 00
F
18
F
mean SD
2.29e + 01(+)/5.49e + 00
2.31e + 01(+)/1.04e + 01
4.01e + 01(+)/5.84e + 00
5.48e + 01(+)/6.62e + 00
2.10e + 01(+)/5.05e + 00
1.85e + 01(≈)/5.64e + 00
1.74
e
+
0
1/
4.65
e
+
0
0
F
19
F
mean SD
4.81e − 01(≈)/1.88e − 01
5.44e − 01(≈)/1.90e − 01
6.22
e
-
02
(
-
)/
2.85
e
-
02
7.17e + 00(+)/1.21e + 00
5.33e − 01(≈)/1.58e − 01
4.61e − 01(≈)/1.28e − 01
4.68e − 01/1.69e − 01
F
20
F
mean SD
1.81
e
+
00
(
-
)/7.77e − 01
2.73e + 00(≈)/7.18e − 01
3.40e + 00(+)/
2.40
e
-
01
3.08e + 00(+)/2.19e − 01
2.10e + 00(−)/5.25e − 01
2.54e + 00(≈)/3.42e − 01
2.72e + 00/5.83e − 01
General merit over contender
3
3
7
19
3
3
“+,” “−,” and “≈” denote that the performance of CPEA (D
=
10) is better than, worse than, and similar to that of other algorithms.Table8 consists of unimodal problems (F
0
1
–
F
0
5) and basic multimodal problems (F
0
6
–
F
20). From the results of Table 8, we can see that PSO, ABC, PSO-FDR, PSO-cf-Local, and CPEA work well for F
0
1, where PSO-FDR, PSO-cf-Local, and CPEA can always achieve the optimal solution in each run. CPEA outperforms other algorithms on F
0
1, F
0
2, F
0
3, F
0
7, F
0
8, F
10, F
12, F
15, F
16, and F
18; the cause of the outstanding performance may be because of the phase transformation mechanism, which leads to faster convergence in the late stage of evolution.


## 3.1. Parameter Study

In this section, the setting basis for initial values of some parameters involved in the cloud particles evolution algorithm model is described. Firstly, we hope that all the cloud particles which are generated by the cloud generator can fall into the search space as much as possible. Secondly, the diversity of the population can be ensured effectively if the cloud particles are distributed in the search space uniformly. Therefore, some experiments are designed to analyze the feasibility of the selected scheme in order to illustrate the rationality of the selected parameters.


### 3.1.1. Definition ofE
n

Let the search rangeS be [
-
r
,
r
] and r is the search radius. Set r
=
2, E
x
1
=
-
2
/
3, E
x
2
=
2
/
3, H
e
=
0.001
E
n. Five hundred cloud particles are generated in every experiment. The experiment is independently simulated 30 times. The cloud particle is called a valid cloud particle if it falls into the search space. Figure 9 is a graph which shows the number of the valid cloud particles generated with different E
n, and Figure 10 is a graph which shows the distribution of valid cloud particles generated with different E
n.Figure 9
Valid cloud particles of differentE
n.Figure 10
Distribution of cloud particles in differentE
n.It can be seen from Figure9 that the number of the valid cloud particles reduces as the E
n increases. The average number of valid cloud particles is basically the same when E
n is (
1
/
3
)
r, (
1
/
2
)
r, and (
2
/
3
)
r, respectively. However, the number of valid cloud particles is relatively fewer when E
n is (
3
/
4
)
r, (
4
/
5
)
r, and r. It can be seen from Figure 10 that all the cloud particles fall into the research space when E
n is (
1
/
3
)
r, (
1
/
2
)
r, and (
2
/
3
)
r, respectively. However, the cloud particles concentrate in a certain area of the search space when E
n is (
1
/
3
)
r or (
1
/
2
)
r. By contrast, the cloud particles are widely distributed throughout all the search space [
-
2
,
2
] when E
n is (
2
/
3
)
r. Therefore, set E
n
=
(
2
/
3
)
r. The dashed rectangle in Figure 10 represents the decision variables space.


### 3.1.2. Definition ofH
e

Ackley Function (multimode) and Sphere Function (unimode) are selected to calculate the number of fitness evaluations (NFES) of the algorithm when the solution reachesf_stop with different H
e. Let the search space be [
-
32,32
] (Ackley) and [
-
100,100
] (Sphere), f
_
stop
=
10
-
3. H
e
1
=
E
n, H
e
2
=
0.1
E
n, H
e
3
=
0.01
E
n, H
e
4
=
0.001
E
n, H
e
5
=
0.0001
E
n, H
e
6
=
0.00001
E
n, and H
e
7
=
0.000001
E
n. The reported values are the average of the results for 50 independent runs. As shown in Figure 11, for Ackley Function, the number of fitness evaluations is least when H
e is equal to 0.001
E
n. For Sphere Function, the number of fitness evaluations is less when H
e is equal to 0.001
E
n. Therefore, H
e is equal to 0.001
E
n in the cloud particles evolution algorithm.NFES of CPEA in differentH
e.
(a)
(b)


### 3.1.3. Definition ofE
x

Ackley Function (multimode) and Sphere Function (unimode) are selected to calculate the number of fitness evaluations (NFES) of the algorithm when the solution reachesf_stop with E
x being equal to r, (
1
/
2
)
r, (
1
/
3
)
r, (
1
/
4
)
r, and (
1
/
5
)
r, respectively. Let the search space be [
-
r
,
r
], f
_
stop
=
10
-
3. As shown in Figure 12, for Ackley Function, the number of fitness evaluations is least when E
x is equal to (
1
/
2
)
r. For Sphere Function, the number of fitness evaluations is less when E
x is equal to (
1
/
2
)
r or (
1
/
3
)
r. Therefore, E
x is equal to (
1
/
2
)
r in the cloud particles evolution algorithm.NFES of CPEA in differentE
x.
(a)
(b)


### 3.1.4. Definition ofc
/
b

Many things in nature have the mathematical proportion relationship between their parts. If the substance is divided into two parts, the ratio of the whole and the larger part is equal to the ratio of the larger part to the smaller part; this ratio is defined as golden ratio, expressed algebraically(12)
a
+
b
a
=
a
b
≝
φ
,
in which φ represents the golden ration. Its value is
(13)
φ
=
1
+
5
2
=
1.6180339887
,
1
φ
=
φ
-
1
=
0.618
⋯
.The golden section has a wide application not only in the fields of painting, sculpture, music, architecture, and management but also in nature. Studies have shown that plant leaves, branches, or petals are distributed in accordance with the golden ratio in order to survive in the wind, frost, rain, and snow. It is the evolutionary result or the best solution which adapts its own growth after hundreds of millions of years of long-term evolutionary process.In this model, competition and reciprocity of the population depend on the fitness and are a dynamic relationship which constantly adjust and change. With competition and reciprocity mechanism, the convergence speed and the population diversity are improved. Premature convergence is avoided, and large search space and complex computational problems are solved. Indirect reciprocity is selected in this paper. Here, setc
/
b
=
(
1
-
0.618
)
≈
0
.
38.


## 3.1.1. Definition ofE
n

Let the search rangeS be [
-
r
,
r
] and r is the search radius. Set r
=
2, E
x
1
=
-
2
/
3, E
x
2
=
2
/
3, H
e
=
0.001
E
n. Five hundred cloud particles are generated in every experiment. The experiment is independently simulated 30 times. The cloud particle is called a valid cloud particle if it falls into the search space. Figure 9 is a graph which shows the number of the valid cloud particles generated with different E
n, and Figure 10 is a graph which shows the distribution of valid cloud particles generated with different E
n.Figure 9
Valid cloud particles of differentE
n.Figure 10
Distribution of cloud particles in differentE
n.It can be seen from Figure9 that the number of the valid cloud particles reduces as the E
n increases. The average number of valid cloud particles is basically the same when E
n is (
1
/
3
)
r, (
1
/
2
)
r, and (
2
/
3
)
r, respectively. However, the number of valid cloud particles is relatively fewer when E
n is (
3
/
4
)
r, (
4
/
5
)
r, and r. It can be seen from Figure 10 that all the cloud particles fall into the research space when E
n is (
1
/
3
)
r, (
1
/
2
)
r, and (
2
/
3
)
r, respectively. However, the cloud particles concentrate in a certain area of the search space when E
n is (
1
/
3
)
r or (
1
/
2
)
r. By contrast, the cloud particles are widely distributed throughout all the search space [
-
2
,
2
] when E
n is (
2
/
3
)
r. Therefore, set E
n
=
(
2
/
3
)
r. The dashed rectangle in Figure 10 represents the decision variables space.


## 3.1.2. Definition ofH
e

Ackley Function (multimode) and Sphere Function (unimode) are selected to calculate the number of fitness evaluations (NFES) of the algorithm when the solution reachesf_stop with different H
e. Let the search space be [
-
32,32
] (Ackley) and [
-
100,100
] (Sphere), f
_
stop
=
10
-
3. H
e
1
=
E
n, H
e
2
=
0.1
E
n, H
e
3
=
0.01
E
n, H
e
4
=
0.001
E
n, H
e
5
=
0.0001
E
n, H
e
6
=
0.00001
E
n, and H
e
7
=
0.000001
E
n. The reported values are the average of the results for 50 independent runs. As shown in Figure 11, for Ackley Function, the number of fitness evaluations is least when H
e is equal to 0.001
E
n. For Sphere Function, the number of fitness evaluations is less when H
e is equal to 0.001
E
n. Therefore, H
e is equal to 0.001
E
n in the cloud particles evolution algorithm.NFES of CPEA in differentH
e.
(a)
(b)


## 3.1.3. Definition ofE
x

Ackley Function (multimode) and Sphere Function (unimode) are selected to calculate the number of fitness evaluations (NFES) of the algorithm when the solution reachesf_stop with E
x being equal to r, (
1
/
2
)
r, (
1
/
3
)
r, (
1
/
4
)
r, and (
1
/
5
)
r, respectively. Let the search space be [
-
r
,
r
], f
_
stop
=
10
-
3. As shown in Figure 12, for Ackley Function, the number of fitness evaluations is least when E
x is equal to (
1
/
2
)
r. For Sphere Function, the number of fitness evaluations is less when E
x is equal to (
1
/
2
)
r or (
1
/
3
)
r. Therefore, E
x is equal to (
1
/
2
)
r in the cloud particles evolution algorithm.NFES of CPEA in differentE
x.
(a)
(b)


## 3.1.4. Definition ofc
/
b

Many things in nature have the mathematical proportion relationship between their parts. If the substance is divided into two parts, the ratio of the whole and the larger part is equal to the ratio of the larger part to the smaller part; this ratio is defined as golden ratio, expressed algebraically(12)
a
+
b
a
=
a
b
≝
φ
,
in which φ represents the golden ration. Its value is
(13)
φ
=
1
+
5
2
=
1.6180339887
,
1
φ
=
φ
-
1
=
0.618
⋯
.The golden section has a wide application not only in the fields of painting, sculpture, music, architecture, and management but also in nature. Studies have shown that plant leaves, branches, or petals are distributed in accordance with the golden ratio in order to survive in the wind, frost, rain, and snow. It is the evolutionary result or the best solution which adapts its own growth after hundreds of millions of years of long-term evolutionary process.In this model, competition and reciprocity of the population depend on the fitness and are a dynamic relationship which constantly adjust and change. With competition and reciprocity mechanism, the convergence speed and the population diversity are improved. Premature convergence is avoided, and large search space and complex computational problems are solved. Indirect reciprocity is selected in this paper. Here, setc
/
b
=
(
1
-
0.618
)
≈
0
.
38.


## 3.2. Experiments Settings

For fair comparison, we set the parameters to be fixed. For CPEA,population
size
=
5, c
h
i
l
d
n
u
m
=
30, E
n
=
(
2
/
3
)
r, (the search space is [
-
r
,
r
]), H
e
=
E
n
/
1000, E
x
=
(
1
/
2
)
r, cd
=
1000, δ
=
10
-
2, and sf
=
1000. For CMA-ES, the parameter values are chosen from [27]. The parameters of JADE come from [28]. For PSO, the parameter values are chosen from [29]. The parameters of DE are set to be F
=
0
.
3 and CR
=
0.9 used in [30]. For ABC, the number of colony size NP
=
20 and the number of food sources FoodNumber
=
NP
/
2 and limit
=
100. A food source will not be exploited anymore and is assumed to be abandoned when limit is exceeded for the source. For CEBA, population
size
=
10, c
h
i
l
d
n
u
m
=
50, E
n
=
1, and H
e
=
E
n
/
6. The parameters of PSO-cf-Local and FDR-PSO come from [31, 32].


## 3.3. Comparison Strategies

Five comparison strategies come from the literature [33] and literature [26] to evaluate the algorithms. The comparison strategies are described as follows.Successful Run. A run during which the algorithm achieves the fixed accuracy level within the max number of function evaluations for the particular dimension is shown as follows:
(14)
SR
=
successful
runs
total
runs
.Error. The error of a solution x is defined as f
(
x
)
-
f
(
x
*
), where x
* is the global minimum of the function. The minimum error is written when the max number of function evaluations is reached in 30 runs. Both the average and standard deviation of the error values are calculated.Convergence Graphs. The convergence graphs show the median performance of the total runs with termination by either the max number of function evaluations or the termination error value.Number of Function Evaluations (NFES). The number of function evaluations is recorded when the value-to-reach is reached. The average and standard deviation of the number of function evaluations are calculated.Acceleration Rate (AR). The acceleration rate is defined as follows, based on the NFES for the two algorithms:
(15)
AR
=
NFE
S
a
1
NFE
S
a
2
,
in which AR
>
1 means Algorithm 2 (a
2) is faster than Algorithm 1 (a
1).The average acceleration rate (A
R
ave) and the average success rate (S
R
ave) over n test functions are calculated as follows:
(16)
A
R
ave
=
1
n
∑
i
=
1
n
A
R
i
,
SR
ave
=
1
n
∑
i
=
1
n
SR
i
.


## 3.4. Comparison of CMA-ES, DE, RES, LOS, ABC, CEBA, and CPEA

This section compares CMA-ES, DE, RES, LOS, ABC, CEBA [34], and CPEA in terms of mean number of function evaluations for successful runs. Results of CMA-ES, DE, RES, and LOS are taken from literature [35]. Each run is stopped and regarded as successful if the function value is smaller than f_stop.As shown in Table1, CPEA can reach f_stop faster than other algorithms with respect to f
07 and f
17. However, CPEA converges slower than CMA-ES with respect to f
18. When the dimension is 30, CPEA can reach f_stop relatively fast with respect to f
33. Nonetheless, CPEA converges slower while the dimension increases.Table 1
Average number of function evaluations to reachf_stop of CPEA versus CMA-ES, DE, RES, LOS, ABC, and CEBA on f
07, f
17, f
18, and f
33 function.
F
f_stop
init
D
CMA-ES
DE
RES
LOS
ABC
CEBA
CPEA
f
07
1
e
-
3
[
-600, 600
]
D
20
3111
8691
—
—
12463
94412
2950
30
4455
11410
2
e
5
—
29357
163545
4270
100
12796
31796
—
—
34777
458102
7260
f
17
0.9
[
-5.12, 5.12]
D
20
68586
12971
—
9.2
e
4
11183
30436
2850
DE:[
-600, 600]
D
30
147416
20150
1.0
e
5
2.3
e
5
26990
44138
6187
100
1010989
73620
—
—
449404
60854
10360
f
18
1
e
-
3
[
-30, 30]
D
20
2667
—
—
6.0
e
4
10183
64400
3330
30
3701
12481
1.1
e
5
9.3
e
4
17810
91749
4100
100
11900
36801
—
—
59423
205470
12248
f
33
0.9
[
-5.12, 5.12]
D
30
152000
1.25
e
6
—
—
—
50578
13100
100
240899
5522
—
—
—
167324
14300
The bold entities indicate the best results obtained by the algorithm.D represents dimension.


## 3.5. Comparison of CMA-ES, DE, PSO, ABC, CEBA, and CPEA

This section compares the error of a solution of CMA-ES, DE, PSO, ABC [36], CEBA, and CPEA under a given NFES. The best solution of every single function in different algorithms is highlighted with boldface.As shown in Table2, CPEA obtains better results in function f
0
1, f
0
2, f
0
4, f
0
5, f
0
7, f
0
9, f
10, f
11, and f
12 while it obtains relatively bad results in other functions. The reason for this phenomenon is that the changing rule of E
n is set according to the priori information and background knowledge instead of adaptive mechanism in CPEA. Therefore, the results of some functions are relatively worse. Besides, E
n which can be optimized in order to get better results is an important controlling parameter in the algorithm.Table 2
Comparison on the ERROR values of CPEA versus CMA-ES, DE, PSO, ABC, and CEBA on functionf
01~f
13.
F
NFES
D
CMA-ESmean/SD
DEmean/SD
PSOmean/SD
ABCmean/SD
CEBAmean/SD
CPEAmean/SD
f
01
30 k
30
5.87e − 29/1.56e − 29
3.20e − 07/2.51e − 07
5.44e − 12/1.56e − 11
9.71e − 16/2.26e − 16
1.74e − 02/1.55e − 02
3.1
9
e
-
33/
6.2
1
e
-
33
30 k
50
1.02e − 28/1.87e − 29
2.95e − 03/2.88e − 03
3.93e − 03/6.62e − 02
2.35e − 08/2.47e − 08
4.07e − 02/2.11e − 02
1.09
e
-
48/
1.41
e
-
48
f
02
10 k
30
1.50e − 05/3.66e − 05
2.30e + 00/1.78e + 00
6.24e − 06/8.74e − 05
0.00
e
+
00/
0.00
e
+
00
6.16e + 00/3.22e + 00
0.00
e
+
00/
0.00
e
+
00
50 k
50
0.00
e
+
00/
0.00
e
+
00
1.65e + 00/2.05e + 00
3.87e + 002.65e + 00
0.00
e
+
00/
0.00
e
+
00
6.08e + 00/5.34e + 00
0.00
e
+
00/
0.00
e
+
00
f
03
10 k
2
2.65e − 30/8.08e − 31
0.00
e
+
00/
0.00
e
+
00
5.01e − 14/8.52e − 14
4.83e − 04/1.40e − 03
5.83e − 04/8.32e − 04
1.52e − 13/2.44e − 13
f
04
100 k
30
8.2
8
e
-
27/
1.6
9
e
-
27
2.47e − 03/2.50e − 03
1.69e − 02/2.57e − 02
1.36e + 00/1.93e + 00
1.59e + 00/7.93e + 00
3.31e − 14/1.42e − 13
100 k
50
2.29e − 06/2.22e − 06
7.37e − 01/1.93e − 01
5.01e + 00/3.88e + 00
1.19e + 02/4.54e + 01
3.13e + 02/3.94e + 02
3.98
e
-
07/
8.79
e
-
07
f
05
30 k
30
8.19e − 09/5.08e − 09
1.27e − 12/3.94e − 12
8.66e − 34/2.88e − 33
8.00e − 13/2.24e − 12
2.54e − 09/4.13e − 09
2.11
e
-
37/
4.20
e
-
36
50 k
50
2.68e − 08/1.67e − 09
1.93e − 13/2.01e − 13
4.37
e
-
54/
7.58
e
-
54
9.54e − 12/1.55e − 11
1.34e − 15/2.01e − 15
1.07e − 23/1.85e − 23
f
06
30 k
2
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
3.60e − 09/1.05e − 08
1.1e − 03/8.2e − 04
7.71e − 09/4.52e − 09
f
07
30 k
30
1.88e − 12/1.40e − 12
3.08e − 06/5.41e − 06
1.64e − 09/2.85e − 08
1.54e − 09/2.78e − 08
2.07e − 06/5.57e − 05
5.32
e
-
16/
7.59
e
-
16
50 k
50
0.00
e
+
00/
0.00
e
+
00
3.28e − 05/5.69e − 05
1.79e − 02/2.49e − 02
1.27e − 09/2.21e − 09
9.49e − 03/4.23e − 03
0.00
e
+
00/
0.00
e
+
00
f
08
30 k
3
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
2.80e − 03/2.13e − 03
3.85e − 06/6.21e − 06
f
09
30 k
2
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
0.00
e
+
00/
0.00
e
+
00
f
10
5 k
2
6.99e − 47/2.21e − 46
1.28e − 26/3.9e − 26
6.59e − 07/7.71e − 07
2.32e − 04/2.60e − 04
8.24e − 08/1.37e − 07
8.29
e
-
59/
3.42
e
-
58
f
11
30 k
30
1.23e − 04/1.82e − 04
3.74e − 01/2.91e − 01
4.78e − 01/1.94e − 01
2.39e + 00/3.27e + 00
1.42e + 00/8.13e + 00
5.70
e
-
40/
6.61
e
-
39
50 k
50
5.58e − 01/1.55e − 01
3.36e + 00/1.55e + 00
3.15e + 01/1.11e + 01
5.71e + 01/6.98e + 00
2.05e + 01/9.40e + 00
2.33
e
-
05/
3.58
e
-
05
f
12
30 k
30
3.19e − 14/1.45e − 14
1.20e − 08/2.24e − 08
2.15e − 05/5.18e − 04
3.56e − 10/1.22e − 10
1.93e − 07/6.45e − 08
2.33
e
-
23/
2.50
e
-
22
50 k
50
1.41e − 14/2.43e − 14
6.77e − 05/1.11e − 05
6.08e − 03/3.35e − 03
6.37e − 10/6.83e − 10
9.76e − 03/3.08e − 03
1.90
e
-
72/
3.29
e
-
72
f
13
30 k
30
3.94e − 11/1.15e − 10
1.34
e
-
31/3.21e − 31
1.34
e
-
31/
0.00
e
+
00
6.57e − 04/3.77e − 03
2.09e − 02/7.34e − 02
1.20e − 12/2.82e − 12
50 k
50
3.08e − 03/3.80e − 03
1.34
e
-
31/
0.00
e
+
00
1.34
e
-
31/
0.00
e
+
00
1.41e − 02/9.29e − 03
4.22e − 02/2.53e − 02
2.36e − 12/1.04e − 12
“SD” stands for standard deviation.


## 3.6. Comparison of JADE, DE, PSO, ABC, CEBA, and CPEA

This section compares JADE, DE, PSO, ABC, CEBA, and CPEA in terms of SR and NFES for most functions atD
=
30 except f
03 (D
=
2). The termination criterion is to find a value smaller than 10
-
8. The best result of different algorithms for each function is highlighted in boldface.As shown in Table3, CPEA can reach the specified precision successfully except for function f
14 which successful run is 96. In addition, CPEA can reach the specified precision faster than other algorithms in function f
0
1, f
0
2, f
0
4, f
11, f
14, and f
17. For other functions, CPEA spends relatively long time to reach the precision. The reason for this phenomenon is that, for some functions, the changing rule of E
n given in the algorithm makes CPEA to spend more time to get rid of the local optima and converge to the specified precision.Table 3
Comparison on the SR and the NFES of CPEA versus JADE, DE, PSO, ABC, and CEBA for each function.
F
f
01  mean/SD
f
02  mean/SD
f
03  mean/SD
f
04  mean/SD
f
07  mean/SD
f
11  mean/SD
f
14  mean/SD
f
16  mean/SD
f
17  mean/SD
JADE
SR
1
1
1
1
1
0.98
0.94
1
1
NFES
2.8e + 4/6.0e + 2
1.0e + 4/3.5e + 2
5.7e + 3/6.9e + 2
6.9e + 4/3.0e + 3
3
.
3
e
+
3/
2.6
e
+
2
6.6e + 4/4.2e + 3
4.0e + 4/2.1e + 3
1
.
3
e
+
3/
3.8
e
+
2
1.3e + 5/1.8e + 3
DE
SR
1
1
1
1
0.96
1
0.94
1
0
NFES
3.3e + 4/6.6e + 2
1.1e + 4/2.3e + 3
1
.
6
e
+
3/
1.1
e
+
2
1.9e + 5/1.2e + 4
3.4e + 4/1.6e + 3
1.4e + 5/1.0e + 4
2.6e + 4/1.6e + 4
2.1e + 4/8.1e + 3
—
PSO
SR
1
1
1
0
0.86
0
0
1
0
NFES
2.9e + 4/7.5e + 2
1.6e + 4/1.4e + 3
7.7e + 3/9.6e + 2
—
2.8e + 4/9.8e + 2
—
—
1.8e + 4/7.2e + 2
—
ABC
SR
1
1
1
0
0.96
0
0
1
1
NFES
1.7e + 4/1.0e + 3
6.3e + 3/1.5e + 3
1.6e + 4/2.6e + 3
—
1.9e + 4/2.3e + 3
—
—
6.5e + 3/3.9e + 2
3.9e + 4/1.5e + 4
CEBA
SR
1
1
1
0
1
1
0.96
1
1
NFES
1.3e + 5/5.8e + 3
5.9e + 4/4.6e + 3
4.2e + 4/6.1e + 3
—
2.2e + 5/1.2e + 4
1.3e + 5/5.8e + 3
2.9e + 5/1.4e + 4
5.4e + 4/3.8e + 3
1.1e + 5/7.3e + 3
CPEA
SR
1
1
1
1
1
1
0.96
1
1
NFES
4.3
e
+
3/
5
.
1
e
+
2
4.5
e
+
3/
2
.
7
e
+
2
5.7e + 3/3.0e + 3
6.1
e
+
4/
2.6
e
+
3
7.7e + 3/3.4e + 2
2.3
e
+
4/
3.4
e
+
3
3.4
e
+
4/
1.6
e
+
3
2.4e + 3/4.3e + 2
1
.
0
e
+
4/
1
.
3
e
+
3As shown in Table4, the solution obtained by CPEA is closer to the best solution compared with other algorithms except for functions f
13 and f
20. The reason is that the phase transformation mechanism introduced in the algorithm can improve the convergence of the population greatly. Besides, the reciprocity mechanism of the algorithm can improve the adaptation and the diversity of the population in case the algorithm falls into the local optimum. Therefore, the algorithm can obtain better results in these functions.Table 4
Comparison on the ERROR values of CPEA versus JADE, DE, PSO, and ABC atD
=
50 except f
21 (D
=
2) and f
22 (D
=
2) for each function.
F
NFES
JADEmean/SD
DEmean/SD
PSOmean/SD
ABCmean/SD
CPEAmean/SD
f
13
40 k
1.34
e
-
31
/
0.00
e
+
00
1.34
e
-
31
/
0.00
e
+
00
1.35e − 31/0.00
e
+
00
3.52e − 03/6.33e − 02
3.64e − 12/4.68e − 12
f
15
40 k
9.18e − 11/7.04e − 11
2.94e − 05/3.66e − 05
1.81e − 09/4.18e − 08
4.11e − 14/2.62e − 14
4.5
9
e
-
53
/
1.36
e
-
52
f
16
40 k
1.15e − 20/1.42e − 20
6.62e − 08/1.24e − 07
6.91e − 19/7.55e − 19
3.41e − 15/5.23e − 15
1.77
e
-
76
/
0.00
e
+
00
f
17
200 k
6.71e − 07/1.32e − 07
2.58e − 01/6.15e − 01
6.73e + 00/2.06e + 00
6.07e − 09/1.85e − 08
0.00
e
+
00
/
0.00
e
+
00
f
18
50 k
1.16e − 07/4.27e − 08
3.89e − 04/2.10e − 04
1.70e − 09/3.13e − 08
2.93e − 09/2.34e − 08
8.88
e
-
16
/
0.00
e
+
00
f
19
40 k
1.14e − 02/1.68e − 02
9.90e − 03/6.66e − 03
9.00e − 06/1.62e − 05
2.22e − 05/1.95e − 04
1.64
e
-
64
/
5.13
e
-
64
f
20
40 k
0.00
e
+
00
/
0.00
e
+
00
0.00
e
+
00
/
0.00
e
+
00
1.40e − 10/2.81e − 10
4.38e − 03/3.33e − 03
1.21e − 10/3.11e − 09
f
21
30 k
0.00
e
+
00
/
0.00
e
+
00
0.00
e
+
00
/
0.00
e
+
00
0.00
e
+
00
/
0.00
e
+
00
1.06e − 17/3.34e − 17
0.00
e
+
00
/
0.00
e
+
00
f
22
30 k
1.28e − 45/1.86e − 45
4.99e − 164/0.00
e
+
00
2.63e − 51/6.82e − 51
6.62e − 18/5.84e − 18
0.00
e
+
00
/
0.00
e
+
00
f
23
30 k
0.00
e
+
00
/
0.00
e
+
00
0.00
e
+
00
/
0.00
e
+
00
0.00
e
+
00
/
0.00
e
+
00
0.00
e
+
00/4.24e − 04
0.00
e
+
00
/
0.00
e
+
00As shown in Table5, CPEA can find the optimal solution for the test functions. However, the stability of CPEA is worse than JADE and DE (f
25, f
26, f
27, f
28 and f
31). This result shows that the setting of H
e needs to be further improved.Table 5
Comparison on the optimum solution of CPEA versus JADE, DE, PSO, and ABC for each function.
F
NFES
JADE mean/SD
DE mean/SD
PSO mean/SD
ABC mean/SD
CPEA mean/SD
f
08
40 k
-
3.86278
/
0.00
e
+
00
-
3.86278
/
0.00
e
+
00
-
3.86278
/
6.20
e
-
17
-
3.86278
/
0.00
e
+
00
-
3.86278
/
0.00
e
+
00
f
24
40 k
-
3.3
1044
/
3.60
e
-
02
-
3.24608
/
5.70
e
-
02
-
3.25800
/
5.90
e
-
02
-
3.33539
/
0.00
e
+
00
-
3.3353
9
/
3.03
e
-
06
f
25
40 k
-
10.1532
/
4.00
e
-
14
-
10.1532
/
4.20
e
-
16
-
5.79196
/
3.30
e
+
00
-
10
.
153
2
/
1.25
e
-
13
-
10
.
153
2
/
1.58
e
-
09
f
26
40 k
-
10.4029
/
9.40
e
-
16
-
10.4029
/
2.50
e
-
16
-
7.14128
/
3.50
e
+
00
-
10
.
4029
/
1.77
e
-
15
-
10
.
402
9
/
1.94
e
-
10
f
27
40 k
-
10.5364
/
8.10
e
-
12
-
10.5364
/
7.10
e
-
16
-
7.02213
/
3.60
e
+
00
-
10
.
536
4
/
4.03
e
-
05
-
10
.
536
4
/
2.02
e
-
09
f
28
40 k
3.00000
/
1.10
e
-
15
3.00000
/
6.4
e
-
16
3.00000
/
1.30
e
-
15
3.0000
0
/
6.17
e
-
05
3.00000
/
2
.
01
e
-
15
f
29
40 k
2.55
e
-
20
/
4.79
e
-
20
0.00
e
+
00
/
0.00
e
+
00
2.10
e
-
03
/
4.28
e
-
03
1.14
e
-
18
/
1.62
e
-
18
0.00
e
+
00
/
0.00
e
+
00
f
30
40 k
-
2.3458
/
0.00
e
+
00
-
2.3458
/
0.00
e
+
00
-
2.3458
/
0.00
e
+
00
-
2.3458
/
0.00
e
+
00
-
2.3458
/
0.00
e
+
00
f
31
40 k
4.78
e
-
31
/
1.07
e
-
30
3.59
e
-
55
/
1.13
e
-
54
1.22
e
-
36
/
3.85
e
-
36
1.40
e
-
11
/
1.38
e
-
11
8.69
e
-
25
/
1.14
e
-
24
f
32
40 k
1
.
49
e
-
02
/
4
.
23
e
-
03
3.46
e
-
02
/
5.14
e
-
03
6.58
e
-
02
/
1.22
e
-
02
2.95
e
-
01
/
1.13
e
-
01
5
.
70
e
-
03
/
3
.
58
e
-
03
Data of functionsf
08, f
24, f
25, f
26, f
27, and f
28 in JADE, DE, and PSO are from literature [28].


## 3.7. Comparison of CMAES, JADE, DE, PSO, ABC, CEBA, and CPEA

The convergence curves of different algorithms on some selected functions are plotted in Figures13 and 14. The solutions obtained by CPEA are better than other algorithms with respect to functions f
0
1, f
11, f
16, f
32, and f
33. However, the solution obtained by CPEA is worse than CMA-ES with respect to function f
0
4. The solution obtained by CPEA is worse than JADE for high dimensional problems such as function f
0
5 and function f
0
7. The reason is that the setting of initial population of CPEA is insufficient for high dimensional problems.Figure 13
Convergence graph for test functionsf
0
1, f
0
4, f
16, and f
32 (D
=
30).Figure 14
Convergence graph for test functionsf
0
5, f
11, f
33 (D
=
50), and f
0
7 (D
=
100).


## 3.8. Comparison of DE, ODE, RDE, PSO, ABC, CEBA, and CPEA

The algorithms DE, ODE, RDE, PSO, ABC, CBEA, and CPEA are compared in terms of convergence speed. The termination criterion is to find error values smaller than10
-
8. The results of solving the benchmark functions are given in Tables 6 and 7. Data of all functions in DE, ODE, and RDE are from literature [26]. The best result of NEFS and the SR for each function are highlighted in boldface. The average success rates and the average acceleration rate on test functions are shown in the last row of the tables.Table 6
Comparison of DE, ODE, RDE, PSO, ABC, CEBA, and CPEA. The best result for each case is highlighted in boldface.
F
D
DE
ODE
ARODE
RDE
ARRDE
PSO
ARPSO
ABC
ARABC
CEBA
ARCEBA
CPEA
ARCPEA
NFES
SR
NFES
SR
NFES
SR
NFES
SR
NFES
SR
NFES
SR
NFES
SR
f
0
1
30
87748
1
47716
1
1.83
115096
1
0.76
20290
1
4.32
18370
1
4.77
133381
1
0.66
5266
1
16.66
f
0
2
30
41588
1
23124
1
1.80
54316
1
0.77
11541
1
3.60
6247
1
6.66
58416
1
0.71
2966
1
14.02
f
0
3
2
4324
1
4776
1
0.90
4904
1
0.88
6110
1
0.71
18050
1
0.24
38278
1
0.11
4350
1
0.99
f
0
4
20
177880
1
168680
1
1.05
231152
1
0.77
—
0
—
—
0
—
—
0
—
11068
1
16.07
f
0
5
30
25140
1
8328
1
3.01
29096
1
0.86
5910
1
4.25
6134
1
4.10
18037
1
1.39
1700
1
14.79
f
0
6
2
8016
1
6608
1
1.21
8248
1
0.97
7185
1
1.12
25407
1
0.32
66319
1
0.12
5033
1
1.59
f
0
7
30
113428
1
69342
0.96
1.64
149744
1
0.76
20700
0.84
5.48
20930
0.96
5.42
209278
1
0.54
9941
1
11.41
f
0
8
3
3376
1
3580
1
0.94
3652
1
0.92
4090
1
0.83
1492
1
2.26
206859
1
0.02
31874
1
0.11
f
0
9
2
5352
1
4468
1
1.20
6292
1
0.85
2830
1
1.89
1077
1
4.97
9200
1
0.58
3856
1
1.39
f
10
2
3608
1
3288
1
1.09
3924
1
0.92
5485
1
0.66
10891
1
0.33
6633
1
0.54
1920
1
1.88
f
11
30
385192
1
369104
1
1.04
498200
1
0.77
—
0
—
—
0
—
142564
1
2.70
12662
1
30.42
f
12
30
187300
1
155636
1
1.20
244396
1
0.76
—
0
—
25787
1
7.23
157139
1
1.19
10183
1
18.39
f
13
30
101460
1
70408
1
1.44
138308
1
0.73
6512
1
15.58
—
0
—
—
0
—
4300
1
23.60
f
14
30
570290
0.28
72250
0.88
7.89
285108
1
2.00
—
0
—
—
0
—
292946
1
1.95
31112
0.98
18.33
f
15
30
96488
1
53304
1
1.81
126780
1
0.76
18757
1
5.14
15150
1
6.37
103475
1
0.93
6466
1
14.92
f
16
2
1016
1
996
1
1.02
1084
1
0.94
428
1
2.37
301
1
3.38
3292
1
0.31
610
1
1.67
f
17
10
328844
1
70389
0.76
4.67
501875
0.96
0.66
—
0
—
7935
1
41.44
48221
1
6.82
3725
1
88.28
f
18
30
169152
1
98296
1
1.72
222784
1
0.76
27250
0.92
6.21
32398
1
5.22
184041
1
0.92
13400
1
12.62
f
19
30
411164
1
337532
1
1.22
927230
0.24
0.44
—
0
—
156871
1
2.62
148014
1
2.78
6233
1
65.97
Ave
0.96
0.98
0.98
0.96
0.86
0.67
4.01
0.79
6.36
0.89
1.31
0.99
18.58Table 7
Comparison of DE, ODE, RDE, PSO, ABC, CEBA, and CPEA. The best result for each case is highlighted in boldface.
F
D
DE
ODE
ARODE
RDE
ARRDE
PSO
ARPSO
ABC
ARABC
CEBA
ARCEBA
CPEA
ARCPEA
NFES
SR
NFES
SR
NFES
SR
NFES
SR
NFES
SR
NFES
SR
NFES
SR
f
20
5
2163
0.88
2024
1
1.07
2904
0.96
0.74
11960
1
0.18
—
0
—
58252
1
0.04
4216
0.92
0.51
f
21
2
7976
1
5092
1
1.56
8332
1
0.96
8372
0.98
0.95
15944
1
0.50
11129
1
0.72
1025
1
7.78
f
22
2
3780
1
3396
1
1.11
3820
1
0.98
5755
1
0.66
2050
1
1.84
6635
1
0.57
855
1
4.42
f
23
10
19528
1
15704
1
1.24
23156
1
0.84
9290
1
2.10
5845
1
3.34
17907
1
1.09
1130
1
17.28
f
24
6
—
0
—
0
—
—
0
—
7966
0.94
—
3018
1
—
—
0
—
42716
0.98
—
f
25
4
—
0
—
0
—
—
0
—
—
0
—
4280
1
—
—
0
—
5471
0.94
—
f
26
4
—
0
—
0
—
—
0
—
—
0
—
8383
1
—
—
0
—
5391
0.96
—
f
27
4
4576
1
4500
1
1.02
5948
1
0.77
—
0
—
4850
1
0.94
—
0
—
5512
0.92
0.83
f
28
2
4780
1
4684
1
1.02
4924
1
0.97
6770
1
0.71
6228
1
0.77
14248
1
0.34
4495
1
1.06
f
29
2
8412
1
5772
1
1.46
10860
1
0.77
—
0
—
1914
1
4.40
37765
1
0.22
4950
1
1.70
f
30
2
3256
1
3120
1
1.04
3308
1
0.98
3275
1
0.99
1428
1
2.28
114261
1
0.03
3025
1
1.08
f
31
4
1108
1
1232
1
0.90
1388
1
0.80
990
1
1.12
9944
1
0.11
131217
1
0.01
2662
1
0.42
f
34
2
5232
1
5956
0.92
0.88
5768
1
0.91
3495
1
1.50
1852
1
2.83
9926
1
0.53
4370
1
1.20
f
35
5
38532
1
16340
1
2.36
46080
1
0.84
11895
1
3.24
25303
1
1.52
35292
1
1.09
2095
1
18.39
f
36
2
7888
1
4272
1
1.85
9148
1
0.86
10910
1
0.72
8359
1
0.94
13529
1
0.58
1340
1
5.89
f
37
2
5284
1
4728
1
1.12
5488
1
0.97
7685
1
0.69
1038
1
5.09
9056
1
0.58
1022
1
5.17
f
38
2
5280
1
4804
1
1.10
5540
1
0.95
7680
1
0.69
1598
1
3.30
8905
1
0.59
990
1
5.33
f
39
10
37824
1
24206
1
1.56
46800
1
0.80
—
0
—
—
0
—
153161
1
0.25
38420
0.94
0.98
f
40
2
2976
1
2872
1
1.03
3200
1
0.93
2790
1
1.07
666
1
4.47
9032
1
0.33
3085
1
0.96
Ave
0.84
0.84
1.27
0.84
0.88
0.73
1.12
0.89
2.31
0.77
0.46
0.98
4.56It can be seen from Tables6 and 7 that CPEA outperforms DE, ODE, RDE, PSO, ABC, and CBEA on 38 of the functions. CPEA performs better than other algorithms on 23 functions. In Table 6, the average AR is 0.98, 0.86, 4.01, 6.36, 1.31, and 18.58 for ODE, RDE, PSO, ABC, CBEA, and CPEA, respectively. In Table 7, the average AR is 1.27, 0.88, 1.12, 2.31, 0.46, and 4.56 for ODE, RDE, PSO, ABC, CBEA, and CPEA, respectively. In Table 6, the average SR is 0.96, 0.98, 0.96, 0.67, 0.79, 0.89, and 0.99 for DE, ODE, RDE, PSO, ABC, CBEA, and CPEA, respectively. In Table 7, the average SR is 0.84, 0.84, 0.84, 0.73, 0.89, 0.77, and 0.98 for DE, ODE, RDE, PSO, ABC, CBEA, and CPEA, respectively. These results clearly demonstrate that CPEA has better performance.


## 3.9. Test on the CEC2013 Benchmark Problems

To further verify the performance of CPEA, a set of recently proposed CEC2013 rotated benchmark problems are used. In the experiment,t-test [37] has been carried out to show the differences between CPEA and the other algorithm. “General merit over contender” displays the difference between the number of better results and the number of worse results, which is used to give an overall comparison between the two algorithms. Table 8 gives the information about the average error, standard deviation, and t-value of 30 runs of 7 algorithms over 300,000 FEs on 20 test functions with D
=
10. The best results among the algorithms are shown in bold.Table 8
Results for CEC2013 benchmark functions over 30 independent runs.
Function
Result
DE
PSO
ABC
CEBA
PSO-FDR
PSO-cf-Local
CPEA
F
0
1
F
mean SD
1.64e − 02(+)/3.98e − 02
1.51e − 14(≈)/5.76e − 14
1.51e − 14(≈)/5.76e − 14
4.99e + 02(+)/1.46e + 02
0
.
00
e
+
00
(
≈
)/
0
.
00
e
+
00
0
.
00
e
+
00
(
≈
)/
0
.
00
e
+
00
0.00e+00/
0
.
00
e
+
00
F
0
2
F
mean SD
1.20e + 05(+)/1.50e + 05
1.73e + 05(+)/2.02e + 05
3.31e + 06(+)/1.13e + 06
3.95e + 06(+)/1.60e + 06
4.37e + 04(+)/3.76e + 04
2.11e + 04(≈)/1.61e + 04
1.82
e
+
0
4/
1.12
e
+
0
4
F
0
3
F
mean SD
1.75e + 06(+)/6.09e + 05
4.63e + 05(+)/8.77e + 05
1.35e + 07(+)/1.30e + 07
6.79e + 08(+)/3.27e + 08
1.38e + 06(+)/1.89e + 06
2.92e + 04(+)/7.10e + 04
4.91
e
+
00/
2.52
e
+
0
1
F
0
4
F
mean SD
1.30e + 03(+)/1.43e + 03
1.97e + 02(−)/
1.11
e
+
0
2
1.16e + 04(+)/3.19e + 03
6.25e + 03(+)/2.25e + 03
1.54
e
+
0
2
(
-
)/1.60e + 02
8.18e + 02(≈)/5.26e + 02
7.26e + 02/3.26e + 02
F
0
5
F
mean SD
1.38e − 01(+)/3.62e − 01
5.30e − 14(−)/5.76e − 14
1.62e − 13(−)/5.72e − 14
1.15e + 02(+)/9.45e + 01
7.57e − 15(−)/2.88e − 14
0
.
00
e
+
00
(
-
)/
0
.
00
e
+
00
7.72e − 05/2.23e − 05
F
0
6
F
mean SD
6.30e + 00(≈)/4.47e + 00
6.21e + 00(≈)/4.36e + 00
1.51
e
+
00
(
-
)/
2.90
e
+
00
5.72e + 01(+)/1.94e + 01
5.28e + 00(≈)/4.92e + 00
9.26e + 00(+)/2.51e + 00
6.21e + 00/4.80e + 00
F
0
7
F
mean SD
7.25e − 01(+)/1.93e + 00
1.63e + 00(+)/1.89e + 00
4.08e + 01(+)/1.18e + 01
4.01e + 01(+)/5.31e + 00
1.89e + 00(+)/1.94e + 00
6.87e − 01(+)/1.23e + 00
1.69
e
-
03/
3.05
e
-
03
F
0
8
F
mean SD
2.02
e
+
01
(
≈
)/6.33e − 02
2.02
e
+
01
(
≈
)/6.88e − 02
2.04e + 01(+)/8.13e − 02
2.02
e
+
0
1
(
≈
)/5.86e − 02
2.02
e
+
0
1
(
≈
)/6.46e − 02
2.02
e
+
0
1
(
≈
)/7.73e − 02
2.02
e
+
0
1/
4.62
e
-
02
F
0
9
F
mean SD
1.03
e
+
00
(
-
)/
8.41
e
-
01
2.97e + 00(≈)/1.26e + 00
5.30e + 00(+)/8.55e − 01
6.83e + 00(+)/5.71e − 01
2.44e + 00(≈)/1.37e + 00
2.26e + 00(≈)/8.62e − 01
2.32e + 00/1.29e + 00
F
10
F
mean SD
5.64e − 02(≈)/4.11e − 02
3.88e − 01(+)/2.90e − 01
1.42e + 00(+)/3.27e − 01
6.24e + 01(+)/1.40e + 01
3.39e − 01(+)/2.25e − 01
1.12e − 01(+)/5.53e − 02
5
.
3
5
e
-
02/
2
.
81
e
-
02
F
1
1
F
mean SD
7.29e − 01(−)/8.63e − 01
5.63e − 01(−)/7.24e − 01
0
.
00
e
+
00
(
-
)/
0
.
00
e
+
00
4.56e + 01(+)/5.89e + 00
6.63e − 01(−)/6.57e − 01
1.22e + 00(−)/1.29e + 00
7.95e + 00/3.30e + 00
F
12
F
mean SD
5.31e + 00(≈)/2.27e + 00
1.38e + 01(+)/6.58e + 00
2.80e + 01(+)/8.22e + 00
4.59e + 01(+)/5.13e + 00
1.38e + 01(+)/4.80e + 00
6.36e + 00(+)/2.19e + 00
5.00
e
+
00/
1.21
e
+
00
F
13
F
mean SD
7.74
e
+
00
(
-
)/4.59e + 00
1.89e + 01(≈)/6.14e + 00
3.39e + 01(+)/9.86e + 00
4.68e + 01(+)/8.05e + 00
1.53e + 01(≈)/6.07e + 00
9.04e + 00(−)/
4.00
e
+
00
1.70e + 01/8.72e + 00
F
14
F
mean SD
2.83e + 00(−)/3.36e + 00
4.82e + 01(−)/5.11e + 01
1.46
e
-
01
(
-
)/
6.09
e
-
02
9.74e + 02(+)/1.12e + 02
5.65e + 01(−)/8.34e + 01
2.00e + 02(−)/1.00e + 02
3.17e + 02/1.53e + 02
F
15
F
mean SD
3.04e + 02(≈)/2.03e + 02
6.88e + 02(+)/2.91e + 02
7.35e + 02(+)/1.94e + 02
7.21e + 02(+)/1.65e + 02
6.38e + 02(+)/1.92e + 02
3.90e + 02(+)/1.24e + 02
3.02
e
+
0
2/
1.15
e
+
0
2
F
16
F
mean SD
9.56e − 01(+)1.42e − 01
7.13e − 01(+)/1.86e − 01
8.93e − 01(+)/1.91e − 01
1.02e + 00(+)/1.14e − 01
7.60e − 01(+)/2.38e − 01
5.88e − 01(+)/2.59e − 01
4.47
e
-
02/
4.10
e
-
02
F
17
F
mean SD
1.45e + 01(≈)/4.11e + 00
8.74e + 00(−)/4.71e + 00
8.55
e
+
00
(
-
)/2.61e + 00
4.64e + 01(+)/7.21e + 00
1.32e + 01(≈)/5.44e + 00
1.37e + 01(≈)/
1.15
e
+
00
1.43e + 01/1.62e + 00
F
18
F
mean SD
2.29e + 01(+)/5.49e + 00
2.31e + 01(+)/1.04e + 01
4.01e + 01(+)/5.84e + 00
5.48e + 01(+)/6.62e + 00
2.10e + 01(+)/5.05e + 00
1.85e + 01(≈)/5.64e + 00
1.74
e
+
0
1/
4.65
e
+
0
0
F
19
F
mean SD
4.81e − 01(≈)/1.88e − 01
5.44e − 01(≈)/1.90e − 01
6.22
e
-
02
(
-
)/
2.85
e
-
02
7.17e + 00(+)/1.21e + 00
5.33e − 01(≈)/1.58e − 01
4.61e − 01(≈)/1.28e − 01
4.68e − 01/1.69e − 01
F
20
F
mean SD
1.81
e
+
00
(
-
)/7.77e − 01
2.73e + 00(≈)/7.18e − 01
3.40e + 00(+)/
2.40
e
-
01
3.08e + 00(+)/2.19e − 01
2.10e + 00(−)/5.25e − 01
2.54e + 00(≈)/3.42e − 01
2.72e + 00/5.83e − 01
General merit over contender
3
3
7
19
3
3
“+,” “−,” and “≈” denote that the performance of CPEA (D
=
10) is better than, worse than, and similar to that of other algorithms.Table8 consists of unimodal problems (F
0
1
–
F
0
5) and basic multimodal problems (F
0
6
–
F
20). From the results of Table 8, we can see that PSO, ABC, PSO-FDR, PSO-cf-Local, and CPEA work well for F
0
1, where PSO-FDR, PSO-cf-Local, and CPEA can always achieve the optimal solution in each run. CPEA outperforms other algorithms on F
0
1, F
0
2, F
0
3, F
0
7, F
0
8, F
10, F
12, F
15, F
16, and F
18; the cause of the outstanding performance may be because of the phase transformation mechanism, which leads to faster convergence in the late stage of evolution.


## 4. Conclusions

This paper presents a new optimization method called cloud particles evolution algorithm (CPEA). The fundamental concepts and ideas come from the formation of the cloud and the phase transformation of the substances in nature. In this paper, CPEA solves 40 benchmark optimization problems which have difficulties in discontinuity, nonconvexity, multimode, and high-dimensionality. The evolutionary process of the cloud particles is divided into three states: the cloud gaseous state, the cloud liquid state, and the cloud solid state. The algorithm realizes the change of the state of the cloud with the phase transformation driving force according to the evolutionary process. The algorithm reduces the population number appropriately and increases the individuals of subpopulations in each phase transformation. Therefore, CPEA can improve the whole exploitation ability of the population greatly and enhance the convergence. There is a larger adjustment forE
n to improve the exploration ability of the population in melting or gasification operation. When the algorithm falls into the local optimum, it improves the diversity of the population through the indirect reciprocity mechanism to solve this problem.As we can see from experiments, it is helpful to improve the performance of the algorithm with appropriate parameter settings and change strategy in the evolutionary process. As we all know, the benefit is proportional to the risk for solving the practical engineering application problems. Generally, the better the benefit is, the higher the risk is. We always hope that the benefit is increased and the risk is reduced or controlled at the same time. In general, the adaptive system adapts to the dynamic characteristics change of the object and disturbance by correcting its own characteristics. Therefore, the benefit is increased. However, the adaptive mechanism is achieved at the expense of computational efficiency. Consequently, the settings of parameter should be based on the actual situation. The settings of parameters should be based on the prior information and background knowledge in some certain situations, or based on adaptive mechanism for solving different problems in other situations. The adaptive mechanism will not be selected if the static settings can ensure to generate more benefit and control the risk. Inversely, if the adaptive mechanism can only ensure to generate more benefit while it cannot control the risk, we will not choose it. The adaptive mechanism should be selected only if it can not only increase the benefit but also reduce the risk. In CPEA,E
n and H
e are statically set according to the priori information and background knowledge. As we can see from experiments, the results from some functions are better while the results from the others are relatively worse among the 40 tested functions.As previously mentioned, parameter setting is an important focus. Then, the next step is the adaptive settings ofE
n and H
e. Through adaptive behavior, these parameters can automatically modify the calculation direction according to the current situation while the calculation is running. Moreover, possible directions for future work include improving the computational efficiency and the quality of solutions of CPEA. Finally, CPEA will be applied to solve large scale optimization problems and engineering design problems.



---

*Source: 434831-2015-02-25.xml*