智能优化算法：骑手优化算法-附代码

智能优化算法：骑手优化算法

摘要：骑手优化算法（Rider Optimization Algorithm,ROA），是于2019年提出的一种新型智能优化算法。该算法主要模拟了在比赛中骑手通过身份划分来引导整个团队走向胜利的过程。具有寻优能力强，收敛速度快等特点。

1.算法原理

在骑手优化算法中，骑手一共分为4类：

(a).绕行手：绕行手通过绕行引导路径到达目标。这意味着其不跟随领先骑手(种群最优个体)。

(b).跟随手：跟随手在大部分区域内依赖或跟随领先骑手。

©.超车手：这类骑手根据自身当前位置及领先骑手的位置来调整下一步行动。

(d).冲刺手：其根据其他骑手的位置以最大速度到达目标点。

1.1 种群位置初始化

将种群G均分为四组：绕行手(B)、跟随手(F)、超车手(O)、冲刺手(A)：
$$R=B+F+O+A\text{\hspace{1em}(1)}$$

$$B=F=O=A=R/4\text{\hspace{1em}(2)}$$

其中$R$为种群规模。

1.2 参数初始化

跟其他优化算法一样，在搜索空间内，随机生成种群T：
T t = { T i , j t } ; 1 ≤ i ≤ R ; 1 ≤ j ≤ Q (3) T_t=\{T_{i,j}^t\};1\leq i\leq R;1\leq j \leq Q \tag{3} Tt​={Ti,jt​};1≤i≤R;1≤j≤Q(3)

