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

Forpersonaluseonlyinstudyandresearch;mercialuse§,通过自己编程实现聚类算法,在此只讨论根据最短距离规则聚类的方法。调用函数:——求矩阵最小值,——比较两数大小,————————调用各子函数,显示聚类结果聚类分析算法假设距离矩阵为vector,a阶,矩阵中最大值为max,令矩阵上三角元素等于max聚类次数=a-1,以下步骤作a-1次循环:求改变后矩阵的阶数,计作c求矩阵最小值,返回最小值所在行e和列f以及值的大小gforl=1:c,为vector(c+1,l)赋值,产生新类令第c+1列元素,第e行和第f行所有元素为,第e列和第f列所有元素为max源程序如下:%,用极差标准化法标准化矩阵functionstd=std1(vector)max=max(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%,用绝对值法求距离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,k));endendendfprintf('绝对值距离矩阵如下:\n');disp(d)%,求矩阵中最小值,并返回行列数及其值function[v1,v2,v3]=min1(vector);%v1为行数,v2为列数,v3为其值[v,v2]=min(min(vector'));[v,v1]=min(min(vector));v3=min(min(vector));%,比较两数大小,返回较小的值functionv1=min(v2,v3);ifv2>v3v1=v3;elsev1=v2;end%,最短距离聚类法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次聚类:\n',k);[e,f,g]=min1(vector);fprintf('最小值=%g,将第%g区和第%g区并为一类,记作G%g\n\n',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