第十章-群集智能与蚁群算法

第十章

群集智能与蚁群算法2013年6月1第1节群体智能概念

1.群体智能(SwarmIntelligence)2研究意义－群体智能概念

群体智能这个概念来自对蜜蜂和蚂蚁的观察。一组相互之间可以进行直接通信或者间接通信(通过改变局部环境)的主体，这组主体能够合作进行分布问题求解。任何启发于群居性昆虫群体和其它动物群体的集体行为而设计的算法和分布式问题解决装置都称为群体智能。群体智能在没有集中控制并且不提供全局模型的前提下，为寻找复杂的分布式问题的解决方案提供了基础。3研究意义－群体智能的特点

分布式：能够适应当前网络环境下的工作状态;鲁棒性：没有中心的控制与数据，个体的故障不影响整个问题的求解;扩充性：个体的增加，系统的通信开销增加小;简单性：个体简单，实现也比较简单。

4研究意义群体智能的研究不仅在多主体仿真、系统复杂性以及NP问题等方面为人工智能、认识科学、计算经济学等领域的基础理论问题的研究开辟了新的研究途径，同时也为诸如组合优化、机器人协作、电信路由控制等实际工程问题提供了新的解决方法。因此，群体智能的研究具有重要意义和广阔的应用前景。5研究现状－国外

美国的SDG组织在系统复杂性方面开展了研究。他们主要通过多主体的仿真来研究系统复杂性。他们开发的SWARM软件包为多学科进行基于多主体的建模提供了一个基础平台；加州工学院专门开设了群体智能的课程；

欧洲联盟资助的SWARM-BOTS项目的主要目标是研究设计和实现自组织和自装配的装置的新途径。它的理论基础是群体智能和蚁群算法的近期研究成果，即对群居性昆虫和其它动物群体的自组织和自装配能力的研究。6研究现状－国内

国家自然科学基金在“十五”期间学科交叉类优先资助领域第一类属于把握科学前沿，推动源头创新的项目。其中第7项认知科学及其信息处理的研究内容就明确列出了群体智能的进化、自适应与现场认知。相关项目还有第9项复杂系统与复杂性。

7主要研究方向蚁群寻食行为研究，相对应组合优化算法和通信网络路由控制算法；群体分工和任务分配行为研究，相对应多主体分工协作算法；巢穴组织和自组织行为及群体分类行为研究，相对应数据分析和图的分割算法；建巢和自装配行为研究，相对应模拟建巢算法；群体合作搬运行为研究，相对应机器人合作搬运算法。8所需解决的关键问题蚁群算法

效率与理论；由于没有标准的测试集，除了寻食模型，蚁卵聚类、蚁群分工和蚁巢自装配等模型都只处于证实阶段理论和实验

；一个多主体自组织模型实验和测试平台

；对于追求效率的实际问题，如何既保持群体智能系统的灵活性和鲁棒性等自组织特征又能保证系统的高效率也是一个关键问题

；群体智能与分布式智能的智能主体研究相结合，将产生新的智能主体协作、建模等算法和机制，提出网络和网格环境的自适应多智能主体系统

。9蚂蚁寻找最短路径原理A）蚁群到达决策点。B）一些蚂蚁选择上方路径，一些蚂蚁选择下方路径。选择是随机的。C）下方短路径蚂蚁到达相反方向的决策点的时间早于选择上方长路径的蚂蚁。D）短路径上外激素以较高的速度积累。。外激素多的短路径将吸收更多的蚂蚁，反过来，更多的蚂蚁在短路径上会留下更多的外激素，加上外激素挥发效应，最后，蚁群都选择了最短路径。

102.蚁群优化算法研究现状1/7

90年代Dorigo最早提出了蚁群优化算法---蚂蚁系统（AntSystem,AS）并将其应用于解决计算机算法学中经典的旅行商问题（TSP）。从蚂蚁系统开始，基本的蚁群算法得到了不断的发展和完善，并在TSP以及许多实际优化问题求解中进一步得到了验证。这些AS改进版本的一个共同点就是增强了蚂蚁搜索过程中对最优解的探索能力，它们之间的差异仅在于搜索控制策略方面。而且，取得了最佳结果的ACO是通过引入局部搜索算法实现的，这实际上是一些结合了标准局域搜索算法的混合型概率搜索算法，有利于提高蚁群各级系统在优化问题中的求解质量。11蚁群优化算法研究现状2/7最初提出的AS有三种版本：Ant-density、Ant-quantity和Ant-cycle。在Ant-density和Ant-quantity中蚂蚁在两个位置节点间每移动一次后即更新信息素，而在Ant-cycle中当所有的蚂蚁都完成了自己的行程后才对信息素进行更新，而且每个蚂蚁所释放的信息素被表达为反映相应行程质量的函数。通过与其它各种通用的启发式算法相比，在不大于75城市的TSP中，这三种基本算法的求解能力还是比较理想的，但是当问题规模扩展时，AS的解题能力大幅度下降。因此，其后的ACO研究工作主要都集中于AS性能的改进方面。较早的一种改进方法是精英策略(ElitistStrategy)，其思想是在算法开始后即对所有已发现的最好路径给予额外的增强，并将随后与之对应的行程记为Tgb(全局最优行程)，当进行信息素更新时，对这些行程予以加权，同时将经过这些行程的蚂蚁记为“精英”，从而增大较好行程的选择机会。这种改进型算法能够以更快的速度获得更好的解。但是若选择的精英过多则算法会由于较早的收敛于局部次优解而导致搜索的过早停滞。

