[求助]viterbi译码源程序

zw123 · 发表于 2006-4-24 17:39

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

各位高手,谁提供一下实现(3,1,3)viterbi译码的源程序,不胜感激!!!

风花雪月 · 发表于 2006-5-15 10:26

/* ** File：Viterbi.cpp ** 功能：给定HMM和观察序列，求最可能的状态 */
//#include "StdAfx.h" #include <math.h> #include "hmm.h" #include "nrutil.h"
#define VITHUGE 100000000000.0
/************************************************************************** ** 函数名称：Viterbi ** 功能：Viterbi算法 ** 参数：phmm：HMM结构指针 ** T：观察值的个数 ** O：观察序列 ** delta，psi为中间变量 ** q：求得的最佳状态序列 ** pprob：概率 **/ void Viterbi(HMM *phmm, int T, int *O, double **delta, int **psi, int *q, double *pprob) { int i, j; /* 状态下标 */ int t; /* 时间下标 */ 
 int maxvalind; double maxval, val;
 /* 1. 初始化 */ for (i = 1; i <= phmm->N; i++) { delta[1] = phmm->pi * (phmm->B[O[1]]); psi[1] = 0; } 
 /* 2. 递归 */ for (t = 2; t <= T; t++) { for (j = 1; j <= phmm->N; j++) { maxval = 0.0; maxvalind = 1; for (i = 1; i <= phmm->N; i++) { val = delta[t-1]*(phmm->A[j]); if (val > maxval) { maxval = val; maxvalind = i; } } delta[t][j] = maxval*(phmm->B[j][O[t]]); psi[t][j] = maxvalind; 
 } }
 /* 3. 终止 */
 *pprob = 0.0; q[T] = 1; for (i = 1; i <= phmm->N; i++) { if (delta[T] > *pprob) { *pprob = delta[T]; q[T] = i; } }
 /* 4. Path (state sequence) backtracking */
 for (t = T - 1; t >= 1; t--) q[t] = psi[t+1][q[t+1]];
}
/************************************************************************** ** 函数名称：ViterbiLog ** 功能：Viterbi算法 ** 参数：phmm：HMM结构指针 ** T：观察值的个数 ** O：观察序列 ** delta，psi为中间变量 ** q：求得的最佳状态序列 ** pprob：概率 **/ void ViterbiLog(HMM *phmm, int T, int *O, double **delta, int **psi, int *q, double *pprob) { int i, j; /* 状态下标 */ int t; /* 时间下标 */ int maxvalind; double maxval, val; double **biot;
 /* 0. 预处理 */
 for (i = 1; i <= phmm->N; i++) phmm->pi = log(phmm->pi); for (i = 1; i <= phmm->N; i++) for (j = 1; j <= phmm->N; j++) { phmm->A[j] = log(phmm->A[j]); }
 biot = dmatrix(1, phmm->N, 1, T); for (i = 1; i <= phmm->N; i++) for (t = 1; t <= T; t++) { biot[t] = log(phmm->B[O[t]]); } /* 1. 初始化 */ for (i = 1; i <= phmm->N; i++) { delta[1] = phmm->pi + biot[1]; psi[1] = 0; } /* 2. 递归 */ for (t = 2; t <= T; t++) { for (j = 1; j <= phmm->N; j++) { maxval = -VITHUGE; maxvalind = 1; for (i = 1; i <= phmm->N; i++) { val = delta[t-1] + (phmm->A[j]); if (val > maxval) { maxval = val; maxvalind = i; } } delta[t][j] = maxval + biot[j][t]; psi[t][j] = maxvalind; } } /* 3. 终止 */ *pprob = -VITHUGE; q[T] = 1; for (i = 1; i <= phmm->N; i++) { if (delta[T] > *pprob) { *pprob = delta[T]; q[T] = i; } } /* 4. 回溯 */
 for (t = T - 1; t >= 1; t--) q[t] = psi[t+1][q[t+1]];
}

zw123 · 发表于 2006-5-15 16:43

太感谢了!! 能否再提供一下hmm.h和nrutil.h这两个头文件? 谢谢!!

账号		自动登录	找回密码
密码			我要加入

[其他相关] [求助]viterbi译码源程序

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（zw123）[求助]viterbi译码源程序

没有头文件