sourcecodes/bnt-master/KPMstats/condgauss_sample.m


1
2
3
4
5
6
7
8
9
10
11

function x = mixgauss_sample(mu, Sigma, labels)
% MIXGAUSS_SAMPLE Sample from a mixture of Gaussians given known mixture labels
% function x = mixgauss_sample(mu, Sigma, labels)

T = length(labels);
[D Q] = size(mu);
x = zeros(D,T);
for q=1:Q
  ndx = find(labels==q);
  x(:,ndx) = gaussian_sample(mu(:,q)', Sigma(:,:,q), length(ndx))';
end