=Paper= {{Paper |id=Vol-3304/paper28 |storemode=property |title=An Adaptive Differential Evolution Based on Multi-subpopulation for Global Optimization |pdfUrl=https://ceur-ws.org/Vol-3304/paper28.pdf |volume=Vol-3304 |authors=Zhoucheng Bao,Kaige Chen,Tingting Pang,Zuling Wang }} ==An Adaptive Differential Evolution Based on Multi-subpopulation for Global Optimization== https://ceur-ws.org/Vol-3304/paper28.pdf
An Adaptive Differential Evolution Based on Multi-
subpopulation for Global Optimization
Zhoucheng Bao1, Kaige Chen1, Tingting Pang1, Zuling Wang2
1
    HangZhou Normal University, Hangzhou, China
2
    Chongqing Vocational Institute of Safety and Technology, Chongqing, China

                 Abstract
                 This paper proposed a multiple-subpopulation based on differential evolution with an adaptive
                 parameter control and population reduction, called EMPADE, for global optimization. In
                 EMPADE, an adaptive parameter control method based on population diversity and fitness
                 value is proposed to set the parameters F. Further, a mutation strategy named DE/current-
                 pbad/1 with archives is designed. In addition, a population size reduction strategy is devised to
                 adaptively eliminate poor performance Individuals. The proposed algorithm has been tested by
                 solving CEC2014 benchmark problems. The experimental results show that the EMPADE
                 algorithm strategy is feasible and effective, and has greater competitiveness compared with
                 other related algorithms.

                 Keywords 1
                 Adaptive, Differential evolution, Population size control, Multiple-population

1. Introduction

    Differential evolution (DE), initially introduced by Storn and Price [1], has become a powerful tool
to solve global optimization problems [2]. Applications of DE for real problems includes pattern
recognition [3], power systems [4] and vehicle routing [5].
    The selection of mutation operator is concerned by many researchers. For example, in JADE [6], a
parameter adaptation strategy along with DE/current-to-pbest mutation strategy is employed to generate
offspring. CoDE CoDE adopts multiple compound mutation operators (DE/rand/1, DE/rand/2 and
DE/current-to-rand/1) and DE variants of parameter control mode [8]. In the above methods, mutation
strategies are performed on a single population, rather than multiple populations.
    Recently, multi-population based DE algorithms have also been proposed. These algorithms generally
divide initial population into multiple sub-populations and share information among different sub-
populations through individual migration. For example, in MPEDE [9] and EDEV[10], the authors tried
to divide the population into a large reward sub-population and three small indicator sub-populations. In
these methods, three different mutation strategies are employed in an ensemble manner. In SAMDE[11],
the entire population is divided into three subpopulations of the same size, each employing a different
mutation strategy. After each generation, the population is randomly reorganized. However, random
recombination of the population could lead to lead to unreasonable allocation of computing resources.
Further, parameter control in these methods does not consider the diversity and performance of sub-
populations, restricting their performance [6].
    In this paper, we propose multiple-population based differential evolution with adaptive parameter
control and population reduction strategies, called EMPADE, for global optimization. In each generation,
three subpopulations with the same size will be divided according to the Euclidean distance, and the


ICBASE2022@3rd International Conference on Big Data & Artificial Intelligence & Software Engineering, October 21-
23, 2022, Guangzhou, China
EMAIL: baozc1031@163.com(Zhoucheng Bao); 2020112011013@stu.hznu.edu.cn(Kaige Chen);
2020112011025@stu.hznu.edu.cn(Tingting Pang); zulingw@163.com(Zuling Wang)
ORCID: 0000-0001-5461-800 (Zhoucheng Bao); 0000-0002-0432-4454(Kaige Chen); 0000-0002-0923-8803(Tingting Pang); 0000-0003-
4850-9976 (Zuling Wang)
              Β© 2022 Copyright for this paper by its authors.
              Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0).
              CEUR Workshop Proceedings (CEUR-WS.org)



                                                                                 225
three subpopulations will use different and fixed mutation operators, that is, similar individuals will form
a new population in each generation, so that the mutation operator can act on individuals. A mutation
strategy named DE/current-pbad/1 has also been devised and performed on the top-ranked
subpopulation. At the same time, a parameter adaptive method based on the diversity and fitness value
of subpopulation is proposed, which can better address different problems. Additionally, a population
reduction strategy is proposed to adaptively remove unpromising individuals. Experiments show that the
performance of EMPADE is better or competitive than related method.
    The rest of this paper is arranged as follows. The details of EMPADE are described in Section 2.
Section 3 gives the related experiments and results of EMAPDE algorithm. Finally, we conclude this
paper in Section 4.

2. Proposed Algorithm

   In this section, we introduce the detailed process and strategies of our algorithm. First, initializing a
population randomly, calculating the Euclidean distance of each individual in the population to the best
individual. Then, dividing the population into three sub-populations of the same size according to the
sort, in which the top-ranked subpopulation has a low diversity, and the bottom-ranked subpopulation
has high diversity. The first third uses a devised mutation strategy called DE/current-to-pbad/1, which
can improve the subpopulation diversity. The middle part uses mutation strategy DE/current-to-rand/1,
which is able to converge fast while maintaining diversity. The last third uses mutation strategy
DE/current-to-pbest/1, which can accelerate convergence. The program runs repeatedly until the
termination conditions are met: three populations are run in parallel to search the space. At the end of
each generation, update the control parameters F and CR and employing the proposed adaptive
elimination operation to update the population size. This is followed by re-sorting and redividing the
population into three sub-populations according to Euclidean distance. The overview of the proposed
algorithm is shown in Algorithm I.
   In the following subsections, we will describe the details of a cooperative mechanism strategy based
on multi group mutation in Section 2.1, a parameter adaptive strategy for multi population in Section 2.2,
and an adaptive elimination operation in Section 2.3.

 Algorithm I EMPADE algorithm.
 1) Set the initial population size 𝑁𝑃 , minimum population 𝑁𝑃          and probability of excellent
     individual(pi).
 2) Calculate the Euclidean distance of each individual to the best individual and sort the
     individuals.
 3) Set the initial parameters of each subpopulation mF, mCR, 𝑁𝑃 = 𝑁𝑃 /3, p = 1, 2,3.
 4) Assign different mutation operators to three subpopulations.
 5) Run the following processes circularly until the termination conditions are met.
     a) Information sharing: the individuals with the best fitness values of three sub populations
          are replaced by the best ones in the whole population
     b) The three sub populations Pi evolve in parallel with different mutation mi operators and
          return to𝑃 , 𝑓𝑖𝑑𝑛𝑒𝑠𝑠 , 𝑁𝑃
     c) Calculate the proportional parameters of control mF and mCr of each sub population (π‘šπΉπ‘
          and π‘šπΆπ‘Ÿ ).
     d) Adjust population size.
     e) Sort the individuals according to Euclidean distance, and redivide three subpopulations of
          the same size.
 6) Output results after calculating resource consumption.




                                                    226
