有没有k-means聚类分析做辐射的？

无敌泡

如题，最近写毕业论文，需要用聚类分析做一下全国辐射的分区，无奈MATLAB水平差，想问下哪位大神可以指点一二？多谢

问题1、源数据该怎么处理？是这个格式么：站点号，辐射值，经度，纬度。

问题2、MATLAB自带的kmeans函数，调用的话，参数怎么设置？怎么体现地理位置信息？
          function [idxbest, Cbest, sumDbest, Dbest] = kmeans(X, k, varargin)
问题3、网上的这个例子，如何修改可以在我这里使用？这个X是随机的，我怎么把我的源数据放进去？

1. %随机获取150个点
   
2. X = [randn(50,2)+ones(50,2);randn(50,2)-ones(50,2);randn(50,2)+[ones(50,1),-ones(50,1)]];
   
3. 
   
4. opts = statset('Display','final');
   
5. 
   
6. %调用Kmeans函数
   
7. %X N*P的数据矩阵
   
8. %Idx N*1的向量,存储的是每个点的聚类标号
   
9. %Ctrs K*P的矩阵,存储的是K个聚类质心位置
   
10. %SumD 1*K的和向量,存储的是类间所有点与该类质心点距离之和
    
11. %D N*K的矩阵，存储的是每个点与所有质心的距离;
    
12. 
    
13. [Idx,Ctrs,SumD,D] = kmeans(X,3,'Replicates',3,'Options',opts);
    
14. 
    
15. %画出聚类为1的点。X(Idx==1,1),为第一类的样本的第一个坐标；X(Idx==1,2)为第二类的样本的第二个坐标
    
16. plot(X(Idx==1,1),X(Idx==1,2),'r.','MarkerSize',14)
    
17. hold on
    
18. plot(X(Idx==2,1),X(Idx==2,2),'b.','MarkerSize',14)
    
19. hold on
    
20. plot(X(Idx==3,1),X(Idx==3,2),'g.','MarkerSize',14)
    
21. 
    
22. %绘出聚类中心点,kx表示是圆形
    
23. plot(Ctrs(:,1),Ctrs(:,2),'kx','MarkerSize',14,'LineWidth',4)
    
24. plot(Ctrs(:,1),Ctrs(:,2),'kx','MarkerSize',14,'LineWidth',4)
    
25. plot(Ctrs(:,1),Ctrs(:,2),'kx','MarkerSize',14,'LineWidth',4)
    
26. 
    
27. legend('Cluster 1','Cluster 2','Cluster 3','Centroids','Location','NW')
    
28. 
    
29. Ctrs
    
30. SumD

复制代码

多谢多谢

无敌泡

木有人来回答，，，，

无敌泡

木有人，是我的问题太不明显了吗

无敌泡

木有人回答，求版主删除此贴，我回头自己再摸索去。等我摸索好了，我再把我所有的MATLAB的代码发上来。@二爷名声在外

无敌泡

楼主当了逃兵，放弃用MATLAB做聚类分析了，改用SPSS，结果凑合

That_gal

我也要研究kme ans，不过刚起步，想先试下spss，但是我没有软件，可不可以给我发一个，万分感谢，113871631@qq.com

无敌泡

That_gal 发表于 2015-5-11 20:08
我也要研究kme ans，不过刚起步，想先试下spss，但是我没有软件，可不可以给我发一个，万分感谢，

只是一个EXE文件，要不我给你共享到网盘吧。不过这个SPSS有点缺陷，就是不太准，还需要自己调整

无敌泡

That_gal 发表于 2015-5-11 20:08
我也要研究kme ans，不过刚起步，想先试下spss，但是我没有软件，可不可以给我发一个，万分感谢，

http://yunpan.cn/cj5EkCBefCqT6  访问密码 c810

That_gal

谢谢啦

guoqu510

正在学习，还是谢谢