%%?input?the?data?隨機200個點
A=rand（10002）;?
x=A（:1）;?
y=A（:2）;
plot（xy‘r.‘）;?
??
%%?計算質(zhì)心??
centroids?=?kMeans（A?3）;??

%%?取得隨機中心??
function?[?centroids?]?=?randCent（?dataSet?k?）??
????[mn]?=?size（dataSet）;%取得列數(shù)??
????centroids?=?zeros（k?n）;??
????for?j?=?1:n??
????????minJ?=?min（dataSet（:j））;??
????????rangeJ?=?max（dataSet（:j））-min（dataSet（:j））;??
????????centroids（:j）?=?minJ+rand（k1）*rangeJ;%產(chǎn)生區(qū)間上的隨機數(shù)??
????end??
end?

function?[?dist?]?=?distence（?vecA?vecB?）??
????dist?=?（vecA-vecB）*（vecA-vecB）‘;%這里取歐式距離的平方??
end

%%?kMeans的核心程序，不斷迭代求解聚類中心??
function?[?centroids?]?=?kMeans（?dataSet?k?）??
????[mn]?=?size（dataSet）;??
????%初始化聚類中心??
????centroids?=?randCent（dataSet?k）;??
????subCenter?=?zeros（m2）;%做一個m*2的矩陣，第一列存儲類別，第二列存儲距離??
????change?=?1;%判斷是否改變??
????while?change?==?1??
????????change?=?0;??
????????%對每一組數(shù)據(jù)計算距離??
????????for?i?=?1:m??
????????????minDist?=?inf;??
????????????minIndex?=?0;??
????????????for?j?=?1:k??
?????????????????dist=?distence（dataSet（i:）?centroids（j:））;??
?????????????????if?dist??????????????????????minDist?=?dist;??