f4.m

K-MEANS聚类算法

function [obj_fcn] =f4(center,cluster_n,data,expo)
dist = distfcm(center, data);       % fill the distance matrix
tmp = dist.^(-2/(expo-1));      % calculate new U, suppose expo != 1
U= tmp./(ones(cluster_n, 1)*sum(tmp));
mf = U.^expo;
obj_fcn = sum(sum((dist.^2).*mf));  % objective function