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

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

/***************************************************************************
**?函數(shù)名稱：Backward
**?功能：后向算法估計(jì)參數(shù)
**?參數(shù)：phmm：指向HMM的指針
**???????T：觀察值序列的長(zhǎng)度
**???????O：觀察值序列
**???????beta：運(yùn)算中用到的臨時(shí)數(shù)組
**???????pprob：返回值，所要求的概率
**/
void?Backward（HMM?*phmm?int?T?int?*O?double?**beta?double?*pprob）
{
	int?????i?j;???/*?狀態(tài)?*/
	int?????t;??????/*?時(shí)間下標(biāo)?*/
	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];
}

/***************************************************************************
**?函數(shù)名稱：BackwardWithScale
**?功能：后向算法估計(jì)參數(shù)（帶比例因子修正）
**?參數(shù)：phmm：指向HMM的指針
**???????T：觀察值序列的長(zhǎng)度
**???????O：觀察值序列
**???????beta：運(yùn)算中用到的臨時(shí)數(shù)組
**???????scale：比例因子數(shù)組
**???????pprob：返回值，所要求的概率
**/
void?BackwardWithScale（HMM?*phmm?int?T?int?*O?double?**beta?
	double?*scale?double?*pprob）
{
	int?????i?j;???/*?狀態(tài)指示?*/
	int?????t;??????/*?時(shí)間下標(biāo)?*/
	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];
		}
	}
}

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????1919??2002-09-26?06:39??C語(yǔ)言的HMM算法庫(kù)\backward.cpp

?????文件???????4045??2002-09-26?05:15??C語(yǔ)言的HMM算法庫(kù)\baum.cpp

?????文件????????237??2005-04-06?20:37??C語(yǔ)言的HMM算法庫(kù)\forword.cpp

?????文件???????2280??2002-09-26?04:37??C語(yǔ)言的HMM算法庫(kù)\hmm.h

?????文件?????????66??2005-04-06?20:33??C語(yǔ)言的HMM算法庫(kù)\hmm.txt

?????文件????????434??2002-09-26?05:27??C語(yǔ)言的HMM算法庫(kù)\hmmrand.cpp

?????文件???????4081??2005-04-06?20:15??C語(yǔ)言的HMM算法庫(kù)\hmmutils.cpp

?????文件??????10836??2005-04-06?20:32??C語(yǔ)言的HMM算法庫(kù)\nrutil.cpp

?????文件???????1448??2002-09-26?03:32??C語(yǔ)言的HMM算法庫(kù)\nrutil.h

?????文件???????1190??2005-04-06?18:36??C語(yǔ)言的HMM算法庫(kù)\ReadMe.txt

?????文件????????290??2005-04-06?18:36??C語(yǔ)言的HMM算法庫(kù)\StdAfx.cpp

?????文件????????667??2005-04-06?18:36??C語(yǔ)言的HMM算法庫(kù)\StdAfx.h

?????文件???????3020??2002-09-26?06:07??C語(yǔ)言的HMM算法庫(kù)\viterbi.cpp

?????目錄??????????0??2005-04-07?21:47??C語(yǔ)言的HMM算法庫(kù)

-----------?---------??----------?-----??----

????????????????30513????????????????????14