ex16.m

花蝴蝶,利用HCM分群,模糊分群,hard-C mean

clc;
clear;
U=[1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 
   0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1];
U1=zeros(2,17);
x=[1 2 3 4 5 6 7 8 15 16 16 17 18 19 20 21 22];
y=[7 3.5 1 2 5 1.5 4 3 3 1.5 4 1.5 3.5 5 1.5 7 4];
num=0;
den=0;
e1=0;
error=100;

plot(x,y,'*')

while error>0.5
U1=zeros(2,17);
% update vi
for i=1:2
    for k=1:17
        num=num+x(k)*U(i,k);
        den=den+U(i,k);
    end
    v(i)=num/den
end

% update U
for i=1:2
    for k=1:17
        e=sqrt((x(k)-v(i))*(x(k)-v(i)))
        c=1:2;
        c(i)=[];
        for j=c
            e1=sqrt((x(k)-v(j))*(x(k)-v(j)));
        end
        min_e=min(e1)
        if e<=min_e
            U1(i,k)=1
        else
            U1(i,k)=0
        end
    end
end
%sprintf('材%d