第5章 Markov链
马尔可夫链-2013

定理2.2 设{Xn,n∈T}为马尔可夫链,则对任意j∈I和n≥ 1,绝对概率pj(n)具有下列性质: (1) pj(n)=pipij(n); (2) pj(n)= pi(n-1)pij ;
iI
iI
(3) PT(n)=PT(0)P(n);
(4) PT(n)=PT(n-1)P. 证明:(1) pj(n)=P{Xn=j}= P{X0=i,Xn=j} =P{Xn=j|X0=i}P{X0=i}= pipij(n).
p0 (0) p
(1) 01
p
(1) 11
p
(1) 11
p
(1) 11
= 1 18 52 52 52 ≈0.28.
26 70 70 70
例2.5 设{Xn,n≥0}是具有3个状态0,1,2的齐次马氏链,一 步转移概率矩阵如右所示: 0 1 2 初始分布pi(0)=P{X0=i}=1/3,i=0,1,2. 0 ¾ ¼ 0 试求(1) P{X0=0,X2=1}; (2) P{X2=1}. 1 ¼ ½ ¼ 解: 先求出二步转移概率矩阵(如右下): 2 0 ¾ ¼ 于是有 0 1 2 (1) P{X0=0,X2=1} 5/16 1/16 0 5/8 =P{X0=0}P{X2=1|X0=0} P2= 1 5/16 1/2 3/16 =p0(0)p01(2)=(1/3)· (5/16)=5/48; 2 3/16 9/16 1/4 (2) p1(2)=P{X2=1} =p0(0)p01(2)+p1(0)p11(2)+p2(0)p21(2) =(1/3)(5/16+1/2+9/16)=11/24.
1,i=j .
例2.1 (一维随机游动) 3 5 1 2 4 设一随机游动的质点, 在如右上图所示的 直线点集I={1,2,3,4,5}作随机游动,并且仅仅在1秒,2秒 …等时刻发生游动.游动的概率规则是:如果Q现在位于点 i(1<i<5), 则下一时刻各以1/3的概率向左或向右移动 一格,或以1/3的概率留在原处; 如果Q现在位于点1(或5) 上,则下一时刻就以概率1移动到点2(或4)上.点1与5称为 反射壁.并称上述这种游动为带有两个反射壁的随机游动. 若以Xn表示时刻n时Q的位置, 不同的位置就是Xn的不同
第5章 Markov过程(2)(使用版)

13/34
三、状态间的关系
1.定义 状态 i可达状态j 简记为 i→j
( n N , 使得pijn) 0
状态 i与状态j互通 i→j 且 2.性质 传递性、对称性 j →i 3.利用首达概率刻画 可达和互通关系
14/34
结论1
解
先按一步转移概率,画出各状态间的传递图
22/34
1/2 1/2 4 1 由图可知 1
1 1/2 2 1 5 1/2 3
状态3为吸收态 故
C1 = {3}为闭集
且
C2 ={1,4} 闭集, C3 ={1,3,4} 闭集,
C4 ={1,2,3,4}闭集, 其中 C1 , C 2 是不可约。
又因状态空间I有闭子集, 故此链为非不可约链。
pij 1, i C
jC
18/34
B. 有关等价类
结论1 等价类若是闭集,则必定是不可约的。 结论2 设C是闭集,当且仅当C中的任何两个状态都互 通时, C是不可约的。 结论3 齐次马氏链不可约的充要条件是它的任意两个 状态均互通。 结论4 包含常返态的等价类是不可约闭集。
19/34
其均值为
1 , 即 1 f ii
平均回到 i 共
1 次 就不再回到 i 了。 1 f ii
也就是说以概率1只有有穷次返回i。
10/34
2.判别
(1)判别是否常返态
fii 1 fii 1
定理3
( piin ) n 1
( piin ) n 1
f ij lim
C1,C2, ,Ck 是互不相交的由正常返态组成的闭集。
第5章 马尔可夫链PPT课件

