搜索
clearclcx=[00;10;01;11;21;12;32;66;76;86;67;77;87;97;78;88;98;89;99];z=zeros(2,2);z1=zeros(2,2);z=x(1:2,1:2);%%寻找聚类中心while1count=zeros(2,1);allsum=zeros(2,2);fori=1:20%对每一个样本,计算到两个聚类中心的距离temp1=sqrt((z(1,1)-x(i,1)).^2+(z(1,2)-x(i,2)).^2);temp2=sqrt((z(2,1)-x(i,1)).^2+(z(2,2)-x(i,2)).^2);if(temp1temp2)count(1)=count(1)+1;allsum(1,1)=allsum(1,1)+x(i,1);allsum(1,2)=allsum(1,2)+x(i,2);elsecount(2)=count(2)+1;allsum(2,1)=allsum(2,1)+x(i,1);allsum(2,2)=allsum(2,2)+x(i,2);endendz1(1,1)=allsum(1,1)/count(1);z1(1,2)=allsum(1,2)/count(1);z1(2,1)=allsum(2,1)/count(2);z1(2,2)=allsum(2,2)/count(2);if(z==z1)breakelsez=z1;endend%%结果显示disp(z1)%%输出聚类中心plot(x(:,1),x(:,2),'b*');holdonplot(z1(:,1),z1(:,2),'ro');title('K均值分类图')xlabel('特征x1')ylabel('特征x2')