12蚁群优化算法研究现状3/7

为了进一步克服AS中暴露出的问题，提出了蚁群系统(AntColonySystem,ACS)。该系统的提出是以Ant-Q算法为基础的。Ant-Q将蚂蚁算法和一种增强型学习算法Q-learning有机的结合了起来。ACS与AS之间存在三方面的主要差异：首先，ACS采用了更为大胆的行为选择规则；其次，只增强属于全局最优解的路径上的信息素。其中，0<ρ<1是信息素挥发参数，是从寻路开始到当前为止全局最优的路径长度。13蚁群优化算法研究现状4/7

再次，还引入了负反馈机制，每当一只蚂蚁由一个节点移动到另一个节点时，该路径上的信息素都按照如下公式被相应的消除一部分，从而实现一种信息素的局部调整，以减小已选择过的路径再次被选择的概率。

14蚁群优化算法研究现状5/7在对AS进行直接完善的方法中，MAX-MINAntSystem是一个典型代表。该算法修改了AS的信息素更新方式，每次迭代之后只有一只蚂蚁能够进行信息素的更新以获取更好的解。为了避免搜索停滞，路径上的信息素浓度被限制在[MAX，MIN]范围内，另外，信息素的初始值被设为其取值上限，这样有助于增加算法初始阶段的搜索能力。15蚁群优化算法研究现状6/7

另一种对AS改进的算法是Rank-basedVersionAS。与“精英策略”相似，在此算法中总是更新更好进程上的信息素，选择的标准是其行程长度决定的排序，且每个蚂蚁放置信息素的强度通过下式中的排序加权处理确定，其中，ｗ为每次迭代后放置信息素的蚂蚁总数。

16蚁群优化算法研究现状7/7这种算法求解TSP的能力与AS、精英策略AS、遗传算法和模拟退火算法进行了比较。在大型TSP问题中（最多包含132座城市），基于AS的算法都显示出了优于GA和SA的特性。而且在Rank-basedAS和精英策略AS均优于基本AS的同时，前者还获得了比精英策略AS更好的解。

173.蚁群优化算法应用现状1/5随着群智能理论和应用算法研究的不断发展，研究者已尝试着将其用于各种工程优化问题，并取得了意想不到的收获。多种研究表明，群智能在离散求解空间和连续求解空间中均表现出良好的搜索效果，并在组合优化问题中表现突出。蚁群优化算法并不是旅行商问题的最佳解决方法，但是它却为解决组合优化问题提供了新思路，并很快被应用到其它组合优化问题中。比较典型的应用研究包括：网络路由优化、数据挖掘以及一些经典的组合优化问题。

18蚁群优化算法应用现状2/5蚁群算法在电信路由优化中已取得了一定的应用成果。HP公司和英国电信公司在90年代中后期都开展了这方面的研究，设计了蚁群路由算法（AntColonyRouting,ACR）。每只蚂蚁就像蚁群优化算法中一样，根据它在网络上的经验与性能，动态更新路由表项。如果一只蚂蚁因为经过了网络中堵塞的路由而导致了比较大的延迟，那么就对该表项做较大的增强。同时根据信息素挥发机制实现系统的信息更新，从而抛弃过期的路由信息。这样，在当前最优路由出现拥堵现象时，ACR算法就能迅速的搜寻另一条可替代的最优路径，从而提高网络的均衡性、负荷量和利用率。目前这方面的应用研究仍在升温，因为通信网络的分布式信息结构、非稳定随机动态特性以及网络状态的异步演化与ACO的算法本质和特性非常相似。19

蚁群优化算法应用现状3/5

基于群智能的聚类算法起源于对蚁群蚁卵的分类研究。Lumer和Faieta将Deneubourg提出将蚁巢分类模型应用于数据聚类分析。其基本思想是将待聚类数据随机地散布到一个二维平面内，然后将虚拟蚂蚁分布到这个空间内，并以随机方式移动，当一只蚂蚁遇到一个待聚类数据时即将之拾起并继续随机运动，若运动路径附近的数据与背负的数据相似性高于设置的标准则将其放置在该位置，然后继续移动，重复上述数据搬运过程。按照这样的方法可实现对相似数据的聚类。20

蚁群优化算法应用现状4/5