2.1.    Cooperative mechanism of adaptive multi group

    When the entire population is divided into three equally sized sub-populations according to their
Euclidean distances to the best individual, the highest ranked population will have the least diversity. If
traditional mutation operators such as DE/rand/1 and DE/rand/2 are used, it cannot help to improve its
diversity, since the Euclidean distance among individuals in the subpopulation is small.
    In view of the fact that and cannot maintain diversity in the subpopulation, a new mutation operator,
called DE/current-to-pbad with archive, is proposed. Compared with the current subpopulation, the
recently explored inferior solutions can improve the diversity of subpopulation. 𝑋      , βˆ’ 𝑋 , is good
for individuals to evolve toward diversity, and random 𝑋 , βˆ’ 𝑋 , makes the individual change value
larger, which is good for jumping out of the local optimum. In DE/current-to-pbad/1, a mutation vector
is generated in the following manner.

 𝑉,    = 𝑋,       +𝐹 βˆ— 𝑋        ,    βˆ’ 𝑋,       + π‘Ÿπ‘Žπ‘›π‘‘ 𝑋     ,    βˆ’ 𝑋    ,     ,     i ∈ 1,2, … , 𝑁𝑃 (1)

   where 𝑋      , represents an individual randomly selected from the last 𝑝𝑖 individuals of the current
