《模糊数学在聚类分析中的作用(matlab代码)》由会员分享，可在线阅读，更多相关《模糊数学在聚类分析中的作用(matlab代码)（8页珍藏版）》请在金锄头文库上搜索。
1、function M,N = Example8_11X=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 4M,N=fuzzy_jlfx(4,5,X);end%function M,N=fuzzy_jlfx(bzh,fa,X)%得到聚类结果X=F_JlSjBzh(bzh,X);%数据标准化R=F_JlR(fa,X);%建立相似矩阵A=fuzzy_
2、cdbb(R); %得到传递闭包矩阵Alamd=fuzzy_lamdjjz(A);%得到 lamdf 截矩阵从而得到聚类结果M,N=F_JlDtjl(R);%动态聚类并画出聚类图%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:1m,n=find(A=U(i);Ni,1=n;Ni,2=m;A(m(1),:)=0;mm=unique(m);Ni,3=mm;len=length(find(m=mm(1);depth=length(find(m=mm(2);index1=f
3、ind(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);endM=1:L;M;ones(1,L);h=(max(U)-min(U)/L;figuretext(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+
4、min(U)/2);text(0,min(U),sprintf( %3.2f ,min(U);hold onfor i=L-1:-1:1m=Ni,2;n=Ni,1;mm=Ni,3;k=find(M(2,:)=mm(1); l=find(M(2,:)=mm(2); x1=M(1,k);y1=M(3,k);x2=M(1,l);y2=M(3,l);x=x1,x1,x2,x2;M(3,k,l)=U(i);M(1,k,l)=sum(M(1,k,l)/length(M(1,k,l);y=y1,U(i),U(i),y2;plot(x,y);text(i,1,sprintf( %d ,M(2,i);text(
5、M(1,k(1),U(i)+h*0.1,sprintf( %3.2f ,U(i);%d ,L-i);%定义函数: X=F_JlSjBzh(cs,X)text(L+1,1-i*h-h,sprintf( endaxis(0 L+1 min(U) max(U) axis off hold offendend %function X=F_JlSjBzh(cs,X) %模糊聚类分析数据标准化变换%数据矩阵 %cs=0, 不变换 ;cs=1, 标准差变换 ;cs=2, 极差变换 ;cs= 其它 , 最大值规格化if (cs=0) return ; endn,m=size(X);%获得矩阵的行列数if (c
6、s=1) %平移 ?标准差变换for (k=1:m)xk=0;for (i=1:n)xk=xk+X(i,k); end xk=xk/n;sk=0;for (i=1:n)sk=sk+(X(i,k)-xkF2;endsk=sqrt(sk/n);for (i=1:n)X(i,k)=(X(i,k)-xk)/sk;endendelseif (cs=2) %平移 ?极差变换for (k=1:m)xmin=X(1,k);xmax=X(1,k);for (i=1:n)if (xminX(i,k)xmin=X(i,k);endif (xmaxX(i,k)xmax=X(i,k);endendfor (i=1:n)
7、X(i,k)=(X(i,k)-xmin)/(xmax-xmin);endendelse %最大值规格化A=max(X); for (i=1:m) X(:,i)=X(:,i)/A(i);endendend%function R=F_JlR(cs,X)%定义函数%模糊聚类分析建立模糊相似矩阵: R=F_JlR(cs,X)%数据矩阵%cs=1, 数量积法%cs=2, 夹角余弦法%cs=3, 相关系数法%cs=4, 指数相似系数法%cs=5, 最大最小法%cs=6, 算术平均最小法%cs=7, 几何平均最小法%cs=8, 一般欧式距离法%cs=9, 一般海明距离法%cs=10, 一般切比雪夫距离法%c
8、s=11, 倒数欧式距离法%cs=12, 倒数海明距离法%cs=13, 倒数切比雪夫距离法%cs=14, 指数欧式距离法%cs=15, 指数海明距离法%cs=16, 指数切比雪夫距离法n,m=size(X);%获得矩阵的行列数R=;if (cs=1)maxM=0;pd=0; %数量积法for (i=1:n) for (j=1:n) if (j=i)x=0;for (k=1:m)x=x+X(i,k)*X(j,k); endif (maxMx)maxM=x;endend; end; endif (maxM0.000001) return ; endfor (i=1:n) for (j=1:n)if
9、 (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)0)pd=1; endend end; end if (pd) for (i=1:n) for (j=1:n)R(i,j)=(R(i,j)+1)/2;end; end; endelseif (cs=2) %夹角余弦法 for (i=1:n) for (j=1:n)xi=0;xj=0;for (k=1:m)xi=xi+X(i,k)A2;xj=xj+X(j,k)A2; s=sqrt(xi*xj);R(
10、i,j)=0;endfor (k=1:m)R(i,j)=R(i,j)+X(i,k)*X(j,k); R(i,j)=R(i,j)/s;end; endendelseif (cs=3) %相关系数法for (i=1:n) for (j=1:n)xi=0;xj=0;for (k=1:m)xi=xi+X(i,k);xj=xj+X(j,k); end xi=xi/m;xj=xj/m;xis=0;xjs=0;for (k=1:m)xis=xis+(X(i,k)-xi)A2;xjs=xjs+(XQ,k)-xjF2;ends=sqrt(xis*xjs);R(i,j)=0;for (k=1:m)R(i,j)=R
11、(i,j)+abs(X(i,k)-xi)*(X(j,k)-xj);endR(i,j)=R(i,j)/s;end; endelseif (cs=4) %指数相似系数法for (i=1:n) for (j=1:n)R(i,j)=0;for (k=1:m)xk=0;for (z=1:n)xk=xk+X(z,k); end xk=xk/n;sk=0;for (z=1:n)sk=sk+(X(z,k)-xkf2;endsk=sk/ n;R(i,j)=R(i,j)+exp(-0.75*(X(i,k)-X(j,k)/sk)A2);endR(i,j)=R(i,j)/m;end; endelseif (cs=7)
12、 %最大最小法 算术平均最小法 几何平均最小法for (i=1:n) for (j=1:n)fz=0;fm=0;for (k=1:m)if (X(j,k)X(i,k)x=X(i,k);else x=X(j,k); endfz=fz+x;endif (cs=5) %最大最小法for (k=1:m) if (X(i,k)X(j,k)x=X(i,k);else x=X(j,k); endfm=fm+x;endelseif (cs=6) for (k=1:m)fm=fm+(X(i,k)+X(j,k)/2;e n d %算术平均最小法else for (k=1:m)fm=fm+sqrt(X(i,k)*X
13、(j,k);end; end%几何平均最小法R(i,j)=fz/fm;end; endelseif (cs=10)C=0; %一般距离法for (i=1:n) for (j=i+1:n)d=0;if (cs=8) for (k=1:m)d=d+(X(i,k)-X(j,k)A2;endd=sqrt(d);%欧式距离elseif (cs=9) for (k=1:m)d=d+abs(X(i,k)-X(j,k);end%海明距离else for (k=1:m) if (dabs(X(i,k)-X(j,k)d=abs(X(i,k)-X(j,k);end; end; end%切比雪夫距离if (Cd)C=
14、d; endend; endC=1/(1+C);for (i=1:n) for (j=1:n)d=0;if (cs=8) for (k=1:m)d=d+(X(i,k)-X(j,k)A2;endd=sqrt(d);%欧式距离elseif (cs=9) for (k=1:m)d=d+abs(X(i,k)-X(j,k);end%海明距离else for (k=1:m) if (dabs(X(i,k)-X(j,k)d=abs(X(i,k)-X(j,k);end; end; end%切比雪夫距离R(i,j)=1-C*d;end; endelseif (cs=13)minM=Inf;%倒数距离法for (i=1:n) for (j=i+1:n)d=0;if (cs=11) for (k=1:m)d=d+(X(i,k)-X(j,k)A2;end
