阿荣带你玩转K均值分类【matlab算法】

[filen, pathn] = uigetfile({'*.data'},'请选择待处理数据');
if isequal(filen,0)
    errordlg('没有选择文件.','温馨提示');
else
    %     [~,tname,~]     = fileparts(filen);
    filename        = strcat(pathn,filen);
    x=load(filename);%从文件中取数据
    x=x';
    [m,n]=size(x);
    k=3;
    G1=[];%暂存G1类
    G2=[];%暂存G2类
    G3=[];%暂存G3类
    g1=[];%G1类
    g2=[];%G2类
    g3=[];%G3类
    while isempty(g1)||(~(isequal(g1,G1))&&~(isequal(g2,G2))&&(~(isequal(g3,G3))))
        G1=g1;
        G2=g2;
        G3=g3;
        g=[];%存放类中心点
        if isempty(g1)
            g=x(1:2,1:3);%初次存放G1类、G2类、G3类的中心点
        else
            g=[g mean(G1,2)];
            g=[g mean(G2,2)];
            g=[g mean(G3,2)];
        end
        g1=[];
        g2=[];
        g3=[];
        for i=1:n %遍历所有数据
            min=(norm(g(1,1)-x(:,i)))^2;
            t=1;%记录数据应该在哪个类里
            for j=1:k %计算每个数据到类中心点的距离
                distant=(norm(g(:,j)-x(:,i)))^2 ;
                if(distant<min)
                    min=distant;
                    t=j;
                end
            end
            if t==1
                g1=[g1 x(:,i)];
            elseif t==2
                g2=[g2 x(:,i)];
            elseif t==3
                g3=[g3 x(:,i)];
            end
        end
    end
    figure('Name','欢迎关注matlab爱好者','Color','w')
    plot(g1(1,:),g1(2,:),'r*');
    hold on;
    plot(g2(1,:),g2(2,:),'b*');
    plot(g3(1,:),g3(2,:),'g*');
    hold off;
    title('K均值算法演示');
end