function Umat = fcm_calcU(mfuz,Dmat);
%   function Umat = fcm_calcU(mfuz,Dmat);
%
%   計(jì)算或者更新模糊劃分矩陣
%
%   mfuz :加權(quán)因子m
%   Dmat ：距離矩陣，大小為 (C,M)，這里使用歐氏距離
%   Umat ：模糊劃分矩陣，大小為(C,M)
%
%   

[C,M] = size(Dmat);
for i=1:M
    temp=0;
     for k=1:C
         if(Dmat(k,i)==0) temp=temp+1; end
     end
     if(temp==0)
          for k=1:C
             nume = Dmat(k,i);
            deno = (nume./Dmat(:,i)).^(2/(mfuz-1));
            Umat(k,i) = 1.0/sum(deno);
         end
      else
             for k=1:C
                 if(Dmat(k,i)==0)
                     Umat(k,i)=1/temp;
                 else
                     Umat(k,i)=0;
        
                 end
             end
         end
     end

              
             
return %% end of function fcm_calcU.m