第四章马氏链(1)

合集下载

4.马尔可夫链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}

马尔可夫链

马尔可夫链
2020年5月21日星期四
例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马氏链

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}为马尔可夫过程。
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 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方程.
相关文档
最新文档