sub population, 𝑝𝑖 ∈ (0,1).
   In summary, we calculate and sort the Euclidean distance for each generation, and divide three
subpopulations. And P1 (”DE/current-pbad/1” for the top ranking), P2(”DE/current-to-rand/1” for the
general ranking) and P3(”DE/current-to-pbest/1” for the worst ranking).

2.2.    Parameter adaptive strategy for multi-subpopulation

   The parameter control of each subpopulation is very meaningful for improving the performance of
the algorithm. In EMPADE, a adaptive parameter control mechanism based on multi subpopulation is
designed that control parameters of each subpopulation can be adaptively changed according to the
population diversity and the best fitness value. The adaptive scheme can coordinate the exploration and
exploitation capabilities of each subpopulation to generate better offspring in the process of evolution.
   The updates of π‘šπΉ and π‘šπΆπ‘Ÿ are shown as follows, respectively. (And in the paper, p represents the
 p subpopulation, and p is 1,2,3)

                                                𝑒
                                      π‘šπΉ =              βˆ— π‘šπ‘’π‘Žπ‘› 𝑆𝐹                                      (2)
                                                2
                                                    βˆ—
                                              𝑒
                                   π‘šπΉ =                βˆ— π‘šπ‘’π‘Žπ‘› 𝑆𝐹                                     (3)
                                              2
   where 𝐷𝑆 ∈ (0,1)is the proportion of each sub population to the diversity of all sub-populations,
and 𝑄𝐹 ∈ (0,1) is the ratio of the best fitness value of the current subpopulation to the sum of the best
fitness values of the three sub-populations. When the diversity of population 𝐷𝑆 is small,             (F)
increased and 𝑄𝐹 βˆ— 𝑄𝐹 (CR) decreased, then mF will increase and mCr will decrease to increase the
diversity of subpopulation.

2.3.    Adaptive elimination operation

   The current multi population algorithms often migrate individuals or copy individuals in different
subpopulations, and the poorly performing individuals in the entire population are not eliminated.
Although the diversity of population can be maintained, the overall performance of the algorithm will be
affected. The linear population size reduction scheme proposed in LSHADE was proven to be excellent
scheme for population size adaptation[12][13]. However, the rapid decline of the population at the
beginning may lead to the loss of diversity of the population. Here, we designed a parabolic population


                                                    227
elimination operation to eliminate under-performing individuals at each generation while maintaining
the diversity of population at the early stage. The adaptive elimination operation is calculated as:

                                          𝑁𝑃   βˆ’ 𝑁𝑃
                        𝑃𝑆     = π‘Ÿπ‘œπ‘’π‘›π‘‘                     βˆ— 𝑛𝑓𝑒𝑠 + 𝑁𝑃                              (4)
                                           max _𝑛𝑓𝑒𝑠

   where 𝑁𝑃       represents the minimum population, 𝑁𝑃       represents the original population. The
poor individuals in the whole population will be randomly selected and removed.

3. Experiments &Results

   All the algorithms in the experiment are run on the same computer, and are aimed at the cec2014
dataset. We will run each function 51 times, record its average value and standard deviation, and mark
the best one as bold. The rank sum test of the algorithm takes 0.05 as the confidence interval. The sign
β€˜+’ means better than the correlation algorithm, β€˜-’ means worse than the correlation algorithm, and
β€˜=’means similar to the correlation algorithm. The error of the experimental results is f-f *, where f *
represents the best fitness value of the function.

3.1.    Parameter Settings

   Table I lists the parameter settings in our algorithm, and the parameters in the comparison algorithm
are all the parameters that get the best results according to the relevant papers.
Table I
Experimental parameter setting of EMPADE
      Parameters             𝑁𝑃                   𝑁𝑃              π‘šπΉ             π‘šπΆπ‘Ÿ              pi
          Value                9 *Dim                15             0.5            0.5          0.11

3.2.    The proposed strategies

    The strategy mechanisms are been explored by comparing EMPADE with EMPADE 1, in which a
