一种基于改进引力搜索算法的多目标优化方法技术

本发明专利技术公开了一种基于改进引力搜索算法的多目标优化方法，其算法将记忆性策略引入万有引力搜索算法，使粒子不仅将群体信息共享，也实现了上下代之间信息的共享，从而平衡粒子的全局与局部搜索能力，克服早熟收敛问题。在此基础上，该算法引入多样性增强机制，即在每次迭代中对粒子的速度和位置加以控制，以减缓多样性的流失，提高粒子的多样性，从而增强非支配解集的多样性与分布性，因此，基于改进引力搜索算法的多目标优化方法能够有效避免多目标优化陷入局部极值，显著改善引力搜索算法运用到多目标优化领域时非支配解的收敛性、多样性与分布性。

【技术实现步骤摘要】

本专利技术涉及。
技术介绍
多目标优化问题在科学研究与工程应用中都处于非常重要的地位。与单目标优化问题的本质区别在于，它的解并不是唯一的，而是存在一组由众多Pareto最优解组成的最优解集合，集合中的各个元素称为Pareto最优解或非支配解。由于在多目标优化算法中，全局最优粒子的计算与搜索(引导粒子的选取)，对多目标优化中解的收敛性与分布性具有重要影响，具备突出的全局寻优能力的进化算法被应用到多目标优化领域中。目前发展比较成熟的多目标进化算法主要包括基于遗传算法GA的多目标优化和基于粒子群算法PSO的多目标优化。在遗传算法方面，主要有非支配排序遗传算法NSGA及其改进算法NSGA-11、小生境非劣解排序遗传算法NPGA，以及加强型非劣解排序遗传算法SPEA2等，其中尤以NSGA-1I搜索寻优效果和收敛速度最为突出。而在粒子群方面，主要就是基于外部存档的粒子群多目标优化算法。由于粒子群算法具备两个特殊的性质:(1)记忆性——用来存储全局最优粒子与个体历史最优值；(2)信息交流性——粒子之间依据记忆特性相互分享最优位置的信息，使得粒子群在多目标优化领域表现出一定的实用性。作为一种新型进化算法，引力搜索算法已经成功应用到单目标优化领域，并在多目标优化领域取得了一定的效果。其根本思想是基于牛顿的万有引力定律:“在宇宙间，每一个粒子由于万有引力的作用而彼此相互吸引，引力的大小与粒子的质量成正比，与他们之间的距离成反比”。所以，通过粒子间的相互吸引，引力搜索算法保证了所有粒子向着质量最大的粒子移动。但是，当引力搜索算法运用到多目标优化时，其自身的一些缺点导致该算法中全局最优粒子的质量较低，多目标优化的效果还有待提高。首先，在引力搜索算法中，只有当前的位置信息在迭代更新过程中起作用，即该算法是一种缺乏记忆性的算法，这就导致种群上下代之间没有信息交流，容易陷入早熟收敛。另一方面，由于引力搜索算法中粒子速度较大，全部都向质量较大的粒子移动，收敛非常迅速，所以种群的多样性降低快速，即多样性迅速流失，不能保证非支配解的多样性与分布性。
技术实现思路
本专利技术要解决的技术问题是提供，它能够有效避免多目标优化陷入局部极值，显著改善引力搜索算法运用到多目标优化领域时非支配解的收敛性、多样性与分布性。为解决上述技术问题，本专利技术采用如下技术方案:包括以下步骤:I)明确要解决的多目标优化问题MOP:该问题的目标函数，整个问题的搜索空间维度n，空间上下限范围ub、lb，以及目标个数m ；2)种群初始化:设定引力场中粒子数目N，最大迭代次数M ;随机初始化种群中各个粒子的位置X、速度V、速度的上下界范围从Vmax到Vmin之间；每个粒子个体历史最优位置Pbesti赋值:将初始每个Xi的位置赋予对应Pbesti ;设置权重系数w的范围；初始化协调因子:变系数C1, C2;3)初始化外部档案NP与全局最优值gbest:分别计算初始X中每个粒子在各个目标下的函数值Π，f2…fm，根据支配关系将初始X划分为两个子种群A和B，分别存储非支配子集A，支配子集B，NP=A,从NP中随机选取一个粒子作为gbest ；4)根据目标问题，计算每个粒子的适应值:首先计算每个粒子在每个目标的函数值，然后按照线性权重的方式，将多个目标的适应度值，规划为一个适应度值:本文档来自技高网...

