资源描述：
《模糊聚类分析步骤.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、原始数据矩阵X=欧式距离相似系数法标准化矩阵明氏距离距离法求分类对象的相似度切比雪夫距离主观评分法模糊相似矩阵R（1）传递闭包法截矩阵直接聚类法布尔矩阵法等价关系矩阵传递闭包法进行聚类（求动态聚类图）根据（0,1）的不同取值分布不同的类。注释（1）：模糊相似矩阵只具有自反性和对称性，不具有传递性，求截矩阵的前提是R是X上的的模糊等价关系。所以要先求得R传递闭包，将模糊相似矩阵转化为模糊等价矩阵。雨量站问题原始数据矩阵：（重要定理：设RF(XX)是相似关系(即R是自反、对称模糊关系)，则e(R)=t(R),即模糊相似关系的传递
2、闭包就是它的等价闭包。）Y的传递闭包（即Y的等价矩阵）：求截矩阵，在程序中我用的k代替了。K=1时，x1,x2,x3,…x11,各成一类，将11个雨量站分成11类。K=0.9095时，将11个雨量站分为10类，X8,X11为一类，其余各自一类。分8类，将x2,x5,x8,x11分一类，其余各自一类分6类，x2x3,x5,x8,x9x11为一类，其余各自一类。分4类，x1，x2,x3,x5,x7，x8,x9x11为一类，其余各自一类。分4类，x1,x3x2x7x8x9x11为一类，x2x4x5为一类，x6一类，x10一类。分3
3、类，x2x4x5x6为一类，x1x3x7x8x9x11一类，x10一类。分2类，x2x4x5x6x10一类，x1x3x7x8x9x11一类分2类，x1x2x4x5x6x10一类，x3x8x9x11一类.分1类。程序一：标准化矩阵：functionY=bzh1(X)[a,b]=size(X);C=max(X);D=min(X);Y=zeros(a,b);fori=1:aforj=1:bY(i,j)=(X(i,j)-D(j))/(C(j)-D(j));%平移极差变化进行数据标准化endendfprintf('标准化矩阵如下：Y=
4、');disp(Y)end程序二：求模糊相似矩阵：functionR=biaod2(Y,c)[a,b]=size(Y);Z=zeros(a);R=zeros(a);fori=1:aforj=1:afork=1:bZ(i,j)=abs(Y(i,k)-Y(j,k))+Z(i,j);R(i,j)=1-c*Z(i,j);%绝对值减数法--欧氏距离求模糊相似矩阵endendendfprintf('模糊相似矩阵如下：R=');disp(R)end程序三：计算传递闭包：functionB=cd3(R)a=size(R);B=zer
5、os(a);flag=0;whileflag==0fori=1:aforj=1:afork=1:aB(i,j)=max(min(R(i,k),R(k,j)),B(i,j));%R与R内积，先取小再取大endendendifB==Rflag=1;elseR=B;%循环计算R传递闭包endend程序四：求截矩阵：function[Dk]=jjz4(B)L=unique(B)';a=size(B);D=zeros(a);form=length(L):-1:1k=L(m);fori=1:aforj=1:aifB(i,j)>=kD(i
6、,j)=1;elseD(i,j)=0;%求?截距阵，当bij≥?时，bij(?)=1；当bij＜?时，bij(?)=0endendendfprintf('当分类系数k=：');disp(L(m));fprintf('所得截距阵为：');disp(D);end
