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