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