模糊聚类matlab源程序代码

%模糊聚类程序 function f=mujl(x,lamda) %输入原始数据以及lamda的值 if lamda>1 disp('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); for i=1:n for j=1:n if dist(i,j)>dmax dmax=dist(i,j); end end end disp('处理后的欧氏距离矩阵，其特点为每项元素均不超过1：'); sdist=dist/dmax %使距离值不超过1 disp('模糊关系矩阵：'); r=ones(n,n)-sdist %计算对应的模糊关系矩阵 t=mhdj(r); le=t-r; while all(all(le==0)==0)==1 %如果t与r相等，则继续求r乘以r r=t; t=mhdj(r); le=t-r; end disp('模糊等价矩阵为：') t for i=1:n k=1; for j=1:n if t(i,j)>=lamda group(i,k)=j; k=k+1; end end end disp('聚类结果如下（数字0为自动填充数据，不是样本序号）：') group(1,:) gru_val=1; for i=2:n k=0; for j=1:i-1 if all(group(i,:)==group(j,:))==1 %两行值完全相等，不输出 k=1;break; end end if k==0 disp('第i类样本序号：'),i gru_val=gru_val+1; disp(group(i,:)) %仅输出不重复的分类 end end gru_val