ACO还在许多经典组合优化问题中获得了成功的应用，如二次规划问题（QAP）、机器人路径规划、作业流程规划、图着色（GraphColoring）等问题。经过多年的发展，ACO已成为能够有效解决实际二次规划问题的几种重要算法之一。AS在作业流程计划（Job-shopScheduling）问题中的应用实例已经出现，这说明了AS在此领域的应用潜力。利用MAX-MINAS解决PAQ也取得了比较理想的效果，并通过实验中的计算数据证明采用该方法处理PAQ比较早的SA算法更好，且与禁忌搜索算法性能相当。利用ACO实现对生产流程和特料管理的综合优化，并通过与遗传、模拟退火和禁忌搜索算法的比较证明了ACO的工程应用价值。21蚁群优化算法应用现状5/5许多研究者将ACO用于了武器攻击目标分配和优化问题、车辆运行路径规划、区域性无线电频率自动分配、Bayesiannetworks的训练和集合覆盖等应用优化问题。Costa和Herz还提出了一种AS在规划问题方面的扩展应用——图着色问题，并取得了可与其他启发式算法相比的效果。22第2节蚁群优化算法概念 2.1蚁群算法原理2.2简化的蚂蚁寻食过程2.3自然蚁群与人工蚁群算法2.4蚁群算法与TSP问题2.5初始的蚁群优化算法—基于图的蚁群系统（GBAS）2.6一般蚁群算法的框架2.7蚁群优化算法—技术问题232.1蚁群算法原理

蚁群算法是对自然界蚂蚁的寻径方式进行模似而得出的一种仿生算法。蚂蚁在运动过程中，能够在它所经过的路径上留下一种称之为外激素(pheromone)的物质进行信息传递，而且蚂蚁在运动过程中能够感知这种物质，并以此指导自己的运动方向，因此由大量蚂蚁组成的蚁群集体行为便表现出一种信息正反馈现象：某一路径上走过的蚂蚁越多，则后来者选择该路径的概率就越大。为了说明蚁群算法的原理，先简要介绍一下蚂蚁搜寻食物的具体过程。在蚁群寻找食物时，它们总能找到一条从食物到巢穴之间的最优路径。这是因为蚂蚁在寻找路径时会在路径上释放出一种特殊的信息素。当它们碰到一个还没有走过的路口时．就随机地挑选一条路径前行。与此同时释放出与路径长度有关的信息素。路径越长，释放的激索浓度越低.当后来的蚂蚁再次碰到这个路口的时候．选择激素浓度较高路径概率就会相对较大。这样形成一个正反馈。最优路径上的激索浓度越来越大．而其它的路径上激素浓度却会随着时间的流逝而消减。最终整个蚁群会找出最优路径。242.2简化的蚂蚁寻食过程1/3蚂蚁从A点出发，速度相同，食物在D点，可能随机选择路线ABD或ACD。假设初始时每条分配路线一只蚂蚁，每个时间单位行走一步，本图为经过9个时间单位时的情形：走ABD的蚂蚁到达终点，而走ACD的蚂蚁刚好走到C点，为一半路程。25简化的蚂蚁寻食过程2/3本图为从开始算起，经过18个时间单位时的情形：走ABD的蚂蚁到达终点后得到食物又返回了起点A，而走ACD的蚂蚁刚好走到D点。26

简化的蚂蚁寻食过程3/3

假设蚂蚁每经过一处所留下的信息素为一个单位，则经过36个时间单位后，所有开始一起出发的蚂蚁都经过不同路径从D点取得了食物，此时ABD的路线往返了2趟，每一处的信息素为4个单位，而ACD的路线往返了一趟，每一处的信息素为2个单位，其比值为2：1。寻找食物的过程继续进行，则按信息素的指导，蚁群在ABD路线上增派一只蚂蚁（共2只），而ACD路线上仍然为一只蚂蚁。再经过36个时间单位后，两条线路上的信息素单位积累为12和4，比值为3：1。若按以上规则继续，蚁群在ABD路线上再增派一只蚂蚁（共3只），而ACD路线上仍然为一只蚂蚁。再经过36个时间单位后，两条线路上的信息素单位积累为24和6，比值为4：1。若继续进行，则按信息素的指导，最终所有的蚂蚁会放弃ACD路线，而都选择ABD路线。这也就是前面所提到的正反馈效应。272.3自然蚁群与人工蚁群算法基于以上蚁群寻找食物时的最优路径选择问题，可以构造人工蚁群，来解决最优化问题，如TSP问题。人工蚁群中把具有简单功能的工作单元看作蚂蚁。二者的相似之处在于都是优先选择信息素浓度大的路径。较短路径的信息素浓度高，所以能够最终被所有蚂蚁选择，也就是最终的优化结果。两者的区别在于人工蚁群有一定的记忆能力，能够记忆已经访问过的节点。同时，人工蚁群再选择下一条路径的时候是按一定算法规律有意识地寻找最短路径，而不是盲目的。例如在TSP问题中，可以预先知道当前城市到下一个目的地的距离。282.4蚁群算法与TSP问题1/3TSP问题表示为一个N个城市的有向图G=（N，A），其中 城市之间距离目标函数为，其中为城市1,2，…n的一个排列，。29

蚁群算法与TSP问题2/3

