%% Script file for examples of Kohonen's Map
load Data6;
X=XX{1};
[mm,nn]=size(X);
clear XX group_num;

NumClusters=10;
C=rand(NumClusters,2);
Corig=C;

erate=0.01;
flag=0;

%With no topology;
for j=1:5*mm
    k=ceil(rand*mm);
    for jj=1:7
      dd(jj)=norm(X(k,:)-C(jj,:));
    end  
      [val,win]=min(dd);
      C(win,:)=C(win,:)+erate*(X(k,:)-C(win,:));
    
    if mod(j,10)==0
        flag=flag+1;
        voronoi(C(:,1),C(:,2));
        hold on
        plot(C(:,1),C(:,2),'r*',X(:,1),X(:,2),'k.');
        axis([0 1 0 1])
        pause(0.01)
        hold off
    end
end