hmm.h

HMM的演示源程序

#ifndef _INC_HMM
#define _INC_HMM

typedef struct {
	int N;
	int M;
	double	**A;
	double	**B;
	double	*pi;
} HMM;

void ReadHMM(FILE *fp, HMM *phmm);
void PrintHMM(FILE *fp, HMM *phmm);
void InitHMM(HMM *phmm, int N, int M);
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 T, int *O, int *q);
int GenInitalState(HMM *phmm);
int GenNextState(HMM *phmm, int q_t);
int GenSymbol(HMM *phmm, int q_t);

void ForwardWithScale(HMM *phmm, int T, int *O, double **alpha, double *scale, 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);
double *** AllocXi(int T, int N);
void FreeXi(double *** xi, int T, int N);
void ComputeGamma(HMM *phmm, int T, double **alpha, double **beta, double **gamma);
void ComputeXi(HMM* phmm, int T, int *O, double **alpha, double **beta, double ***xi);


#define MAX(x,y)        ((x) > (y) ? (x) : (y))
#define MIN(x,y)        ((x) < (y) ? (x) : (y))

#endif