TSP问题的人工蚁群算法中，假设m只蚂蚁在图的相邻节点间移动，从而协作异步地得到问题的解。每只蚂蚁的一步转移概率由图中的每条边上的两类参数决定：1信息素值也称信息素痕迹。2可见度，即先验值。信息素的更新方式有2种，一是挥发，也就是所有路径上的信息素以一定的比率进行减少，模拟自然蚁群的信息素随时间挥发的过程；二是增强，给评价值“好”(有蚂蚁走过)的边增加信息素。30蚁群算法与TSP问题3/3蚂蚁向下一个目标的运动是通过一个随机原则来实现的，也就是运用当前所在节点存储的信息，计算出下一步可达节点的概率，并按此概率实现一步移动，逐此往复，越来越接近最优解。蚂蚁在寻找过程中，或者找到一个解后，会评估该解或解的一部分的优化程度，并把评价信息保存在相关连接的信息素中。312.5初始的蚁群优化算法—基于图的蚁群系统（GBAS）1/12初始的蚁群算法是基于图的蚁群算法，graph-basedantsystem,简称为GBAS，是由GutjahrWJ在2000年的FutureGenerationComputingSystems提出的，课本的参考文献2。算法步骤如下：STEP0

对n个城市的TSP问题，城市间的距离矩阵为，给TSP图中的每一条弧赋信息素初值，假设m只蚂蚁在工作，所有蚂蚁都从同一城市出发。当前最好解是 。32初始的蚁群优化算法—基于图的蚁群系统（GBAS）2/12STEP1

（外循环）如果满足算法的停止规则，则停止计算并输出计算得到的最好解。否则使蚂蚁s从起点出发，用表示蚂蚁s行走的城市集合，初始为空集，。STEP2(内循环)按蚂蚁的顺序分别计算。当蚂蚁在城市i，若 完成第s只蚂蚁的计算。否则，若，则以概率 ， 到达j， ；若则到达 重复STEP2。33初始的蚁群优化算法—基于图的蚁群系统（GBAS）3/12STRP3

对 ，若，按中城市的顺序计算路径程度；若，路径长度置为一个无穷大值（即不可达）。比较m只蚂蚁中的路径长度，记走最短路径的蚂蚁为t。若，则。用如下公式对W路径上的信息素痕迹加强，对其他路径上的信息素进行挥发。得到新的，重复步骤STEP1。34初始的蚁群优化算法—基于图的蚁群系统（GBAS）4/12在STEP3

中，挥发因子对于一个固定的，满足并且

经过k次挥发，非最优路径的信息素逐渐减少至消失。35初始的蚁群优化算法—基于图的蚁群系统（GBAS）5/12以上算法中，在蚂蚁的搜寻过程中，以信息素的概率分布来决定从城市i到城市j的转移。算法中包括信息素更新的过程

1信息素挥发（evaporation）信息素痕迹的挥发过程是每个连接上的信息素痕迹的浓度自动逐渐减弱的过程，由表示，这个挥发过程主要用于避免算法过快地向局部最优区域集中，有助于搜索区域的扩展。

2信息素增强（reinforcement）增强过程是蚁群优化算法中可选的部分，称为离线更新方式（还有在线更新方式）。这种方式可以实现由单个蚂蚁无法实现的集中行动。也就是说，增强过程体现在观察蚁群（m只蚂蚁）中每只蚂蚁所找到的路径，并选择其中最优路径上的弧进行信息素的增强，挥发过程是所有弧都进行的，不于蚂蚁数量相关。这种增强过程中进行的信息素更新称为离线的信息素更新。在STEP3中，蚁群永远记忆到目前为止的最优解。36图的蚁群系统（GBAS）6/12可以验证，下式满足：即是一个随机矩阵。四个城市的非对称TSP问题，距离矩阵和城市图示如下：37

初始的蚁群优化算法—基于图的蚁群系统（GBAS）7/12假设共4只蚂蚁，所有蚂蚁都从城市A出发，挥发因子。此时，观察GBAS的计算过程。矩阵共有12条弧，初始信息素记忆矩阵为：38初始的蚁群优化算法—基于图的蚁群系统（GBAS）8/12执行GBAS算法的步骤2，假设蚂蚁的行走路线分别为：当前最优解为，这个解是截止到当前的最优解，碰巧是实际最优解39初始的蚁群优化算法—基于图的蚁群系统（GBAS）9/12按算法步骤3的信息素更新规则，得到更新矩阵这是第一次外循环结束的状态。40初始的蚁群优化算法—基于图的蚁群系统（GBAS）10/12重复外循环，由于上一次得到的W2已经是全局最优解，因此按算法步骤3的信息素更新规则，无论蚂蚁如何行走，都只是对W2路线上的城市信息素进行增强，其他的城市信息素进行挥发。得到更新矩阵这是第一次外循环结束的状态。41初始的蚁群优化算法—基于图的蚁群系统（GBAS）11/12重复外循环，由于W2全局最优解，GBAS只记录第一个最优解，因此一但得到了全局最优解，信息素的更新将不再依赖于以群的行走路线，而只是不断增强最优路线的信息素，同时进行挥发。第三次外循环后得到的信息素矩阵为：42初始的蚁群优化算法—基于图的蚁群系统（GBAS）12/12蚂蚁以一定的概率从城市i到城市j进行转移，信息素的更新在STEP3完成，并随K而变化。假设第K次外循环后得到信息素矩阵，得到当前最优解。第K次循环前的信息素和最优解为，经过第K次外循环后，得到。由于蚂蚁的一步转移概率是随机的，从到也是随机的，是一个马尔可夫过程。432.6一般蚁群算法的框架一般蚁群算法的框架和GBAS基本相同，有三个组成部分：蚁群的活动；信息素的挥发；信息素的增强；主要体现在前面的算法中步骤2和步骤3中的转移概率公式和信息素更新公式。442.7蚁群优化算法—技术问题（1）解的表达形式与算法的实现（2）每一节点的记忆信息和系数的确定（3）蚁群的规模和停止规则（4）信息素的更改45（1）