cooperative mechanism strategy of adaptive multi group mutation operator is employed. EMPADE 2
uses an adaptive elimination operation based on EMPADE 1 and EMPADE uses a parameter adaptive
strategy for multi population based on EMPADE 2. The results on 30D are shown in Table II. Compared
with other algorithms, 25 of the 30 benchmark functions have the best results. The EMPADE could be
significantly better than three algorithms to be compared, especially on unimodal and hybrid functions.
The results therefore indicated the strategies proposed are effective.

3.3.    Comparisons with related algorithms

    EMPADE is compared to four DE variants: JADE[6], EPSDE[7], MPEDE[9], EDEV[10]. And Table
II show the comparison results of the four algorithms. Looking at the experimental results on 30D,
                                                             βˆ’
EMPADE achieves the best mean values on 22 functions      βˆ—
                                                            out of 30 functions. For unimodal functions
F1–F3, EMPADE shows the best performance, the best mean fitness value can be found on all functions.
EPSDE, MPEDE and EDEV can find the best mean fitness value on F2 and F3. And JADE can also find
the best mean fitness value on F2. For the multi-modal benchmark functions F4-F16, compared with
other algorithms, EMPADE still has the best performance and can be more competitive on eleven
functions (F4, F5, F6, F7, F8, F9, F11, F12, F14, F15, F16), and finds the best fitness values on F4, F7
and F8. JADE also has a good performance and can achieve good results on ten functions (F4, F5, F7,
F8, F10, F11, F12, F13, F14, F16), and finds the best fitness values on F4 and F8. EPSDE, MPEDE and
EDEV have competitive performance on one function (F8), three functions (F6, F7, F8) and three
functions (F7, F8, F13), respectively, and they can all find the best mean fitness value on F8. For the
hybrid functions F17-F22, EMPADE is the most competitive, better than other related algorithms in all

                                                  228
functions. For the complex composition functions F23-F30, EDEV achieves the best results, with
excellent performance on 7 functions (F23, F25, F26, F27, F28, F29 and F30). JADE, EPSDE, MPEDE
and EMPADE perform well on one function (F27), three functions (F23, F25, F29), one function (F26)
and two functions (F24, F26), respectively.
  Based on the above experiments and results, it can be concluded that EMPADE is better or
competitive than related algorithms. especially in low dimensional unimodal functions, multimodal
functions and hybrid functions.

Table II
Comparison results within 30D
                 EMPADE 1                     EMPADE 2                     EMPADE
