%隐马尔可夫模型维特比算法：求最优路径I%输入模型（A,B,pai），观测O=（o1,o2,...oT）%A状态转移概率矩阵%B观测概率矩阵%pai初始状态概率矩阵%T总时刻,N状态数%Q={1，2，3}状态集合，V={红，白}观测集合%课本p213，习题10.3A=[0.5 0.2 0.3;0.3 0.5 0.2;0.2 0.3 0.5];B=[0.5 0.5;0.4 0.6;0.7 0.3];pai=[0.2;0.4;0.4];T=4;N=3;O=[1;2;1;2];%1代表红球，2代表白球delta=zeros(T,N);pesai=zeros(T,N);I=zeros(T,1);%初始化 此时t=1for i=1:Ndelta(1,i)=pai(i)*B(i,O(1));end%递推for t=2:Tfor i=1:Ndelta(t,i)=max(delta(t-1,:).*A(:,i)')*B(i,O(t));[u,pesai(t,i)]=max(delta(t-1,:).*A(:,i)');%输出使得delta最大的结点iendend%终止P=max(delta(T,:));[r,I(T)]=max(delta(T,:));P%最优路径回溯for t=T-1:-1:1I(t)=pesai(t+1,I(t+1));endI %输出最优路径
输出结果如下
>> ma_viterbiP =0.0030I =3222