状态.
精选PPT课件
18
马尔可夫链
一般,一个特定的参保人年理赔要求的次数是参数为λ 的泊松随机变量,那么此参保人相继的状态将构成一个马 尔可夫链,并具有转移概率
但昨天没下雨,那么明天下雨的概率为0.5;如果昨天下雨
但今天没下雨,那么明天下雨的概率为0.4;如果昨、今两
天都没下雨,那么明天下雨的概率为0.2.
假设在时间n的状态只依赖于在时间n-1是否下雨,那么
上述模型就不是一个马尔可夫链.
但是,当假定在任意时间的状态是由这天与前一天两者
的天气条件所决定时,上面的模型就可以转变为一个马尔
令Xn为第n天结束时的存货量,则
XSX-nYn-nY++n1+1=,1,
若Xn≥s, 若Xn<s.
构成的{Xn,n≥1}是Markov链.
例5.11 以Sn表示保险公司在时刻n的盈余,这里的时间以
适当的单位来计算(如天,月等), 初始盈余S0=x显然为
已知,但未来的盈余S1,S2,…却必须视为随机变量,增量
参保人的状态随着参保人要求理赔的次数而一年一年
地变化.低的状态对应于低的年保险金. 如果参保人在上
一年没有理赔要求,他的状态就将降低; 如果参保人在上
一年至少有一次理赔要求,他的状态一般会增加(可见,无
理赔是好的,并且会导致低保险金;而要求理赔是坏的,一
般会导致更高的保险金).
对于给定的一个好-坏系统, 以si(k)记一个在上一年 处在状态i,且在该年有k次理赔要求的参保人在下一年的
矩阵为
p11 p12 p13 p14
P=
p21 p22 p23 p24 0010
0001
例5.5(赌徒的破产或称带吸收壁的随机游动)系统的状态
随机过程 研究生 课程介绍

第0章 课程介绍及课时安排 授课人:刘玉婷 ytliu@ 理学院数学系
提纲
教材及参考书目 主要内容 考试安排
教材及参考书目
教材
《随机过程及其在金融领域中的应用》王军 王 娟 清华大学出版社 北京交通大学出版社
参考书目
《应用随机过程》 林元烈 清华大学出版社 《应用随机过程》柳金甫 李学伟 中国铁道出版 社
第4章 Poisson过程
第6课:3.5 + 4.1 第7课:4.1
复习:第15课 答疑:第16课 – 机械楼N201
考核方式
平时作业 10%
每章之后留习题若干,下次课上交 作业纸作答(不返回) ( )
期末考试 90%
闭卷 仅考所学内容
主要学习内容
第2章 概率空间
第1课:2.1 + 2.2 第2课:2.3 第3课:2.4 arkov链
第9课:5.1 + 5.2 第10课:5.2 第11课:5.3 第12课:5.3 第13课:5.4 第14课:5.5
第3章 随机过程
第4课:3.1 + 3.2 +3.3 第5课:3.4 + 3.6
Markov链(讲义)

^‡:
G
˜m• E ;
i0 , i1 , · · · , in+1 , k (0.1.1)
(ii) é?¿ n 9G
P(Xn+1 = in+1 |X0 = i0 , · · · , Xn = in ) = P(Xn+1 = in+1 |Xn = in ).
E kš²…
48 {x} gñ. ±Ï, P d(x) •8Ü {n ≥ 1 :
½Â 0.2.7 (±Ï, š±Ï) é x ∈ E , ½Â x
px,x > 0}
(n)
•ŒúÏê. (e8Ü´˜8, d(x) := 0.) e d(x) = 1, ¡ x ´š±Ï
, XJÙ¤kG
( aperiodic). ¡ Markov ó´š±Ï
6
Ún 0.2.6 X ØŒ
…=
E عý4f8. x ∈ C … y ∈ C k px,y = 0, n ≥ 1, ù
(n)
y².7‡5(‡y) e C † x Œˆ y gñ.
E •48, K
¿©5 (‡y)e•3 x, y ∈ E , k x → y , Kd·K0.2.5•, y ∈ {x}, ùÒ`²
py,z = 0, l y → z, q x → y, ? G •¹ x x → z ù† z ∈ F gñ. •
`² F ´•¹ x
• 48, =
48, K7k G ⊃ F . eØ,, •3 y ∈ F … y ∈ G,
q G ´48, dÚn0.2.4•, x → y . ù† y ∈ F gñ.
Markov链

