模糊数学教学课件2.ppt

模糊聚类程序 function julei(data) %%%%%%%%%%%%%%%模糊聚类%%%%%%%%%%%%%%%%%%%%%%% DATAFORCLUS=data; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%基于模糊等价关系的模糊聚类%%%%%%%%%%%%%%%%%%%%%% %----------构造相似关系-----------% numrows=size(DATAFORCLUS,1); numcols=size(DATAFORCLUS,2); disp(请选择对象之间相似性统计量的方式: ); disp(1-相关系数法|2-夹角余弦法); wayforr_ij=input(请输入: ); switch wayforr_ij case 1, %-----------------------------------相关系数法 for i=1:numrows, for j=1:numrows, meani=mean(DATAFORCLUS(i,:));meanj=mean(DATAFORCLUS(j,:)); simiR(i,j)=sum((DATAFORCLUS(i,:)-meani).*(DATAFORCLUS(j,:)-meanj))/... (sqrt(sum((DATAFORCLUS(i,:)-meani).^2))*sqrt(sum((DATAFORCLUS(j,:)-meanj).^2))); end end case 2, %-----------------------------------夹角余弦法 for i=1:numrows, for j=1:numrows, simiR(i,j)=sum(DATAFORCLUS(i,:).*DATAFORCLUS(j,:))/... (sqrt(sum(DATAFORCLUS(i,:).*DATAFORCLUS(i,:)))*sqrt(sum(DATAFORCLUS(j,:).*DATAFORCLUS(j,:)))); end end end 模糊聚类程序 %-------改造成等价关系----------% sign=0; numselfmul=1; simiRk=eye(numrows); equi_tem=simiR; while sign==0, for i=1:numrows, for j=1:numrows, for c=1:numrows, rij_temp(c)=min([equi_tem(i,c) equi_tem(c,j)]); end simiRk(i,j)=max(rij_temp); end end %--------------% if sum(sum(simiRk-equi_tem,1))~=0, numselfmul=numselfmul+1; equi_tem=simiRk; else sign=1; break end %--------------% end if sign==1, disp(从相似矩阵到等价矩阵改造成功!!!); else disp(从相似矩阵到等价矩阵改造失败!!!); end equiR=simiRk; numclass=input(请输入聚类数: ); 模糊聚类程序 %---------在不同的截集水平进行聚类--------------% clasc=0; comp_vec(1,1:numrows)=0; index=0; clasc=0; tip=0; alpha=0; temnumeachclass=0; while (tip==0), %alpha=input(请输入进行分类的截集水平λ: ); %alpha=0.5; %调试 if (alpha0 || alpha1), error(您输入的截集水平λ不符合分类要求-执行结束!!!); end comp_arr=o