随机过程-第四章 更新过程
随机过程第4章

2、更新函数 令m(t)=E[N(t)],称m(t)为更新函数。显然m(t)是单调递 增的,因而其反函数m-1(t)存在 Theorem: m t Fn t
n 1
Pr oof : m t E N t nP N t n
0 t
是方程 最后证明解的惟一性,设 K
K (t ) H (t ) K (t s)dF (s)
0 t
的解,且满足有界性条件,则
H F K K
连续代换有
H F (H F K ) H F H F (F K ) K H F H F H F K H F H F K
定理二:如果对于计数过程任意相继出现的两个 质点的点间间距 Xn是相互独立,且服从同一个指数分布:
et t 0 f t t0 0
则计数过程构成强度为λ 的泊松过程
自然延伸—更新过程
令{N (t ), t 0}是一个计数过程,而以 X n记这个过程的第 n 1个和第 n个事件(质点)之间的 时间,n 1.
0 t t t
定义(更新方程)如下形式的积分方程称 为更新方程
K (t ) H (t ) K (t s)dF (s)
0
t
其中H(t),F(t)为已知,且当t<0时, H(t),F(t) 均为0,当H(t)在任何区间上有界时称此方程为 适定更新方程,简称更新方程。
二、更新方程的解
n2
H (t ) F H (t ) ( ( Fn 1 F )) H (t )
随机过程习题集-第四章马尔可夫过程

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 } 为此马氏链的平稳分布。
随机过程-第四章 更新过程

N (t ) 的情况。 t
为考虑 N (t ) 的发散速度,我们先考虑到达时刻 TN (t ) ( TN (t ) 表示在时刻 t 或时刻 t 之前 最后一次更新发生的时刻,以此类推,则 TN (t )1 表示在时刻 t 之后第一次更新发生的时刻) 。 利用 TN (t ) 和 TN (t )1 ,我们提出并证明以下命题。
X
i 1
n
i
n
Tn n
以概率 1 成立。因 0 ,这意味着当 n 时,Tn ,这即是说无穷多次更新只可能 在无限长的时间内发生,因此在有限时间内至多只能发生有限次更新。因此,更新过程亦可 写成
N (t ) max n, Tn t
4.2 N (t ) 的分布
这其中利用了 X n , n 1, 2, 的独立同分布性质,这里 [1 F (b)] (0,1) 。又因为
k
Tmk t Tk T0 t , T2k Tk t ,, Tmk T( m1)k t
而且更新区间(相当于时间间隔)服从独立同分布,即
N (t ) 的分布至少在理论上能够得到,首先我们注意这样一个重要的关系:到时刻 t 为
止的更新次数大于或等于 n 当且仅当在 t 之前或在时刻 t 发生第 n 次更新,即
N (t ) n Tn t
所以
P N (t ) n P N (t ) n P N (t ) n 1 P Tn t P Tn1 t
N (t ) sup n, Tn t
定义 4.1 更新过程:计数过程 N (t ), t 0 称为更新过程。
在更新过程中我们将事件发生一次叫做一次更新, 从而 X n 就是第 n 1 次与第 n 次更新 相距的时间,Tn 表示第 n 次更新发生的时刻, 而 N (t ) 就是 t 时刻或 t 时刻之前发生的总的更 新次数。更新过程一个典型的例子是机器零件的更换。 我们首先要回答是第一个问题是在有限时间内是否会有无限多次更新发生。答案是不 会发生这种情况的概率为 1。由强大数定律可知
随机过程第4章Markov过程(PDF)

