coordinate_file = 'taylor.xlsx';
[num,txt] = xlsread(coordinate_file);
popularity=num(1:end,17);
acousticness=num(1:end,8);
energy=num(1:end,10);
Max_lat=max(popularity);
Min_lat=min(popularity);
Max_lng=max(acousticness);
Min_lng=min(acousticness);
Max_jd=max(energy);
Min_jd=min(energy);
figure(1);
plot3(popularity,acousticness,energy,'b.','markersize',7);
% 聚类分析
[type,center]=kmeans(num,8);
rows=size(center,1);
for i=1:rows
fprintf('聚类号为%f 的聚类中心坐标为:%f,%f,为%f\n',i,center(i,1),center(i,2),center(i,3));
end
h=figure(2);
hold on;
plot3(popularity(type==1),acousticness(type==1),energy(type==1),'r+','markersize',5);
plot3(popularity(type==2),acousticness(type==2),energy(type==2),'bx','markersize',5);
plot3(popularity(type==3),acousticness(type==3),energy(type==3),'g*','markersize',5);
plot3(popularity(type==4),acousticness(type==4),energy(type==4),'m.','markersize',5);
plot3(popularity(type==5),acousticness(type==5),energy(type==5),'c+','markersize',5);
plot3(popularity(type==6),acousticness(type==6),energy(type==6),'m*','markersize',5);
plot3(popularity(type==7),acousticness(type==7),energy(type==7),'c+','markersize',5);
plot3(popularity(type==8),acousticness(type==8),energy(type==8),'y*','markersize',5);
plot3(center(:,2),center(:,1),center(:,3),'bo','markersize',16) ;
hold off
大概是什么问题呢?