第三章-马尔可夫链.doc

  • 格式:docx
  • 大小:71.43 KB
  • 文档页数:19

第三章马尔可夫链

、马尔可夫链的概念

马尔可夫过程是一类有重要应用意义的随机过程,它具有如下特征:随机过 程‘将来’所处的状态仅与‘现在’所处的状态有关,而与‘过去’曾处于什么 状态无关。

马尔可夫过程按其状态和时间参数是离散还是连续的可以分成三类

(1) 时间和状态都是离散的马尔可夫过程,称为马尔可夫链。

(2) 时间连续、状态离散的马尔可夫过程,称为连续时间的马尔可夫链

(3) 时间和状态都连续的马尔可夫过程。

本章介绍马尔可夫链 定义1设{Xn,n 0}为随机序列,其状态空间为I {i0,i1,i2,},如果对任意正

整数n及任意n+2个状态i0,i1,i2, ,in 1 I,有

P{Xn 1 in 1X0 io,X1 i1, ,Xn *}

P{Xn1 in1 Xn 叨 则称此随机序列{Xn,n 0}为马尔可夫链。

若将时刻n称为‘现在’,将时刻n+1称为‘将来’,而把0,1, 2,……,n-1 称为‘过去’。定义中的等式便可通俗解释为:在已知 {Xn,n 0} ‘现在’所

处的状态条件下,‘将来’所要达到的状态与‘过去’所经历的状态无关,这一 特性常称为马尔可夫的无后效性。

例1.一个n级数字传输系统,每一级的输入和输出信号只取 0或1两个值,每 一级的输出是下一级的输入;并假定当一级输入为0时,其输出为0和为1的概 率分别为p和1-p;当输入为1时,其输出为1和0的概率分别为p和1-p(见图)令Xn表示第n级输出,则{ Xn,n > 0}便为一个马尔可夫链。

例2•从1, 2,……,N数字中任取一个数,记为 X0;再从1, 2,……,X0数 字中任取一个数,记为X1;再从1, 2,……,X1中任取一个数,记为X2;依此 类推,在1, 2, , Xn-1中任取一个数,记为Xn。可以证明{ Xn,n > 0}为马 尔可夫链。

事实上,{ Xn,n >0}的状态空间为I={1 , 2,……,N},对任意正整数n,取n+1 个状态

i0,i1,i2, ,i n I , 由题意可知 故{ Xn,n >0}为马尔可夫链。

二、转移概率

由马尔可夫链的无后效性和乘法公式有

P{Xo io,Xi ii, ,Xn in}

P{Xn inXnl in l} ?P{ Xo io,Xi ii, ,Xn 1 in l}

P{Xn inXnl in l}P{ Xn 1 ini Xn 2 in 2} P{ Xi hX° io}P{X。 i。}

由此可见,马尔可夫链的统计特性完全由条件概率

P{Xn 1 iniXn in}所确定,所以如何确定这个条件概率就显得非常重要,我

们把这个条件概率称为一步转移概率。一般一步转移概率为 P{Xn 1 j Xn i}, 它表示系统在时刻n处于状态i的条件下,到时刻n+1转移到状态j的概率,记 为 Pj(n)。

定义2称条件概率Pij(n) P{Xn 1 j|Xn i}为马尔可夫链{Xn,n 0}在时刻n 的一步转移概率。

一般,转移概率Pij(n)不仅与状态i, j有关,而且与时刻n有关,但当它与时刻n 无关时,表示马尔可夫链具有平稳转移概率, 即与起点无关,此时我们称马尔可

夫链是齐次的。

定义3如果对任意的i,j I,马尔可夫链的转移概率Pj(n)与n无关,则称 {Xn,n 0}为齐次马尔可夫链,并记Pjj(n)卩耳。

下面我们只讨论齐次马尔可夫链。

设P为一步转移概率Pj所组成的矩阵,状态空间I {1,2, },称

Pl1 P12 P1n

P P21 P22 P2n P{Xn inXo io,Xi i1 , , X n 1 in i}?P{Xo io,Xi ini} 为马尔可夫链的一步转移概率矩阵。

转移概率矩阵具有下面性质

(1)Pij 0,i,j I (2) Pij 1 , i I

j i

称具有上面两条性质的矩阵为随机矩阵。

下面给出n步转移概率的概念

定义4称条件概率

P(n) P{Xm n jXm i},',j 5 0, n 1

为马尔可夫链的n步转移概率,并称

P(n) (P『)

为马尔可夫链的n步转移概率矩阵。其中p(n) 0, p(n) 1

j i

当n=1时,P⑴P,规定

切普曼--柯尔莫哥洛夫方程

定理1设{Xn, n 0}为马尔可夫链,则对任意正整数n,0 l n,和状态i, j I,

n步转移概率具有下列性质

(1) (n) Pij Pik) pkjn "(切普曼一柯尔莫哥洛夫方程)

