基于动态局部进化的高斯粒子群优化算法

1.本发明涉及粒子群算法技术领域，具体为基于动态局部进化的高斯粒子群优化算法。


背景技术：

2.在过去的十年中，进化计算(ec)开始在不同的领域受到极大的关注。大量的研究人员从事进化计算的理论和应用研究。不同类型的进化算法(ea)，包括遗传算法(ga)、差分进化(de)等，已经针对不同的需求被开发出来。作为进化论的一个重要分支，群集智能(si)在优化领域表现突出。
3.大多数si算法的灵感来自于自然界的生物，如蚁群优化(aco)、粒子群优化(pso)等。pso的灵感来自于鸟群和鱼群。粒子在pso中被用来模拟鱼或鸟的团队行为。pso算法在经济调度、电力系统和其他许多领域表现出显著的性能。然而随着技术的进步，传统pso算法不能很好地解决日益复杂的优化问题。同时，pso算法的性能大大依赖于各个参数的设定，在面对全新问题是很容易陷入局部最优，为此，我们提出基于动态局部进化的高斯粒子群优化算法。


技术实现要素：

4.本发明的目的在于提供基于动态局部进化的高斯粒子群优化算法，以解决上述背景技术中提出的问题。
5.为实现上述目的，本发明提供如下技术方案：基于动态局部进化的高斯粒子群优化算法，具体的步骤为下：
6.s1、对标准粒子群算法进行优化，选择粒子群优化算法的宏观方向参数，表示粒子构造及种群规模，加快算法寻优的速度，避免算法陷入局部极值；
7.s2、设置粒子群优化算法的微观方向参数，根据社会搜索和认知搜索更新粒子自身运动，对粒子群算法的权重参数按照构建的非线性函数进行更新优化，以此加快收敛的速度，并且保证较高的精度；
8.s3、对粒子群算法的两个学习因子更新公式进行改进，保证粒子在运算前期具有较好的全局搜索能力，同时保证粒子在运算后期有较强的局部收敛能力，使收敛速度加快，结合当前参数，使用粒子群优化算法进行训练学习，以使粒子群中的粒子不断靠拢最优粒子；
9.s4、对粒子群算法的速度更新公式进行改进，使得如果粒子在上一次迭代中朝着全局最优方向飞行，那么在下一次迭代中的粒子飞行速度不变，否则按照速度更新公式进行更新速度；
10.s5、通过优化后的粒子群算法中粒子个体之间的协作，最终达到群体最优，通过个体与群体最佳位置的多次迭代，通过训练完成的动态局部进化模型对渲染数据进行时间预估，制定集群作业调度策略，以达到减少真实感渲染时间的目的；
11.s6、利用粒子群算法对所述函数模型进行全局搜索，同时，运行辅助算法对粒子群算法中的粒子以及粒子的个体最优位置pbest进行选择性更新，当所述粒子群算法达到最大的迭代次数时，根据最终代粒子群搜索到的最优解的相应变量值对所述待优化目标问题的参数进行对应设置，以获得所述待优化目标问题的解决方案。
12.优选的，所述s4中用构成径向基函数神经网络所需要的参数构成改进粒子群算法中的粒子，将构成的粒子代入目标函数中，得到初始适应度值。
13.优选的，所述s5中通过优化后的粒子群算法进行粒子位置及速度的更新，并计算更新之后的适应度值。
14.优选的，所述s6中判断适应度值是否达到设定的目标值，或者迭代次数到达最大值，满足其中一项条件即可停止迭代，取当前迭代出来的粒子作为最佳粒子，用构成最佳粒子的参数作为最优径向基函数神经网络的参数；若不满足条件，则继续执行s2，直到满足结束迭代条件得出最优径向基函数神经网络参数。
15.优选的，所述用于根据粒子群算法的粒子位置更新公式进行计算得到每个粒子在当前运行代数下的新位置矢量，对所述新位置矢量的目标函数值进行计算，并根据目标函数值的优劣以对每个粒子的个体最优位置和种群最优位置进行更新；
16.优选的，所述设ppbest表示当前粒子本身搜索到的最佳方案，ppbest＝(ppbest1,ppbest2,...,ppbestn)，其中，ppbestm表示当前粒子本身在第m维搜索到的最佳方案，m∈[1，n]，n为单个粒子的维度，粒子自身运动更新为：
[0017]
vi(t 1)＝ωvi(t) c1r1(pi(t) rg1gi(t)-xi(t)) c2r2(pg(t)-xi(t))
[0018]
xi(t 1)＝xi(t) vi(t 1)
[0019]
优选的，所述s4中利用辅助算法对第t代粒子以及第t代之前的h-1代粒子进行聚类并生成新个体，得到第t 1代备选粒子，并计算粒子的目标函数值，用于根据预设的第一概率在所述第t 1代备选粒子中选择出替换粒子对第t 1代粒子进行替换，以对所述第t 1代粒子进行更新，用于根据预设的第二概率在所述第t 1代备选粒子中选择出若干比较粒子，根据目标函数值依次判断所述比较粒子是否优于所述第t 1代粒子的个体最优位置pbest；
[0020]
优选的，所述vi(t 1)是粒子i第t 1次迭代的速度，xi(t 1)是粒子i第t 1次迭代的位置，pi(t)是粒子i第t次迭代的极值，pg(t)是第t次迭代粒子群全局最优解，t是当前时刻迭代次数，r1、r2是(0,1)之间的任一随机数，c1、c2是大于零的学习因子，ω是惯性权重，rg1是(0,1)之间的任一随机数，gi(t)是粒子i在第t次迭代时产生的高斯扰动，基于更新的粒子速度和位置继续在粒子群个体和蝙蝠个体中进行搜索，直至达到最大迭代次数，输出最优个体。
[0021]
优选的，所述惯性权重计算如下：
[0022]
其中，ωmin是最小惯性权重，ωmax是最大惯性权重，t是当前时刻迭代次数；tmax是最大迭代次数。
[0023]
高斯扰动计算如下：
[0024]
gi(t)＝rg2*gaussian(μ,σ2)
[0025]
其中，rg2是(0,1)之间的任一随机数，μ是均值，σ2是方差。
[0026]
与现有技术相比，本发明的有益效果是：
[0027]
1.该发明解决了传统粒子群算法搜索精度低，容易陷入局部最优的问题，为单目标优化问题提供高精度搜索结果。
具体实施方式
[0028]
下面将结合本发明实施例，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0029]
本发明提供一种技术方案：基于动态局部进化的高斯粒子群优化算法，具体的步骤为下：
[0030]
s1、对标准粒子群算法进行优化，选择粒子群优化算法的宏观方向参数，表示粒子构造及种群规模，加快算法寻优的速度，避免算法陷入局部极值；
[0031]
s2、设置粒子群优化算法的微观方向参数，根据社会搜索和认知搜索更新粒子自身运动，对粒子群算法的权重参数按照构建的非线性函数进行更新优化，以此加快收敛的速度，并且保证较高的精度；
[0032]
s3、对粒子群算法的两个学习因子更新公式进行改进，保证粒子在运算前期具有较好的全局搜索能力，同时保证粒子在运算后期有较强的局部收敛能力，使收敛速度加快，结合当前参数，使用粒子群优化算法进行训练学习，以使粒子群中的粒子不断靠拢最优粒子；
[0033]
s4、对粒子群算法的速度更新公式进行改进，使得如果粒子在上一次迭代中朝着全局最优方向飞行，那么在下一次迭代中的粒子飞行速度不变，否则按照速度更新公式进行更新速度；
[0034]
s5、通过优化后的粒子群算法中粒子个体之间的协作，最终达到群体最优，通过个体与群体最佳位置的多次迭代，通过训练完成的动态局部进化模型对渲染数据进行时间预估，制定集群作业调度策略，以达到减少真实感渲染时间的目的；
[0035]
s6、利用粒子群算法对所述函数模型进行全局搜索，同时，运行辅助算法对粒子群算法中的粒子以及粒子的个体最优位置pbest进行选择性更新，当所述粒子群算法达到最大的迭代次数时，根据最终代粒子群搜索到的最优解的相应变量值对所述待优化目标问题的参数进行对应设置，以获得所述待优化目标问题的解决方案。
[0036]
优选的，所述s4中用构成径向基函数神经网络所需要的参数构成改进粒子群算法中的粒子，将构成的粒子代入目标函数中，得到初始适应度值。
[0037]
优选的，所述s5中通过优化后的粒子群算法进行粒子位置及速度的更新，并计算更新之后的适应度值。
[0038]
优选的，所述s6中判断适应度值是否达到设定的目标值，或者迭代次数到达最大值，满足其中一项条件即可停止迭代，取当前迭代出来的粒子作为最佳粒子，用构成最佳粒子的参数作为最优径向基函数神经网络的参数；若不满足条件，则继续执行s2，直到满足结束迭代条件得出最优径向基函数神经网络参数。
[0039]
优选的，所述用于根据粒子群算法的粒子位置更新公式进行计算得到每个粒子在当前运行代数下的新位置矢量，对所述新位置矢量的目标函数值进行计算，并根据目标函
数值的优劣以对每个粒子的个体最优位置和种群最优位置进行更新；
[0040]
优选的，所述设ppbest表示当前粒子本身搜索到的最佳方案，ppbest＝(ppbest1,ppbest2,...,ppbestn)，其中，ppbestm表示当前粒子本身在第m维搜索到的最佳方案，m∈[1，n]，n为单个粒子的维度，粒子自身运动更新为：
[0041]
vi(t 1)＝ωvi(t) c1r1(pi(t) rg1gi(t)-xi(t)) c2r2(pg(t)-xi(t))
[0042]
xi(t 1)＝xi(t) vi(t 1)
[0043]
优选的，所述s4中利用辅助算法对第t代粒子以及第t代之前的h-1代粒子进行聚类并生成新个体，得到第t 1代备选粒子，并计算粒子的目标函数值，用于根据预设的第一概率在所述第t 1代备选粒子中选择出替换粒子对第t 1代粒子进行替换，以对所述第t 1代粒子进行更新，用于根据预设的第二概率在所述第t 1代备选粒子中选择出若干比较粒子，根据目标函数值依次判断所述比较粒子是否优于所述第t 1代粒子的个体最优位置pbest；
[0044]
优选的，所述vi(t 1)是粒子i第t 1次迭代的速度，xi(t 1)是粒子i第t 1次迭代的位置，pi(t)是粒子i第t次迭代的极值，pg(t)是第t次迭代粒子群全局最优解，t是当前时刻迭代次数，r1、r2是(0,1)之间的任一随机数，c1、c2是大于零的学习因子，ω是惯性权重，rg1是(0,1)之间的任一随机数，gi(t)是粒子i在第t次迭代时产生的高斯扰动，基于更新的粒子速度和位置继续在粒子群个体和蝙蝠个体中进行搜索，直至达到最大迭代次数，输出最优个体。
[0045]
优选的，所述惯性权重计算如下：
[0046]
其中，ωmin是最小惯性权重，ωmax是最大惯性权重，t是当前时刻迭代次数；tmax是最大迭代次数。
[0047]
高斯扰动计算如下：
[0048]
gi(t)＝rg2*gaussian(μ,σ2)
[0049]
其中，rg2是(0,1)之间的任一随机数，μ是均值，σ2是方差。
[0050]
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。