function CollectAnaly(X)
format long;
sz = size(X);
N = sz(1);         %樣本個數
n = sz(2);
D = zeros(n,n);

totalClass = N;
RecordClass = zeros(N,N+1);
RecordClass(:,1) = ones(N,1);
RecordClass(:,2) = 1:N;
disp('聚類前的N個類：');
disp(RecordClass);

while totalClass > 1
    minClaDist = inf;
    for i=1:totalClass
        for j=i+1:totalClass
            distClass = DistClass(X,RecordClass,i,j,RecordClass(i,1),RecordClass(j,1));
            if distClass < minClaDist
                minClaDist = distClass;
                i1 = i;
                j1 = j;
            end
        end
    end
    t1 = RecordClass(i1,1);
    t2 = RecordClass(j1,1);
    RecordClass(i1,(t1+2):(t1+t2+1)) = RecordClass(j1,2:(t2+1));
    RecordClass(i1,1) = RecordClass(i1,1) + RecordClass(j1,1);
    RecordClass(j1:(totalClass-1),:) = RecordClass((j1+1):totalClass,:);
    RecordClass(totalClass:N,:) = zeros(N-totalClass+1,N+1);
    totalClass = totalClass - 1;
    str1 = strcat('第',num2str(N - totalClass));
    str1 = strcat(str1, '次聚類，');
    str1 = strcat(str1, '第');
    str1 = strcat(str1, num2str(i1));
    str1 = strcat(str1, '類和第');
    str1 = strcat(str1, num2str(j1));
    str1 = strcat(str1, '類合并：');
    disp(str1);
    disp(RecordClass(1:totalClass,:))
end


function d = DistSamp(X1,X2)
format long;
d = sqrt(dot( X1-X2,X1-X2));
format short;

function d = DistClass(X,C1,l1,l2,n,m)
format long;
d = inf;
for i=1:n
    for j=1:m
        dc = DistSamp(X(C1(l1,i+1),:),X(C1(l2,j+1),:));
        if dc < d
            d = dc;
        end
    end
end
format short;