第四章马氏链(1)
4.马尔可夫链1

qr p
i 1
p,
pij
r, q,
0,
i i 1
j i 1 j i j i 1 其它
其一步转移概率矩阵为
i 2 i 1 i i 1... ... ...
... ... ... ... ... ... ... ... ... ... q r p 0 0 ... 0 ... i 1 ... 0 q r p 0 ... 0 ... i ... 0 0 q r p ... 0 ... i 1 ... ... ... ... ... ... ... ... ...
(时间离散、状态连续的马尔可夫过程,通常用泛函中 二元函数的范数进行研究)
例1 独立过程 X t,t T 是马尔可夫过程
证 设Xt,t T是独立过程,对于t1 t2 ... tn T,
X t1, X t2 ,..., X tn 相互独立,因此
P X tn xn X t1 x1, X t2 x2,..., X tn1 xn1 =P X tn xn =P X tn xn X tn1 xn1
率r原地不动。若以 X n 表示时刻 n 时质点的位置,
则X n,n 0,1, 2,...是一个随机过程。而且当
X n i 时,X n+1,X n+2,...,X n+k,...等 n时刻后质点所处的状态,只与X n i 有关,而与
质点在n以前是如何到达i的完全无关。所以它是一
个齐次马尔可夫链,其状态空间为I= ,-2,-1,0,1,2,
则称 Xn,nT 为马尔可夫链.
定义4.1 设有随机过程Xn ,n T,若
对于任意的整数n T和任意的 i0,i1,...,in1 I 条件概率满足
随机过程习题集-第四章马尔可夫过程

1第四章 马尔可夫过程内容提要1. 马尔可夫过程的概念 (1)马尔可夫过程给定随机过程{}(),X t t T ∈,如果对122,∀≥∀<<<∈n n t t t T ,有11221111{()|(),(),,()}{()|()}n n n n n n n n P X t x X t x X t x X t x P X t x X t x ----<====<=则称{}(),X t t T ∈为马尔可夫过程。
称(){}:,==∈E x X t x t T 为状态空间。
参数集和状态空间都是离散的马尔可夫过程称为离散参数马氏链. 参数连续、状态空间离散的马尔可夫过程称为连续参数马氏链. (2)k 步转移概率设{}(),0,1,2,=X n n 为离散参数马氏链,称()(),(,){|},0,1=+==≥≥i j p n k P X n k j X n i n k为{}(),0,1,2,=X n n 在时刻n 的k 步转移概率,称(),(,)((,)),P =∈i j n k p n k i j E为{}(),0,1,2,=X n n 在时刻n 的k 步转移概率矩阵. 特别地,当1k =时,在时刻n 的一步转移概率和一步转移概率矩阵分别简记为()ij p n 和()n P . (3)初始分布、绝对分布称((0)),,==∈i p P X i i E 为离散参数马氏链{}(),0,1,2,=X n n 的初始分布,记为0P ,称()(){},,==∈j p n P X n j j E 为马尔可夫链{}0n X n ≥的绝对分布,记为P n . (4)离散参数齐次马氏链设{}(),0,1,2,=X n n 是一离散参数马氏链,如果其一步转移概率()ij p n 恒与起始时刻n 无关,记为ij p ,则称{}(),0,1,2,=X n n 为离散参数齐次马氏链。
若{}(),0,1,2,=X n n2是离散参数齐次马氏链,则其k 步转移概率记为(),i j p k ,一步转移概率矩阵和k 转移概率矩阵分别记为P 和().P k(5) 离散参数齐次马氏链的遍历性离散参数齐次马氏链{X (n ) ,n=0,1,2… },若对一切状态i ,j ,存在与i 无关的极限()()lim 0,ij j n p n i j E →+∞=π>∈则称此马氏链具有遍历性.0,1j j j Ej E ππ∈>∈=∑若且则称{},j j E π∈为离散参数齐次马氏链{X (n ) ,n=0,1,2… }的极限分布,或称为最终分布,记为{},j j E ∏=∈π(6)离散参数齐次马氏链的平稳分布离散参数齐次马氏链{X (n ) ,n=0,1,2… },若存在{v j , j ∈E } 满足条件:1)0,2)13)j jj Ej i iji Ev j E vv v p ∈∈≥∈==∑∑则称此马氏链是平稳的,称 { v j , j ∈E } 为此马氏链的平稳分布。
第四章 马尔可夫链

