function dan_new2(train_file,cha_num,lei_num)      %train_file 是打開的輸入文件名 lie_num是類數


train=train_file(:,1:(2+cha_num));
%train=xlsread(train_file);

yang_num=size(train,1);

%########################### 一、準備工作 ##############################

%初始化聚類,第一維與第二維記錄樣本特征,第三維記錄類別號,初始化時共60個類別
for ii=1:yang_num
    julei(1,:,ii)=train(ii,:);
    for jj=2:yang_num
        julei(jj,:,ii)=0;
    end    
end
kk=size(julei,2);
%%%%%%%%%%%%%%%%%% 初始化聚類結束 %%%%%%%%%%%%%%%%%%%%%%%%

% dd用于存放任意兩個類中各個樣本間的距離
dd=9999; 

%%%%%%%%%%%%%%%%%%%%%%%%% 初始化距離矩陣 %%%%%%%%%%%%%%%%%%%%%%%%%%%%   
% d矩陣中存放所有各個樣本間的距離，如第一行第二列，是第一類到第二類的距離
% 對角線及對角線以下的無用，全部置為9999，那么選擇類與類的距離的時候便不會選到他們    

for ii=1:yang_num
    for jj=1:yang_num
        d(ii,jj)=9999;
    end
end
%%%%%%%%%%%%%%%%%%%%%%%%% 初始化距離矩陣結束 %%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%   有效聚類數      %%%%%%%%%%%%%%%%%%%%%%%%
%有效聚類是指內容不為空的類，比如第一類與第二類合并后，將第二類的樣本放入第一類
%而將第二類清零，此時第二類為無效聚類。
n=0;
for ii=1:yang_num
    if julei(1,1,ii)~=0
        n=n+1;
    end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 準備工作結束 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%############################ 二、開始聚類 ###############################
while n~=lei_num
    %算距離
    for ii=1:yang_num
        for jj=1:yang_num
            if julei(1,1,ii)~=0&julei(1,1,jj)~=0&ii<jj
                for mm=1:yang_num
                    for nn=1:yang_num
                        if julei(mm,1,ii)~=0&julei(nn,1,jj)~=0
                            dd(mm,nn)=sum((julei(mm,3:(2+cha_num),ii)-julei(nn,3:(2+cha_num),jj)).^2).^0.5;
                        end
                    end
                end
                d(ii,jj)=min(min(dd));
            end 
            dd=zeros(1);
            dd=9999;
        end
    end
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    %合并
    for ii=1:yang_num
        for jj=1:yang_num
            if d(ii,jj)==min(min(d))
                ii
                jj
                num1=0;num2=0;
                for mm=1:yang_num
                    if julei(mm,1,ii)~=0
                        num1=num1+1;
                    end
                    if julei(mm,1,jj)~=0
                        num2=num2+1;
                    end                        
                end
                number=1;
                for nn=(num1+1):(num1+num2)
                    julei(nn,:,ii)=julei(number,:,jj);
                    number=number+1;
                end
                julei(:,:,jj)=0;
            end
        end
    end
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    %初始化
    for ii=1:yang_num
        for jj=1:yang_num
            d(ii,jj)=9999;
        end
    end

    %計算有效聚類
    n=0;
    for ii=1:yang_num
        if julei(1,1,ii)~=0
            n=n+1;
        end
    end 
    
end
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   
%#############################              ###############################
%############################# 聚類全部結束  ###############################

%############################  三、后期處理 ############################
%挑出有效類
a=1;
for ii=1:yang_num
    
        if julei(1,1,ii)~=0
            num1=0;
                for mm=1:yang_num
                    if julei(mm,1,ii)~=0
                        num1=num1+1;
                    end                                       
                end
            Final_ju(1:num1,:,a)=[julei(1:num1,:,ii)];
            a=a+1;
            
        end
    
end



disp('聚類結果')
Final_ju

aa=size(Final_ju,1);
%當只有三個類的時候
%區聚類所屬并計算模板
if lei_num==3
    num1=0;num2=0;num3=0;
    for ii=1:aa
        if Final_ju(ii,1,1)~=0
            num1=num1+1;
        end
        if Final_ju(ii,1,2)~=0
            num2=num2+1;
        end
        if Final_ju(ii,1,3)~=0
            num3=num3+1;
        end
    end
end
    
if lei_num==3&cha_num==4
    figure
plot(Final_ju(1:num1,5,1),Final_ju(1:num1,4,1),'r*',Final_ju(1:num2,5,2),Final_ju(1:num2,4,2),'m*',Final_ju(1:num3,5,3),Final_ju(1:num3,4,3),'g*')

figure
colordef black;
plot3(Final_ju(1:num1,3,1),Final_ju(1:num1,4,1),Final_ju(1:num1,5,1),'r*',Final_ju(1:num2,3,2),Final_ju(1:num2,4,2),Final_ju(1:num2,4,2),'m*',Final_ju(1:num3,3,3),Final_ju(1:num3,4,3),Final_ju(1:num3,5,3),'g*')
grid
end