Bayesian and HSMM

全面的HSMM 代码，包括详细例子。运行DEMO便可

%% set up parameters
% prior stuff for synthetic data
N = 4;
T = 500;
obs_dim = 2;
durparams = [12182530];
dur_hypparams = {85};
obs_hypparams = {zeros（obs_dim1）eye（obs_dim）obs_dim+30.1};

% inference
Nmax = 10; % weak limit approximation parameter
Niter = 100;
plot_results = 1;
plot_results_every = 10;

%% create truth model
for state=1:N
    truth_obs_distns{state} = observations.gaussian（obs_hypparams{:}）;
    truth_dur_distns{state} = durations.poisson（dur_hypparams{:}durparams（state））;
end

truthmodel = hsmm（Ttruth_obs_distnstruth_dur_distns）;

%% generate data from truth model
[data labels] = truthmodel.generate（）;

if plot_results
    fig = figure（）;
    truthmodel.plot（datafig）;
    title（‘TRUTH‘）;
end

%% create posterior model
for state=1:Nmax
    obs_distns{state} =

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\
     目录           0  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+durations\
     目录           0  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+durations\.svn\
     文件         392  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+durations\.svn\entries
     目录           0  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+durations\.svn\text-base\
     文件        1178  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+durations\.svn\text-base\poisson.m.svn-base
     文件        1178  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+durations\poisson.m
     目录           0  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+observations\
     目录           0  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+observations\.svn\
     文件         396  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+observations\.svn\entries
     目录           0  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+observations\.svn\text-base\
     文件        2641  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+observations\.svn\text-base\gaussian.m.svn-base
     文件        2643  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+observations\gaussian.m
     目录           0  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+util\
     目录           0  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+util\.svn\
     文件         845  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+util\.svn\entries
     目录           0  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+util\.svn\text-base\
     文件         186  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+util\.svn\text-base\logsumexp.m.svn-base
     文件         849  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+util\.svn\text-base\plot_gaussian_2D.m.svn-base
     文件         212  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+util\.svn\text-base\print_dot.m.svn-base
     文件         298  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+util\.svn\text-base\sample_discrete.m.svn-base
     文件         186  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+util\logsumexp.m
     文件         850  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+util\plot_gaussian_2D.m
     文件         212  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+util\print_dot.m
     文件         298  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\+util\sample_discrete.m
     文件         966  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\README.mkd
     文件        1288  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\demo.m
     文件        3341  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\hsmm.m
     文件         767  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\initial_state.m
     文件        7038  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\states.m
     文件        3135  2012-04-10 14:36  mattjj-matlab-hsmm-40a597a\transitions.m
............此处省略0个文件信息