clusternormalize.m

混合高斯模型 对于给定的数据

function mixture1=ClusterNormalize(mixture)

   s = sum([mixture.cluster(:).pb]);
   for k=1:mixture.K
      mixture.cluster(k).pb = mixture.cluster(k).pb/s;     % normalize the probility
      mixture.cluster(k).invR = inv(mixture.cluster(k).R); % the inverse of covariance matrix R
      mixture.cluster(k).const = -(mixture.M*log(2*pi) +log(det(mixture.cluster(k).R)))/2; % ??
   end
   mixture1=mixture;

end