数据聚类实验报告记录附代码文档格式.docx

1、本实验采用的是k-means聚类算法，类中心的选择基于簇中对象的平均值。（1）获取用户的类中心数目k和具有n个对象的数据集；（2）任意选取k个对象作为初始的簇中心；（3）根据簇中对象的平均值，将每个对象重新赋给最类似的簇；（4）更新簇的平均值，即类中心，重新计算每个簇中对象的平均值；（5）如果新的平均值发生改变，转至（2）步；（6）新的平均值不在发生改变，算法聚类结束。2.2 算法流程图 开 始 选择k值，选取k个 对象作为平均值 聚 类 对新生成的簇重新 计算平均值 Y 新的平均 值改变 N 结 束 图1 k-means算法流程图 k-means算法流程图，如图1所示。 k-means算法中

2、的k，由用户输入，最终得到的类别数即为用户输入的数目。聚类过程中，涉及到初始类中心的选择。在程序中，对于类中心，是选择前k个作为初始类中心，对于数据的组织，前k个数据有较大差别，可以提高程序的运行效率和分类结果的准确率。3 实验结果分析 在实验中，利用k-means聚类算法对“ch7 iris.txt”数据集进行聚类实验。当k=3时，k-means算法聚类效果如图2所示： 图2 k=3时聚类效果 当k=4时，k-means算法聚类效果如图3所示： 图3 k=4时聚类效果 K-means聚类算法的收敛性和初值的选取有关。初始的聚类中心的不同，对聚类结果没有很大的影响，而对迭代次数有显著的影响。数

3、据的输入顺序不同，同样影响迭代次数，而对聚类结果没有太大的影响。4 实验结论 K-means聚类算法对于类别数的选择k值有较高的要求，如果类别数较少，则不能区分数据。K-means聚类算法找出平均误差最小的k个划分。当结果簇是密集的，而簇与簇之间的区别明显时，它的效果较好。该算法只有在簇的平均值被定义的情况下才能使用。对于初始类中心的选择，特别重要。对于分类的准确度和距离影响明显。而且该算法对孤立点是敏感的。所以如果数据集中存在有极大值的对象，应该消除这种敏感性。5 实验心得体会1、初始值可的选取 K-means聚类算法对于类别数目的选择，需要使用该算法的人员对于数据分类有一定的了解，并且可以

4、根据观察部分原始抽样数据，得出该样本数据的大致类别数目，否则，应用该方法的聚类可能会出现较大的错误率。2、初始类中心的选取初始类中心的选择对聚类的准确度有较大的影响。在初始类中心的选择时，最好选择两两距离较大，且能代表不同数据样本类别的点作为初始的类中心点。参考文献1 数据挖掘：概念与技术/（加）韩家炜，（加）坎伯（Kamber,M.）著；范明等译.-北京：机械工业出版社，2001.8 .2 效琴，戴汝源数据挖掘中聚类分析的技术方法J微计算机信息，2003，19（1）3 贺玲，吴玲达，蔡益朝数据挖掘中的聚类算法综述J计算机应用研究，2007，24（1）:10-134 孙吉贵，刘杰，赵连宇聚类算

5、法研究J软件学报，2008，19（1）：48-615 冯晓蒲，张铁峰四种聚类方法之比较J微型机与应用，2010，16附录（源代码）Matlab%K-means算法主程序 k=3;x= 0.224 0.624 0.067 0.043 0.749 0.502 0.627 0.541 0.557 0.541 0.847 1.000 0.110 0.502 0.051 0.043 0.722 0.459 0.663 0.584 0.776 0.416 0.831 0.831 0.196 0.667 0.067 0.043 0.612 0.333 0.612 0.584 0.612 0.416 0.81

