- [/font]%----------------------------------------------------------------------
- % starting value for P
- %----------------------------------------------------------------------
- T = n;
- sm=[S(1:T-1) S(2:T)]; % switching number, sm: (T-1)*2
- tm=zeros(sn,sn); % total number
- for im=1:sn
- for jm=1:sn
- tm(im,jm)=size(sm(sm(:,1)==im & sm(:,2)==jm,:),1); % total number of each transition from i state to j state
- end
- rr = pi_pri(im,:)+tm(im,:);
- P_draw = chi2rnd(2*rr);
- P(im,:)=P_draw./repmat(sum(P_draw,2),1,sn);
- % P(im,:)=draw_dirichlet(1,sn,pi_pri(im,:)+tm(im,:),0); % rm: prior for number of state transition (ij)
- end
前面设定的sn=2


雷达卡



京公网安备 11010802022788号