只要知道初始概率和一步转移概率就可以描述 马尔可夫链的统计特性
例4.1:无限制随机游动。 设质点在数轴上移动,每次移动一格,向右 移动的概率为p,向左移动的概率为q 1- p,这种 运动称为无限制随机游动,以X n 表示时刻n质点 所处的位置,则{ X n , n T }是一个齐次马尔可夫 链,试写出它的一步和k步转移概率。
定义4.8:若i , 则称常返状态i 为正常返的。 若i=, 则称常返状态i 为零常返的。非周期的 正常返态称为遍历状态。
( f ij( n )与pijn )关系如下: ( 定理4.4: pijn )= fij( k ) p (jjn k )= f ij( n k ) p (jjk ) k 1 k 0 n n
i常返的等价定义: 定义4.7:i常返 fii 1 定理4.6:i常返 gii 1; i 非常返则gii 0 定理4.5:i常返 n 0 pii
(n)
定理4.7 设 i 常返且有周期d,则 lim pii
n
( nd )
d
i
推论 设i常返,则
( (1) i 零常返 lim piin ) 0 n
bi , pij ri a i
j i 1 j i j i 1
4.2 状态的分类
一、状态的分类
依概率性质对状态进行分类。 自状态1出发再返回状态1的可能步数(时 刻)为T={4, 6, 8, 10.……},最大公约数 为2.但2不属于T,好由1出发经两步不能返 回1.把2定义为状态1的周期。
x y k 从而, x y j i k ( j i) k ( j i) x , y 2 2 由于x, y都只能取整数,所以k ( j i ) 必须是偶数。又在k步中哪x步向右,哪y步 向左是任意的,选取的方法有ckx 种。于是 P
5.4.15.4Markov链状态分类准备知识
主要研究马氏链的状态按其概率特性进行分类, 并研究这些分类的判断准则.
2
可达
i, j E, n 0,使得Pij(n) 0,称状态i可达状态j,记为i j; 若不可达,记为i j.
3
互通
若i j,且j i,则称i, j互通,记为i j.
4
互通属性
命题5.1:可达关系与互通关系都具有传递性.
命题5.2
互通关系是一等价关系,即它满足:
(1)自反性: i i; (2)对称性:若i j,则j i; (3)传递性:若i j, j k,则i k;
注 : 利用等价关系,可将马氏链的状态空间分为若干等价类; 每一个等价类内状态彼此互通,不同类间的状态不可能彼此互通; 然而某一类出发以正的概率到达另一类的情形是可能的.
nP
k 1
X n j,Tij k X 0 i
nP
k 1
Xn
j X 0 i,Tij k
P
Tij k
X0 i
nP
k 1
Xn
jXk
j
f (k)1 ij jj
9
可达等价性
命题5.3 i, j E, fij 0 i j.
证明:
0.
P (n) ij
f P n
(k) (nk)
k 1 ij jj
f P ( n) (0) ij jj
0
i j.
10
目录
5.4 Markov链状态分类——准备知识
互通、互通等价类、首达时间、首达分解定理
1
状态分类的意义
1
1/3
1
1/3
2
1/3
1/3
4 1/3 3
1
1/3
《应用随机过程》-课程教学大纲
《应用随机过程》课程教学大纲一、课程基本信息课程代码:16055502课程名称:应用随机过程英文名称:Applied Stochastic Processes课程类别:专业课学时:32学分: 2适用对象:财经类专业本科生考核方式:考试先修课程:微积分、线性代数、概率论二、课程简介中文简介紧抓课程改革核心环节,不断提升教学质量,将“课程思政”作为融合德育与智育的融合主渠道,是逐步实现“立德树人”的综合教育理念的前进方向。
《应用随机过程》是面向经济统计专业三年级学生开设的一门必修课,随机过程通常被视为概率论的动态部分,即研究的是随机现象的动态特征,着重对随时间和空间变化的随机现象提出各种不同的模型并研究其内在的性质与相互联系。
具有较强的理论性。
该学科在社会科学、自然科学、经济和管理等各个领域中都有广泛的应用,培养学生的科学精神,探索自然和人类的奥秘。
英文简介The course Applied Stochastic Processes is one of the compulsory courses for the junior undergraduates majoring in Economic Statistics,which is usually viewed as the dynamic part of probability theories. It focuses on the dynamic feature of stochastic phenomena and emphasizes modeling the stochastic phenomena varying with time and space .Moreover,it explores the inner property and relationship among various models and it is quite theoretical and widely used in social science,natural science,Economic and management science etc.三、课程性质与教学目的本课程是经济统计专业一门应用性很强的专业课。
北京交通大学硕士研究生课程《随机过程》5.3-1
5.3.1 分类的准备知识
一、分类的准备知识 1)可达:
i , j E, n 0,使得p
( n) ij
0,称状态i可达
状态j,记为i j;若不可达,记为i j.
2)互通:
若i j,且j i,则称 i , j互通,记为 i j.
定理5-4:可达关系与互通关系都具有传递性.
5.3.2 状态的类别
例子: 1
7 8
1
9
1 2/3
1
1/3 1
2
1
3
1
6
4,8,12,16, 点1:D 6,12,18,24, d 1 2, 10,14,16,20,
1
5
4
1
p
( 2) 11
0
4k 6l : (6) 点2:D d 2 2 , p 22 0 k 1,2; l 0,1, 4k 6l : (4) 点7:D d 7 2 , p 77 0 k 0,1,; l 1,2
5.3.4 互通等价类
补例2:
1 2 1 2 1 4 1 2 1 2 1 4
设E 1,2,3,4
1/2
1
0 0 0 0 P 1 1 4 4 0 0 0 1 解:
分类的准备知识 状态类别 分类属性 互通等价类 常返性的判定
5.3 状态的分类
0、前言 1)研究内容 主要研究马氏链的状态按其概率特性进行分类,并 研究这些分类的判断准则. 2)例证—例5-5 1
1
1/3
2
1/3 1/3
3
1/3
4
1/3 1
1/3
观察得到: 状态 2 、 3 、 4 有进有出,而 且经过有限次转移都能到 达状态 1 ;而一旦到达 1 后 ,则会永远停止在 1 上,不 再转移出来;由此可见各 状态的概率性质是不同的.
No.12-第5章-马尔可夫预测的基本原理
P(2)
0.76 0.72
p (1) 11
p (1) 21
p (1) 11
p (1) 11
p p (1) (1) 12 21
p p (1) (1) 22 21
p p (1) (1) 11 12
p p (1) (1) 21 12
p (1) 12
p (1) 22
p2 (1)
p2 (2) ?
p22
p2 (2) p1(1) p12 p2 (1) p22
P(2) ( p1(2)
p2 (2)) ( p1(1)
p2
(1))
p11 p21
p12
p22
P(1)P P(0)PP P(0)P2
[0.8
0.2]
0.76 0.72
以一个月为单位,经观察统计,知其从某个月份到下月份, 机床出现故障的概率为0.3。在这一段时间内,故障机床经维修 恢复到正常状态的概率为0.9。
p12=0.3
1
2
p21=0.9
0.7 P 0.9
两步状态转移概率:
0.3 0.1
即有 p11 0.7, p12 0.3, p21 0.9, p22 0.1
状态转移概率矩阵
状态
1
P
1 p11 2 p21
2
状态
1
p12
p22
1 0.8 2 0.6
2
0.2
0.4
例5.1 考察一台机床的运行状态。机床的运行存在正常和故
障两种状态 S 1, 2 。机床在运行中出现故障:1->2;处于故
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
������, ������ ∈ ������; ∀������ ∈ ������ 为随机矩阵 ,若 ������������������ ≥ 0(������, ������ ∈ ������) ,且对
定义 5.1.4 称矩阵 ������ = ������������������ ∀ ������ ∈ ������,有
0 0 0 0
11
5.1 基本概念
例 5.1.8 (Wright-Fisher 遗传模型)基因控制着生物的特征,它们是成 对出现的.控制同一特征的不同基因称为等位基因,记这对等位基因为 ������和������ , 分别称为显性的与隐性的.在一个总体中基因������和������ 出现的频率称为基因频率, 分别记为������和1 − ������. 设总体中的个体数为2������,每个个体的基因按基因������的基因频率的大小,在下 一代中转移成为基因 ������.即如果在第 ������ 代母体中基因������出现了 ������ 次,基因 ������ 出现了
6
5.1 基本概念
例 5.1.3 在任意给定的一天,加里的心情或者是快乐的(cheerful,C),或 者是一般的(so-so,S),或者是忧郁的(glum,G). 如果今天他是快乐的,则明天 他分别以概率 0.5,0.4,0.1 是 C,S,G.如果今天他感觉一般,则明天他分别 以概率 0.3,0.4,0.3 为 C,S,G.如果今天他是忧郁的,则明天他分别以概率 0.2,0.3,0.5 为 C,S,G. 以 ������������ 记加里在第 ������ 天的心情, 则 ������������ , ������ ≥ 0 是一个三个状态的马尔可夫链 (状态 0=C,状态 1=S,状态 2=G),具有转移概率矩阵 0.5 ������ = 0.3 0.2 0.4 0.4 0.3 0.1 0.3 0.5
3
5.1 基本概念
定义 5.1.1 随机过程{������������ , ������ = 0, 1, 2, ⋯ }称为 Markov 链,若它的状态空间 可列,并且对任意的������ ≥ 0,及任意状态������ , ������, ������0 , ������1 , ⋯ , ������������−1 ,有 ������ ������������ +1 = ������ ������������ = ������ , ������������−1 = ������������−1 , ⋯ , ������1 = ������1 , ⋯ , ������0 = ������0 = ������ ������������ +1 = ������ ������������ = ������ 定义 5.1.2 移概率. 定义 5.1.3 若对 ∀ ������, ������ ∈ ������, ������������������ ������ ≡ ������������������ , 即 ������������������ 与 ������ 无关, 则称 ������������ , ������ ≥ 0 为 时齐马尔可夫链;否则,就称为非时齐的. ∀ ������, ������ ∈ ������ ,称 ������ ������������ +1 = ������|������������ = ������ ≜ ������������������ ������ 为 ������ 时刻的一步转
随机过程 ������ ������ , ������ ∈ ������ , 若对任意������1 < ������2 < ⋯ < ������������ , ������������ , 1 ≤ ������ ≤ ������, 及 ������ ⊂ ℝ , 总有 ������ ������ ������ ∈ ������|������ ������1 = ������1 , ������ ������2 = ������2 , ⋯ , ������ ������������ = ������������ = ������ ������ ������ ∈ ������|������ ������������ = ������������ 则称此过程为马尔可夫过程(Markov process),简称马氏过程. 称 ������ ������, ������ ; ������, ������ = ������ ������ ������ ∈ ������|������ ������ = ������ ������ < ������ 为转移概率函数. ������ ������ 的取值全体构成的集合记为 ������,称为状态空间. 离散时间参数 ������ = 0, 1,2, ⋯ = ℕ0 ,状态空间 ������ = 0, 1,2, ⋯ 可列的马尔可 夫过程,通常称为马尔可夫链(Markov chain),简称马氏链.
������ ������ ������2������ ������������ (1 ������ 2������−������
������ ∈������ ������������������
������ × ������
= 1.
5
5.1 基本概念
例 5.1.1(天气预报)假设明天下雨的机会只依赖于前一天的天气条件,即 今天是否下雨,而不依赖过去的天气条件.再假设如果今天下雨,那么明天下雨 的概率为 ������ ;如果今天没有下雨,那么明天下雨的概率为 ������ . 如果下雨, 我们假定过程在状态 0; 如果不下雨, 我们假定过程在状态 1. 那 么,上面的内容是一个两个状态的马尔可夫链,其转移概率矩阵给定为 ������ = ������ ������ 1 − ������ 1 − ������
7
5.1 基本概念
例 5.1.4(赌徒的破产或称带吸收壁的随机游动)系统的状态是0~������,反映 赌博者在赌博期间拥有的金钱数额,当他输光或拥有钱数为������时,赌博停止,否 则他将持续赌博.每次以概率������赢得 1,以概率������ = 1 − ������输掉 1.这个系统的转 移矩阵为: 1 ������ ������ = ⋮ 0 0 ������������ ,������ +1 = ������, 0 0 0 0 ������ 0 ⋮ ⋮ ⋮ 0 0 0 0 0 0 ⋯ 0 0 0 ⋯ 0 0 0 ⋮ ⋮ ⋮ ⋯ ������ 0 ������ ⋯ 0 0 1
例 5.1.2(一个通信系统)考察一个传送数字 0 和 1 的通信系统.每个数字 的传送必须经过几个阶段,在每个阶段有一个概率������使进入的数字在离开时不改 变.以 ������������ 记第 ������ 个阶段进入的数字,则 ������������ , ������ = 0,1,2, ⋯ 是一个两个状态的马 尔可夫链,具有转移概率矩阵 ������ ������ = 1 − ������ 1 − ������ ������
例 5.1.6(自由随机游动)设一个球在全直线上做无限制的随机游动,它的 状态为0, ±1, ±2, ⋯.它仍是一个 Markov 链,转移矩阵为: ⋯ ⋯ ⋯ ⋯ ������������ ,������ +1 = ������, ⋮ ⋮ ⋮ ⋮ ⋮ ������ 0 ������ 0 ⋯ 0 0 ������ 0 ������ ⋯ 0 ⋮ ⋮ ⋮ ⋮ ⋮ 0 0 0 0 ⋯ ������ 0 0 0 0 ⋯ 0 ⋮ ⋮ ⋮ ⋮ ⋮ ������������ ,������−1 = 1 − ������, ⋮ ⋮ 0 0 0 0 ⋮ ⋮ 0 ������ ������ 0 ⋮ ⋮ ⋮ 0 0 ⋮ 0 ������ ⋮ ⋯ ⋯ ⋯ ⋯
4
5.1 基本概念
本课程只讨论时齐 Markov 链,并且简称为 Markov 链. 记 ������ = ������������������ ,称 ������ 为 ������������ , ������ ≥ 0 ������00 ������10 ������ = ������������������ = ⋮ ������������0 ⋮ 易见,������������������ (������ , ������ ∈ ������)有性质: 1 2
������ ∈������
的一步转移概率矩阵: ������01 ������02 ⋯ ������11 ������12 ⋯ ⋮ ⋮ ⋮ ������������1 ������������3 ⋯ ⋮ ⋮ ⋮
(5.1.1)
������������������ ≥ 0, ������������������ = 1,
应用随机过程(Applied Stochastic Processes)
第5章 Markov链
1
第5章 Markov链
5.1 基本概念 5.2 状态的分类及性质 5.3 极限定理及平稳分布 5.4 Markov链的应用 5.5 连续时间Markov链
2
5.1 基本概念
5.1.1 Markov 链的定义及一些例子
������ =
������ = 0, ±1, ±2, ⋯பைடு நூலகம்
10
5.1 基本概念
例 5.1.7 (图上的简单随机游动)设有一蚂蚁在如图 5-1 所示的架子上爬 行,当两个节点相邻时,蚂蚁将爬向它邻近的一点,并且爬向任何一个邻近节点 的概率是相同的. 则此 Markov 链的转移矩阵为: 0 1 2 1 ������ = 4 0 1 2 0 1 4 0 1 0 2 1 0 2 1 0 4 1 0 1 0 2 0 0 0 0 1 4 0 0 1 0 0 0 0 1 2 0 图 5-1 例 5.1.5 图示
������+������ ×(������+������)
������������ ,������−1 = 1 − ������,