6、2 0.875 0.055 0.584 0.067 0.082 0.557 0.541 0.627 0.624 0.165 0.208 0.592 0.667 0.027 0.376 0.067 0.043 0.639 0.376 0.612 0.498 0.667 0.208 0.812 0.710 0.306 0.710 0.086 0.043 0.196 0.000 0.424 0.376 0.612 0.502 0.694 0.792 0.137 0.416 0.067 0.000 0.471 0.082 0.510 0.376 0.694 0.416 0.761 0.831 0.41

7、6 0.831 0.035 0.043 0.361 0.376 0.439 0.498 0.416 0.333 0.694 0.957 0.306 0.792 0.051 0.125 0.361 0.416 0.592 0.584 0.612 0.416 0.761 0.710 0.388 0.749 0.118 0.082 0.529 0.082 0.592 0.584 0.945 0.251 1.000 0.918 0.306 0.584 0.118 0.043 0.443 0.502 0.643 0.459 0.722 0.502 0.796 0.918 0.082 0.667 0.00

8、0 0.043 0.557 0.208 0.663 0.584 0.945 0.333 0.965 0.792 0.137 0.584 0.153 0.043 0.584 0.376 0.561 0.498 0.667 0.541 0.796 0.831 0.196 0.584 0.102 0.125 0.694 0.333 0.643 0.541 0.529 0.333 0.643 0.710 0.251 0.584 0.067 0.043 0.471 0.376 0.592 0.584 0.584 0.333 0.780 0.831 0.137 0.459 0.102 0.043 0.33

9、3 0.169 0.475 0.416 0.863 0.333 0.863 0.749 0.251 0.875 0.086 0.000 0.416 0.290 0.490 0.459 0.584 0.333 0.780 0.875 0.165 0.459 0.086 0.043 0.306 0.416 0.592 0.584 0.498 0.251 0.780 0.541 0.333 0.624 0.051 0.043 0.667 0.459 0.627 0.584 0.557 0.584 0.780 0.957 0.027 0.416 0.051 0.043 0.361 0.416 0.52

10、5 0.498 0.471 0.416 0.643 0.710 0.196 0.624 0.051 0.082 0.333 0.251 0.576 0.459 0.667 0.459 0.780 0.957 0.027 0.502 0.051 0.043 0.416 0.251 0.510 0.459 0.416 0.290 0.694 0.749 0.224 0.749 0.153 0.125 0.361 0.290 0.541 0.498 0.667 0.541 0.796 1.000 0.224 0.749 0.102 0.043 0.388 0.376 0.541 0.498 0.55

11、7 0.208 0.678 0.749 0.278 0.710 0.086 0.043 0.224 0.208 0.337 0.416 0.529 0.584 0.745 0.918 0.165 0.416 0.067 0.043 0.584 0.502 0.592 0.584 0.082 0.459 0.086 0.043 0.333 0.125 0.510 0.498 0.557 0.376 0.780 0.710 0.306 0.792 0.118 0.125 0.388 0.333 0.592 0.498 0.918 0.416 0.949 0.831 0.196 0.584 0.08

12、6 0.043 0.165 0.169 0.388 0.376 0.835 0.376 0.898 0.710 0.165 0.459 0.086 0.000 0.251 0.290 0.490 0.541 0.804 0.667 0.863 1.000 0.137 0.584 0.102 0.043 0.443 0.416 0.541 0.584 0.584 0.290 0.729 0.749 0.000 0.416 0.016 0.000 0.498 0.376 0.627 0.541 0.388 0.208 0.678 0.792 0.388 1.000 0.086 0.125 0.66

13、7 0.459 0.576 0.541 0.584 0.502 0.729 0.918 0.224 0.624 0.067 0.082 0.416 0.290 0.525 0.376 0.945 0.749 0.965 0.875 0.224 0.749 0.086 0.082 0.361 0.208 0.490 0.416 0.471 0.082 0.678 0.584 0.224 0.710 0.086 0.125 0.498 0.333 0.510 0.498 0.361 0.333 0.663 0.792 0.224 0.541 0.118 0.165 0.498 0.333 0.62

