nodepos=xlsread('nodes');
id=nodepos(:,1);
x=nodepos(:,2);
y=nodepos(:,3);
pos=nodepos(:,2:3);
% Plot nodes
scatter(x(id== 0), y(id == 0), 100, 'k', '*')
hold on
scatter(x(id== 1), y(id == 1), 100, 'g', 's')
hold on
scatter(x(id== 2), y(id == 2), 120, 'r', '>')
hold on
scatter(x(id== 3), y(id == 3), 130, 'g', 'o')
hold on
scatter(x(id== 4), y(id == 4), 140, 'c', 'p')
hold on
scatter(x(id== 5), y(id == 5), 150, 'm', 'x')
hold on
scatter(x(id== 6), y(id == 6), 160, 'y', '+')
hold on
scatter(x(id== 7), y(id == 7), 170, 'b', 'h')
hold on
scatter(x(id== 8), y(id ==8), 180, 'r', '^')
hold on
scatter(x(id== 9), y(id == 9), 190, 'b', 'd')
hold on
grid on
opts = statset('Display','final');
[idx,C1] =
kmeans(pos,2,'Distance','cityblock','Replicates',10,'Options',opts);
figure;
plot(pos(idx==1,1),pos(idx==1,2),'rp','MarkerSize',12)
hold on
plot(pos(idx==2,1),pos(idx==2,2),'bh','MarkerSize',12)
plot(C1(:,1),C1(:,2),'kx','MarkerSize',15,'LineWidth',3)
legend('Cluster 1','Cluster 2','Centroids','Location','NE')
title 'Clusters'
hold off