# 从数据估计隐马尔可夫过程

 In[1]:= Xdata = TemporalData[Automatic, {CompressedData[" 1:eJxNk9t1wzAMQwNwhg7QFTpKR8gC3f+vIi7k1ic5lig+AAj+fP98v/16vb7O /+P8zSPr/OzJY81oo7vR9PCsNrRnRE3O1s92mM05z5C1tbvbjF15GHX2e5KS tHR67SlNZ9Mn3RZM5tyahbjHG9i4xAymczhZpMckpJYrjdySTRAT1JH/ktJO YZSCkCUO5YbCL1AFViEMnAbOISJQgQWRjYx/+ClEnUcXRoZ0ssbNhHBguEKF k5CviIVw8BggyKA1bEoZ8W7fueow+d4FEAJYIo86UOKGqSLgq6y6ja5psrhO CVSnN+bDAu3sDvUVO15FJnijQm4Px7pG5B56xdKDcLj6egkV7atKd3M7U1Qz 9u15uNnFpfo708u9KgY0KdjLc0XBgyFZY/aTofTKz91W3+cjmqLmVf16l3Vf ptRFV2rfu43UJVMFT+wX64AKKw== "], {{0, 20, 1}}, 50, {"Discrete", 50}, { "Discrete", 1}, 1, {ValueDimensions -> 1, MetaInformation -> {}}}, False, 10.];
 In[2]:= Xhmm = EstimatedProcess[data, HiddenMarkovProcess[2, 3]]
 Out[2]=

 In[3]:= XLogLikelihood[hmm, data]
 Out[3]=

 In[5]:= XHistogram[data, Automatic, "PDF", PlotTheme -> "Monochrome"]
 Out[5]=

 In[6]:= XeprocBW = EstimatedProcess[data, HiddenMarkovProcess[2, "Gaussian"]]
 Out[6]=
 In[7]:= XeprocVT = EstimatedProcess[data, HiddenMarkovProcess[2, "Gaussian"], ProcessEstimator -> "ViterbiTraining"]
 Out[7]=

BaumWelch 估计过程中数据有更高的对数似然值.

 In[8]:= XLogLikelihood[#, data] & /@ {eprocBW, eprocVT}
 Out[8]=

## Mathematica

Questions? Comments? Contact a Wolfram expert »