第四章 Markov 过程本章我们先讨论一类特殊的参数离散状态空间离散的随机过程,参数为0{0,1,2,}T N ==L ,状态空间为可列{1,2,}I =L 或有限{1,2,,}I n =L 的情况,即讨论的过程为Markov 链。
Markov 链最初由Markov 于1906年引入,至今它在自然科学、工程技术、生命科学及管理科学等诸多领域中都有广泛的应用。
之后我们将讨论另一类参数连续状态空间离散的随机过程,即纯不连续Markov 过程。
§4.1 Markov 链的定义与性质一、Markov 链的定义定义 4.1设随机序列{;0}n X n ≥的状态空间为I ,如果对0n N ∀∈,及0110011,,,,,{,,,}0n n n n i i i i I P X i X i X i +∈===>L L ,有:11001111{,,,}{}n n n n n n n n P X i X i X i X i P X i X i ++++=======L (4.1.0)则称{;0}n X n ≥为Markov 链。
注1:等式(4.1.0)刻画了Markov 链的特性,称此特性为Markov 性或无后效性,简称为马氏性。
Markov 链也称为马氏链。
定义4.2 设{;0}n X n ≥为马氏链,状态空间为I ,对于,i j I ∀∈,称1{}()ˆn n i j P X j X i p n +===为马氏链{;0}n X n ≥在n 时刻的一步转移概率。
注2:一步转移概率满足:()0,,()1,i j i jj Ip n i j Ipn i I ∈≥∈=∈∑若对于,i j I ∀∈,有1{}()ˆn n i j i j P X j X i p n p +===≡即上面式子的右边与时刻n 无关,则称此马氏链为齐次(或时齐的)马氏链。
设{}0()(0),p i P X i i I ==∈,如果对一切i I ∈都有00()0,()1i Ip i p i ∈≥=∑,称0()p i 为马氏链的初始分布。
第四章-马尔可夫链-随机过程

计算 n 步转移概率的方法。
切普曼一柯尔莫哥格夫方程:对一切n,m 0,一切 i,j,有(4.2.1)
P nm ij
Pikn Pkmj
k0
证明:
P nm ij
P{ X nm
j|
X0
i}
P{Xn k | X0 i}P{Xnm j | Xn k, X0 i}
顾客数构成一个泊松过程。所以,
Pi, j
e t (t )i1 j dG(t ), j 1,
0
(i 1 j)!
i 1
这是因为若一个来客发现有 i 个人在系统中,那么下一个来客将
发现人数为 i+1 减去已服务完毕的人数,易知有 i+1-j 个人被服
务完毕的概率(对相继来到之间的时间取条件)等于上式的右端。
0
0
0 P43
例 4.1(b) G/M/1 排队系统。假设顾客 依照一个任意的更新过
程来到一个单服务台的服务中心,来到间隔分布为 G。进一步
假设服务分布是指数分布,参数为。若以 Xn 记第 n 个顾客来
到时见到系统中的顾客数,以 Yn 记第 n 个顾客与第(n+1)个顾客
不可被 d 整除的 n 有 Piin 0,且 d 是具有此性质的最大整数(d 是
{n : Piin 0}的最大公约数)。(若对一切 n>0, Piin 0,则定义 i 的周 期是无穷大。)具有周期 1 的状态称为非周期的(aperiodic)。以 d(i)记 i 的周期。
例设马尔可夫链的状态空间I={1,2,,9}, 转移概率如下图
P nm ij
随机过程-第四章 更新过程

N n ,则在观察 X1 ,, X n 之后与观察 X n1 , X n 2 , 之前我们停止观察。
例 4.2(a)掷硬币试验的停时:设 X n , n 1, 2, 相互独立且使得
1 P X n 0 P X n 1 , 0表示反面,1表示正面 2
所以
mk j
t Tmk t
( m 1) k 1 n mk
P T
n
t kP Tmk t
综合以上得
M (t ) Fn (t ) P Tn t
n 1 n 1
P Tn t P Tn t
n 1 nk
验,设成功的概率为 P ,失败的概率为 1 P 。以试验成功作为事件(更新) ,则此过程是 更新过程。求 P{N (t ) n} 和更新函数 M (t ) 。 解:依题意易知,过程的时间间隔 X i 服从独立的同几何分布,即
P{X i n} P(1 P)n1 , i 1, 2,, n 1, 2,
则第 k 次成功(更新)发生的时刻 Tk
X
i 1
k
i
具有负二项分布
C k 1P k (1 p)n k , n k P{Tk n} n 1 0, n k
上式表示:在第 n 次贝努利试验取得第 k 次成功(更新)的概率。 因此
P{N (t ) k} Fk (t ) Fk 1 (t ) P{Tk t} P{Tk 1 t}
命题 4.3 当 t 时,以概率 1 保证
证明:因为 TN (t ) t TN (t ) 1 ,于是有
随机过程-第四章

