www.gusucode.com > 用粒子滤波算法进行跟踪的matlab代码 > gmm_utilities/gmm_samples_old.m
function s = gmm_samples(g, N) % only works for equi-weight gaussians k = ceil(stratified_random(N)*size(g.x,2)); i=1; while i<=N ki = k(i); j = i+1; while j<=N & k(j)==ki j = j+1; end s(:,i:j-1) = gauss_samples(g.x(:,ki), g.P(:,:,ki), j-i); i = j; end