智能优化算法：人工电场优化算法-附代码

本文链接：https://blog.csdn.net/u011835903/article/details/118929142

人工电场优化算法（AEFA）是一种基于库仑定律的元启发式优化算法，通过模拟电荷间的吸引力寻找全局最优解。算法在非线性优化问题中表现出强搜索能力和高精度。在每次迭代中，电荷的运动由库仑力、加速度和速度更新规则决定，电荷量根据适应度函数动态调整，从而逐步收敛到最佳解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

智能优化算法：人工电场优化算法

文章目录

智能优化算法：人工电场优化算法

摘要：人工电场算法（Artificial Electric Field Algorithm，AEFA）是由 Anita 和 Anupam Yadav于 2019 年提出的元启发式优化算法。AEFA 算法是一个杰出的非线性优化算法，具有全局搜索能力强、优化精度高和适应范围广等特点。

1.算法原理

人工电场算法（AEFA）受库仑定律的启发，通过模拟带电粒子在静电场的运动，将其演化成随机搜索最优解的过程。电荷在静电力的作用下相互吸引或排斥，使电荷能在搜索空间中移动。在 AEFA 算法中，仅考虑电荷的吸引力，而忽略电荷的排斥力，因而电荷量大的带电粒子能吸引所有其他电荷量较低的粒子向其靠近。

在搜索空间中，每一个电荷代表一个可行解，其强度由他们的电荷量来衡量，电荷的电荷量越大，表明该电荷越接近理论最优解。如图1所示，电荷面积的大小表示该电荷量的大小，电荷 $Q_1$ 分别受到其它三个电荷的吸引力，根据运动定律，形成一个合力 $F$ 和该方向的加速度。由于电荷 $Q_4$ 的电荷量最大，其吸引力也越大，电荷 $Q_1$ 合力 $F$ 的方向更接近 $Q_1$ 与 $Q_4$ 之间的中心连线。因此，AEFA 算法通过模拟电荷间的相互作用力，当搜索空间存在电荷量大的电荷时，其它电荷都向电荷量大的方向靠近，使算法收敛到最优解。

在这里插入图片描述

图1.库仑力的作用。

AEFA 算法可看作是一个遵循库仑静电力定律和运动定律孤立的电荷系统。在此基础上，定义了算法的物理性质。

设 $d$ 维搜索空间中，第 $i$ 个电荷的位置为：
$X_i=(x_i^1,x_i^2,...,x_i^d);i=1,2,...,N\tag{1}$
式中， $x_i^d$ 为电荷 $i$ 在第 $d$ 维上的位置， $N$ 为电荷的总数。

电荷 $i$ 在时刻 $t$ 获得的最佳适应度值的位置由下式确定：
$p_i^d(t+1)=\begin{cases}p_i^d(t),if\,f(P_i(t))<f(X_i(t+1))\\x_i^d(t+1),else \end{cases}\tag{2}$
所有电荷的全局最佳适应度的位置用 $P_{best}=X_{best}$ 表示。

在时刻 $t$ ，电荷 $j$ 在第 $d$ 维上受到电荷 $i$ 的库仑力如公式所示：
$F_{ij}^d = K(t)\frac{Q_i(t)Q_j(t)(p_j^d(t)-X_i^d(t))}{R_{ij}(t)+\xi} \tag{3}$
式中， $Q_i$ 和 $Q_j$ 分别为作用电荷 $i$ 和被作用电荷 $j$ 的电荷量， $\xi$ 表示一个极小的常量。 $R_{ij}$ 为电荷 $i$ 与电荷 $j$ 之间的欧式距离，由下式得出：
$R_{ij}(t)=||X_i^t(t),X_j(t)||_2\tag{4}$
库仑常数 $K (t)$ 为当前迭代数和系统迭代数的函数，为了控制算法的搜索精度，呈指数递减。时刻 $t$ 的库仑常数 $K (t)$ 可由下式计算：
$K(t)=K_0*exp(-\alpha \frac{iter}{maxiter}) \tag{5}$
式中， $K_0$ 为初值， $\alpha$ 为常数， $i t e r$ 是当前迭代的次数， $m a x i t e r$ 是系统迭代的次数。在算法伊始，库仑常数常被初始化为一个较高的值，以便于算法的初期探索，然后逐次迭代递减以控制算法的搜索精度。

在 $d$ 维搜索空间中，作用在电荷 $i$ 总的作用力等于来自其他所有电荷作用力的总和，其大小为：
$F_i^d(t)=\sum_{j=1,j\ne i}^N rand_jF_{ij}^d(t)\tag{6}$
式中， $rand_j$ 为[0,1]之间的随机数， $F_{ij}^d$ 为电荷 $j$ 作用在电荷 $i$ 上的库仑力。在任意时刻 $t$ ，电荷 $i$ 位于第 $d$ 维时的电场强度由下式给出：
$E_i^d(t)=\frac{F_i^d(t)}{Q_j(t)}\tag{7}$
根据牛顿第二运动定律，得出电荷 $i$ 在时刻 $t$ 的加速度：
$a_i^d(t)=Q_i(t)E_i^d(t)/M_i(t)\tag{8}$
在每一次迭代过程中，电荷 $i$ 根据计算得到的加速度来更新电荷的速度和位置，更新方式如下式所示：
$V_i^d(t+1)=rand_i*V_i^d(t) + a_i^d(t) \tag{9}$

$X_i^d(t+1)=X_i^d(t)+V_i^d(t+1)\tag{10}$

式中， $V_i^d(t),X_i^d(t)$ 分别为电荷 $i$ 在时刻 $t$ 的速度和位置。

电荷的电荷量通过适应度函数计算得出，假设初始时每个电荷的电荷量相等。
$Q_i(t)=Q_j(t);i,j=1,2,...,N \tag{11}$
采取以下式子更新电荷 $i$ 的电荷量：
$q_i(t)=exp(\frac{fit_{pi}(t) -worst(t)}{best(t)-worst(t)})\tag{12}$