%------------------------------------------------------------
% K- Means Clustering Algorithm
%-----------------------------------------------------------
clear all
%生成數(shù)據(jù)集，共生成三個(gè)類
%--------------------------------------------------------
mu1=[0 10];
sigma1=[1 0;0 1];
cluster1 = mvnrnd(mu1,sigma1,200);
mu2=[10 0];
sigma2=[1 .5;.5 2];
cluster2 = mvnrnd(mu2,sigma2,200);
mu3=[0 0];
sigma3=[3 1;1 2];
cluster3 = mvnrnd(mu3,sigma3,200);
data=[cluster1;cluster2;cluster3];
plot(data(:,1),data(:,2),'r.')
hold on;
%---------------------------------------------------------
%數(shù)據(jù)的初始化
%------------------------------------------------------
k=3;
dist=zeros(k,1);
[dataNumber,cols]=size(data);
indexdata=zeros(dataNumber,2);
for i=1:dataNumber
    indexdata(i,1)=i;
end

%生成質(zhì)心集和質(zhì)心集副本并初始化
center=zeros(k,cols);
centercopy=zeros(k,cols);
flag=50;

for i=1:k
    centercopy(i,1)=data(i,1);
    centercopy(i,2)=data(i,2);
end
plot(centercopy(:,1),centercopy(:,2),'go','MarkerSize',10);



%計(jì)算每個(gè)數(shù)據(jù)到各質(zhì)心的距離,劃分聚類
while(flag~=0)
    clusterNumber=ones(k,1);
    for i=1:dataNumber
                winner=1;
        dist(winner)=inf; 
            for j=1:k
                dist(j)=(centercopy(j,1)-data(i,1))*(centercopy(j,1)-data(i,1))+(centercopy(j,2)-data(i,2))*(centercopy(j,2)-data(i,2));      
            end
            
            for j=2:k
                if(dist(j)<dist(winner))
                    winner=j;
                end
            end
            
        indexdata(i,2)=winner;
        center(winner,1)=center(winner,1)+data(i,1);
        center(winner,2)=center(winner,2)+data(i,2);
        clusterNumber(winner)=clusterNumber(winner)+1;
    end

    %計(jì)算更新的簇均值
    for j=1:k
        center(j,1)=center(j,1)/clusterNumber(j);
        center(j,2)=center(j,2)/clusterNumber(j);
    end
    
    flag=flag-1;
    centercopy=center;
%     if(sum(sum(center==centercopy))==1)
%         flag=1;
%     else
%         centercopy=center;
%     end
end

%draw center and cluster
for i=1:dataNumber
    switch  indexdata(i,2)
        case 1
            plot(data(i,1),data(i,2),'r.');
            hold on;
        case 2
            plot(data(i,1),data(i,2),'g.');
            hold on;
        case 3
            plot(data(i,1),data(i,2),'b.');
            hold on;
    end
end
plot(center(:,1),center(:,2),'ko','MarkerSize',10);