14、7 0.459 0.557 0.290 0.663 0.710 0.196 0.416 0.102 0.043 0.639 0.416 0.576 0.541 0.804 0.502 0.847 0.710 0.251 0.624 0.086 0.043 0.667 0.416 0.678 0.667 0.498 0.416 0.510 0.710 0.110 0.502 0.102 0.043 0.388 0.251 0.424 0.376 0.804 0.416 0.812 0.624 0.306 0.584 0.086 0.125 0.333 0.169 0.459 0.376 1.00

15、0 0.749 0.914 0.792 0.333 0.918 0.067 0.043 0.471 0.290 0.694 0.624 0.557 0.333 0.694 0.584 0.196 0.502 0.035 0.043 0.471 0.584 0.592 0.624 0.945 0.416 0.863 0.918 0.165 0.667 0.067 0.000 0.557 0.125 0.576 0.498 0.584 0.459 0.761 0.710 0.224 0.584 0.086 0.043 0.333 0.208 0.510 0.498 0.722 0.459 0.74

16、5 0.831 0.055 0.125 0.051 0.082 0.498 0.416 0.612 0.541 0.722 0.459 0.694 0.918 0.196 0.624 0.102 0.208 0.196 0.125 0.388 0.376 0.694 0.502 0.831 0.918 0.137 0.416 0.067 0.082 0.388 0.416 0.541 0.459 0.667 0.416 0.714 0.918 0.082 0.502 0.067 0.043 0.529 0.376 0.561 0.498 0.612 0.416 0.714 0.792 0.19

17、6 0.541 0.067 0.043 0.388 0.333 0.525 0.498 0.443 0.416 0.694 0.710 ;n,d = size（x）;bn=round（n/k*rand）;%第一个随机数在前1/K的范围内 nc=x（bn,:）;x（2*bn,:x（3*bn,:;%初始聚类中心 %nc=x（bn,:x（4*bn,:% 4类cid,nr,centers = kmeans（x,k,nc）%调用kmeans函数 for i=1:150 if cid（i）=1, plot（x（i,1）,x（i,2）,r*） % 显示第一类 hold on else if cid（i）=2

18、, b*） %显示第二类 if cid（i）=3, g*） %显示第三类 %else %if cid（i）=4, %plot（x（i,1）,x（i,2）,k*） %显示第四类 %hold on %end end strt=红色*为第一类；蓝色*为第二类；绿色*为第三类；黑色*为第四类 ;text（-4,-3.6,strt）; %kmeans.m主类function cid,nr,centers = kmeans（x,k,nc） % 设置cid为分类结果显示矩阵cid = zeros（1,n）;oldcid = ones（1,n）;nr = zeros（1,k）;maxgn= 100;iter

19、= 1;while iter maxgn %计算每个数据到聚类中心的距离 for i = 1:n dist = sum（repmat（x（i,:）,k,1）-nc）.2,2）;m,ind = min（dist）; % 将当前聚类结果存入cid中 cid（i） = ind;k %找到每一类的所有数据，计算他们的平均值，作为下次计算的聚类中心ind = find（cid=i）;nc（i,:） = mean（x（ind,:）;% 统计每一类的数据个数 nr（i） = length（ind）;iter = iter + 1;maxiter = 2;move = 1; maxiter & move =

20、0 move = 0; % 对所有的数据进行再次判断，寻求最佳聚类结果 for i = 1: dist = sum（repmat（x（i,: r = cid（i）; % 将当前数据属于的类给r dadj = nr./（nr+1）.*dist; % 计算调整后的距离 m,ind = min（dadj）; % 早到该数据距哪个聚类中心最近 if ind = r % 如果不等则聚类中心移动 cid（i） = ind;%将新的聚类结果送给cid ic = find（cid = ind）;%重新计算调整当前类别的聚类中心 nc（ind,:） = mean（x（ic,: move = 1; end iter = iter+1;endcenters = nc;if move = 0 disp（No points were moved after the initial clustering procedure.） elseSome points were moved after the initial clustering procedure.