|
设B(k)是一个m×m矩阵,其元素为sbki,j=P(Xt=i,Xt+k=j),i,j∈ 十、K∈ Zand B(0)=diag(η,…,ηm)。然后b(k)=pXg=1λgB(k- g) Qg。证据首先考虑k=1的情况,p、 LetYkt={Xt+k-g:g=1,Pg6=k},thenbki,j=P(Xt=i,Xt+k=j)=XYktP(Xt=i,Xt+k=j|Ykt)P(Ykt)=XYktP(Xt+k=j | Xt=i,Ykt)P(Xt=i|Ykt)P(Ykt)=XYktpXg=1,g6=kλgqgXt+k-g、 jP(Xt=i | Ykt)P(Ykt)+XYktλkqki,jP(Xt=i | Ykt)P(Ykt)=pXg=1,g6=kλgXXt+k-gqgXt+k-g、 jP(Xt=i | Xt+k-g) P(Xt+k)-g) +λk^ηiqki,j=pXg=1,g6=kλgmXh=1bk-gi,hqgh,j+λk^ηiqki,jxg=1λgB(k)的第(i,j)个元素- g) 需要Qgas。C一类一般的MTDg模型let B(k)是一个m×m矩阵,其元素为bki,j=P(Xt=i,Xt+k=j),i,j=1,m、 k∈ Z、 其中B(0)=diag(η,…,ηm)。矩阵B(k)代表随机变量Xt的二元分布。然后,我们得到了b(k)=bk1,1··bk1,m-1^η-下午-1i=1bk1,i。。。。。。。。。。。。bkm-1,1··bkm-1米-1^ηm-1.-下午-1i=1bkm-1,i^η-下午-1i=1bki,1·^ηm-1.-下午-1i=1bki,m-12^ηm- 下午1点多-1i,j=1bki,j,其中,独立元素的总数为m- 每k为2m+1。p阶MTD模型的参数由向量λ=(λ,…,λp)和矩阵Qg组成,例如Qg=1T^η+~Qg,~Qg=qg1,1····qg1,m-1.-下午-1i=1qg1,i。。。。。。。。。。。。~qgm-1,1··qgm-1米-1.-下午-1i=1qgm-1.我-下午-1i=1ci~qgi,1··-下午-1i=1ciqgi,m-下午1点-1i,j=1ciqgi,j,式中,ηQg=0,g和ci=^ηi/^ηm。与这些定义一致,p阶马尔可夫链readP(Xt=i | Xt)的条件概率-1=我,Xt-p=ip)=ηi+pXg=1agig,i,其中agig,i≡ λgqgig,i.在此框架内,二元分布和矩阵qg满足以下矩阵方程组b(k)- ηTη=pXg=1B(k- g) Ag,Ag在哪里≡ λgQg。利用经验二元分布,可以对上述线性系统进行反演,以确定模型的参数。
|