一步转移概率
定义4.2 称条件概率pij(n)= P{Xn+1=j|Xn=i} 为马尔 可夫链{Xn,nT }在时刻n的一步转移概率, 简称转移概率,其中i,jI。 定义4.3 若对任意的i,jI,马尔可夫链{Xn,nT } 的转移概率pij(n)与n无关,则称马尔可夫链是 齐次的,并记pij(n)为pij。
如果d>1,就称i为周期的, 如果d=1,就称i为非周期的。
引理4.1 如果i的周期为d,则存在正整数M,对一切 ( nd ) n≥M ,有 p ii 0。
例4.6
设马尔可夫链的状态空间I={1,2,,9},转移概率如 下图所示。
1 1
8
9
1
1 3
2
1
7
1
1 6
1
3
1
5
2 3
4
1
从状态1出发再返回状态1的可能步数为T={4,6,8,10, },T的最大公约数为2,从而状态1的周期为2。
P{ X n j | X n1 i}P{ X n1 i} pi (n 1) pij
iI iI iI
(3)(4)为(1)(2)的矩阵表示。
定理4.3 设{Xn,nT }为马尔可夫链,则对任意 整数i1, i2,,inI和n1 ,有性质
P{ X1 i1 ,, X n in } pi pii1 pi1i2 pin1in
证明: (1) p j (n) P{ X n j} P{ X 0 i , X n j}
P{ X n j | X 0 i}P{ X 0 i}
iI ( p i p ijn ) iI iI
(2) p j (n) P{ X n j} P{ X n1 i , X n j}
马尔可夫链

例7 设马氏链{Xn}的状态空间为 I={1, 2, 3, 4, 5}, 转移概率矩阵为
1 2
1
2
0 0
0
1 2
1 2
0
0
0
P 0 0 1 0 0
3 / 16 . 1/ 4
于是: (1) P{X0 0, X2 1}
P{ X0 0}P{ X2 1 | X0 0} 1 5 5 ;
3 16 48
2020年5月21日星期四
(2)P{X2 1}
P{X0 0}P{X2 1 | X0 0} P{X0 1}P{X2 1 | X0 1}
显然有
p(n) 11
p(n) 21
P(n)
p(n j1
)
L
p(n) 12
p(n) 22
p(n) 1j
L
p(n) 2j
L
p(n) j2
p(n) jj
L
LL
L
(1)
0
p(n) ij
1
(2)
p(n) ij
1,
i
1,
2,L
j
2020年5月21日星期四
切普曼-柯尔莫哥洛夫方程(C-K方程): 对任意的m,n≥0,有
的矩阵
p11 p21
P
L
pj1 L
p12 L p22 L LL pj2 L LL
p1 j L
p2 j L
L
L
p jj L
L L
称为一步转移概率矩阵. 显然有
(1) 0 pij 1
(2)
pij 1, i 1, 2,L
j
2020年5月21日星期四
3、马尔可夫链举例
第四章马尔可夫链

i1
Pi , j 0
j . i 1 ,i-1 , i 1
1 0 0 0 0 . .
q
0
p
0
0
.
.
0 q 0 p 0 . .
P
0
0
q
0
p
.
.
0 0 0 q 0 . . . . . . . . .
.
例题:带2个吸收壁的随机游动
质点在数轴上移动,规律同上例。随机游动的状态 空间I={0,1,2…a}, 其中0和a为吸收态 。求一步转移p12 p1n Pp21 p22 p2n
称为系统状态的一步转移概率矩阵,它具有 如下性质:
1. pij 0, i, jI
2. pij 1, i, jI jI
满足上述两个性质的矩阵成为随机矩阵
.
定义4.4
称条件概率 p i(n ) j P { X m n j|X m i}i,j I,m 0 ,n 1 为马尔可夫链{Xn,n∈T}的n步转移概率,并称
0 1 1
.
马尔可夫链的状态分类
周期、非周期 常返、非常返
其中,常返分为正常返、零常返 非周期的正常返称为遍历状态
到达和互通
.
设马尔可夫链的状态空间I={1,2,3,4,5,6,7,8,9}, 状态转移图如下图
8
9
2
7
1
3
6
5
4
观察状态1
.
定义4.6 如集合{n: n≥1,pii(n)>0}非空,则称该集合的 最大公约数d=d(i)=G.C.D{n:pii(n)>0}为状态i 的周期。如d>1就称i为周期的,如d=1就称i 为非周期的。
.
马氏链

