wolffd@0: function [obs, hidden] = dhmm_sample(initial_prob, transmat, obsmat, numex, len)
wolffd@0: % SAMPLE_DHMM Generate random sequences from a Hidden Markov Model with discrete outputs.
wolffd@0: %
wolffd@0: % [obs, hidden] = sample_dhmm(initial_prob, transmat, obsmat, numex, len)
wolffd@0: % Each row of obs is an observation sequence of length len.
wolffd@0: 
wolffd@0: hidden = mc_sample(initial_prob, transmat, len, numex);
wolffd@0: obs = multinomial_sample(hidden, obsmat);
wolffd@0: 
wolffd@0: