clear?all
clc
close?all
mdist=input（‘輸入坐標文件名字\n‘）;
disp（‘讀取數據坐標‘）
%獲取坐標
%文件格式為數據的二維的坐標，
%共兩列第一列為x軸坐標，第二列為y軸坐標
xx=load（mdist）;
%獲取數據的個數存至number
[number?row]=size（xx）;
%獲取距離矩陣，第二參數指定距離計算方法
%可選參數有‘euclidean‘：歐氏距離（默認）；‘seuclidean‘：標準化歐氏距離；
%‘mahalanobis‘：馬氏距離；‘cityblock‘：布洛克距離；‘minkowski‘：明可夫斯基距離；
%‘cosine‘：余弦距離?‘correlation‘：相關性?‘hamming‘：漢明距離?‘jaccard‘：Jaccard相似度
%‘chebychev‘：Chebychev距離。
yy=pdist（xx‘euclidean‘）;
%獲取聚類，第二參數指定層次聚類方式
%‘single‘：單連通，最短距離法（默認）；‘complete‘：全連通，最長距離法；‘average‘：未加權平均距離法；?
%‘weighted‘：?加權平均法；‘centroid‘：?質心距離法；‘median‘：加權質心距離法；‘ward‘：內平方距離法（最小方差算法）
zz=linkage（yy‘single‘）;
%指定獲取簇類個數
Ncluster=input（‘輸入類個數\n‘）;
%獲取指定Ncluster個數的聚類結果
c?=?cluster（?zz‘maxclust‘?Ncluster?）;?
?
%獲取繪圖所需顏色
%需要用到linspecer.m文件
%下載地址
%http://www.mathworks.com/matlabcentral/fileexchange/42673-beautiful-and-distinguishable-line-colors-+-colormap
?if（Ncluster>12）
?Color?=?linspecer（?Ncluster?）;
?else
?????Color?=?linspecer（?Ncluster?‘qualitative‘?）;
?????end
?????
?
?for?i=1:Ncluster
?????for?j?=?1:number
?????if（c（j）?==?i）????
?????hold?on
?????plot（xx（j1）xx（j2）‘o‘‘MarkerFaceColor‘Color（i:）‘MarkerEdgeColor‘Color（i:））
?????end
?????end

?end