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