main.m

matlab实现的k均值聚类算法

clc
clear all;
close all;
load x_data
Data=x_data';
% Data=[0 0 1 1 1 2 2 3 6 6 7 7 7 8 8 8 8 9 9 9
%       0 1 0 1 2 1 2 2 6 7 6 7 8 6 7 8 9 7 8 9;];
[DataRow,DataColumn]=size(Data);
% model=pca(Data,2)
Step=20;
NumKind=5;
Center=Data(:,200:NumKind+199);
[KindData,KindNum,E]=Clustering(Center,Data);
NewCenter=CaculateCenter(KindData,KindNum,DataRow);
while (sum(sum(NewCenter~=Center))) & Step
    Center=NewCenter;
    [KindData,KindNum,E]=Clustering(Center,Data);
    NewCenter=CaculateCenter(KindData,KindNum,DataRow);
    Step=Step-1;
end

for i=1:NumKind
    KindData(:,1:KindNum(i),i) 
    pause
    clc
end
% 
% if DataRow~=NumKind
if (sum(sum(NewCenter~=Center)))
    disp('The result can not display on the two-dimensional figure ')
else
%     XMin=min(Data(1,:));
%     YMin=min(Data(2,:));
%     YMax=max(Data(2,:));
%     XMax=max(Data(1,:));
%     KindDisplay(KindData,KindNum,XMin,XMax,YMin,YMax)
end
E