wolffd@0: function bnet = hodbn_to_bnet(dbn, T)
wolffd@0: % DBN_TO_BNET Convert a DBN to a static network by unroll for T slices
wolffd@0: % bnet = dbn_to_bnet(dbn, T)
wolffd@0: ss = length(dbn.intra);
wolffd@0: [row,order] = size(dbn.equiv_class);
wolffd@0: eclass = [];
wolffd@0: for i = 1:min(order,T)
wolffd@0:      eclass = [eclass ; dbn.equiv_class(:,i)];
wolffd@0: end
wolffd@0: if T > order
wolffd@0:     eclass = [eclass ; repmat(dbn.equiv_class(:,order),T-order,1)]; 
wolffd@0: end
wolffd@0: 
wolffd@0: dnodes = unroll_set(dbn.dnodes_slice, ss, T);
wolffd@0: ns = repmat(dbn.node_sizes_slice(:), 1, T);
wolffd@0: dag = unroll_higher_order_topology(dbn.intra, dbn.inter, T, dbn.intra1);
wolffd@0: onodes = unroll_set(dbn.observed(:), ss, T);
wolffd@0: bnet = mk_bnet(dag, ns(:), 'discrete', dnodes(:), 'equiv_class', eclass(:), 'observed', onodes(:));
wolffd@0: bnet.CPD = dbn.CPD;
wolffd@0:  
wolffd@0: