gmm_samples_old.m

来自「高斯滤波器 matlab toolbox For GMMs and Gauss」· M 代码 · 共 17 行

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