
Exposure to respirable fractions of crystalline silica quartz dust particles is associated with silicosis, cancer and the development of autoimmune conditions. Early cellular interactions are not well understood, partly due to a lack of suitable technological methods. Improved techniques are needed to better quantify and study high-level respirable crystalline silica exposure in human populations. Techniques that can be applied to complex biological matrices are pivotal to understanding particle-cell interactions and the impact of particles within real, biologically complex environments. In this study, we investigated whether imaging flow cytometry could be used to assess the interactions between cells and crystalline silica when present within complex biological matrices. Using the respirable-size fine quartz crystalline silica dust Min-u-sil® 5, we first validated previous reports that, whilst associating with cells, crystalline silica particles can be detected solely through their differential light scattering profile using conventional flow cytometry. This same property reliably identified crystalline silica in association with primary monocytic cells in vitro using an imaging flow cytometry assay, where darkfield intensity measurements were able to detect crystalline silica concentrations as low as 2.5 µg/mL. Finally, we ultilised fresh whole blood as an exemplary complex biological matrix to test the technique. Even after the increased sample processing required to analyse cells within whole blood, imaging flow cytometry was capable of detecting and assessing silica-association to cells. As expected, in fresh whole blood exposed to crystalline silica, neutrophils and cells of the monocyte/macrophage lineage phagocytosed the particles. In addition to the use of this technique in in vitro exposure models, this method has the potential to be applied directly to ex vivo diagnostic studies and research models, where the identification of crystalline silica association with cells in complex biological matrices such as bronchial lavage fluids, alongside additional functional and phenotypic cellular readouts, is required.
Citation: Bradley Vis, Jonathan J. Powell, Rachel E. Hewitt. Imaging flow cytometry methods for quantitative analysis of label-free crystalline silica particle interactions with immune cells[J]. AIMS Biophysics, 2020, 7(3): 144-166. doi: 10.3934/biophy.2020012
[1] | Qing Wu, Chunjiang Zhang, Mengya Zhang, Fajun Yang, Liang Gao . A modified comprehensive learning particle swarm optimizer and its application in cylindricity error evaluation problem. Mathematical Biosciences and Engineering, 2019, 16(3): 1190-1209. doi: 10.3934/mbe.2019057 |
[2] | Wenqiang Zhang, Chen Li, Mitsuo Gen, Weidong Yang, Zhongwei Zhang, Guohui Zhang . Multiobjective particle swarm optimization with direction search and differential evolution for distributed flow-shop scheduling problem. Mathematical Biosciences and Engineering, 2022, 19(9): 8833-8865. doi: 10.3934/mbe.2022410 |
[3] | Shangbo Zhou, Yuxiao Han, Long Sha, Shufang Zhu . A multi-sample particle swarm optimization algorithm based on electric field force. Mathematical Biosciences and Engineering, 2021, 18(6): 7464-7489. doi: 10.3934/mbe.2021369 |
[4] | Rongmei Geng, Renxin Ji, Shuanjin Zi . Research on task allocation of UAV cluster based on particle swarm quantization algorithm. Mathematical Biosciences and Engineering, 2023, 20(1): 18-33. doi: 10.3934/mbe.2023002 |
[5] | Xin Zhou, Shangbo Zhou, Yuxiao Han, Shufang Zhu . Lévy flight-based inverse adaptive comprehensive learning particle swarm optimization. Mathematical Biosciences and Engineering, 2022, 19(5): 5241-5268. doi: 10.3934/mbe.2022246 |
[6] | Xianzi Zhang, Yanmin Liu, Jie Yang, Jun Liu, Xiaoli Shu . Handling multi-objective optimization problems with a comprehensive indicator and layered particle swarm optimizer. Mathematical Biosciences and Engineering, 2023, 20(8): 14866-14898. doi: 10.3934/mbe.2023666 |
[7] | Dashe Li, Xueying Wang, Jiajun Sun, Huanhai Yang . AI-HydSu: An advanced hybrid approach using support vector regression and particle swarm optimization for dissolved oxygen forecasting. Mathematical Biosciences and Engineering, 2021, 18(4): 3646-3666. doi: 10.3934/mbe.2021182 |
[8] | Xiangyang Ren, Shuai Chen, Kunyuan Wang, Juan Tan . Design and application of improved sparrow search algorithm based on sine cosine and firefly perturbation. Mathematical Biosciences and Engineering, 2022, 19(11): 11422-11452. doi: 10.3934/mbe.2022533 |
[9] | Jinyun Jiang, Jianchen Cai, Qile Zhang, Kun Lan, Xiaoliang Jiang, Jun Wu . Group theoretic particle swarm optimization for gray-level medical image enhancement. Mathematical Biosciences and Engineering, 2023, 20(6): 10479-10494. doi: 10.3934/mbe.2023462 |
[10] | Li-zhen Du, Shanfu Ke, Zhen Wang, Jing Tao, Lianqing Yu, Hongjun Li . Research on multi-load AGV path planning of weaving workshop based on time priority. Mathematical Biosciences and Engineering, 2019, 16(4): 2277-2292. doi: 10.3934/mbe.2019113 |
Exposure to respirable fractions of crystalline silica quartz dust particles is associated with silicosis, cancer and the development of autoimmune conditions. Early cellular interactions are not well understood, partly due to a lack of suitable technological methods. Improved techniques are needed to better quantify and study high-level respirable crystalline silica exposure in human populations. Techniques that can be applied to complex biological matrices are pivotal to understanding particle-cell interactions and the impact of particles within real, biologically complex environments. In this study, we investigated whether imaging flow cytometry could be used to assess the interactions between cells and crystalline silica when present within complex biological matrices. Using the respirable-size fine quartz crystalline silica dust Min-u-sil® 5, we first validated previous reports that, whilst associating with cells, crystalline silica particles can be detected solely through their differential light scattering profile using conventional flow cytometry. This same property reliably identified crystalline silica in association with primary monocytic cells in vitro using an imaging flow cytometry assay, where darkfield intensity measurements were able to detect crystalline silica concentrations as low as 2.5 µg/mL. Finally, we ultilised fresh whole blood as an exemplary complex biological matrix to test the technique. Even after the increased sample processing required to analyse cells within whole blood, imaging flow cytometry was capable of detecting and assessing silica-association to cells. As expected, in fresh whole blood exposed to crystalline silica, neutrophils and cells of the monocyte/macrophage lineage phagocytosed the particles. In addition to the use of this technique in in vitro exposure models, this method has the potential to be applied directly to ex vivo diagnostic studies and research models, where the identification of crystalline silica association with cells in complex biological matrices such as bronchial lavage fluids, alongside additional functional and phenotypic cellular readouts, is required.
In the past few decades, many meta-heuristic algorithms have been proposed to handle complex optimization problems, which are usually difficult to solve with traditional optimization methods relying on gradients. Meta-heuristic algorithms are problem-solving methods that mimic natural or abstract processes to efficiently explore solution spaces and find near-optimal solutions for complex optimization problems. Kennedy and Eberhart [1] developed particle swarm optimization algorithm (PSO), a population-based optimization algorithm inspired by social behavior and collective intelligence, where particles iteratively update their positions based on their own experience and the information from the best-performing particles in the swarm. Holland [2] proposed genetic algorithm (GA) that imitates the process of natural selection and genetic evolution to iteratively search for optimal solutions to complex problems. Wang et al. [3] proposed a monarch butterfly optimization algorithm (MBO) that imitates the migration behavior of monarch butterflies, where the search process is guided by a combination of exploration and exploitation strategies to find optimal solutions. Li et al. [4] presented a slime mold algorithm (SMA) that simulates the foraging behavior of slime mold to efficiently explore and exploit search spaces by forming networks of interconnected pathways. Wang [5] presented a moth search algorithm (MSA) that mimics the behavior of moths to navigate toward the optimal solution by adjusting their flight direction based on attraction and repulsion forces. Yang et al. [6] proposed a hunger games search (HGS) algorithm, where individuals compete for resources and survival. Ahmadianfar et al. [7] presented Runge-Kutta optimizer (RUN) that utilizes the logic of slope variations computed by the Runge-Kutta method to perform active exploration and exploitation for global optimization with an enhanced solution quality mechanism to avoid local optima and improve convergence speed. Tu et al. [8] proposed colony predation algorithm (CPA) that mimics the collective hunting behavior of predators to explore and exploit the solution space for optimization efficiently. Ahmadianfar et al. [9] presented INFO which utilizes the weighted mean concept, enhanced updating rules, vector combining and local search to optimize various problems. Heidari et al. [10] proposed Harris hawks optimization (HHO), which is inspired by the hunting behavior of Harris hawks, utilizing different strategies such as exploration, exploitation and cooperative hunting to solve optimization problems. Su et al. [11] proposed a RIME optimization algorithm, a nature-inspired metaheuristic algorithm based on the freezing process of water droplets in rime formation, utilizing dynamic freezing and melting mechanisms for efficient optimization. Compared with conventional algorithms, these meta-heuristic algorithms can handle non-differentiable, multimodal, hybrid and high-dimensional problems without requiring explicit gradient information.
After they are proposed, many researchers successfully solve complex optimization problems by using them in intelligent systems. Li et al. [12] introduced ACO-S, a bionic optimization algorithm-based dimension reduction method specifically designed for high-dimensional datasets like microarray data, demonstrating its ability to generate compact and informative gene subsets with high classification accuracy compared to existing bionic optimization algorithms. Thawkar et al. [13] presented a hybrid feature selection method, combining the butterfly optimization algorithm (BOA) and the ant lion optimizer (ALO) to effectively predict the benign or malignant status of breast tissue using mammogram images. Chakraborty et al. [14] developed a computational tool using a modified whale optimization method (WOA) to quickly and accurately determine the severity of COVID-19 illness through chest X-ray image analysis. Gao et al. [15] introduced a multi-objective optimization model and an improved genetic algorithm for cooperative mission assignment of heterogeneous UAVs, considering the balance between mission gains and UAV losses. Yu et al. [16] presented a multi-objective optimization problem and an improved genetic algorithm for cooperative mission planning of heterogeneous UAVs in cross-regional joint operations, considering makespan minimization, value expectation maximization, flexible base return and ammunition inventory constraints. Li et al. [17] investigated the scheduling problem of maintenance and repair tasks for carrier-based aircraft fleets in hangar bays, proposing a model and an improved teaching-learning-based optimization algorithm to enhance maintenance efficiency and reduce manual scheduling burden.
Among the algorithms mentioned above, PSO is popular because of its fast convergence, few parameters and easy programming implementation. However, the performance of the standard PSO depends mainly on its parameter settings [18] and different parameter settings may lead to different convergence outcomes. In addition, the update strategy of particles during the search process is mainly aimed at the global best, which leads the population to lose diversity in the population and increases the risk of being trapped in local optima [19,20,21]. To address these drawbacks and improve the performance of PSO, many PSO variants have been developed. Generally, research on PSO variants can be classified into three parts: parameter modification, population topology structure modification and novel learning strategy.
1) Parameter modification. PSO has some parameters such as inertia weight and acceleration coefficients. They have been modified in many ways. Shi and Eberhart [22] introduced the concept of inertia weight in PSO to improve velocity control and search effectiveness. Afterward, various adjustment strategies were advocated, including defining it as a time-varying function called global PSO (GPSO) [23]. Liu et al. [24] proposed a PSO algorithm with chaos inertia weight to enlarge the search range and increase population diversity. Chen et al. [25] proposed a hybrid PSO with sine cosine acceleration coefficients where the performance of PSO is improved by introducing sine and cosine acceleration coefficients. By using these strategies, the performance of PSO can be improved to some extent.
2) Population topology structure modification. In the standard PSO, the population size is fixed. A larger population size may increase computation costs, resulting in slow convergence, while a smaller size may converge in local optima. Some researchers have improved the performance of PSO by modifying the population topology structure. Kennedy and Mendes [26] proposed the local version of PSO (LPSO) with ring topology, in which different topologies based on neighborhood connections are used to maintain the diversity of the population. Liang and Suganthan [27] presented a dynamic multi-swarm PSO (DMS-PSO), in which the population comprises many small subgroups and a dynamic strategy is used to maintain population diversity.
3) Novel learning strategy. A well-designed learning strategy can enhance the efficiency and effectiveness of the PSO algorithm. Liang et al. [28] proposed the CLPSO, in which each particle updates its velocity by incorporating historical best information from different particles. This approach enables particles to acquire valuable information and improves their search behavior. Wang et al. [29] proposed a PSO-HLM algorithm that balances exploration and exploitation by utilizing a hybrid learning model and a multi-pools fusion strategy, improving convergence speed and avoiding local optima. Zhou et al. [30] proposed a LFIACL-PSO algorithm that enhances PSO by incorporating Lévy flight-based inverse learning, comprehensive learning strategy, ring-type topology and adaptive update of acceleration coefficients.
In this paper, an improved PSO combined with double-chaos search (DCS-PSO) is proposed. Chaotic dynamics is a branch of nonlinear dynamical systems. Chaos motion exhibits ergodicity, randomness and regularity properties and can traverse all states without repetition. For optimization problems, these properties of chaos can be used as a mechanism to avoid getting trapped in local optima during the search process. In DCS-PSO, we utilize two distinct chaotic systems to explore the search space globally and narrow it down, enabling PSO to focus on the vicinity of the optimal solution and alleviate its premature convergence caused by getting trapped in local optima. Additionally, we incorporate a chaotic motion search in parallel with the PSO search. If the solution obtained through the chaotic search outperforms the PSO's global best, it will be replaced with the solution from the chaotic motion to enhance population diversity. This strategy prevents all particles in the PSO from relying solely on the global optimal particle for learning, thereby facilitating escape from local optima.
The main contributions of this paper are described as follows.
1) Introducing a novel approach, DCS-PSO, which combines the strengths of PSO and chaotic dynamics.
2) Addressing the issue of premature convergence by narrowing down the search space using double-chaos search and enhancing population diversity through the logistic map, resulting in improved convergence accuracy and speed.
3) Conducting extensive experimental evaluations to demonstrate the effectiveness and superiority of DCS-PSO compared to standard PSO and other optimization algorithms commonly used in the field.
The rest of this paper is organized as follows. Section 2 introduces the related work. Section 3 describes the proposed DCS-PSO. Section 4 is the experimental studies. Conclusion and future work are provided in Section 5.
PSO is an optimization algorithm based on swarm intelligence inspired by collective behavior in ecosystems such as bird foraging. In PSO, the problem is treated as an optimization problem of finding the optimal solution in a D-dimensional search space, where each particle i in the population represents a possible solution and has two vectors: a position vector Xi=[xi,1,xi,2,⋅⋅⋅,xi,D] and a velocity vector Vi=[vi,1,vi,2,⋅⋅⋅,vi,D]. In searching for the optimal solution, particles move and interact with each other in the search space. The particle i adjusts its position based on the optimal position vector pbesti=[pbesti,1,pbesti,2,⋅⋅⋅pbesti,D] it has reached in the past and the globally optimal position vector gbest=[gbest1,gbest2,⋅⋅⋅,gbestD] discovered by the entire swarm. PSO attempts to find a balance between local search and global search in this way. The update rules for the velocity and position of particle i in the d th dimension are as follows:
v(n+1)i,d=v(n)i,d+c1×r1×(pbesti,d−xi,d)+c2×r2×(gbestd−xi,d), | (2.1) |
x(n+1)i,d=x(n)i,d+v(n+1)i,d | (2.2) |
where c1 and c2 denote the acceleration coefficients, usually set to 2, r1 and r2 are two uniformly distributed values in the range [0, 1]. The parameters Vmin and Vmax can be used to determine an appropriate range of velocity values for each particle to prevent it from wandering too far outside the search space.
The global minimum optimization problem (2.3) for a D-dimensional continuous object can be described as:
min f(X), X=[x1,⋅⋅⋅,xD], s.t. xd∈[ad,bd], d=1,⋅⋅⋅,D | (2.3) |
The standard PSO can be described in the following steps:
Step 1: Initialize N particles. For each particle i in the population:
Step 1.1: Initialize Xi randomly in the range [a,b].
Step 1.2: Initialize Vi randomly in the range [Vmin,Vmax].
Step 1.3: Evaluate the fitness value fi.
Step 1.4: Set pbesti as the particle's current position and gbest as the particle with the best fitness value in the population.
Step 2: Loop until the stop condition is satisfied:
Step 2.1: Update the velocity Vi and position Xi for each particle i according to Eqs (2.1) and (2.2).
Step 2.2: Evaluate the fitness value fi of each particle i.
Step 2.3: Update pbesti if fi<fpbesti,∀i⩽N for each particle i.
Step 2.4: Update the global optimal position gbest fgbest⩽fi, ∀i⩽N.
Chaos is one of the most exciting properties exhibited by the brain [31]. It is the phenomenon of complex, unpredictable and random-like behavior arising from simple deterministic nonlinear systems rather than a state of disorder [32]. Researchers have recently developed many optimization algorithms based on chaos theory to address complex optimization problems [33,34,35].
COA was first proposed by Li and Jiang [32,35] in 1997. Its core idea is to introduce the chaotic state into optimization variables via the carrier wave method, map the traversal range of chaos motion to the value range of optimization variables and then search for the global optimal solution using chaotic dynamics instead of random search. The chaotic dynamics equation chosen in COA is the Logistic map [36], which is defined as follows:
x(n+1)=μx(n)(1−x(n)), 0 ⩽x(0)⩽1 | (2.4) |
where μ is the control parameter, n=0,1,2,.... Although the above equation is definite, it exhibits chaotic dynamics when μ=4 and x(0)∉{0,0.25,0.5,0.75,1}, the property that a minute difference in the initial value of the chaotic variable would lead to significant differences in output. This is the sensitive dependence on initial conditions, the basic characteristic of chaos. This property enables the trajectory of chaotic variables to traverse the entire search space. The process of COA can be defined through the following equation.
cx(k+1)d=4cx(k)d(1−cx(k)d), d=1,2,⋅⋅⋅,D, | (2.5) |
where cxd is the d th chaotic variable and k represents the number of iterations. Obviously, cx(k)d is distributed in the range (0, 1) when the initial cx(0)d∈(0,1) and cx(0)d∉{0.25,0.5,0.75}.
For the optimization problem (2.3), the procedures of COA are as follows.
Step 1: Set k=0, and initialize D chaotic variables cx(0)d in the range (0, 1) with very small differences. The current optimal function value f∗ is initialized to a large value, with X∗=[x∗1,x∗2,⋅⋅⋅x∗D] representing the corresponding optimal variables.
Step 2: Map the D chaotic variables cxd to the D optimization variables xd of the optimization problem (2.3) through the carrier wave method of Eq (2.6).
x(k)d=cx(k)d(bd−ad)+ad | (2.6) |
Step 3: Iterative search using chaotic variables.
If f(X(k))⩽f∗, then f∗=f(X(k)), X∗=X(k). Else if f(X(k))>f∗, then give up X(k).
Step 4: Let k=k+1 and continue the chaotic traversal through Eq (2.5).
Step 5: Repeat steps 2 to 4. If the best function value f∗ remains unchanged after several iterations, the second carrier wave is applied according to Eq (2.7).
x(k+1)d=x∗d+α⋅cx(k+1)d | (2.7) |
where x∗d is the current optimal solution, α is an adjusting constant which can be less than 1 and α⋅cx(k+1)d generates i chaotic states with small ergodic ranges around x∗d.
Step 6: Continue the iterative search using the chaotic variables after the second carrier wave until the stop criterion is satisfied.
As a simple and efficient optimization algorithm, COA has garnered the attention of many researchers for further improvement. Its primary disadvantage is that the sensitive dependence of chaotic motion on the initial values will cause the algorithm's instability and a large number of traversal searches are required before the search space is reduced. Moreover, the lack of an explicit search stop condition results in the need for parameter tuning according to specific problems, limiting the algorithm's generality. Therefore, many variants of COA aimed at these problems have emerged.
Xiu et al. [37] proposed a double-chaos optimization algorithm (DCOA) that uses two different chaotic systems to explore the search space independently and in parallel and narrows down the search space when the distance between the optimal positions of the two systems satisfies specific criteria. Therefore, DCOA overcomes some of the deficiencies of COA. Additionally, DCOA provides a condition for narrowing the search space, which improves the algorithm's generality, avoids multiple blind searches and reduces its running cost. Liang and Gu [38] developed a chaos optimization algorithm based on parallel computing (PCOA), which performs parallel searches using several different sets of chaotic variables, effectively reducing the sensitive dependence of chaotic motion on initial values.
Based on in-depth research of the chaos algorithms and the PSO algorithm, we combined the advantages of the two types of algorithms and proposed DCS-PSO. The details are as follows:
Although a single chaos search mechanism is easy to implement and can avoid being trapped in local optima, it often requires a large number of iterations to narrow down the search space to ensure that the optimal solution is within the narrowed area.
Additionally, the stop condition of the search needs to be adjusted based on the specific problem, limiting the algorithm's generality. In contrast, the DCOA determines the requirement for narrowing the space based on the distance between the optimal positions found by each chaotic system, which improves the algorithm's generality. However, there is a risk that the algorithm may continuously narrow down the search space, making it difficult to converge accurately to a point and potentially excluding the global optimal solution from the search space. Nonetheless, experimental verification (see Section 4.1) suggests that the double-chaos search mechanism can narrow the search space to the neighborhood of the global optimum after only one sufficient search.
In the double-chaos search mechanism, we select two chaotic systems: The logistic map and the tent map. The orbit points of the logistic map are distributed near the edges, which to some extent can explain why chaotic motion has the advantage of escaping from local optima [39]. However, several breakpoints {0,0.25,0.5,0.75,1} in the logistic map make it difficult to traverse specific chaos states during the search process. The orbit points of the tent map are distributed more evenly. Its equation is defined in Eq (3.1). Figure 1(a) and (b) show the scatter plots of the two chaos maps for 2000 iterations, respectively and (c) and (d) show their chaotic dynamics.
x(n+1)={x(n)/β, 0<x(n)⩽β(1−x(n))/(1−β), β<x(n)⩽1 | (3.1) |
where β=0.4.
The double-chaos search mechanism is illustrated in Figure 2. In the search space A, two chaotic systems cx and cy are used to perform independent parallel searches. When the distance between the optimal solutions cx∗ and cy∗ found by cx and cy is small enough, namely,
‖X∗−Y∗‖2<γ‖a−b‖2 | (3.2) |
where X∗ and Y∗ represent the corresponding optimization variables, γ∈(0,0.25] and is set to 0.15 in this study, the search space is narrowed down from A to B according to Eqs (3.3) and (3.4).
ad=max(ad,(min(x∗d,y∗d)−ξ⋅γ⋅‖X∗−Y∗‖2)) | (3.3) |
bd=min(bd,(max(x∗d,y∗d)+ξ⋅γ⋅‖X∗−Y∗‖2)) | (3.4) |
where ξ∈[1,2] and it is set to 1.5 in this study.
Based on the chaos search mechanism, we propose a two-stage iteration strategy named DCS-PSO. The algorithm employs the double-chaos search mechanism for global exploration and narrowing the search space and then uses the improved PSO for further fine search in the narrowed search space.
For the D-dimensional optimization problem (2.3), the detailed steps of the first stage of DCS-PSO are listed in Algorithm 1:
Algorithm 1 Space reduction based on double-chaos search mechanism. |
1. Initialize two chaotic systems cx(0) and cy(0) in the range (0, 1); |
2. Initialize the optimal function values of the two chaotic systems f∗x and f∗y with large values; |
3. Initialize the corresponding optimization variables X∗ and Y∗; |
4. set k=0; |
5. while True do: |
6. Map cx(k) and cy(k) to optimization variables X(k) and Y(k); |
7. Update f∗x=f(X(k)) and X∗=X(k), if f(X(k))<f∗x; |
8. Update f∗y=f(Y(k)) and Y∗=Y(k), if f(Y(k))<f∗y; |
9. if k>h then |
10. if ‖X∗−Y∗‖2<γ‖a−b‖2 then |
11. Narrow down the search space [a,b] according to Eqs (3.3) and (3.4); |
12. break |
13. end if |
14. end if |
15. Update cx(k) and cy(k) according to Eqs (2.5) and (3.1), respectively; |
16. k++; |
17. end while |
18. Output the narrowed search space a=[a1,a2,⋅⋅⋅ad],b=[b1,b2,⋅⋅⋅,bd] |
19. Output cbest, one of X∗ and Y∗ that has the smallest fitness value |
The flowchart of the first stage of DCS-PSO is shown in Figure 3(a).
In conventional PSO, all particles learn from the global optimal particle in the population to update their positions and velocities. This mechanism is advantageous in achieving fast convergence, but it also leads to the need for more diversity and often causes premature convergence due to falling into local optimum [40]. The conventional PSO may require longer to find the optimal solution when dealing with multimodal optimization problems with a large search space.
Although the risk of being trapped in local optima can be reduced by narrowing the search space, PSO may still fall into a local optimum in the reduced search space. To enhance the population's diversity, in the narrowed search space, we perform a chaos search on the best solution found by the two chaotic systems using the logistic map. The best solution found by the logistic map and the population guides the algorithm to converge toward the global optimum. The velocity of particle i on dimension d is updated according to the following rules.
best={cbest, if f(cbest)⩽f(gbest)gbest, else | (3.5) |
v(n+1)i,d=w×v(n)i,d+c1×r1×(pbesti,d−xi,d)+c2×r2×(bestd−xi,d) | (3.6) |
where best denotes the best solution between chaos search and population search, f(cbest) is the optimal fitness value found by chaos search and cbest is the corresponding solution.
The detailed steps of the second stage of DCS-PSO are listed in Algorithm 2:
Algorithm 2 Improved PSO combined with chaos |
1. Input the narrowed search space [a,b] and cbest |
2. Transform cbest to chaotic variables cx; |
3. Initialize N particles' position and velocity within the new search space; |
4. Evaluate the fitness value of every particle; |
5. Update pbesti and gbest; |
6. for i=1: Max_Iter do |
7. Implement chaos search globally and update cbest via logistic map; |
8. if f(cbest)⩽f(gbest) then |
9. best=cbest; |
10. else best=gbest; |
11. end if |
12. Update particles' velocity and position according to Eqs (3.6) and (2.2), respectively; |
13. Evaluate the fitness value of every particle; |
14. Update pbesti and gbest; |
15. end for |
The flowchart of the second stage of DCS-PSO is shown in Figure 3(b).
It is worth noting that many researchers have attempted to introduce the chaos search mechanism into intelligent optimization algorithms [39,41]. The most common strategy is to use chaotic local search (CLS) after each iteration to conduct a fine search near the current optimum. Although this approach improves the convergence accuracy of the algorithm, it needs to address the problem of the population falling into local optima. Therefore, in addition to CLS, other strategies are often required to maintain the diversity of the population. However, some of these strategies may sacrifice the performance of convergence speed.
In this chapter, simulation experiments are conducted to validate the performance of DCS-PSO. We first conduct experiments using six 2-dimensional benchmark functions to validate the effectiveness of the DCS-PSO search space narrowing mechanism. After that, the performance of the DCS-PSO algorithm is compared with five typical PSO variants placed on ten benchmark functions in terms of solution accuracy, statistical significance test and convergence speed. Finally, to verify the performance of DCS-PSO in complex and multimodal conditions, an authoritative test suite CEC2017 is selected in this paper.
All simulation experiments were conducted on the same PC. The configuration information of the PC is shown in Table 1.
PC configuration | Information |
CPU | Inter(R) Core(TM) i7-8700 |
Frequency | 3.2GHz CPU |
RAM | 16.0GB |
Operating system | Windows 10(64Bit) |
Language | Python 3.9.13 |
To validate the effectiveness of the search space narrowing mechanism in DCS-PSO, we conducted experiments using six 2-dimensional benchmark functions [38,42]. DCS-PSO was compared with other chaos-based algorithms, namely COA, DCOA and PCOA. The narrowed search space can be visually observed through the contour map of the 2-dimensional function. The six benchmark functions are listed below.
(1) Jong function:
fJ=100⋅(x21−x2)2+(1−x1)2, −2.048⩽x1,x2⩽2.048. | (4.1) |
The function has a global minimum of 0 at (1, 1).
(2) Camel function:
fC=(4−2.1⋅x21+x413)⋅x21+x1⋅x2+(−4+4⋅x22)⋅x22, −2⩽x1,x2⩽2. | (4.2) |
The function has a global minimum of −1.031628 at two points (−0.0898, 0.7126) and (0.0898, −0.7126).
(3) GP-Goldstein-Price:
fG(x)=[1+(x1+x2+1)2(19−14x1+3x21−14x2+6x1x2+3x22)] ×[30+(2x1−3x2)2(18−32x1+12x21+48x2−36x1x2+27x22)], −2⩽x1,x2⩽2. | (4.3) |
The function has a global minimum of 3 at (0, −1) and there are four local minima in the minimization region.
(4) BR-Branin:
fB(x)=(x2−5.14π2x21+5πx1−6)2+10(1−18π)cosx1+10, −5⩽x1⩽10, 0⩽x2⩽15. | (4.4) |
The global minimum of this function is approximately equal to 0.398 at three points (−3.142, 12.275), (3.142, 2.275) and (9.425, 2.425).
(5) RA-Rastrigin:
fR(x)=x21+x22−cos18x1−cos18x2, −1⩽x1,x2⩽1 | (4.5) |
The function has a global minimum of −2 at (0, 0) and its lattice structure has about 50 local minima.
(6) SH-Shuber:
fs(x)={5∑i=1icos((i+1)x1+i)}{5∑i=1icos((i+1)x2+i)}, −10⩽x1,x2⩽10. | (4.6) |
The function has 760 local minima, 18 of which are global with a value of −186.7309.
The RA-Rastrigin function and the SH-Shuber function are relatively complex and have a large number of local minima, which can effectively test the algorithm's ability to escape from the local optima. Their 3D stereograms are shown in Figure 4.
To ensure the fairness of the experimental results, each algorithm will run on each benchmark function independently 50 times and the mean execution time, the mean optimal values and the standard deviations will be provided. The specific parameter settings of all algorithms are listed in Table 2.
Algorithm | Parameters |
COA | max_iter1=105, max_iter2=1.5×105, α=0.3 |
DCOA | h=3000, γ=0.25, ξ = 1.5 |
PCOA | n_group=3, max_iter1=105, max_iter2=1.5×105, α=0.3 |
DCS-PSO | N=20, w=0.4, c1=c2=2, Vmax=0.2×Range, max_iter=1000 |
Table 3 shows the execution time (E.T.), the mean optimal values (Mean) and standard deviations (S.D.) of COA, DCOA, PCOA and DCS-PSO. The best results of each algorithm on each test function are indicated in bold and the second-best results are marked with an underline.
Function | Criteria | COA | DCOA | PCOA | DCS-PSO |
fJ | E.T. | 30.8898 | 0.5753 | 17.2234 | 0.2793 |
Mean | 0.0000 | 0.0005 | 0.0000 | 0.0000 | |
S.D. | 0.0000 | 0.0014 | 0.0000 | 0.0000 | |
fC | E.T. | 7.6564 | 0.6478 | 10.3990 | 0.3610 |
Mean | −1.0316 | −1.0312 | −1.0316 | −1.0316 | |
S.D. | 0.0000 | 0.0012 | 0.0000 | 0.0000 | |
fG | E.T. | 14.1867 | 1.0656 | 14.0814 | 0.6097 |
Mean | 3.0022 | 3.0100 | 3.0024 | 3.0000 | |
S.D. | 0.0022 | 0.0126 | 0.0022 | 0.0000 | |
fB | E.T. | 1.5663 | 0.8653 | 5.7455 | 0.4513 |
Mean | 0.3983 | 0.4003 | 0.3980 | 0.3979 | |
S.D. | 0.0003 | 0.0039 | 0.0001 | 0.0000 | |
fR | E.T. | 4.4150 | 0.6652 | 5.8342 | 0.4366 |
Mean | −1.9993 | −1.9652 | −1.9992 | −2.0000 | |
S.D. | 0.0006 | 0.0558 | 0.0008 | 0.0000 | |
fS | E.T. | 5.8609 | 2.1743 | 13.8397 | 2.0615 |
Mean | −186.6428 | -186.4517 | −186.6861 | −186.7309 | |
S.D. | 0.0711 | 0.4203 | 0.0680 | 0.0000 |
Based on the mean values in Table 3, the proposed DCS-PSO outperforms the other three chaos-based methods. Compared with COA, DCOA and PCOA, DCS-PSO achieves superior results on five functions (fJ, fC, fG, fR and fS). COA performs well in solving simple optimization problems and acquires the theoretically optimal solutions on both functions fJ and fC. However, it performs poorly on function fS, which has multiple local minima, possibly due to the challenge of attaining certain chaos states in a large search space. The convergence accuracy of DCOA is relatively poor. As mentioned earlier, DCOA may exclude the optimal solution while continuously narrowing the search space. However, compared with COA and PCOA, DCOA shows a shorter convergence time, which reflects its efficiency in narrowing the search space. The overall performance of PCOA is comparable to COA and it is the only algorithm that converges to the global optimum on function fB. DCS-PSO has a faster convergence speed due to the combination of double-chaos search and PSO. The S.D. values in Table 3 show that DCS-PSO has the smallest S.D. values among the six benchmark functions, indicating that DCS-PSO is more stable and reliable.
Figure 5(a) and (b) show the contour plots of functions fR and fS projected onto the XOY plane, respectively. The red rectangle in each plot represents the narrowed search space, demonstrating the ability of DCS-PSO to reduce the search space of complex optimization problems. These plots visually display that DCS-PSO can accurately converge to the vicinity of the global optimum, validating the effectiveness of the proposed method.
In this subsection, five representative PSO variants are selected for comparison. The first one is the PSO with inertia weight (PSO-w) [22], which can better control the search process and is characteristic of fast convergence when dealing with unimodal problems. The second one is the global version of PSO (GPSO) [23], whose inertia weight is decreased linearly from 0.9 to 0.4, improving the population's global search ability. The third one is the local version of PSO (LPSO) [26] with ring topology, in which different topologies based on neighborhood connections are used to maintain the diversity of the population. The fourth one is the Comprehensive Learning PSO (CLPSO) [28], in which each particle updates its velocity by learning the historical best information of different particles. The fifth is the dynamic multi-swarm PSO (DMS-PSO) [27], in which the population comprises many small subgroups and a dynamic strategy is used to maintain population diversity. CLPSO and DMS-PSO are designed to improve the performance of PSO on multimodal problems.
To investigate the performance of DCS-PSO on different optimization problems, ten benchmark functions are tested. Based on the properties of these functions, they are divided into two groups. The first group consists of 5 unimodal benchmark functions with only one global optimum value, mainly used to test the convergence accuracy and speed. The second group consists of 5 multimodal benchmark functions used to test the algorithm's global search ability and its ability to escape from local optima. These benchmark functions are widely used in the literature [25,43,44]. The details of these benchmarks are shown in Table 4.
Name | Functions | Range | Optimum |
Sphere | f1(x)=D∑i=1x2i | [−100, 100] | 0 |
Schwefel's 2.22 | f2(x)=D∑i=1|xi|+D∏i=1|xi| | [−10, 10] | 0 |
Schwefel's 1.2 | f3(x)=D∑i=1(i∑j=1xj)2 | [−100, 100] | 0 |
Rosenbrock | f4(x)=D−1∑i=1[100(xi+1−x2i)2+(xi−1)2] | [−10, 10] | 0 |
Noise | f5(x)=D∑i=1ix4i+random[0,1) | [−1.28, 1.28] | 0 |
Rastrigin | f6(x)=D∑i=1(x2i−10cos(2πxi)+10) | [−5.12, 5.12] | 0 |
Ackley | f7(x)=−20exp(−0.2√1DD∑i=1x2i)−exp(1DD∑i=1cos(2πxi))+20+e | [−32, 32] | 0 |
Griewank | f8(x)=14000D∑i=1x2i−D∏i=1cos(xi√i)+1 | [−600, 600] | 0 |
Penalized 1 | f9(x)=πD{10sin2(πy1)+D−1∑i=1(yi−1)2[1+10sin2(πyi+1)]+(yD−1)2}+D∑i=1u(xi,10,100,4) | [−50, 50] | 0 |
Penalized 2 | f10(x)=110{sin2(3πx1)+D∑i=1(xi−1)2[1+sin2(3πxi+1)]+(xD−1)2[1+sin2(2πxD)]}+D∑i=1u(xi,5,100,4) | [−50, 50] | 0 |
In the first group, the global optimal position of function f4 is [1]D.The global optimal position of the remaining four functions is [0]D and the theoretical optimal values of the five functions are 0. f1 is the sphere function and is easy to solve. f2 and f3 are the Schwefel functions. f4 is the Rosenbrock function, it becomes a simple multimodal problem when its dimension is more than three and has a narrow gap between the perceived local optima and the global optimum. f5 is the noise function, which is challenging to find the global optimum because there is a perturbation factor.
In the second group, the global optimal positions of all functions are [0]D and the global optimal values are 0. Rastrigin's function (f6) has a large number of local optima. f7 has a narrow global basin and a large number of local optima. f8 is the Griewank's function. f9 and f10 are the penalized functions.
In f9, the yi is updated by yi=1+14(xi+1).
In f9 and f10, u(xi,a,k,m)={k(xi−a)m, xi>a 0, −a⩽xi⩽ak(−xi−a)m, xi<−a.
In the experiments, a population size of 40 and a maximum number of iterations of 2000 were used for all algorithms. The parameters for the five PSO variants were kept consistent with their respective literature. The benchmark functions were tested with dimensionalities of 10, 30 and 50 and each benchmark function was independently run 50 times in different dimensions to ensure the fairness of the experimental results. The error between the best optimum found by the algorithm and the actual global optimum was recorded and the mean error and standard deviation were displayed. The specific parameter settings for all algorithms can be found in Table 5.
Algorithm | Parameters |
PSO-w | w:0.4, c1=c2=2.0, Vmax=0.2×Range |
GPSO | w:0.9∼0.4, c1=c2=2.0, Vmax=0.2×Range |
LPSO | χ=0.7298, c1=c2=1.49445, Vmax=0.2×Range |
CLPSO | w:0.9∼0.4, c=1.49445, Vmax=0.2×Range, m=7 |
DMS-PSO | χ=0.7298, c1=c2=1.49445, Vmax=0.2×Range, M=4, R=10 |
DCS-PSO | w:0.4,c1=c2=2.0,Vmax=0.2×Range, h=3000,γ=0.15,ξ=1.5 |
Furthermore, the nonparametric Wilcoxon rank-sum tests were executed for each benchmark function to determine whether the results obtained by DCS-PSO are significantly different from the best results achieved by the other algorithms
Tables 6, 7 and 8 display the means (Mean) and standard deviations (S.D.) of six algorithms on 10 test functions with dimensions of 10, 30 and 50, respectively. The k value presented in the second column of each table is the outcome of the nonparametric Wilcoxon rank-sum tests on our proposed algorithm and the second-best algorithms (indicated with underlining in tables). A k value of 1 indicates a significant difference in the performances of the two algorithms with 95% certainty, while a k value of 0 suggests no significant difference. Moreover, Table 7 ranks the algorithms according to their mean solution accuracy. The best results for each benchmark function are denoted in bold and the second-best results are underlined.
Function | k | Criteria | PSO-w | GPSO | LPSO | CLPSO | DMS-PSO | DCS-PSO |
f1 | 1 | Mean | 1.84 × 10−144 | 1.47 × 10−73 | 2.78 × 10−169 | 2.18 × 10−21 | 6.30 × 10−62 | 2.35 × 10−151 |
S.D. | 9.13 × 10−144 | 1.04 × 10−72 | 0.00 | 5.44 × 10−21 | 4.25 × 10−61 | 7.75 × 10−151 | ||
f2 | 1 | Mean | 1.89 × 10−81 | 6.24 × 10−38 | 7.59 × 10−09 | 4.30 × 10−14 | 3.45 × 10−20 | 9.65 × 10−84 |
S.D. | 1.26 × 10−80 | 2.99 × 10−37 | 3.80 × 10−08 | 1.39 × 10−13 | 5.86 × 10−20 | 4.81 × 10−83 | ||
f3 | 1 | Mean | 1.22 × 10−44 | 9.93 × 10−17 | 6.11 × 10−13 | 2.18 × 10−20 | 6.13 × 10−10 | 1.21 × 10−48 |
S.D. | 7.42 × 10−44 | 3.36 × 10−16 | 4.32 × 10−12 | 3.25 × 10−20 | 3.33 × 10−09 | 2.38 × 10−48 | ||
f4 | 1 | Mean | 2.79 | 2.80 | 4.00 | 2.49 | 1.60 | 9.11 |
S.D. | 1.39 | 1.16 | 2.56 | 4.95 × 10−01 | 9.79 × 10−01 | 1.89 | ||
f5 | 1 | Mean | 7.23 × 10−04 | 1.02 × 10−03 | 3.25 × 10−03 | 1.49 × 10−03 | 8.80 × 10−04 | 2.15 × 10−04 |
S.D. | 3.90 × 10−04 | 5.46 × 10−04 | 2.43 × 10−03 | 7.09 × 10−04 | 4.49 × 10−04 | 1.58 × 10−04 | ||
f6 | 0 | Mean | 5.99 | 2.13 | 9.35 | 0.00 | 2.67 | 0.00 |
S.D. | 2.74 | 1.15 | 3.53 | 0.00 | 1.36 | 0.00 | ||
f7 | 0 | Mean | 4.14 × 10−15 | 4.07 × 10−15 | 4.67 × 10−15 | 3.16 × 10−14 | 4.00 × 10−15 | 3.91 × 10−15 |
S.D. | 7.03 × 10−16 | 5.02 × 10−16 | 1.62 × 10−15 | 5.29 × 10−14 | 0.00 | 5.62 × 10−16 | ||
f8 | 1 | Mean | 7.44 × 10−02 | 6.88 × 10−02 | 1.06 × 10−01 | 1.64 × 10−03 | 7.12 × 10−02 | 6.65 × 10−02 |
S.D. | 3.41 × 10−02 | 3.13 × 10−02 | 5.81 × 10−02 | 7.61 × 10−04 | 3.55 × 10−02 | 3.75 × 10−02 | ||
f9 | 1 | Mean | 1.24 × 10−02 | 4.71 × 10−32 | 6.22 × 10−03 | 2.99 × 10−30 | 1.38 × 10−28 | 1.78 × 10−01 |
S.D. | 6.16 × 10−02 | 1.11 × 10−47 | 4.40 × 10−02 | 3.07 × 10−30 | 2.65 × 10−28 | 6.11 × 10−02 | ||
f10 | 1 | Mean | 3.56 × 10−04 | 9.33 × 10−20 | 9.80 × 10−04 | 3.16 × 10−27 | 2.57 × 10−13 | 1.35 × 10−32 |
S.D. | 2.52 × 10−03 | 1.90 × 10−19 | 6.86 × 10−03 | 6.79 × 10−27 | 4.33 × 10−13 | 8.29 × 10−48 |
Function | k | Criteria | PSO-w | GPSO | LPSO | CLPSO | DMS-PSO | DCS-PSO |
f1 | 1 | Mean | 1.28 × 10−39 | 3.42 × 10−14 | 3.82 | 3.83 × 10−13 | 6.76 × 10−28 | 1.05 × 10−40 |
S.D. | 2.48 × 10−39 | 1.20 × 10−13 | 2.18× 10+01 | 9.48 × 10−13 | 5.49 × 10−28 | 2.01 × 10−40 | ||
Rank | 2 | 4 | 6 | 5 | 3 | 1 | ||
f2 | 1 | Mean | 1.67 × 10−09 | 1.27 × 10−08 | 1.83 | 4.45 × 10−08 | 3.74 × 10−13 | 2.49 × 10−23 |
S.D. | 1.15 × 10−08 | 2.46 × 10−08 | 1.83 | 2.84 × 10−08 | 2.33 × 10−13 | 5.70 × 10−23 | ||
Rank | 3 | 4 | 6 | 5 | 2 | 1 | ||
f3 | 1 | Mean | 6.50 × 10−01 | 9.58 × 10+01 | 4.79× 10+02 | 4.37 × 10−04 | 1.16 × 10+01 | 5.90 |
S.D. | 5.22 × 10−01 | 5.10 × 10+01 | 3.08× 10+02 | 3.04 × 10−04 | 5.95 | 4.45 | ||
Rank | 2 | 5 | 6 | 1 | 4 | 3 | ||
f4 | 1 | Mean | 2.98× 10+01 | 4.17 × 10+01 | 9.42× 10+01 | 1.81 × 10 +01 | 2.26 × 10+01 | 2.21 × 10+01 |
S.D. | 2.16× 10+01 | 2.87 × 10+01 | 5.46× 10+01 | 3.80 | 2.29 | 6.32 × 10−01 | ||
Rank | 4 | 5 | 6 | 1 | 3 | 2 | ||
f5 | 1 | Mean | 8.74 × 10−03 | 1.89 × 10−02 | 8.00 × 10−02 | 3.83 × 10−02 | 9.99 × 10−03 | 5.12 × 10−03 |
S.D. | 3.28 × 10−03 | 6.58 × 10−03 | 3.51 × 10−02 | 1.47 × 10−02 | 3.98 × 10−03 | 1.69 × 10−03 | ||
Rank | 2 | 4 | 6 | 5 | 3 | 1 | ||
f6 | 1 | Mean | 3.88 × 10+01 | 3.34 × 10+01 | 5.12 × 10+01 | 4.40 × 10−11 | 2.73 × 10+01 | 1.31 × 10+02 |
S.D. | 8.06 | 8.19 | 1.62 × 10+01 | 6.58 × 10−11 | 8.88 × 10+00 | 1.01 × 10+01 | ||
Rank | 4 | 3 | 5 | 1 | 2 | 6 | ||
f7 | 1 | Mean | 3.81 × 10−01 | 1.14 × 10−06 | 4.32 | 3.15 × 10−03 | 4.37 × 10−14 | 1.20 × 10−14 |
S.D. | 6.47 × 10−01 | 1.04 × 10−06 | 1.16 | 2.18 × 10−02 | 4.86 × 10−14 | 3.90 × 10−15 | ||
Rank | 5 | 3 | 6 | 4 | 2 | 1 | ||
f8 | 0 | Mean | 1.46 × 10−02 | 1.38 × 10−02 | 4.14 × 10−01 | 1.41 × 10−12 | 4.44 × 10−03 | 1.15 × 10−02 |
S.D. | 1.70 × 10−02 | 1.55 × 10−02 | 3.12 × 10−01 | 2.52 × 10−12 | 4.99 × 10−03 | 1.40 × 10−02 | ||
Rank | 5 | 4 | 6 | 1 | 2 | 3 | ||
f9 | 0 | Mean | 1.39 × 10−01 | 1.04 × 10−02 | 3.87 | 2.32 × 10−14 | 1.08 × 10−26 | 2.96 × 10−02 |
S.D. | 2.62 × 10−01 | 3.14 × 10−02 | 2.66 | 5.72 × 10−14 | 2.24 × 10−26 | 4.74 × 10−02 | ||
Rank | 4 | 3 | 6 | 2 | 1 | 4 | ||
f10 | 1 | Mean | 2.28 × 10−01 | 1.10 × 10−02 | 2.16× 10+01 | 1.03 × 10−11 | 5.31 × 10−03 | 6.22 × 10−32 |
S.D. | 6.63 × 10−01 | 4.43 × 10−02 | 1.12× 10+01 | 1.45 × 10−11 | 4.13 × 10−04 | 7.45 × 10−32 | ||
Rank | 5 | 4 | 6 | 2 | 3 | 1 | ||
Ave. rank | 3.6 | 3.9 | 5.9 | 2.7 | 2.5 | 2.3 | ||
Final rank | 4 | 5 | 6 | 3 | 2 | 1 |
Function | k | Criteria | PSO-w | GPSO | LPSO | CLPSO | DMS-PSO | DCS-PSO |
f1 | 1 | Mean | 5.20 × 10−17 | 3.68 × 10−04 | 1.30 × 10+03 | 5.70 × 10−05 | 1.26 × 10−06 | 3.93 × 10−18 |
S.D. | 1.91 × 10−16 | 4.22 × 10−04 | 4.91 × 10+02 | 6.50 × 10−05 | 3.14 × 10−06 | 1.05 × 10−17 | ||
f2 | 0 | Mean | 9.06 × 10−06 | 2.01 × 10−01 | 2.26 × 10+01 | 2.00 × 10−01 | 1.48 × 10−04 | 5.87 × 10−07 |
S.D. | 5.83 × 10−05 | 1.41 | 4.26 | 1.41 | 6.99 × 10−04 | 2.82 × 10−06 | ||
f3 | 1 | Mean | 4.52 × 10+02 | 4.38 × 10+03 | 7.71 × 10+03 | 1.55 | 1.90 × 10+03 | 1.33 × 10+03 |
S.D. | 1.68 × 10+02 | 1.88× 10+03 | 2.11 × 10+03 | 9.59 × 10−01 | 2.30 × 10+03 | 5.06× 10+02 | ||
f4 | 1 | Mean | 7.59 × 10+01 | 1.34× 10+02 | 3.19 × 10+03 | 1.69 × 10+02 | 1.46 × 10+02 | 4.52 × 10−02 |
S.D. | 4.33 × 10+01 | 1.52× 10+02 | 1.65 × 10+03 | 2.75 × 10+02 | 1.41 × 10+02 | 2.98 × 10−02 | ||
f5 | 1 | Mean | 3.75 × 10−02 | 7.50 × 10−02 | 9.04 × 10−01 | 1.72 × 10−01 | 4.52 × 10−02 | 1.85 × 10−02 |
S.D. | 1.12 × 10−02 | 2.09 × 10−02 | 3.67 × 10−01 | 4.36 × 10−02 | 1.67 × 10−02 | 6.22 × 10−03 | ||
f6 | 1 | Mean | 7.48 × 10+01 | 8.78× 10+01 | 1.16 × 10+02 | 1.32 × 10−03 | 7.67 × 10+01 | 4.79 × 10+01 |
S.D. | 1.62 × 10+01 | 1.83× 10+01 | 2.04 × 10+01 | 6.87 × 10−03 | 1.60 × 10+01 | 3.18 | ||
f7 | 0 | Mean | 1.38 | 1.40 × 10−01 | 9.00 | 1.22 | 1.14 | 9.05 × 10−01 |
S.D. | 8.89 × 10−01 | 3.52 × 10−01 | 1.08 | 5.83 × 10−01 | 7.77 × 10−01 | 8.47 × 10−01 | ||
f8 | 0 | Mean | 1.36 × 10−02 | 8.18 × 10−03 | 1.22 × 10+01 | 3.18 × 10−04 | 3.79 × 10−02 | 1.43 × 10−02 |
S.D. | 2.61 × 10−02 | 1.02 × 10−02 | 4.71 | 1.41 × 10−03 | 7.54 × 10−02 | 1.76 × 10−02 | ||
f9 | 0 | Mean | 1.60 × 10−01 | 1.04 × 10−01 | 2.08 × 10+01 | 8.95 × 10−02 | 7.84 × 10−08 | 2.99 × 10−02 |
S.D. | 2.65 × 10−01 | 1.96 × 10−01 | 8.57 | 1.44 × 10−01 | 4.22 × 10−07 | 4.91 × 10−02 | ||
f10 | 1 | Mean | 4.07 | 9.60 × 10−01 | 4.06 × 10+02 | 1.35 | 8.78 | 3.53 × 10−01 |
S.D. | 3.36 | 1.12 | 3.10 × 10+02 | 1.87 | 6.03 | 3.69 × 10−01 |
1) Results for the 10-D problems: For f1, the performance of DCS-PSO is second only to LPSO, which maintains population diversity based on different neighborhood topologies, so its performance varies greatly on problems with different topologies. For f2 and f3, our DCS-PSO very clearly outperforms all other algorithms, in solution accuracy. DCS-PSO achieves comparable result to other algorithms on f4 with CLPSO and DMS-PSO exceeding it. f5 is a noise function that is difficult to solve due to the perturbation factor. However, DCS-PSO still gets the best result on this problem.
For multimodal problems, comparing DCS-PSO with CLPSO and DMS-PSO, which are algorithms proposed to improve the performance of PSO on multimodal problems, can demonstrate the competitiveness of DCS-PSO. DCS-PSO and CLPSO perform equally well with zero error on Rastrigin's function (f6), which has many local optima. For f7, DCS-PSO achieves the best result, but according to the statistical test result, its performance is comparable to that of DMS-PSO. For f8, PSO-w, GPSO, DMS-PSO and DCS-PSO perform similarly, but CLPSO performs relatively better. For two penalized functions f9 and f10, although the performance of DCS-PSO on f9 is not as good as that of GPSO, CLPSO and DMS-PSO, DCS-PSO achieves the lowest error on f10.
Overall, DCS-PSO obtains the best results in solving three unimodal functions (f2, f3 and f5) and three multimodal functions (f6, f7 and f10). Additionally, Table 6 shows that the S.D. values of DCS-PSO on most test functions are much smaller than those of other algorithms, indicating that DCS-PSO is more stable during the convergence process.
2) Results for the 30-D and 50-D problems: The experiments conducted on 10-D problems are repeated for 30-D and 50-D problems and the results are presented in Tables 7 and 8. All algorithms exhibited similar characteristics as on the 10-D problems. However, high dimensionality leads to increased complexity, resulting in a decline in the performance of optimization algorithms on most problems. As shown in Tables 7 and 8, for f1, f2 and f5 with 30 and 50 dimensions, DCS-PSO consistently outperforms all other algorithms, particularly on f2. Notably, for the Rosenbrock function (f4), DCS-PSO performs similarly to other algorithms in 10 dimensionalities. However, its performance does not deteriorate as the dimension increases, illustrating that high dimensionality does not significantly decrease the performance of DCS-PSO on this problem.
For f6 with 30 and 50 dimensions, CLPSO always performs best among all algorithms. DCS-PSO is still best on f7 with 30 dimensionalities, while GPSO performs better on this problem with 50 dimensions. CLPSO and DMS-PSO retain their superiority on f8 and f9 with 30 and 50 dimensionalities, respectively. However, the nonparametric Wilcoxon rank-sum test results indicate that DCS-PSO is not statistically significantly poorer than CLPSO and DMS-PSO. DCS-PSO remains the best-performing algorithm on f10, which suggests that its performance is not affected by the increased dimensionality of this function.
3) Comparison of convergence speed: Figure 6 shows the convergence characteristics in terms of the error between the best fitness value found and the actual best function value of the median run of each algorithm for each test function with 30 dimensions. Since the convergence characteristics for the 10-D and 50-D cases are similar to those in the 30-D case, they are not presented.
Combining the numerical results and convergence plots, it can be concluded that all algorithms perform well in solving f1 and f2 with three different dimensionalities. However, DCS-PSO is faster than other algorithms. PSO exhibits similar convergence characteristics as DCS-PSO on f2. The convergence characteristics of each algorithm on f1 and f2 are shown in Figure 6(a) and (b). Figure 6(c) indicates that CLPSO has better convergence speed and accuracy than other algorithms on f3. Although DCS-PSO has slightly lower convergence accuracy than CLPSO on f4, it has the fastest convergence speed as illustrated in Figure 6(d). All algorithms perform similarly on the noise problem (f5). Figure 6(e) displays that these algorithms show a slow decrease in error, but DCS-PSO has the fastest convergence speed and the lowest error.
The advantages of CLPSO and DMS-PSO are revealed through some multimodal problems. As shown in Figure 6(f), although the convergence speed is slow, CLPSO has the best search accuracy on f6. In contrast, other algorithms cannot continue to search effectively after reaching a certain accuracy. From the convergence characteristics in Figure 6(g), the performances of DCS-PSO and DMS-PSO are also excellent on f7, but DCS-PSO has the fastest convergence speed. Although DCS-PSO is not as good as CLPSO and DMS-PSO in terms of average errors on f8 and f9, there is no statistically significant difference according to the Wilcoxon tests, as verified in Figure 6(h, i). It can be seen from Figure 6(j) that DCS-PSO has the best search efficiency on f10. Considering the numerical results and the convergence speeds, DCS-PSO performs best among all tested algorithms.
4) Discussion: The nonparametric Wilcoxon rank-sum tests indicate that DCS-PSO is significantly better than or equivalent to the other algorithms on most problems with three different dimensionality cases. When the problems are 10-dimensional, DCS-PSO is significantly better than other algorithms on f2, f3, f5 and f10, and is equivalent to CLPSO and DMS-PSO on f6 and f7, respectively. For f8 and f9, with 30 and 50 dimensionalities, DCS-PSO performs equivalently to CLPSO and DMS-PSO, respectively. All statistical results demonstrate that DCS-PSO performs very well overall.
In the above experiments, DCS-PSO shows better convergence speed and accuracy on most unimodal problems, indicating that it is more effective in solving such problems. Simultaneously,
the performance of DCS-PSO on most multimodal problems is comparable to that of CLPSO and DMS-PSO, demonstrating its ideal global search performance and ability to escape local optima.
This section further uses a widely recognized CEC2017 test suite to confirm the performance of DCS-PSO with PSO variants. In CEC2017, there are 29 functions (F2 has been excluded because it shows unstable behavior, especially for higher dimensions) which can be classified into four categories: unimodal functions (f1, f3), multimodal functions (f4 – f10), hybrid functions (f11 – f20) and composition functions (f21 – f30). Due to space constraints, we only select D = 30 for analysis and discussion. The best value for each function is shown in bold.
Table 9 shows that CLPSO and DMS-PSO obtain the best result in unimodal functions, followed by DCS-PSO. In multimodal functions, DCS-PSO achieves the best performance on f4, f5, f8 and f9, 4 of 7 multimodal functions in terms of mean value, followed by GPSO, CLPSO, and DMS-PSO, which all obtain the best solutions on one function.
Function | Criteria | PSO-w | GPSO | LPSO | CLPSO | DMS-PSO | DCS-PSO |
F1 | Mean | 5.73 × 10+08 | 9.15 × 10+09 | 5.24 × 10+09 | 1.01 × 10+02 | 3.26 × 10+03 | 6.44 × 10+03 |
S.D. | 2.36 × 10+09 | 1.29 × 10+10 | 9.65 × 10+09 | 5.12 × 10−01 | 4.53 × 10+03 | 1.93 × 10+03 | |
F3 | Mean | 2.52 × 10+03 | 5.01 × 10+03 | 9.58 × 10+03 | 4.90 × 10+02 | 2.16 × 10+02 | 5.45 × 10+02 |
S.D. | 1.23 × 10+03 | 2.38 × 10+03 | 5.11 × 10+03 | 1.39 × 10+01 | 3.43 × 10+02 | 3.90 × 10−13 | |
F4 | Mean | 1.03 × 10+02 | 1.73 × 10+02 | 2.15 × 10+02 | 5.53 × 10+02 | 7.39 × 10+01 | 1.22 |
S.D. | 3.19 × 10+01 | 1.02 × 10+02 | 1.92 × 10+02 | 7.54 | 2.39 × 10+01 | 1.75 | |
F5 | Mean | 1.43 × 10+02 | 7.61 × 10+01 | 1.42 × 10+02 | 6.00 × 10+02 | 5.88 × 10+01 | 2.09 × 10+01 |
S.D. | 3.11 × 10+01 | 2.08 × 10+01 | 3.78 × 10+01 | 0.00 | 1.30 × 10+01 | 7.01 | |
F6 | Mean | 4.01 × 10+01 | 1.43 × 10+01 | 5.14 × 10+01 | 7.87 × 10+02 | 9.86 × 10−04 | 3.75 |
S.D. | 1.49 × 10+01 | 5.47 | 1.39 × 10+01 | 6.43 | 2.47 × 10−03 | 3.17 | |
F7 | Mean | 1.29 × 10+02 | 1.12 × 10+01 | 1.91 × 10+02 | 8.49 × 10+02 | 9.76 × 10+01 | 1.28 × 10+02 |
S.D. | 3.46 × 10+01 | 1.03 × 10+01 | 4.86 × 10+01 | 1.03 × 10+01 | 1.99 × 10+01 | 9.23 | |
F8 | Mean | 1.20 × 10+02 | 6.83 × 10+01 | 1.24 × 10+02 | 9.09 × 10+02 | 5.94 × 10+01 | 1.62 × 10+01 |
S.D. | 3.18 × 10+01 | 1.66 × 10+01 | 2.95 × 10+01 | 3.05 | 1.68 × 10+01 | 5.52 | |
F9 | Mean | 1.38 × 10+03 | 1.53 × 10+02 | 2.26 × 10+03 | 3.27 × 10+03 | 2.33 × 10+01 | 6.09 × 10−01 |
S.D. | 1.18 × 10+03 | 1.64 × 10+02 | 9.36 × 10+02 | 2.89 × 10+02 | 5.63 × 10+01 | 1.83 | |
F10 | Mean | 3.61 × 10+03 | 3.19 × 10+03 | 3.91 × 10+03 | 1.18 × 10+03 | 2.97 × 10+03 | 1.99 × 10+03 |
S.D. | 5.67 × 10+02 | 6.28 × 10+02 | 7.81 × 10+02 | 4.51 | 5.89 × 10+02 | 2.87 × 10+02 | |
F11 | Mean | 1.30 × 10+02 | 1.71 × 10+02 | 1.76 × 10+02 | 5.97 × 10+05 | 4.74 × 10+01 | 1.12 × 10+02 |
S.D. | 4.18 × 10+01 | 5.41 × 10+01 | 5.11 × 10+01 | 3.96 × 10+05 | 3.12 × 10+01 | 2.44 × 10+01 | |
F12 | Mean | 4.56 × 10+07 | 6.40 × 10+07 | 1.76 × 10+08 | 2.62 × 10+03 | 8.25 × 10+04 | 2.79 × 10+03 |
S.D. | 2.83 × 10+08 | 1.12 × 10+08 | 5.00 × 10+08 | 1.55 × 10+03 | 2.16 × 10+04 | 2.07 × 10+03 | |
F13 | Mean | 1.77 × 10+06 | 5.01 × 10+07 | 1.88 × 10+07 | 1.75 × 10+04 | 1.55 × 10+04 | 1.27 × 10+04 |
S.D. | 8.67 × 10+06 | 1.71 × 10+08 | 1.02 × 10+08 | 9.65 × 10+03 | 1.65 × 10+04 | 3.86 × 10+03 | |
F14 | Mean | 1.80 × 10+04 | 3.70 × 10+04 | 1.10 × 10+04 | 1.61 × 10+03 | 1.32 × 10+04 | 1.02 × 10+02 |
S.D. | 2.36 × 10+04 | 4.40 × 10+04 | 1.30 × 10+04 | 2.70 × 10+01 | 1.23 × 10+04 | 2.45 × 10+01 | |
F15 | Mean | 3.99 × 10+03 | 1.97 × 10+04 | 1.07 × 10+04 | 2.16 × 10+03 | 6.90 × 10+03 | 1.09 × 10+02 |
S.D. | 3.96 × 10+03 | 3.05 × 10+04 | 1.25 × 10+04 | 1.32 × 10+02 | 9.16 × 10+03 | 1.14 × 10+02 | |
F16 | Mean | 1.04 × 10+03 | 7.61 × 10+02 | 1.13 × 10+03 | 1.84 × 10+03 | 7.95 × 10+02 | 1.79 × 10+03 |
S.D. | 2.73 × 10+02 | 2.47 × 10+02 | 2.95 × 10+02 | 4.18 × 10+01 | 2.10 × 10+02 | 1.25 × 10+02 | |
F17 | Mean | 4.00 × 10+04 | 1.57 × 10+04 | 1.21 × 10+04 | 2.16 × 10+05 | 1.89 × 10+02 | 4.62 × 10+03 |
S.D. | 6.54 × 10+02 | 1.54 × 10+03 | 2.56 × 10+03 | 1.61 × 10+05 | 1.14 × 10+02 | 6.54 × 10+02 | |
F18 | Mean | 1.41 × 10+05 | 4.45 × 10+05 | 1.12 × 10+05 | 1.95 × 10+03 | 1.27 × 10+05 | 5.56 × 10+02 |
S.D. | 1.04 × 10+05 | 3.62 × 10+05 | 1.15 × 10+05 | 1.88 × 10+01 | 1.06 × 10+05 | 8.59 × 10+02 | |
F19 | Mean | 6.20 × 10+03 | 4.47 × 10+05 | 6.10 × 10+03 | 6.30 × 10+03 | 9.39 × 10+03 | 7.78 × 10+03 |
S.D. | 6.09 × 10+03 | 2.71 × 10+06 | 6.37 × 10+03 | 2.64 × 10+01 | 1.28 × 10+04 | 1.08 × 10+03 | |
F20 | Mean | 1.70 × 10+04 | 1.45 × 10+04 | 1.92 × 10+04 | 2.36 × 10+03 | 2.73 × 10+02 | 6.00 × 10+03 |
S.D. | 3.26 × 10+03 | 3.44 × 10+02 | 2.75 × 10+03 | 5.04 | 1.41 × 10+02 | 5.32 × 10+02 | |
F21 | Mean | 3.39 × 10+02 | 2.80 × 10+02 | 3.34 × 10+02 | 2.65 × 10+03 | 2.62 × 10+02 | 1.03 × 10+02 |
S.D. | 5.45 × 10+01 | 4.25 × 10+01 | 4.90 × 10+01 | 5.40 × 10+02 | 2.11 × 10+01 | 8.18 × 10−01 | |
F22 | Mean | 4.15 × 10+02 | 6.72 × 10+02 | 1.00 × 10+03 | 2.72 × 10+03 | 2.38 × 10+02 | 1.09 × 10+02 |
S.D. | 9.94 × 10+02 | 7.48 × 10+02 | 1.50 × 10+03 | 1.15 × 10+01 | 5.92 × 10+02 | 1.44 | |
F23 | Mean | 6.42 × 10+02 | 5.82 × 10+02 | 6.73 × 10+02 | 2.92 × 10+03 | 4.30 × 10+02 | 3.25 × 10+02 |
S.D. | 1.01 × 10+02 | 9.36 × 10+01 | 1.10 × 10+02 | 1.21 × 10+01 | 2.10 × 10+01 | 7.91 | |
F24 | Mean | 6.61 × 10+02 | 6.34 × 10+02 | 7.14 × 10+02 | 2.89 × 10+03 | 4.94 × 10+02 | 3.53 × 10+02 |
S.D. | 7.19 × 10+01 | 8.33 × 10+01 | 8.07 × 10+01 | 2.23 × 10−01 | 2.60 × 10+01 | 1.09 × 10+01 | |
F25 | Mean | 3.95 × 10+02 | 4.13 × 10+02 | 4.47 × 10+02 | 3.47 × 10+03 | 3.97 × 10+02 | 4.26 × 10+02 |
S.D. | 1.50 × 10+01 | 3.64 × 10+01 | 3.05 × 10+01 | 7.89 × 10+02 | 1.53 | 2.34 × 10+01 | |
F26 | Mean | 1.24 × 10+03 | 1.51 × 10+03 | 2.45 × 10+03 | 3.21 × 10+03 | 1.49 × 10+03 | 4.02 × 10+02 |
S.D. | 1.29 × 10+03 | 7.71 × 10+02 | 1.36 × 10+03 | 6.70 | 5.18 × 10+02 | 5.94 × 10+01 | |
F27 | Mean | 5.88 × 10+02 | 5.80 × 10+02 | 6.21 × 10+02 | 3.21 × 10+03 | 5.38 × 10+02 | 4.03 × 10+02 |
S.D. | 5.87 × 10+01 | 5.54 × 10+01 | 5.36 × 10+01 | 5.65 | 1.25 × 10+01 | 4.45 | |
F28 | Mean | 4.34 × 10+02 | 5.05 × 10+02 | 5.49 × 10+02 | 3.45 × 10+03 | 3.40 × 10+02 | 3.88 × 10+02 |
S.D. | 4.76 × 10+01 | 7.23 × 10+01 | 1.46 × 10+02 | 4.85 | 5.79 × 10+01 | 8.36 | |
F29 | Mean | 6.25 × 10+02 | 8.71 × 10+02 | 6.24 × 10+02 | 8.91 × 10+02 | 6.03 × 10+02 | 6.38 × 10+02 |
S.D. | 4.87 × 10+01 | 6.31 × 10+01 | 2.17 × 10+02 | 7.46 × 10+01 | 1.55 × 10+02 | 7.88 × 10+01 | |
F30 | Mean | 4.52 × 10+04 | 4.50 × 10+05 | 3.78 × 10+05 | 1.12 × 10+03 | 5.81 × 10+03 | 3.14 × 10+05 |
S.D. | 1.10 × 10+05 | 1.41 × 10+06 | 4.93 × 10+05 | 2.22 × 10+03 | 2.34 × 10+03 | 5.14 × 10+05 |
For the 10 hybrid functions, DCS-PSO reaches the best solutions on f13, f14, f15 and f18, 4 of 10 hybrid functions. DMS-PSO also performs remarkably since it reaches the best solutions on 3 out of 10 functions, followed by GPSO, LPSO and CLPSO, which all obtain the best solutions on one function. In the ten composition functions, DMS-PSO attains remarkable results on f21, f22, f23, f24, f26 and f27, 6 out of 10 composition functions. The comparison results indicate that DCS-PSO works well in complex and multimodal situations.
Friedman test is used in this section to compare the performance of DCS-PSO and other PSO variants. The results of the Friedman test are shown in Table 10. It can be seen from Table 10 that DCS-PSO achieves a remarkable result on the CEC2017, followed by DMS-PSO, PSO-w, GPSO, CLPSO and LPSO. The overall performance in CEC2017 reveals that DCS-PSO has a distinct advantage over other compared algorithms.
Average rank | Algorithm | Ranking |
1 | DCS-PSO | 2.03 |
2 | DMS-PSO | 2.17 |
3 | PSO-w | 3.69 |
4 | GPSO | 4.07 |
5 | CLPSO | 4.45 |
6 | LPSO | 4.59 |
In this paper, we propose an improved particle swarm optimization combined with double-chaos search mechanism, namely DCS-PSO. The standard PSO is simple and efficient, but it is difficult to maintain population diversity when dealing with complex problems with a large search space, leading to premature convergence and getting trapped in local optima. The chaos search mechanism has the advantages of global traversal and avoiding falling into local optima.
In the first stage of DCS-PSO, we adopt the double-chaos search mechanism to narrow the search space to the vicinity of the optimal solution, effectively reducing the risk of PSO getting trapped in local optima. In the second stage, to enhance the population diversity, the logistic map is utilized to perform a global search in the narrowed search space and the best solution obtained from both chaos search and population search will guide the population to converge. Experimental studies show that the algorithm can effectively narrow the search space and has better convergence accuracy and speed for most functions.
Of course, according to the no free lunch theorem, DCS-PSO cannot optimally solve every kind of global optimization problem. Although we use the chaos search mechanism to improve the convergence speed and ability to escape the local optima of DCS-PSO, it also brings about a problem worth further investigation: how to ensure that the optimal solution must exist in the narrowed search space. In the future, we will continue conducting in-depth research on chaos theory, optimize our DCS-PSO algorithm and try to apply DCS-PSO to solving these real-world problems.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
This research was supported by the National Natural Science Foundation of China (Grant No 82260849); the National Natural Science Foundation of China (Grant No 61562045); and Jiangxi University of Chinese Medicine Science and Technology Innovation Team Development Program (Grant No CXTD22015).
[1] | (1997) World Health OrganizationSilica, some silicates, coal dust and para-aramid fibrils. IARC Monographs on the Evaluation of Carcinogenic Risks to Humans. Lyon: IARC Press. |
[2] | Health and Safety Executive Research Report RR878 - Levels of respirable dust and respirable crystalline silica at construction sites.Available from: https://www.hse.gov.uk/research/rrhtm/rr878.htm. |
[3] | World Health Organization (2000) Crystalline Silica, Quartz (Concise International Chemical Assessment Document 24), ISBN-10: 9241530243.Available from: http://www.inchem.org/documents/cicads/cicads/cicad24.htm. |
[4] |
Pollard KM (2016) Silica, Silicosis, and Autoimmunity. Front Immunol 7: 97. doi: 10.3389/fimmu.2016.00097
![]() |
[5] |
Miller FW, Alfredsson L, Costenbader KH, et al. (2012) Epidemiology of environmental exposures and human autoimmune diseases: findings from a National Institute of Environmental Health Sciences Expert Panel Workshop. J Autoimmun 39: 259-271. doi: 10.1016/j.jaut.2012.05.002
![]() |
[6] |
Bartůňková J, Pelclova D, Fenclová Z, et al. (2006) Exposure to silica and risk of ANCA-associated vasculitis. Am J Ind Med 49: 569-576. doi: 10.1002/ajim.20327
![]() |
[7] |
Pollard KM, Hultman P, Kono DH (2010) Toxicology of autoimmune diseases. Chem Res Toxicol 23: 455-466. doi: 10.1021/tx9003787
![]() |
[8] |
Lee S, Matsuzaki H, Kumagai-Takei N, et al. (2014) Silica exposure and altered regulation of autoimmunity. Environ Health Prev Med 19: 322-329. doi: 10.1007/s12199-014-0403-9
![]() |
[9] |
Van der Woude FJ, Lobatto S, Permin H, et al. (1985) Autoantibodies against neutrophils and monocytes. Tool for diagnosis and marker of disease activity in Wegener's granulomatosis. Lancet 325: 425-429. doi: 10.1016/S0140-6736(85)91147-X
![]() |
[10] |
Germolec D, Kono DH, Pfau JC, et al. (2012) Animal models used to examine the role of the environment in the development of autoimmune disease: findings from an NIEHS Expert Panel Workshop. J Autoimmun 39: 285-293. doi: 10.1016/j.jaut.2012.05.020
![]() |
[11] |
Riediker M, Zink D, Kreyling W, et al. (2019) Particle toxicology and health-where are we? Part Fibre Toxicol 16: 19. doi: 10.1186/s12989-019-0302-8
![]() |
[12] |
Murphy CJ, Vartanian AM, Geiger FM, et al. (2015) Biological responses to engineered nanomaterials: Needs for the next decade. ACS Central Sci 1: 117-123. doi: 10.1021/acscentsci.5b00182
![]() |
[13] |
Fröhlich E (2015) Value of phagocyte function screening for immunotoxicity of nanoparticles in vivo. Int J Nanomed 10: 3761-3778. doi: 10.2147/IJN.S83068
![]() |
[14] |
Jochums A, Friehs E, Sambale F, et al. (2017) Revelation of Different Nanoparticle-Uptake Behavior in Two Standard Cell Lines NIH/3T3 and A549 by Flow Cytometry and Time-Lapse Imaging. Toxics 5: 15. doi: 10.3390/toxics5030015
![]() |
[15] |
Parks CG, Miller FW, Pollard KM, et al. (2014) Expert panel workshop consensus statement on the role of the environment in the development of autoimmune disease. Int J Mol Sci 15: 14269-14297. doi: 10.3390/ijms150814269
![]() |
[16] |
Baumann D, Hofmann D, Nullmeier S, et al. (2013) Complex encounters: nanoparticles in whole blood and their uptake into different types of white blood cells. Nanomedicine 8: 699-713. doi: 10.2217/nnm.12.111
![]() |
[17] |
Vranic S, Boggetto N, Contremoulins V, et al. (2013) Deciphering the mechanisms of cellular uptake of engineered nanoparticles by accurate evaluation of internalization using imaging flow cytometry. Part Fibre Toxicol 10: 2. doi: 10.1186/1743-8977-10-2
![]() |
[18] | Phanse Y, Ramer-Tait AE, Friend SL, et al. (2012) Analyzing cellular internalization of nanoparticles and bacteria by multi-spectral imaging flow cytometry. J Vis Exp 64: e3884. |
[19] |
Rieger AM, Hall BE, Barreda DR (2010) Macrophage activation differentially modulates particle binding, phagocytosis and downstream antimicrobial mechanisms. Dev Comp Immunol 34: 1144-1159. doi: 10.1016/j.dci.2010.06.006
![]() |
[20] | Marangon I, Boggetto N, Ménard-Moyon C, et al. (2013) Localization and relative quantification of carbon nanotubes in cells with multispectral imaging flow cytometry. J Vis Exp 82: e50566. |
[21] |
Smirnov A, Solga MD, Lannigan J, et al. (2015) An improved method for differentiating cell-bound from internalized particles by imaging flow cytometry. J Immunol Methods 423: 60-69. doi: 10.1016/j.jim.2015.04.028
![]() |
[22] |
Tian L, Dai S, Wang J, et al. (2008) Nanoquartz in Late Permian C1 coal and the high incidence of female lung cancer in the Pearl River Origin area: a retrospective cohort study. BMC Public Health 8: 398. doi: 10.1186/1471-2458-8-398
![]() |
[23] |
Hornung V, Bauernfeind F, Halle A, et al. (2008) Silica crystals and aluminum salts activate the NALP3 inflammasome through phagosomal destabilization. Nat Immunol 9: 847-856. doi: 10.1038/ni.1631
![]() |
[24] |
Beamer CA, Holian A (2005) Scavenger receptor class A type I/II (CD204) null mice fail to develop fibrosis following silica exposure. Am J Physiol Lung Cell Mol Physiol 289: L186-L195. doi: 10.1152/ajplung.00474.2004
![]() |
[25] |
Beamer GL, Seaver BP, Jessop F, et al. (2016) Acute exposure to crystalline silica reduces macrophage activation in response to bacterial lipoproteins. Front Immunol 7: 49. doi: 10.3389/fimmu.2016.00049
![]() |
[26] | Schorn C, Janko C, Latzko M, et al. (2012) Monosodium urate crystals induce extracellular DNA traps in neutrophils, eosinophils, and basophils but not in mononuclear cells. Front Immunol 3: 277. |
[27] |
Satpathy SR, Jala VR, Bodduluri SR, et al. (2015) Crystalline silica-induced leukotriene B4-dependent inflammation promotes lung tumour growth. Nat Commun 6: 7064. doi: 10.1038/ncomms8064
![]() |
[28] |
Chubb LG, Cauda EG (2017) Characterizing particle size distributions of crystalline silica in gold mine dust. Aerosol Air Qual Res 17: 24-33. doi: 10.4209/aaqr.2016.05.0179
![]() |
[29] |
Dominical V, Samsel L, McCoy JP (2017) Masks in imaging flow cytometry. Methods 112: 9-17. doi: 10.1016/j.ymeth.2016.07.013
![]() |
[30] |
Boyette LB, Macedo C, Hadi K, et al. (2017) Phenotype, function, and differentiation potential of human monocyte subsets. PLoS One 12: e0176460. doi: 10.1371/journal.pone.0176460
![]() |
[31] |
Rosales C (2018) Neutrophil: A cell with many roles in inflammation or several cell types? Front Physiol 9: 113. doi: 10.3389/fphys.2018.00113
![]() |
[32] |
Hewitt RE, Vis B, Pele LC, et al. (2017) Imaging flow cytometry assays for quantifying pigment grade titanium dioxide particle internalization and interactions with immune cells in whole blood. Cytom Part A 91: 1009-1020. doi: 10.1002/cyto.a.23245
![]() |
[33] |
Sellamuthu R, Umbright C, Roberts JR, et al. (2011) Blood gene expression profiling detects silica exposure and toxicity. Toxicol Sci 122: 253-264. doi: 10.1093/toxsci/kfr125
![]() |
[34] |
Ishihara Y, Yasuhara T, Ishiyama S, et al. (2001) The role of leukocytes during acute phase inflammation in crystalline silica-induced lung injury. Exp Lung Res 27: 589-603. doi: 10.1080/019021401753181845
![]() |
[35] |
Baldwin PEJ, Yates T, Beattie H, et al. (2019) Exposure to Respirable Crystalline Silica in the GB Brick Manufacturing and Stone Working Industries. Ann Work Expos Health 63: 184-196. doi: 10.1093/annweh/wxy103
![]() |
[36] | Künzli N, Perez L, Rapp R (2010) Air quality and health. Lausanne: European Respiratory Society. World Health Organ 89. |
[37] |
Schmid O, Cassee FR (2017) On the pivotal role of dose for particle toxicology and risk assessment: exposure is a poor surrogate for delivered dose. Part Fibre Toxicol 14: 52. doi: 10.1186/s12989-017-0233-1
![]() |
[38] |
Fröhlich E, Mercuri A, Wu S, et al. (2016) Measurements of deposition, lung surface area and lung fluid for simulation of inhaled compounds. Front Pharmacol 7: 181. doi: 10.3389/fphar.2016.00181
![]() |
[39] |
Friedetzky A, Garn H, Kirchner A, et al. (1998) Histopathological changes in enlarged thoracic lymph nodes during the development of silicosis in rats. Immunobiology 199: 119-132. doi: 10.1016/S0171-2985(98)80068-5
![]() |
[40] |
Xu M, Qing M, Peng D (2014) Silicon dioxide particles deposited in vessels and cartilage of the femoral head. Yonsei Med J 55: 1447-1449. doi: 10.3349/ymj.2014.55.5.1447
![]() |
[41] |
Slavin RE, Swedo JL, Brandes D, et al. (1985) Extrapulmonary silicosis: a clinical, morphologic, and ultrastructural study. Hum Pathol 16: 393-412. doi: 10.1016/S0046-8177(85)80233-1
![]() |
[42] |
Chu Z, Huang Y, Li L, et al. (2012) Physiological pathway of human cell damage induced by genotoxic crystalline silica nanoparticles. Biomaterials 33: 7540-7546. doi: 10.1016/j.biomaterials.2012.06.073
![]() |
[43] | Ueki A, Yamaguchi M, Ueki H, et al. (1994) Polyclonal human T-cell activation by silicate in vitro. Immunology 82: 332. |
[44] |
Wu P, Hyodoh F, Hatayama T, et al. (2005) Induction of CD69 antigen expression in peripheral blood mononuclear cells on exposure to silica, but not by asbestos/chrysotile-A. Immunol Lett 98: 145-152. doi: 10.1016/j.imlet.2004.11.005
![]() |
[45] |
Braakhuis HM, Park MVDZ, Gosens I, et al. (2014) Physicochemical characteristics of nanomaterials that affect pulmonary inflammation. Part Fibre Toxicol 11: 18. doi: 10.1186/1743-8977-11-18
![]() |
[46] |
Vis B, Hewitt RE, Faria N, et al. (2018) Non-functionalized ultrasmall silica nanoparticles directly and size-selectively activate T cells. ACS Nano 12: 10843-10854. doi: 10.1021/acsnano.8b03363
![]() |
[47] |
Vis B, Hewitt RE, Monie TP, et al. (2020) Ultrasmall silica nanoparticles directly ligate the T cell receptor complex. P Natl Acad Sci 117: 285-291. doi: 10.1073/pnas.1911360117
![]() |
[48] |
Eleftheriadis T, Pissas G, Zarogiannis S, et al. (2019) Crystalline silica activates the T-cell and the B-cell antigen receptor complexes and induces T-cell and B-cell proliferation. Autoimmunity 52: 136-143. doi: 10.1080/08916934.2019.1614171
![]() |
[49] |
Hewitt RE, Chappell HF, Powell JJ (2020) Small and dangerous? Potential toxicity mechanisms of common exposure particles and nanoparticles. Curr Opin Toxicol 19: 93-98. doi: 10.1016/j.cotox.2020.01.006
![]() |
[50] |
Gilberti RM, Joshi GN, Knecht DA (2008) The phagocytosis of crystalline silica particles by macrophages. Am J Resp Cell Mol Biol 39: 619-627. doi: 10.1165/rcmb.2008-0046OC
![]() |
[51] |
Pacheco P, White D, Sulchek T (2013) Effects of microparticle size and Fc density on macrophage phagocytosis. PLoS One 8: e60989. doi: 10.1371/journal.pone.0060989
![]() |
[52] | Dale DC, Boxer L, Liles WC (2008) The phagocytes: neutrophils and monocytes. Blood, J Am Soc Hematol 112: 935-945. |
[53] |
Brinkmann V, Zychlinsky A (2012) Neutrophil extracellular traps: is immunity the second function of chromatin? J Cell Biol 198: 773-783. doi: 10.1083/jcb.201203170
![]() |
[54] |
Desai J, Foresto-Neto O, Honarpisheh M, et al. (2017) Particles of different sizes and shapes induce neutrophil necroptosis followed by the release of neutrophil extracellular trap-like chromatin. Sci Rep 7: 15003. doi: 10.1038/s41598-017-15106-0
![]() |
![]() |
![]() |
1. | Santuan Qin, Huadie Zeng, Wei Sun, Jin Wu, Junhua Yang, Multi-Strategy Improved Particle Swarm Optimization Algorithm and Gazelle Optimization Algorithm and Application, 2024, 13, 2079-9292, 1580, 10.3390/electronics13081580 | |
2. | Guangyun Cui, Zhen Qi, Huaqing Zhao, Ranhang Zhao, Haofang Wang, Jiaxing Zhao, Application of F-HGAPSO Algorithm in Reservoir Flood Control Optimal Operation, 2024, 0920-4741, 10.1007/s11269-024-04045-x |
PC configuration | Information |
CPU | Inter(R) Core(TM) i7-8700 |
Frequency | 3.2GHz CPU |
RAM | 16.0GB |
Operating system | Windows 10(64Bit) |
Language | Python 3.9.13 |
Algorithm | Parameters |
COA | max_iter1=105, max_iter2=1.5×105, α=0.3 |
DCOA | h=3000, γ=0.25, ξ = 1.5 |
PCOA | n_group=3, max_iter1=105, max_iter2=1.5×105, α=0.3 |
DCS-PSO | N=20, w=0.4, c1=c2=2, Vmax=0.2×Range, max_iter=1000 |
Function | Criteria | COA | DCOA | PCOA | DCS-PSO |
fJ | E.T. | 30.8898 | 0.5753 | 17.2234 | 0.2793 |
Mean | 0.0000 | 0.0005 | 0.0000 | 0.0000 | |
S.D. | 0.0000 | 0.0014 | 0.0000 | 0.0000 | |
fC | E.T. | 7.6564 | 0.6478 | 10.3990 | 0.3610 |
Mean | −1.0316 | −1.0312 | −1.0316 | −1.0316 | |
S.D. | 0.0000 | 0.0012 | 0.0000 | 0.0000 | |
fG | E.T. | 14.1867 | 1.0656 | 14.0814 | 0.6097 |
Mean | 3.0022 | 3.0100 | 3.0024 | 3.0000 | |
S.D. | 0.0022 | 0.0126 | 0.0022 | 0.0000 | |
fB | E.T. | 1.5663 | 0.8653 | 5.7455 | 0.4513 |
Mean | 0.3983 | 0.4003 | 0.3980 | 0.3979 | |
S.D. | 0.0003 | 0.0039 | 0.0001 | 0.0000 | |
fR | E.T. | 4.4150 | 0.6652 | 5.8342 | 0.4366 |
Mean | −1.9993 | −1.9652 | −1.9992 | −2.0000 | |
S.D. | 0.0006 | 0.0558 | 0.0008 | 0.0000 | |
fS | E.T. | 5.8609 | 2.1743 | 13.8397 | 2.0615 |
Mean | −186.6428 | -186.4517 | −186.6861 | −186.7309 | |
S.D. | 0.0711 | 0.4203 | 0.0680 | 0.0000 |
Name | Functions | Range | Optimum |
Sphere | f1(x)=D∑i=1x2i | [−100, 100] | 0 |
Schwefel's 2.22 | f2(x)=D∑i=1|xi|+D∏i=1|xi| | [−10, 10] | 0 |
Schwefel's 1.2 | f3(x)=D∑i=1(i∑j=1xj)2 | [−100, 100] | 0 |
Rosenbrock | f4(x)=D−1∑i=1[100(xi+1−x2i)2+(xi−1)2] | [−10, 10] | 0 |
Noise | f5(x)=D∑i=1ix4i+random[0,1) | [−1.28, 1.28] | 0 |
Rastrigin | f6(x)=D∑i=1(x2i−10cos(2πxi)+10) | [−5.12, 5.12] | 0 |
Ackley | f7(x)=−20exp(−0.2√1DD∑i=1x2i)−exp(1DD∑i=1cos(2πxi))+20+e | [−32, 32] | 0 |
Griewank | f8(x)=14000D∑i=1x2i−D∏i=1cos(xi√i)+1 | [−600, 600] | 0 |
Penalized 1 | f9(x)=πD{10sin2(πy1)+D−1∑i=1(yi−1)2[1+10sin2(πyi+1)]+(yD−1)2}+D∑i=1u(xi,10,100,4) | [−50, 50] | 0 |
Penalized 2 | f10(x)=110{sin2(3πx1)+D∑i=1(xi−1)2[1+sin2(3πxi+1)]+(xD−1)2[1+sin2(2πxD)]}+D∑i=1u(xi,5,100,4) | [−50, 50] | 0 |
Algorithm | Parameters |
PSO-w | w:0.4, c1=c2=2.0, Vmax=0.2×Range |
GPSO | w:0.9∼0.4, c1=c2=2.0, Vmax=0.2×Range |
LPSO | χ=0.7298, c1=c2=1.49445, Vmax=0.2×Range |
CLPSO | w:0.9∼0.4, c=1.49445, Vmax=0.2×Range, m=7 |
DMS-PSO | χ=0.7298, c1=c2=1.49445, Vmax=0.2×Range, M=4, R=10 |
DCS-PSO | w:0.4,c1=c2=2.0,Vmax=0.2×Range, h=3000,γ=0.15,ξ=1.5 |
Function | k | Criteria | PSO-w | GPSO | LPSO | CLPSO | DMS-PSO | DCS-PSO |
f1 | 1 | Mean | 1.84 × 10−144 | 1.47 × 10−73 | 2.78 × 10−169 | 2.18 × 10−21 | 6.30 × 10−62 | 2.35 × 10−151 |
S.D. | 9.13 × 10−144 | 1.04 × 10−72 | 0.00 | 5.44 × 10−21 | 4.25 × 10−61 | 7.75 × 10−151 | ||
f2 | 1 | Mean | 1.89 × 10−81 | 6.24 × 10−38 | 7.59 × 10−09 | 4.30 × 10−14 | 3.45 × 10−20 | 9.65 × 10−84 |
S.D. | 1.26 × 10−80 | 2.99 × 10−37 | 3.80 × 10−08 | 1.39 × 10−13 | 5.86 × 10−20 | 4.81 × 10−83 | ||
f3 | 1 | Mean | 1.22 × 10−44 | 9.93 × 10−17 | 6.11 × 10−13 | 2.18 × 10−20 | 6.13 × 10−10 | 1.21 × 10−48 |
S.D. | 7.42 × 10−44 | 3.36 × 10−16 | 4.32 × 10−12 | 3.25 × 10−20 | 3.33 × 10−09 | 2.38 × 10−48 | ||
f4 | 1 | Mean | 2.79 | 2.80 | 4.00 | 2.49 | 1.60 | 9.11 |
S.D. | 1.39 | 1.16 | 2.56 | 4.95 × 10−01 | 9.79 × 10−01 | 1.89 | ||
f5 | 1 | Mean | 7.23 × 10−04 | 1.02 × 10−03 | 3.25 × 10−03 | 1.49 × 10−03 | 8.80 × 10−04 | 2.15 × 10−04 |
S.D. | 3.90 × 10−04 | 5.46 × 10−04 | 2.43 × 10−03 | 7.09 × 10−04 | 4.49 × 10−04 | 1.58 × 10−04 | ||
f6 | 0 | Mean | 5.99 | 2.13 | 9.35 | 0.00 | 2.67 | 0.00 |
S.D. | 2.74 | 1.15 | 3.53 | 0.00 | 1.36 | 0.00 | ||
f7 | 0 | Mean | 4.14 × 10−15 | 4.07 × 10−15 | 4.67 × 10−15 | 3.16 × 10−14 | 4.00 × 10−15 | 3.91 × 10−15 |
S.D. | 7.03 × 10−16 | 5.02 × 10−16 | 1.62 × 10−15 | 5.29 × 10−14 | 0.00 | 5.62 × 10−16 | ||
f8 | 1 | Mean | 7.44 × 10−02 | 6.88 × 10−02 | 1.06 × 10−01 | 1.64 × 10−03 | 7.12 × 10−02 | 6.65 × 10−02 |
S.D. | 3.41 × 10−02 | 3.13 × 10−02 | 5.81 × 10−02 | 7.61 × 10−04 | 3.55 × 10−02 | 3.75 × 10−02 | ||
f9 | 1 | Mean | 1.24 × 10−02 | 4.71 × 10−32 | 6.22 × 10−03 | 2.99 × 10−30 | 1.38 × 10−28 | 1.78 × 10−01 |
S.D. | 6.16 × 10−02 | 1.11 × 10−47 | 4.40 × 10−02 | 3.07 × 10−30 | 2.65 × 10−28 | 6.11 × 10−02 | ||
f10 | 1 | Mean | 3.56 × 10−04 | 9.33 × 10−20 | 9.80 × 10−04 | 3.16 × 10−27 | 2.57 × 10−13 | 1.35 × 10−32 |
S.D. | 2.52 × 10−03 | 1.90 × 10−19 | 6.86 × 10−03 | 6.79 × 10−27 | 4.33 × 10−13 | 8.29 × 10−48 |
Function | k | Criteria | PSO-w | GPSO | LPSO | CLPSO | DMS-PSO | DCS-PSO |
f1 | 1 | Mean | 1.28 × 10−39 | 3.42 × 10−14 | 3.82 | 3.83 × 10−13 | 6.76 × 10−28 | 1.05 × 10−40 |
S.D. | 2.48 × 10−39 | 1.20 × 10−13 | 2.18× 10+01 | 9.48 × 10−13 | 5.49 × 10−28 | 2.01 × 10−40 | ||
Rank | 2 | 4 | 6 | 5 | 3 | 1 | ||
f2 | 1 | Mean | 1.67 × 10−09 | 1.27 × 10−08 | 1.83 | 4.45 × 10−08 | 3.74 × 10−13 | 2.49 × 10−23 |
S.D. | 1.15 × 10−08 | 2.46 × 10−08 | 1.83 | 2.84 × 10−08 | 2.33 × 10−13 | 5.70 × 10−23 | ||
Rank | 3 | 4 | 6 | 5 | 2 | 1 | ||
f3 | 1 | Mean | 6.50 × 10−01 | 9.58 × 10+01 | 4.79× 10+02 | 4.37 × 10−04 | 1.16 × 10+01 | 5.90 |
S.D. | 5.22 × 10−01 | 5.10 × 10+01 | 3.08× 10+02 | 3.04 × 10−04 | 5.95 | 4.45 | ||
Rank | 2 | 5 | 6 | 1 | 4 | 3 | ||
f4 | 1 | Mean | 2.98× 10+01 | 4.17 × 10+01 | 9.42× 10+01 | 1.81 × 10 +01 | 2.26 × 10+01 | 2.21 × 10+01 |
S.D. | 2.16× 10+01 | 2.87 × 10+01 | 5.46× 10+01 | 3.80 | 2.29 | 6.32 × 10−01 | ||
Rank | 4 | 5 | 6 | 1 | 3 | 2 | ||
f5 | 1 | Mean | 8.74 × 10−03 | 1.89 × 10−02 | 8.00 × 10−02 | 3.83 × 10−02 | 9.99 × 10−03 | 5.12 × 10−03 |
S.D. | 3.28 × 10−03 | 6.58 × 10−03 | 3.51 × 10−02 | 1.47 × 10−02 | 3.98 × 10−03 | 1.69 × 10−03 | ||
Rank | 2 | 4 | 6 | 5 | 3 | 1 | ||
f6 | 1 | Mean | 3.88 × 10+01 | 3.34 × 10+01 | 5.12 × 10+01 | 4.40 × 10−11 | 2.73 × 10+01 | 1.31 × 10+02 |
S.D. | 8.06 | 8.19 | 1.62 × 10+01 | 6.58 × 10−11 | 8.88 × 10+00 | 1.01 × 10+01 | ||
Rank | 4 | 3 | 5 | 1 | 2 | 6 | ||
f7 | 1 | Mean | 3.81 × 10−01 | 1.14 × 10−06 | 4.32 | 3.15 × 10−03 | 4.37 × 10−14 | 1.20 × 10−14 |
S.D. | 6.47 × 10−01 | 1.04 × 10−06 | 1.16 | 2.18 × 10−02 | 4.86 × 10−14 | 3.90 × 10−15 | ||
Rank | 5 | 3 | 6 | 4 | 2 | 1 | ||
f8 | 0 | Mean | 1.46 × 10−02 | 1.38 × 10−02 | 4.14 × 10−01 | 1.41 × 10−12 | 4.44 × 10−03 | 1.15 × 10−02 |
S.D. | 1.70 × 10−02 | 1.55 × 10−02 | 3.12 × 10−01 | 2.52 × 10−12 | 4.99 × 10−03 | 1.40 × 10−02 | ||
Rank | 5 | 4 | 6 | 1 | 2 | 3 | ||
f9 | 0 | Mean | 1.39 × 10−01 | 1.04 × 10−02 | 3.87 | 2.32 × 10−14 | 1.08 × 10−26 | 2.96 × 10−02 |
S.D. | 2.62 × 10−01 | 3.14 × 10−02 | 2.66 | 5.72 × 10−14 | 2.24 × 10−26 | 4.74 × 10−02 | ||
Rank | 4 | 3 | 6 | 2 | 1 | 4 | ||
f10 | 1 | Mean | 2.28 × 10−01 | 1.10 × 10−02 | 2.16× 10+01 | 1.03 × 10−11 | 5.31 × 10−03 | 6.22 × 10−32 |
S.D. | 6.63 × 10−01 | 4.43 × 10−02 | 1.12× 10+01 | 1.45 × 10−11 | 4.13 × 10−04 | 7.45 × 10−32 | ||
Rank | 5 | 4 | 6 | 2 | 3 | 1 | ||
Ave. rank | 3.6 | 3.9 | 5.9 | 2.7 | 2.5 | 2.3 | ||
Final rank | 4 | 5 | 6 | 3 | 2 | 1 |
Function | k | Criteria | PSO-w | GPSO | LPSO | CLPSO | DMS-PSO | DCS-PSO |
f1 | 1 | Mean | 5.20 × 10−17 | 3.68 × 10−04 | 1.30 × 10+03 | 5.70 × 10−05 | 1.26 × 10−06 | 3.93 × 10−18 |
S.D. | 1.91 × 10−16 | 4.22 × 10−04 | 4.91 × 10+02 | 6.50 × 10−05 | 3.14 × 10−06 | 1.05 × 10−17 | ||
f2 | 0 | Mean | 9.06 × 10−06 | 2.01 × 10−01 | 2.26 × 10+01 | 2.00 × 10−01 | 1.48 × 10−04 | 5.87 × 10−07 |
S.D. | 5.83 × 10−05 | 1.41 | 4.26 | 1.41 | 6.99 × 10−04 | 2.82 × 10−06 | ||
f3 | 1 | Mean | 4.52 × 10+02 | 4.38 × 10+03 | 7.71 × 10+03 | 1.55 | 1.90 × 10+03 | 1.33 × 10+03 |
S.D. | 1.68 × 10+02 | 1.88× 10+03 | 2.11 × 10+03 | 9.59 × 10−01 | 2.30 × 10+03 | 5.06× 10+02 | ||
f4 | 1 | Mean | 7.59 × 10+01 | 1.34× 10+02 | 3.19 × 10+03 | 1.69 × 10+02 | 1.46 × 10+02 | 4.52 × 10−02 |
S.D. | 4.33 × 10+01 | 1.52× 10+02 | 1.65 × 10+03 | 2.75 × 10+02 | 1.41 × 10+02 | 2.98 × 10−02 | ||
f5 | 1 | Mean | 3.75 × 10−02 | 7.50 × 10−02 | 9.04 × 10−01 | 1.72 × 10−01 | 4.52 × 10−02 | 1.85 × 10−02 |
S.D. | 1.12 × 10−02 | 2.09 × 10−02 | 3.67 × 10−01 | 4.36 × 10−02 | 1.67 × 10−02 | 6.22 × 10−03 | ||
f6 | 1 | Mean | 7.48 × 10+01 | 8.78× 10+01 | 1.16 × 10+02 | 1.32 × 10−03 | 7.67 × 10+01 | 4.79 × 10+01 |
S.D. | 1.62 × 10+01 | 1.83× 10+01 | 2.04 × 10+01 | 6.87 × 10−03 | 1.60 × 10+01 | 3.18 | ||
f7 | 0 | Mean | 1.38 | 1.40 × 10−01 | 9.00 | 1.22 | 1.14 | 9.05 × 10−01 |
S.D. | 8.89 × 10−01 | 3.52 × 10−01 | 1.08 | 5.83 × 10−01 | 7.77 × 10−01 | 8.47 × 10−01 | ||
f8 | 0 | Mean | 1.36 × 10−02 | 8.18 × 10−03 | 1.22 × 10+01 | 3.18 × 10−04 | 3.79 × 10−02 | 1.43 × 10−02 |
S.D. | 2.61 × 10−02 | 1.02 × 10−02 | 4.71 | 1.41 × 10−03 | 7.54 × 10−02 | 1.76 × 10−02 | ||
f9 | 0 | Mean | 1.60 × 10−01 | 1.04 × 10−01 | 2.08 × 10+01 | 8.95 × 10−02 | 7.84 × 10−08 | 2.99 × 10−02 |
S.D. | 2.65 × 10−01 | 1.96 × 10−01 | 8.57 | 1.44 × 10−01 | 4.22 × 10−07 | 4.91 × 10−02 | ||
f10 | 1 | Mean | 4.07 | 9.60 × 10−01 | 4.06 × 10+02 | 1.35 | 8.78 | 3.53 × 10−01 |
S.D. | 3.36 | 1.12 | 3.10 × 10+02 | 1.87 | 6.03 | 3.69 × 10−01 |
Function | Criteria | PSO-w | GPSO | LPSO | CLPSO | DMS-PSO | DCS-PSO |
F1 | Mean | 5.73 × 10+08 | 9.15 × 10+09 | 5.24 × 10+09 | 1.01 × 10+02 | 3.26 × 10+03 | 6.44 × 10+03 |
S.D. | 2.36 × 10+09 | 1.29 × 10+10 | 9.65 × 10+09 | 5.12 × 10−01 | 4.53 × 10+03 | 1.93 × 10+03 | |
F3 | Mean | 2.52 × 10+03 | 5.01 × 10+03 | 9.58 × 10+03 | 4.90 × 10+02 | 2.16 × 10+02 | 5.45 × 10+02 |
S.D. | 1.23 × 10+03 | 2.38 × 10+03 | 5.11 × 10+03 | 1.39 × 10+01 | 3.43 × 10+02 | 3.90 × 10−13 | |
F4 | Mean | 1.03 × 10+02 | 1.73 × 10+02 | 2.15 × 10+02 | 5.53 × 10+02 | 7.39 × 10+01 | 1.22 |
S.D. | 3.19 × 10+01 | 1.02 × 10+02 | 1.92 × 10+02 | 7.54 | 2.39 × 10+01 | 1.75 | |
F5 | Mean | 1.43 × 10+02 | 7.61 × 10+01 | 1.42 × 10+02 | 6.00 × 10+02 | 5.88 × 10+01 | 2.09 × 10+01 |
S.D. | 3.11 × 10+01 | 2.08 × 10+01 | 3.78 × 10+01 | 0.00 | 1.30 × 10+01 | 7.01 | |
F6 | Mean | 4.01 × 10+01 | 1.43 × 10+01 | 5.14 × 10+01 | 7.87 × 10+02 | 9.86 × 10−04 | 3.75 |
S.D. | 1.49 × 10+01 | 5.47 | 1.39 × 10+01 | 6.43 | 2.47 × 10−03 | 3.17 | |
F7 | Mean | 1.29 × 10+02 | 1.12 × 10+01 | 1.91 × 10+02 | 8.49 × 10+02 | 9.76 × 10+01 | 1.28 × 10+02 |
S.D. | 3.46 × 10+01 | 1.03 × 10+01 | 4.86 × 10+01 | 1.03 × 10+01 | 1.99 × 10+01 | 9.23 | |
F8 | Mean | 1.20 × 10+02 | 6.83 × 10+01 | 1.24 × 10+02 | 9.09 × 10+02 | 5.94 × 10+01 | 1.62 × 10+01 |
S.D. | 3.18 × 10+01 | 1.66 × 10+01 | 2.95 × 10+01 | 3.05 | 1.68 × 10+01 | 5.52 | |
F9 | Mean | 1.38 × 10+03 | 1.53 × 10+02 | 2.26 × 10+03 | 3.27 × 10+03 | 2.33 × 10+01 | 6.09 × 10−01 |
S.D. | 1.18 × 10+03 | 1.64 × 10+02 | 9.36 × 10+02 | 2.89 × 10+02 | 5.63 × 10+01 | 1.83 | |
F10 | Mean | 3.61 × 10+03 | 3.19 × 10+03 | 3.91 × 10+03 | 1.18 × 10+03 | 2.97 × 10+03 | 1.99 × 10+03 |
S.D. | 5.67 × 10+02 | 6.28 × 10+02 | 7.81 × 10+02 | 4.51 | 5.89 × 10+02 | 2.87 × 10+02 | |
F11 | Mean | 1.30 × 10+02 | 1.71 × 10+02 | 1.76 × 10+02 | 5.97 × 10+05 | 4.74 × 10+01 | 1.12 × 10+02 |
S.D. | 4.18 × 10+01 | 5.41 × 10+01 | 5.11 × 10+01 | 3.96 × 10+05 | 3.12 × 10+01 | 2.44 × 10+01 | |
F12 | Mean | 4.56 × 10+07 | 6.40 × 10+07 | 1.76 × 10+08 | 2.62 × 10+03 | 8.25 × 10+04 | 2.79 × 10+03 |
S.D. | 2.83 × 10+08 | 1.12 × 10+08 | 5.00 × 10+08 | 1.55 × 10+03 | 2.16 × 10+04 | 2.07 × 10+03 | |
F13 | Mean | 1.77 × 10+06 | 5.01 × 10+07 | 1.88 × 10+07 | 1.75 × 10+04 | 1.55 × 10+04 | 1.27 × 10+04 |
S.D. | 8.67 × 10+06 | 1.71 × 10+08 | 1.02 × 10+08 | 9.65 × 10+03 | 1.65 × 10+04 | 3.86 × 10+03 | |
F14 | Mean | 1.80 × 10+04 | 3.70 × 10+04 | 1.10 × 10+04 | 1.61 × 10+03 | 1.32 × 10+04 | 1.02 × 10+02 |
S.D. | 2.36 × 10+04 | 4.40 × 10+04 | 1.30 × 10+04 | 2.70 × 10+01 | 1.23 × 10+04 | 2.45 × 10+01 | |
F15 | Mean | 3.99 × 10+03 | 1.97 × 10+04 | 1.07 × 10+04 | 2.16 × 10+03 | 6.90 × 10+03 | 1.09 × 10+02 |
S.D. | 3.96 × 10+03 | 3.05 × 10+04 | 1.25 × 10+04 | 1.32 × 10+02 | 9.16 × 10+03 | 1.14 × 10+02 | |
F16 | Mean | 1.04 × 10+03 | 7.61 × 10+02 | 1.13 × 10+03 | 1.84 × 10+03 | 7.95 × 10+02 | 1.79 × 10+03 |
S.D. | 2.73 × 10+02 | 2.47 × 10+02 | 2.95 × 10+02 | 4.18 × 10+01 | 2.10 × 10+02 | 1.25 × 10+02 | |
F17 | Mean | 4.00 × 10+04 | 1.57 × 10+04 | 1.21 × 10+04 | 2.16 × 10+05 | 1.89 × 10+02 | 4.62 × 10+03 |
S.D. | 6.54 × 10+02 | 1.54 × 10+03 | 2.56 × 10+03 | 1.61 × 10+05 | 1.14 × 10+02 | 6.54 × 10+02 | |
F18 | Mean | 1.41 × 10+05 | 4.45 × 10+05 | 1.12 × 10+05 | 1.95 × 10+03 | 1.27 × 10+05 | 5.56 × 10+02 |
S.D. | 1.04 × 10+05 | 3.62 × 10+05 | 1.15 × 10+05 | 1.88 × 10+01 | 1.06 × 10+05 | 8.59 × 10+02 | |
F19 | Mean | 6.20 × 10+03 | 4.47 × 10+05 | 6.10 × 10+03 | 6.30 × 10+03 | 9.39 × 10+03 | 7.78 × 10+03 |
S.D. | 6.09 × 10+03 | 2.71 × 10+06 | 6.37 × 10+03 | 2.64 × 10+01 | 1.28 × 10+04 | 1.08 × 10+03 | |
F20 | Mean | 1.70 × 10+04 | 1.45 × 10+04 | 1.92 × 10+04 | 2.36 × 10+03 | 2.73 × 10+02 | 6.00 × 10+03 |
S.D. | 3.26 × 10+03 | 3.44 × 10+02 | 2.75 × 10+03 | 5.04 | 1.41 × 10+02 | 5.32 × 10+02 | |
F21 | Mean | 3.39 × 10+02 | 2.80 × 10+02 | 3.34 × 10+02 | 2.65 × 10+03 | 2.62 × 10+02 | 1.03 × 10+02 |
S.D. | 5.45 × 10+01 | 4.25 × 10+01 | 4.90 × 10+01 | 5.40 × 10+02 | 2.11 × 10+01 | 8.18 × 10−01 | |
F22 | Mean | 4.15 × 10+02 | 6.72 × 10+02 | 1.00 × 10+03 | 2.72 × 10+03 | 2.38 × 10+02 | 1.09 × 10+02 |
S.D. | 9.94 × 10+02 | 7.48 × 10+02 | 1.50 × 10+03 | 1.15 × 10+01 | 5.92 × 10+02 | 1.44 | |
F23 | Mean | 6.42 × 10+02 | 5.82 × 10+02 | 6.73 × 10+02 | 2.92 × 10+03 | 4.30 × 10+02 | 3.25 × 10+02 |
S.D. | 1.01 × 10+02 | 9.36 × 10+01 | 1.10 × 10+02 | 1.21 × 10+01 | 2.10 × 10+01 | 7.91 | |
F24 | Mean | 6.61 × 10+02 | 6.34 × 10+02 | 7.14 × 10+02 | 2.89 × 10+03 | 4.94 × 10+02 | 3.53 × 10+02 |
S.D. | 7.19 × 10+01 | 8.33 × 10+01 | 8.07 × 10+01 | 2.23 × 10−01 | 2.60 × 10+01 | 1.09 × 10+01 | |
F25 | Mean | 3.95 × 10+02 | 4.13 × 10+02 | 4.47 × 10+02 | 3.47 × 10+03 | 3.97 × 10+02 | 4.26 × 10+02 |
S.D. | 1.50 × 10+01 | 3.64 × 10+01 | 3.05 × 10+01 | 7.89 × 10+02 | 1.53 | 2.34 × 10+01 | |
F26 | Mean | 1.24 × 10+03 | 1.51 × 10+03 | 2.45 × 10+03 | 3.21 × 10+03 | 1.49 × 10+03 | 4.02 × 10+02 |
S.D. | 1.29 × 10+03 | 7.71 × 10+02 | 1.36 × 10+03 | 6.70 | 5.18 × 10+02 | 5.94 × 10+01 | |
F27 | Mean | 5.88 × 10+02 | 5.80 × 10+02 | 6.21 × 10+02 | 3.21 × 10+03 | 5.38 × 10+02 | 4.03 × 10+02 |
S.D. | 5.87 × 10+01 | 5.54 × 10+01 | 5.36 × 10+01 | 5.65 | 1.25 × 10+01 | 4.45 | |
F28 | Mean | 4.34 × 10+02 | 5.05 × 10+02 | 5.49 × 10+02 | 3.45 × 10+03 | 3.40 × 10+02 | 3.88 × 10+02 |
S.D. | 4.76 × 10+01 | 7.23 × 10+01 | 1.46 × 10+02 | 4.85 | 5.79 × 10+01 | 8.36 | |
F29 | Mean | 6.25 × 10+02 | 8.71 × 10+02 | 6.24 × 10+02 | 8.91 × 10+02 | 6.03 × 10+02 | 6.38 × 10+02 |
S.D. | 4.87 × 10+01 | 6.31 × 10+01 | 2.17 × 10+02 | 7.46 × 10+01 | 1.55 × 10+02 | 7.88 × 10+01 | |
F30 | Mean | 4.52 × 10+04 | 4.50 × 10+05 | 3.78 × 10+05 | 1.12 × 10+03 | 5.81 × 10+03 | 3.14 × 10+05 |
S.D. | 1.10 × 10+05 | 1.41 × 10+06 | 4.93 × 10+05 | 2.22 × 10+03 | 2.34 × 10+03 | 5.14 × 10+05 |
Average rank | Algorithm | Ranking |
1 | DCS-PSO | 2.03 |
2 | DMS-PSO | 2.17 |
3 | PSO-w | 3.69 |
4 | GPSO | 4.07 |
5 | CLPSO | 4.45 |
6 | LPSO | 4.59 |
PC configuration | Information |
CPU | Inter(R) Core(TM) i7-8700 |
Frequency | 3.2GHz CPU |
RAM | 16.0GB |
Operating system | Windows 10(64Bit) |
Language | Python 3.9.13 |
Algorithm | Parameters |
COA | max_iter1=105, max_iter2=1.5×105, α=0.3 |
DCOA | h=3000, γ=0.25, ξ = 1.5 |
PCOA | n_group=3, max_iter1=105, max_iter2=1.5×105, α=0.3 |
DCS-PSO | N=20, w=0.4, c1=c2=2, Vmax=0.2×Range, max_iter=1000 |
Function | Criteria | COA | DCOA | PCOA | DCS-PSO |
fJ | E.T. | 30.8898 | 0.5753 | 17.2234 | 0.2793 |
Mean | 0.0000 | 0.0005 | 0.0000 | 0.0000 | |
S.D. | 0.0000 | 0.0014 | 0.0000 | 0.0000 | |
fC | E.T. | 7.6564 | 0.6478 | 10.3990 | 0.3610 |
Mean | −1.0316 | −1.0312 | −1.0316 | −1.0316 | |
S.D. | 0.0000 | 0.0012 | 0.0000 | 0.0000 | |
fG | E.T. | 14.1867 | 1.0656 | 14.0814 | 0.6097 |
Mean | 3.0022 | 3.0100 | 3.0024 | 3.0000 | |
S.D. | 0.0022 | 0.0126 | 0.0022 | 0.0000 | |
fB | E.T. | 1.5663 | 0.8653 | 5.7455 | 0.4513 |
Mean | 0.3983 | 0.4003 | 0.3980 | 0.3979 | |
S.D. | 0.0003 | 0.0039 | 0.0001 | 0.0000 | |
fR | E.T. | 4.4150 | 0.6652 | 5.8342 | 0.4366 |
Mean | −1.9993 | −1.9652 | −1.9992 | −2.0000 | |
S.D. | 0.0006 | 0.0558 | 0.0008 | 0.0000 | |
fS | E.T. | 5.8609 | 2.1743 | 13.8397 | 2.0615 |
Mean | −186.6428 | -186.4517 | −186.6861 | −186.7309 | |
S.D. | 0.0711 | 0.4203 | 0.0680 | 0.0000 |
Name | Functions | Range | Optimum |
Sphere | f1(x)=D∑i=1x2i | [−100, 100] | 0 |
Schwefel's 2.22 | f2(x)=D∑i=1|xi|+D∏i=1|xi| | [−10, 10] | 0 |
Schwefel's 1.2 | f3(x)=D∑i=1(i∑j=1xj)2 | [−100, 100] | 0 |
Rosenbrock | f4(x)=D−1∑i=1[100(xi+1−x2i)2+(xi−1)2] | [−10, 10] | 0 |
Noise | f5(x)=D∑i=1ix4i+random[0,1) | [−1.28, 1.28] | 0 |
Rastrigin | f6(x)=D∑i=1(x2i−10cos(2πxi)+10) | [−5.12, 5.12] | 0 |
Ackley | f7(x)=−20exp(−0.2√1DD∑i=1x2i)−exp(1DD∑i=1cos(2πxi))+20+e | [−32, 32] | 0 |
Griewank | f8(x)=14000D∑i=1x2i−D∏i=1cos(xi√i)+1 | [−600, 600] | 0 |
Penalized 1 | f9(x)=πD{10sin2(πy1)+D−1∑i=1(yi−1)2[1+10sin2(πyi+1)]+(yD−1)2}+D∑i=1u(xi,10,100,4) | [−50, 50] | 0 |
Penalized 2 | f10(x)=110{sin2(3πx1)+D∑i=1(xi−1)2[1+sin2(3πxi+1)]+(xD−1)2[1+sin2(2πxD)]}+D∑i=1u(xi,5,100,4) | [−50, 50] | 0 |
Algorithm | Parameters |
PSO-w | w:0.4, c1=c2=2.0, Vmax=0.2×Range |
GPSO | w:0.9∼0.4, c1=c2=2.0, Vmax=0.2×Range |
LPSO | χ=0.7298, c1=c2=1.49445, Vmax=0.2×Range |
CLPSO | w:0.9∼0.4, c=1.49445, Vmax=0.2×Range, m=7 |
DMS-PSO | χ=0.7298, c1=c2=1.49445, Vmax=0.2×Range, M=4, R=10 |
DCS-PSO | w:0.4,c1=c2=2.0,Vmax=0.2×Range, h=3000,γ=0.15,ξ=1.5 |
Function | k | Criteria | PSO-w | GPSO | LPSO | CLPSO | DMS-PSO | DCS-PSO |
f1 | 1 | Mean | 1.84 × 10−144 | 1.47 × 10−73 | 2.78 × 10−169 | 2.18 × 10−21 | 6.30 × 10−62 | 2.35 × 10−151 |
S.D. | 9.13 × 10−144 | 1.04 × 10−72 | 0.00 | 5.44 × 10−21 | 4.25 × 10−61 | 7.75 × 10−151 | ||
f2 | 1 | Mean | 1.89 × 10−81 | 6.24 × 10−38 | 7.59 × 10−09 | 4.30 × 10−14 | 3.45 × 10−20 | 9.65 × 10−84 |
S.D. | 1.26 × 10−80 | 2.99 × 10−37 | 3.80 × 10−08 | 1.39 × 10−13 | 5.86 × 10−20 | 4.81 × 10−83 | ||
f3 | 1 | Mean | 1.22 × 10−44 | 9.93 × 10−17 | 6.11 × 10−13 | 2.18 × 10−20 | 6.13 × 10−10 | 1.21 × 10−48 |
S.D. | 7.42 × 10−44 | 3.36 × 10−16 | 4.32 × 10−12 | 3.25 × 10−20 | 3.33 × 10−09 | 2.38 × 10−48 | ||
f4 | 1 | Mean | 2.79 | 2.80 | 4.00 | 2.49 | 1.60 | 9.11 |
S.D. | 1.39 | 1.16 | 2.56 | 4.95 × 10−01 | 9.79 × 10−01 | 1.89 | ||
f5 | 1 | Mean | 7.23 × 10−04 | 1.02 × 10−03 | 3.25 × 10−03 | 1.49 × 10−03 | 8.80 × 10−04 | 2.15 × 10−04 |
S.D. | 3.90 × 10−04 | 5.46 × 10−04 | 2.43 × 10−03 | 7.09 × 10−04 | 4.49 × 10−04 | 1.58 × 10−04 | ||
f6 | 0 | Mean | 5.99 | 2.13 | 9.35 | 0.00 | 2.67 | 0.00 |
S.D. | 2.74 | 1.15 | 3.53 | 0.00 | 1.36 | 0.00 | ||
f7 | 0 | Mean | 4.14 × 10−15 | 4.07 × 10−15 | 4.67 × 10−15 | 3.16 × 10−14 | 4.00 × 10−15 | 3.91 × 10−15 |
S.D. | 7.03 × 10−16 | 5.02 × 10−16 | 1.62 × 10−15 | 5.29 × 10−14 | 0.00 | 5.62 × 10−16 | ||
f8 | 1 | Mean | 7.44 × 10−02 | 6.88 × 10−02 | 1.06 × 10−01 | 1.64 × 10−03 | 7.12 × 10−02 | 6.65 × 10−02 |
S.D. | 3.41 × 10−02 | 3.13 × 10−02 | 5.81 × 10−02 | 7.61 × 10−04 | 3.55 × 10−02 | 3.75 × 10−02 | ||
f9 | 1 | Mean | 1.24 × 10−02 | 4.71 × 10−32 | 6.22 × 10−03 | 2.99 × 10−30 | 1.38 × 10−28 | 1.78 × 10−01 |
S.D. | 6.16 × 10−02 | 1.11 × 10−47 | 4.40 × 10−02 | 3.07 × 10−30 | 2.65 × 10−28 | 6.11 × 10−02 | ||
f10 | 1 | Mean | 3.56 × 10−04 | 9.33 × 10−20 | 9.80 × 10−04 | 3.16 × 10−27 | 2.57 × 10−13 | 1.35 × 10−32 |
S.D. | 2.52 × 10−03 | 1.90 × 10−19 | 6.86 × 10−03 | 6.79 × 10−27 | 4.33 × 10−13 | 8.29 × 10−48 |
Function | k | Criteria | PSO-w | GPSO | LPSO | CLPSO | DMS-PSO | DCS-PSO |
f1 | 1 | Mean | 1.28 × 10−39 | 3.42 × 10−14 | 3.82 | 3.83 × 10−13 | 6.76 × 10−28 | 1.05 × 10−40 |
S.D. | 2.48 × 10−39 | 1.20 × 10−13 | 2.18× 10+01 | 9.48 × 10−13 | 5.49 × 10−28 | 2.01 × 10−40 | ||
Rank | 2 | 4 | 6 | 5 | 3 | 1 | ||
f2 | 1 | Mean | 1.67 × 10−09 | 1.27 × 10−08 | 1.83 | 4.45 × 10−08 | 3.74 × 10−13 | 2.49 × 10−23 |
S.D. | 1.15 × 10−08 | 2.46 × 10−08 | 1.83 | 2.84 × 10−08 | 2.33 × 10−13 | 5.70 × 10−23 | ||
Rank | 3 | 4 | 6 | 5 | 2 | 1 | ||
f3 | 1 | Mean | 6.50 × 10−01 | 9.58 × 10+01 | 4.79× 10+02 | 4.37 × 10−04 | 1.16 × 10+01 | 5.90 |
S.D. | 5.22 × 10−01 | 5.10 × 10+01 | 3.08× 10+02 | 3.04 × 10−04 | 5.95 | 4.45 | ||
Rank | 2 | 5 | 6 | 1 | 4 | 3 | ||
f4 | 1 | Mean | 2.98× 10+01 | 4.17 × 10+01 | 9.42× 10+01 | 1.81 × 10 +01 | 2.26 × 10+01 | 2.21 × 10+01 |
S.D. | 2.16× 10+01 | 2.87 × 10+01 | 5.46× 10+01 | 3.80 | 2.29 | 6.32 × 10−01 | ||
Rank | 4 | 5 | 6 | 1 | 3 | 2 | ||
f5 | 1 | Mean | 8.74 × 10−03 | 1.89 × 10−02 | 8.00 × 10−02 | 3.83 × 10−02 | 9.99 × 10−03 | 5.12 × 10−03 |
S.D. | 3.28 × 10−03 | 6.58 × 10−03 | 3.51 × 10−02 | 1.47 × 10−02 | 3.98 × 10−03 | 1.69 × 10−03 | ||
Rank | 2 | 4 | 6 | 5 | 3 | 1 | ||
f6 | 1 | Mean | 3.88 × 10+01 | 3.34 × 10+01 | 5.12 × 10+01 | 4.40 × 10−11 | 2.73 × 10+01 | 1.31 × 10+02 |
S.D. | 8.06 | 8.19 | 1.62 × 10+01 | 6.58 × 10−11 | 8.88 × 10+00 | 1.01 × 10+01 | ||
Rank | 4 | 3 | 5 | 1 | 2 | 6 | ||
f7 | 1 | Mean | 3.81 × 10−01 | 1.14 × 10−06 | 4.32 | 3.15 × 10−03 | 4.37 × 10−14 | 1.20 × 10−14 |
S.D. | 6.47 × 10−01 | 1.04 × 10−06 | 1.16 | 2.18 × 10−02 | 4.86 × 10−14 | 3.90 × 10−15 | ||
Rank | 5 | 3 | 6 | 4 | 2 | 1 | ||
f8 | 0 | Mean | 1.46 × 10−02 | 1.38 × 10−02 | 4.14 × 10−01 | 1.41 × 10−12 | 4.44 × 10−03 | 1.15 × 10−02 |
S.D. | 1.70 × 10−02 | 1.55 × 10−02 | 3.12 × 10−01 | 2.52 × 10−12 | 4.99 × 10−03 | 1.40 × 10−02 | ||
Rank | 5 | 4 | 6 | 1 | 2 | 3 | ||
f9 | 0 | Mean | 1.39 × 10−01 | 1.04 × 10−02 | 3.87 | 2.32 × 10−14 | 1.08 × 10−26 | 2.96 × 10−02 |
S.D. | 2.62 × 10−01 | 3.14 × 10−02 | 2.66 | 5.72 × 10−14 | 2.24 × 10−26 | 4.74 × 10−02 | ||
Rank | 4 | 3 | 6 | 2 | 1 | 4 | ||
f10 | 1 | Mean | 2.28 × 10−01 | 1.10 × 10−02 | 2.16× 10+01 | 1.03 × 10−11 | 5.31 × 10−03 | 6.22 × 10−32 |
S.D. | 6.63 × 10−01 | 4.43 × 10−02 | 1.12× 10+01 | 1.45 × 10−11 | 4.13 × 10−04 | 7.45 × 10−32 | ||
Rank | 5 | 4 | 6 | 2 | 3 | 1 | ||
Ave. rank | 3.6 | 3.9 | 5.9 | 2.7 | 2.5 | 2.3 | ||
Final rank | 4 | 5 | 6 | 3 | 2 | 1 |
Function | k | Criteria | PSO-w | GPSO | LPSO | CLPSO | DMS-PSO | DCS-PSO |
f1 | 1 | Mean | 5.20 × 10−17 | 3.68 × 10−04 | 1.30 × 10+03 | 5.70 × 10−05 | 1.26 × 10−06 | 3.93 × 10−18 |
S.D. | 1.91 × 10−16 | 4.22 × 10−04 | 4.91 × 10+02 | 6.50 × 10−05 | 3.14 × 10−06 | 1.05 × 10−17 | ||
f2 | 0 | Mean | 9.06 × 10−06 | 2.01 × 10−01 | 2.26 × 10+01 | 2.00 × 10−01 | 1.48 × 10−04 | 5.87 × 10−07 |
S.D. | 5.83 × 10−05 | 1.41 | 4.26 | 1.41 | 6.99 × 10−04 | 2.82 × 10−06 | ||
f3 | 1 | Mean | 4.52 × 10+02 | 4.38 × 10+03 | 7.71 × 10+03 | 1.55 | 1.90 × 10+03 | 1.33 × 10+03 |
S.D. | 1.68 × 10+02 | 1.88× 10+03 | 2.11 × 10+03 | 9.59 × 10−01 | 2.30 × 10+03 | 5.06× 10+02 | ||
f4 | 1 | Mean | 7.59 × 10+01 | 1.34× 10+02 | 3.19 × 10+03 | 1.69 × 10+02 | 1.46 × 10+02 | 4.52 × 10−02 |
S.D. | 4.33 × 10+01 | 1.52× 10+02 | 1.65 × 10+03 | 2.75 × 10+02 | 1.41 × 10+02 | 2.98 × 10−02 | ||
f5 | 1 | Mean | 3.75 × 10−02 | 7.50 × 10−02 | 9.04 × 10−01 | 1.72 × 10−01 | 4.52 × 10−02 | 1.85 × 10−02 |
S.D. | 1.12 × 10−02 | 2.09 × 10−02 | 3.67 × 10−01 | 4.36 × 10−02 | 1.67 × 10−02 | 6.22 × 10−03 | ||
f6 | 1 | Mean | 7.48 × 10+01 | 8.78× 10+01 | 1.16 × 10+02 | 1.32 × 10−03 | 7.67 × 10+01 | 4.79 × 10+01 |
S.D. | 1.62 × 10+01 | 1.83× 10+01 | 2.04 × 10+01 | 6.87 × 10−03 | 1.60 × 10+01 | 3.18 | ||
f7 | 0 | Mean | 1.38 | 1.40 × 10−01 | 9.00 | 1.22 | 1.14 | 9.05 × 10−01 |
S.D. | 8.89 × 10−01 | 3.52 × 10−01 | 1.08 | 5.83 × 10−01 | 7.77 × 10−01 | 8.47 × 10−01 | ||
f8 | 0 | Mean | 1.36 × 10−02 | 8.18 × 10−03 | 1.22 × 10+01 | 3.18 × 10−04 | 3.79 × 10−02 | 1.43 × 10−02 |
S.D. | 2.61 × 10−02 | 1.02 × 10−02 | 4.71 | 1.41 × 10−03 | 7.54 × 10−02 | 1.76 × 10−02 | ||
f9 | 0 | Mean | 1.60 × 10−01 | 1.04 × 10−01 | 2.08 × 10+01 | 8.95 × 10−02 | 7.84 × 10−08 | 2.99 × 10−02 |
S.D. | 2.65 × 10−01 | 1.96 × 10−01 | 8.57 | 1.44 × 10−01 | 4.22 × 10−07 | 4.91 × 10−02 | ||
f10 | 1 | Mean | 4.07 | 9.60 × 10−01 | 4.06 × 10+02 | 1.35 | 8.78 | 3.53 × 10−01 |
S.D. | 3.36 | 1.12 | 3.10 × 10+02 | 1.87 | 6.03 | 3.69 × 10−01 |
Function | Criteria | PSO-w | GPSO | LPSO | CLPSO | DMS-PSO | DCS-PSO |
F1 | Mean | 5.73 × 10+08 | 9.15 × 10+09 | 5.24 × 10+09 | 1.01 × 10+02 | 3.26 × 10+03 | 6.44 × 10+03 |
S.D. | 2.36 × 10+09 | 1.29 × 10+10 | 9.65 × 10+09 | 5.12 × 10−01 | 4.53 × 10+03 | 1.93 × 10+03 | |
F3 | Mean | 2.52 × 10+03 | 5.01 × 10+03 | 9.58 × 10+03 | 4.90 × 10+02 | 2.16 × 10+02 | 5.45 × 10+02 |
S.D. | 1.23 × 10+03 | 2.38 × 10+03 | 5.11 × 10+03 | 1.39 × 10+01 | 3.43 × 10+02 | 3.90 × 10−13 | |
F4 | Mean | 1.03 × 10+02 | 1.73 × 10+02 | 2.15 × 10+02 | 5.53 × 10+02 | 7.39 × 10+01 | 1.22 |
S.D. | 3.19 × 10+01 | 1.02 × 10+02 | 1.92 × 10+02 | 7.54 | 2.39 × 10+01 | 1.75 | |
F5 | Mean | 1.43 × 10+02 | 7.61 × 10+01 | 1.42 × 10+02 | 6.00 × 10+02 | 5.88 × 10+01 | 2.09 × 10+01 |
S.D. | 3.11 × 10+01 | 2.08 × 10+01 | 3.78 × 10+01 | 0.00 | 1.30 × 10+01 | 7.01 | |
F6 | Mean | 4.01 × 10+01 | 1.43 × 10+01 | 5.14 × 10+01 | 7.87 × 10+02 | 9.86 × 10−04 | 3.75 |
S.D. | 1.49 × 10+01 | 5.47 | 1.39 × 10+01 | 6.43 | 2.47 × 10−03 | 3.17 | |
F7 | Mean | 1.29 × 10+02 | 1.12 × 10+01 | 1.91 × 10+02 | 8.49 × 10+02 | 9.76 × 10+01 | 1.28 × 10+02 |
S.D. | 3.46 × 10+01 | 1.03 × 10+01 | 4.86 × 10+01 | 1.03 × 10+01 | 1.99 × 10+01 | 9.23 | |
F8 | Mean | 1.20 × 10+02 | 6.83 × 10+01 | 1.24 × 10+02 | 9.09 × 10+02 | 5.94 × 10+01 | 1.62 × 10+01 |
S.D. | 3.18 × 10+01 | 1.66 × 10+01 | 2.95 × 10+01 | 3.05 | 1.68 × 10+01 | 5.52 | |
F9 | Mean | 1.38 × 10+03 | 1.53 × 10+02 | 2.26 × 10+03 | 3.27 × 10+03 | 2.33 × 10+01 | 6.09 × 10−01 |
S.D. | 1.18 × 10+03 | 1.64 × 10+02 | 9.36 × 10+02 | 2.89 × 10+02 | 5.63 × 10+01 | 1.83 | |
F10 | Mean | 3.61 × 10+03 | 3.19 × 10+03 | 3.91 × 10+03 | 1.18 × 10+03 | 2.97 × 10+03 | 1.99 × 10+03 |
S.D. | 5.67 × 10+02 | 6.28 × 10+02 | 7.81 × 10+02 | 4.51 | 5.89 × 10+02 | 2.87 × 10+02 | |
F11 | Mean | 1.30 × 10+02 | 1.71 × 10+02 | 1.76 × 10+02 | 5.97 × 10+05 | 4.74 × 10+01 | 1.12 × 10+02 |
S.D. | 4.18 × 10+01 | 5.41 × 10+01 | 5.11 × 10+01 | 3.96 × 10+05 | 3.12 × 10+01 | 2.44 × 10+01 | |
F12 | Mean | 4.56 × 10+07 | 6.40 × 10+07 | 1.76 × 10+08 | 2.62 × 10+03 | 8.25 × 10+04 | 2.79 × 10+03 |
S.D. | 2.83 × 10+08 | 1.12 × 10+08 | 5.00 × 10+08 | 1.55 × 10+03 | 2.16 × 10+04 | 2.07 × 10+03 | |
F13 | Mean | 1.77 × 10+06 | 5.01 × 10+07 | 1.88 × 10+07 | 1.75 × 10+04 | 1.55 × 10+04 | 1.27 × 10+04 |
S.D. | 8.67 × 10+06 | 1.71 × 10+08 | 1.02 × 10+08 | 9.65 × 10+03 | 1.65 × 10+04 | 3.86 × 10+03 | |
F14 | Mean | 1.80 × 10+04 | 3.70 × 10+04 | 1.10 × 10+04 | 1.61 × 10+03 | 1.32 × 10+04 | 1.02 × 10+02 |
S.D. | 2.36 × 10+04 | 4.40 × 10+04 | 1.30 × 10+04 | 2.70 × 10+01 | 1.23 × 10+04 | 2.45 × 10+01 | |
F15 | Mean | 3.99 × 10+03 | 1.97 × 10+04 | 1.07 × 10+04 | 2.16 × 10+03 | 6.90 × 10+03 | 1.09 × 10+02 |
S.D. | 3.96 × 10+03 | 3.05 × 10+04 | 1.25 × 10+04 | 1.32 × 10+02 | 9.16 × 10+03 | 1.14 × 10+02 | |
F16 | Mean | 1.04 × 10+03 | 7.61 × 10+02 | 1.13 × 10+03 | 1.84 × 10+03 | 7.95 × 10+02 | 1.79 × 10+03 |
S.D. | 2.73 × 10+02 | 2.47 × 10+02 | 2.95 × 10+02 | 4.18 × 10+01 | 2.10 × 10+02 | 1.25 × 10+02 | |
F17 | Mean | 4.00 × 10+04 | 1.57 × 10+04 | 1.21 × 10+04 | 2.16 × 10+05 | 1.89 × 10+02 | 4.62 × 10+03 |
S.D. | 6.54 × 10+02 | 1.54 × 10+03 | 2.56 × 10+03 | 1.61 × 10+05 | 1.14 × 10+02 | 6.54 × 10+02 | |
F18 | Mean | 1.41 × 10+05 | 4.45 × 10+05 | 1.12 × 10+05 | 1.95 × 10+03 | 1.27 × 10+05 | 5.56 × 10+02 |
S.D. | 1.04 × 10+05 | 3.62 × 10+05 | 1.15 × 10+05 | 1.88 × 10+01 | 1.06 × 10+05 | 8.59 × 10+02 | |
F19 | Mean | 6.20 × 10+03 | 4.47 × 10+05 | 6.10 × 10+03 | 6.30 × 10+03 | 9.39 × 10+03 | 7.78 × 10+03 |
S.D. | 6.09 × 10+03 | 2.71 × 10+06 | 6.37 × 10+03 | 2.64 × 10+01 | 1.28 × 10+04 | 1.08 × 10+03 | |
F20 | Mean | 1.70 × 10+04 | 1.45 × 10+04 | 1.92 × 10+04 | 2.36 × 10+03 | 2.73 × 10+02 | 6.00 × 10+03 |
S.D. | 3.26 × 10+03 | 3.44 × 10+02 | 2.75 × 10+03 | 5.04 | 1.41 × 10+02 | 5.32 × 10+02 | |
F21 | Mean | 3.39 × 10+02 | 2.80 × 10+02 | 3.34 × 10+02 | 2.65 × 10+03 | 2.62 × 10+02 | 1.03 × 10+02 |
S.D. | 5.45 × 10+01 | 4.25 × 10+01 | 4.90 × 10+01 | 5.40 × 10+02 | 2.11 × 10+01 | 8.18 × 10−01 | |
F22 | Mean | 4.15 × 10+02 | 6.72 × 10+02 | 1.00 × 10+03 | 2.72 × 10+03 | 2.38 × 10+02 | 1.09 × 10+02 |
S.D. | 9.94 × 10+02 | 7.48 × 10+02 | 1.50 × 10+03 | 1.15 × 10+01 | 5.92 × 10+02 | 1.44 | |
F23 | Mean | 6.42 × 10+02 | 5.82 × 10+02 | 6.73 × 10+02 | 2.92 × 10+03 | 4.30 × 10+02 | 3.25 × 10+02 |
S.D. | 1.01 × 10+02 | 9.36 × 10+01 | 1.10 × 10+02 | 1.21 × 10+01 | 2.10 × 10+01 | 7.91 | |
F24 | Mean | 6.61 × 10+02 | 6.34 × 10+02 | 7.14 × 10+02 | 2.89 × 10+03 | 4.94 × 10+02 | 3.53 × 10+02 |
S.D. | 7.19 × 10+01 | 8.33 × 10+01 | 8.07 × 10+01 | 2.23 × 10−01 | 2.60 × 10+01 | 1.09 × 10+01 | |
F25 | Mean | 3.95 × 10+02 | 4.13 × 10+02 | 4.47 × 10+02 | 3.47 × 10+03 | 3.97 × 10+02 | 4.26 × 10+02 |
S.D. | 1.50 × 10+01 | 3.64 × 10+01 | 3.05 × 10+01 | 7.89 × 10+02 | 1.53 | 2.34 × 10+01 | |
F26 | Mean | 1.24 × 10+03 | 1.51 × 10+03 | 2.45 × 10+03 | 3.21 × 10+03 | 1.49 × 10+03 | 4.02 × 10+02 |
S.D. | 1.29 × 10+03 | 7.71 × 10+02 | 1.36 × 10+03 | 6.70 | 5.18 × 10+02 | 5.94 × 10+01 | |
F27 | Mean | 5.88 × 10+02 | 5.80 × 10+02 | 6.21 × 10+02 | 3.21 × 10+03 | 5.38 × 10+02 | 4.03 × 10+02 |
S.D. | 5.87 × 10+01 | 5.54 × 10+01 | 5.36 × 10+01 | 5.65 | 1.25 × 10+01 | 4.45 | |
F28 | Mean | 4.34 × 10+02 | 5.05 × 10+02 | 5.49 × 10+02 | 3.45 × 10+03 | 3.40 × 10+02 | 3.88 × 10+02 |
S.D. | 4.76 × 10+01 | 7.23 × 10+01 | 1.46 × 10+02 | 4.85 | 5.79 × 10+01 | 8.36 | |
F29 | Mean | 6.25 × 10+02 | 8.71 × 10+02 | 6.24 × 10+02 | 8.91 × 10+02 | 6.03 × 10+02 | 6.38 × 10+02 |
S.D. | 4.87 × 10+01 | 6.31 × 10+01 | 2.17 × 10+02 | 7.46 × 10+01 | 1.55 × 10+02 | 7.88 × 10+01 | |
F30 | Mean | 4.52 × 10+04 | 4.50 × 10+05 | 3.78 × 10+05 | 1.12 × 10+03 | 5.81 × 10+03 | 3.14 × 10+05 |
S.D. | 1.10 × 10+05 | 1.41 × 10+06 | 4.93 × 10+05 | 2.22 × 10+03 | 2.34 × 10+03 | 5.14 × 10+05 |
Average rank | Algorithm | Ranking |
1 | DCS-PSO | 2.03 |
2 | DMS-PSO | 2.17 |
3 | PSO-w | 3.69 |
4 | GPSO | 4.07 |
5 | CLPSO | 4.45 |
6 | LPSO | 4.59 |