利用Matlab软件实现聚类分析.docxVIP

利用Matlab软件实现聚类分析 §8.利用matlab和spss软件实现聚类分析 1.用MATLAB编程 运用matlab中的一些基本矩阵计算方法，通过自己编程实现聚类算法，在此只讨论根据最短距离规则聚类的方法。 min1。M——找到矩阵的最小值，返回最小值的行和列以及值Min2的大小M——比较两个数字的大小并返回较小的值 std1.m——用极差标准化法标准化矩阵 ds1。M——用绝对值距离法计算距离矩阵 cluster.m——应用最短距离聚类法进行聚类分析 打印1。M——调用每个子函数以显示聚类结果 假设距离矩阵为vector，a阶，矩阵中最大值为max，令矩阵上三角元素等于max 聚类时间=A-1，以下步骤为A-1周期： 求改变后矩阵的阶数,计作c 找到矩阵的最小值，并返回最小值所在的行E和列f以及值的大小g forl=1:c,为vector(c+1,l)赋值，产生新类 设C+1列中的元素、E行和f行中的所有元素以及E列和f列中的所有元素为 %std1.m,用极差标准化法标准化矩阵 函数std=std1（向量） max=max(vector);%对列求最大值 min=min（向量）； [a,b]=size(vector);%矩阵大小,a为行数,b为列数 fori=1:a forj=1:b std（i，j）=（向量（i，j）-min（j））/（最大（j）-min（j））； %ds1.m,用绝对值法求距离 函数D=ds1（向量）； [a,b]=size(vector); d=零（a）； fori=1:a forj=1:a fork=1:b d（i，j）=d（i，j）+abs（向量（i，k）-向量（j，k））； fprintf(绝对值距离矩阵如下：\n); 高级惩教主任（丁） %min1.m,求矩阵中最小值，并返回行列数及其值 函数[v1，v2，v3]=min1（向量）；%V1是行数，V2是列数，V3是其值 [v,v2]=min(min(vector)); [v，v1]=min（min（向量））； v3=min(min(vector)); %min2。m、 比较两个数字的大小并返回较小的值 functionv1=min(v2,v3); ifv2v3 %cluster.m,最短距离聚类法 functionresult=聚类（向量）； [a,b]=size(vector); max=max（max（向量））； fori=1:a forj=i:b vector(i,j)=max; fork=1：（b-1） [c,d]=size(vector); Fprintf（群集%G:\n，K）； [e,f,g]=min1(vector); Fprintf（最小值=%G，将区域%G和区域%G合并为一个类，并记录为G%G\n\n，G，e，F，C+1）； forl=1:c 向量（c+1，l）=min2（向量（e，l），向量（f，l））； vector(c+1,l)=min2(vector(l,e),vector(l,f)); 向量（1:c+1，c+1）=最大值； vector(1:c+1,e)=max; 向量（1:c+1，f）=最大值； vector(e,1:c+1)=max; 向量（f，1:c+1）=最大值； %print1,调用各子函数 functionprint=print1（文件名，a，b）；%A是区域数，B是指标数 fid=fopen(filename,r) 向量=fscanf（fid，%g，[ab]）； fprintf(标准化结果如下：\n) v1=std1（向量） v2=ds1(v1); 集群（v2）； print1(fname,9,7) 2.直接调用matlab函数 2.1调用函数 层次聚类的计算步骤： ①计算n个样本两两间的距离{dij}，记d ② 构造n个类，每个类只包含一个样本； ③合并距离最近的两类为一新类； ④ 计算新类与当前类之间的距离；如果类数等于1，则转到5）；否则返回3）；⑤ 绘制聚类图； ⑥决定类的个数和类；