lei.m

用于 分层聚类法分类；原码为单连接法分类

%%单连接法 --最近邻近点分类  （可改 完全连接法    平均连接法）
function le=lei(x)
%用x为要分类的欧式距离（或其他）矩阵，
m=length(x);
%%找最短值前，先排除0
for i=1:m
    for j=1:m
        if x(i,j)==0
            x(i,j)=inf;
        end
    end
end
q=0;k=m;
for p=(m+1):2*m
%%保存分类矩阵   i,j表示国家代号；每一行表示被分为一组；k表示被分为一组后的新的国家代号，l表i，j的距离

%通过k的连接，可实现把几个国家分为一类。用l评定分类的距离
l=min(min(x));
if l~=inf
    [i,j]=find(x==l);
end
k=k+1;

q=q+1;
le(q,:)=[i,j ,k,l];
for n=1:size(x,2)
    if  n~=i|n~=j
       b(n)=min(x(i,n),x(j,n));  %%%%%改为b（n）=max(x(i,n),x(j,n))为完全连接法     
                                %%%%      b(n)=(x(i,n)+x(j,n))/2为平均连接法
  elseif j==p
      b(j)=inf;
       end
       
end
b(i)=inf;
b(j)=inf;
x(i,:)=inf;
x(j,:)=inf;
x(:,i)=inf;
x(:,j)=inf;
x(:,p)=b';
x(p,:)=inf;%%%%
end