极光优化算法:图像分割和特征选择中的应用

1.摘要
本文介绍了一种基于极光现象的算法——极光优化算法(Polar Lights Optimizer, PLO)。极光是一种独特的自然奇观,当来自太阳风的高能粒子在地磁场和地球大气层的影响下汇聚在地球两极时,就会发生极光。PLO集成了旋转运动和极光椭圆步行,前者有助于局部开发,而后者则促进全局探索。通过这两种策略的协同结合,PLO实现了局部开发与全局探索的平衡。
2.算法原理
旋转运动
PLO受到高能粒子向地球旅行过程中旋转运动的启发。这些高能粒子,主要是电子和质子,从大约1.5亿公里外的太阳发射向地球。地球被其磁场所包围,这个磁场向外延伸约50,000至65,000公里。当这些粒子接近地球时,会受到地球磁场的阻力,并在磁场的影响下沿不同方向辐射。这些带电粒子在接近地球时会经历沿磁场线的旋转运动,这种运动可以通过洛伦兹力来解释:
F L = q ν B F_L=q\nu B FL=qνB
FL使带电粒子受到向心力的作用,导致带电粒子在磁场中沿磁力线旋转运动。带电粒子的方程可以用牛顿第二定律来描述:
F L = m d ν d t F_L=m\frac{d\nu}{dt} FL=mdtdν
因此,整理可得:
m d ν d t = q ν B m\frac{d\nu}{dt}=q\nu B mdtdν=qνB
在理想状态下,可以使用一个简单的方程来描述带电粒子的运动,但在实际大气中,空气分子的阻力会使得粒子的圆周运动变得不平滑,这主要是由于大气的阻尼作用。因此,文中引入了一个阻尼因子α来修改原有的速度方程:
m d ν d t = q ν B − α ν m\frac{d\nu}{dt}=q\nu B-\alpha\nu mdtdν=qνB−αν
求解一阶微分方程可得:
ν ( t ) = C e q B − a m t \nu(t)=Ce^{\frac{qB-a}mt} ν(t)=CemqB−at
极光椭圆步行

极光椭圆步行灵感来源于对极光的深入研究,特别是极光沿着名为极光椭圆的椭圆带形成的现象。极光椭圆的大小与行星际磁场的南北分量有关,其边界随地磁活动的变化而变化。极光椭圆步行的复杂波动将显著影响全局搜索,它是不可预测的混沌,满足了PLO对快速全局搜索解空间的需求。此外,Levy飞行(LF)经常用于元启发式算法中以增强全局探索,其本质上是一种随机非高斯游走:
L e v y ( d ) ∼ ∣ d ∣ − 1 − β , 0 < β ≤ 2 Levy(d)\sim|d|^{-1-\beta},0<\beta\leq2 Levy(d)∼∣d∣−1−β,0<β≤2
在极光椭圆步行中,LF模拟的高能粒子受地磁活动和大气的影响,极光椭圆边界在极方向上收缩,在赤道方向上膨胀:
A o = L e v y ( d ) × ( X a ν g ( j ) − X ( i , j ) ) + L B + r 1 × ( U B − L B ) / 2 X a ν g = 1 N × ∑ i = 1 N X ( i ) \begin{aligned}&Ao=Levy(d)\times\left(X_{a\nu g}(j)-X(i,j)\right)+LB+r_1\times(UB-LB)/2\\\\&X_{a\nu g}=\frac1N\times\sum_{i=1}^NX(i)\end{aligned} Ao=Levy(d)×(Xaνg(j)−X(i,j))+LB+r1×(UB−LB)/2Xaνg=N1×i=1∑NX(i)
PLO包含了两种主要搜索策略:旋转运动和极光椭圆步行。这两种运动模式分别代表不同的搜索策略和行为特征,共同构成了粒子从太阳到地球极点的飞行轨迹:
X n e w ( i , j ) = X ( i , j ) + r 2 × ( W 1 × ν ( t ) + W 2 × A o ) X_{new}(i,j)=X(i,j)+r_2\times(W_1\times\nu(t)+W_2\times Ao) Xnew(i,j)=X(i,j)+r2×(W1×ν(t)+W2×Ao)
其中,自适应权值W1和W2:
W 1 = 2 ( 1 + e − 2 ( t / T ) 4 ) − 1 W 2 = e − ( 2 t / T ) 3 \begin{aligned}&W_1=\frac2{(1+e^{-2(t/T)^4})}-1\\\\&W_2=e^{-(2t/T)^3}\end{aligned} W1=(1+e−2(t/T)4)2−1W2=e−(2t/T)3
粒子碰撞

粒子碰撞策略基于太阳风中的带电粒子(如电子和质子)在进入地球大气层时的相互碰撞现象。这些碰撞不仅引发了能量的转移和形式的转换,增强了极光的形成过程,还能导致粒子速度和方向的改变,从而增加了极光的亮度和复杂性。
在PLO算法中,粒子碰撞策略被用来帮助算法从局部最优解中跳出,通过模拟粒子间混乱的碰撞来创建新的解空间位置,这种策略使得粒子在解空间中可能以不可预测的方式移动。通过这种模拟,PLO算法能够有效地进行全局搜索,从而探索更广阔的解空间:
X n e w ( i , j ) = X ( i , j ) + sin ( r 3 × π ) × ( X ( i , j ) − X ( a , j ) ) , r 4 < K a n d r 5 < 0.05 X_{new}(i,j)=X(i,j)+\sin(r_3\times\pi)\times(X(i,j)-X(a,j)),r_4<Kandr_5<0.05 Xnew(i,j)=X(i,j)+sin(r3×π)×(X(i,j)−X(a,j)),r4<Kandr5<0.05
其中,参数w表述为:
K = ( t / T ) K=\sqrt{(t/T)} K=(t/T)
流程图

伪代码

3.结果展示


4.参考文献
[1] Yuan C, Zhao D, Heidari A A, et al. Polar Lights Optimizer: Algorithm and Applications in Image Segmentation and Feature Selection[J]. Neurocomputing, 2024: 128427.
5.代码获取
DAMO开发者矩阵,由阿里巴巴达摩院和中国互联网协会联合发起,致力于探讨最前沿的技术趋势与应用成果,搭建高质量的交流与分享平台,推动技术创新与产业应用链接,围绕“人工智能与新型计算”构建开放共享的开发者生态。
更多推荐

所有评论(0)