解的表达形式与算法的实现1/4

----解的表达形式解的表达形式基于TSP问题的蚁群优化算法，其解的形式是所有城市的一个排列（闭圈，这种情况下谁在第一并不重要），信息素痕迹按每个弧记录。而对于一般以顺序作为解的优化问题，谁在第一是很重要的。蚁群算法在解决这类问题时，只需要建立一个虚拟的始终点，就可以把TSP问题的解法推广，用于诸多的优化问题。诸如车间作业及下料等问题，他们的共同特点是解以一个顺序表示。TSP问题寻找的是最短回路，而一般优化问题中，STEP3中的判断条件需要根据实际问题进行修改。46解的表达形式与算法的实现2/4

----算法的实现例：0-1背包问题的解顺序表达形式与算法实现。设有一个容积为b的背包，n个尺寸分别为，价值分别为的物品，0-1背包问题的数学模型为：假设其解的顺序表达形式为，其中为的一个排列。47

解的表达形式与算法的实现3/4

----算法的实现建立有向图，其中A中共有条弧。初始信息素痕迹定义为。设第s只蚂蚁第k步所走的路线为，表示蚂蚁从0点出发，顺序到达。第步按TSP算法的转移概率公式行走选择。若则，否则，此蚂蚁不再继续行走，退回起点。48解的表达形式与算法的实现4/4----算法的实现对蚁群重复以上过程，比较m只蚂蚁的装包值 并记忆具有最大装包值的蚂蚁为t。把GBAS算法中步骤3中的改为，若满足此条件则替换当前最好解为，对W上的弧进行信息素的加强，其他弧进行信息素的挥发。算法中记录了三个信息：信息素痕迹；行走路线；和问题的约束条件 ，以确定是否将加入。49（2）每一节点的记忆信息和系数的确定----需要记忆的信息1/3算法中需要记忆的信息有三部分。第一部分信息是存在每个节点的路由表数据结构，由此决定的的转移概率为其中T可以看成节点i的邻域。50每一节点的记忆信息和系数的确定----需要记忆的信息2/3第二部分需要记忆的信息是每个蚂蚁的记忆表中存储着的自身的历史信息，这一部分主要由算法的中的记忆，表示蚂蚁已经行走过的节点。第三部分为问题的约束条件。在GBAS中，T集合表示满足约束条件的候选集，在背包问题的蚁群算法中由判别条件， 来实现记

忆功能。51每一节点的记忆信息和系数的确定----系数的确定3/3

残留信息的相对重要程度和预见值的相对重要程度体现了相关信息痕迹和预见度对蚂蚁决策的相对影响。Dorigo在求解TSP问题时，推荐参数的最佳设置为：。52（3）蚁群的规模和停止规则一、蚁群大小一般情况下蚁群中蚂蚁的个数不超过TSP图中节点的个数。二、终止条件

1给定一个外循环的最大数目，表明已经有足够的蚂蚁工作；

2当前最优解连续K次相同而停止，其中K是一个给定的整数，表示算法已经收敛，不再需要继续；

3目标值控制规则，给定优化问题（目标最小化）的一个下界和一个误差值，当算法得到的目标值同下界之差小于给定的误差值时，算法终止。53（4）信息素的更改1/6信息素的更新分为离线和在线两种方式。离线方式（同步更新方式）的主要思想是在若干只蚂蚁完成n个城市的访问后，统一对残留信息进行更新处理。信息素的在线更新（异步更新方式）即蚂蚁每行走一步，立即回溯并且更新行走路径上的信息素。54信息素的更改2/6离线方式的信息素更新可以进一步分为单蚂蚁离线更新和蚁群离线更新。蚁群离线更新方式是在蚁群中的m只蚂蚁全部完成n城市的访问（第k-1次蚁群循环）后，统一对残留信息进行更新处理。其中，为第k-1次循环后的的信息素的痕迹值。单蚂蚁离线更新是在第s只蚂蚁完成对所有n个城市的访问后，进行路径回溯，更新行走路径上的信息素，同时释放分配给它的资源。更新公式为第s+1只蚂蚁根据重新计算路由表。

55

