Mercurial > hg > camir-aes2014
comparison toolboxes/FullBNT-1.0.7/bnt/examples/dynamic/HHMM/Map/sample_from_map.m @ 0:e9a9cd732c1e tip
first hg version after svn
author | wolffd |
---|---|
date | Tue, 10 Feb 2015 15:05:51 +0000 |
parents | |
children |
comparison
equal
deleted
inserted
replaced
-1:000000000000 | 0:e9a9cd732c1e |
---|---|
1 if 0 | |
2 % Generate some sample paths | |
3 | |
4 bnet = mk_map_hhmm('p', 1); | |
5 % assign numbers to the nodes in topological order | |
6 U = 1; A = 2; C = 3; F = 4; O = 5; | |
7 | |
8 | |
9 seed = 0; | |
10 rand('state', seed); | |
11 randn('state', seed); | |
12 | |
13 % control policy = sweep right then left | |
14 T = 10; | |
15 ss = 5; | |
16 ev = cell(ss, T); | |
17 ev(U,:) = num2cell([R*ones(1,5) L*ones(1,5)]); | |
18 | |
19 % fix initial conditions to be in left most state | |
20 ev{A,1} = 1; | |
21 ev{C,1} = 1; | |
22 evidence = sample_dbn(bnet, 'length', T, 'evidence', ev) | |
23 | |
24 | |
25 % Now do same but with noisy actuators | |
26 | |
27 bnet = mk_map_hhmm('p', 0.8); | |
28 evidence = sample_dbn(bnet, 'length', T, 'evidence', ev) | |
29 | |
30 end | |
31 | |
32 % Now do same but with 4 observations per slice | |
33 | |
34 bnet = mk_map_hhmm('p', 0.8, 'obs_model', 'four'); | |
35 ss = bnet.nnodes_per_slice; | |
36 | |
37 ev = cell(ss, T); | |
38 ev(U,:) = num2cell([R*ones(1,5) L*ones(1,5)]); | |
39 ev{A,1} = 1; | |
40 ev{C,1} = 1; | |
41 evidence = sample_dbn(bnet, 'length', T, 'evidence', ev) |