%%單連接法 --最近鄰近點分類  （可改 完全連接法    平均連接法）
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