模糊數(shù)學(xué)在聚類(lèi)分析中的作用(matlab代碼)

1、function [M,N] = Example8_11 X=[1.8 2.1 3.2 2.2 2.5 2.8 1.9 2.0; 95 99 101 103 98 102 120 130; 0.15 0.21 0.18 0.17 0.16 0.20 0.09 0.11]; X=X' %X=[80 10 6 2;50 1 6 4;90 6 4 6;40 5 7 3;10 1 2 4

2、][M,N]=fuzzy_jlfx(4,5,X); end %% function [M,N]=fuzzy_jlfx(bzh,fa,X)%得到聚類(lèi)結(jié)果 [X]=F_JlSjBzh(bzh,X);%數(shù)據(jù)標(biāo)準(zhǔn)化 [R]=F_JlR(fa,X);%建立相似矩陣 [A]=fuzzy_cdbb(R);%得到傳遞閉包矩陣 [Alamd]=fuzzy_lamdjjz(A);%得到lamdf截矩陣從而得到聚類(lèi)結(jié)果 [M,N]=F_JlDtjl(R);%

3、動(dòng)態(tài)聚類(lèi)并畫(huà)出聚類(lèi)圖 %% function [M,N]=F_JlDtjl(R) %clc; [A]=fuzzy_cdbb(R); U=unique(A); L=length(U); M=1:L; for i=L-1:-1:1[m,n]=find(A==U(i));N{i,1}=n;N{i,2}=m;A(m(1),:)=0;mm=unique(m);N{i,3}=mm;len=length(find(m==mm(1)));depth=l

4、ength(find(m==mm(2)));index1=find(M==mm(1));MM=[M(1:index1-1),M(index1+depth:L)]; %index2=find(MM==mm(2));M=M(index1:index1+depth-1);M=[MM(1:index2-1),M,MM(index2:end)]; end M=[1:L;M;ones(1,L)]; h=(max(U)-min(U))/L; fig

5、ure text(L,1,sprintf('%d',M(2,L))); text(L+1,1-h,sprintf('%d',L)); text(0,1,sprintf('%3.2f',1)); text(0,(1+min(U))/2,sprintf('%3.2f',(1+min(U))/2)); text(0,min(U),sprintf('%3.2f',

6、min(U))); hold on for i=L-1:-1:1m=N{i,2};end end end %% function[R]=F_JlR(cs,X)%定義函數(shù) %模糊聚類(lèi)分析建立模糊相似矩陣: [R]=F_JlR(cs,X) %X,數(shù)據(jù)矩陣 %cs=1,數(shù)量積法 %cs=2,夾角余弦法 %cs=3,相關(guān)系數(shù)法 %cs=4,指數(shù)相似系數(shù)法 %cs=5,最大最小法 %cs=6,算術(shù)平均最小法 %cs=7,幾何平均最小法 %cs=8

7、,一般歐式距離法 %cs=9,一般海明距離法 %cs=10,一般切比雪夫距離法 %cs=11,倒數(shù)歐式距離法 %cs=12,倒數(shù)海明距離法 %cs=13,倒數(shù)切比雪夫距離法 %cs=14,指數(shù)歐式距離法 %cs=15,指數(shù)海明距離法 %cs=16,指數(shù)切比雪夫距離法 [n,m]=size(X);%獲得矩陣的行列數(shù) R=[]; if(cs==1)maxM=0;pd=0;%數(shù)量積法for(i=1:n)for(j=1:n)if(j~=i)x=

8、0;for(k=1:m)x=x+X(i,k)*X(j,k);endif(maxM<x)maxM=x;endend;end;endif(maxM<0.000001)return;endfor(i=1:n)for(j=1:n)if(i==j)R(i,j)=1;else R(i,j)=0;for(k=1:m)R(i,j)=R(i,j)+X(i,k)*X(j,k);endR(i,j)=R(i,j)/maxM;if(R(i,j)<

9、;0)pd=1;endendend;endif(pd)for(i=1:n)for(j=1:n)R(i,j)=(R(i,j)+1)/2;end;end;end elseif(cs==2)%夾角余弦法for(i=1:n)for(j=1:n)xi=0;xj=0;for(k=1:m)xi=xi+X(i,k)^2;xj=xj+X(j,k)^2;ends=sqrt(xi*xj);R(i,j)=0;for(k=1:m)R(i,j)=R(i,j)+X(