目录
- 1.UAV路径规划数学模型
- 2.改进点
- 3.结果展示
- 4.参考文献
- 5.代码获取
1.UAV路径规划数学模型
【智能算法应用】蜣螂优化算法DBO求解UAV路径规划
2.改进点
Logistics混沌映射
X n + 1 = μ X n ( 1 − X n ) , X n ∈ ( 0 , 1 ) (1) X_{_{n+1}} = \mu X_{_n}( 1 - X_{_n} ) ,\quad X_{_n} \in ( 0 ,1 )\tag{1} Xn+1=μXn(1−Xn),Xn∈(0,1)(1)
动态自适应权重
论文使用Gauss变异替换原有的更新方式,利用边界最大最小值两者作差与高斯分布算子进行交互形成高斯分布数。 这种变异不仅可以跳出局部最优并且可以优化算法的寻优能力,同时,少数麻雀种群会扩大搜索范围,增强种群多样性,有利于搜索潜在区域,从而加快算法的收敛速度。改进后的发现者位置更新公式为:
X i , j t + 1 = { [ X i , j t + ω ( f i , j t − X i , j t ) ] γ , R 2 < S T G ( ε ) g + X i , j t ( 1 − g ) , R 2 ⩾ S T (2) \boldsymbol{X}_{i,j}^{t+1} = \left\{\begin{array}{ll}\left[ \boldsymbol{X}_{i,j}^{t} + \boldsymbol{\omega}(\boldsymbol{f}_{i,j}^{t} - \boldsymbol{X}_{i,j}^{t} ) \right]\boldsymbol{\gamma} ,& R_{2} < \mathrm{ST}\\[2ex]\boldsymbol{G}( \boldsymbol{\varepsilon} ) \boldsymbol{g} + \boldsymbol{X}_{i,j}^{t} ( 1 - g ) ,& R_{2} \geqslant \mathrm{ST}\end{array}\right.\tag{2} Xi,jt+1=⎩ ⎨ ⎧[Xi,jt+ω(fi,jt−Xi,jt)]γ,G(ε)g+Xi,jt(1−g),R2<STR2⩾ST(2)
参数表述为:
ω = e 2 ( 1 − t / i t e r m a x ) − e − 2 ( 1 − t / i t e r m a x ) e 2 ( 1 − t / i t e r m a x ) + e − 2 ( 1 − t / i t e r m a x ) G ( ε ) = X i , j t + γ ( d 1 − d 2 ) g (3) \omega = \frac{\mathrm{e}^{2( 1-t/\mathrm{iter}_{\mathrm{max}})} - \mathrm{e}^{-2( 1-t/\mathrm{iter}_{\mathrm{max}})}}{\mathrm{e}^{2( 1-t/\mathrm{iter}_{\mathrm{max}})} + \mathrm{e}^{-2( 1-t/\mathrm{iter}_{\mathrm{max}})}}\\G( \varepsilon ) = X_{i,j}^{t} + \gamma( d_{1} - d_{2} ) g\tag{3} ω=e2(1−t/itermax)+e−2(1−t/itermax)e2(1−t/itermax)−e−2(1−t/itermax)G(ε)=Xi,jt+γ(d1−d2)g(3)
量子粒子群
量子粒子群算法能够在迭代过程中覆盖整个搜索空间,需要获得一个自适应的势能场,在分析粒子群演化时,粒子群在迭代过程中会收敛至一个点,这个点作为吸引子不断吸引粒子,并且有一定的概率能使得粒子有机会摆脱吸引子的束缚。本文中将采取量子粒子群生成新解的方式替换原有的追随者更新方式:
X i , j t + 1 = { P i , j t ± 0.5 L i , j t ln ( 1 μ i , j t + 1 ) , i > n 2 X p t + ∣ X i , j t − X p t ∣ A + L , i ⩽ n 2 (4) \boldsymbol{X}_{i,j}^{t+1} = \begin{cases}\boldsymbol{P}_{i,j}^{t} \pm 0. 5\boldsymbol{L}_{i,j}^{t}\ln\!\left(\frac{1}{\boldsymbol{\mu}_{i,j}^{t+1}}\right),&i > \frac{n}{2}\\\\\boldsymbol{X}_{p}^{t} + | X_{i,j}^{t} -\boldsymbol{X}_{p}^{t} |\boldsymbol{A}^{+} \boldsymbol{L},&i \leqslant \frac{n}{2}\end{cases}\tag{4} Xi,jt+1=⎩ ⎨ ⎧Pi,jt±0.5Li,jtln(μi,jt+11),Xpt+∣Xi,jt−Xpt∣A+L,i>2ni⩽2n(4)
其中,参数表述为:
P i , j t = φ i , j t X i , j t + ( 1 − φ i , j t ) G i t L i , j ′ = 2 a t ∣ C t − X i , j ι ∣ a t = a 0 + ( a 1 − a 0 ) ( T − t ) / T C ι = ∑ 1 M P i , j t (5) \begin{aligned} &\boldsymbol{P}_{i,j}^{t} = \boldsymbol{\varphi}_{i,j}^{t}\boldsymbol{X}_{i,j}^{t} + \left( 1 - \boldsymbol{\varphi}_{i,j}^{t} \right)\boldsymbol{G}_{i}^{t} \\ &L_{i,j}^{\prime} = 2a_{_t}\mid C_{_t} - X_{_{i,j}}^{\iota}\mid \\ &a_{_t} = a_{_0} + \left( a_{_1} - a_{_0} \right) ( T - t )/T \\ &C_{\iota} = \sum_{1}^{M}P_{i,j}^{t} \end{aligned}\tag{5} Pi,jt=φi,jtXi,jt+(1−φi,jt)GitLi,j′=2at∣Ct−Xi,jι∣at=a0+(a1−a0)(T−t)/TCι=1∑MPi,jt(5)
其中,Pti,j 为吸引子位置; Lti,j 为粒子当前位置与个体平均最优位置的距离; μti,j 均匀分布并且 μti,j ∈(0,1]; Gti 为 t 时刻麻雀种群得全局最优位置; φti,j
为一个随机数∈ (0,1]; Ct 为个体平均最好位置;a1和 a0 分别为收缩扩张因子取1和0. 5。
柯西变异
在种群的每一代更新结束之前将柯西变异算子扰动策略加入到麻雀算法:
X i , j t + 1 = X b e s t ( t ) + c a u c h y ( 0 , 1 ) X b e s t ( t ) (6) \boldsymbol{X}_{i,j}^{t+1} = \boldsymbol{X}_{\mathrm{best}}( t ) + \mathrm{cauchy}( 0 ,1 )\boldsymbol{X}_{\mathrm{best}}( t )\tag{6} Xi,jt+1=Xbest(t)+cauchy(0,1)Xbest(t)(6)
3.结果展示
CEC2005测试
无人机三维路径规划
4.参考文献
[1] 吴学礼,王超,赵俊棋,等.改进麻雀算法的无人机三维路径规划[J].科学技术与工程,2024,24(15):6534-6542.