annotate toolboxes/FullBNT-1.0.7/bnt/examples/dynamic/HHMM/Map/sample_from_map.m @ 0:cc4b1211e677 tip

initial commit to HG from Changeset: 646 (e263d8a21543) added further path and more save "camirversion.m"
author Daniel Wolff
date Fri, 19 Aug 2016 13:07:06 +0200
parents
children
rev   line source
Daniel@0 1 if 0
Daniel@0 2 % Generate some sample paths
Daniel@0 3
Daniel@0 4 bnet = mk_map_hhmm('p', 1);
Daniel@0 5 % assign numbers to the nodes in topological order
Daniel@0 6 U = 1; A = 2; C = 3; F = 4; O = 5;
Daniel@0 7
Daniel@0 8
Daniel@0 9 seed = 0;
Daniel@0 10 rand('state', seed);
Daniel@0 11 randn('state', seed);
Daniel@0 12
Daniel@0 13 % control policy = sweep right then left
Daniel@0 14 T = 10;
Daniel@0 15 ss = 5;
Daniel@0 16 ev = cell(ss, T);
Daniel@0 17 ev(U,:) = num2cell([R*ones(1,5) L*ones(1,5)]);
Daniel@0 18
Daniel@0 19 % fix initial conditions to be in left most state
Daniel@0 20 ev{A,1} = 1;
Daniel@0 21 ev{C,1} = 1;
Daniel@0 22 evidence = sample_dbn(bnet, 'length', T, 'evidence', ev)
Daniel@0 23
Daniel@0 24
Daniel@0 25 % Now do same but with noisy actuators
Daniel@0 26
Daniel@0 27 bnet = mk_map_hhmm('p', 0.8);
Daniel@0 28 evidence = sample_dbn(bnet, 'length', T, 'evidence', ev)
Daniel@0 29
Daniel@0 30 end
Daniel@0 31
Daniel@0 32 % Now do same but with 4 observations per slice
Daniel@0 33
Daniel@0 34 bnet = mk_map_hhmm('p', 0.8, 'obs_model', 'four');
Daniel@0 35 ss = bnet.nnodes_per_slice;
Daniel@0 36
Daniel@0 37 ev = cell(ss, T);
Daniel@0 38 ev(U,:) = num2cell([R*ones(1,5) L*ones(1,5)]);
Daniel@0 39 ev{A,1} = 1;
Daniel@0 40 ev{C,1} = 1;
Daniel@0 41 evidence = sample_dbn(bnet, 'length', T, 'evidence', ev)