matlab、lingo程序代码14-模糊聚类(聚类分析)

模糊聚类functionc=fuz_hc(a,b)%模糊矩阵的合成运算程序%输入模糊矩阵a,b,输出合成运算结果cm=size(a,1);n=size(b,2);p=size(a,2);%错误排除ifsize(a,2)~=size(b,1)disp('输入数据错误！');return;end%合成运算fori=1:mforj=1:nfork=1:ptemp(k)=min(a(i,k),b(k,j));endc(i,j)=max(temp);endenddisp('模糊矩阵a与b作合成运算后结果矩阵c为：');c%求模糊等价矩阵functionr_d=mhdj(r)[m,n]=size(r);fori=1:nforj=1:nfork=1:nr1(i,j,k)=min(r(i,k),r(k,j));endr1max(i,j)=r1(i,j,1);endendfori=1:nforj=1:nfork=1:nifr1(i,j,k)>r1max(i,j)r1max(i,j)=r1(i,j,k);endendr_d(i,j)=r1max(i,j);endend%模糊聚类程序functionf=mujl(x,lamda)%输入原始数据以及lamda的值iflamda>1disp('error!')%错误处理end[n,m]=size(x);y=pdist(x);disp('欧式距离矩阵：');dist=squareform(y)%欧氏距离矩阵dmax=dist(1,1);fori=1:nforj=1:nifdist(i,j)>dmaxdmax=dist(i,j);endendenddisp('处理后的欧氏距离矩阵，其特点为每项元素均不超过1：');sdist=dist/dmax%使距离值不超过1disp('模糊关系矩阵：');r=ones(n,n)-sdist%计算对应的模糊关系矩阵t=mhdj(r);le=t-r;whileall(all(le==0)==0)==1%如果t与r相等，则继续求r乘以rr=t;t=mhdj(r);le=t-r;enddisp('模糊等价矩阵为：')tfori=1:nk=1;forj=1:nift(i,j)>=lamdagroup(i,k)=j;k=k+1;endendenddisp('聚类结果如下(数字0为自动填充数据，不是样本序号)：')group(1,:)fori=2:nk=0;forj=1:i-1ifall(group(i,:)==group(j,:))==1%两行值完全相等，不输出k=1;break;endendifk==0disp(group(i,:))%仅输出不重复的分类endend%模糊聚类程序functionf=mujl(x,lamda)%输入原始数据以及lamda的值iflamda>1disp('error!')%错误处理end[n,m]=size(x);xmax=max(x);xmin=min(x);x=(x-xmin(ones(n,1),:))./(xmax(ones(n,1),:)-xmin(ones(n,1),:))y=pdist(x);disp('欧式距离矩阵：');dist=squareform(y)%欧氏距离矩阵dmax=dist(1,1);fori=1:nforj=1:nifdist(i,j)>dmaxdmax=dist(i,j);endendenddisp('处理后的欧氏距离矩阵，其特点为每项元素均不超过1：');sdist=dist/dmax%使距离值不超过1disp('模糊关系矩阵：');r=ones(n,n)-sdist%计算对应的模糊关系矩阵t=mhdj(r);le=t-r;whileall(all(le==0)==0)==1%如果t与r相等，则继续求r乘以rr=t;t=mhdj(r);le=t-r;enddisp('模糊等价矩阵为：')tfori=1:nk=1;forj=1:nift(i,j)>=lamdagroup(i,k)=j;k=k+1;endendenddisp('聚类结果如下(数字0为自动填充数据，不是样本序号)：')group(1,:)gru_val=1;fori=2:nk=0;forj=1:i-1ifall(group(i,:)==group(j,:))==1%两行值完全相等，不输出k=