Kmean.rar_kmean_kmean 二维_matlab 随机点_matlab随机点_二维 高斯分布

function [cluster,u,sum_old]=kmfunc2(r,k) load('Data.mat'); x=dataset2; sum_old=inf; L=size(x); l=L(1); zc=zeros(r,10); for yy=1:r u=rand(k,2); for zz=1:10 cluster=zeros(l,1); count=zeros(k,1); u_new=zeros(k,2); for ii=1:l maxi=inf; for jj=1:k; dist=norm(x(ii,:)-u(jj,:)); if dist<maxi maxi=dist; cluster(ii)=jj; end end count(cluster(ii))=count(cluster(ii))+1; u_new(cluster(ii),:)=u_new(cluster(ii),:)+x(ii,:); end for ii=1:k; if count(ii)~=0 u_new(ii,:)=u_new(ii,:)/count(ii); end end u=u_new; sum=0; for ii=1:l sum=sum+norm(x(ii,:)-u(cluster(ii),:))^2; end zc(yy,zz)=sum; end if sum_old>sum sum_old=sum; mii=yy; end end figure (2); subplot(1,2,1); for ii=1:l if cluster(ii)==1 plot (x(ii,1),x(ii,2),'.r'); hold on end if cluster(ii)==2 plot (x(ii,1),x(ii,2),'.g'); hold on end if cluster(ii)==3 plot (x(ii,1),x(ii,2),'.b'); hold on end end hold off; subplot(1,2,2) plot(1:10,zc(mii,:)/l); xlabel('iteration times') ylabel('sum-squared-error (divided by the sample size n)') end