function [loglike, pmat] = loglikeGMM(X, mus, sigmas, weights);

% 措阿 啊快矯救狼 雇膠媚俊 措茄 X狼 阿 單撈磐 氦磐狼 肺弊 快檔甫 府畔茄促. 
% 肚茄 阿 雇膠媚 己盒苞 單撈磐 痢俊 措茄 俺喊 膠內(nèi)絹牢 ( p(x_n|m)c_m )甫
% 器竊竅絆 樂綽 M*numpoints 青紡闌 府畔茄促.
%
% 澇仿 牢磊 :
% X : 青撈 d-瞞盔牢 漂隆 氦磐甸牢 單撈磐 青紡撈促.
% mus : d*M 青紡肺 弊 青籃 雇膠媚 己盒俊 措茄 乞閉 氦磐.
% sigmas : d*M 青紡肺 弊 青籃 雇膠媚 去欽 己盒俊 措茄 盒魂 氦磐.
% weights : M-瞞盔 氦磐肺 雇膠媚 啊吝摹.
%
% 免仿 牢磊 :
% loglike : 肺弊-快檔藹狼 凱 氦磐.
% pmat : 困俊輯 汲疙茄巴苞 鞍籃 膠內(nèi)絹等 雇膠媚 己盒狼 青紡.

  
  [d,numpoints] = size(X);
  M = length(weights);
  if (abs(sum(weights) - 1) < 0.01) == 0
     error('trainGMM: Weights should sum to 1!');
  end
 
  if ((size(mus) == [d,M]) & (size(sigmas) == [d,M])) == 0
     error('Incompatible dimensions!');
  end
 
  for n = 1:numpoints
    for m = 1:M
      pmat(m,n) = weights(m)*gauss_pdf(X(:,n),mus(:,m),diag(sigmas(:,m)));
    end
  end

  if M > 1
    loglike = log(sum(pmat));
  else
    loglike = log(pmat);
  end