x1=[0 0];x2=[0 1];x3=[1 0];x4=[1 1];x5=[2 1];
x6=[1 2];x7=[2 2];x8=[3 2];x9=[6 6];x10=[7 6];
x11=[8 6];x12=[6 7];x13=[7 7];x14=[8 7];x15=[9 7];
x16=[7 8];x17=[8 8];x18=[9 8];x19=[8 9];x20=[9 9];
X=[x1;x2;x3;x4;x5;x6;x7;x8;x9;x10;
    x11;x12;x13;x14;x15;x16;x17;x18;x19;x20];
c1=X(1,:);z1=X(1,:);c2=X(2,:);z2=X(2,:);s1=[];s2=[];  %初始化
for j=1:10
    for i=1:20
        if norm(X(i,:)-z1)<norm(X(i,:)-z2) %判別條件  
            s1=[s1;X(i,:)];
            z1=mean(s1);
        else 
            s2=[s2;X(i,:)];
            z2=mean(s2);
        end
    end
    if c1==z1&c2==z2
        break        %當滿足分類條件后跳出循環，分類結束
    else
        c1=z1;c2=z2;
    end
    s1=[];s2=[];   %若分類未結束，則清空分類情況 
end
plot(s1(:,1),s1(:,2),'ro',s2(:,1),s2(:,2),'bo');
xlabel('x1');ylabel('x2');grid;