www.pudn.com > gmm_utilities.zip > 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