2025年KBS新算法 SCI1区TOP:长颖燕麦优化算法AOO,深度解析+性能实测
目录
- 1.摘要
- 2.算法原理
- 3.结果展示
- 4.参考文献
- 5.文章&代码获取
1.摘要
本文提出了一种新颖的元启发式算法——长颖燕麦优化算法(AOO),该算法灵感来自动画燕麦在环境中的自然行为。AOO模拟了长颖燕麦的三种独特行为:(i) 通过自然元素如风、水和动物进行种子传播;(ii) 在吸湿运动的影响下,长颖燕麦种子的主要芒刺发生变形和旋转,使整个种子能够滚动并传播;(iii) 在滚动传播过程中,当种子遇到障碍物时,会储存能量,并在特定条件下触发推进机制,进一步传播种子。
2.算法原理
参数计算
长颖燕麦种子传播过程的特征与种子主芒的长度、质量和滚动过程中的偏心系数有关。
{ m = 0.5 × r d i m L = N × r d i m e = 0.5 × r d i m c = 1 − ( t T ) 3 \left\{ \begin{array}{l} m = 0.5 \times \frac{r}{dim} \\ L = N \times \frac{r}{dim} \\ e = 0.5 \times \frac{r}{dim} \\ c = 1 - \left(\frac{t}{T}\right)^3 \end{array} \right. ⎩ ⎨ ⎧m=0.5×dimrL=N×dimre=0.5×dimrc=1−(Tt)3
其中, m m m表示长颖燕麦种子质量, L L L为长颖燕麦主芒长度, e e e为滚动种子时的偏心旋转系数。
探索阶段
在部分长颖燕麦从植物上脱落后,它们的传播主要是通过风、水或动物的影响或作用来实现的。这种分散模式表现出显著的随机性,允许在广泛的解空间中进行探索。
W = c π × ( 2 × r dim − 1 ) ⊗ U B W=\frac{c}{\pi}\times(2\times r_{\dim}-1)\otimes UB W=πc×(2×rdim−1)⊗UB
{ X t + 1 ( i ) = 1 N × ∑ i = 1 N X t ( i ) + W , if m o d ( i , N / 10 ) = 0 , X t + 1 ( i ) = X b e s t + W , if m o d ( i , N / 10 ) = 1 , X t + 1 ( i ) = X t ( i ) + W , else . \left\{ \begin{array}{l} X_{t+1}(i) = \frac{1}{N} \times \sum_{i=1}^{N} X_t(i) + W, \text{ if } \mod(i, N/10) = 0, \\ X_{t+1}(i) = X_{best} + W, \text{ if } \mod(i, N/10) = 1, \\ X_{t+1}(i) = X_t(i) + W, \text{ else}. \end{array} \right. ⎩ ⎨ ⎧Xt+1(i)=N1×∑i=1NXt(i)+W, if mod(i,N/10)=0,Xt+1(i)=Xbest+W, if mod(i,N/10)=1,Xt+1(i)=Xt(i)+W, else.
开发阶段
在此阶段,剩余的长颖燕麦种子根据是否遇到障碍物被分为两种传播方式,并假设这两种情况的概率相等。在没有障碍物的情况下,种子的位置变化通过湿度引起的应力梯度驱动吸湿滚动。受Lindtner等人研究的启发,他们表明纤维素微纤丝的取向决定了各向异性的膨胀,论文通过临界曲率模型来描述这种运动,采用了快折屈曲的方式。滚动机制通过偏心旋转和扭矩公式:
A = U B − ∣ U B × t × sin ( 2 × π × r ) T ∣ A=UB-\left|\frac{UB\times t\times\sin(2\times\pi\times r)}{T}\right| A=UB− TUB×t×sin(2×π×r)
R = ( m × e + L 2 ) × r dim ( − A , A ) dim R= \begin{pmatrix} m\times e+L^2 \end{pmatrix}\times\frac{r_{\dim}(-A,A)}{\dim} R=(m×e+L2)×dimrdim(−A,A)
X t ( i ) = X b e s t + R + c × L e v y ( dim ) ⊗ X b e s t X_t(i)=X_{best}+R+c\times Levy(\dim)\otimes X_{best} Xt(i)=Xbest+R+c×Levy(dim)⊗Xbest
当种子在传播过程中遇到障碍时,假设主芒在储存能量的驱动下进行抛射。整个种子传播过程以抛射运动为特征,AOO使用简单弹丸运动模型进行位置更新:
B = U B − ∣ U B × t × cos ( 2 × π × r ) T ∣ B=UB-\left|\frac{UB\times t\times\cos(2\times\pi\times r)}{T}\right| B=UB− TUB×t×cos(2×π×r)
{ k = 0.5 + 0.5 × r x = 3 × r d i m θ = π × r α = 1 π × e r T \left\{ \begin{array}{l} k = 0.5 + 0.5 \times r \\ x = 3 \times \frac{r}{dim} \\ \theta = \pi \times r \\ \alpha = \frac{1}{\pi} \times e^{\frac{r}{T}} \end{array} \right. ⎩ ⎨ ⎧k=0.5+0.5×rx=3×dimrθ=π×rα=π1×eTr
J = 2 × k × x 2 × sin ( 2 θ ) m g × r dim ( − B , B ) dim × ( 1 − α ) J=\frac{2\times k\times x^2\times\sin(2\theta)}{mg}\times\frac{r_{\dim}(-B,B)}{\dim}\times(1-\alpha) J=mg2×k×x2×sin(2θ)×dimrdim(−B,B)×(1−α)
X t ( i ) = X b e s t + J + c × L e v y ( dim ) ⊗ X b e s t X_t(i)=X_{best}+J+c\times Levy(\dim)\otimes X_{best} Xt(i)=Xbest+J+c×Levy(dim)⊗Xbest
流程图
伪代码
3.结果展示
4.参考文献
[1] Wang R B, Hu R B, Geng F D, et al. The Animated Oat Optimization Algorithm: A Nature-Inspired Metaheuristic for Engineering Optimization and a Case Study on Wireless Sensor Networks[J]. Knowledge-Based Systems, 2025: 113589.