数学建模常用模型---马尔可夫预测方法

浅夏110

马尔可夫预测的性质及运用
8 L: e- {" o. L& Y
       对事件的全面预测，不仅要能够指出事件发生的各种可能结果，而且还必须给出每一种结果出现的概率，说明被预测的事件在预测期内出现每一种结果的可能性程度。这就是关于事件发生的概率预测。
  t4 Y& B& g) I) b
, y6 v: N$ [! k) ^1 f( r       马尔可夫（Markov)预测法，就是一种关于事件发生的概率预测方法。它是根据事件的目前状况来预测其将来各个时刻（或时期)变动状况的一种预测方法。马尔可夫预测法是地理预测研究中重要的预测方法之一。
0 a# X  ]$ w- z% c3 Q0 e  V7 x; q
3 O; J! z7 d$ ]; p7 z基本概念
6 n0 G& a9 O+ c; p. n4 ?# H' H+ z- u
（一)状态、状态转移过程与马尔可夫过程

1.状态  在马尔可夫预测中，“状态”是一个重要的术语。所谓状态，就是指某一事件在某个时刻（或时期)出现的某种结果。一般而言，随着所研究的事件及其预测的目标不同，状态可以有不同的划分方式。譬如，在商品销售预测中，有“畅销”、“一般”、“滞销”等状态；在农业收成预测中，有“丰收”、“平收”、“欠收”等状态；在人口构成预测中，有“婴儿”、“儿童”、“少年”、“青年”、“中年”、“老年”等状态；等等。
$ ^( O! p* O% p3 o8 [+ C2 u5 Q2 X
( K8 E; {* g+ y2.状态转移过程  在事件的发展过程中，从一种状态转变为另一种状态，就称为状态转移。事件的发展，随着时间的变化而变化所作的状态转移，或者说状态转移与时间的关系，就称为状态转移过程，简称过程。



' B& q' R" w+ K0 R! K（二)状态转移概率与状态转移概率矩阵

* g  I8 ?2 f9 ^1.状态转移概率  在事件的发展变化过程中，从某一种状态出发，下一时刻转移到其它状态的可能性，称为状态转移概率。根据条件概率的定义，由状态Ei转为状态Ej的状态转移概率P（Ei→Ej)就是条件概率P（Ej/Ei)，即
4 }: t5 n& E0 G! u$ |
" _7 k' @) F9 b# Q                             
& i3 y; ]3 y) a3 B
- n" ~7 r$ e* Z. \2.状态转移概率矩阵  假定某一种被预测的事件有E1，E2，…，En，共n个可能的状态。记Pij为从状态Ei转为状态Ej的状态转移概率，作矩阵
# o0 F; |& C, S% x1 C: T
                        
" e) \+ X- |' b8 Q) l" E7 f
0 a" P0 }6 ~6 U! }则称P为状态转移概率矩阵。

如果被预测的某一事件目前处于状态Ei，那么在下一个时刻，它可能由状态Ei转向E1，E2，…Ei…En中的任一个状态。所以Pij满足条件：　

& u7 t5 i9 L6 @) S4 d                         

一般地，我们将满足上面条件的任何矩阵都称为随机矩阵，或概率矩阵。不难证明，如果P为概率矩阵，则对任何数m＞0，矩阵Pm都是概率矩阵。
/ {% }* {  g, W1 r3 V
2 }( P4 u$ j+ C* M如果P为概率矩阵，而且存在整数m＞0，使得概率矩阵Pm中诸元素皆非零，则称P为标准概率矩阵。可以证明，如果P为标准概率矩阵，则存在非零向量，而且满足，使得：
; r1 G( ~1 k! j
( y. C7 ^6 l% p4 |5 M                              ap=a

这样的向量α称为平衡向量，或终极向量。
: w$ A% m- i3 ?5 ~  ^
3.状态转移概率矩阵的计算  计算状态转移概率矩阵P，就是要求每个状态转移到其它任何一个状态的转移概率Pij（i，j=1，2，…，n)。为了求出每一个Pij，我们采用频率近似概率的思想来加以计算。
: B' ]; C& V6 q4 S3 Y/ |$ y
举例如下：
6 ^, f" p% @# ^( H2 |* i
考虑某地区农业收成变化的三个状态，即“丰收”、“平收”和“欠收”。记E1为“丰收”状态，E2为“平收”状态，E3为“欠收”状态。下表给出了该地区1950—1989年期间农业收成的情况以及状态变化：

1950	1951	1952	1953	1954	1955	1956	1957	1958	1959
1	2	3	4	5	6	7	8	9	10
E1	E1	E2	E3	E2	E1	E3	E2	E1	E2
1960	1961	1962	1963	1964	1965	1966	1967	1968	1969
11	12	13	14	15	16	17	18	19	20
E3	E1	E2	E3	E1	E2	E1	E3	E3	E1
1970	1971	1972	1973	1974	1975	1976	1977	1978	1979
21	22	23	24	25	26	27	28	29	30
E3	E3	E2	E1	E1	E3	E2	E2	E1	E2
1980	1981	1982	1983	1984	1985	1986	1987	1988	1989
31	32	33	34	35	36	37	38	39	40
E1	E3	E2	E1	E1	E2	E2	E3	E1	E2

以下，我们来计算该地区农业收成变化的状态转移概率矩阵。

从表2-18中可知，在15个从E1出发（转移出去)的状态转移中，有3个是从E1转移到E1的（即1→2，24→25，34→35)，有7个是从E1转移到E2的（即2→3，9→10，12→13，15→16，29→30，35→36，39→40)，有5个是从E1转移到E3的（即6→7，17→18，20→21，25→26，31→32)。

