Quadratic Surrogate Attractor for Particle Swarm Optimization
Organizations: Center for Automotive Research, The Ohio State University, Columbus, OH 43212 USA
Abstract
This paper presents a particle swarm optimization algorithm that leverages surrogate modeling to replace the conventional global best solution with the minimum of an n-dimensional quadratic form, providing a better-conditioned dynamic attractor for the swarm. This refined convergence target, informed by the local landscape, enhances global convergence behavior and increases robustness against premature convergence and noise, while incurring only minimal computational overhead. The surrogate-augmented approach is evaluated against the standard algorithm through a numerical study on a set of benchmark optimization functions that exhibit diverse landscapes. To ensure statistical significance, 400 independent runs are conducted for each function and algorithm, and the results are analyzed based on their statistical characteristics and corresponding distributions. The quadratic surrogate attractor consistently outperforms the conventional algorithm across all tested functions. The improvement is particularly pronounced for quasi-convex functions, where the surrogate model can exploit the underlying convex-like structure of the landscape.
Figures & tables
| 0.72984 | 2.8 | 2.05 | 2 | 200 | 52 | 1.2 |
| Function | N p | Bounds | Mean Q.S. | Mean Std. | Rel. Diff. | Time Q.S. [ ] | Time Std. [ ] | Rel. Diff. | IQR 25-75 Q.S. | IQR 25-75 Std. |
|---|---|---|---|---|---|---|---|---|---|---|
| Ackley 2D | 6 | [[-32.768, 32.768],[-32.768, 32.768]] | 1.307e-2 | 3.294e-1 | -96.03% | 91.27 | 84.41 | + 8.12 % | (1.112e-11)–(3.568e-7) | (1.037e-7)–(1.478e-4) |
| Griewank 2D | 6 | [[-600, 600],[-600, 600]] | 2.119e-1 | 7.927e-1 | -73.27% | 117.55 | 102.41 | + 14.78 % | (1.265e-12)–(1.300e-1) | (1.010e-8)–(2.610e-1) |
| Sphere 2D | 6 | [[-10, 10],[-10, 10]] | 1.182e-7 | 1.235e-7 | -4.292% | 67.19 | 60.32 | + 11.38 % | (1.791e-18)–(5.523e-11) | (2.838e-14)–(7.071e-10) |
| Sphere 3D | 10 | [[-10, 10],[-10, 10],[-10, 10]] | 3.122e-7 | 2.334e-6 | -86.62% | 75.09 | 67.59 | + 11.10 % | (3.306e-17)–(1.637e-10) | (1.525e-11)–(4.652e-8) |
| Flower 2D | 6 | [[-100, 100],[-100, 100]] | 8.045e-4 | 4.776e-1 | -99.83% | 76.43 | 69.21 | + 10.43 % | (1.501e-12)–(7.327e-8) | (1.132e-6)–(7.560e-4) |
| Flower 3D | 10 | [[-100, 100],[-100, 100],[-100, 100]] | 4.565e-4 | 5.374e-1 | -99.92% | 96.25 | 87.39 | + 10.14 % | (1.221e-9)–(2.168e-6) | (9.853e-5)–(3.932e-2) |
Particle Swarm Optimization
Particle Swarm Optimization (PSO) uses interacting candidate solutions to search complex spaces and tune models, designs, and control systems.
See all 21 papers