主要内容

HMM密码

隐马尔可夫模型后验状态概率

语法

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。

R2006a之前引入