www.gusucode.com > 用粒子滤波算法进行跟踪的matlab代码 > gmm_utilities/kernel_distance_KLD.m
function K = kernel_distance_KLD(g1, g2, N) s = kernel_samples(g1, N); w1 = kernel_evaluate(g1, s); w2 = kernel_evaluate(g2, s); %K = sum(log(w1) - log(w2)) / N; ii = find(w1 ~= 0 & w2 ~= 0); K = sum(log(w1(ii)) - log(w2(ii))) / length(ii);