%%?K-mens方法的matlab實現
%%?數據準備和初始化
clc
clear?
x=[164?62;156?50;168?86;160?60;162?54;187?75;162?55;167?75;160.5?56;160?53;158?55;164?60;165?50;174?64;166?55;158?47.5;162?60;175?62;170?65;161?49;169?55;161?46;160?45;167?44;176?73;169?58;178?54;165?66;155?49;183?68;171?61;179?64;154?57;172?60;172?52;173?59;172?58;175?62];
z=zeros（22）;
z1=zeros（22）;
z=x（1:21:2）;
%%?尋找聚類中心
while?1
????count=zeros（21）;
????allsum=zeros（22）;
????for?i=1:37%對每一個樣本i，計算到2個聚類中心的距離
???????temp1=sqrt（（z（11）-x（i1））.^2+（z（12）-x（i2））.^2）;
???????temp2=sqrt（（z（21）-x（i1））.^2+（z（22）-x（i2））.^2）;
????????if（temp1????????????count（1）=count（1）+1;
????????????allsum（11）=allsum（11）+x（i1）;
????????????allsum（12）=allsum（12）+x（i2）;
????????else
????????????count（2）=count（2）