基于改进人工蜂群算法的 K 均值聚类算法(Matlab代码实现)

💥💥💞💞欢迎来到本博客❤️❤️💥💥

****🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。

⛳️**座右铭:**行百里者,半于九十。

📋📋📋++本文目录如下:++🎁🎁🎁

目录

[💥1 概述](#💥1 概述)

[📚2 运行结果](#📚2 运行结果)

[🎉3 参考文献](#🎉3 参考文献)

[🌈4 Matlab代码、数据、文章讲解](#🌈4 Matlab代码、数据、文章讲解)


💥1 概述

基于以上三点改进提出了 IABC-Kmeans 算法。该算法的基本思想是: 通过 IABC 算法进行一次迭代,将迭代得到的新位置作为 KMC 的初始点并进行一次 K 均值聚类,再用聚类获得的新的中心点更新蜂群; 如此交替执行 IABC 算法和 K均值聚类,直到算法结束。

算法基本步骤描述如下:

步骤 1 设置引领蜂、跟随蜂和侦察蜂的数量( 一般有引领蜂数量 = 跟随蜂数量) ; 最大迭代次数 MCN 以及控制参数Limit; 当前迭代次数Cycle,初始值为1; 聚类类别数 k; 利用最大最小距离积法初始化蜂群,产生 { Z1,Z2,...,ZN } 个初始蜂 群。

步骤 2 对初始蜂群进行一次聚类划分,根据式( 5) 计算每只蜜蜂的适应度,按照适应度大小排序,将前一半作为引领蜂,后一半作为跟随蜂。

步骤 3 引领蜂利用式( 6) 对其邻域进行搜索,得到新的位置,按照贪婪选择原则,如果新的位置的适应度大于原先位置的适应度,则用新的位置更新原位置; 否则,保持原位置不变。当所有引领蜂完成邻域搜索后,根据式( 4) 计算概率 Pi。

步骤 4 跟随蜂利用算得的概率 Pi 并基于轮盘赌原则选择引领蜂,原则上,Pi 越大,表明引领蜂 i 的适应度值越大,被跟随蜂选中的概率也越大。当跟随蜂完成引领蜂选择后,利用式( 7) 对邻域搜索,同样按照贪婪选择原则选择适应度高的位置。

步骤 5 在所有跟随蜂完成搜索后,将得到的位置作为聚类中心,对数据集进行一次 K 均值迭代聚类,根据聚类划分,用每一类的新的聚类中心更新蜂群。

步骤 6 如果某引领蜂在 Limit 次迭代后,结果都没有改变,则由引领蜂变为侦察蜂,并随机产生一个新的位置取代原位置。

步骤 7 如果当前迭代次数大于最大次数 MCN,则迭代结束,算法结束; 否则转向步骤 2,Cycle = Cycle + 1。

📚 2 运行结果

复现结果:

部分代码:

function [Colony Obj Fit oBas]=GreedySelection(Colony1,Colony2,ObjEmp,ObjEmp2,FitEmp,FitEmp2,fbas,ABCOpts,i)

oBas=fbas;

Obj=ObjEmp;

Fit=FitEmp;

Colony=Colony1;

if (nargin==8)%Inside the body of a user-defined function, NARGIN returns the number of input arguments that were used to call the function.

for ind=1:size(Colony1,1)%ind=1:5,对所有食物源进行贪婪选择

if (FitEmp2(ind)>FitEmp(ind))%如果Vi的适应度值大于Xi的,替换,,

oBas(ind)=0;

%zj因为这是已经被新的位置更新了,所以其开采度应该置为零,表示这是第一次,没有被开采过

Fit(ind)=FitEmp2(ind);

Obj(ind)=ObjEmp2(ind);

Colony(ind,:)=Colony2(ind,:);

else%否则不变,并且计数器bas+1

oBas(ind)=fbas(ind)+1;

%zj因为新的位置的适应度没有当前的好(大),所以在当前位置上仍保留当前解,表示当前又被开采了一次

Fit(ind)=FitEmp(ind);

Obj(ind)=ObjEmp(ind);

Colony(ind,:)=Colony1(ind,:);

end;

end; %for

end; %if

if(nargin==9)%第i个引领蜂被跟随,只对第i个食物源进行贪婪选择

ind=i;

if (FitEmp2(ind)>FitEmp(ind))

oBas(ind)=0;

Fit(ind)=FitEmp2(ind);

Obj(ind)=ObjEmp2(ind);

Colony(ind,:)=Colony2(ind,:);

else

oBas(ind)=fbas(ind)+1;

Fit(ind)=FitEmp(ind);

Obj(ind)=ObjEmp(ind);

Colony(ind,:)=Colony1(ind,:);

end;

end;

🎉3 参考文献

文章中一些内容引自网络,会注明出处或引用为参考文献,难免有未尽之处,如有不妥,请随时联系删除。

[1]喻金平,郑杰,梅宏标.基于改进人工蜂群算法的K均值聚类算法[J].计算机应用,2014,34(04):1065-1069+1088.

🌈4 Matlab代码、数据、文章讲解

相关推荐
APP 肖提莫1 分钟前
MyBatis-Plus分页拦截器,源码的重构(重构total总数的计算逻辑)
java·前端·算法
OTWOL9 分钟前
两道数组有关的OJ练习题
c语言·开发语言·数据结构·c++·算法
地球资源数据云21 分钟前
全国30米分辨率逐年植被覆盖度(FVC)数据集
大数据·运维·服务器·数据库·均值算法
qq_4335545436 分钟前
C++ 面向对象编程:递增重载
开发语言·c++·算法
带多刺的玫瑰1 小时前
Leecode刷题C语言之切蛋糕的最小总开销①
java·数据结构·算法
巫师不要去魔法部乱说1 小时前
PyCharm专项训练5 最短路径算法
python·算法·pycharm
dundunmm2 小时前
论文阅读:Deep Fusion Clustering Network With Reliable Structure Preservation
论文阅读·人工智能·数据挖掘·聚类·深度聚类·图聚类
qystca2 小时前
洛谷 P11242 碧树 C语言
数据结构·算法
汤姆和佩琦2 小时前
2024-12-25-sklearn学习(20)无监督学习-双聚类 料峭春风吹酒醒,微冷,山头斜照却相迎。
学习·聚类·sklearn
冠位观测者2 小时前
【Leetcode 热题 100】124. 二叉树中的最大路径和
数据结构·算法·leetcode