wolffd@0: function x = mixgauss_sample(mu, Sigma, labels)
wolffd@0: % MIXGAUSS_SAMPLE Sample from a mixture of Gaussians given known mixture labels
wolffd@0: % function x = mixgauss_sample(mu, Sigma, labels)
wolffd@0: 
wolffd@0: T = length(labels);
wolffd@0: [D Q] = size(mu);
wolffd@0: x = zeros(D,T);
wolffd@0: for q=1:Q
wolffd@0:   ndx = find(labels==q);
wolffd@0:   x(:,ndx) = gaussian_sample(mu(:,q)', Sigma(:,:,q), length(ndx))';
wolffd@0: end