gusucode.com > 用粒子滤波算法进行跟踪的matlab代码 > gmm_utilities/gmm_counting_algorithm.m
function w = gmm_counting_algorithm(g, prob, N) %function w = gmm_counting_algorithm(g, prob, N) % % INPUTS: % g - Gaussian mixture % prob - probability mass to enclose [0, 1] % N - number of samples to count % % OUTPUT: % w - the (approximate) weight of g that marks the level-set enclosing % mass, prob. % % Tim Bailey 2006. if prob<=0 | prob>=1 error('Probability must be in range (0,1)') end s = gmm_samples(g, N); w = gmm_evaluate(g, s); [ws, i] = sort(w); M = N - ceil(prob*N); M(M==0) = 1; w = ws(M);