信息素的更改3/6TSP问题中，蚁群优化算法根据信息素痕迹更新方式不同可以分为不同的算法，采用离线方式，并且时，其中W为t循环中m只蚂蚁所行走的最佳路线或第t只蚂蚁所行走的一条路径。Q为一个常数，该算法名为蚁环算法（ant-cyclealgotithm）,特点是行走的路径越短对应保存的信息素的值就越大。56信息素的更改4/6

GBAS算法是典型的离线信息素更新方式。该算法中，蚁群中蚂蚁的先后出行顺序没有相关性，但是每次循环需要记忆m只蚂蚁的行走路径，以进行比较选择最优路径。相对而言，单蚂蚁离线更新方式记忆信息少，只需要记忆第s只蚂蚁的路径，并通过信息素更新后，释放该蚂蚁的所有记录信息。实际上这种方式等价于蚁群离线方式中只有一只蚂蚁。57

信息素的更改5/6与单蚂蚁离线更新方式相比，信息量记忆更小的是信息素在线更新方式，即蚂蚁每走一步，马上回溯并且更新刚刚走过的路径上的信息素，其规则为其中，k为蚂蚁行走的第k步。58信息素的更改6/6

蚁量算法（ant-quantityalgorithm）的信息素更新为,Q为常量，表示i到j的距离，这样信息浓度会随城市距离的减小而加大。蚁密算法（ant-densityalgorithm）信息素更新为。以上三种算法中，蚁环算法效果最好，因为他用的是全局信息，而其余两种算法用的是局部信息。蚁环离线更新方法很好地保证了残留信息不至于无限积累，非最优路径会逐渐随时间推移被忘记。59

蚁群优化算法—参考书1智能蚁群算法及应用吴启迪上海科技出版社从基本结构、算法特点、改进方法、突破途径、实现模式及应用模式等方面进行了论述。主要内容有蚁群算法的由来、研究成果、应用综述、算法的具体描述及改进、算法的典型优化问题求解模式、算法的典型应用及拓展应用。60蚁群优化算法—参考书2蚁群算法及其应用李士勇哈工大出版社国内首部蚁群算法的专著，系统地阐述蚁群算法的基本原理、基本蚁群算法及改进算法，蚁群算法与遗传、免疫算法的融合，自适应蚁群算法，并行蚁群算法，蚁群算法的收敛性与理论模型及其在优化问题中的应用。本书可供人工智能、计算机科学、信息科学、控制工程、管理工程、交通工程、网络工程、智能优化算法及智能自动化等领域的广大师生和科技人员学习及参考。61

蚁群优化算法—参考文献题目：群智能理论及应用电子学报,2003年S1期

【作者】彭喜元彭宇戴毓丰【关键词】群智能微粒群算法蚁群算法优化算法62基于群体智能的聚类算法CSI的研究CSI聚类算法主要步骤；基本模型简化：概率转换公式；实验结果

。63基于蚁群算法的聚类算法主要步骤：随机分布待聚类模式；每只蚂蚁计算当前对象在局部环境的群体相似度，并通过概率转换函数得到拾起或放下对象的概率，以这个概率行动；经过群体大量的相互作用，最终得到若干聚类中心；最后收集聚类结果。

64概率转换公式的简化基本模型简化模型65实验结果

66电信消费数据聚类分析实验结果比较

67基于群体智能的文档聚类算法CSIM的研究为了处理聚类过程中出现的散点以及克服算法的一些随机因素，更是为了提高算法的效率，我们将基于群体智能的文档聚类算法与经典的K均值算法相结合，对算法进行了改进。混合算法的过程是这样的：首先采用基于群体智能文档聚类算法对聚类文档进行处理，得到初始的聚类中心个数和聚类中心模板，然后运用K均值算法再次聚类。这样，既保留了群体智能算法的自组织特征，又结合了K均值算法的高效率，同时也克服了两种算法的弱点，如群体智能算法的随机性和K均值算法的聚类中心个数的参数预定及输入顺序敏感。我们将算法缩写为CSIM。68第3节改进的蚁群算法MacroDorigoGambardella69

带精英策略的蚂蚁系统带精英策略的蚂蚁系统（AntSystemwithelitiststrategy,ASelite）是最早的改进蚂蚁系统遗传算法中的精英策略传统的遗传算法可能会导致最适应个体的遗传信息丢失精英策略的思想是保留住一代中的最适应个体蚂蚁系统中的精英策略每次循环之后给予最优解以额外的信息素量这样的解被称为全局最优解（global-bestsolution）找出这个解的蚂蚁被称为精英蚂蚁(elitistants)70

带精英策略的蚂蚁系统信息素根据下式进行更新其中71

带精英策略的蚂蚁系统上式中表示精英蚂蚁引起的路径(i,j)上的信息素量的增加特点：可以使蚂蚁系统找出更优的解找到这些解的时间更短精英蚂蚁过多会导致搜索早熟收敛是精英蚂蚁的个数是所找出的最优解的路径长度72