k I

(2) (n) Pij Pi/ Pg Pkn 和

k1 1 kn 1 1

(3) p(n) pp(n 1)

(4) p(n) Pn

证明(1)利用全概率公式和马尔可夫性,有 (0)

Pij 0, i

1, i j,即P(0)为单位矩阵

j Pi(n) P{Xmn jXm i} P{ (Xml k), X m n jXm i}

k I

P{ (Xm1

k I k, X m n j)Xm i}

k I P{Xm 1 k, X m n jXm i}

P{Xm i,Xm 1 k, X m n j}

k I

P{Xm i}

P{Xm i}P{Xmi kXm i}P{Xmn jXm Mmi k}

k I P{Xm i}

P{Xm 1 kXm i}P{Xmn jXml k}

k I

(1) (n 1) (1) (n 1) Pik (m)Pkj (m 1) Pik Pkj

k I k I

(1)式是关于转移概率的一个重要结果,切普曼 -柯尔莫哥洛夫方程(简称为

C-K方程),直观上可以作如下解释:马尔可夫链 { Xn,n >0}在时刻m处于状态

i,经过n步,即在时刻m+n转移到状态j的过程可以视为它在时刻 m处于状态i ,

先经过I步,即在时刻m I遍历所有状态k (k 1,2,),然后再经过n I步,

即在时刻m+n转到状态j的转移过程(见下图)

[J] ” mH ,nr+n

C-K 方程的矩阵形式为 P(n) P(l)P(n l)

当 l 1 时,即为( 3) P(n) PP(n 1) ,再利用归纳法可证( 4) 在(1)中令丨1,k ,得p(n) piklpk;j1)这是一个递推公式,逐步递推可

kI

证( 2)

例 3 (随机游动)

设质点在线段上做随机游动。 (见图)。每隔一秒钟移动一步。当质点处于‘ O' 点时,必然要以概率 1 向右移动一步至‘ 1'点;当质点处于‘ 4'点时,下一步 必然以概率 1 向左移动一步至‘ 3'点;当质点处于其它点时,下一步便均分别

以概率 向左、向右或停留在原地不动。

令Xn表示n次移动后质点所处的位置。显然,{ Xn,n >0}为一齐次马尔可夫链, 其状态空间为 I={0 , 1 , 2, 3, 4}

试求{ Xn,n >0}的一步和二步转移概率矩阵:

解:按题意可知

同样可求得其它转移概率

于是便得一步转移概率矩阵 步转移概率矩阵便为

1. 一维分布

定义5设{Xn,n 0}为齐次马尔可夫链,

Pj(O) P{X。 j},j I

为{Xn,n 0}的初始分布,Pj(0)称为初始概率。将其写成向量形式为

PtO) (5(0门2(0), ,PN(0),),称为初始概率向量

f 0 1 0 0 0 1

r 0 1 0 0 0 1

1 1 1 1 1 1

---00 --------- 0 0

3 3 3 3 3 3

^111 1 1 1

0 ---------- 0 0 ---------- 0

3 3 3 3 3 3

…1 1 1 1 1 1

0 0——- 0 0 -----------

3 3 3 3 3 3

0 0 0 10 Q 0 0 1 0

\

i

齐次马尔可夫链的有限维分布 r 1 i 1 0

0 3 3 3

[ 5 2 [ a g g 9 g

i 2 3 2 1

9 9 9 9 9

0 1 2 5 1

9 9

1 1 ]

为I,称下列一组概率 定义6设{Xn,n 0}为齐次马尔可夫链,其状态空间为I,称下列一组概率

Pj(n) Pg j},j I

为{Xn,n 0}的绝对分布,Pj(n)称为绝对概率。将其写成向量形式为

pT(n) (Pi(n), P2(n), ,PN(n),),称为绝对概率向量。

定理2设{Xn, n 0}为齐次马尔可夫链,则对任意j I和n 1,绝对概率具有

下列性质

(1) Pj( n) Pi(0) Pi(n)

i I

(2) Pj( n) Pi(n 1)Pij

i I

(3) PT( n) PT (0)P(n)

(4) PT( n) PT(n 1)P

(3)( 4)式是(1)( 2)的矩阵形式。

2. n维分布

定理3设{Xn,n 0}为齐次马尔可夫链,对任意 m, ,in I和n 1,则马尔可 证明(1) Pj(n) P{Xn j} P{Xo i,Xn j}

i I

P{Xo i}P{Xn

i I Xo i} Pi (0)P(n)

i I

(2)Pj( n) P{Xn j} P{Xni i,Xn j}

i I

P{Xn i i}P{Xn

i I j Xn i i} Pi(n 1) Pj

i I