马尔科夫预测方法

马尔可夫预测方法本节主要内容：几个基本概念状态；状态转移过程；马尔科夫过程；状态转移概率；状态转移概率矩阵。马尔可夫预测法状态转移概率；状态转移概率矩阵。对事件的全面预测，不仅要能够指出事件发生的各种可能结果，而且还必须给出每一种结果出现的概率。马尔可夫（Markov）预测法，就是一种预测事件发生的概率的方法。它是基于马尔可夫链，根据事件的目前状况预测其将来各个时刻（或时期）变动状况的一种预测方法。马尔可夫预测法是对地理事件进行预测的基本方法，它是地理预测中常用的重要方法之一。状态。指某一事件在某个时刻（或时期）出现的某种结果。状态转移过程。事件的发展，从一种状态转变为另一种状态，称为状态转移。马尔可夫过程。在事件的发展过程中，若每次状态的转移都仅与前一时刻的状态有关，而与过去的状态无关，或者说状态转移过程是无后效性的，则这样的状态转移过程就称为马尔可夫过程。几个基本概念状态转移概率。在事件的发展变化过程中，从某一种状态出发，下一时刻转移到其它状态的可能性，称为状态转移概率。由状态Ei转为状态Ej的状态转移概率是)(jiEEPijijjiPEEPEEP)/()(（3.7.1）状态转移概率矩阵。假定某一个事件的发展过程有n个可能的状态，即E1，E2，…，En。记为从状态Ei转变为状态Ej的状态转移概率，则矩阵几个基本概念)(jiEEP称为状态转移概率矩阵。概率矩阵。一般地，将满足条件（3.7.3）的任何矩阵都称为随机矩阵，或概率矩阵。nnnnnnPPPPPPPPPP212222111211（3.7.2）),,2,1(1),,2,1,(101niPnjiPnjijij（3.7.3）几个基本概念不难证明，如果P为概率矩阵，则对于任何整数m0，矩阵都是概率矩阵。mP标准概率矩阵、平衡向量。如果P为概率矩阵，而且存在整数m0，使得概率矩阵中诸元素皆非零，则称P为标准概率矩阵。可以证明，如果P为标准概率矩阵，则存在非零向量，而且满足，使得：（3.7.4）这样的向量α称为平衡向量，或终极向量。这就是说，标准概率矩阵一定存在平衡向量。mP],,,[21nxxxix10ixniix11几个基本概念P状态转移概率矩阵的计算。计算状态转移概率矩阵P，就是求从每个状态转移到其它任何一个状态的状态转移概率。为了求出每一个，一般采用频率近似概率的思想进行计算。)n,1,2,(i,jPij几个基本概念例题1：考虑某地区农业收成变化的三个状态，即“丰收”、“平收”和“欠收”。记E1为“丰收”状态，E2为“平收”状态，E3为“欠收”状态。表3.7.1给出了该地区1960～1999年期间农业收成的状态变化情况。试计算该地区农业收成变化的状态转移概率矩阵。表3.7.1某地区农业收成变化的状态转移情况年份1960196119621963196419651966196719681969序号状态年份序号状态年份序号状态年份序号状态1E1197011E3198021E3199031E12E1197112E1198122E3199132E33E2197213E2198223E2199233E24E3197314E3198324E1199334E15E2197415E1198425E1199435E16E1197516E2198526E3199536E27E3197617E1198627E2199637E28E2197718E3198728E2199738E39E1197819E3198829E1199839E110E2197920E1198930E2199940E2从表3.7.1中可以知道，在15个从E1出发（转移出去）的状态中，(1)有3个是从E1转移到E1的（即1→2，24→25，34→35）(2)有7个是从E1转移到E2的（即2→3，9→10，12→13，15→16，29→30，35→36，39→40）(3)有5个是从E1转移到E3的（即6→7，17→18，20→21，25→26，31→32）①计算：2000.0153)()(111111EEPEEPP4667.0157)()(122112EEPEEPP3333.0155)()(133113EEPEEPP所以同理可得：5385.0137)()(211221EEPEEPP1538.0132)()(222222EEPEEPP3077.0134)()(233223EEPEEPP3636.0114)()(311331EEPEEPP4545.0115)()(322332EEPEEPP1818.0112)()(333333EEPEEPP②结论：该地区农业收成变化的状态转移概率矩阵为1818.04545.03636.03077.01538.05385.03333.04667.02000.0P（3.6.5）状态概率及其计算状态概率：表示事件在初始（k＝0）状态为已知的条件下，经过k次状态转移后，在第k个时刻（时期）处于状态的概率。且：根据马尔可夫过程的无后效性及Bayes条件概率公式，有)(kπj1)(1njjk（3.7.6）),,2,1()1()(1njPkkniijjj（3.7.7）jE记行向量，则由（3.7.7）式可以得到逐次计算状态概率的递推公式：)](,),(),([)(21kkkknkPPkkPPP)0()1()()0()1()2()0()1(1（3.7.8）式中，为初始状态概率向量。)]0(,),0(),0([)0(21n第k个时刻（时期）的状态概率预测如果某一事件在第0个时刻（或时期）的初始状态已知，即已知，则利用递推公式(3.7.8)式，就可以求得它经过k次状态转移后，在第k个时刻（时期）处于各种可能的状态的概率，即，从而就得到该事件在第k个时刻（时期）的状态概率预测。)0()(k马尔可夫预测法例题2：将例题1中1999年的农业收成状态记为=[0,1,0]，将状态转移概率矩阵（3.7.5）式及代入递推公式（3.7.8）式，可求得2000——2010年可能出现的各种状态的概率（见表3.7.2）。)0(表3.7.2某地区1990—2000年农业收成状态概率预测值年份2000200120022003状态概率E10.5385E20.1528E30.3077E10.3024E20.414E30.2837E10.3867E20.3334E30.2799E10.3587E20.3589E30.2779年份2004200520062007状态概率E10.3677E20.3509E30.2799E10.3647E20.3532E30.2799E10.3656E20.3524E30.2799E10.3653E20.3526E30.2799年份200820092010状态概率E10.3653E20.3525E30.2799E10.3653E20.3525E30.2799E10.3653E20.3525E30.2799终极状态概率预测①定义：经过无穷多次状态转移后所得到的状态概率称为终极状态概率，即：②终极状态概率应满足的条件：)(lim)](lim,),(lim),(lim[21kkkkknkkkP),,2,1(10niinii11马尔可夫预测法③例题：在例1中，设终极状态的状态概率为则],,[3211818.04545.03636.03077.01538.05385.03333.04667.02000.0],,[],,[321321即：求解该方程组得：＝0.3653，＝0.3525，＝0.2799。这说明，该地区农业收成的变化过程，在无穷多次状态转移后，“丰收”和“平收”状态出现的概率都将大于“欠收”状态出现的概率。3213331232111818.03077.03333.04545.01538.04667.03636.05385.02000.0123在地理事件的预测中，被预测对象所经历的过程中各个阶段（或时点）的状态和状态之间的转移概率是最为关键的。马尔可夫预测的基本方法就是利用状态之间的转移概率矩阵预测事件发生的状态及其发展变化趋势。马尔可夫预测法的基本要求是状态转移概率矩阵必须具有一定的稳定性。因此，必须具有足够的统计数据，才能保证预测的精度与准确性。换句话说，马尔可夫预测模型必须建立在大量的统计数据的基础之上。这一点也是运用马尔可夫预测方法预测地理事件的一个最为基本的条件。