Functions        Mean (Std)                   Mean (Std)                   Mean (Std)
F1               2.28E-05(1.34E-04)           6.67E-07(2.26E-06)           0.00E+00(0.00E+00)
F2               0.00E+00(0.00E+00)           0.00E+00(0.00E+00)           0.00E+00(0.00E+00)
F3               0.00E+00(0.00E+00)           0.00E+00(0.00E+00)           0.00E+00(0.00E+00)
F4               1.26E+00(8.88E+00)           0.00E+00(0.00E+00)           0.00E+00(0.00E+00)
F5               2.04E+01(5.41E-02)           2.02E+01(6.66E-0)            2.01E+01(1.23E-01)
F6               3.04E+00(2.05E+00)           2.47E+00(1.33E+00)           1.17E+00(1.32E+00)
F7               0.00E+00(0.00E+00)           0.00E+00(0.00E+00)           0.00E+00(0.00E+00)
F8               3.22E-05(5.83E-05)           0.00E+00(0.00E+00)           0.00E+00(0.00E+00)
F9               2.22E+01(6.12E+00)           2.19E+01(5.23E+00)           2.34E+01(6.24E+00)
F10              1.45E+01(4.09E+00)           2.45E-01(5.85E-02)           8.04E-02(4.02E-02)
F11              2.25E+03(3.85E+02)           1.70E+03(3.57E+02)           1.72E+03(3.86E+02)
F12              4.65E-01(6.78E-02)           1.69E-01(5.88E-02)           1.61E-01(8.10E-02)
F13              2.38E-01(3.30E-02)           2.06E-01(3.95E-02)           2.18E-01(5.09E-02)
F14              2.22E-01(2.52E-02)           2.14E-01(2.44E-02)           2.17E-01(3.54E-02)
F15              3.75E+00(5.78E-01)           2.54E+00(5.97E-01)           2.58E+00(6.51E-01)
F16              9.86E+00(4.12E-01)           9.30E+00(5.48E-01)           9.21E+00(6.16E-01)
F17              2.64E+02(1.55E+02)           2.64E+02(1.46E+02)           1.81E+02(1.53E+02)
F18              1.08E+01(3.90E+00)           1.09E+01(3.88E+00)           9.73E+00(5.20E+00)
F19              4.03E+00(4.68E-01)           3.25E+00(4.67E-01)           2.95E+00(4.19E-01)
F20              1.15E+01(4.02E+00)           9.57E+00(3.46E+00)           5.81E+00(2.03E+00)
F21              1.19E+02(8.70E+01)           1.04E+02(8.73E+01)           8.55E+01(7.92E+01)
F22              1.02E+02(5.88E+01)           7.08E+01(6.28E+01)           6.36E+01(6.24E+01)
F23              3.15E+02(2.76E-12)           3.15E+02(2.76E-12)           3.15E+02(2.76E-12)
F24              2.24E+02(7.36E-01)           2.24E+02(7.15E-01)           2.24E+02(2.01E+00)
F25              2.03E+02(1.53E-01)           2.03E+02(1.41E-01)           2.03E+02(1.14E-01)
F26              1.00E+02(3.24E-02)           1.00E+02(4.12E-02)           1.00E+02(5.22E-02)
F27              3.59E+02(3.57E+01)           3.59E+02(3.57E+01)           3.81E+02(3.07E+01)
F28              7.61E+02(4.09E+01)           7.57E+02(4.32E+01)           8.03E+02(4.98E+01)
F29              6.87E+02(1.20E+02)           6.97E+02(9.17E+01)           7.15E+02(7.85E-01)
F30              5.72E+02(1.78E+02)           5.63E+02(1.36E+02)           7.41E+02(3.75E+02)
                 20/7/3                       10/17/3                      +/-/=


4. Conclusions

   This paper proposes an adaptive DE based on multi-subpopulation. In the proposed method, we
design an adaptive parameter control scheme, which can comprehensively consider the optimal solution
and diversity of each sub-population. This strategy can improve the robustness of multiple population
based on DE. Further, a mechanism of eliminating poor individuals is proposed to preserve the
population diversity at the early stage as well as improve the convergence speed of the population. In
addition, we propose a mutation strategy to solve the problem of insufficient diversity of the top-ranked
sub-populations when obtaining sub-populations based on Euclidean distance. The experimental results


                                                  229
clearly show the significance of the proposed strategy, and the obtained algorithm can achieve better or
comparable performance than related method.

