一种基于增强局部搜索能力的人工蜂群优化方法技术

本发明专利技术的一种基于增强局部搜索能力的人工蜂群优化方法，通过采用一种遍历性好、分布均匀的高维混沌序列产生方式，在一定程度上避免了ABC中随机初始化的盲目性，采用基于对数的适应度评价方式，扩大了种群个体适应值的差异性，减小了选择压力，并采用一种增强局部搜索能力的搜索策略，在一定程度上有效地避免了陷入局部最优，从而提高了收敛速度，改善算法性能。

【技术实现步骤摘要】
一种基于增强局部搜索能力的人工蜂群优化方法
本专利技术涉及优化领域，尤其涉及一种基于增强局部搜索能力的人工蜂群优化方法。
技术介绍
优化是许多科学和工程领域中一个重要的问题，在生产与生活中，人们经常会遇到如何安排生产资源、销售计划，或者怎样合理选取参数、路径，才能使总利润最大、总成本最小、总时间最短，或使相应产品指标达到最优、资源利用率迖到最经济等诸如此类的问题，这些问题都可以归结为最优化问题。传统的优化方法具有对初值的强依赖性、容易陷入局部极值点、可能得不到全局最优值、一般要求解析的目标函数是连续且可导的，变量也是连续的、而对于一些复杂的非线性的、离散的或多目标值的优化问题，用传统的优化方法很难得到理想的结果等缺点。为了解决这些问题，研究者们开始将目光转向对大自然中生物所展示出的智能行为的模拟和研究，从而开启了一个重要的研究领域仿生智能计算来解决以往难以实现的优化问题。人工蜂群算法(ABC)是一种模拟蜜蜂群体寻找优良蜜源的仿生智能计算方法，与遗传算法、粒子群算法、蚁群算法等智能计算方法相比，该算法的突出优点是在每次迭代中都进行全局和局部搜索、控制参数少、易于实现、计算方便等。但是，人工蜂群算法仍存在一些问题，种群初始化盲目性、局部搜索能力弱和易陷入局部最优的缺点。因此，本专利技术人对其进一步的探索和研究，提出一种基于增强局部搜索能力的人工蜂群优化方法。
技术实现思路
本专利技术的技术目的在于提出一种基于增强局部搜索能力的人工蜂群优化方法，以增强局部搜索能力，与ABC的较强全局搜索能力平衡，提高收敛速度，改善算法性能。为了解决上述技术问题，本专利技术的技术方案如下：一种基于增强局部搜索能力的人工蜂群优化方法，包括以下步骤：步骤1、设置ESABC的参数，所述参数包括种群大小SN、采蜜蜂个数NE、跟随蜂个数NO、最大迭代次数MCN、个体维数D、阈值limit；步骤2、按公式(1)进行初始化，得到初始种群：其中，i＝1,2,…,SN,j＝1,2,…,D；SN表示食物源(解)的数量；D表示解的维数；由公式(2)所得：所述公式(2)为洛伦兹混沌系统，通过该公式产生三组分布均匀、遍历性好的混沌迭代序列；其中，取x(0),y(0),z(0)为初始值，δ，γ，β为洛伦兹系统的参数，取值分别为δ＝10，β＝8/3,γ＞24.74，从产生的三组混沌序列中随机选一维，记作步骤3、根据公式(3)分别计算各个解的适应值fit，并将适应值排名前SN/2的解作为初始的采蜜蜂种群，剩余的蜜蜂则为观察蜂种群：fiti＝0.1/(0.1+1/|lgfi|),0≤fi≤10-λ(3)其中，fiti表示第i个个体的适应值，fi表示第i个个体的函数值，λ则由计算机的计算精度决定；步骤4、在由步骤3选择出来的采蜜蜂附近进行邻域搜索，搜索策略如公式(4)所示：其中，r为{1,2,…,SN}内的随机整数，且r≠i，i和j的含义同公式(1)，为[-1,1]内的随机数，随机为第r个蜜蜂的第j维，代表当代种群中最优个体的第j维，即当前种群最优解引导下一代个体的进化方向，可以达到增强局部搜索能力和提高解精度的效果；步骤5、按照公式(3)计算其适应值fit_new，若fit<fit_new，更新当前采蜜蜂，trial(i)＝0,否则，trial(i)++；步骤6、对于跟随蜂，首先根据公式(5)计算选择概率P，以概率P(i)搜索新的食物源，然后转化为采蜜蜂进行邻域搜索，根据公式(3)计算适应值，比较并判断是否更新，同时修改标识向量trial(i)；其中概率计算如公式(5)所示：步骤7、判断trial(i)>limit是否成立，若成立，则放弃该食物源，进入侦查蜂阶段，根据公式(1)生成新的食物源，若不成立，转至步骤8；步骤8、记录最优解；步骤9、判断是否满足终止条件iter≥MCN，若满足，则输出最优解，否则，返回步骤4。采用上述方案后，本专利技术具有以下特点：一、初始化种群时采用一种遍历性好、分布均匀的高维混沌序列产生方式，增强了初始种群的多样性和分布性，在一定程度上避免了ABC中随机初始化的盲目性；二、采用基于对数的适应度评价方式划分初始的采蜜蜂种群和观察蜂种群，通过该方法把个体间差异明显化，扩大了种群个体适应值的差异性，减小了选择压力，进而把函数值相似但不同的种群个体区分开，使得优秀个体有更大的概率被跟随开采；三、本专利技术并采用一种增强局部搜索能力的搜索策略，在一定程度上有效地避免了陷入局部最优，提高了解的精度。下面结合附图和具体实施方式对本专利技术的技术方案进行详细说明。附图说明图1是本专利技术一种基于增强局部搜索能力的人工蜂群优化方法的流程简图；图2-12是本专利技术分别对于11个测试函数在D＝30的情况下的进化曲线图；图13-22是本专利技术分别对于11个测试函数在D＝60的情况下的进化曲线图。具体实施方式本专利技术公开的一种基于增强局部搜索能力的人工蜂群优化方法(ESABC)，能够达到搜索某一问题最优解的目的。下面是本专利技术的一个实施例，具体说明如何应用本专利技术ESABC。实施例：通过MATLAB(商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境)对12个典型的函数最小化问题的仿真寻优来测试本专利技术ESABC的性能，并与ABC(传统人工蜂群算法)、GABC(用于与本专利技术ESABC性能比较，参考ZhuG,KwongS.Gbest-guidedartificialbeecolonyalgorithmfornumericalfunctionoptimization[J].AppliedMathematics&Computation,2010,217(7):3166-3173.)进行比较。测试函数介绍如表1所示。表1.测试函数表1中F01-F04为单峰函数，F05-F12为多峰函数。ESABC，ABC和GABC被用于优化上述函数，三个算法在相同的实验背景下运行，且每个测试函数独立运行10次以避免偶然性，并记录最优值、最差值、平均值和方差。如图1所示，是采用本专利技术的一种基于增强局部搜索能力的人工蜂群优化方法对上述12个典型的函数最小化问题的仿真寻优过程，包括：S1、设置初始参数。本例设初始参数为：维数D＝30&60，种群个数SN＝150，最大迭代次数MCN＝1000，阈值limit＝100；S2、按公式(1)进行初始化，得到初始种群：其中，i＝1,2,…,SN,j＝1,2,…,D；SN表示食物源(解)的数量；D表示解的维数；由公式(2)所得：所述公式(2)为洛伦兹混沌系统，通过该公式产生三组分布均匀、遍历性好的混沌迭代序列，在一定程度上避免了ABC中随机初始化的盲目性；其中，取x(0),y(0),z(0)为初始值，δ，γ，β为洛伦兹系统的参数，取值分别为δ＝10，β＝8/3,γ＞24.74，从产生的三组混沌序列中随机选一维，记作S3、根据公式(3)分别计算各个解的适应值fit，即将步骤S2得到的初始个体向量(解)代入测试函数，即得到目标函数值fi，如ellipitic函数，目标函数值然后将fi代入公式(3)计算适应值fit，并将适应值排名前75的解作为初始的采蜜蜂种群，剩余的蜜蜂则为观察蜂种群：fiti＝0.1/(0.1+1/本文档来自技高网...

【技术保护点】
一种基于增强局部搜索能力的人工蜂群优化方法，其特征在于，包括以下步骤：步骤1、设置ESABC的参数，所述参数包括种群大小SN、采蜜蜂个数NE、跟随蜂个数NO、最大迭代次数MCN、个体维数D、阈值limit；步骤2、按公式(1)进行初始化，得到初始种群：

【技术特征摘要】
1.一种基于增强局部搜索能力的人工蜂群优化方法，其特征在于，包括以下步骤：步骤1、设置ESABC的参数，所述参数包括种群大小SN、采蜜蜂个数NE、跟随蜂个数NO、最大迭代次数MCN、个体维数D、阈值limit；步骤2、按公式(1)进行初始化，得到初始种群：其中，i＝1,2,…,SN,j＝1,2,…,D；SN表示食物源的数量；D表示解的维数；由公式(2)所得：所述公式(2)为洛伦兹混沌系统，通过该公式产生三组混沌迭代序列；其中，取x(0),y(0),z(0)为初始值，δ，γ，β为洛伦兹系统的参数，取值分别为δ＝10，β＝83,γ＞24.74，从产生的三组混沌序列中随机选一维，记作步骤3、根据公式(3)分别计算各个解的适应值fit，并将适应值排名前SN/2的解作为初始的采蜜蜂种群，剩余的蜜蜂则为观察蜂种群：fiti＝0.1/(0.1+1/|lgfi|),0≤fi≤10-λ(3)其中，fiti表示第i个个体的适应值，fi表示第i个个体的函数值，λ则由计算机的计算精度决定；步骤4、在由步骤3选择出来的采蜜蜂附...

【专利技术属性】
技术研发人员：柳培忠，刘晓芳，顾培婷，骆炎民，汪鸿翔，陈智，范宇凌，
申请(专利权)人：华侨大学，泉州市华工智能技术有限公司，
类型：发明
国别省市：福建,35