【技术保护点】
一种基于改进引力搜索算法的多目标优化方法，其特征在于：包括以下步骤：1）明确要解决的多目标优化问题MOP：该问题的目标函数，整个问题的搜索空间维度n，空间上下限范围ub、lb，以及目标个数m；2）种群初始化：设定引力场中粒子数目N，最大迭代次数M；随机初始化种群中各个粒子的位置X、速度V、速度的上下界范围从Vmax到Vmin之间；每个粒子个体历史最优位置pbesti赋值：将初始每个Xi的位置赋予对应pbesti；设置权重系数w的范围；初始化协调因子：变系数c1,c2；3）初始化外部档案NP与全局最优值gbest：分别计算初始X中每个粒子在各个目标下的函数值f1,f2…fm，根据支配关系将初始X划分为两个子种群A和B，分别存储非支配子集A,支配子集B，NP=A，从NP中随机选取一个粒子作为gbest；4）根据目标问题，计算每个粒子的适应值：首先计算每个粒子在每个目标的函数值，然后按照线性权重的方式，将多个目标的适应度值，规划为一个适应度值：fitness＝f1*(1?t/M)+f2*t/M其中，f1与f2表示两个目标问题的目标函数值，t表示当前迭代次数，M表示最大迭代次数；5）计算引力加速度：根据引力搜索算法的介绍，加速度的计算公式为：aijd(t)=Fijd(t)Mi(t),其中，是两个粒子在d维度上的相互作用力：Fijd(t)=G(t)Mi(t)×Mj(t)Rij(t)+ϵ(xjd(t)-xid(t)),Mi(t)是粒子i的质量：Mi(t)=mi(t)Σj=1Nmj(t),而粒子质量m的计算，是根据适应度值，通过归一化实现的：mi(t)=fiti(t)-worst(t)best(t)-worst(t),对于不同的多目标问题，best与worst的标准不同，若求最小值，则：best(t)=minj∈{1,...,N}fitj(t),若求最大值，则：best(t)=maxj∈{1,...,N}fitj(t),6）更新粒子的速度与位置：在这个更新过程中，粒子个体最优和全局最优被引入，修正了粒子的速度更新公式；同时，引入协调因子与权重系数的概念，用于调整所记忆的历史信息在优化过程中的影响比例，以控制粒子的全局与局部搜索的平衡，所以粒子i在t时刻的速度和位置的进化公式为：vid(t+1)=w·vid(t)+c1·r1(t)·(pbestid(t)-xid(t))+c2·r2(t)·(gbestid(t)-xid(t))+aid(t)xid(t+1)=xid(t)+vid(t+1)其中，r1,r2是从0到1间均匀分布的相互独立的随机数序列，c1,c2作为协调因子，设定为数值不断变化的系数，w是变权重系数；在粒子速度与位置更新过程中，若出现速度或位置超出搜索限制的情况，需要加入多样性增强机制来控制粒子的速度与位置；7）更新外部存档：由于外部存档用来存放非支配解子集，所以需要将每一代种群产生的新的粒子与NP内的粒子进行比较，保留非支配解，删除支配解；另外，在每一次更新过程中，都根据密度评估技术、分散度技术等删除最劣解，其中解的密度越低，分散度越好，表明解越占优；8）更新全局最优与个体最优粒子：全局最优gbest，是从外部存档NP中随机选择的，而每个粒子的个体最优pbest是通过判断新的粒子与个体历史最优pbest的支配关系来更新的，如果新的粒子被pbest支配，则pbest的值不变；如果新的粒子支配pbest，则pbest被替换；9）不断搜索迭代，直至达到迭代次数或满足要求精度为止。FDA0000440962460000012.jpg...

【技术特征摘要】
1.一种基于改进引力搜索算法的多目标优化方法，其特征在于:包括以下步骤: O明确要解决的多目标优化问题MOP:该问题的目标函数，整个问题的搜索空间维度η,空间上下限范围ub、lb,以及目标个数m ； 2)种群初始化:设定引力场中粒子数目N，最大迭代次数M;随机初始化种群中各个粒子的位置X、速度V、速度的上下界范围从Vmax到Vmin之间；每个粒子个体历史最优位置Pbesti赋值:将初始每个Xi的位置赋予对应Pbesti ;设置权重系数w的范围；初始化协调因子:变系数C1, C2 ； 3)初始化外部档案NP与全局最优值gbest:分别计算初始X中每个粒子在各个目标下的函数值Π，f2…fm，根据支配关系将初始X划分为两个子种群A和B，分别存储非支配子集A，支配子集B，NP=A,从NP中随机选取一个粒子作为gbest ； 4)根据目标问题，计算每个粒子的适应值:首先计算每个粒子在每个目标的函数值，然后按照线性权重的方式，将多个目标的适应度值，规划为一个适应度值:fitness = fI*(1-t/M)+f2*t/M 其中，fl与f2表示两个目标问题的目标函数值，t表示当前迭代次数，M表示最...

【专利技术属性】
技术研发人员：孙根云，张爱竹，王振杰，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：