2 b: X  J$ P7 Y5 s( J# V        故
: }! W1 i, ~4 J" J            
按照上述同样的办法计算可以得到
6 q$ s# t7 ~$ a' e9 _% V/ H" N                 
! {6 ~5 s9 ]+ {: R2 a所以，该地区农业收成变化的状态转移概率矩阵为
                  
                                                                                  （1）
5 D7 v! E# J0 r0 v+ \% N3 k马尔可夫预测法

" `4 N3 k" Z( M$ ~' M& V1 I& i为了运用马尔可夫预测法对事件发展过程中状态出现的概率进行预测，还需要再介绍一个名词：状态概率πj（k)。πj（k)表示事件在初始（k=0)时状态为已知的条件下，经过k次状态转移后，第k个时刻（时期)处于状态Ej的概率。根据概率的性质，显然有：
  K; I7 J( ?8 S' t" d$ J; Y                     
                                                        （2）
* F1 a) p+ `' X4 h3 J' g
从初始状态开始，经过k次状态转移后到达状态Ej这一状态转移过程，可以看作是首先经过（k-1)次状态转移后到达状态Ei（i=1，2，…，n)，然后再由Ei经过一次状态转移到达状态Ej。根据马尔可夫过程的无后效性及Bayes条件概率公式，有
                     
1 d. ^6 h. V" n3 h                                                                                                (3)

若记行向量π(k)=[π1(k)，π2(k)，…，πn(k)]，则由(3)式可得逐次计算状态概率的递推公式：
1 n' ~2 Y( b  A0 Q) {                     
                                                                                            （4）
, Q, ?$ H- y) T& F, v( z- m, N
- v( T8 J+ h7 [4 B0 {$ i2 ?! g在（4）式中，π（0)=[π1（0)，π2（0)，…，πn（0)]为初始状态概率向量。

6 ]" [2 o+ M/ b. J由上述分析可知，如果某一事件在第0个时刻（或时期)的初始状态已知（即π（0)已知)，则利用递推公式（4)式，就可以求得它经过k次状态转移后，在第k个时刻（时期)处于各种可能的状态的概率（即π（k))，从而得到该事件在第k个时刻（时期)的状态概率预测。

3 h4 i+ u; C+ t8 o1 k; N. v在前例中，如果将1989年的农业收成状态记为π（0)=[0，1，0]（因为1989年处于“平收”状态)，则将状态转移概率矩阵（1)式及π（0)代入递推公式（4)式，就可以求得1990—2000年可能出现的各种状态的概率（见下表)。
; M7 f1 l8 `1 K  A' j1 @
) P6 E1 e: r$ Q

0 i8 X( J( E/ F4 y（二)终极状态概率预测

经过无穷多次状态转移后所得到的状态概率称为终极状态概率，或称平衡状态概率。如果记终极状态概率向量为π=[π1，π2，…，πn]，则
( _! z1 k0 P& `# Y+ P# z3 I0 J                     
$ Q- y2 v. |2 h$ ^                                                                                                               （5）
, n2 u) I7 \, `' v9 G! n
+ b( @* L7 K7 m3 A6 {- E即有：
# U6 D+ y$ T+ ]! P* \                     
                                           （6）
7 C$ e- \6 x3 m! }/ D, y/ t
将（6)式代入马尔可夫预测模型的递推公式（4)式得
/ Z& q. g% g8 J* k8 z- A/ J# t
2 F! s; t# b7 ~% {            π=πP
0 b( T0 K. p- d: o; z5 y; z2 v
2 c, ]; E% h/ K. }8 ]这样，就得到了终极状态概率应满足的条件：

7 [4 o4 y. t' _  i$ s（1)π=πP
5 _# s1 @8 V5 Z' d
  X# p8 @3 l) o! Q% F& o8 h# P（2)

9 a  O9 X3 s9 `! n8 t( L$ X% @（3）
7 i& X1 t0 e" n* Z' [+ k: K5 |
以上条件（2)与（3)是状态概率的要求，其中，条件（2)表示，在无穷多次状态转移后，事件必处在n个状态中的任意一个；条件（1)就是用来计算终极状态概率的公式。终极状态概率是用来预测马尔可夫过程在遥远的未来会出现什么趋势的重要信息。

* f; B$ i  @+ J& F- [在前例关于某地区农业收成状态概率的预测中，设终极状态的概率为π=[π1，π2，π3]，则
6 p# _$ y5 M$ t! t6 X4 w( `
                     
即：
                        
: G9 _5 c- N' ^, C) m                                                                                                          （7）

7 f; G3 e' W8 B) p# o) h4 X# J* D求解方程组（7)式得：π1=0.3653，π2=0.3525，π3=0.2799。这说明，该地区农业收成的变化，在无穷多次状态转移后，“丰收”和“平收”状态出现的概率都将大于“欠收”状态出现的概率。

7 D3 t9 U+ |$ r5 ?8 k9 a$ p9 {$ M, s在地理事件的预测中，被预测对象所经历的过程中各个阶段（或时点)的状态和状态之间的转移概率是最为关键的。马尔可夫预测的基本方法就是利用状态之间的转移概率矩阵预测事件发生的状态及其发展变化趋势。马尔可夫预测法的基本要求是状态转移概率矩阵必须具有一定的稳定性。因此，必须具有足够多的统计数据，才能保证预测的精度与准确性。

" x3 }- P( P# e8 p3 l
0 d' _+ j: w: ^+ K
# j4 Z  a# z/ J7 ?& @" C7 p