$$T_{i,j}=\{\begin{array}{l}{\theta }_i;ifj=1\\ T_{i,j-1}+\phi ;ifj\ne 1\&T_{i,j-1}+\phi \le 360\\ T_{i,j-1}+\phi -360;otherwise\end{array}\text{\hspace{1em}(4)}$$

$${\theta }_i=i\ast \frac{360^o}{R}\text{\hspace{1em}(5)}$$

$$\phi =\frac{360}{Q}\text{\hspace{1em}(6)}$$

其中Q为空间维度。自行车的齿轮E用式（7）表示：
E = { E i } ; 1 ≤ i ≤ R (7) E=\{E_i\};1\leq i \leq R \tag{7} E={Ei​};1≤i≤R(7)
初始时所有个体的齿轮参数均设为，迭代时其值从{0,1,2,3,4}中选取。

自行车的加度度e和制动器K：
e = { e i } ; 1 ≤ i ≤ R (8) e=\{e_i\};1\leq i \leq R \tag{8} e={ei​};1≤i≤R(8)

K = { K i } ; 1 ≤ i ≤ R (9) K=\{K_i\}; 1\leq i\leq R \tag{9} K={Ki​};1≤i≤R(9)

其值均处于[0,1]内，但前者初始为0，后者为1。

1.3 速度限制

该算法选择了一种动态速度边界：
$$V_{max}^i=\frac{X_U^i-X_L^i}{T_{OFF}}\text{\hspace{1em}(10)}$$
其中$X_U$、$X_L$分别为骑手$i$ 在此次迭代时的位置上、下限，$T_{OFF}$即为停止时间，对应最大迭代次数。因此，骑手$i$对速度的限制为;
$$V_i^E=\frac{V_{max}^i}{|E|}\text{\hspace{1em}(11)}$$
这里|E|为齿轮数,一般取5。

1.4 寻找成功率

成功率是根据骑手位置和目标位置之间的距离计算得出的：
$$r_i=\frac{1}{||X_i-L_T||}\text{\hspace{1em}(12)}$$
这里$X_i$为当前骑手i 的位置，LT为目标位置。

1.5 绕行手位置更新

绕行手绕过公共路径，因此位置更新是随机的：
$$X_{t+1}^B=\delta [X_t(\eta ,j)\ast \beta (j)+X_t(\xi ,j)\ast [1-\beta (j)]]\text{\hspace{1em}(13)}$$
其中 $\delta ,\beta (j)$均为[0,1]内一随机数，$\eta ,\xi$均为[1,R]内一随机整数，表示随机选择出的一个骑手。

1.6 跟随手位置更新

跟随手通过跟随领先者来更新他的位置：
$$X_{t+1}^F(i,k)=X^L(L,k)+cos(T_{i,k}^t\ast X^L(L,k)\ast d_i^t)\text{\hspace{1em}(14)}$$
式中$k$为当前空间维度，$X^L$为领导者$L$的位置，$d$ 的计算方式为：
$$d_i^t=v_i^t(1/T_{OFF})\text{\hspace{1em}(15)}$$

$$v_i^t=\frac{1}{3}[E_i^t\ast V_i^E+V_{max}^i\ast e_i^t+(1-K_i^t)\ast V_{max}^i]\text{\hspace{1em}(16)}$$

并不是每个跟随手的每个维度都要更新！是否更新取决于坐标选择器$P_{ON}^t$
$$P_{ON}^t=(\frac{t}{T_{OFF}})\ast Q\text{\hspace{1em}(17)}$$
这里t 为当前迭代次数。将其与相应的维度指数相乘，并将结果与空间维度Q进行比较，以确定位置变化需要考虑的坐标：
k < < { ( j ∗ P O N t ) ; i f   ( j ∗ P O N t ) < Q , ∀ j (18) k<<\{(j*P_{ON}^t);if\,(j*P_{ON}^t)<Q,\forall j \tag{18} k<<{(j∗PONt​);if(j∗PONt​)<Q,∀j(18)

1.7 超车手位置更新

超车手的位置更新如下式:
$$X_{t+1}^O(i,k)=X_t(i,k)+[D_t^I(i)\ast X^L(L,k)]\text{\hspace{1em}(19)}$$
中的k 为当前维度，D为当前个体的方向指示器:
$$D_t^I(i)=[\frac{2}{1-log(S_t^R(i))}]-1\text{\hspace{1em}(20)}$$
其中$S_t^R(i)$为相对成功率，超车手的相对成功率由其适应度值与种群内最差个体的适应度值之比确定：
$$S_t^R(i)=\frac{r_t(i)}{ma{x}_{r=i}^R{r}_t(i)}\text{\hspace{1em}(21)}$$
同跟随手，为了确定要更新的维度，首先要确定当前超车手与最优个体的间距：
$$l(i,j)=|X_t(i,j)-X^L(L,j)|\text{\hspace{1em}(22)}$$
之后要通过对比当前超车手各维度的平均位置$\mu$来确定:
k < < { j ; i f   l ( i , j ) < μ i   ∀ j (23) k<<\{j; if\,l(i,j)<\mu_i\, \forall j \tag{23} k<<{j;ifl(i,j)<μi​∀j(23)

1.8 冲刺手位置更新

试图超越领导者的冲刺手遵循与跟随手类似的更新过程。但是，冲刺手会更新坐标中的所有值，而不仅仅是选定维度中的值：
$$X_{t+1}^A(i,j)=X^L(L,j)+[cos(T_{i,j}^t)\ast X(L,j{)}^L]+d_i^t\text{\hspace{1em}(24)}$$

1.9 参数更新

活动计数A更新：
$$A_c^{t+1}(i)=\{\begin{array}{l}1;if{r}_{t+1}(i)>r_t(i)\\ 0;otherwise\end{array}\text{\hspace{1em}(25)}$$
转向角更新：
$$T_{i,j}^{t+1}=\{\begin{array}{l}{T}_{i+1,j}^t;if{A}_c^{t+1}(i)=1\\ T_{i-1,j}^t;if{A}_c^{t+1}(i)=0\end{array}\text{\hspace{1em}(26)}$$
齿轮更新：
$$E_i^{t+1}=\{\begin{array}{l}{E}_i^t+1;if{A}_c^{t+1}(i)=1\&E_i^t\ne |E|\\ E_i^t-1;if{A}_c^{t+1}=0\&E_i^t\ne 0\\ E_i^t;otherwise\end{array}\text{\hspace{1em}(27)}$$
加速度的更新方式为：
$$e_i^{t+1}=\frac{E_i^{t+1}}{|E|}\text{\hspace{1em}(28)}$$
制动器的更新方式为：
$$K_i^{t+1}=[1-\frac{E_i^{t+1}}{|E|}]\text{\hspace{1em}(29)}$$
算法伪代码如下：

1.Input: Random positions of theriders Xt
2.Output:Leading rider XL
3.Begin
4.	Initialize the population
5.	initialize the rider parameters: Steering angle T,Gear E,Accelerator e,and Brake K.
6.	Find the success rate rt
7.	while t<TOFF
8.		for i=1 to R
9.			update the positon of the bypass rider using Eq.(13)
10.			update the positon of the follower using Eq.(14)
11.			update the positon of the overtaker using Eq.(19)
12.			update the positon of the attacker using Eq.(24)
13.			Rank the riders based on rt
14.			Select the rider having the maximum rt as the leading rider
15.			update T, E, e,and K
16.			Return XL
17.			t=t+1
18.		end for
19.	end while
20.Terminate

2.实验结果

3.参考文献

[1] Binu D , Kariyappa B S . RideNN: A New Rider Optimization Algorithm-Based Neural Network for Fault Diagnosis in Analog Circuits[J]. Instrumentation and Measurement, IEEE Transactions on, 2019, 68(1):2-26.

4.Matlab代码