资源描述：
《利用matlab和spss软件实现聚类分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、《计量地理学》（徐建华，高等教育出版社，2005）配套实习指导§8.利用Matlab和SPSS软件实现聚类分析1.用Matlab编程实现运用Matlab中的一些基本矩阵计算方法，通过自己编程实现聚类算法，在此只讨论根据最短距离规则聚类的方法。调用函数：min1.m——求矩阵最小值，返回最小值所在行和列以及值的大小min2.m——比较两数大小，返回较小值std1.m——用极差标准化法标准化矩阵ds1.m——用绝对值距离法求距离矩阵cluster.m——应用最短距离聚类法进行聚类分析print1.m——调用各子函数，显示聚类结
2、果聚类分析算法假设距离矩阵为vector，a阶，矩阵中最大值为max，令矩阵上三角元素等于max聚类次数=a-1,以下步骤作a-1次循环:求改变后矩阵的阶数,计作c求矩阵最小值，返回最小值所在行e和列f以及值的大小gforl=1:c,为vector(c+1,l)赋值，产生新类50《计量地理学》（徐建华，高等教育出版社，2005）配套实习指导令第c+1列元素，第e行和第f行所有元素为，第e列和第f列所有元素为max源程序如下：%std1.m,用极差标准化法标准化矩阵functionstd=std1(vector)max=ma
3、x(vector);%对列求最大值min=min(vector);[a,b]=size(vector);%矩阵大小,a为行数,b为列数fori=1:aforj=1:bstd(i,j)=(vector(i,j)-min(j))/(max(j)-min(j));endend%ds1.m,用绝对值法求距离functiond=ds1(vector);[a,b]=size(vector);d=zeros(a);fori=1:aforj=1:afork=1:bd(i,j)=d(i,j)+abs(vector(i,k)-vector(j
4、,k));endendendfprintf('绝对值距离矩阵如下：');disp(d)%min1.m,求矩阵中最小值，并返回行列数及其值function[v1,v2,v3]=min1(vector);%v1为行数，v2为列数，v3为其值[v,v2]=min(min(vector'));[v,v1]=min(min(vector));v3=min(min(vector));%min2.m，比较两数大小，返回较小的值functionv1=min(v2,v3);ifv2>v3v1=v3;else50《计量地理学》（徐建华，高
5、等教育出版社，2005）配套实习指导v1=v2;end%cluster.m,最短距离聚类法functionresult=cluster(vector);[a,b]=size(vector);max=max(max(vector));fori=1:aforj=i:bvector(i,j)=max;endend;fork=1:(b-1)[c,d]=size(vector);fprintf('第%g次聚类:',k);[e,f,g]=min1(vector);fprintf('最小值=%g,将第%g区和第%g区并为一类，记作G
6、%g',g,e,f,c+1);forl=1:cifl<=min2(e,f)vector(c+1,l)=min2(vector(e,l),vector(f,l));elsevector(c+1,l)=min2(vector(l,e),vector(l,f));endend;vector(1:c+1,c+1)=max;vector(1:c+1,e)=max;vector(1:c+1,f)=max;vector(e,1:c+1)=max;vector(f,1:c+1)=max;end%print1,调用各子函数funct
7、ionprint=print1(filename,a,b);%a为地区个数，b为指标数fid=fopen(filename,'r')vector=fscanf(fid,'%g',[ab]);fprintf('标准化结果如下：')v1=std1(vector)v2=ds1(v1);cluster(v2);%输出结果print1('fname',9,7)50《计量地理学》（徐建华，高等教育出版社，2005）配套实习指导2.直接调用Matlab函数实现2.1调用函数层次聚类法（HierarchicalClustering）的
8、计算步骤：①计算n个样本两两间的距离{dij}，记D②构造n个类，每个类只包含一个样本；③合并距离最近的两类为一新类；④计算新类与当前各类的距离；若类的个数等于1，转到5）；否则回3）；⑤画聚类图；⑥决定类的个数和类；Matlab软件对系统聚类法的实现（调用函数说明）：cluster从连接输出(link
