统计学习方法第10章隐马尔可夫模型后向算法例题10.2代码实践

统计学习方法第10章隐马尔可夫模型后向算法例题10.2代码实践

from numpy import *

Q=[1,2,3]
V=['红','白']
A=mat([[0.5,0.2,0.3],[0.3,0.5,0.2],[0.2,0.3,0.5]])
B=mat([[0.5,0.5],[0.4,0.6],[0.7,0.3]])
Pi=mat([0.2,0.4,0.4])
O=['红','白','红']

m=len(O)
backAlpha=mat(ones((m,1))).T
for i in range(m-1,0,-1):
    backAlpha=multiply(B[:,V.index(O[i])].T*A,backAlpha)

result=multiply(Pi,B[:,V.index(O[0])].T)*backAlpha.T
print(result)

结果：

[[ 0.133182]]

似乎结果跟前向算法有一点点不同

made by zcl at CUMT

I know I can because I have a heart that beats