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)