基本蚁群算法及其改进算法.ppt

《基本蚁群算法及其改进算法.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、蚁群算法的基本原理及其改进算法专业：控制工程年级：2009级姓名：胡训智学号：30956060指导老师：周润景教授算法的提出算法的基本原理模型建立算法的实现算法改进结论参考文献蚁群算法的提出蚁群算法(antcolonyoptimization,ACO)，又称蚂蚁算法，是一种用来寻找优化路径的机率型算法。它由MarcoDorigo于1992年在他的博士论文中提出，其灵感来源于蚂蚁在寻找食物过程中发现路径的行为。MacroDorigo基本原理NestFoodObstacle图1蚂蚁正常行进，突然环境改变，增加了障碍

2、物基本原理NestFoodObstacle图2蚂蚁以等同概率选择各条路径较短路径信息素浓度高，选择该路径的蚂蚁增多基本原理图3蚂蚁选路过程示例EABDHCEABDHCd=0.5d=0.5d=1d=130ants30ants15ants15ants15ants15antst=0EABDHC30ants30ants20ants20ants10ants10antst=1基本原理NestFoodObstacle图4蚂蚁最终绕过障碍物找到最优路径模型建立基于蚂蚁构造墓地和分类幼体的聚类分析模型基于蚂蚁觅食行为和信息素的聚

3、类分析模型基于蚂蚁构造墓地和分类幼体的聚类分析模型蚁群构造墓地行为和分类幼体行为统称之为蚁群聚类行为。生物学家经过长期的观察发现，在蚂蚁群体中存在一种本能的聚集行为。蚂蚁往往能在没有关于蚂蚁整体的任何指导性信息情况下，将其死去的同伴的尸体安放在一个固定的场所。真实蚁群的聚类行为DeneubougJL等人也用pheidolepallidula蚂蚁做了实验。发现蚁群会根据蚂蚁幼体的大小将其放置在不同的位置，分别把其堆放在蚁穴周围和中央的位置。真实的蚁群聚类行为的实验结果右图，四张照片分别对应为实验初始状态、3小时、

4、6小时和36小时的蚁群聚类情况。基于蚂蚁构造墓地和分类幼体的聚类分析模型基本模型经过利用个体与个体和个体与环境之间的交互作用，实现了自组织聚类，并成功的应用于机器人的控制中(一群类似于蚂蚁的机器人在二维网格中随意移动并可以搬运基本物体，最终把它们聚集在一起)。该模型成功的应用引起了各国学者的广泛关注和研究的热潮。LumerE和FaietaB通过在Denurbourg的基本分类模型中引入数据对象之间相似度的概念，提出了LF聚类分析算法，并成功的将其应用到数据分析中。基于蚂蚁觅食行为和信息素的聚类分析模型蚂蚁在觅食

5、的过程中，能够分为搜索食物和搬运食物两个环节。每个蚂蚁在运动过程中都将会在其所经过的路径上留下信息素，而且能够感知到信息素的存在及其强度，比较倾向于向信息素强度高的方向移动，同样信息素自身也会随着时间的流逝而挥发，显然某一路径上经过的蚂蚁数目越多，那么其信息素就越强，以后的蚂蚁选择该路径的可能性就比较高，整个蚁群的行为表现出了信息正反馈现象。程序流程图程序初始化X=load('data.txt');[N,n]=size(X);%N=测试样本数;n=测试样本的属性数;K=4;%K=组数;R=100;%R=蚂蚁数;

6、t_max=1000;%t_max=最大迭代次数;best_solution_function_value=inf;信息素矩阵初始化信息素矩阵维数为N*K（样本数*聚类数）初始值为0.01。c=10^-2;tau=ones(N,K)*c;%信息素矩阵，初始值为0.01的N*K矩阵（样本数*聚类数）蚂蚁路径的选择及标识定义标识字符矩阵solution_string，维数为R*N+1，初始值都为0，以信息矩阵中信息素的值确定路径（即确定分到哪一组），具体方法如下：如果该样本各信息素的值都小于信息素阈值q，则取信息素

7、最大的为作为路径。若最大值有多个，则从相同的最大值中随机取一个，作为路径。若信息数大于阈值q，则求出各路径信息素占该样本总信息素的比例，以概率确定路径。聚类中心选择聚类中心为该类所有样本的各属性值的平均值。偏离误差计算偏离误差的计算，即各样本到其对应的聚类中心的欧式距离之和MIN。MIN越小，聚类效果越好。计算各只蚂蚁的MIN值，找到最小的MIN值，该值对应的路径为本次迭代的最佳路径。信息素更新对信息素矩阵进行更新，更新方法为新值为原信息素值乘以(1-rho)，rho为信息素蒸发率，在加上最小偏差值的倒数。fo

8、ri=1:Ntau(i,best_solution(1,i))=(1-rho)*tau(i,best_solution(1,i))+1/tau_F;信息数更新之后，再根据新的信息数矩阵，判断路径。进行迭代运算。直到达到最大迭代次数，或偏离误差达到要求值。参数调试（蚂蚁数R最大迭代次数t_max）Rt_maxMINTime（s）10100504312.1007101000450457.2