function Cmat = fcm_calcC(Xdat,C,mfuz,Umat);
%   function Cmat = fcm_calcC(Xdat,K,mfuz,Umat);
%
%   計算或者更新聚類中心矩陣（聚類原型模式矩陣）
%   
%   Xdat ：原數(shù)據(jù)集，大小為(M,N)
%   C    ：聚類數(shù)目
%   M:數(shù)據(jù)的個數(shù)，N：數(shù)據(jù)的維數(shù)
%   mfuz ：加權指數(shù)m
%   Umat ：模糊劃分矩陣，大小為(C,M)
%   Cmat ：聚類中心矩陣，大小為(C,N)
%
%   
%Umat=[1,2,3,4;5,6,7,8;9,10,11,12]
%Xdat=[1,0;0,1;2,2;1,1]
%C=3;
%mfuz=2;
[M,N] = size(Xdat);
for c=1:C
    tmp(c,:) = (Umat(c,:).^mfuz);
end

for c=1:C
    for m=1:M
    nume(m,:) = tmp(c,m)*Xdat(m,:);
    end
    deno = sum(tmp(c,:));
    snume=sum(nume);
    Cmat(c,:) = snume/deno;
end

return %% end of function fcm_calcC()