蚁群系统蚁群系统(AntColonySystem,ACS)是由Dorigo和Gambardella在1996年提出的蚁群系统做了三个方面的改进：状态转移规则为更好更合理地利用新路径和利用关于问题的先验知识提供了方法全局更新规则只应用于最优的蚂蚁路径上在建立问题解决方案的过程中，应用局部信息素更新规则73蚁群系统状态转移规则一只位于节点r的蚂蚁通过应用下式给出的规则选择下一个将要移动到的城市s其中，S根据下列公式得到74蚁群系统状态转移规则q是在[0,1]区间均匀分布的随机数q0的大小决定了利用先验知识与探索新路径之间的相对重要性。上述状态转移规则被称为伪随机比例规则特点：倾向于选择短的且有着大量信息素的边作为移动方向75蚁群系统全局更新规则只有全局最优的蚂蚁才被允许释放信息素目的：使蚂蚁的搜索主要集中在当前循环为止所找出的最好路径的领域内全局更新在所有蚂蚁都完成它们的路径之后执行，使用下式对所建立的路径进行更新76蚁群系统全局更新规则为信息素挥发参数，0<<1为到目前为止找出的全局最优路径全局更新规则的另一个类型称为迭代最优区别:使用代替，为当前迭代(循环)中的最优路径长度这两种类型对蚁群系统性能的影响差别很小，全局最优的性能要稍微好一些77蚁群系统局部更新规则类似于蚁密和蚁量模型中的更新规则蚂蚁应用下列局部更新规则对它们所经过的边进行激素更新实验发现，可以产生好的结果，其中n是城市的数量，是由最近的邻域启发产生的一个路径长度局部更新规则可以有效地避免蚂蚁收敛到同一路径78最大-最小蚂蚁系统蚁群算法将蚂蚁的搜索行为集中到最优解的附近可以提高解的质量和收敛速度，从而改进算法的性能。但这种搜索方式会使早熟收敛行为更容易发生最大-最小蚂蚁系统(Max-MinAntSystem,MMAS)能将这种搜索方式和一种能够有效避免早熟收敛的机制结合在一起，从而使算法获得最优的性能79最大-最小蚂蚁系统MMAS和AS主要有三个方面不同：为了充分利用循环最优解和到目前为止找出的最优解，在每次循环之后，只有一只蚂蚁进行信息素更新。这只蚂蚁可能是找出当前循环中最优解的蚂蚁，也可能是找出从实验开始以来最优解的蚂蚁为避免搜索的停滞，在每个解的元素上的的信息素轨迹量的值域范围被限制在区间内将信息素轨迹初始化为80信息素轨迹更新在MMAS中，只有一只蚂蚁用于在每次循环后更新信息轨迹经修改的轨迹更新规则如下：表示迭代最优解或全局最优解的值在蚁群算法中主要使用全局最优解，而在MMAS中则主要使用迭代最优解81信息素轨迹的限制不管是选择迭代最优还是全局最优蚂蚁来进行信息素更新，都可能导致搜索的停滞。停滞现象发生的原因：在每个选择点上一个选择的信息素轨迹量明显高于其他的选择。避免停滞状态发生的方法：影响用来选择下一解元素的概率，它直接依赖于信息素轨迹和启发信息。通过限制信息素轨迹的影响，可以很容易地避免各信息素轨迹之间的差异过大。82信息素轨迹的限制MMAS对信息素轨迹的最小值和最大值分别施加了和的限制，从而使得对所有信息素轨迹，有MMAS收敛：在每个选择点上，其中一个解元素上的轨迹量为，而所有其他可选择的解元素上的轨迹量为。若MMAS收敛，通过始终选择信息素量最大的解元素所构造的解将与算法找出的最优解相一致83信息素轨迹的限制的选取的选取要基于两点假设最优解在搜索停滞发生之前不久被找出对解构造的主要影响是由信息素轨迹的上限与下限之间的相对差异决定84信息素轨迹的限制在一个选择点上选择相应解元素的概率Pdec直接取决于和在每个选择点上蚂蚁需在avg=n/2个解元素中选择蚂蚁构造最优解，需作n次正确的决策85信息素轨迹的初始化在第一次循环后所有信息素轨迹与相一致通过选择对这种类型的轨迹初始化来增加在算法的第一次循环期间对新解的探索当将信息素轨迹初始化为时，选择概率将增加得更加缓慢实验表明，将初始值设为可以改善最大-最小蚂蚁系统的性能86信息素轨迹的平滑化基本思想：通过增加选择有着低强度信息素轨迹量解元素的概率以提高探索新解的能力平滑机制有助于对搜索空间进行更有效的探索87蚁群算法的应用88混流装配线调度混流装配线(sequencingmixedmodelsonanassemblyline,SMMAL)是指一定时间内，在一条生产线上生产出多种不同型号的产品，产品的品种可以随顾客需求的变化而变化。SMMAL是车间作业调度问题(job-shopschedulingproblem,JSP)的具体应用之一。89问题描述以汽车组装为例，即在组装所有车辆的过程中，所确定的组装顺序应使各零部件的使用速率均匀化。如果不同型号的汽车消耗零部件的种类大致相同，那么原问题可简化为单级SMMAL调度问题。90问题描述i表示车型数的标号n表示需要装配的车型数m表示装配线上需要的零部件种类总数p表示生产调度中子装配的标号表示零部件p的理想使用速率j表示车型调度结果(即排序位置)的标号D表示在一个生产循环中需要组装的各种车型的总和91问题描述di表示在一个生产循环中车型i的数量bip表示生产每辆i车型需要零部件p的数量表示在组装线调度中前j-1台车消耗零部件p的数量和92蚁群算法在SMMAL中的应用假设有3种车型A、B、C排序，每个生产循环需A型车3辆，B型车2辆，C型车1辆，则每个循环共需生产6辆车。采用下图的搜索空间定义，列表示6个排序阶段，行表示有3种车型可以选择。蚁群算法就是不断改变圆圈的大小，最终寻找到满意的可行解。搜索的初始状态93简单SMMAL排序的搜索空间举例经过若干次迭代之后，搜索空间变化，此时最可能的可行解为B-A-C-A-B-A若干次迭代后的状态94局部搜索()的计算局部搜索采用的是贪心策略基本思路：每一步均从当前可选择策略中选取，使目标函数值增加最少的策略，即在确定第j个位置组装的车型时，如果有多种车型可供选择，则从中选择一种车型i，使第j个位置组装车型i时各零部件的使用速率最为均匀。95状态转移概率状态转移概率公式如下96信息素更新规则LB表示目标函数的下限值表示当前目标函数的平均值Zcutr表示当前的目标函数值这种动态标记的方法可在搜索过程中加大可行解间信息素的差别，避免算法早熟97实验数据98实验参数设置蚂蚁系统蚂蚁数量N_ant=5最大循环周期Ncmax=400=0.2Q=20000=0.9LB=0.0蚁群系统q0=0.5全局更新规则中的

