1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
% make the structure of an embedded HMM with 2 rows and 3 columns
% 1------------>2
% |\ \ | \ \
% 3->4->5 6->7->8
n = 8;
dag = zeros(n);
dag(1,[2 3 4 5])=1;
dag(2,[6 7 8])=1;
for i=3:4
dag(i,i+1)=1;
end
for i=6:7
dag(i,i+1)=1;
end
ns = 2*ones(1,n);
bnet = mk_bnet(dag,ns);
for i=1:n
bnet.CPD{i}=tabular_CPD(bnet,i);
end
[jtree, root, cliques] = graph_to_jtree(moralize(bnet.dag), ones(1,n), {}, {});
%[jtree, root, cliques, B, w, elim_order, moral_edges, fill_in_edges] = dag_to_jtree(bnet);
|