统计学习方法第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