/*
**??????File:???backward.cpp
**??????功能：給定觀察值序列和HMM模型，利用前向后向算法
**????????????求取其概率
*/

//#include?“StdAfx.h“
#include?
#include?“hmm.h“

/***************************************************************************
**?函數名稱：Backward
**?功能：后向算法估計參數
**?參數：phmm：指向HMM的指針
**???????T：觀察值序列的長度
**???????O：觀察值序列
**???????beta：運算中用到的臨時數組
**???????pprob：返回值，所要求的概率
**/
void?Backward（HMM?*phmm?int?T?int?*O?double?**beta?double?*pprob）
{
	int?????i?j;???/*?狀態?*/
	int?????t;??????/*?時間下標?*/
	double?sum;
?
?
	/*?1.?初始化?*/
?
	for?（i?=?1;?i?<=?phmm->N;?i++）
		beta[T][i]?=?1.0;
?
	/*?2.?遞歸?*/
?
	for?（t?=?T?-?1;?t?>=?1;?t--）?
	{
		for?（i?=?1;?i?<=?phmm->N;?i++）?
		{
			sum?=?0.0;
			for?（j?=?1;?j?<=?phmm->N;?j++）
				sum?+=?phmm->A[i][j]?*?（phmm->B[j][O[t+1]]）*beta[t+1][j];
			beta[t][i]?=?sum;
?
		}
	}
?
	/*?3.?終止?*/
	*pprob?=?0.0;
	for?（i?=?1;?i?<=?phmm->N;?i++）
		*pprob?+=?beta[1][i];
}

/***************************************************************************
**?函數名稱：BackwardWithScale
**?功能：后向算法估計參數（帶比例因子修正）
**?參數：phmm：指向HMM的指針
**???????T：觀察值序列的長度
**???????O：觀察值序列
**???????beta：運算中用到的臨時數組
**???????scale：比例因子數組
**???????pprob：返回值，所要求的概率
**/
void?BackwardWithScale（HMM?*phmm?int?T?int?*O?double?**beta?
	double?*scale?double?*pprob）
{
	int?????i?j;???/*?狀態指示?*/
	int?????t;??????/*?時間下標?*/
	double?sum;
?
?
	/*?1.?初始化?*/
	for?（i?=?1;?i?<=?phmm->N;?i++）
		beta[T][i]?=?1.0/scale[T];?
?
	/*?2.?遞歸?*/
	for?（t?=?T?-?1;?t?>=?1;?t--）?
	{
		for?（i?=?1;?i?<=?phmm->N;?i++）?
		{
			sum?=?0.0;
			for?（j?=?1;?j?<=?phmm->N;?j++）
				sum?+=?phmm->A[i][j]?*?（phmm->B[j][O[t+1]]）*beta[t+1][j];
			beta[t][i]?=?sum/scale[t];
		}
	}
}

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????目錄??????????0??2012-02-16?15:12??HMM論文

?????目錄??????????0??2012-02-16?15:12??HMM論文\hmm2

?????文件???????1921??2005-04-02?19:11??HMM論文\hmm2\backward.cpp

?????文件???????4048??2005-04-02?19:10??HMM論文\hmm2\baum.cpp

?????目錄??????????0??2012-02-16?15:12??HMM論文\hmm2\Debug

?????文件???????3880??2012-02-13?15:05??HMM論文\hmm2\Debug\backward.obj

?????文件???????8030??2012-02-13?15:05??HMM論文\hmm2\Debug\baum.obj

?????文件???????4185??2012-02-13?15:05??HMM論文\hmm2\Debug\forward.obj

?????文件?????237607??2012-02-13?15:05??HMM論文\hmm2\Debug\hmm2.exe

?????文件?????234644??2012-02-13?15:05??HMM論文\hmm2\Debug\hmm2.ilk

?????文件?????218400??2012-02-13?15:05??HMM論文\hmm2\Debug\hmm2.pch

?????文件?????418816??2012-02-13?15:05??HMM論文\hmm2\Debug\hmm2.pdb

?????文件???????3435??2012-02-13?15:05??HMM論文\hmm2\Debug\hmmrand.obj

?????文件???????9296??2012-02-13?15:05??HMM論文\hmm2\Debug\hmmutils.obj

?????文件???????2049??2012-02-13?15:05??HMM論文\hmm2\Debug\main.obj

?????文件??????21056??2012-02-13?15:05??HMM論文\hmm2\Debug\nrutil.obj

?????文件??????41984??2012-02-15?13:42??HMM論文\hmm2\Debug\vc60.idb

?????文件??????53248??2012-02-13?15:05??HMM論文\hmm2\Debug\vc60.pdb

?????文件???????4997??2012-02-13?15:05??HMM論文\hmm2\Debug\viterbi.obj

?????文件???????2244??2005-04-02?19:12??HMM論文\hmm2\forward.cpp

?????文件???????2278??2005-04-02?19:11??HMM論文\hmm2\hmm.h

?????文件???????4825??2012-02-13?20:19??HMM論文\hmm2\hmm2.dsp

?????文件????????531??2012-02-13?20:19??HMM論文\hmm2\hmm2.dsw

?????文件??????50176??2012-02-15?20:00??HMM論文\hmm2\hmm2.ncb

?????文件??????50688??2012-02-15?20:00??HMM論文\hmm2\hmm2.opt

?????文件????????242??2012-02-15?13:42??HMM論文\hmm2\hmm2.plg

?????文件????????436??2005-04-02?19:11??HMM論文\hmm2\hmmrand.cpp

?????文件???????4081??2005-04-02?19:11??HMM論文\hmm2\hmmutils.cpp

?????文件?????????37??2005-04-02?19:15??HMM論文\hmm2\main.cpp

?????文件??????10823??2005-04-02?19:11??HMM論文\hmm2\nrutil.cpp

............此處省略11個文件信息