MATLAB----C均值聚类分析

clear all;
data=[
1702.8	1639.79	2068.74
1877.93	1860.96	1975.3
867.81	2334.68	2535.1
1831.49	1713.11	1604.68
460.69	3274.77	2172.99
2374.98	3346.98	975.31
2271.89	3482.97	946.7
1783.64	1597.99	2261.31
198.83	3250.45	2445.08
1494.63	2072.59	2550.51
1597.03	1921.52	2126.76
1598.93	1921.08	1623.33
1243.13	1814.07	3441.07
2336.31	2640.26	1599.63
354	3300.12	2373.61 
2144.47	2501.62	591.51 
426.31	3105.29	2057.8   
1507.13	1556.89	1954.51 
343.07	3271.72	2036.94   
2201.94	3196.22	935.53   
2232.43	3077.87	1298.87   
1580.1	1752.07	2463.04  
1962.4	1594.97	1835.95  
1495.18	1957.44	3498.02  
1125.17	1594.39	2937.73   
24.22	3447.31	2145.01  
1269.07	1910.72	2701.97  
1802.07	1725.81	1966.35   
1817.36	1927.4	2328.79  
1860.45	1782.88	1875.13   
];
 [IDX,C,SUMD,D] = kmeans(data,4);
plot3(data(:,1),data(:,2),data(:,3),'*');
grid;
D=D'
minD=min(D);
index1 = find(D(1,:) ==min(D))
index2 = find(D(2,:) ==min(D))
index3 = find(D(3,:) ==min(D))
index4 = find(D(4,:) ==min(D))
line(data(index1,1),data(index1,2),data(index1,3),'linestyle', 'none','marker','*','color','g');
line(data(index2,1),data(index2,2),data(index2,3),'linestyle', 'none','marker','*','color','r');
line(data(index3,1),data(index3,2),data(index3,3),'linestyle', 'none','marker','+','color','b');
line(data(index4,1),data(index4,2),data(index4,3),'linestyle', 'none','marker','+','color','y');
title('C均值聚类分析图');
xlabel('第一特征坐标');
ylabel('第二特征坐标');
zlabel('第三特征坐标');
C