Table III
Comparison results with related methods in 30D
                    JADE                 EPSDE                 MPEDE                EDEV                EMPADE
  Functions       Mean (Std)           Mean (Std)           Mean (Std)            Mean (Std)           Mean (Std)
     F1       2.36E+03(2.89E+03)   8.12E+04(4.66E+05)   6.19E-06(4.42E-05)    6.33E+03(1.71E+04)   0.00E+00(0.00E+00)
     F2       0.00E+00(0.00E+00)   0.00E+00(0.00E+00)   0.00E+00(0.00E+00)    0.00E+00(0.00E+00)   0.00E+00(0.00E+00)
     F3       2.16E-11(8.08E-11)   0.00E+00(0.00E+00)   0.00E+00(0.00E+00)    0.00E+00(0.00E+00)   0.00E+00(0.00E+00)
     F4       0.00E+00(0.00E+00)   3.03E+00(2.21E+00)   5.67E-09(3.40E-08)    2.23E-01(8.07E-01)   0.00E+00(0.00E+00)
     F5       2.03E+01(3.44E-02)   2.04E+01(3.42E-02)   2.04E+01(5.41E-02)    2.04E+01(5.06E-02)   2.01E+01(1.23E-01)
     F6       9.15E+00(2.59E+00)   1.89E+01(1.67E+00)   7.53E-01(1.02E+00)    6.01E+00(4.15E+00)   1.17E+00(1.32E+00)
     F7       3.98E-12(2.80E-11)   1.35E-03(3.95E-03)   1.45E-04(1.04E-03)    0.00E+00(0.00E+00    0.00E+00(0.00E+00)
     F8       0.00E+00(0.00E+00)   0.00E+00(0.00E+00)   0.00E+00(0.00E+00)    0.00E+00(0.00E+00)   0.00E+00(0.00E+00)
     F9       2.61E+01(4.09E+00)   4.29E+01(6.61E+00)   2.69E+01(6.51E+00)    3.44E+01(4.62E+00)   2.34E+01(6.24E+00)
    F10       7.35E-03(1.43E-02)   2.27E-01(1.83E-01)   1.32E+00(7.00E-01)    1.61E+00(1.14E+01)   8.04E-02(4.02E-02)
    F11       1.73E+03(2.11E+02)   3.49E+03(4.45E+02)   2.43E+03(3.37E+02)    2.95E+03(7.29E+02)   1.72E+03(3.86E+02)
    F12       2.71E-01(2.87E-02)   5.12E-01(5.41E-02)   4.74E-01(8.78E-02)    6.71E-01(1.71E-01)   1.61E-01(8.10E-02)
    F13       2.03E-01(3.62E-02)   2.58E-01(4.38E-02)   2.11E-01(3.85E-02)    2.06E-01(3.19E-02)   2.18E-01(5.09E-02)
    F14       2.28E-01(3.38E-02)   3.00E-01(8.39E-02)   2.37E-01(2.77E-02)    2.31E-01(3.34E-02)   2.17E-01(3.54E-02)
    F15       3.16E+00(3.27E-01)   5.37E+00(7.27E-01)   4.04E+00(5.31E-01)    4.13E+00(4.98E-01)   2.58E+00(6.51E-01)
    F16       9.40E+00(3.83E-01)   1.12E+01(4.09E-01)   1.00E+01(5.37E-01)    1.01E+01(3.76E-01)   9.21E+00(6.16E-01)
    F17       9.97E+02(3.81E+02)   5.00E+04(5.06E+04)   2.13E+02(1.47E+02)    9.09E+02(6.23E+02)   1.81E+02(1.53E+02)
    F18       6.03E+01(2.99E+01)   2.02E+02(3.62E+02)   1.54E+01(7.11E+00)    2.55E+01(1.74E+01)   9.73E+00(5.20E+00)
    F19       4.58E+00(7.63E-01)   1.31E+01(1.24E+00)   3.78E+00(6.85E-01)    4.70E+00(1.56E+00)   2.95E+00(4.19E-01)
    F20       2.20E+03(2.98E+03)   9.07E+01(1.55E+02)   9.06E+00(2.96E+00)    1.50E+01(3.61E+00)   5.81E+00(2.03E+00)
    F21       2.79E+02(1.53E+02)   1.53E+04(3.15E+04)   1.04E+02(9.60E+01)    3.55E+02(1.79E+02)   8.55E+01(7.92E+01)
    F22       1.33E+02(7.57E+01)   2.06E+02(8.50E+01)   1.03E+02(7.17E+01)    7.48E+01(5.44E+01)   6.36E+01(6.24E+01)
    F23       3.15E+02(3.21E-12)   3.14E+02(1.38E-12)   3.15E+02(3.21E-12))   3.14E+02(0.00E+00)   3.15E+02(2.76E-12)
    F24       2.25E+02(3.05E+00)   2.29E+02(5.66E+00)   2.25E+02(2.54E+00)    2.25E+02(2.61E+00)   2.24E+02(2.01E+00)
    F25       2.05E+02(1.70E+00)   2.00E+02(4.42E-01)   2.03E+02(3.01E-01)    2.00E+02(1.38E+00)   2.03E+02(1.14E-01)
    F26       1.02E+02(1.40E+01)   1.00E+02(4.82E-02)   1.00E+02(2.84E-02)    1.00E+02(3.77E-02)   1.00E+02(5.22E-02)
    F27       3.42E+02(4.92E+01)   8.24E+02(1.05E+02)   3.73E+02(4.50E+01))   3.50E+02(4.47E+01)   3.81E+02(3.07E+01)
    F28       7.89E+02(4.09E+01)   3.93E+02(1.19E+01)   8.34E+02(3.28E+01)    3.80E+02(6.16E+00)   8.03E+02(4.98E+01)
    F29       7.77E+02(1.97E+02)   2.14E+02(1.26E+00)   6.82E+02(1.33E+02)    2.14E+02(1.11E+00)   7.15E+02(7.85E-01)
    F30       1.36E+03(4.76E+02)   5.91E+02(1.38E+02)   6.08E+02(1.85E+02)    4.23E+02(1.41E+02)   7.41E+02(3.75E+02)
                   19/8/3               22/4/4                 19/11/0             17/6/7                +/-/=


