k均值聚类算方法MATLAB程序

function kmeans
clear all
clc
x=[0 1 0 1 2 1 2 3 6 7 8 6 7 8 9 7 8 9 8 9;
   0 0 1 1 1 2 2 2 6 6 6 7 7 7 7 8 8 8 9 9];
figure（1）
plot（x（1:）x（2:）‘r*‘）
%%第一步选取聚类中心，即令K=2
Z1=[x（11）;x（21）];
Z2=[x（110）;x（210）];
% K=1;%记录迭代的次数
dif1=inf;
dif2=inf;
%%第二步计算各点与聚类中心的距离
while（dif1>eps&dif2>eps）
   R1=[];
   R2=[]; 
for i=1:20
    dist1=sqrt（（x（1i）-Z1（1））.^2+（x（2i）-Z1（2））.^2）;
    dist2=sqrt（（x（1i）-Z2（1））.^2+（x（2i）-Z2（2））.^2）;
    temp=[x（1i）x（2i）]‘;
    if dist1