www.pudn.com > gmm_utilities.zip > gmm_distance_bayes.m

function D = gmm_distance_bayes(g1, g2) 
 
% Normalising constant from Gmm multiplication 
 
D = 0; 
for i=1:size(g1.x,2) 
    for j=1:size(g2.x,2) 
        wij = gauss_likelihood(g1.x(:,i)-g2.x(:,j), g1.P(:,:,i)+g2.P(:,:,j)); 
        D = D + g1.w(i) * g2.w(j) * wij; 
    end 
end