clusterdemo.m

DNa聚类分析软件

clc,clear;
path(path,'c:\');
load demo.txt;
sum=0;
cofi=0;
cofj=0;
for i=1:1:20
  for j=1:1:20
     avei=mean(demo(i,:));
     avej=mean(demo(j,:));
     sum=0;
     cofi=0;
     cofj=0;
     for k=1:1:4
         sum=sum+abs(demo(i,k)-avei)*abs(demo(j,k)-avej);
         cofi=cofi+(demo(i,k)-avei)^2;
         cofj=cofj+(demo(j,k)-avej)^2;
     end
     cof(i,j)=sum/(sqrt(cofi)*sqrt(cofj));
 end
end

for i=1:1:20
    for j=1:1:20
        cof(i,j)=round(cof(i,j)*1000)/1000;
    end
end

for i=1:1:20
    for j=1:1:20
        cof(i,j)=round(cof(i,j)*100)/100;
    end
end

%传递闭包
R=0;  
result=cof;
for k=1:1:5
  flag=0;
  last=result;
  for i=1:1:20
    for j=1:1:20
         f=max(min(result(i,:),result(:,j)'));         
         result(i,j)=f; 
    end
  end
  for i=1:1:20
      for j=1:1:20
         if result(i,j)~=last(i,j)
             flag=flag+1;
         end
     end
  end

  k
  flag 
end

temp=result;
for thre=0:0.01:0.05
    thre=1-thre;
    result=temp;
     for j=1:1:20        
        for k=1:1:20
            if result(j,k)>=thre 
                result(j,k)=1;
            else
                result(j,k)=0;
            end
        end
    end    
end