[发明专利]一种基于策略适应性差分进化的全局优化方法

摘要

一种基于策略适应性差分进化的全局优化方法，首先，计算出种群中各个体与当前种群中最优个体之间的距离，并分别根据距离和目标函数值对整个种群进行排名；然后，根据距离排名和目标函数值排名的平均误差值判断当前种群中个体的分布情况，进而判断算法搜索所处的状态，即全局探测和局部搜索，同时对每种搜索状态设置多个不同的变异策略；最后，对于种群中的每个个体从其对应的状态策略池中随机选择一个变异策略来产生新个体，从而达到平衡算法全局探测能力和局部增强能力的效果。本发明有效避免策略选择不当而影响算法的性能、提升优化性能。

主权项

一种基于策略适应性差分进化的全局优化方法，其特征在于：所述优化方法包括以下步骤：1)初始化：设置种群规模N_P，初始交叉概率C_R，初始增益常数F；2)随机生成初始种群P＝{x^1,g,x^2,g,...,x^Np,g}，并计算出各个体的目标函数值，其中，g为进化代数，x^i,g,i＝1,2,…,Np表示第g代种群中的第i个个体，若g＝0，则表示初始种群；3)根据各个体x^i,g的目标函数值f(x^i,g)对各个体进行降序排列，并记下各个体的排名F^i,g，并找出当前种群中的最优个体x^best,g，其中，F^i,g表示第g代种群中第i个个体的目标函数值排名；4)根据公式(1)计算出初始种群中各个体与最优个体x^best,g之间的距离d^i,g；$d^{i,g}=\left(\underset{i=1}{\overset{N_p}{\Σ}}\underset{k=i+1}{\overset{N_P}{\Σ}}\sqrt{\underset{j=1}{\overset{N}{\Σ}}{\left(x_j^{i,g}-x_j^{best,g}\right)}^2}\right)/\left(N_p\right(N_p-1)/2)---\left(1\right)$其中，d^i,g表示第g代种群中第i个个体与最优个体x^best,g之间的距离，表示第g代种群中第i个个体x^i,g的第j维元素，表示第g代种群中第最优个体x^best,g的第j维元素，N为问题维数，N_P为种群规模；5)根据各个体与最优个体之间的距离d^i,g进行降序排列，并记下各个体的排名D^i,g，D^i,g表示第g代种群中第i个个体的距离排名；6)根据式(2)计算每一代中目标值排名和距离排名的平均误差值E^g：$E^g=\frac{1}{N_p}\underset{i=1}{\overset{N_p}{\Σ}}\left|F^{i,g}-D^{i,g}\right|---\left(2\right)$其中E^g表示第g代种群的平均误差值；7)根据公式(3)将平均误差值E^g进行归一化处理：${\stackrel{\‾}{E}}^g=\frac{E^g-E_{min}^g}{E_{max}^g-E_{min}^g}---\left(3\right)$其中，表示平均误差值E^g的归一化值，平均误差值表示E^g的最小值，取值始终为0，为E^g的最大值，当种群规模N_p为偶数时，当种群规模N_p为奇数时，$E_{max}^g=({N_p}^2-1)/2N_p;$8)判断进化过程所处的状态，对种群中的每个个体随机选择变异策略进行变异：8.1)如果则算法处于全局探测阶段，根据式(4)进行变异：$v_j^{i,g}=\left\{\begin{array}{cc}{x}_j^{pbest,g}+F\·(x_j^{a,g}-x_j^{b,g}),& ifrandn(1,3)=1\\ x_j^{a,g}+F\·(x_j^{pbest,g}-x_j^{a,g})+F\·(x_j^{b,g}-x_j^{c,g}),& ifrandn(1,3)=2\\ x_j^{i,g}+F\·(x_j^{pbest,g}-x_j^{i,g})+F\·(x_j^{a,g}-x_j^{b,g}),& otherwise\end{array}\right.---\left(4\right)$8.2)如果则算法处在局部搜索阶段，则根据式(5)进行变异：$v_j^{i,g}=\left\{\begin{array}{cc}{x}_j^{best,g}+F\·(x_j^{a,g}-x_j^{b,g}),& ifrandn(1,3)=1\\ x_j^{a,g}+F\·(x_j^{best,g}-x_j^{a,g})+F\·(x_j^{b,g}-x_j^{c,g}),& ifrandn(1,3)=2\\ x_j^{i,g}+F\·(x_j^{best,g}-x_j^{i,g})+F\·(x_j^{a,g}-x_j^{b,g}),& otherwise\end{array}\right.---\left(5\right)$在8.1)和8.2)中，rand(0,1)表示在区间[0,1]之间随机生成一个小数，j＝1,2,…,N，N为问题维数，g为进化代数，randn(1,3)表示在区间[1,3]之间随机生成一个整数，a,b,c∈{1,2,...,N_p}，a≠b≠c≠i，i为当前目标个体的索引，为第g代种群中第i个目标个体的变异个体的第j维元素，分别为第g代种群中第a、b、c个个体的第j维元素，为从0.5N_prandb(0,1)个个体中随机选取的最优个体的第j维元素，randb(0,1)表示随机产生0到1之间的小数，为当前第g代种群中的最优个体的第j维元素，F表示增益常数；9)根据公式(6)对每个变异个体进行交叉生成新个体trial^i,g：${\mathrm{trial}}_j^{i,g}=\left\{\begin{array}{cc}{v}_j^{i,g}& \begin{array}{ccc}if\left(randb\right(0,1)\≤C_R& or& j=rnbr\left(j\right)\end{array}\\ x_j^{i,g}& otherwise\end{array}\right.---\left(6\right)$其中，j＝1,2,…,N，表示第g代种群中第i个目标个体对应的新个体trial^i,g的第j维元素，randb(0,1)表示为随机产生0到1之间的小数，rnbr(j)表示随机产生1到N之间的整数，C_R表示交叉概率；10)根据公式(7)对每个新个体进行种群更新：$x^{i,g+1}=\left\{\begin{array}{cc}{\mathrm{trial}}^{i,g},& iff\left({\mathrm{trial}}^{i,g}\right)\≤f\left(x^{i,g}\right)\\ x^{i,g},& otherwise\end{array}\right.---\left(7\right)$其中，${\mathrm{trial}}^{i,g}=({\mathrm{trial}}_1^{i,g},{\mathrm{trial}}_2^{i,g},...,{\mathrm{trial}}_N^{i,g}),$$x^{i,g+1}=(x_1^{i,g+1},x_2^{i,g+1},...,x_N^{i,g+1}),$公式(7)表明，如果新个体优于目标个体，则新个体替换目标个体，否则保持目标个体不变；11)判断是否满足终止条件，如果满足，则保存结果并退出，否则返回步骤3)。