5. References

[1] R. Storn and K. Price, β€œDifferential evolution: A simple and efficient adaptive scheme for global
    optimization over continuous spaces,” Journal of Global Optimization, vol. 23, no. pp. 1–8,1, 1995.
[2] X. Li, L. Wang, Q. Jiang, and N. Li, β€œDifferential evolution algorithm with multi-population
    cooperation and multi-strategy integration,” Neu-rocomputing, vol. 421, no. 1, pp. 285–302, 2021.
[3] U. Maulik and I. Saha, β€œModified differential evolution based fuzzy clustering for pixel
    classification in remote sensing imagery,” Pattern Recognition, vol. 42, no. 9, pp. 2135–2149,
    2009.
[4] G. Y. Yang, Y. D. Zhao, and K. P. Wong, β€œA modified differential evolution algorithm with fitness
    sharing for power system planning,” IEEE Transactions on Power Systems, vol. 23, no. 2, pp. 514–
    522, 2008.


                                                         230
[5] L. Song and D. Yun, β€œAn improved differential evolution algorithm with local search for
     capacitated vehicle routing problem,” in 2018 Tenth International Conference on Advanced
     Computational Intelligence (ICACI), 2018.
[6] J. Zhang and A. C. Sanderson, β€œJade: adaptive differential evolution with optional external archive,”
     IEEE Transactions on evolutionary computation, vol. 13, no. 5, pp. 945–958, 2009.
[7] R. Mallipeddi and P. N. Suganthan, β€œDifferential evolution algorithm with ensemble of parameters
     and mutation and crossover strategies,” in Swarm, Evolutionary, and Memetic Computing - First
     International Conference on Swarm, Evolutionary, and Memetic Computing, SEM- CCO 2010,
     Chennai, India, December 16-18, 2010. Proceedings, 2010.
[8] W. Yong, Z. Cai, and Q. Zhang, β€œDifferential evolution with composite trial vector generation
     strategies and control parameters,” IEEE Trans- actions on Evolutionary Computation, vol. 15, no.
     1, pp. 55–66, 2011.
[9] G. Wu, R. Mallipeddi, P. N. Suganthan, W. Rui, and H. Chen, β€œDifferential evolution with multi-
     population based ensemble of mutation strategies,” Information Sciences An International Journal,
     vol. 329, no. C, pp. 329–345, 2016.
[10] Guohua, Wu, Xin, Shen, Haifeng, Li, Huangke, Chen, Anping, and Lin, β€œEnsemble of differential
     evolution variants,” Information Sciences: An International Journal, vol. 423, pp. 172–186, 2018.
[11] L. Zhu, Y. Ma, and Y. Bai, β€œA self-adaptive multi-population differential evolution algorithm,”
     Natural Computing, vol. 19, no. 1, pp. 211–235, 2020.
[12] Z. Meng, J. S. Pan, and K. K. Tseng, β€œPade: An enhanced differential evolution algorithm with
     novel control parameter adaptation schemes for numerical optimization,” Knowledge-Based
     Systems, 2019.
[13] J. Wei, Z. Wang, Y. Xu, and Z. Chen, β€œA differential evolution with multi-factor ranking based
     parameter adaptation for global optimization,” in 2021 IEEE Congress on Evolutionary
     Computation (CEC),2021




                                                  231