隐马尔可夫模型与异常检测的论文

/*
**      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论文hmm2ackward.cpp

     文件       4048  2005-04-02 19:10  HMM论文hmm2aum.cpp

     目录          0  2012-02-16 15:12  HMM论文hmm2Debug

     文件       3880  2012-02-13 15:05  HMM论文hmm2Debugackward.obj

     文件       8030  2012-02-13 15:05  HMM论文hmm2Debugaum.obj

     文件       4185  2012-02-13 15:05  HMM论文hmm2Debugforward.obj

     文件     237607  2012-02-13 15:05  HMM论文hmm2Debughmm2.exe

     文件     234644  2012-02-13 15:05  HMM论文hmm2Debughmm2.ilk

     文件     218400  2012-02-13 15:05  HMM论文hmm2Debughmm2.pch

     文件     418816  2012-02-13 15:05  HMM论文hmm2Debughmm2.pdb

     文件       3435  2012-02-13 15:05  HMM论文hmm2Debughmmrand.obj

     文件       9296  2012-02-13 15:05  HMM论文hmm2Debughmmutils.obj

     文件       2049  2012-02-13 15:05  HMM论文hmm2Debugmain.obj

     文件      21056  2012-02-13 15:05  HMM论文hmm2Debug
rutil.obj

     文件      41984  2012-02-15 13:42  HMM论文hmm2Debugvc60.idb

     文件      53248  2012-02-13 15:05  HMM论文hmm2Debugvc60.pdb

     文件       4997  2012-02-13 15:05  HMM论文hmm2Debugviterbi.obj

     文件       2244  2005-04-02 19:12  HMM论文hmm2forward.cpp

     文件       2278  2005-04-02 19:11  HMM论文hmm2hmm.h

     文件       4825  2012-02-13 20:19  HMM论文hmm2hmm2.dsp

     文件        531  2012-02-13 20:19  HMM论文hmm2hmm2.dsw

     文件      50176  2012-02-15 20:00  HMM论文hmm2hmm2.ncb

     文件      50688  2012-02-15 20:00  HMM论文hmm2hmm2.opt

     文件        242  2012-02-15 13:42  HMM论文hmm2hmm2.plg

     文件        436  2005-04-02 19:11  HMM论文hmm2hmmrand.cpp

     文件       4081  2005-04-02 19:11  HMM论文hmm2hmmutils.cpp

     文件         37  2005-04-02 19:15  HMM论文hmm2main.cpp

     文件      10823  2005-04-02 19:11  HMM论文hmm2
rutil.cpp

............此处省略11个文件信息