www.gusucode.com > 用粒子滤波算法进行跟踪的matlab代码 > gmm_utilities/kernel_distance_bhattacharyya.m

    function [Bc,Bd] = kernel_distance_bhattacharyya(g1, g2, N)

s = kernel_samples(g1, N); % use g1 as proposal 
w1 = kernel_evaluate(g1, s);
w2 = kernel_evaluate(g2, s);

Bc = sum(sqrt(w1.*w2)./w1) / N; % Bhattacharyya coefficient
Bd = -log(Bc);