clear all;
clc;
%可以輸入任意的一個(gè)樣本集，并制定要分多少類(lèi)k
filename=input('請(qǐng)輸入識(shí)別數(shù)據(jù)文件名：','s');
k=input('請(qǐng)輸入想要分的類(lèi)數(shù)：');
%打開(kāi)制定的數(shù)據(jù)集，讀到矩陣Temp_Data中
Temp_Data=xlsread(filename);
%得到數(shù)據(jù)的總數(shù)量
[tk,temp]=size(Temp_Data);
%創(chuàng)建連個(gè)空的矩陣用于存儲(chǔ)k類(lèi)，和其他的欲分類(lèi)的數(shù)據(jù)
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; %講樣本的分類(lèi)標(biāo)志全部置零
%一下開(kāi)始分類(lèi)
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
%得到總的分類(lèi)結(jié)果，根據(jù)第一列的分類(lèi)標(biāo)志
last_data=[std_data;other_data];