三、分类 马氏过程{X(t),t∈T }按参数 T 和状态空间 E 的情况一般分三类
(1)T离散 如例2
E离散的马氏过程,称为马氏链
(2 )T 连续,E 离散称为马氏过程 如:例 1 ,电话…
它表示,已知 n 时刻处于状态 i,经 k 个单位时间 后处于(转移到)状态 j 的概率(条件概率)
一般 pij ( n, n k ) 与 n 有关,如果不依赖于 n,则称过 程{X(n),n=0,1,2… }为时齐(齐次)马氏链,即 有 pij ( n, n k ) = pij ( k ) , k≥1 的马氏链是时齐马氏链,
或说:如果过程{ X(t),t∈T }的 n 维联合分布函 数可表示为
Fn ( x1 , x 2 , , x n ; t1 , t 2 , , t n )
n
F ( x k , t k ) ,n=1,2… =k 1
则称{ X(t),t∈T }为独立随机过程。
证明:设 0≤t1<t2<… tn<t∈T 由条件 X(t1),X(t2),…X(tn)相互独立 故事件 ( X ( t 1 )
一般规定
1, i j pij (0) ij 0, i j
说明:k 步转移概率,可由一步转移概率矩阵获 得, 这说明一步转移概率矩阵是马氏过程最基本 的,它完全确定了链的状态转移的统计规律。
花粉位置,用平面直角坐标系描述,t 时刻花粉的位 置。X(t)( 模标) Y(t)(纵标)t≥0,都是马氏过程。
二、马氏过程。“无后效性”的特点在数学上的定义
第四章随机过程

设状态空间为 S ,一步转移概率为 P ,初始分布为 p i = P( X 0 = i ), i ∈ S 的齐次
Markov 链 {X n , n ≥ 0},令 Pij( n ) = P ( X n + m = j X m = i ) = P (X n = j X 0 = i ), n ≥ 2 ,表示
证:
5
Pij( n ) = P( X n = j X 0 = i ) = ∑ P(τ ij = l , X n = j X 0 = i )
l =1 n n
= ∑ P(τ ij = l X 0 = i ) P( X n = j τ ij = l , X 0 = i )
l =1 n
= ∑ P(τ ij = l X 0 = i ) P( X n = j X 0 = i, X 1 ≠ j , L X l −1 ≠ j , X l = j )
4.4 常返与瞬过
在事件 {X 0 = i}上引入一个重要的概率 f ij( n ) ,表示从 i 出发在 n 步转移时首次 到达 j 的概率。用式子表示即是
f ij( 0) = 0, f ij( n ) = P( X n = j , X k ≠ j , k = 1, L n − 1 X 0 = i ) 。
i =1
(m) 定理 4.3.3 的一个直接推论是: 若 Pji > 0 ,存在正整数 N 使得对所有的 ( m + nd ( i )) n > N 恒有 Pji > 0。
定理 4.3.4:设 P 为不可约、非周期、有限状态 Markov 链的一步转移概率矩阵, 则存在正整数 N 使得当 n > N 时, n 步转移概率矩阵 P ( n ) 的所有元素都大于 0。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
4.1 更新过程定义
上一章我们看到泊松过程的到达时间间隔是服从独立同分布的指数随机变量。现将其 进行推广,考虑到达时间间隔服从独立同分布,但分布函数任意,这样得到的计数过程称为 更新过程。 设 X n , n 1, 2, 是一列服从独立同分布的非负随机变量,分布函数为 F ( x) ,为避 免显而易见的平凡情形, 假设 F (0) P X n 0 1 。 将 X n 解释为第 n 1 个与第 n 个事件 之间相距的时间,记 E ( X n ) 有 0 。令 Tn
这其中利用了 X n , n 1, 2, 的独立同分布性质,这里 [1 F (b)] (0,1) 。又因为
k
Tmk t Tk T0 t , T2k Tk t ,, Tmk T( m1)k t
而且更新区间(相当于时间间隔)服从独立同分布,即
P 1 因 此 存 在 a 0 , 使 得 P Xn a 0 , 从 而 由 于 F( 0 ) X n 0 , P X n a 1 。而 F (a) P X n a P X n a P X n a
为 避 免 因 可 能 的
N (t ) sup n, Tn t
定义 4.1 更新过程:计数过程 N (t ), t 0 称为更新过程。
在更新过程中我们将事件发生一次叫做一次更新, 从而 X n 就是第 n 1 次与第 n 次更新 相距的时间,Tn 表示第 n 次更新发生的时刻, 而 N (t ) 就是 t 时刻或 t 时刻之前发生的总的更 新次数。更新过程一个典型的例子是机器零件的更换。 我们首先要回答是第一个问题是在有限时间内是否会有无限多次更新发生。答案是不 会发生这种情况的概率为 1。由强大数定律可知
定义 4.2 的理解:我们依次观察诸 X n ,以 N 记在停止观察之前所观察的次数。若
N n ,则在观察 X1 ,, X n 之后与观察 X n1 , X n 2 , 之前我们停止观察。
例 4.2(a)掷硬币试验的停时:设 X n , n 1, 2, 相互独立且使得
1 P X n 0 P X n 1 , 0表示反面,1表示正面 2
且因为随机变量 X n , n 1, 2, 服从独立同分布且分布函数为 F ( x) ,记 Fn 为 Tn 的分 布函数,则 Fn 是 F 自身的 n 次卷积。因此可得
P N (t ) n Fn (t ) Fn1 (t )
令 M (t ) E[ N (t )] ,称 M (t ) 为更新函数。
TN (t ) N (t )
N (t ) 时,
TN (t ) N (t )
。但由于 t 时 N (t ) ,所以当 t 时,
TN (t ) N (t )
。
又
TN (t )1 N (t )
TN (t )1 N (t ) 1 。 ,类似地可推得当 t 时, N (t ) 1 N (t ) N (t ) TN (t ) 1
P Tk T0 t , T2k Tk t ,, Tmk T( m1) k t [ P Tk T0 t]m [ P Tk t]m
所以
P Tmk t [ P Tk t]m (1 )m
对任意整数 j 0 ,有
T
t
使所发生的更新总数 N () 有限的唯一方法是有一个到达时间间隔为无大,即
P N () P X n , n P X n P X n 0 ?(对比有 n1 n1
限 t 的情形) 于是当 t 趋于无穷时 N (t ) 趋于无穷。接下来,我们还需要进一步考虑的是 N (t ) 趋于无 穷的速度,即要考虑 lim
如果令
N min n, X1 X n 10
则 N 是一个停时,即在连续掷硬币试验过程中,当出现硬币正面次数达到 10 次时停止 试验。
定理 4.1 瓦尔德等式 :若 X1 , X 2 , 是独立同分布的随机变量,期望有限,且 N 是
X1 , X 2 , 的停时,使得 E[ N ] ,则
n
0
xdF ( x) 为相继发生的两事件的时间间隔的均值,且
X , n 1, T
i 1 i
0
0 ,显然, Tn 表示第 n 个事件发生的时刻。因为至
时刻 t 为止已发生的事件个数等于使第 n 个事件在时间 t 或 t 之前发生的最大的 n 值, 所以到 时刻 t 时刻已发生的事件的个数 N (t ) 为
t
N (t ) 的情况。 t
为考虑 N (t ) 的发散速度,我们先考虑到达时刻 TN (t ) ( TN (t ) 表示在时刻 t 或时刻 t 之前 最后一次更新发生的时刻,以此类推,则 TN (t )1 表示在时刻 t 之后第一次更新发生的时刻) 。 利用 TN (t ) 和 TN (t )1 ,我们提出并证明以下命题。
P Tn t kP Tmk t
m 1
P Tn t k (1 ) m
n 1 m 1
k 1
P Tn t
n 1
k 1
k
因此,命题 4.2 得证。
例 4.1 考虑一个时间离散的计数过程 {N (t ), t 1, 2,} ,在败的概率为 1 P 。以试验成功作为事件(更新) ,则此过程是 更新过程。求 P{N (t ) n} 和更新函数 M (t ) 。 解:依题意易知,过程的时间间隔 X i 服从独立的同几何分布,即
P{X i n} P(1 P)n1 , i 1, 2,, n 1, 2,
利用极限的夹逼定理可知,当 t 时,
t 。因此命题 4.3 得证。 N (t ) 1 1 。因此, 称为
命题 4.3 的解释:以概率 1 保证,长时间后更新发生的速率将等于 更新过程的速率。
接下来我们感兴趣的是更新过程平均速度的期望
M (t ) 1 是否也同样收敛于速率 。然 t
k 1 k n k Cn 1 P (1 p ) nk t n k 1
C
t
k n 1
P k 1 (1 p)n k 1
更新函数
M (k ) E[ N (t )] rP{N (t ) r}
n 1
k
4.3 更新定理
在讨论更新定理之前, 我们先讨论若干极限定理, 这对于我们更好地理解更新定理有所 帮助。 若以 N () lim N (t ) 记所发生的更新总数,易知以概率 1 保证 N () 。这是因为
P X n a P X n a
造
成
的
F (a) P X n a P X n a P X a 1 的情况,不妨取 0 b a ,则有 F (b) P X n b P X n a 1
又对任意固定的 t ,总能找到一正整数 k ,使得 kb t ,所以
k 1
P Tn t [ P Tn t P Tn t ]
n 1 nk n2k
k 1
2 k 1
3 k 1
P Tn t
n 1 k 1 n 1
k 1
( m 1) k 1 n mk
m 1
P Tn t
命题 4.3 当 t 时,以概率 1 保证
证明:因为 TN (t ) t TN (t ) 1 ,于是有
N (t ) 1 , ( EX n ) 。 t
TN (t ) N (t )
其中,
T t N (t ) 1 N (t ) N (t )
表示前 N (t ) 个事件(或更新)到达时间间隔的均值,由强大数定律可得,当
E[ X n ] E[ N ]E[ X ]
n 1 N
证明:令
1, N n In 0, N n
则有
N
X
n 1
n
X n In
n 1
因此
E[ X n ] E[ X n I n ] E[ X n I n ]
X
i 1
n
i
n
Tn n
以概率 1 成立。因 0 ,这意味着当 n 时,Tn ,这即是说无穷多次更新只可能 在无限长的时间内发生,因此在有限时间内至多只能发生有限次更新。因此,更新过程亦可 写成
N (t ) max n, Tn t
4.2 N (t ) 的分布
n 1 n 1 n 1 n 1 n 1
这其中由于 I n 非负,求期望与求和顺序交换是合理的。 注:对于泊松过程, M (t ) E[ N (t )] t 。 (请同学们自证)
命题 4.2 M (t ) 是不减函数,且对一切 0 t , M (t ) 。
证明:因为 N (t ) 是不减函数,所以 M (t ) 也是不减的。接下来证明 M (t ) 的有限性。
Tk t Tk kb X1 b, X 2 b, , X k b
c
(
思
考
如
果
是
X1 b, X 2 b, , X k b 成立吗?)
于是
P Tk t 1 P X1 b, X 2 b, , X k b 1 [1 F (b)]k 1
N (t ) 的分布至少在理论上能够得到,首先我们注意这样一个重要的关系:到时刻 t 为
止的更新次数大于或等于 n 当且仅当在 t 之前或在时刻 t 发生第 n 次更新,即
N (t ) n Tn t