隐马尔可夫模型后验状态概率
PSTATES=hmmdecode(seq,TRANS,emi)
[PSTATES,logpseq]=hmmdecode(…)
[PSTATES,logpseq,FORWARD,BACKWARD,S]=hmmdecode(…)
hmmdecode(…,'Symbols',符号)
PSTATES=hmmdecode(seq,TRANS,emi)
计算后验状态概率,PSTATES公司
,序列的序号
,从隐马尔可夫模型。后验状态概率是处于状态的条件概率K步进我,给定观察到的符号序列,符号
. 你用转移概率矩阵来指定模型,反式
,以及排放概率矩阵,电磁干扰
.反式(i,j)
是从状态转换的概率我
陈述J
.电磁干扰(k,seq)
是符号序号
从状态发出K
.
PSTATES公司
数组的长度与序号
模型中每个状态对应一行。那个(我,J)的th元素PSTATES公司
给出模型处于状态的概率我在J第四步,给定顺序序号
.
注意
函数HMM密码
在第一次发射之前,从步骤0处处于状态1的模型开始。HMM密码
计算中的概率PSTATES公司
基于模型从状态1开始的事实。
[PSTATES,logpseq]=hmmdecode(…)
退货logpseq公司
,序列概率的对数序号
,给定转移矩阵反式
和发射矩阵电磁干扰
.
[PSTATES,logpseq,FORWARD,BACKWARD,S]=hmmdecode(…)
返回按比例缩放的序列的前向概率和后向概率S
.
hmmdecode(…,'Symbols',符号)
指定发射的符号。符号
可以是数字数组、字符串数组或符号名称的单元格数组。默认符号是整数1.
通过N
,在哪里N
是可能排放的数量。
反式=[0.95,0.05;0.10,0.90]; emis=[1/6 1/6 1/6 1/6 1/6 1/6 1/6;1/10 1/10 1/10 1/10 1/10 1/2]; [seq,states]=hmmgenerate(100,trans,emis);pStates=hmmdecode(seq,trans,emis)[序号,状态]=hmmgenerate(100,trans,emis,…'符号',{'1','2','3','4','5','6'})pStates=hmmdecode(seq,trans,emis,…'符号“,{1”,“2”,“3”,“4”,“5”,“6});
[1] 杜宾,R.,S。艾迪,A。克罗和G。米奇森。生物序列分析. 英国剑桥:剑桥大学出版社,1998。