第12卷第5期 智能系统学报 Vol.12 No.5 2017年10月 CAAI Transactions on Intelligent Systems 0ct.2017 D0I:10.11992/is.201706014 网络出版地址：http:/kns.cmki.ne/kcms/detail/23.1538.TP.20170831.1058.006.html 一种精英反向学习的萤火虫优化算法 魏伟一，文雅宏 (西北师范大学计算机科学与工程学院，甘肃兰州730070) 摘要：为了提高传统莹火虫算法的收敛速度和求解精度，提出了一种精英反向学习的萤火虫优化算法。通过反向 学习策略构造精英群体，在精英群体构成的区间上求普通群体的反向解，增加了群体的多样性，提高了算法的收敛 速度：同时，为了避免最优个体陷入局部最优，使整个群体在搜索过程中出现停滞，提出了差分演化变异策略：最后， 提出了一种线性递减的自适应步长来平衡算法的开发能力。实验结果表明，算法在收敛速度和收敛精度上有更好 的效果。 关键词：萤火虫算法：精英反向学习：优化算法：精英群体：反向解：反向学习策略：差分演化变异：自适应步长 中图分类号：TP309.2文献标志码：A文章编号：1673-4785(2017)05-0710-07 中文引用格式：魏伟一，文雅宏.一种精英反向学习的萤火虫优化算法[J】.智能系统学报，2017,12(5)：710-716， 英文引用格式：WEI Weiyi,WEN Yahong.Firefly optimization algorithm utilizing elite opposition-based learning[J].CAAI transactions on intelligent systems,2017,12(5):710-716. Firefly optimization algorithm utilizing elite opposition-based learning WEI Weiyi,WEN Yahong (College of Computer Science and Engineering,Northwest Normal University,Lanzhou 730070,China) Abstract:To increase the convergence speed and solution accuracy of the traditional firefly algorithm,in this paper,we propose a firefly optimization algorithm that utilizes elite opposition-based learning.Using an opposition- based learning strategy,we constructed an elite group and,in the interval of the elite group,we solved the opposite solutions of the ordinary groups.This strategy could increase group diversity and improve the convergence speed of the algorithm.To prevent the optimal individual from falling into the local optimum,which could cause stagnation of the whole group during the search process,we introduce a differential evolutionary mutation strategy.Finally,we propose an adaptive step size with a linear decrease to balance the development ability of the algorithm. Experimental results show that the proposed algorithm can increase convergence speed and accuracy. Keywords:firefly algorithm;elite opposition-based learning;optimized algorithm;elite group;opposite solutions; opposition-based learning strategy;differential evolutionary mutation;adaptive step size 萤火虫算法(firefly algorithm,FA)是受自然界问题的优化上，萤火虫算法依然存在收敛速度慢、 中萤火虫发光特性的启发，由剑桥学者Yang于 解的精度不高、容易陷入局部最优等不足。近年 2008年提出的一种群体智能随机优化算法[1-1。在来，很多学者已经进行了多角度的改进。文献[10] 多个科学与工程领域中，萤火虫算法已得到成功的 为了解决萤火虫算法过早地收敛和陷入局部最优 应用[6)，虽然FA表现出了良好的性能，但在一些 的不足，利用广义反向学习策略来优化萤火虫算 法。文献[11]采用正交学习策略改进FA算法，利 收稿日期：2017-06-07.网络出版日期：2017-08-31. 基金项目：甘肃省科技计划资助项目(1506RZA130):甘肃省高等学校 用精英萤火虫来构造指导向量，通过指导向量引导 科研项目(2014B-018). 通信作者：文雅宏.E-mail:wwyahong@126.com. 群体向全局最优区域移动。文献「12]提出基于蛙
第 １２ 卷第 ５ 期 智 能 系 统 学 报 Ｖｏｌ．１２ №．５ ２０１７ 年 １０ 月 ＣＡＡＩ Ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ Ｉｎｔｅｌｌｉｇｅｎｔ Ｓｙｓｔｅｍｓ Ｏｃｔ． ２０１７ ＤＯＩ：１０．１１９９２ ／ ｔｉｓ．２０１７０６０１４ 网络出版地址：ｈｔｔｐ： ／ ／ ｋｎｓ．ｃｎｋｉ．ｎｅｔ ／ ｋｃｍｓ／ ｄｅｔａｉｌ ／ ２３．１５３８．ＴＰ．２０１７０８３１．１０５８．００６．ｈｔｍｌ 一种精英反向学习的萤火虫优化算法 魏伟一，文雅宏 （西北师范大学 计算机科学与工程学院，甘肃 兰州 ７３００７０） 摘 要：为了提高传统萤火虫算法的收敛速度和求解精度，提出了一种精英反向学习的萤火虫优化算法。 通过反向 学习策略构造精英群体，在精英群体构成的区间上求普通群体的反向解，增加了群体的多样性，提高了算法的收敛 速度；同时，为了避免最优个体陷入局部最优，使整个群体在搜索过程中出现停滞，提出了差分演化变异策略；最后， 提出了一种线性递减的自适应步长来平衡算法的开发能力。 实验结果表明，算法在收敛速度和收敛精度上有更好 的效果。 关键词：萤火虫算法；精英反向学习；优化算法；精英群体；反向解；反向学习策略；差分演化变异；自适应步长 中图分类号：ＴＰ３０９．２ 文献标志码：Ａ 文章编号：１６７３－４７８５（２０１７）０５－０７１０－０７ 中文引用格式：魏伟一，文雅宏．一种精英反向学习的萤火虫优化算法［Ｊ］． 智能系统学报， ２０１７， １２（５）： ７１０－７１６． 英文引用 格 式： ＷＥＩ Ｗｅｉｙｉ， ＷＥＮ Ｙａｈｏｎｇ． Ｆｉｒｅｆｌｙ ｏｐｔｉｍｉｚａｔｉｏｎ ａｌｇｏｒｉｔｈｍ ｕｔｉｌｉｚｉｎｇ ｅｌｉｔｅ ｏｐｐｏｓｉｔｉｏｎ⁃ｂａｓｅｄ ｌｅａｒｎｉｎｇ ［ Ｊ ］． ＣＡＡＩ ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ ｉｎｔｅｌｌｉｇｅｎｔ ｓｙｓｔｅｍｓ， ２０１７， １２（５）： ７１０－７１６． Ｆｉｒｅｆｌｙ ｏｐｔｉｍｉｚａｔｉｏｎ ａｌｇｏｒｉｔｈｍ ｕｔｉｌｉｚｉｎｇ ｅｌｉｔｅ ｏｐｐｏｓｉｔｉｏｎ⁃ｂａｓｅｄ ｌｅａｒｎｉｎｇ ＷＥＩ Ｗｅｉｙｉ， ＷＥＮ Ｙａｈｏｎｇ （Ｃｏｌｌｅｇｅ ｏｆ Ｃｏｍｐｕｔｅｒ Ｓｃｉｅｎｃｅ ａｎｄ Ｅｎｇｉｎｅｅｒｉｎｇ， Ｎｏｒｔｈｗｅｓｔ Ｎｏｒｍａｌ Ｕｎｉｖｅｒｓｉｔｙ， Ｌａｎｚｈｏｕ ７３００７０，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｔｏ ｉｎｃｒｅａｓｅ ｔｈｅ ｃｏｎｖｅｒｇｅｎｃｅ ｓｐｅｅｄ ａｎｄ ｓｏｌｕｔｉｏｎ ａｃｃｕｒａｃｙ ｏｆ ｔｈｅ ｔｒａｄｉｔｉｏｎａｌ ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ， ｉｎ ｔｈｉｓ ｐａｐｅｒ， ｗｅ ｐｒｏｐｏｓｅ ａ ｆｉｒｅｆｌｙ ｏｐｔｉｍｉｚａｔｉｏｎ ａｌｇｏｒｉｔｈｍ ｔｈａｔ ｕｔｉｌｉｚｅｓ ｅｌｉｔｅ ｏｐｐｏｓｉｔｉｏｎ⁃ｂａｓｅｄ ｌｅａｒｎｉｎｇ． Ｕｓｉｎｇ ａｎ ｏｐｐｏｓｉｔｉｏｎ⁃ ｂａｓｅｄ ｌｅａｒｎｉｎｇ ｓｔｒａｔｅｇｙ， ｗｅ ｃｏｎｓｔｒｕｃｔｅｄ ａｎ ｅｌｉｔｅ ｇｒｏｕｐ ａｎｄ， ｉｎ ｔｈｅ ｉｎｔｅｒｖａｌ ｏｆ ｔｈｅ ｅｌｉｔｅ ｇｒｏｕｐ， ｗｅ ｓｏｌｖｅｄ ｔｈｅ ｏｐｐｏｓｉｔｅ ｓｏｌｕｔｉｏｎｓ ｏｆ ｔｈｅ ｏｒｄｉｎａｒｙ ｇｒｏｕｐｓ． Ｔｈｉｓ ｓｔｒａｔｅｇｙ ｃｏｕｌｄ ｉｎｃｒｅａｓｅ ｇｒｏｕｐ ｄｉｖｅｒｓｉｔｙ ａｎｄ ｉｍｐｒｏｖｅ ｔｈｅ ｃｏｎｖｅｒｇｅｎｃｅ ｓｐｅｅｄ ｏｆ ｔｈｅ ａｌｇｏｒｉｔｈｍ． Ｔｏ ｐｒｅｖｅｎｔ ｔｈｅ ｏｐｔｉｍａｌ ｉｎｄｉｖｉｄｕａｌ ｆｒｏｍ ｆａｌｌｉｎｇ ｉｎｔｏ ｔｈｅ ｌｏｃａｌ ｏｐｔｉｍｕｍ， ｗｈｉｃｈ ｃｏｕｌｄ ｃａｕｓｅ ｓｔａｇｎａｔｉｏｎ ｏｆ ｔｈｅ ｗｈｏｌｅ ｇｒｏｕｐ ｄｕｒｉｎｇ ｔｈｅ ｓｅａｒｃｈ ｐｒｏｃｅｓｓ， ｗｅ ｉｎｔｒｏｄｕｃｅ ａ ｄｉｆｆｅｒｅｎｔｉａｌ ｅｖｏｌｕｔｉｏｎａｒｙ ｍｕｔａｔｉｏｎ ｓｔｒａｔｅｇｙ． Ｆｉｎａｌｌｙ， ｗｅ ｐｒｏｐｏｓｅ ａｎ ａｄａｐｔｉｖｅ ｓｔｅｐ ｓｉｚｅ ｗｉｔｈ ａ ｌｉｎｅａｒ ｄｅｃｒｅａｓｅ ｔｏ ｂａｌａｎｃｅ ｔｈｅ ｄｅｖｅｌｏｐｍｅｎｔ ａｂｉｌｉｔｙ ｏｆ ｔｈｅ ａｌｇｏｒｉｔｈｍ． Ｅｘｐｅｒｉｍｅｎｔａｌ ｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ ｔｈｅ ｐｒｏｐｏｓｅｄ ａｌｇｏｒｉｔｈｍ ｃａｎ ｉｎｃｒｅａｓｅ ｃｏｎｖｅｒｇｅｎｃｅ ｓｐｅｅｄ ａｎｄ ａｃｃｕｒａｃｙ． Ｋｅｙｗｏｒｄｓ：ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ； ｅｌｉｔｅ ｏｐｐｏｓｉｔｉｏｎ⁃ｂａｓｅｄ ｌｅａｒｎｉｎｇ； ｏｐｔｉｍｉｚｅｄ ａｌｇｏｒｉｔｈｍ； ｅｌｉｔｅ ｇｒｏｕｐ； ｏｐｐｏｓｉｔｅ ｓｏｌｕｔｉｏｎｓ； ｏｐｐｏｓｉｔｉｏｎ⁃ｂａｓｅｄ ｌｅａｒｎｉｎｇ ｓｔｒａｔｅｇｙ； ｄｉｆｆｅｒｅｎｔｉａｌ ｅｖｏｌｕｔｉｏｎａｒｙ ｍｕｔａｔｉｏｎ； ａｄａｐｔｉｖｅ ｓｔｅｐ ｓｉｚｅ 收稿日期：２０１７－０６－０７． 网络出版日期：２０１７－０８－３１． 基金项目：甘肃省科技计划资助项目（１５０６ＲＪＺＡ１３０）；甘肃省高等学校 科研项目（２０１４Ｂ⁃０１８）． 通信作者：文雅宏．Ｅ⁃ｍａｉｌ：ｗｗｙａｈｏｎｇ＠ １２６．ｃｏｍ． 萤火虫算法（ ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ，ＦＡ） 是受自然界 中萤火虫发光特性的启发， 由剑桥学者 Ｙａｎｇ 于 ２００８ 年提出的一种群体智能随机优化算法［１－５］ 。 在 多个科学与工程领域中，萤火虫算法已得到成功的 应用［６－９］ ，虽然 ＦＡ 表现出了良好的性能，但在一些 问题的优化上，萤火虫算法依然存在收敛速度慢、 解的精度不高、容易陷入局部最优等不足。 近年 来，很多学者已经进行了多角度的改进。 文献［１０］ 为了解决萤火虫算法过早地收敛和陷入局部最优 的不足，利用广义反向学习策略来优化萤火虫算 法。 文献［１１］采用正交学习策略改进 ＦＡ 算法，利 用精英萤火虫来构造指导向量，通过指导向量引导 群体向全局最优区域移动。 文献［１２］ 提出基于蛙
第5期 魏伟一，等：一种精英反向学习的萤火虫优化算法 ·711. 跳的萤火虫算法，在原始的FA算法中引入蛙跳算 有两个重要的影响因子，即亮度I和吸引度B。亮 法中的分群。同时，为了加强算法的局部开发能 度高说明其所处位置好，并吸引亮度低的个体向其 力，引入了模拟退火的思想。该算法对于高维多模 靠近。吸引度高则萤火虫移动的距离大。从FA开 态函数的优化问题，表现得还不够理想。文献[13] 始，萤火虫的个体随机地分布在指定的局域内，个 提出了一种基于多种群学习机制的萤火虫优化算 体的亮度由目标函数决定。 法，把萤火虫分为不同的子群，同时，子群建立学习 设1。表示萤火虫个体的固有亮度，Y为介质的 机制，实现不同子群间的信息交流，完成局部和全 光亮度吸收系数，「：为任意两个个体i和j的相对距 局的寻优。文献[14]引入模式搜索思想，把FA算 离（一般使用欧氏距离），B。为萤火虫个体固有吸引 法与模式搜索相结合，FA算法具有较强的全局搜索 度，随距离r变化的个体光强度1表示为 能力，模式搜索具有较好的局部搜索能力，利用两 I=loe-v2 (1) 者的优势来提高FA算法的性能。文献[15]针对高 则萤火虫i与萤火虫j之间的相互吸引力计算公 维问题，提出了多维反向学习的萤火虫算法，用反 式为 向学习策略初始化萤火虫种群。同时，用基于多维 B=B。×e-"g (2) 的方法更新不同维度上萤火虫的位置。算法在收 设x,(t)和x,(t)分别表示萤火虫i和j在t时刻 敛速度和精度上比原始萤火虫算法更优。 的位置，则两者之间的距离计算公式为 以上文献虽然对FA算法做了很好的改进，但 是在收敛速度和精度上还不够理想，为了更好地提 r=x(t)x(t)2 m)2 高FA算法的收敛速度和收敛精度，本文基于文献 (3) [16]利用精英反向学习策略来改进差分演化算法 萤火虫i向萤火虫移动，其位置更新方程为 的思想，提出了一种精英反向学习的萤火虫算法， 在文献[16]中，通过设置一个参数来选取精英个 多nu+)=x0+a0-0)+md-》 体，而本文根据原解和反向解适应度值的大小选取 (4) 精英个体，这样能更充分地利用精英群体的良好信 式中：x4(t+1)表示萤火虫i在t+1时刻的位置： 息，提高算法的收敛速度。同时，本文采用了差分 a∈[0,1]，表示步长因子。 演化策略(differential evolutionary mutation)来增强 2精英反向学习的萤火虫算法 算法的局部搜索能力。最后，为了增强和平衡算法 的开发能力，本文提出了一种线性递减的自适应步 2.1 反向学习策略 长。在5个标准测试函数上进行实验，并和多个改 反向学习策略是近年来计算智能领域出现的 进的FA算法进行实验对比。结果表明，本文算法 新概念-]，其主要思想是对一个问题的可行解， 在收敛速度和收敛精度上更好。 求其反向解，并对原解和反向解进行评估，从中选 出较优的解作为下一代个体。其中反向点和反向 1萤火虫算法(FA) 解的定义如下。 FA是受自然界中萤火虫个体通过发光来吸引 定义1反向点(opposite point,OP)[8)。设x= 同伴求偶或觅食行为的启发而提出的一种元启发 (x1,x2,…,xD)是D维空间中的一个点，且x1,x2,… 式算法[1-】，萤火虫之间相互吸引以及位置迭代更 xo∈R,x∈[a,b],则x对应的反向点x·=(x, 新的过程是搜索和优化的过程。寻找最亮萤火虫 x2,…,x0)定义为 的问题是求解最优值的问题，不断用最好的位置替 xi=ai+bi-x (5) 换较差的位置来完成整个搜索过程。在一定的搜 定义2动态反向学习策略(dynamic opposition 索区域内所有发光弱的萤火虫向发光强的萤火虫 based leaming,DOBL)。]i设x,(t)=((t),xa(t),…, 移动，从而实现位置寻优[]。每个萤火虫被看作一 xD(t))为问题第t代的一个可行解，x(t)是其j维 个个体，个体主要有“位置、亮度、吸引度”等属性， 上的分量，x(t)是x(t)对应的反向解，则x(t)是
跳的萤火虫算法，在原始的 ＦＡ 算法中引入蛙跳算 法中的分群。 同时，为了加强算法的局部开发能 力，引入了模拟退火的思想。 该算法对于高维多模 态函数的优化问题，表现得还不够理想。 文献［１３］ 提出了一种基于多种群学习机制的萤火虫优化算 法，把萤火虫分为不同的子群，同时，子群建立学习 机制，实现不同子群间的信息交流，完成局部和全 局的寻优。 文献［１４］引入模式搜索思想，把 ＦＡ 算 法与模式搜索相结合，ＦＡ 算法具有较强的全局搜索 能力，模式搜索具有较好的局部搜索能力，利用两 者的优势来提高 ＦＡ 算法的性能。 文献［１５］针对高 维问题，提出了多维反向学习的萤火虫算法，用反 向学习策略初始化萤火虫种群。 同时，用基于多维 的方法更新不同维度上萤火虫的位置。 算法在收 敛速度和精度上比原始萤火虫算法更优。 以上文献虽然对 ＦＡ 算法做了很好的改进，但 是在收敛速度和精度上还不够理想，为了更好地提 高 ＦＡ 算法的收敛速度和收敛精度，本文基于文献 ［１６］利用精英反向学习策略来改进差分演化算法 的思想，提出了一种精英反向学习的萤火虫算法， 在文献［１６］ 中，通过设置一个参数来选取精英个 体，而本文根据原解和反向解适应度值的大小选取 精英个体，这样能更充分地利用精英群体的良好信 息，提高算法的收敛速度。 同时，本文采用了差分 演化策略（ ｄｉｆｆｅｒｅｎｔｉａｌ ｅｖｏｌｕｔｉｏｎａｒｙ ｍｕｔａｔｉｏｎ） 来增强 算法的局部搜索能力。 最后，为了增强和平衡算法 的开发能力，本文提出了一种线性递减的自适应步 长。 在 ５ 个标准测试函数上进行实验，并和多个改 进的 ＦＡ 算法进行实验对比。 结果表明，本文算法 在收敛速度和收敛精度上更好。 １ 萤火虫算法（ＦＡ） ＦＡ 是受自然界中萤火虫个体通过发光来吸引 同伴求偶或觅食行为的启发而提出的一种元启发 式算法［１－２］ ，萤火虫之间相互吸引以及位置迭代更 新的过程是搜索和优化的过程。 寻找最亮萤火虫 的问题是求解最优值的问题，不断用最好的位置替 换较差的位置来完成整个搜索过程。 在一定的搜 索区域内所有发光弱的萤火虫向发光强的萤火虫 移动，从而实现位置寻优［１７］ 。 每个萤火虫被看作一 个个体，个体主要有“位置、亮度、吸引度” 等属性， 有两个重要的影响因子，即亮度 Ｉ 和吸引度 β。 亮 度高说明其所处位置好，并吸引亮度低的个体向其 靠近。 吸引度高则萤火虫移动的距离大。 从 ＦＡ 开 始，萤火虫的个体随机地分布在指定的局域内，个 体的亮度由目标函数决定。 设 Ｉ０ 表示萤火虫个体的固有亮度，γ 为介质的 光亮度吸收系数，ｒｉｊ为任意两个个体 ｉ 和 ｊ 的相对距 离（一般使用欧氏距离），β０ 为萤火虫个体固有吸引 度，随距离 ｒ 变化的个体光强度 Ｉ 表示为 Ｉ ＝ Ｉ０ ｅ －γｒ ２ （１） 则萤火虫 ｉ 与萤火虫 ｊ 之间的相互吸引力计算公 式为 β ＝ β０ × ｅ －γｒ ｉｊ （２） 设 ｘｉ（ｔ）和 ｘｊ（ｔ）分别表示萤火虫 ｉ 和 ｊ 在 ｔ 时刻 的位置，则两者之间的距离计算公式为 ｒｉｊ ＝ ‖ｘｉ（ｔ） － ｘｊ（ｔ）‖２ ＝ ∑ ｄ ｍ ＝ １ （ｘｉｍ － ｘｊｍ ） ２ （３） 萤火虫 ｉ 向萤火虫 ｊ 移动，其位置更新方程为 ｘｉ＋１ (ｔ ＋ １) ＝ ｘｉ (ｔ) ＋ β ｘｊ (ｔ) － ｘｉ ( (ｔ) ) ＋ α ｒａｎｄ － １ ２ æ è ç ö ø ÷ （４） 式中：ｘｉ＋１（ ｔ＋１） 表示萤火虫 ｉ 在 ｔ ＋ １ 时刻的位置； α∈［０，１］，表示步长因子。 ２ 精英反向学习的萤火虫算法 ２．１ 反向学习策略 反向学习策略是近年来计算智能领域出现的 新概念［１８－１９］ ，其主要思想是对一个问题的可行解， 求其反向解，并对原解和反向解进行评估，从中选 出较优的解作为下一代个体。 其中反向点和反向 解的定义如下。 定义 １ 反向点（ｏｐｐｏｓｉｔｅ ｐｏｉｎｔ，ＯＰ） ［１８］ 。 设ｘ＝ （ｘ１，ｘ２，…，ｘＤ ）是 Ｄ 维空间中的一个点，且 ｘ１，ｘ２，…， ｘＤ∈Ｒ，ｘｊ∈［ ａｊ，ｂｊ ］，则 ｘ 对应的反向点 ｘ ∗ ＝ （ ｘ ∗ １ ， ｘ ∗ ２ ，…，ｘ ∗ Ｄ ）定义为 ｘ ∗ ｊ ＝ ａｊ ＋ ｂｊ － ｘｊ （５） 定义 ２ 动态反向学习策略（ｄｙｎａｍｉｃ ｏｐｐｏｓｉｔｉｏｎ ｂａｓｅｄ ｌｅａｒｎｉｎｇ，ＤＯＢＬ）。 ［１８］设 ｘｉ（ｔ）＝ （ｘｉ１（ｔ），ｘｉ２（ｔ），…， ｘｉＤ（ｔ））为问题第 ｔ 代的一个可行解，ｘｉｊ（ ｔ）是其 ｊ 维 上的分量，ｘ ∗ ｉｊ （ｔ）是 ｘｉｊ（ｔ）对应的反向解，则 ｘ ∗ ｉ （ｔ）是 第 ５ 期 魏伟一，等：一种精英反向学习的萤火虫优化算法 ·７１１·
·712· 智能系统学报 第12卷 x:(t)对应的反向解，其中 N(t))。[a(t),b(t)]为精英群体所构造的区间， xi(t)=k(a;(t)+b(t))-x(t) (6) 当反向解越过边界[a(t),b,(t)]时，可以用下列方 式中：a(t)=min(x,(t)),b(t)=max(x(t)为当前 式进行重置： 搜索区域的最小值和最大值，其随着迭代的改变， xi=a,(t),x写b(t) D是解空间的维数：k是介于0~1的随机数。 2.3差分演化变异策略 2.2精英反向学习 在FA中，群体中的最优个体x✉引导群体向最 反向解的引入，可以扩大算法的搜索区域，但 优方向移动，如果x陷人局部最优，群体的移动终 对那些原解适应度值大于反向解适应度值的个体， 止，即收敛到局部最优，则群体无法到达全局最优。 对其进行反向区域的搜索，浪费时间，则应加强其 因此，本文为了求得全局最优解，引入差分变异策 领域搜索。而对原解适应度值小于反向解适应度 略，对x进行变异操作，使其陷入局部最优的概率 值的个体，对其进行反向区域的搜素价值要高于其 减小。 领域的开发价值。因此，本文将原解适应度值小于 本文要对最优个体进行变异操作，使其跳出局 反向解适应度值的个体作为研究对象，求其反向 部最优的概率增大，因此选择“DE/best/1”作为变 解，既可以扩大搜素区域，也能有效避免盲目搜索 异操作，公式为 带来的时间浪费。 Cg=xe,+F·(xn1J-xa2） (9) 同时，本文为了提高算法的收敛速度，首先在当前 式中：x,为最优个体的第j维；F是缩放系数；n1、 解所构造的空间中，求所有当前解的反向解：然后，通 n2是[1，n]上两个互不相同的随机整数，代表不同 过比较适应度值，选出那些原解适应度值大于反向解 个体的下标；j是维度：c是变异后的值。 适应度值的个体组成精英群体：最后，在精英群体构造 将变异后的个体和父代个体进行如下交叉 的新的搜索空间上，再求原解适应度值小于反向解适 操作： 应度值的个体的反向解。如果算法能收敛到全局最优 Ci' rand[0,1]≤CR或者j=rand(1,D) 解，则精英群体所形成的搜索区间必将收敛到最优解 .j' 其他 所在的区域，这样充分利用了精英群体的有效信息， (10) 在精英群体所构成的动态定义区间上生成反向解，引 式中：x为新生成个体的第j维；rand[0,1]是 导搜索向最优解靠近。 [0,1]的随机数：CR是交叉概率；rand(1,D)是[1， 定义3精英(elite)。设x,(t)=(xa,x2,…, D]的一个随机整数；参数F、CR分别设为1和0.1。 xD)是第1次迭代的一个解，其反向解为x(t), 2.4自适应步长 f(x)为目标函数。当f(x,(t)≥f八x:(t))时，称 在原始萤火虫算法FA中，步长因子α在每次 x,(t)为第t次迭代的精英个体，记为N,(t);当 迭代时保持不变。但是当α取较大的值时，增强了 f八x:(t))<f八x(t))时，称x,(t)为第t次迭代的普 算法的全局搜索能力，降低了算法的收敛速度和搜 索的精度：当α取较小的值时，有利于算法的局部 通个体，记为Q(t)。若精英群体的规模为p(1< p≤n,n为解的总个数)时，则p个精英个体可表示 搜索，提高了搜索精度和算法的收敛速度。在算法 为{N(t),N2(t),…,Nn(t)}C{x(t),x2(t),…, 迭代前期，较大的α有利于算法的全局搜索：在后 期，较小的α显得更有利。因此本文对α采用动态 x(t)}。 递减的方式，计算公式为 定义4精英反向解(elite opposite solution)i。 t 设x为普通个体x:在j维上的值，则其反向解可定 T- 10 义为 a,1=a·（T),t=1,2,…,T(11) xi(t)=k(a;(t)+b(t))-x(t) (7)》 式中t为迭代次数。 式中：k是介于0~l的随机数；a(t)=min(V,(t), 2.5E0FA算法描述 N2(t),…,N(t)）;b(t)=max(N(t),N2(t),…, EOFA算法流程如下
ｘｉ（ｔ）对应的反向解，其中 ｘ ∗ ｉｊ （ｔ） ＝ ｋ（ａｊ（ｔ） ＋ ｂｊ（ｔ）） － ｘｉｊ（ｔ） （６） 式中： ａ（ｔ） ＝ ｍｉｎ（ｘｉｊ（ｔ）），ｂ（ｔ） ＝ ｍａｘ（ｘｉｊ（ｔ）） 为当前 搜索区域的最小值和最大值，其随着迭代的改变， 而发生变化； ｉ ∈ ［１，ｎ］，ｊ ∈ ［１，Ｄ］；ｎ 是种群大小； Ｄ 是解空间的维数；ｋ 是介于 ０～１ 的随机数。 ２．２ 精英反向学习 反向解的引入，可以扩大算法的搜索区域，但 对那些原解适应度值大于反向解适应度值的个体， 对其进行反向区域的搜索，浪费时间，则应加强其 领域搜索。 而对原解适应度值小于反向解适应度 值的个体，对其进行反向区域的搜素价值要高于其 领域的开发价值。 因此，本文将原解适应度值小于 反向解适应度值的个体作为研究对象，求其反向 解，既可以扩大搜素区域，也能有效避免盲目搜索 带来的时间浪费。 同时，本文为了提高算法的收敛速度，首先在当前 解所构造的空间中，求所有当前解的反向解；然后，通 过比较适应度值，选出那些原解适应度值大于反向解 适应度值的个体组成精英群体；最后，在精英群体构造 的新的搜索空间上，再求原解适应度值小于反向解适 应度值的个体的反向解。 如果算法能收敛到全局最优 解，则精英群体所形成的搜索区间必将收敛到最优解 所在的区域［１６］ ，这样充分利用了精英群体的有效信息， 在精英群体所构成的动态定义区间上生成反向解，引 导搜索向最优解靠近。 定义 ３ 精英（ ｅｌｉｔｅ）。 设 ｘｉ（ｔ） ＝ （ｘｉ１ ，ｘｉ２ ，…， ｘｉＤ） 是第 ｔ 次迭代的一个解，其反向解为 ｘ ∗ ｉ （ｔ）， ｆ（ｘ） 为目标函数。 当 ｆ（ｘｉ（ｔ）） ≥ ｆ（ｘ ∗ ｉ （ｔ）） 时，称 ｘｉ（ｔ） 为第 ｔ 次迭代的精英个体， 记为 Ｎｉ（ｔ） ； 当 ｆ（ｘｉ（ｔ）） ＜ ｆ（ｘ ∗ ｉ （ｔ）） 时，称 ｘｔ（ｔ） 为第 ｔ 次迭代的普 通个体，记为 Ｑｉ（ｔ）。 若精英群体的规模为 ｐ（１ ＜ ｐ ≤ｎ，ｎ 为解的总个数） 时，则 ｐ 个精英个体可表示 为 ｛Ｎ１（ｔ），Ｎ２（ｔ），…，Ｎｐ（ｔ）｝ ⊆ ｛ｘ１（ｔ），ｘ２（ｔ），…， ｘｎ（ｔ）｝。 定义 ４ 精英反向解（ｅｌｉｔｅ ｏｐｐｏｓｉｔｅ ｓｏｌｕｔｉｏｎ） ［１８］ 。 设 ｘｉｊ 为普通个体 ｘｉ 在 ｊ 维上的值，则其反向解可定 义为 ｘ ∗ ｉｊ （ｔ） ＝ ｋ（ａｊ（ｔ） ＋ ｂｊ（ｔ）） － ｘｉｊ（ｔ） （７） 式中：ｋ 是介于 ０ ～ １ 的随机数； ａｊ（ｔ） ＝ ｍｉｎ（Ｎ１ｊ（ｔ）， Ｎ２ｊ（ｔ），…，Ｎｐｊ（ｔ））；ｂｊ（ｔ） ＝ ｍａｘ（Ｎ１ｊ（ｔ），Ｎ２ｊ（ｔ），…， Ｎｐｊ（ｔ））。 ［ａｊ（ｔ），ｂｊ（ｔ）］ 为精英群体所构造的区间， 当反向解越过边界 ［ａｊ（ｔ），ｂｊ（ｔ）］ 时，可以用下列方 式进行重置： ｘ ∗ ｉｊ ＝ ａｊ（ｔ）， ｘ ∗ ｉｊ ＜ ａｊ（ｔ） ｘ ∗ ｉｊ ＝ ｂｊ（ｔ）， ｘ ∗ ｉｊ ＞ ｂｊ（ｔ） { （８） ２．３ 差分演化变异策略 在 ＦＡ 中，群体中的最优个体 ｘｂｅｓｔ引导群体向最 优方向移动，如果 ｘｂｅｓｔ 陷入局部最优，群体的移动终 止，即收敛到局部最优，则群体无法到达全局最优。 因此，本文为了求得全局最优解，引入差分变异策 略，对 ｘｂｅｓｔ 进行变异操作，使其陷入局部最优的概率 减小。 本文要对最优个体进行变异操作，使其跳出局 部最优的概率增大，因此选择“ＤＥ ／ ｂｅｓｔ ／ １” 作为变 异操作，公式为 ｃｉｊ ＝ ｘｂｅｓｔ， ｊ ＋ Ｆ· ｘｎ１ ，ｊ － ｘｎ２ ，ｊ ( ) （９） 式中： ｘｂｅｓｔ，ｊ 为最优个体的第 ｊ 维； Ｆ 是缩放系数； ｎ１ 、 ｎ２ 是 ［１，ｎ］ 上两个互不相同的随机整数，代表不同 个体的下标； ｊ 是维度；ｃｉｊ是变异后的值。 将变异后的个体和父代个体进行如下交叉 操作： ｘ － ｂｅｓｔ， ｊ ＝ ｃｉｊ， ｒａｎｄ［０，１］ ≤ ＣＲ 或者 ｊ ＝ ｒａｎｄ（１，Ｄ） ｘｂｅｓｔ， ｊ， 其他 { （１０） 式中：ｘ － ｂｅｓｔ，ｊ 为新生成个体的第 ｊ 维； ｒａｎｄ［０，１］ 是 ［０，１］的随机数；ＣＲ 是交叉概率；ｒａｎｄ（１，Ｄ）是［１， Ｄ］的一个随机整数；参数 Ｆ、ＣＲ 分别设为 １ 和 ０．１。 ２．４ 自适应步长 在原始萤火虫算法 ＦＡ 中，步长因子 α 在每次 迭代时保持不变。 但是当 α 取较大的值时，增强了 算法的全局搜索能力，降低了算法的收敛速度和搜 索的精度；当 α 取较小的值时，有利于算法的局部 搜索，提高了搜索精度和算法的收敛速度。 在算法 迭代前期，较大的 α 有利于算法的全局搜索；在后 期，较小的 α 显得更有利。 因此本文对 α 采用动态 递减的方式，计算公式为 αｔ＋１ ＝ αｔ·（ Ｔ － ｔ １０ Ｔ ）， ｔ ＝ １，２，…，Ｔ （１１） 式中 ｔ 为迭代次数。 ２．５ ＥＯＦＡ 算法描述 ＥＯＦＡ 算法流程如下。 ·７１２· 智 能 系 统 学 报 第 １２ 卷
第5期 魏伟一，等：一种精英反向学习的萤火虫优化算法 ·713· 输入目标函数和搜索空间： l)Sphere函数 输出全局最优解和最优位置。 1)初始化参数m,n,T,a,B,y,在[m,n]上生成 )=2.5e（-512,5.12 i=1 初始种群x,(i=1,2,…,n)。 Sphere函数为多维单峰值函数，在点x=(0,0,, 2)执行EOFA算法搜索。 0)处取得极小值0。 3)把x(t)(i=1,2,…,n)带入目标函数，计算 2)Rosenbrock函数 函数值，把目标函数的值作为每个个体的亮度 值1.(t)。 f=[100(,-》+(k-1)1. 4)在[m,n]生成x,(t)的反向解x,(t),并计算 x:∈(-2.048,2.048) 亮度I(t)。对I(t)和(t)进行比较，fI(t)>I Rosenbrock函数为多维病态二次函数，在点x= (t),则x,(t)为精英个体N(t),记精英群体大小为 (1,1,…,1)处取得全局极小值0。 p(p>1,i=1,2,…,p)。设精英个体的区间范围为 3)Ackley函数 [a(t),b(t)](若精英群体的规模小于2，则在x (t)(i=1,2,…,n)构成的区间上求x(i=1,2,…, -e点m+2.7128, f3(x)=-20e√名 n)的反向解)。elsex,(t)是普通个体，普通群体大小 x∈（-32.7,32.7) 为n-po Ackley函数为多维多峰值函数，在点x=(0,0,…, 5)用式(7)在精英个体构成的区间[a,(t),b,(t)] 0)处取得全局极小值0。 上计算普通群体的反向解x:'(t)(i=1,…,n-p)。 4)Griewank函数 6)精英群体和普通群体的反向解群体构成当 f(x)= os()+1， 前新种群，计算新种群的亮度，并进行排序，选出最 优的个体x(t)。 x:∈(-600,600) 7)用式(3)计算每个个体i和最优个体x(t) Griewank函数为多维多峰值函数，在点x=(0,0, 之间的距离R(t)。 …,0)处取得全局极小值0。 8)用式(2)计算吸引力B(t)。个体i向最优个 5)Rastrigin函数 体(t)移动，用公式(6)更新位置x“(t)。 9)用式(11)计算a(t),并用式(9)、(10)对最 f(x)=10d+ A-10e(2小， 优个体进行位置扰动。 x:∈(-5.12,5.12) 10)算法搜索结束，输出全局最优解和最优 Rastrigin函数为多维多峰值函数，在点x=(0,0,…, 位置。 0)处取得全局极小值0。 若种群的规模为n,空间维度为D,则种群初始 3.2EOFA算法的测试结果 化的时间复杂度为O(nD);从迭代开始到结束的整 实验环境为：Inter Core(TM)i5-2450MCPU@ 个过程中，迭代的次数为t,其中3)是计算种群的亮 2.50GHz,内存4GB,Window7操作系统，MATLAB 度，复杂度为O(nDt):4)~9)是建立新的种群，并进 7.8.0版本。分别选取标准的FA算法[)，LFA算 行位置的更新，复杂度为O((6n-p)·D),p(p≤n) 法[20]，MFA算法[2)与本文提出的E0FA算法在5 是精英群体的规模。因此本文算法的时间复杂度 为O(nDr)。 种标准的测试函数上进行实验比较，种群规模n取 40,初始a取值为0.98。维度D取10和30，y=1, 3实验仿真及分析 MFA算法中方向向量的个数m取30，其他参数分 3.1测试函数 别取T=1000,B=1。分别记录4种算法迭代1000 在仿真实验中，本文采用下列5个常用的标准 次并在测试函数上独立运行40次的最优值、最差值 测试函数对算法进行测试。 和平均值，结果如表1所示
输入 目标函数和搜索空间； 输出 全局最优解和最优位置 。 １）初始化参数 ｍ，ｎ，Ｔ，α，β，γ，在［ｍ，ｎ］上生成 初始种群 ｘｉ（ｉ ＝ １，２，…，ｎ）。 ２）执行 ＥＯＦＡ 算法搜索。 ３）把 ｘｉ（ｔ） （ｉ ＝ １，２，…，ｎ）带入目标函数，计算 函数值， 把 目 标 函 数 的 值 作 为 每 个 个 体 的 亮 度 值Ｉｉ（ｔ）。 ４）在［ｍ，ｎ］生成 ｘｉ（ｔ）的反向解 ｘ ∗ ｉ （ｔ），并计算 亮度 Ｉ ∗ ｉ （ｔ）。 对 Ｉｉ（ｔ）和 Ｉ ∗ ｉ （ｔ）进行比较，ｉｆ Ｉｉ（ｔ） ＞Ｉ ∗ ｉ （ｔ），则 ｘｉ（ｔ）为精英个体 Ｎｉ（ｔ），记精英群体大小为 ｐ（ｐ＞１，ｉ ＝ １，２，…，ｐ）。 设精英个体的区间范围为 ［ａｊ（ｔ），ｂｊ（ ｔ）］ （若精英群体的规模小于 ２，则在 ｘｉ （ｔ）（ｉ ＝ １，２，…，ｎ）构成的区间上求 ｘｉ（ ｉ ＝ １，２，…， ｎ）的反向解）。 ｅｌｓｅｘｉ（ｔ）是普通个体，普通群体大小 为ｎ－ｐ。 ５）用式（７）在精英个体构成的区间［ａｊ（ｔ），ｂｊ（ｔ）］ 上计算普通群体的反向解 ｘｉ ′（ｔ）（ｉ ＝１，…，ｎ－ｐ）。 ６）精英群体和普通群体的反向解群体构成当 前新种群，计算新种群的亮度，并进行排序，选出最 优的个体 ｘｂｅｓｔ（ｔ） 。 ７）用式（３）计算每个个体 ｉ 和最优个体 ｘｂｅｓｔ（ｔ） 之间的距离 Ｒｉｊ（ｔ） 。 ８）用式（２）计算吸引力 βｉｊ（ｔ）。 个体 ｉ 向最优个 体 ｘｂｅｓｔ（ｔ） 移动，用公式（６）更新位置 ｘ ｎｅｗ ｉ （ｔ） 。 ９）用式（１１）计算 α（ ｔ），并用式（９）、（１０）对最 优个体进行位置扰动。 １０）算法搜索结束，输出全局最优解和最优 位置。 若种群的规模为 ｎ，空间维度为 Ｄ，则种群初始 化的时间复杂度为 Ｏ（ｎＤ）；从迭代开始到结束的整 个过程中，迭代的次数为 ｔ，其中 ３）是计算种群的亮 度，复杂度为 Ｏ（ｎＤｔ）；４） ～９）是建立新的种群，并进 行位置的更新，复杂度为 Ｏ（（６ｎ－ｐ）·Ｄｔ），ｐ（ｐ≤ｎ） 是精英群体的规模。 因此本文算法的时间复杂度 为 Ｏ（ｎＤｔ）。 ３ 实验仿真及分析 ３．１ 测试函数 在仿真实验中，本文采用下列 ５ 个常用的标准 测试函数对算法进行测试。 １） Ｓｐｈｅｒｅ 函数 ｆ １（ｘ） ＝ ∑ ｄ ｉ ＝ １ ｘ ２ ｉ ， ｘｉ ∈ （ － ５．１２，５．１２） Ｓｐｈｅｒｅ 函数为多维单峰值函数，在点 ｘ ＝ （０，０，…， ０）处取得极小值 ０。 ２） Ｒｏｓｅｎｂｒｏｃｋ 函数 ｆ ２（ｘ） ＝ ∑ ｄ－１ ｉ ＝ １ ［１００ （ｘｉ＋１ － ｘ ２ ｉ ） ２ ＋ （ｘｉ － １） ２ ］， ｘｉ ∈ （ － ２．０４８，２．０４８） Ｒｏｓｅｎｂｒｏｃｋ 函数为多维病态二 次 函 数， 在 点 ｘ ＝ （１，１，…，１）处取得全局极小值 ０。 ３） Ａｃｋｌｅｙ 函数 ｆ ３（ｘ） ＝ － ２０ｅ －０．２ １ ｄ ∑ ｄ ｉ ＝ １ ｘ ２ ｉ － ｅ １ ｄ ∑ ｄ ｉ ＝ １ ｃｏｓ（２πｘ ｉ ） ＋ ２２．７１２ ８， ｘｉ ∈ （ － ３２．７，３２．７） Ａｃｋｌｅｙ 函数为多维多峰值函数，在点 ｘ ＝ （０，０，…， ０）处取得全局极小值 ０。 ４） Ｇｒｉｅｗａｎｋ 函数 ｆ ４（ｘ） ＝ ∑ ｄ ｉ ＝ １ ｘ ２ ｉ ４ ０００ － ∏ ｄ ｉ ＝ １ ｃｏｓ（ ｘ ｉ ｉ ） ＋ １， ｘｉ ∈ （ － ６００，６００） Ｇｒｉｅｗａｎｋ 函数为多维多峰值函数，在点 ｘ ＝ （ ０，０， …，０）处取得全局极小值 ０。 ５） Ｒａｓｔｒｉｇｉｎ 函数 ｆ ５（ｘ） ＝ １０ｄ ＋ ∑ ｄ ｉ ＝ １ ［ｘ ２ ｉ － １０ｃｏｓ（２πｘｉ）］， ｘｉ ∈ （ － ５．１２，５．１２） Ｒａｓｔｒｉｇｉｎ 函数为多维多峰值函数，在点 ｘ ＝ （０，０，…， ０）处取得全局极小值 ０。 ３．２ ＥＯＦＡ 算法的测试结果 实验环境为：Ｉｎｔｅｒ Ｃｏｒｅ（ ＴＭ） ｉ５⁃２４５０Ｍ ＣＰＵ＠ ２．５０ ＧＨｚ，内存 ４ ＧＢ，Ｗｉｎｄｏｗ７ 操作系统，ＭＡＴＬＡＢ ７．８． ０ 版本。 分别选取标准的 ＦＡ 算法［１］ ，ＬＦＡ 算 法［２０］ ，ＭＦＡ 算法［２１］ 与本文提出的 ＥＯＦＡ 算法在 ５ 种标准的测试函数上进行实验比较，种群规模 ｎ 取 ４０，初始 α 取值为 ０．９８。 维度 Ｄ 取 １０ 和 ３０，γ ＝ １， ＭＦＡ 算法中方向向量的个数 ｍ 取 ３０，其他参数分 别取Ｔ ＝ １ ０００，β ＝ １。 分别记录 ４ 种算法迭代 １ ０００ 次并在测试函数上独立运行 ４０ 次的最优值、最差值 和平均值，结果如表 １ 所示。 第 ５ 期 魏伟一，等：一种精英反向学习的萤火虫优化算法 ·７１３·
.714. 智能系统学报 第12卷 表14种算法的实验结果 Table 1 Experimental results of four algorithms 函数 10维 30维 算法 最优值 最差值 平均值 最优值 最差值 平均值 FA 2.7502×103 6.2938×103 4.2120×10 3.7861×10 8.9832×10 5.0921×10 MFA 3.811×10- 4.5317×10-6 3.6714×10-6 4.2306×102 7.4644×102 5.4182×102 LFA 1.7142×10-3 7.9325×10- 4.5214×103 1.7412×10- 1.2901 7.1094×10 EOFA 3.7941×108 2.3972×10-2 2.0069x10B 1.0092×10-9 3.1434×10-9 2.0041×10-9 FA 7.7421×10- 8.1837×10- 7.9224×10- 2.9013 5.3953 3.0391 MFA 4.0928×102 1.0928×10 13420×10 1.0085 4.9032 2.6059 LFA 2.0492×10- 8.4952×10 4.1182×10 1.1196 8.9013 4.9081 EOFA 3.2920×10-4 5.3898×104 4.2539×10 1.8046×10-2 2.2333×102 2.0110×102 FA 2.395×10- 5.0919×10- 3.3183×10- 4.3942×10- 7.0392×10 5.3092×10 MFA 7.3478×10-6 8.8721×10-6 7.8134×10-6 3.7314×10-2 5.0930×10-2 4.2091×102 LFA 7.3517×10- 2.9205×102 1.4518×102 2.3971×10-1 6.7836×10- 3.0240×10~1 EOFA 2.9766×10-8 5.4730×10 3.4613×10-8 1.3410×10-6 1.7892×10-6 1.5753×10-6 FA 4.4720×103 3.0951×102 1.9771×102 3.0501×10- 5.9081×10 4.3018×10 MFA 1.0978×10 4.0928×10+ 2.4980×10 3.0958×102 5.0968×10-2 4.4103×102 LFA 3.1192×10- 8.0937×10-3 5.4283×10- 2.1976×10 5.8921×10 3.8349×10 EOFA 1.0411×109 6.1995×10-9 3.6203×10-9 1.2169×10-6 3.9001×10-6 2.7101×10-6 FA 3.0967×10- 6.3187×101 4.8301×10 5.5430 7.1025 6.5018 MFA 4.8241×10-2 2.8710×10- 1.3027×10- 3.9816 5.1183 4.1510 LFA 6.3387×10- 7.2014×10 6.5211×10 5.9082 9.2205 7.2098 EOFA 2.0689×10-0 4.0900×10-10 3.1325×1010 5.8798e-08 8.3268×10-6 4.5036×10-6 FA是原始的萤火虫算法。LFA是根据Levy分 FA 布来设置一种随机步长对传统萤火虫算法进行改 LFA 10° -MFA --EOFA 进，其主要优点是算法收敛到局部最优的概率降 10 低。MFA算法是从随机生成的方向向量中选择使 种群进化到最优的方向向量，方向向量的个数对算 10W 法的性能有很大的影响，数量越大，算法收敛性越 10 ×10 00.10.20.30.40.50.60.70.80.91.0 好。由表1可知，EOFA、LFA和MFA算法在10维 迭代次数 和30维函数上都优于FA算法。本文提出的EOFA (a)f:Sphere函数 算法在5种测试函数上的函数值都小于FA、LFA、 10㎡ MFA算法在测试函数上的值，即EOFA算法的收敛 10 性更好，在每个测试函数上EOFA算法的求解精度 EOFA 10 比其他3种算法都高。 为了更好地验证EOFA算法的有效性，本文用 10 图描述4种算法的收敛性，由于受篇幅的限制，仅给 10 ×10 00.10.20.30.40.50.60.70.80.91.0 出4个代表性的函数收敛曲线图，结果如图1、2 迭代次数 所示。 (b)f5:Rosenbrock函数
表 １ ４ 种算法的实验结果 Ｔａｂｌｅ １ Ｅｘｐｅｒｉｍｅｎｔａｌ ｒｅｓｕｌｔｓ ｏｆ ｆｏｕｒ ａｌｇｏｒｉｔｈｍｓ 函数 算法 １０ 维 ３０ 维 最优值 最差值 平均值 最优值 最差值 平均值 ｆ １ ＦＡ ２．７５０ ２×１０ －３ ６．２９３ ８×１０ －３ ４．２１２ ０×１０ －３ ３．７８６ １×１０ －１ ８．９８３ ２×１０ －１ ５．０９２ １×１０ －１ ＭＦＡ ３．８１１×１０ －６ ４．５３１ ７×１０ －６ ３．６７１ ４×１０ －６ ４．２３０ ６×１０ －２ ７．４６４ ４×１０ －２ ５．４１８ ２×１０ －２ ＬＦＡ １．７１４ ２×１０ －３ ７．９３２ ５×１０ －３ ４．５２１ ４×１０ －３ １．７４１ ２×１０ －１ １．２９０ １ ７．１０９ ４×１０ －１ ＥＯＦＡ ３．７９４ １×１０ －１３ ２．３９７ ２×１０ －１２ ２．００６ ９×１０ －１３ １．００９ ２×１０ －９ ３．１４３ ４×１０ －９ ２．００４ １×１０ －９ ｆ ２ ＦＡ ７．７４２ １×１０ －１ ８．１８３ ７×１０ －１ ７．９２２ ４×１０ －１ ２．９０１ ３ ５．３９５ ３ ３．０３９ １ ＭＦＡ ４．０９２ ８×１０ －２ １．０９２ ８×１０ －１ １ ３４２ ０×１０ －１ １．００８ ５ ４．９０３ ２ ２．６０５ ９ ＬＦＡ ２．０４９ ２×１０ －１ ８．４９５ ２×１０ －１ ４．１１８ ２×１０ －１ １．１１９ ６ ８．９０１ ３ ４．９０８ １ ＥＯＦＡ ３．２９２ ０×１０ －４ ５．３８９ ８×１０ －４ ４．２５３ ９×１０ －４ １．８０４ ６×１０ －２ ２．２３３ ３×１０ －２ ２．０１１ ０×１０ －２ ｆ ３ ＦＡ ２．３９５×１０ －３ ５．０９１ ９×１０ －３ ３．３１８ ３×１０ －３ ４．３９４ ２×１０ －１ ７．０３９ ２×１０ －１ ５．３０９ ２×１０ －１ ＭＦＡ ７．３４７ ８×１０ －６ ８．８７２ １×１０ －６ ７．８１３ ４×１０ －６ ３．７３１ ４×１０ －２ ５．０９３ ０×１０ －２ ４．２０９ １×１０ －２ ＬＦＡ ７．３５１ ７×１０ －３ ２．９２０ ５×１０ －２ １．４５１ ８×１０ －２ ２．３９７ １×１０ －１ ６．７８３ ６×１０ －１ ３．０２４ ０×１０ －１ ＥＯＦＡ ２．９７６ ６×１０ －８ ５．４７３ ０×１０ －８ ３．４６１ ３×１０ －８ １．３４１ ０×１０ －６ １．７８９ ２×１０ －６ １．５７５ ３×１０ －６ ｆ ４ ＦＡ ４．４７２ ０×１０ －３ ３．０９５ １×１０ －２ １．９７７ １×１０ －２ ３．０５０ １×１０ －１ ５．９０８ １×１０ －１ ４．３０１ ８×１０ －１ ＭＦＡ １．０９７ ８×１０ －４ ４．０９２ ８×１０ －４ ２．４９８ ０×１０ －４ ３．０９５ ８×１０ －２ ５．０９６ ８×１０ －２ ４．４１０ ３×１０ －２ ＬＦＡ ３．１１９ ２×１０ －３ ８．０９３ ７×１０ －３ ５．４２８ ３×１０ －３ ２．１９７ ６×１０ －１ ５．８９２ １×１０ －１ ３．８３４ ９×１０ －１ ＥＯＦＡ １．０４１ １×１０ －９ ６．１９９ ５×１０ －９ ３．６２０ ３×１０ －９ １．２１６ ９×１０ －６ ３．９００ １×１０ －６ ２．７１０ １×１０ －６ ｆ ５ ＦＡ ３．０９６ ７×１０ －１ ６．３１８ ７×１０ －１ ４．８３０ １×１０ －１ ５．５４３ ０ ７．１０２ ５ ６．５０１ ８ ＭＦＡ ４．８２４ １×１０ －２ ２．８７１ ０×１０ －１ １．３０２ ７×１０ －１ ３．９８１ ６ ５．１１８ ３ ４．１５１ ０ ＬＦＡ ６．３３８ ７×１０ －１ ７．２０１ ４×１０ －１ ６．５２１ １×１０ －１ ５．９０８ ２ ９．２２０ ５ ７．２０９ ８ ＥＯＦＡ ２．０６８ ９×１０ －１０ ４．０９０ ０×１０ －１０ ３．１３２ ５×１０ －１０ ５．８７９ ８ｅ⁃０８ ８．３２６ ８×１０ －６ ４．５０３ ６×１０ －６ ＦＡ 是原始的萤火虫算法。 ＬＦＡ 是根据 Ｌｅｖｙ 分 布来设置一种随机步长对传统萤火虫算法进行改 进，其主要优点是算法收敛到局部最优的概率降 低。 ＭＦＡ 算法是从随机生成的方向向量中选择使 种群进化到最优的方向向量，方向向量的个数对算 法的性能有很大的影响，数量越大，算法收敛性越 好。 由表 １ 可知，ＥＯＦＡ、ＬＦＡ 和 ＭＦＡ 算法在 １０ 维 和 ３０ 维函数上都优于 ＦＡ 算法。 本文提出的 ＥＯＦＡ 算法在 ５ 种测试函数上的函数值都小于 ＦＡ、ＬＦＡ、 ＭＦＡ 算法在测试函数上的值，即 ＥＯＦＡ 算法的收敛 性更好，在每个测试函数上 ＥＯＦＡ 算法的求解精度 比其他 ３ 种算法都高。 为了更好地验证 ＥＯＦＡ 算法的有效性，本文用 图描述 ４ 种算法的收敛性，由于受篇幅的限制，仅给 出 ４ 个代表性的函数收敛曲线图，结果如图 １、２ 所示。 （ａ） ｆ １ ：Ｓｐｈｅｒｅ 函数 （ｂ） ｆ ２ ：Ｒｏｓｅｎｂｒｏｃｋ 函数 ·７１４· 智 能 系 统 学 报 第 １２ 卷
第5期 魏伟一，等：一种精英反向学习的萤火虫优化算法 .715· 10 10 -FA 10 …LFA --…MfA 109 EOFA ---EOFA 10 10- 100 ×10 10 00.10.20.30.40.50.60.70.80.91.0 00.10.20.30.40.50.60.70.80.91.0 迭代次数 迭代次数 (c)f方：Ackley函数 (c)方：Ackley函数 10 -EA S10 LFA -MFA ---EOFA 10 EOFA 10 10 0 0.102030.40$0.6070809i0*10 迭代次数 10 ×10 00.10.20.30.40.50.60.70.80.91.0 (d)f:Griewank函数 迭代次数 图2维度为30时算法收敛曲线对比 Fig.2 Comparison of convergent graphs for a (d)f:Griewank函数 dimensionality of 30 图1维度为10时算法收敛曲线对比 由图1、2可以看出，对于每一个测试函数， Fig.1 Comparison of convergent graphs for a EOFA算法总比其他3种算法表现出更好的收敛 dimensionality of 10 性，因为它构建了动态的精英反向解区间，同时，精 英群体的规模自适应的改变，使普通个体向最优个 体移动的速度加快。对于函数f、f,当函数维度为 FA 10和30时，FA、LFA、MFA这3种算法出现了早熟， A -…MFA 而EOFA算法继续收敛，且收敛速度比其他3种算 109 ---EOFA 法都快。对于函数方、5虽然4种算法都出现了早 熟，但EOFA算法解的精确度比其他算法更好。当 10 维数从10增加到30时，4种算法的性能都有所下 降，但EOFA算法的性能优于其他3种算法。EOFA 算法具有较优性能的原因是：首先EOFA算法采用 100 ×103 00.10.20.30.40.50.60.70.80.91.0 反向学习策略，构造精英群体和普通群体，扩大了 迭代次数 搜索范围，通过生成每个个体的反向解，增加解的 (a)f:Sphere函数 多样性来提高种群多样性。同时加入了差分演化 变异策略，使其跳出局部最优的概率增大。最后， 103 为了增强算法开发能力，采用递减的自适应步长。 FA a -MEA 4结束语 EOFA 本文提出的精英反向学习萤火虫算法 (EOFA),通过精英反向学习策略生成当前解的反 向解，评估当前解和反向解，构建精英群体和普通 群体，增加了群体的多样性：在动态的精英区间上 102 二×10 0.10.2030.40.50.60.70.80.91.0 求普通群体的反向解，提高算法的收敛速度。差分 迭代次数 演化变异策略对最优个体进行变异操作，对其领域 (b)f:Rosenbrock函数 空间进行搜索，增强了EOFA的局部开采能力。同 时，采用自适应步长，提高和平衡算法的开发能力
（ｃ） ｆ ３ ：Ａｃｋｌｅｙ 函数 （ｄ） ｆ ４ ：Ｇｒｉｅｗａｎｋ 函数 图 １ 维度为 １０ 时算法收敛曲线对比 Ｆｉｇ． １ Ｃｏｍｐａｒｉｓｏｎ ｏｆ ｃｏｎｖｅｒｇｅｎｔ ｇｒａｐｈｓ ｆｏｒ ａ ｄｉｍｅｎｓｉｏｎａｌｉｔｙ ｏｆ １０ （ａ） ｆ １ ：Ｓｐｈｅｒｅ 函数 （ｂ） ｆ ２ ：Ｒｏｓｅｎｂｒｏｃｋ 函数 （ｃ） ｆ ３ ：Ａｃｋｌｅｙ 函数 （ｄ） ｆ ４ ：Ｇｒｉｅｗａｎｋ 函数 图 ２ 维度为 ３０ 时算法收敛曲线对比 Ｆｉｇ．２ Ｃｏｍｐａｒｉｓｏｎ ｏｆ ｃｏｎｖｅｒｇｅｎｔ ｇｒａｐｈｓ ｆｏｒ ａ ｄｉｍｅｎｓｉｏｎａｌｉｔｙ ｏｆ ３０ 由图 １、２ 可以看出，对于每一个测试函数， ＥＯＦＡ 算法总比其他 ３ 种算法表现出更好的收敛 性，因为它构建了动态的精英反向解区间，同时，精 英群体的规模自适应的改变，使普通个体向最优个 体移动的速度加快。 对于函数 ｆ １ 、 ｆ ４ ，当函数维度为 １０ 和 ３０ 时，ＦＡ、ＬＦＡ、ＭＦＡ 这 ３ 种算法出现了早熟， 而 ＥＯＦＡ 算法继续收敛，且收敛速度比其他 ３ 种算 法都快。 对于函数 ｆ ２ 、 ｆ ３ 虽然 ４ 种算法都出现了早 熟，但 ＥＯＦＡ 算法解的精确度比其他算法更好。 当 维数从 １０ 增加到 ３０ 时，４ 种算法的性能都有所下 降，但 ＥＯＦＡ 算法的性能优于其他 ３ 种算法。 ＥＯＦＡ 算法具有较优性能的原因是：首先 ＥＯＦＡ 算法采用 反向学习策略，构造精英群体和普通群体，扩大了 搜索范围，通过生成每个个体的反向解，增加解的 多样性来提高种群多样性。 同时加入了差分演化 变异策略，使其跳出局部最优的概率增大。 最后， 为了增强算法开发能力，采用递减的自适应步长。 ４ 结束语 本文 提 出 的 精 英 反 向 学 习 萤 火 虫 算 法 （ＥＯＦＡ），通过精英反向学习策略生成当前解的反 向解，评估当前解和反向解，构建精英群体和普通 群体，增加了群体的多样性；在动态的精英区间上 求普通群体的反向解，提高算法的收敛速度。 差分 演化变异策略对最优个体进行变异操作，对其领域 空间进行搜索，增强了 ＥＯＦＡ 的局部开采能力。 同 时，采用自适应步长，提高和平衡算法的开发能力。 第 ５ 期 魏伟一，等：一种精英反向学习的萤火虫优化算法 ·７１５·
·716 智能系统学报 第12卷 通过实验结果得出，EOFA算法在解的精度和收敛 LI Yang.Leapfrog firefly algorithm and application in di- 速度上都表现出更好的性能。本文只考虑了最优 spatch of power system containing wind farm[D].Shanghai: 个体对每个个体的影响，下一步工作是将个体邻域 East China University of Science and Technol ogy,2013. [13]符强，童楠，赵一鸣.一种基于多种群学习机制的萤火 的信息加入，进一步提高算法性能。 虫优化算法[J].计算机应用研究，2013,30(12)： 参考文献： 3600-3603 FU Qiang,TONG Nan,ZHAO Yiming.Firefly algorithm [1]YANG X S.Firefly algorithms for multimodal optimiz-ation based on multi-grouplearning mechanism[J].Application [C]//Stochastic Algorithms:Foundations and Applic- research of computers,2013,30(12):3600-3603. ations.Sapporo,Japan,2009,5792:169-178. [14]FISTER I,JR F,YANG X S,et al.A comprehensivereview of [2]朱书伟，周治平，张道文.融合并行混沌萤火虫算法的 firefly algorithms [J].Swarm and evolutionary computation, K-调和均值聚类[J].智能系统学报，2015,10(6)： 2013,13(1)：34-46. 872-880. [15]VERMA O P,AGGARWAL D,PATODI T.Opposition ZHU Shuxin,ZHOU Zhiping,ZHANG Daowen.Kharm- and dimensional based modified firefly algorithm[J].Ex- onic means clustering merged with parallel chaotic firefly pert-systems with applications an international journal, algorithm[J].CAAI transactions on intelligent systems, 2016,44(C):168-176. 2015,37(2):342-347 [16]汪慎文，丁立新，谢大同.应用精英反向学习策略的混合 [3]赵杰，雷秀娟，吴振强.基于最优类中心扰动的萤火虫聚 差分演化算法[J].武汉大学学报：理学版，2013,59 类算法[J].计算机工程与科学，2015,37(2)： (3):111-16. 342-347. WANG Shenwen,DING Lixin,XIE Datong.A hybrid ZHAO Jie,LEI Xiujuan,WU Zhengiang.A clustering al- differential evolution with elite opposition-based learning[]] gorithm for Fireflies based on optimal class center pert- Journal of Wuhan university:natural science edition,2013. urbation[J].Computer engineering science,2015,37 59(3)：111-116. (2):342-347. [17]程美英，倪志伟，朱旭辉.萤火虫优化算法理论研究综 「4]莫愿斌，马彦追，郑巧燕.单纯形法的改进萤火虫算法及 述[J].计算机科学，2015,42(4)：19-24. 其在非线性方程组求解中的应用[J].智能系统学报， CHENG Meiying,NI Zhiwei,ZHU Xuhui.Over-view on 2014,9(6):747-755. glow-worm swarm optimization or firefly algorithm J]. MO Yuanbin,MA Yanzhui,ZHENG Qiaoyan.Improved Computer science,2015,42(4):19-24. firefly algorithm based on simplex method and its appli- [18]汪慎文，丁立新，谢大同.应用反向学习策略的群搜索 cation in solving non-near equation groups J].CAAI tr- 优化算法[J].计算机科学，2012,39(9)：183-187. ansactions on intelligent systems,2014,9(6):747-755. WANG Shenwen,DING Lixin,XIE Datong.Group search [5 HORNG M H.Vector quantization using the firefly al- optimizer applying opposition-based learning[J].Computer gorithm for image compression J].Expert systems with science,2012,39(9):183-187. applications,2012,39(1):1078-1091. [19]周新宇，吴志健，王晖，等.一种精英反向学习的粒子 [6]MARICHELVAM M K,PRABAHARAN T,YANG XS.A 群优化算法[J].电子学报，2013,41(8)：1647-1652. discrete firefly algorithm for the multi-objective hy-brid flow ZHOU Xinyu,WU Zhijian,WANG Hui,et al.A particle shop scheduling problems [J ]IEEE transactio-ns on swarm optimization algorithm based on elite reverse learning evolutionary computation,2014,18(2):301-305. [J].Sinica Acta electronica,2013,41(8):1647-1652. [7]YANG X S,HOSSEINI SSS,GANDOMI A.Firefly [20]YANG XS.Firefly Algorithm Levy flights and global optimization algorithm for solving non-convex economic dispatch pr- [C]//Research and Development in Intelligent Systems oblems with valve loading effect J].Applied soft com XXVI.London,Springer,2010:209-218. puting,2012,12(3):1180-1186. [21]TILAHUN S L.HONG C O.Modified firefly algorithm J]. [8 SENTHILNATH J,OMKAR S,MANI V.Clustering using Joural of applied mathematics,2012,2012(12):2428-2439 firefly algorithm:performance study J].Swarm and 作者简介： evolutionary computation,2011,1(3):164-171. 魏伟一，男，1976年生，博士.副教 9]FALCON R.ALMEID M.NAYAK A.Fault identification with 授，CC℉会员，主要研究方向为智能信 binary adaptive fireflies in parallel and distrib-uted systems 息处理数字图像处理。 [C1//2011 IEEE Congress on Evolutionary Computation (CEC).New Orleans.USA.2011:1359-1366. [10]YU S,ZHU S,MA Y,et al.Enhancing firefly algorit- hmusing generalized opposition based learning[J].Com- puting,2015,97(7):741-754. [11]周凌云，丁立新，何进荣.精英正交学习萤火虫算法 文雅宏，男，1993年生，硕土研究 [J].计算机科学，2015.42(10)：211-216. 生，主要研究方向为数字图像处理、智 ZHOU Lingyun,DING Lixin,HE Jinrong.Elite-orth-ogonal 能计算。 learning firefly algorithm[J.Computer science,2015,42 (10):211-216. [12]李洋.蛙跳萤火虫算法及其在含风电场的电力系统调度 中的应用[D].上海：华东理工大学，2013
通过实验结果得出，ＥＯＦＡ 算法在解的精度和收敛 速度上都表现出更好的性能。 本文只考虑了最优 个体对每个个体的影响，下一步工作是将个体邻域 的信息加入，进一步提高算法性能。 参考文献： ［１］ＹＡＮＧ Ｘ Ｓ． Ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍｓ ｆｏｒ ｍｕｌｔｉｍｏｄａｌ ｏｐｔｉｍｉｚ⁃ａｔｉｏｎ ［ Ｃ ］ ／ ／ Ｓｔｏｃｈａｓｔｉｃ Ａｌｇｏｒｉｔｈｍｓ： Ｆｏｕｎｄａｔｉｏｎｓ ａｎｄ Ａｐｐｌｉｃ⁃ ａｔｉｏｎｓ． Ｓａｐｐｏｒｏ， Ｊａｐａｎ， ２００９， ５７９２： １６９－１７８． ［２］朱书伟，周治平，张道文． 融合并行混沌萤火虫算法的 Ｋ－调和均值聚类 ［ Ｊ］． 智能系统学报， ２０１５， １０ （ ６）： ８７２－８８０． ＺＨＵ Ｓｈｕｘｉｎ， ＺＨＯＵ Ｚｈｉｐｉｎｇ， ＺＨＡＮＧ Ｄａｏｗｅｎ． Ｋｈａｒｍ － ｏｎｉｃ ｍｅａｎｓ ｃｌｕｓｔｅｒｉｎｇ ｍｅｒｇｅｄ ｗｉｔｈ ｐａｒａｌｌｅｌ ｃｈａｏｔｉｃ ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ ［ Ｊ ］． ＣＡＡＩ ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ ｉｎｔｅｌｌｉｇｅｎｔ ｓｙｓｔｅｍｓ， ２０１５，３７（２）：３４２－３４７． ［３］赵杰， 雷秀娟，吴振强． 基于最优类中心扰动的萤火虫聚 类算 法 ［ Ｊ ］． 计 算 机 工 程 与 科 学， ２０１５， ３７ （ ２ ）： ３４２－３４７． ＺＨＡＯ Ｊｉｅ， ＬＥＩ Ｘｉｕｊｕａｎ， ＷＵ Ｚｈｅｎｑｉａｎｇ． Ａ ｃｌｕｓｔｅｒｉｎｇ ａｌ⁃ ｇｏｒｉｔｈｍ ｆｏｒ Ｆｉｒｅｆｌｉｅｓ ｂａｓｅｄ ｏｎ ｏｐｔｉｍａｌ ｃｌａｓｓ ｃｅｎｔｅｒ ｐｅｒｔ⁃ ｕｒｂａｔｉｏｎ［ Ｊ］． Ｃｏｍｐｕｔｅｒ ｅｎｇｉｎｅｅｒｉｎｇ ＆ ｓｃｉｅｎｃｅ， ２０１５， ３７ （２）： ３４２－３４７． ［４］莫愿斌， 马彦追， 郑巧燕．单纯形法的改进萤火虫算法及 其在非线性方程组求解中的应用［ Ｊ］． 智能系统学报， ２０１４， ９（６）： ７４７－７５５． ＭＯ Ｙｕａｎｂｉｎ， ＭＡ Ｙａｎｚｈｕｉ， ＺＨＥＮＧ Ｑｉａｏｙａｎ． Ｉｍｐｒｏｖｅｄ ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ ｂａｓｅｄ ｏｎ ｓｉｍｐｌｅｘ ｍｅｔｈｏｄ ａｎｄ ｉｔｓ ａｐｐｌｉ⁃ ｃａｔｉｏｎ ｉｎ ｓｏｌｖｉｎｇ ｎｏｎ⁃ｎｅａｒ ｅｑｕａｔｉｏｎ ｇｒｏｕｐｓ ［ Ｊ］． ＣＡＡＩ ｔｒ⁃ ａｎｓａｃｔｉｏｎｓ ｏｎ ｉｎｔｅｌｌｉｇｅｎｔ ｓｙｓｔｅｍｓ，２０１４， ９（６）： ７４７－７５５． ［５ ］ ＨＯＲＮＧ Ｍ Ｈ． Ｖｅｃｔｏｒ ｑｕａｎｔｉｚａｔｉｏｎ ｕｓｉｎｇ ｔｈｅ ｆｉｒｅｆｌｙ ａｌ⁃ ｇｏｒｉｔｈｍ ｆｏｒ ｉｍａｇｅ ｃｏｍｐｒｅｓｓｉｏｎ ［ Ｊ］． Ｅｘｐｅｒｔ ｓｙｓｔｅｍｓ ｗｉｔｈ ａｐｐｌｉｃａｔｉｏｎｓ， ２０１２， ３９（１）： １０７８－１０９１． ［６］ ＭＡＲＩＣＨＥＬＶＡＭ Ｍ Ｋ， ＰＲＡＢＡＨＡＲＡＮ Ｔ， ＹＡＮＧ ＸＳ． Ａ ｄｉｓｃｒｅｔｅ ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ ｆｏｒ ｔｈｅ ｍｕｌｔｉ⁃ｏｂｊｅｃｔｉｖｅ ｈｙ⁃ｂｒｉｄ ｆｌｏｗ ｓｈｏｐ ｓｃｈｅｄｕｌｉｎｇ ｐｒｏｂｌｅｍｓ ［ Ｊ ］． ＩＥＥＥ ｔｒａｎｓａｃｔｉｏ⁃ｎｓ ｏｎ ｅｖｏｌｕｔｉｏｎａｒｙ ｃｏｍｐｕｔａｔｉｏｎ， ２０１４， １８（２）： ３０１－３０５． ［７］ ＹＡＮＧ Ｘ Ｓ， ＨＯＳＳＥＩＮＩ Ｓ Ｓ Ｓ， ＧＡＮＤＯＭＩ Ａ． Ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ ｆｏｒ ｓｏｌｖｉｎｇ ｎｏｎ⁃ｃｏｎｖｅｘ ｅｃｏｎｏｍｉｃ ｄｉｓｐａｔｃｈ ｐｒ⁃ ｏｂｌｅｍｓ ｗｉｔｈ ｖａｌｖｅ ｌｏａｄｉｎｇ ｅｆｆｅｃｔ ［ Ｊ ］． Ａｐｐｌｉｅｄ ｓｏｆｔ ｃｏｍ⁃ ｐｕｔｉｎｇ， ２０１２， １２（３）： １１８０－１１８６． ［８］ ＳＥＮＴＨＩＬＮＡＴＨ Ｊ， ＯＭＫＡＲ Ｓ， ＭＡＮＩ Ｖ． Ｃｌｕｓｔｅｒｉｎｇ ｕｓｉｎｇ ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ： ｐｅｒｆｏｒｍａｎｃｅ ｓｔｕｄｙ ［ Ｊ ］． Ｓｗａｒｍ ａｎｄ ｅｖｏｌｕｔｉｏｎａｒｙ ｃｏｍｐｕｔａｔｉｏｎ， ２０１１， １（３）： １６４－１７１． ［９］ＦＡＬＣＯＮ Ｒ， ＡＬＭＥＩＤ Ｍ， ＮＡＹＡＫ Ａ． Ｆａｕｌｔ ｉｄｅｎｔｉｆｉｃａｔｉｏｎ ｗｉｔｈ ｂｉｎａｒｙ ａｄａｐｔｉｖｅ ｆｉｒｅｆｌｉｅｓ ｉｎ ｐａｒａｌｌｅｌ ａｎｄ ｄｉｓｔｒｉｂ⁃ｕｔｅｄ ｓｙｓｔｅｍｓ ［ Ｃ ］ ／ ／ ２０１１ ＩＥＥＥ Ｃｏｎｇｒｅｓｓ ｏｎ Ｅｖｏｌｕｔｉｏｎａｒｙ Ｃｏｍｐｕｔａｔｉｏｎ （ＣＥＣ）．Ｎｅｗ Ｏｒｌｅａｎｓ， ＵＳＡ， ２０１１： １３５９－１３６６． ［１０］ ＹＵ Ｓ， ＺＨＵ Ｓ， ＭＡ Ｙ， ｅｔ ａｌ． Ｅｎｈａｎｃｉｎｇ ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔ⁃ ｈｍｕｓｉｎｇ ｇｅｎｅｒａｌｉｚｅｄ ｏｐｐｏｓｉｔｉｏｎ ｂａｓｅｄ ｌｅａｒｎｉｎｇ ［ Ｊ］． Ｃｏｍ⁃ ｐｕｔｉｎｇ， ２０１５， ９７（７）： ７４１－７５４． ［１１］周凌云， 丁立新， 何进荣． 精英正交学习萤火虫算法 ［Ｊ］． 计算机科学， ２０１５， ４２（１０）： ２１１－２１６． ＺＨＯＵ Ｌｉｎｇｙｕｎ，ＤＩＮＧ Ｌｉｘｉｎ， ＨＥ Ｊｉｎｒｏｎｇ． Ｅｌｉｔｅ⁃ｏｒｔｈ⁃ｏｇｏｎａｌ ｌｅａｒｎｉｎｇ ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ［ Ｊ］．Ｃｏｍｐｕｔｅｒ ｓｃｉｅｎｃｅ， ２０１５， ４２ （１０）： ２１１－２１６． ［１２］李洋．蛙跳萤火虫算法及其在含风电场的电力系统调度 中的应用［Ｄ］． 上海： 华东理工大学，２０１３． ＬＩ Ｙａｎｇ． Ｌｅａｐｆｒｏｇ ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ ａｎｄ ａｐｐｌｉｃａｔｉｏｎ ｉｎ ｄｉ⁃ ｓｐａｔｃｈ ｏｆ ｐｏｗｅｒ ｓｙｓｔｅｍ ｃｏｎｔａｉｎｉｎｇ ｗｉｎｄ ｆａｒｍ［Ｄ］． Ｓｈａｎｇｈａｉ： Ｅａｓｔ Ｃｈｉｎａ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｓｃｉｅｎｃｅ ａｎｄ Ｔｅｃｈｎｏｌ ｏｇｙ，２０１３． ［１３］符强，童楠， 赵一鸣．一种基于多种群学习机制的萤火 虫优化算法［ Ｊ］． 计算机应用研究， ２０１３， ３０ （ １２）： ３６００－３６０３． ＦＵ Ｑｉａｎｇ， ＴＯＮＧ Ｎａｎ， ＺＨＡＯ Ｙｉｍｉｎｇ． Ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ ｂａｓｅｄ ｏｎ ｍｕｌｔｉ⁃ｇｒｏｕｐｌｅａｒｎｉｎｇ ｍｅｃｈａｎｉｓｍ［ Ｊ］． Ａｐｐｌｉｃａｔｉｏｎ ｒｅｓｅａｒｃｈ ｏｆ ｃｏｍｐｕｔｅｒｓ， ２０１３， ３０（１２）： ３６００－３６０３． ［１４］ＦＩＳＴＥＲ Ｉ， ＪＲ Ｆ， ＹＡＮＧ Ｘ Ｓ， ｅｔ ａｌ． Ａ ｃｏｍｐｒｅｈｅｎｓｉｖｅｒｅｖｉｅｗ ｏｆ ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍｓ ［ Ｊ］． Ｓｗａｒｍ ａｎｄ ｅｖｏｌｕｔｉｏｎａｒｙ ｃｏｍｐｕｔａｔｉｏｎ， ２０１３， １３（１）： ３４－４６． ［１５］ ＶＥＲＭＡ Ｏ Ｐ， ＡＧＧＡＲＷＡＬ Ｄ， ＰＡＴＯＤＩ Ｔ． Ｏｐｐｏｓｉｔｉｏｎ ａｎｄ ｄｉｍｅｎｓｉｏｎａｌ ｂａｓｅｄ ｍｏｄｉｆｉｅｄ ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ［ Ｊ］． Ｅｘ⁃ ｐｅｒｔ⁃ｓｙｓｔｅｍｓ ｗｉｔｈ ａｐｐｌｉｃａｔｉｏｎｓ ａｎ ｉｎｔｅｒｎａｔｉｏｎａｌ ｊｏｕｒｎａｌ， ２０１６， ４４（Ｃ）： １６８－１７６． ［１６］汪慎文，丁立新，谢大同．应用精英反向学习策略的混合 差分演化算法［ Ｊ］．武汉大学学报： 理学版， ２０１３， ５９ （３）： １１１－１６． ＷＡＮＧ Ｓｈｅｎｗｅｎ， ＤＩＮＧ Ｌｉｘｉｎ， ＸＩＥ Ｄａｔｏｎｇ． Ａ ｈｙｂｒｉｄ ｄｉｆｆｅｒｅｎｔｉａｌ ｅｖｏｌｕｔｉｏｎ ｗｉｔｈ ｅｌｉｔｅ ｏｐｐｏｓｉｔｉｏｎ⁃ｂａｓｅｄ ｌｅａｒｎｉｎｇ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｗｕｈａｎ ｕｎｉｖｅｒｓｉｔｙ： ｎａｔｕｒａｌ ｓｃｉｅｎｃｅ ｅｄｉｔｉｏｎ， ２０１３， ５９（３）： １１１－１１６． ［１７］程美英， 倪志伟， 朱旭辉． 萤火虫优化算法理论研究综 述［Ｊ］． 计算机科学， ２０１５， ４２（４）： １９－２４． ＣＨＥＮＧ Ｍｅｉｙｉｎｇ， ＮＩ Ｚｈｉｗｅｉ， ＺＨＵ Ｘｕｈｕｉ． Ｏｖｅｒ⁃ｖｉｅｗ ｏｎ ｇｌｏｗ⁃ｗｏｒｍ ｓｗａｒｍ ｏｐｔｉｍｉｚａｔｉｏｎ ｏｒ ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ ［ Ｊ ］． Ｃｏｍｐｕｔｅｒ ｓｃｉｅｎｃｅ， ２０１５， ４２（４）： １９－２４． ［１８］汪慎文， 丁立新， 谢大同． 应用反向学习策略的群搜索 优化算法［Ｊ］． 计算机科学， ２０１２， ３９（９）： １８３－１８７． ＷＡＮＧ Ｓｈｅｎｗｅｎ，ＤＩＮＧ Ｌｉｘｉｎ， ＸＩＥ Ｄａｔｏｎｇ． Ｇｒｏｕｐ ｓｅａｒｃｈ ｏｐｔｉｍｉｚｅｒ ａｐｐｌｙｉｎｇ ｏｐｐｏｓｉｔｉｏｎ⁃ｂａｓｅｄ ｌｅａｒｎｉｎｇ［Ｊ］． Ｃｏｍｐｕｔｅｒ ｓｃｉｅｎｃｅ， ２０１２， ３９（９）： １８３－１８７． ［１９］周新宇， 吴志健， 王晖，等． 一种精英反向学习的粒子 群优化算法［Ｊ］． 电子学报， ２０１３， ４１（８）： １６４７－１６５２． ＺＨＯＵ Ｘｉｎｙｕ， ＷＵ Ｚｈｉｊｉａｎ， ＷＡＮＧ Ｈｕｉ， ｅｔ ａｌ． Ａ ｐａｒｔｉｃｌｅ ｓｗａｒｍ ｏｐｔｉｍｉｚａｔｉｏｎ ａｌｇｏｒｉｔｈｍ ｂａｓｅｄ ｏｎ ｅｌｉｔｅ ｒｅｖｅｒｓｅ ｌｅａｒｎｉｎｇ ［Ｊ］．Ｓｉｎｉｃａ Ａｃｔａ ｅｌｅｃｔｒｏｎｉｃａ， ２０１３， ４１（８）：１６４７－１６５２． ［２０］ＹＡＮＧ ＸＳ． Ｆｉｒｅｆｌｙ Ａｌｇｏｒｉｔｈｍ Ｌéｖｙ ｆｌｉｇｈｔｓ ａｎｄ ｇｌｏｂａｌ ｏｐｔｉｍｉｚａｔｉｏｎ ［Ｃ］ ／ ／ Ｒｅｓｅａｒｃｈ ａｎｄ Ｄｅｖｅｌｏｐｍｅｎｔ ｉｎ Ｉｎｔｅｌｌｉｇｅｎｔ Ｓｙｓｔｅｍｓ ＸＸＶＩ． Ｌｏｎｄｏｎ， Ｓｐｒｉｎｇｅｒ， ２０１０：２０９－２１８． ［２１］ ＴＩＬＡＨＵＮ Ｓ Ｌ， ＨＯＮＧ Ｃ Ｏ． Ｍｏｄｉｆｉｅｄ ｆｉｒｅｆｌｙ ａｌｇｏｒｉｔｈｍ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ ａｐｐｌｉｅｄ ｍａｔｈｅｍａｔｉｃｓ， ２０１２， ２０１２（１２）： ２４２８－２４３９． 作者简介： 魏伟一 ，男，１９７６ 年生，博士，副教 授，ＣＣＦ 会员，主要研究方向为智能信 息处理、数字图像处理。 文雅宏 ，男，１９９３ 年生，硕士研究 生，主要研究方向为数字图像处理、智 能计算。 ·７１６· 智 能 系 统 学 报 第 １２ 卷