Mercurial > hg > camir-aes2014
comparison toolboxes/FullBNT-1.0.7/bnt/examples/static/Brutti/Belief_IOhmm.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 % Sigmoid Belief IOHMM | |
2 % Here is the model | |
3 % | |
4 % X \ X \ | |
5 % | | | | | |
6 % Q-|->Q-|-> ... | |
7 % | / | / | |
8 % Y Y | |
9 % | |
10 clear all; | |
11 clc; | |
12 rand('state',0); randn('state',0); | |
13 X = 1; Q = 2; Y = 3; | |
14 % intra time-slice graph | |
15 intra=zeros(3); | |
16 intra(X,[Q Y])=1; | |
17 intra(Q,Y)=1; | |
18 % inter time-slice graph | |
19 inter=zeros(3); | |
20 inter(Q,Q)=1; | |
21 | |
22 ns = [1 3 1]; | |
23 dnodes = [2]; | |
24 eclass1 = [1 2 3]; | |
25 eclass2 = [1 4 3]; | |
26 bnet = mk_dbn(intra, inter, ns, dnodes, eclass1, eclass2); | |
27 bnet.CPD{1} = root_CPD(bnet, 1); | |
28 % ========================================================== | |
29 bnet.CPD{2} = softmax_CPD(bnet, 2); | |
30 bnet.CPD{4} = softmax_CPD(bnet, 5, 'discrete', [2]); | |
31 % ========================================================== | |
32 bnet.CPD{3} = gaussian_CPD(bnet, 3); | |
33 | |
34 % make some data | |
35 T=20; | |
36 cases = cell(3, T); | |
37 cases(1,:)=num2cell(round(rand(1,T)*2)+1); | |
38 %cases(2,:)=num2cell(round(rand(1,T))+1); | |
39 cases(3,:)=num2cell(rand(1,T)); | |
40 | |
41 engine = bk_inf_engine(bnet, 'exact', [1 2 3]); | |
42 | |
43 % log lik before learning | |
44 [engine, loglik] = enter_evidence(engine, cases); | |
45 | |
46 % do learning | |
47 ev=cell(1,1); | |
48 ev{1}=cases; | |
49 [bnet2, LL2] = learn_params_dbn_em(engine, ev, 3); |