k均值聚类分析matlab代码

clear;
clc;
x=[0321546567];
y=[0821383445];
    ca1（1）=x（4）;%center of first cluster
    cb1（1）=y（4）;
    ca2（1）=x（2）;%center of second cluster
    cb2（1）=y（2）;
    ca3（1）=x（9）;%center of third cluster
    cb3（1）=y（9）;
    d1=0;d2=0;d3=0;Z=zeros（13）;d=0;
    suma1=ca1（1）;sumb1=cb1（1）;
    suma2=ca2（1）;sumb2=cb2（1）;
    suma3=ca3（1）;sumb3=cb3（1）;
    n1=1;n2=1;n3=1;
    k=1;
    t=0;
    cluster=zeros（110）;
    ra1=0.0;
    rb1=0.0;
    ra2=0.0;
    rb2=0.0;
    ra3=0.0;
    rb3=0.0;
while t==0
    for i=1:10
    d1=（x（i）-ca1（k））^2+（y（i）-cb1（k））^2;
    d2=（x（i）-ca2（k））^2+（y（i）-cb2（k））^2;
    d3=（x（i）-ca3（k））^2+（y（i）-cb3（k））^2;
    Z=[d1d2d3];
    d=min（Z）;
    if d==d1
        suma1=suma1+x（i）;
        sumb1=sumb1+y（i）;
        n1=n1+1;
        cluster（i）=1;
    elseif d==d2
        suma2=suma2+x（i）;