和局部更新规则中的均取0.199实验参数设置最大-最小蚂蚁系统选取全局最优解带有精英策略的蚂蚁系统精英蚂蚁数量:1只100实验结果101实验结果分析直接用贪心策略求解结果：3293.4375蚂蚁系统求解SMMAL问题的性能较差对于这个具体的问题，带精英策略的蚂蚁系统的求解性能并不好于蚂蚁系统蚁群系统的性能相对于前两者而言，有了很大幅度的提高最大-最小蚂蚁系统的性能最好，大多数情况下的求解结果已达到实际的最优解102实验界面103实验界面104蚁群系统在TSP问题中的应用10城市TSP问题20城市TSP问题105蚁群系统在TSP问题中的应用30城市TSP问题48城市TSP问题106Questions?107第4节粒子群优化粒子群优化(ParticleSwarmOptimization,PSO)，又称微粒群算法，是由J.Kennedy和RCEberhart等于1995年开发的一种演化机制。

粒子(particle)”是一个折衷的选择,因为既需要将群体中的成员描述为没有质量、没有体积的,同时也需要描述它的速度和加速状态。ParticleSwarmOptimization(PSO)appliestoconceptofsocialinteractiontoproblemsolving.Itwasdevelopedin1995byJamesKennedyandRussEberhart[Kennedy,J.andEberhart,R.(1995).“ParticleSwarmOptimization”,Proceedingsofthe1995IEEEInternationalConferenceonNeuralNetworks,pp.1942-1948,IEEEPress.](/members/payman/swarm/kennedy95-ijcnn.pdf)

108粒子群优化

SwarmTopologyInPSO,therehavebeentwobasictopologiesusedintheliteratureRingTopology(neighborhoodof3)StarTopology(globalneighborhood)

I4

I0

I1

I2

I3

I4

I0

I1

I2

I3109特点分布式搜寻

具记忆性组件较少，容易实现

适合在连续性的范围内搜寻

110演算法介绍每个寻优的问题解都被想象成一只鸟，我们也称为“Particle”。所有的Particle都有一个fitnessfunction以判断目前的位置之好坏，

每一个Particle必须赋予记忆性，能记得所搜寻到最佳位置。每一个Particle还有一个速度以决定飞行的距离与方向。

111粒子群优化:

TheAnatomyofaParticleAparticle(individual)iscomposedof:Threevectors:Thex-vectorrecordsthecurrentposition(location)oftheparticleinthesearchspace,Thep-vectorrecordsthelocationofthebestsolutionfoundsofarbytheparticle,andThev-vectorcontainsagradient(direction)forwhichparticlewilltravelinifundisturbed.IkX=<xk0,xk1,…,xkn-1>P=<pk0,pk1,…,pkn-1>V=<vk0,vk1,…,vkn-1>x_fitness=?p_fitness=?112粒子群优化:

TheAnatomyofaParticleAparticle(individual)iscomposedof:Twofitnessvalues:Thex-fitnessrecordsthefitnessofthex-vector,andThep-fitnessrecordsthefitnessofthep-vector.IkX=<xk0,xk1,…,xkn-1>P=<pk0,pk1,…,pkn-1>V=<vk0,vk1,…,vkn-1>x_fitness=?p_fitness=?113速度更新Vid

：每一Particle在第d维之速度i：Particle