随机过程-第四章 更新过程
随机过程第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。
随机过程-课件-第四章
随机过程-课件-第四章第四章Poion过程4.1齐次Poion过程到达时间间隔与等待时间的分布1、定理4-1强度为的齐次Poion过程{Nt,t0}的到达时间间隔序列某n,n1,2,是独立1同分布的随机变量序列,且是具有相同均值证:事件的指数分布。
即事件某1t等某1t发生当且仅当Poion过程在区间0,t内没有事件发生,价于{Nt0},所以有P(某tt)P(Nt0)et因此,某1具有均值为1的指数分布,再求已知某1的条件下,某2的分布。
P(某2t|某1)P(在,+t内没有事件发生|某1)(由独立增量性)(由平稳增量性)et上式表明P(在,+t内没有事件发生)P(在0,t内没有事件发生)某2与某1相互独立,而且某2也是一个具有均值为1的指数分布的随机变量,重复同样的推导可以证明定理4-1的结论。
2、定理4-2等待时间Sn服从参数为n,的分布,即分布密度为f(t)et证:(t)n1,t0(n1)!因为第n个事件在时刻t或之前发生当且仅当到时间t已发生的事件数目至少是n,即事件NtnSnt是等价的,因此P(Snt)P(Ntn)ejnt(t)jj!上式两边对t求导得Sn的分布密度为f(t)etjnj1(t)j(t)etj!(j1)!jnet(t),t0(n1)!n1注:定理4-2又给出了定义Poion过程的另一种方法。
从一列均值为1/的独立同分布的指数随机变量序列某n,n1出发,定义第n个事件发生的时刻为Sn,则Sn某1某2某n这样就定义了一个计数过程,且所得计数过程Nt,t0就是参数为的Poion过程。
3、定理4-3条件随机变量(某1证:对|Nt1)U(0,t),即在区间0,t内为均匀分布。
t,(某1|Nt1)的分布函数为P(某1,Nt1)P(某1|Nt1)P(Nt1)P(在0,内有一个事件发生,在,t内没有事件发生)P(Nt1)P(在0,内有一个事件发生)P(在,t内没有事件发生)P(Nt1)P(N1)P(Nt0)P(N1)ee(t)tett这说明(某14、顺序统计量|Nt1)在0,t上服从均匀分布。
《应用随机过程》-课程教学大纲
《应用随机过程》课程教学大纲一、课程基本信息课程代码: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.三、课程性质与教学目的本课程是经济统计专业一门应用性很强的专业课。
随机过程_课件---第四章
随机过程_课件---第四章第四章 Poisson 过程4.1 齐次Poisson 过程到达时间间隔与等待时间的分布1、定理4-1强度为λ的齐次Poisson 过程{,0}t N t≥的到达时间间隔序列{},1,2,n X n = 是独立同分布的随机变量序列,且是具有相同均值1λ的指数分布。
证:事件{}1X t >发生当且仅当Poisson 过程在区间[]0,t 内没有事件发生,即事件{}1X t >等价于{0}tN =,所以有()(0)t t t P X t P N e λ->===因此,1X 具有均值为1λ的指数分布,再求已知1X 的条件下,2X 的分布。
(](](]211(|)(|)((0tP X t X s P X s P P e λ->====在s,s+t 内没有事件发生(由独立增量性)在s,s+t 内没有事件发生)(由平稳增量性)在,t 内没有事件发生)上式表明2X 与1X 相互独立,而且2X 也是一个具有均值为1λ的指数分布的随机变量,重复同样的推导可以证明定理4-1的结论。
2、定理4-2等待时间n S 服从参数为n ,λ的Γ分布,即分布密度为1()(),(1)!n tt f t e n λλλ--=- 0t ≥证:因为第n 个事件在时刻t 或之前发生当且仅当到时间t 已发生的事件数目至少是n ,即事件{}{}t n N n S t ≥?≤是等价的,因此()()()!j tn t j nt P S t P N n ej λλ∞-=≤=≥=∑上式两边对t 求导得n S 的分布密度为11()()()!(1)!(),0(1)!j j tt j nj nn tt t f t e e j j t et n λλλλλλλλλ-∞∞--==--=-+-=≥-∑∑注:定理4-2又给出了定义Poisson 过程的另一种方法。
从一列均值为1/λ的独立同分布的指数随机变量序列{},1n X n ≥出发,定义第n 个事件发生的时刻为n S ,则12n n S X X X =+++这样就定义了一个计数过程,且所得计数过程{},0t N t ≥就是参数为λ的Poisson 过程。
刘次华随机过程 第四章马尔可夫(Markov)链
p
-1
0
1
i-1
i
i+1
一步转移概率:
pi,i+1 = p pi,i−1 = q = 1− p pii = 0
4.1 马尔可夫链与转移概率
k步转移概率:
i经过k步进入j,向右移了x步,向左移了y步
则
⎧x
⎨ ⎩
x
+ −
y y
= =
k j
−
i
⇒
⎧ ⎪⎪
x
⎨
⎪ ⎪⎩
y
= =
k k
+ −
(j 2 (j 2
定义4.1:若随机过程{Xn,n∈T },对任意 n∈T和i0, i1, …, in+1 ∈I,条件概率 P{Xn+1=in+1|X0=i0,X1=i1,…,Xn=in} = P{Xn+1=in+1|Xn=in}, 则称{Xn,n∈T }为马尔可夫链,简称马氏 链。
4.1 马尔可夫链与转移概率
定义4.2:称条件概率pij(n)= P{Xn+1=j|Xn=i} 为马尔可夫链{Xn,n∈T }在时刻n的一 步转移概率,简称转移概率,其中i,j∈I。 I j
)
=
⎧0 , i ≠ ⎩⎨1 , i =
j j
4.1 马尔可夫链与转移概率
定义4.5:
初始概率 绝对概率
p j = P{X 0 = j} pj( n ) = P{ Xn = j }
{ } 初始分布
pj , j∈I
{ } 绝对分布
pj (n) , j ∈ I
初始概率向量
pT (0) = ( p1, p2 , )
− i) − i)
应用随机过程4-更新过程
N (t ) k 1
X
k
, t 0
假设2
c (1 )
其中 0 称为 相 对安 全 负载 。
U (t ) , a.s. {ct S (t ), t 0}为齐次的独立增量过程。盈余过程 lim t 当盈余过程取负值时,称保险公司“破产”。T inf{t : U (t ) 0}
2010-9-2
定理4.3.2
(Blackwell更新定理)
记 E( X n ) ,
(1). 若 F 不是格点的,则对一切 a 0 ,当 t 时 a M (t a ) M (t )
(2). 若 F 是格点的,周期为 d,则当 n 时 d P{在nd处发生更新}
E[TN ( t ) 1 ] E[ X 1 X 2 X N ( t ) 1 ] E ( X 1 ) E ( N (t ) 1)
二* 、更新方程在人口学中的一个应用
设 B (t ) 为 t 时刻女婴的出生率,已知过去的 B (t ),t 0 ,要预测未 来的 B (t ),t 0 。
注: Feller初等更新定理是Blackwell更新定理的特殊情形。
2010-9-2
理学院 施三支
定理4.3.3
(关键更新定理)
记 E ( X n ) ,设函数 h (t ), t [0, ] ,满足
(1). h(t ) 非负不增;(2).
0
h(t ) dt 。 H (t ) 是更新方程
2010-9-2 理学院 施三支
例4.3.1
(剩余寿命与年龄的极限分布)
以 r (t ) TN ( t ) 1 t 表示时刻 t 的剩余寿命,即从 t 开始到下 次更新的时间,s (t ) t TN ( t ) 为 t 时刻的年龄。 求 r (t ) 和 s (t ) 的 极限分布。
应用随机过程-更新过程(PDF)
4.1 更新过程的定义及若干分布 4.2 更新方程及其应用 4.3 更新定理 4.4 Lundberger-Cramer 破产论 4.5 更新过程的推广
2010-9-2
理学院 施三支
4.1 更新过程的定义及若干分布
一、更新过程的定义
定义4.1.1 设{ X n , n 1,2, }是独立同分布的非负随机变量列,
t
其中f
(t)
m(t) F (t )。
f
(t) 0 m(t s) f
(s)ds
t
定义4.2.1 称积分方程 K (t) H (t) K (t s)dF (s) 为更新方程 0
其中 H (t),F (t) 为已知,且当 t 0时,H (t) F (t) 0。
当 H (t) 有上界时,称之为适定的。
假定1 {X k , k 1}是恒正的、独立同分布的随机变量列,
F (x) 是 X1 分布函数, 是 X1 的期望;{N (t), t 0} 是参数为 ,
且与{X k , k 1}独立的泊松过程。
2010-9-2
理学院 施三支
N (t)
到t时刻为止的索赔总额: S(t) X k ,t 0 k 1
则它表示初始盈余为u时,保险公司永不破产的概率,称为生存概率
2010-9-2
理学院 施三支
4.5 更新过程的推广
一、延迟更新过程
更新过程要求时间间隔是独立同分布的序列,如果放宽第一个时 间间隔X1,允许其分布不同,则由X1 ,X2 ,…确定的计数过程为 延迟更新过程。
二、更新回报过程
N (t)
设 R(t) Ri , 其 中 {N (t), t 0} 是 一 个 更 新 过 程 , i 1
4-4 更新过程
f
0 n 1
sn
(t)e stdt
φ(s) 1 φ(s)
f
n 1 0
sn
(t)e stdt
= φ S (s)
n 1
n
n 1
[φ x (s)]n
:
φ(s) Λ(s) 1-φ(s)
Λ(s) φ(s) 1 Λ(s)
更新过程的强度
λ(t)表示某时刻发生更新的强度:
P N ( t ) k P N ( t ) k
k 1 n 1 n 1 k n k
k 1
P N ( t ) n P S n t Fsn ( t )
n 1 n 1 n 1
更新过程的强度
该时间均匀分布于0到1小时之间
问:长时间工作情况下,更换电池的速率? 解 设第 i 次电池的使用时间是 xi,购买电池花费时间 ui
它们都是随机变量,且平均更新间隔为:
μ E x i E ui
电池更新的速率是
45 t dt 45.5
0
平均更新次数: m( t ) Fs ( t )
n 1
n
令λ(t)表示某时刻发生更新的强度:
即在 [t,t+dt) 里发生更新事件的平均次数
λ (t )dt dm (t )
λ(t)
f
n 1
sn
(t)dt
拉氏 Λ(s) 变换
λ (t)e st dt
0
故不等式左右均收敛于 E{x}
于是
t lim E{x} t N(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
为 避 免 因 可 能 的
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
且因为随机变量 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 ) 为更新函数。
t
N (t ) 的情况。 t
为考虑 N (t ) 的发散速度,我们先考虑到达时刻 TN (t ) ( TN (t ) 表示在时刻 t 或时刻 t 之前 最后一次更新发生的时刻,以此类推,则 TN (t )1 表示在时刻 t 之后第一次更新发生的时刻) 。 利用 TN (t ) 和 TN (t )1 ,我们提出并证明以下命题。
命题 4.3 当 t 时,以概率 1 保证
证明:因为 TN (t ) t TN (t ) 1 ,于是有
N (t ) 1 , ( EX n ) 。 t
[物理]第4章 更新过程
F H ( t ) F K (t )] k n
k 1
H (t ) m H (t )
三、瓦尔德(Wald)等式
1,停时:设{Xn,n≥1}为随机变量序列,N为取非负整数的随 机变量。若对一切的n=0,1,2, …,事件{N=n}仅依赖于X1, X2, …, Xn,而与Xn+1 Xn+2,…独立,则称N关于{Xn,n≥1}为停时(Stopping time),或称马尔可夫时(Markov time)。 直观意义:当我们观察诸Xn ,以N表示停止观察前所观察的次数, 如果N=n ,那么,我们是在已经观察X1, X2, …, Xn后,还未观察 Xn+1 Xn+2,…前停止观察。 2,Wald等式 Theorem:设{Xn,n≥1}为独立同分布随机变量序列 EX n N是关于{Xn,n≥1}的停时,
定理:设更新方程中H(t)为有界函数,则
方程存在惟一的在有限区间内有界的解
K (t ) H (t ) H (t s)dm(s)
0
t
证明:先证K(t)有界,因为H(t)有界, m(t) 有界不减,所以,对任何T>0有
0t T T
sup K (t ) sup H (t ) sup H (t s)dm( s)
n2
H (t ) F H (t ) ( ( Fn 1 F )) H (t )
n2
H (t ) F [ H (t ) ( Fn ) H (t )]
n 1
H (t ) F K (t ) H (t ) K (t s )dF ( s )
0 t
是方程 最后证明解的惟一性,设 K
- 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