Mercurial > hg > camir-aes2014
comparison toolboxes/FullBNT-1.0.7/bnt/examples/dynamic/HHMM/mk_hhmm_topo_F1.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 function [intra, inter, Qnodes, Fnodes, Onode] = mk_hhmm_topo_F1(D, all_Q_to_Qs, Ops) | |
2 % MK_HHMM_TOPO Make Hierarchical HMM topology assuming level 1 can finish | |
3 % function [intra, inter, Qnodes, Fnodes, Onode] = mk_hhmm_topo(D, all_Q_to_Qs, Ops, F1) | |
4 % | |
5 % D is the depth of the hierarchy | |
6 % If all_Q_to_Qs = 1, level i connects to all levels below, else just to i+1 [0] | |
7 % Ops are the Q parents of the observed node [Qnodes(end)] | |
8 % If F1=1, level 1 can finish (restart), else there is no F1->Q1 arc [0] | |
9 | |
10 Qnodes = 1:D; | |
11 | |
12 if nargin < 2, all_Q_to_Qs = 1; end | |
13 if nargin < 3, Ops = Qnodes(D); end | |
14 if nargin < 4, F1 = 0; end | |
15 | |
16 if F1 | |
17 Fnodes = 2*D:-1:D+1; % must number from bottom to top | |
18 Onode = 2*D+1; | |
19 ss = 2*D+1; | |
20 else | |
21 Fnodes = (2*D)-1:-1:D+1; | |
22 Onode = 2*D; | |
23 ss = 2*D; | |
24 end | |
25 | |
26 intra = zeros(ss); | |
27 intra(Ops, Onode) = 1; | |
28 for d=1:D-1 | |
29 if all_Q_to_Qs | |
30 intra(Qnodes(d), Qnodes(d+1:end)) = 1; | |
31 else | |
32 intra(Qnodes(d), Qnodes(d+1)) = 1; | |
33 end | |
34 end | |
35 for d=D:-1:3 | |
36 intra(Fnodes(d), Fnodes(d-1)) = 1; | |
37 end | |
38 if F1 | |
39 intra(Fnodes(2), Fnodes(1)) = 1; | |
40 end | |
41 if all_Q_to_Qs | |
42 for d=1:D | |
43 intra(Qnodes(d), Fnodes(d:end)) = 1; | |
44 end | |
45 else | |
46 for d=1:D | |
47 if d < D | |
48 intra(Qnodes(d), Fnodes([d d+1])) = 1; | |
49 else | |
50 intra(Qnodes(d), Fnodes(d)) = 1; | |
51 end | |
52 end | |
53 end | |
54 | |
55 inter = zeros(ss); | |
56 for d=1:D | |
57 inter(Qnodes(d), Qnodes(d)) = 1; | |
58 end | |
59 for d=1:D | |
60 if d==1 | |
61 inter(Fnodes(d), Qnodes(d)) = 1; | |
62 else | |
63 inter(Fnodes(d), Qnodes([d-1 d])) = 1; | |
64 end | |
65 end |