智能优化算法：原子轨道搜索算法-附代码

智能优化算法：原子轨道搜索算法

摘要：原子轨道搜索算法（Atomic orbital search, AOS），是于2021年提出的一种新型智能优化算法，该算法主要是基于量子力学的一些原理以及量子的原子模型行为，来进行寻优，具有寻优能力强，收敛速度快等特点。

1.算法原理

原子轨道搜索算法中，每个候选解的初始化，采用在搜索空间内随机初始化的方式。
$$x_i^j(0)=x_{i,\min }^j+\mathrm{rand}\cdot \left(x_{i,\max }^j-x_{i,\min }^j\right),\{\begin{array}{l}i=1,2,\cdots ,m\\ j=1,2,\cdots ,d\end{array}\text{\hspace{1em}(1)}$$
其中，$x_i^j(0)$表示候选解的初始位置；$x_{i,\min }^j$和$x_{i,\max }^j$是第$i$个候选解的第$j$个决策变量的最小和最大界限；$rand$是[0,1]范围内的均匀分布随机向量。

根据提供的基于量子的原子模型的详细信息，每个电子都有一个能量状态，在数学模型中被视为候选解的目标函数值。具有更好目标函数值的候选解代表具有更低能级的电子，而具有更高能级的电子在具有更差目标函数值的候选解的数学模型中被考虑。向量方程用于包含不同候选解(电子)的目标函数值(能级)，如下所示：
$$E=\left[\begin{array}{c}{E}_1\\ E_2\\ ⋮\\ E_i\\ ⋮\\ E_m\end{array}\right],i=1,2,\cdots ,m\text{\hspace{1em}(2)}$$
其中，$E$是目标函数值的向量，$E_i$是第$i$个候选解的能级，$m$是搜索空间(电子云)内候选解(电子)的数量。

因此每个候选解，以及每个候选解对应的能量密度（即适应度值），可以用如下表达式表达：
$$X^k=\left[\begin{array}{c}{X}_1^k\\ X_2^k\\ ⋮\\ X_i^k\\ ⋮\\ X_p^k\end{array}\right]=\left[\begin{array}{cccccc}{x}_1^1& x_1^2& \cdots & x_1^j& \cdots & x_1^d\\ x_2^1& x_2^2& \cdots & x_2^j& \cdots & x_2^d\\ ⋮& ⋮& & ⋮& \ddots & ⋮\\ x_i^1& x_i^2& \cdots & x_i^j& \cdots & x_i^d\\ ⋮& ⋮& & ⋮& \ddots & ⋮\\ x_p^1& x_p^2& \cdots & x_p^j& \cdots & x_p^d\end{array}\right],\{\begin{array}{l}i=1,2,\cdots ,p\\ j=1,2,\cdots ,d\\ k=1,2,\cdots ,n\end{array}\text{\hspace{1em}(3)}$$

$$E^k=\left[\begin{array}{c}{E}_1^k\\ E_2^k\\ ⋮\\ E_i^k\\ ⋮\\ E_p^k\end{array}\right],\{\begin{array}{l}i=1,2,\cdots ,p\\ k=1,2,\cdots ,n\end{array}\text{\hspace{1em}(4)}$$

其中，$X_i^k$是第$k$个虚拟层中的第$i$个候选解，$n$是虚拟创建的最大层数，$p$是第$k$个虚拟层中的候选解总数，$d$是问题维度，$E_i^k$是第$k$个虚拟层中第$i$个候选解的目标函数值。

在每个虚层中具有最佳目标函数值的候选解被视为每个虚层中具有最低能级的电子($L{E}^k$)。此外，在所有候选解之间具有最佳目标函数值的候选解被认为是原子中具有最低能级($LE$)的电子。通过考虑所选择的层中所有候选解的位置和目标函数值的平均值，来确定每个所考虑的虚层中候选解的结合状态和结合能。该目的的数学方程如下：
$$B{S}^k=\frac{{\sum }_{i=1}^p{X}_i^k}{p},\{\begin{array}{l}i=1,2,\cdots ,p\\ k=1,2,\cdots ,n\end{array}\text{\hspace{1em}(5)}$$

$$BE=\frac{{\sum }_{i=1}^m{E}_i}{m},i=1,2,\cdots ,m\text{\hspace{1em}(6)}$$

其中， $BS$和$BE$是原子的结合态和结合能；$X_i$和$E_i$是第$i$个候选解在原子中的位置和目标函数值。

为了从数学上表示光子对原子核周围电子的作用，在 (0,1)范围内为每个电子生成一个均匀分布的随机数($\varphi$)。此外，光子速率($PR$)被确定为代表考虑光子对电子作用的概率的参数。如果每个电子的随机生成数($\varphi$)大于$PR$，光子对电子的作用是可能的，因此根据光子的发射和吸收来考虑原子核周围不同层之间的电子运动。在这方面，将每个假想层中的每个候选解($X_i^k$)的能级($E_i^k$)与该层的结合能($B{E}^k$)进行比较，以决定光子的发射和吸收。如果特定层中的候选解高于该层的结合能($E_i^k\ge B{E}^k$)，考虑了光子的发射。在这个过程中，候选解倾向于发射具有$\beta$和$\gamma$能量的光子，以便同时达到原子的结合态($BE$)和原子中具有最低能级($LE$)的电子态。该过程中候选解的位置更新过程的数学方程如下：
$$X_{i+1}^k=X_i^k+\frac{{\alpha }_i\times \left({\beta }_i\times LE-{\gamma }_i\times BS\right)}{k},\{\begin{array}{l}i=1,2,\cdots ,p\\ k=1,2,\cdots ,n\end{array}\text{\hspace{1em}(7)}$$
其中，$X_i^k$和$X_{i+1}^k$是第$k$层第$i$个候选解的当前和未来位置；$LE$是原子中能级最低的候选解；$BS$是原子的结合态；${\alpha }_i,{\beta }_i,{\gamma }_i$ 是包含随机生成的数的向量，这些数均匀分布在 (0,1)中，用于确定释放的能量。

如果特定层中候选解的能级低于该层的结合能($E_i^k<B{E}^k$)，则考虑光子的吸收。在此过程中，候选解倾向于吸收具有$\beta$和$\gamma$能量的光子，以便同时达到层的结合态($B{S}^k$)和所考虑层内具有最低能级($L{E}^k$)的电子态。该过程中候选解的位置更新过程的数学方程如下：
$$X_{i+1}^k=X_i^k+{\alpha }_i\times \left({\beta }_i\times L{E}^k-{\gamma }_i\times B{S}^k\right),\{\begin{array}{l}i=1,2,\cdots ,p\\ k=1,2,\cdots ,n\end{array}\text{\hspace{1em}(8)}$$
其中，$X_i^k$和$X_{i+1}^k$是第$k$层第$i$个候选解的当前和未来位置；$L{E}^k$是第$k$层最低能级的候选解；$B{S}^k$是第$k$层的结合态；${\alpha }_i,{\beta }_i,{\gamma }_i$ 是包含随机生成的数的向量，这些数均匀分布在 (0,1)中，用于确定释放的能量。如果每个电子的随机生成数($\varphi$)小于$PR$($\varphi <PR$)，光子对电子的作用是不可能的，因此电子在原子核周围不同层之间的运动是基于一些其他作用来考虑的，例如与粒子或磁场的相互作用，这也会导致能量的吸收或发射。在这种情况下，基于这些影响的候选解的位置更新过程考虑如下：
$$X_{i+1}^k=X_i^k+r_i,\{\begin{array}{l}i=1,2,\cdots ,p\\ k=1,2,\cdots ,n\end{array}\text{\hspace{1em}(9)}$$
其中，$X_i^k$和$X_{i+1}^k$是第$k$层第$i$个候选解的当前和未来位置；$r_1$是包含随机生成的、均匀分布在(0,1)范围内的数的向量。

整个算法的流程图如下：

2.实验结果

3.参考文献

[1] Azizi M . Atomic Orbital Search: A Novel Metaheuristic Algorithm[J]. Applied Mathematical Modelling, 2021, 93(1).

4.Matlab