matlab實現FCM算法

% author:wangjunzuo
% date:2019/5/21
% fuction:fcm algrithmn

load data 
load label
maxgen = 100;  %??????
m = 2;         %2??
threshold = 10e-1000;   %?????
cluster_n = 3;       %????

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%???

[data_row,data_col] =size(data);
U=rand(cluster_n,data_row); %????01?????
col_sum=sum(U);    %???????
U=U./repmat(col_sum,cluster_n,1);  %????????????1
J(1)=0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i = 1:maxgen
    mf = U.^m;
    center = mf*data./(repmat(sum(mf,2),1,data_col)); %??????
    dist= zeros(size(center, 1), size(data, 1));
    for k=1:size(center,1)
        dist(k,:) = sqrt(sum(((data-repmat(center(k,:),data_row,1)).^2),2)');
    end
     J(i)= sum(sum((dist.^2).*mf));  %??????
     tmp = dist.^(-2/(m-1)); 
     U=tmp./repmat(sum(tmp),cluster_n,1);  %???????
     if i > 1               %????????????????
		if abs(J(i) - J(i-1)) < threshold 
            break; 
        end,
     end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?????
[max_vluae,index]=max(U);
index=index';
count = 0;
for i=1:data_row
    if index(i)==label(i)
        count = count + 1;
    end
end
accurcy = count / data_row;