模型假设
钢琴每周需求量服从波松分布,平均每周 架 钢琴每周需求量服从波松分布,平均每周1架. 存贮策略:当周末库存量为零时,订购 架 存贮策略:当周末库存量为零时,订购3架,周初 到货;否则,不订购. 到货;否则,不订购 以每周初的库存量作为状态变量, 以每周初的库存量作为状态变量,状态转移具有 无后效性. 无后效性 在稳态情况下计算失去销售机会的概率和每周的 平均销售量, 作为该存贮策略的评价指标. 平均销售量 作为该存贮策略的评价指标
11.1 健康与疾病
通过有实际背景的例子介绍马氏链的基本概念和性质. 通过有实际背景的例子介绍马氏链的基本概念和性质 人的健康状态随着时间的推移会随机地发生转变. 人的健康状态随着时间的推移会随机地发生转变 保险公司要对投保人未来的健康状态作出估计, 保险公司要对投保人未来的健康状态作出估计 以制 订保险金和理赔金的数额 . 人的健康状况分为健康和疾病两种状态, 例1. 人的健康状况分为健康和疾病两种状态,设对特 定年龄段的人,今年健康、 定年龄段的人,今年健康、明年保持健康状态的概率 而今年患病、明年转为健康状态的概率为0.7. 为0.8, 而今年患病、明年转为健康状态的概率为 若某人投保时健康, 年后他仍处于健康状态的概率. 若某人投保时健康 问10年后他仍处于健康状态的概率 年后他仍处于健康状态的概率
模型建立
Dn P 0
Dn~第n周需求量,均值为 的波松分布 周需求量, 第 周需求量 均值为1的波松分布
P( Dn = k ) = e / k! (k = 0,1,2L)
1 0.368 2 0.184 3 0.061 >3
−1
0.368
0.019
Sn~第n周初库存量 状态变量 ) Sn ∈{1,2,3} 状态转移阵 周初库存量(状态变量 第 周初库存量 p11 p12 p13 Sn − Dn , Dn < Sn 状态转 S = n+1 P = p21 p22 p23 Dn ≥ Sn 移规律 3,
马氏链简介

显然,cij 0
n
n
cij (ai1b1 j ai2b2 j ainbnj ), (i 1,2,, n)
j 1
j 1
n
n
cij (ai1b1 j ai2b2 j ainbnj ), (i 1,2,, n)
j 1
j 1
n
用随机变量 X n 表示第 n 个月的经营状况
Xn 1 Xn 2
表示销路好; 表示销路坏;
n 0,1,2
X n 称为这个经营系统的状态。
用ai n表示第 n 月处于状态 i 的概率, i 1,2,
即
ai n PXn i ai n 称为状态概率。
pij 表示已知这月处于状态 i ,下月处于状态 j 的概率,i, j 1,2, 即
每次传播消息的失真率为 p, 0 p 1,
即 ai 将消息传给 ai1, 时,传错的概率为 p
这样经过长时间传播第n个人得知消息时,消息 的真实程度如何?
第n个人知道消息可能是真,也可能是假, 有两种状态,记为
Xn 1 Xn 2
表示消息假; 表示消息真;
n 0,1,2
用ai n表示第 n 个人处于状态 i 的概率, i 1,2,
n01 2 3
4
ห้องสมุดไป่ตู้
1 0.5 0.45 0.445 0.4445 ?
0 0.5 0.55 0.555 0.5555 ?
a1(n)
5 10
5 102
a2 (n)
4 9
5 10n
5 10
1 (
1 10n1
4马氏链

可见,{Xn,n=0,1,2,…}是一个马氏链。
Pij ( m , m + n)∆ P { X n + m = j | X m = i}
称为马氏链在时刻 m 系统处于状态 i 的条件下,在时 刻 m+n 转移到状态 j 的转移概率。
2.转移概率的性质
(1) Pij≥0;
(2)
∑ P (m , m + n) = 1, i = 0,1, 2,⋯
对任意的 n 及 i 0 , i 1 , ⋯ , i n , i n + 1 ∈ x ,
P {X n +1 = i n +1 X 0 = i 0 , X 1 = i1 , ⋯ , X n = i n } 0 i n+1 > i n =1 = P{X n+1 = i n+1 | X n = i n } i n +1 ≤ i n in
马尔可夫链及其概率分布 引言
直观上,过程(或系统)在时刻t0所处的状态为已 知的条件下,过程在时刻t>t0所处状态的条件分布与过 程在时刻t0之前所处的状态无关。 用分布函数表达此性质,设随机过程{X(t),t∈T}, 状态空间为χ,若对于t 的任意n个值t1<t2<…<tn,n≥3, 有
P {X ( t n ) ≤ xn X ( t1 ) = x1 , X ( t 2 ) = x 2 , ⋯ , X ( t n−1 ) = x n−1 }
条件下, 即在 X ( t i ) = x i , i = 1,2,⋯ , n − 1条件下,X ( t n )的条件分 布函数等于在条件 X ( t n−1 ) = x n−1下X ( t n )的条件分布函 数。
则称过程{X(t),t∈T}具有马尔可夫性,或称 {X(t),t∈T}为马尔可夫过程。
第6讲 第4章马尔科夫链(1)

P X n j P X 0 k P X n j | X 0 k k P X 0 k pkjn j I k
若记: pj P X 0 j, n 1,2
P0 p1, p2 ,
202022128是具有三个状态的齐次马氏链初始分20202212920202213016162411202022131161616162020221321221202022133于多少日为雨天的概率各等又已知的一步转移概率矩阵试写出马氏链天状态表示第表示雨天状态为逆事件任一天晴或雨是互晴天转雨天的概率为雨天转晴天的概率为设任意相继的两天中例48先求出转移概率矩阵121323202022134181118121213231204167040050599503997060030105995
Xn 表示比赛 n 局后,甲的累积得分,n 1, 2, .
( 1) 写出状态空间; ( 2) 求2步转移概率; ( 3) 问在甲获得1分的情况下, 最多再赛2局可以 结束的概率.
2020/1/22
20
解 (1) S 2, 1, 0, 1, 2
(2)
-2 -1 0 1 2
2 1 0 0 0 0
1 q r p 0 0
P
0 0
1
0
q 0
r q
p r
0
p
2 0 0 0 0 1
2020/1/22
21
-2 -1 0
12
2 1
0
0
0
0
P 2 1
0
q rp q2
r2 pq 2rq
2 pr r2 2 pq
p2 2 pr
随机数学 第5讲 第4章马尔科夫链(1)

等时刻发生游动 .
分析:状态空间I={1,2,3,4,5} 这是一个马氏链, 且是齐次的.
1
注意: 以后,我们所讨论的马氏链都是齐次的,并简 称马氏链。
我们把各状态之间的转移概率排列起来,形成矩阵, 为转移概率矩阵。
一步转移概率阵:
P = P(1) = ( pij ) =
n步转移概率矩阵:
a1 a2 a 1 ⎡ p ( n) p ( n) 11 12 P ( n) = ( pij ( n)) = a ⎢ p22 ( n ) 2 ⎢ p21 ( n )
ai
⎢ ⎢ ⎢ pi 1 ( n ) ⎢ ⎢ ⎣ p21 ⎢ ⎢ ⎢ a i ⎢ pi 1 ⎢ ⎣
a1
a2
p12 p22 pi 2
p1 j p1 j pij
aj
⎤ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦
0⎤ 0⎥ ⎥ ⎥ ⎥ ⎦
p1 j ( n ) p1 j ( n ) pij ( n )
= ∑ P ( X ( 0 ) = k , X ( n1 ) = i1 , X ( n2 ) = i2 ,
k
定理4.3 齐次马尔科夫链的有限维分布满足:
(1)P ( X ( 0 ) = i0 , X ( n1 ) = i1 , X ( n2 ) = i2 ,
(2)P ( X ( n1 ) = i1 , X ( n2 ) = i2 ,
m m−1
马尔可夫链的无后效性,用条件分布律表 示,就有了马尔可夫链的如下定义:
二、 转移概率、转移概率矩阵
定义4.2
称 P{ X m+ n = j | X m = i}
为马氏链在时刻m处于状态i条件下, 在时刻 m + n
例1 某学生在高中中每个学期期中、期末六次考试及
第四章 马尔可夫链

第四章 马尔可夫链随机过程在不同时刻下的状态之间一般具有某种关系,马尔可夫(Markov )过程就是描述一类状态之间具有某种特殊统计联系的随机过程.Markov 过程在近代物理学、生物学、管理科学、信息处理与数字计算方法等领域都有重要的应用.按其状态和时间参数是连续的或离散的,它可分为三类:(1)时间、状态都是离散的Markov 过程,称为Markov 链;(2)时间连续、状态离散的Markov 过程,称为连续时间的Markov 链;(3)时间、状态都连续的Markov 过程.本章主要讨论Markov 链,有关连续时间的Markov 链的相关理论将在下章讨论.4.1 马尔可夫链的概念和例子独立随机试验模型最直接的推广就是Markov 链模型,早在1906年俄国数学家Markov 对它进行研究而得名,以后Kolmogorov 、Feller 、Doob 等数学家发展了这一理论.4.1 .1 Markov 链的定义假设Markov 过程{,}n X n T ∈的参数集T 是离散时间集合,即{0,1,2,}T =,相应n X 可能取值的全体组成的状态空间是离散状态集012{,,,}I i i i =.定义 4.1 设有一随机过程{,}n X n T ∈,若对于任意整数n T ∈和任意011,,,n i i i I +∈,条件概率满足11001111{|,,,}{|}n n n n n n n n P X i X i X i X i P X i X i ++++=======则称{,}n X n T ∈为离散时间的Markov 链,简称Markov 链(Markov chains )或马氏链.从定义可以看出:Markov 链具有Markov 性(即无后效性),如果把时刻n 看作现在,那么,1n +是将来的时刻,而0,1,2,,1n -是过去的时刻.Markov 性表示在确切知道系统现在状态的条件下,系统将来的状况与过去的状况无关,而且Markov 链的统计特征完全由条件概率11{|}n n n n P X i X i ++==所决定. 因此,如何确定这个条件概率,是研究Markov 链理论和应用中十分重要的问题之一. 4.1.2 转移概率定义 4.2 称条件概率1(){|}ij n n p n P X j X i +=== (4.1)为Markov 链{,}n X n T ∈在时刻n 的一步转移概率,其中,i j I ∈,简称转移概率(transition probability ).一般地,转移概率()ij p n 不仅仅与状态,i j 有关,而且与时刻n 有关,如果()ij p n 不依赖时刻n 时,则称Markov 链具有平稳转移概率.定义 4.3 若对任意,i j I ∈,Markov 链{,}n X n T ∈的转移概率()ij p n 与n 无关,则称Markov 链是齐次的(或称时齐的)(time homogeneous -),并记()ij p n 为ij p . 下面只讨论齐次Markov 链,并且通常将“齐次”两字省去.定义 4.4 设P 表示一步转移概率ij p 所组成的矩阵,且状态空间{1,2,}I =,则1112121222...........................n n p p p P p p p ⎛⎫ ⎪= ⎪ ⎪⎝⎭称为系统状态的一步转移概率矩阵(transition probability matrix ),它具有性质: (1)0,,ij p i j I ≥∈; (2)1,ijj Ipi I ∈=∈∑.(2)式说明一步转移概率矩阵中任一行元素之和为1,通常称满足性质(1)(2)的矩阵为随机矩阵.定义 4.5 称条件概率(){|},n ij m n m p P X j X i +=== ,,0,1i j I m n ∈≥≥ (4.2)为Markov 链{,}n X n T ∈的n 步转移概率,并称()()()n n ij P p =为Markov 链{,}n X n T ∈的n 步转移矩阵.其中()()0,1n n ij ij j Ip p ∈≥=∑,即()n P 也是一个随机矩阵.特别地,当1n =时,(1)ij ij p p =,此时,一步转移矩阵(1)P P =.我们还规定(0)0,1,iji jpi j ≠⎧=⎨=⎩Markov 链n 步转移概率满足重要的Chapman Kolmogorov -方程(简称C K -方程)。
随机过程(七)-马氏链

第四章Markov过程主要内容⏹离散时间Markov链⏹转移概率⏹平稳分布⏹状态分类⏹极限定理⏹连续时间Markov链⏹Kolomogrov微分方程⏹连续时间马氏过程第一节 离散时间Markov 链一、Markov 链的定义⏹ 直观含义:要确定过程将来的状态,只需知道过程现在的状态就足够了,并不需要知道过程以往的状态。
⏹ 定义:随机过程{,0,1,2,}n X n =⋅⋅⋅称为马氏链(Markov 链),若它只取有限或可列个值E 0, E 1,E 2,…,且对任意的n ≥0及状态011,,,,,n i j i i i -⋅⋅⋅有10011111{|,,,,}{|}n n n n n n P X j X i X i X i X i P X j X i +--+===⋅⋅⋅=====用条件概率的语言来说11011{,,|,,,}{,,|}n n k k n k k n P X j X j X X X P X j X j X ++==⋅⋅⋅===注:1、E 0, E 1,E 2,…称为Markov 链的状态,通常用0,1,2,…来标记E 0, E 1,E 2,…。
{0,1,2,…}称为过程的状态空间,记为S 。
2、若Markov 链的状态是有限的,则称为有限链,否则称为无限链。
2、条件概率11{|}n n n n P X i X i --==,n =1,2,……称为Markov 链的一步转移概率。
3、若转移概率1{|}n n P X j X i -==只与状态,i j 有关,而与时间n 无关,则称该Markov 链是时齐Markov 链,并记1{|}ij n n p P X j X i -===,否则称Markov 链是非时齐的。
矩阵000102101112012()ij ij Sn n n p p p p p p P p p p p ∈⎛⎫ ⎪ ⎪ ⎪== ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭称为转移矩阵。
4、(){|}k ij n k n p P X j X i +===称为k 步转移概率,()k P 称为k 步转移矩阵。
第四章马氏链(1)

P { X m + n = j | X m = i}
=
j1 , j2 ,, jn-1
P { X m+1 = j1 ,, X n+ m-1 = jn-1 , X n+ m = j | X m = i}
而P { Xm+1 = j1 ,, Xn+m-1 = jn-1 , Xn+m = j | Xm = i}
=
=
P { X m = i , X m +1 = j1 , , X n+ m -1 = jn-1 , X n+ m = j |}
P { X m = i} pij1 p j1 j2 p jn-1 j P { X m = i}
P { X m = i}
= pij1 p j1 j2 p jn-1 j
r
注释:如果把转移概率写成矩阵的形式,那么C-K方程
具有以下简单的形式 P(m+k)=P(m)P(k) 步转移概率完全决定。 m, k≧0
特别地,对齐次马氏链有P(n)=Pn, n步转移概率由一
证:
Pij (m + k ) = P { Xn+m+k = j | X n = i}
= P { X n+ m = r , X n+ m + k = j | X n = i}
=
r
= pir (m) prj (k )
r
例4.1.4 求带有两个反射壁的一维随机游动的两步转移 概率矩阵。 1 0 0 0 0 0 1 / 3 1 / 3 1 / 3 0 解 :P ( 2) = P 2 = 0 1 / 3 1 / 3 1 / 3 0 0 1 / 3 1 / 3 1 / 3 0 0 0 0 1 0
数学模型4-1马氏模型

马氏链的遍历极限(II)
若马氏链{ξn: n≥ 0}的状态空间S为有限 集(不妨设S={1,2,L, N}),且转移矩阵 矩阵的每个元素为正,则它存在唯一不 变概率分布π=(π1,π2,L, πN), 满足如下 (指数)遍历性
4
马氏链的遍历极限(III)
令Ti(ω)是(ξ1(ω),L,ξn(ω))首次出现状态 i 的时间.那么μi=E(Ti(ω) | ξ0(ω)=i)就是 一个平均返回(状态i)时间.有结论如下
也就是说随机过程下一时间的发展只和包括当 前时间在内的最近的k个时间的状态有关 而和 这k个时间之前的历史没有关系, (其中k=0, 1, 2, L) ,我们把这样的随机过程叫做k-阶马氏 链.
关于名称的一点说明
参考书中,看到马氏链(过程)的时候要根据上下 文进行判断.有的时候是指普遍的马氏链(包括 高阶,一阶,零阶),有时候特指一阶马氏链. 在大多数情况下,如不特别说明,通常是特指 一阶时齐的马氏链. 如果将一个 k-阶马氏链的相邻 k 个时间的状态 合为一个新的状态: yn=(xn,xn-1,L,xn-k+1) ,则 {yn} 是一个 1-阶马氏链.程
3
时齐马氏链性质 (I)
时齐马氏链由转移概率矩阵和初分布完 全确定,设转移概率矩阵为P=(pij),初 始分布: ,则
时齐马氏链性质 (II)
若记μi(n)=P(ξn=i), μ(n)=(μi(n): i∈ S),即 所谓绝对概率,则:
马氏链的不变分布
记之为pi,j(n,n+k) 矩阵P(n,n+k)=( pi,j(n,n+k) )称为从n出 发的k步转移概率矩阵
高阶马氏过程
若一个随机过程满足:
零(1)阶马氏过程
马氏链

M ( I Q) 1 Q s
s 0
0 R有非 Q 零元素 y ( y1 , y2 , yk r ) Me
e (1,1,,1)T
yi ~ 从第 i 个非吸收状态出发,被某个 吸收状态吸收前的平均转移次数。
1.2 钢琴销售的存贮策略
背景与问题
钢琴销售量很小,商店的库存量不大以免积压资金
年总人数at等级结构等级i12k如助教讲师教授数量分布ntn年属于等级i的人数t01比例分布转移矩阵qpijij是每年从i转至j的比例调入比例年退出总人数退出的比例每年从的人数年调入年调入总人数每年调入i的比例在总调入人数中每年从i转至j的比例人数等级基本模型ij调入退出转移总人数分布可预测已知ijij等级结构一致与马氏链基本方程等级结构at状态概率p转移概率矩阵应满足用调入比例进行稳定控制ijij问题
P(DDH)=P(dddd,dr)=P(ddd)P(ddr) =11/2=1/2
P(RHH)=P(rrdr,dr)=P(rdr)P(rdr) =1/21/2=1/4
i 1 j 1
需求不超过存量,销售需求
3 i
需求超过存量,销售存量
[ jP( Dn j S n i) iP( Dn i Sn i)]P( S n i)
i 1 j 1
n充分大时 P(Sn i) wi
0.632 0.285 0.896 0.263 0.977 0.452 0.857
0.9 0.089 1.0 0.105 1.1 0.122 1.2 0.139
k
状态转移阵
第n周(n充分大)失去销售机会的概率 P P( Dn Sn )
P
0.8 0.073
4 马氏链

P X1 1, X 2 1, X 3 1| X 0 0
7.2 马尔可夫链的状态分类
设 { Xn , n >0 } 是齐次马尔可夫链,其状态空间
I = { 0, 1, 2, … },转移概率是 pij , i , j I
1
7 1
8
1
9
1
1 1/3 121源自3615
2/3
4
1
(1)状态的周期性
2
设初始分布p1 0 P X 0 1 , p0 0 P X 0 0 1 ,
若系统经n级传输后输出为 1,求原发字符也是 1的概率。
0.9 0.1 0.82 0.18 2 解:(1)P , P 2 P , 0.1 0.9 0.18 0.82 0.756 0.244 3 P 3 P 0.244 0.756
马尔可夫链
内容提要
马尔可夫链的概念及转移概率 马尔可夫链的状态分类 状态空间的分解
pij(n) 的渐近性质与平稳分布
马尔可夫过程的四种类型
马尔可夫链(马氏链)
时间、状态都离散 时间连续、状态离散
连续时间马氏链
马尔可夫序列
时间离散、状态连续
时间、状态都连续
连续时间马尔可夫过程(或扩散过程)
[定义] 如集合 { n : n 1, pii(n) > 0 } 非空,则称该集合
的最大公约数 d = d(i) = G.C.D{ n : pii(n) > 0 }为状态 i
则称 { Xn , n T } 为马尔可夫链,简称马氏链。
马氏性 (无后效性)
刘次华版随机4-1离散参数马尔可夫链(1)

成立,称{X(n), n≥0} 为离散参数马氏链.
解放军电子技术学院
卢
定义2 设{X(n):n≥0} 为马氏链,状态空间为 E={0,1,2,…},称条件概率
(n) pij (m) P( X (m n) j | X (m) i)
为马氏链在m 时刻的n 步转移概率. 称
P pij (m)
卢
第m步状态
第m+k步状态
第m+k+s步状态
pij
( s k )
pir
rE
(k )
prj
( s)
解放军电子技术学院
卢
要证
pij ( s k ) pir ( k ) prj ( s ) .
rE
pij
rE
( s k )
P{ X (m k s) j X (m) i}
2) p j (n) pi (n 1) pij
iE
iE
令 PT (n) ( p1 (n), p1(n), )
PT (n) PT (0)P( n) PT (n) PT (n 1) P
卢
解放军电子技术学院
(n) 1) p j (n) pi pij iE
p j (n) PX n j P X 0 i, X n j
rE
pir
rE
(k )
prj .
解放军电子技术学院
卢
(s)
已证
pij ( s k ) pir ( k ) prj ( s ) .
rE
称k步转移矩阵
若记 P
(k )
( pij
(k )
),
第4章马尔可夫链1-2

定义 1 设有随机过程{ X n , n T } ,若对于任意的整数 n T 和任意的 i0 , i1 , , in1 I ,条件概率满足
转移概率矩阵为
q 0 p 0 P 0 q 0 p
设在第k步转移中向右移了x步,向左移了y步,且 经过k步转移状态从i进入j,则
x y k x y j i
从而
k ( j i) k ( n 和 i , j I ,n 步转移概率 ij 具有下列
性质
( n) ( l ) ( n l ) (1) pij pik pkj ; k I
(2) p
( n) ij
k1I
kn1I
pik1 pk1k2 pkn1 j ;
(3) P ( n ) PP ( n1) ; (4) P ( n ) P n .
第4章 马尔可夫链
定义 2.9 设 X t , t T 为随机过程,若对任意正 整数 n 及 t1 t2 , tn , P X (t1 ) x1 , , X t n1 xn1 0 ,且其 条件分布
P X (tn ) xn | X t1 x1 ,, X t n1 xn1 P X ( t n ) xn | X t n 1 x n 1
定义 2 称条件概率
pij (n) P{ X n1 j | X n i }
为马尔可夫链 { X n , n T } 在时刻 n 的一步转移概率,其 中 i , j I ,简称为转移概率。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
pij pij 1 = P { X m+1 = j | X m = i}
称为齐次马氏链的一步转移概率;
P P(1) = pij (1)
a1 a2 P (1) = ai a1 p11 p 21 pi1 a2 p12 p 22 pi 2 L L L L L L aj p1 j p2 j p ij L L L L L L
所以 P { X m+n = j | X m = i} 与m无关。
因此,马氏链的齐次性可写为
P X m1 + n = j | X m1 = i = P X m2 + n = j | X m2 = i
{
}
{
}
定义4.1.4
称条件概率
Pij ( n) P { X m + n = j | X m = i} ,
Pij n, n + 1 = Pij m, m + 1 即马氏链{Xn,n0}的转移概率Pij(n,n+1)与n无关,
则称转移概率具有平稳性,这时,马尔可夫链称为是
若对任意的正整数m,n及任意的ai,aj,有
齐次的。
定理4.1.1:若{Xn}为齐次马氏链,则对任意正整数n,及任意 的i,j,有 P { X m+n = j | X m = i} 与m无关。 证明:
r
注释:如果把转移概率写成矩阵的形式,那么C-K方程
具有以下简单的形式 P(m+k)=P(m)P(k) 步转移概率完全决定。 m, k≧0
特别地,对齐次马氏链有P(n)=Pn, n步转移概率由一
证:
Pij (m + k ) = P { Xn+m+k = j | X n = i}
= P { X n+ m = r , X n+ m + k = j | X n = i}
关的,所以{ Xn,n=0 , 1,2,… }是一马氏链,且是齐次
的。它的一步转移概率和一步转移概率矩阵分别为
1 3 , j = i - 1, i , i + 1, 1 < i < 5 pij = P{Xn+1 = j | Xn = i} = 1, i = 1, j = 2 或 i = 5, j = 4 - 0, j i 2.
j =0 ij j =0
m+n
= j | X m = i}
= P { X m + n = j} | X m = i = 1. j
3.齐次马尔可夫链及一步转移概率
定义4.1.3 若对任意的i,j,有
Pij n,n + 1 = Pij m ,m + 1
即马氏链{Xn,n0}的转移概率Pij(n,n+1)与n无关, 则称转移概率具有平稳性,这时,马尔可夫链称为是 齐次的。
对任意的 n 及 i 0 , i 1 , L , i n , i n + 1 S ,
P {X n +1 = i n +1 X 0 = i 0 , X 1 = i1 , L , X n = i n } 0 i n +1 > i n =1 = P{X n +1 = i n +1 | X n = i n } i n +1 i n in
设订货和进货不需要时间,每天的需求量 立同分布且 P{Yn = j} = a j ( j = 0,1, 2,...) 。
独
4. n步转移概率及C-K方程
称条件概率Pij (m, m + n) P { X n+m = j | X m = i} 为马尔 可夫链在时刻m处于状态i的条件下,在时刻m+n步转移到 状态j的n步转移概率。
在原处;如果Q现在位于1(或5)这点上,则下一时刻就 以概率1移动到2(或4)点上。1和5这两点称为反射壁。 上面这种游动称为带有两个反射壁的随机游动。
1
2
3
4Hale Waihona Puke 5若以 Xn 表示时刻 n 时 Q 的位置,不同的位置就是 Xn 的不同 状态,那么{Xn,n=0,1,2,…}是一随机过程,状态空间 就是 I ,而且当 Xn=i,iI 为已知时 ,Xn+1 所处的状态的概率 分布只与Xn=i有关,而与Q在时刻n以前如何到达i是完全无
=
r
= pir (m) prj (k )
r
例4.1.4 求带有两个反射壁的一维随机游动的两步转移 概率矩阵。 1 0 0 0 0 0 1 / 3 1 / 3 1 / 3 0 解 :P ( 2) = P 2 = 0 1 / 3 1 / 3 1 / 3 0 0 1 / 3 1 / 3 1 / 3 0 0 0 0 1 0
则称{Xn,n=0,1,2,…}为马氏链。
{
}
Pij (m, m + n) P { X n+m = j | X m = i}
称为马氏链在时刻 m 系统处于状态 i 的条件下,经过
n步转移到状态 j 的转移概率。
设{Xn,n0},其状态空间为S,若对于任意的 正整数n和任意的 i0 , i1 , , in+1 , 定义4.1.2
1/ 3 1/ 3
0
5/ 9 2/ 9 1/ 9 2/9 3/9 2/9 1/ 9 2/ 9 5/ 9 0 1/ 3 1/ 3
0 0 1 / 9 1 / 9 1 / 3
5、有限维分布
1.有限维分布 设马氏链{Xn,n≥0},状态空间S,n步转移概率矩阵P(n).
(1)一维分布
有
P { X n+1 = in+1 X 0 = i0 , X 1 = i1 , , X n = in } = P { X n + 1 = in + 1 | X n = i n }
则称{Xn,n0}为马氏链。 注:定义4.1.2与定义4.1.1是等价的。
例4.1.1:记从数1,2, …,N中任取一数为X0,当n1时, 记从数1,2, …,Xn-1中任取一数为Xn,问{Xn,n=0,1, 2,…}是马氏链吗? 证:{Xn,n=0,1,2,…}的状态空间S={i,1iN},
P { X m + n = j | X m = i}
=
j1 , j2 ,, jn-1
P { X m+1 = j1 ,, X n+ m-1 = jn-1 , X n+ m = j | X m = i}
而P { Xm+1 = j1 ,, Xn+m-1 = jn-1 , Xn+m = j | Xm = i}
X n +1 = X n - (X n ) + U n +1
1, x 0 ( x) = 0, x = 0
其中
可证{Xn, n=0,1,2,…}是齐次Markov链,其一步转移概率为 pij = P(Un = j - i + (i)), n 1
14
例:(订货问题)设某商店使用(s,S)订货策略,每 天早上检查某商品的剩余量,设为x,则订购额为:
r
既:“从Xn= i 出发,经时刻m转移到中间状态r,再从 r经k时段转移到 j 状态”这样一些事件的和事件。
=
r
P { X n = i , X n + m = r , X n + m + k = j}
P { X n = i} pir ( m ) prj (k ) P{ X n = i }
P{ X n = i }
i , j S , m 0, n 1
为齐次马氏链{Xn,n≥0}的n步转移概率,并称由pij(n) 组成的矩阵
p11(n ) p12(n ) p 21(n ) p 22(n ) P(n ) = pij (n ) = pi 1(n ) pi 2(n )
=
=
P { X m = i , X m +1 = j1 , , X n+ m -1 = jn-1 , X n+ m = j |}
P { X m = i} pij1 p j1 j2 p jn-1 j P { X m = i}
P { X m = i}
= pij1 p j1 j2 p jn-1 j
称X0的分布 q j (0) = P{ X0 = j}, j = 0,1, 2,
i S
由于Xn, n=0,1,2,…独立同分布,因而
P {X n +1 = j | X n = i } = P {X n +1 = j }
= q j = P{X m +1 = j | X m = i}
所以{Xn}为齐次马氏链。其一步转移概率P:
pij = qj ,
i,j S .
例: M/G/1 排队系统 假设顾客依参数为λ的Poisson过程来到只有一个服务员的服 务站,若服务员空闲来客就立刻得到服务,否则排队等待直 至轮到他。设每名顾客接受服务的时间独立同分布,分布函 数为G(x),且与顾客到达过程相互独立。这个系统称为M/G/1 排队系统. (M--到达的时间间隔服从指数分布, G--服务时间 的分布,1--单个服务员)。 令Xn--第n个顾客结束服务时剩下的顾客数, Un--第n个顾客接受服务的时间内来到服务机构的顾客数,则
L p1j(n ) L L p 2j(n ) L L L L pij (n ) L L L
为齐次马尔可夫链的n步转移概率矩阵。
其中 p ij ( n) 0, p ij ( n) = 1.
a j x
定理4.1.2 设{Xn,n=0,1,…}为齐次马氏链,则对于任 意的正整数k,m,有Pij (m + k ) = Pir (m) Prj (k ) 此方程称为Chapman-kolmogorov(切普曼-柯尔莫哥 洛夫)方程,简称C-K方程.