clear all;
clc;
%可以輸入任意的一個樣本集，并制定要分多少類k
filename=input('請輸入識別數(shù)據(jù)文件名：','s');
k=input('請輸入想要分的類數(shù)：');
%打開制定的數(shù)據(jù)集，讀到矩陣Temp_Data中
Temp_Data=xlsread(filename);
%得到數(shù)據(jù)的總數(shù)量
[tk,temp]=size(Temp_Data);
%創(chuàng)建連個空的舉證
std_data=[];
other_data=[];
for Ans=0:(k-1)
    std_data=[std_data;Temp_Data(tk*Ans/k+1,:)];
    other_data=[other_data;Temp_Data(tk*Ans/k+2:tk*(Ans+1)/k,:)];
end

std_data(:,1)=[1:k];
other_data(:,1)=0;

Sum=zeros(k,1);
for ii=1:(tk-k)
    for jj=1:k
        Sum(jj)=sqrtm(sum((other_data(ii,:)-std_data(jj,:)).^2));
    end
    [Result,Min]=min(Sum);
    other_data(ii,1)=Min;
    std_data(Min)=(std_data(Min)+other_data(ii))/2;
end
last_data=[std_data;other_data];