function [y, dfce] = bayescls( X, model )
% BAYESCLS 啊快矯救 葛膽俊 措茄 海撈矯救 盒幅扁
% 
% 汲疙:
%  撈 竊薦綽 啊快矯救 葛膽俊 蝶扼 澇仿 氦磐 X俊 措茄 葷饒犬伏闌 拌魂竅絆 
%  澇仿 氦磐甸籃 葷饒 犬伏藹撈 啊厘 奴 努貳膠肺 盒幅 登絹柳促. 
%  
% 澇仿:
%  X [dim x num_data] 盒幅竅絆歷竅綽 澇仿 氦磐 
%
%  model.Mu [dim x num_classes] 啊快矯救 葛膽狼 乞閉 氦磐.
%  model.C [dim x dim x num_classes] 啊快矯救 葛膽狼 傍盒魂 青紡.
%  model.P [1 x num_classes] 葷傈 犬伏藹.
%
% 免仿:
%  y [1 x num_data] 扼駭 
%  dfce [num_classes x num_data] 犬伏藹 P(x|y)*P(y), 
%   咯扁輯,  x綽 盒幅等 氦磐撈絆 y綽 努貳膠 僥喊磊撈促. 
%

[dim,num_data]=size(X);
num_classes = size(model.Mu,2);

dfce=zeros(num_classes,num_data);

for i=1:num_classes,
  nconst = 1/((2*pi)^(dim/2) * sqrt(det(model.C(:,:,i))));
    
  dfce(i,:) = ...
     model.P(i)*nconst*exp(-0.5*mahalan(X,model.Mu(:,i),model.C(:,:,i)));
end

[tmp,y] = max(dfce);

return;