1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
| opts = statset('Display','final');
[idx,ctrs] = kmeans(Y,4,...
'Distance','city',...
'Replicates',9,...
'Options',opts);
plot(Y(idx==1,1),Y(idx==1,2),Y(idx==1,3),Y(idx==1,4),'r.','MarkerSize',12)
hold on
plot(Y(idx==2,1),Y(idx==2,2),Y(idx==2,3),Y(idx==2,4),'b.','MarkerSize',12)
hold on
plot(Y(idx==3,1),Y(idx==3,2),Y(idx==3,3),Y(idx==3,4),'y.','MarkerSize',12)
hold on
plot(Y(idx==4,1),Y(idx==4,2),Y(idx==4,3),Y(idx==4,4),'g.','MarkerSize',12)
plot(ctrs(:,1),ctrs(:,2),ctrs(:,3),ctrs(:,4),'kx',...
'MarkerSize',12,'LineWidth',2)
plot(ctrs(:,1),ctrs(:,2),ctrs(:,3),ctrs(:,4),'ko',...
'MarkerSize',12,'LineWidth',2)
legend('Cluster 1','Cluster 2','Cluster 3','Cluster 4','Centroids',...
'Location','NW') |
Partager