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('輸入數(shù)據(jù)錯誤！');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)%輸入原始數(shù)據(jù)以及l(fā)amda的值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相等，則繼續(xù)求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('聚類結果如下(數(shù)字0為自動填充數(shù)據(jù)，不是樣本序號)：')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)%輸入原始數(shù)據(jù)以及l(fā)amda的值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相等，則繼續(xù)求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('聚類結果如下(數(shù)字0為自動填充數(shù)據(jù)，不是樣本序號)：')group(1,:)gru_val=1;fori=2:nk=0;forj=1:i-1ifall(group(i,:)==group(j,:))==1%兩行值完全相等，不輸出k=