/**************************************************************************      File:   hmm.h**      功能：提供HMM的數據結構和定義*/#include "StdAfx.h"#include <stdlib.h>#include <stdio.h>#include <math.h>typedef struct {	int N;		/* 狀態的個數；  Q={1,2,...,N} */	int M; 		/* 觀察值得個數; V={1,2,...,M}*/	double	**A;	/* A[1..N][1..N]. a[i][j] 是從t時刻的狀態i轉移到t+1時刻
					的狀態j的轉移概率*/	double	**B;	/* B[1..N][1..M]. b[j][k] 是在狀態j觀察到觀察值k的概率 */	double	*pi;	/* pi[1..N] pi[i] 是進入Markov鏈的狀態分布 */} HMM;

//關于HMM的文件存取、初始化、拷貝函數void ReadHMM(FILE *fp, HMM *phmm);void PrintHMM(FILE *fp, HMM *phmm);void InitHMM(HMM *phmm, int N, int M, int seed);void CopyHMM(HMM *phmm1, HMM *phmm2);void FreeHMM(HMM *phmm);
//關于序列的存取，狀態的計算等等void ReadSequence(FILE *fp, int *pT, int **pO);void PrintSequence(FILE *fp, int T, int *O);void GenSequenceArray(HMM *phmm, int seed, int T, int *O, int *q);int GenInitalState(HMM *phmm);int GenNextState(HMM *phmm, int q_t);int GenSymbol(HMM *phmm, int q_t);//前向后向算法、BaumWelch算法void Forward(HMM *phmm, int T, int *O, double **alpha, double *pprob);void ForwardWithScale(HMM *phmm, int T, int *O, double **alpha,        double *scale, double *pprob);void Backward(HMM *phmm, int T, int *O, double **beta, double *pprob);void BackwardWithScale(HMM *phmm, int T, int *O, double **beta,        double *scale, double *pprob);void BaumWelch(HMM *phmm, int T, int *O, double **alpha, double **beta,        double **gamma, int *niter, 	double *plogprobinit, double *plogprobfinal);
//Viterbi算法                                  double *** AllocXi(int T, int N);void FreeXi(double *** xi, int T, int N);void ComputeGamma(HMM *phmm, int T,double **gamma,double ***xi);void ComputeXi(HMM* phmm, int T, int *O, double **alpha, double **beta,        double ***xi);void Viterbi(HMM *phmm, int T, int *O, double **delta, int **psi,        int *q, double *pprob);void ViterbiLog(HMM *phmm, int T, int *O, double **delta, int **psi,        int *q, double *pprob);//隨機數相關函數int hmmgetseed(void);void hmmsetseed(int seed);double hmmgetrand(void); #define MAX(x,y)        ((x) > (y) ? (x) : (y))#define MIN(x,y)        ((x) < (y) ? (x) : (y))