gusucode.com > 用粒子滤波算法进行跟踪的matlab代码 > gmm_utilities/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