南邮 随机过程课件 4
合集下载
随机过程课件打印版

当An An 1 , n 1
当An An 1 , n 1
9
A1 A2
连续性定理
A1 A2
则称P为(Ω,F)上的概率,(Ω,F,P)称 为概率空间,P(A)为事件A的概率。
An Ai 新事件:lim n i 1
lim An Ai
n i 1
3 对于R n中的任意区域, a1 , b1; a2 , b2 ;;a n , bn ,其中 ai bi , i 1,, n
F b1 , b2 ,, bn F b1 ,, bi 1 , ai , bi 1 ,bn F b1 ,, bi 1 , ai , bi 1 ,, b j 1 , a j , b j 1 ,, bn ,
d P({e : g( X ) y, e X }) dy
如果上式右端概率的导数对于y处处存在,那么这 个导数就给出了随机变量Y的概率密度
fY ( y)
19
20
n维联合分布函数F x1 , x 2 , x n 具有下列性质 :
三、边缘分布
若二维联合分布函数中有一个变元趋于无 穷,则其极限函数便是一维分布函数,对于这 种特殊性质,我们称其为边缘分布。 对于任意两个随机变量X,Y,其联合分布函数为: F ( x, y ) 则: FX ( x ) P ( X x ) P ( X x , Y ) F ( x , )
P( X x,Y y) P((X x) (Y y)) P( X x)P(Y y)kFra biblioteknpkq
nk
, k 0 ,1 , 2 n
p
P(X k)
k
k!
当An An 1 , n 1
9
A1 A2
连续性定理
A1 A2
则称P为(Ω,F)上的概率,(Ω,F,P)称 为概率空间,P(A)为事件A的概率。
An Ai 新事件:lim n i 1
lim An Ai
n i 1
3 对于R n中的任意区域, a1 , b1; a2 , b2 ;;a n , bn ,其中 ai bi , i 1,, n
F b1 , b2 ,, bn F b1 ,, bi 1 , ai , bi 1 ,bn F b1 ,, bi 1 , ai , bi 1 ,, b j 1 , a j , b j 1 ,, bn ,
d P({e : g( X ) y, e X }) dy
如果上式右端概率的导数对于y处处存在,那么这 个导数就给出了随机变量Y的概率密度
fY ( y)
19
20
n维联合分布函数F x1 , x 2 , x n 具有下列性质 :
三、边缘分布
若二维联合分布函数中有一个变元趋于无 穷,则其极限函数便是一维分布函数,对于这 种特殊性质,我们称其为边缘分布。 对于任意两个随机变量X,Y,其联合分布函数为: F ( x, y ) 则: FX ( x ) P ( X x ) P ( X x , Y ) F ( x , )
P( X x,Y y) P((X x) (Y y)) P( X x)P(Y y)kFra biblioteknpkq
nk
, k 0 ,1 , 2 n
p
P(X k)
k
k!
随机过程精品课件 (4)

续EX.3 醉汉问题 酒吧
1 2 3 4 5
家
醉汉在街上徘徊, 在每一个街口以1/3的概 率停下, 以1/3的概率向前或向后. 若他又返回酒吧或到家门, 不再游动 . 状态转移图为
湖南大学
1 1 1/3
1/3 2
1/3 1/3
1/3 3
1/3 1/3
1/3 4
1/3
1 5
分析状态“ 2”的类型很困难 .
1 ]2 [ 1 2 3 2 1, 3 3 1 1 3 湖南大学
1 n 2( ) , 3
n 2,3,
状态3是非常返的.
1 因 , d 3 1, 3 由于 2 3, 3 4,
(1) p3的,且 E={1}∪{5}∪{2, 3, 4}. 关于常返状态有下结论: 定理4.3.6 设i , j∈E, 且i ≠j,则
湖南大学
1)若状态i 和j 互通, i 是常返态,则j 也是 常返态; 2)若状态i 和j 互通, 且 j 是常返态, 则 f ij 1; 3)若状态i 是常返的, 且 i j , 则 j i .
( i j )
四、闭集、状态空间的分解
一步转 移概率
定义4.3.9 设E 是状态空间, C E , 若对
其中H=C(j ).
证明 f kj
1, k H ; 0, k H 且 k N .
f ij P {经有限步到达 j X ( 0) i }
k E
P {经有限步到达 j X (0) i , X (1) k }
P { X (1) k X ( 0 ) i }
湖南大学
定理4.3.8 分解定理 齐次马氏链的状态空间可唯一地分解 成有限个或可列多个不相交的状态子集之 并. E=N∪C1∪C2∪… 其中 1)N是所有非常返态所成之集; 2)每个Cn,(n=1,2,…)均为常返状态 组成的不可约闭集.
南邮系统工程课件第四章

3
1. 数学模型建模程序
含义: 将现实世界中的原型概括形成数学模型的过程。 四个环节:
数 学 模 型 数 学 结 论
现实世 界的分 析、预 测、决 策或控 制
现实 世界 的原 型
翻译Biblioteka 实验分析抽象比较检验
4
(1)变量类型
领域 计量经济 控制理论 类型1
对模型起作用但 是不属于模型描 述范围的因素
模型所研究的 因素 类型2 内生变量 输出变量
试计算公司在第t年的总价值。
31
请用矩阵形式表示状态方程与输出方程。
x1 (t ) : 第t年交换设备的全部价值 x2 (t ) : 第t年电缆的全部价值 y (t ):第t年公司的总价值,y (t ) x1 (t ) x2 (t )
x1 (t 1) 0.8x1 (t ) 0.15x2 (t ) 0.75u(t ) x2 (t 1) x2 (t ) 0.25u(t )
有向图 权重有向图 邻接矩阵 可达矩阵
7
(1)能源需求权重有向图模型
1
-1.8
2
2.6
1.2
0.7
-0.4
1 能源供给 3 能源需求
2 价格 4 人口
3
5.8
4
有向边上数字代表起始节点变化单位值对于相邻 节点的影响 1和2节点之间的权重为-1.8,表示能源供给增加 1个单位,价格下跌1.8个单位
28
29
(2)系统的状态空间描述
系统工程中所涉及的系统分为:连续系统和离 散系统两类。
连续系统的系统方程
X (t ) f ( X (t ), u (t ), t ) Y (t ) g ( X (t ), u (t ), t )
《随机过程——计算与应用》课件随机过程引论课件4

称 Z (t) E Zt 2
为复随机过程Z的均方值函数.
对任意的s,t T , 称 RZ (s,t) E[ZsZt ]
为复随机过程Z的相关函数.
称 CZ (s,t) Cov(Zs , Zt ) E[(Zs mZ (s))(Zt mZ (t))]
为复随机过程Z的协方差函数.
由以上定义可得
则称之为该二维随机过程的互相关函数,记 RXY (s, t).
若 Cov( X s ,Yt ) E[(X s mX (s))(Yt mY (t))] 存在
则称之为该二维随机过程的互协方差函数,记 CXY (s, t).
互协方差函数可以定义两个过程的相关性 设有随机过程X={Xt, t∈T}和Y={Yt, t∈T}, 对任意 的 s,t∈T, 若
则称之为随机过程X的协方差函数.记为CX (s, t).
4. 相关函数
设 X {Xt , t T} 是一实值随机过程,对任意s,t∈T,
若
E[X s Xt ] 存在
则称之为随机过程X的(自)相关函数.记为 RX (s,t).
5. 均方值函数
设 X {Xt , t T} 是一实值随机过程,对任意3.10 设 X {Xt ,t [a,b]} 是正交增量过程, 且Xa 0 则
(1) RX (s,t) X (min( s,t)) s,t [a,b]
(2) X(t )是单调不减函数
两个随机过程的互相关函数与互协方差函数
设{Xt ,Yt ,t T}为二维随机过程,对任意s,t T,若 E[XsYt ] 存在
2
s,t
CX (s,t) RX (s,t) mX (s)mX (t)
RX (s,t)
a2 cos(t s),
随机过程课件PPT资料(正式版)

应怎样分才合理呢➢?」
☞随机事件:样本空间的子集,常记为 A ,B ,…它是满足某些条件的样本点所组成的集合.
排队和服务系统 ◙A∩勤B 奋⇔、A刻B :苦A、与合➢B作的、积探事索件;; 更新过程 为从事科学研究打下坚实的基础;
☞抽取的是精装中➢文版数学书 ⇒
➢ 时间序列分析
➢ 鞅过程
绪论
《随机过程》基础
概率(或然率或几率) ——随机事件出现的可能 性的量度;
概率论其起源与博弈、 、天气预报等问题有 关
⊕16世纪意大利学者开始研究掷骰子等赌博 中的一些问题;
⊕17世纪中叶,「现有两个赌徒相约赌若干 局,谁先赢S局就算赢了,当赌徒A赢K局(K<S), 而赌徒B赢L局(L<S)时,赌博中止,赌资应怎 样分才合理呢?」
随机过程课件
《随机过程》
➢ 教材: ◙ 张卓奎,陈慧婵,随机过程.西安电子科技大 学.2003.
➢ 主要参考文献: ◙ 胡奇英编著,随机过程.西安电子科技大学.1998. ◙ 周荫清 ,随机过程习题集. 清华大学出版社, 2004. ◙ 林元纟金烈 ,应用随机过程. 清华大学出版社, 2002.
……
➢ 随机过程理论在社会科学中例如在社会统计, 学、经 济、金融工程、管理中也得到极其广泛的应用。
➢ 为从事科学研究打下坚实的基础;
绪论
教学目标
➢ 充分理解、熟练掌握教材的内容 ◙ 熟练掌握基本的数学概念和定理;
◙ 熟练掌握随机过程研究对象的数学描述;
Hale Waihona Puke ➢ 通过学习和练习,具备一定的分析、解决本专业具体 问题的能力;
☞拉普拉斯曾说:“生活中最重要的问题,其中 绝大多数在实质上只是概率的问题”。
☞概率论是研究随机现象数量规律的数学分支。 在实际中,人们往往还需要研究在时间推进中某 一特定随机现象的演变情况,描述这种演变的就 是概率论中的随机过程。
☞随机事件:样本空间的子集,常记为 A ,B ,…它是满足某些条件的样本点所组成的集合.
排队和服务系统 ◙A∩勤B 奋⇔、A刻B :苦A、与合➢B作的、积探事索件;; 更新过程 为从事科学研究打下坚实的基础;
☞抽取的是精装中➢文版数学书 ⇒
➢ 时间序列分析
➢ 鞅过程
绪论
《随机过程》基础
概率(或然率或几率) ——随机事件出现的可能 性的量度;
概率论其起源与博弈、 、天气预报等问题有 关
⊕16世纪意大利学者开始研究掷骰子等赌博 中的一些问题;
⊕17世纪中叶,「现有两个赌徒相约赌若干 局,谁先赢S局就算赢了,当赌徒A赢K局(K<S), 而赌徒B赢L局(L<S)时,赌博中止,赌资应怎 样分才合理呢?」
随机过程课件
《随机过程》
➢ 教材: ◙ 张卓奎,陈慧婵,随机过程.西安电子科技大 学.2003.
➢ 主要参考文献: ◙ 胡奇英编著,随机过程.西安电子科技大学.1998. ◙ 周荫清 ,随机过程习题集. 清华大学出版社, 2004. ◙ 林元纟金烈 ,应用随机过程. 清华大学出版社, 2002.
……
➢ 随机过程理论在社会科学中例如在社会统计, 学、经 济、金融工程、管理中也得到极其广泛的应用。
➢ 为从事科学研究打下坚实的基础;
绪论
教学目标
➢ 充分理解、熟练掌握教材的内容 ◙ 熟练掌握基本的数学概念和定理;
◙ 熟练掌握随机过程研究对象的数学描述;
Hale Waihona Puke ➢ 通过学习和练习,具备一定的分析、解决本专业具体 问题的能力;
☞拉普拉斯曾说:“生活中最重要的问题,其中 绝大多数在实质上只是概率的问题”。
☞概率论是研究随机现象数量规律的数学分支。 在实际中,人们往往还需要研究在时间推进中某 一特定随机现象的演变情况,描述这种演变的就 是概率论中的随机过程。
随机过程4(34)精品PPT课件

1. 功率谱密度的概念
●工程实际中, 能量有限的信号x(t)称为 能量型信号, 可以定义它的总能量:
x2 (t)dt
当时间趋于无穷时,它的平均功率趋于零.
●另一类信号x(t),其能量是无限的,但平均功率有限.即
P lim 1 T x2(t)dt
T 2T T
称为 功率型信号.周期信号就是常见的功率信号.
说明信号的总能量等于能谱密度在全频域上的积分. 右式也是总能量的谱表达式.
由于实际中很多信号(函数)的总能量是无限的, 不满足绝对可积的条件,所以通常研究x(t)在 (-∞,+ ∞)上的平均功率,即
lim 1 T x2 (t)dt
T 2T T
为了能利用Fouier变换给出平均功率的谱表达式, 构造一个截尾函数:
§4 平稳过程的功率谱密度
● 之前对平稳过程的讨论都是在进行的. 在时域上描述了平稳过程的统计特征.
● 但对许多物理和工程领域中问题,不仅要研究其 在时域上的特性,还要研究其在频域内的特征,即 从频率的角度来研究随机过程的统计特征. 例如对信号处理、线性系统分析以及随机振动的 研究. 其中广泛采用的方法是频率域分析方法.
T e jt X (t)dt
T
称 lim E[ 1 T X 2 (t) dt] 为过程的平均功率.
T 2T T
定理1 设{X(t), -∞<t<+ ∞}是平稳过程,若RX(τ) 绝对可积,则有
S X
()
lim
T t
T
lim 1
T 2T
E[ FX (,T ) 2 ]=
即
S(x )
lim 1 T 2T
Fx (,T )
2
lim
随机过程 第4章 Markov过程

(C-K 方程)
证明:由全概率公式,有:
( m+r ) pij (n) = P{ X n + m + r = j X n = i}
= ∑ P{ X n + m + r = j, X n + m = k X n = i} = ∑ P{ X n + m + r = j X n + m = k , X n = i} ⋅P{ X n + m = k X n = i} = ∑ P{ X n + m + r = j X n + m = k}P{ X n + m = k X n = i}
为齐次马氏链的 m 步转移(概率)矩阵。 显然有:
(m) pij (n) ≥ 0, i , j ∈ I
∑p
j∈I
(m) ij
( n) = 1 , i ∈ I
m = 1 时,即为一步转移矩阵。
规定:
⎧1, (0) pij (n) = δ ij = ⎨ ⎩0,
二、切普曼-柯尔莫哥洛夫(C-K)方程
i= j i≠ j
由01nll相互独立111100nnnnpxixixixi??l111100111100nnnnnnpxixixixipxixixi????ll第四章markov过程611121100011211000?nnnnnnnnnpiiiiiiipiiiii?????????l??l1nnnpii??11??nnnnixixp故210lnxn满足markov性且1100nnijnnkkkkppxjxipji??10nnkknji?pji?ipji?q?二随机游动1无限制的随机游动
性质 5 设{ X n , n ≥ 0 }为马氏链, 其状态空间为 I, 则对任意给定的 n 个整数,
随机过程4-3

pi(0) pii1 (n1 ) pi1i2 (n2 n1 ) pi2i3 (n3 n2 ) pim1im (nm nm 1 )
i
遍历性
lim pij ( n) p j , i , j E
n
(2.15)
如果 { p j , j 1, 2,...} 满足 率的极限分布。
第17页
二、柯尔莫哥洛夫向前和向后方程
设 { X (t ), t [0, )} 是状态有限(即具有有限多个状态) 的马尔科夫过程,E {0,1, 2,..., N }
定义 设状态有限的马尔科夫过程 X(t) 的转移概率函数 为 pij (t ) ,若 1, i j (3.10) lim pij (t ) ij t 0 0, i j
个时刻 t1 , t2 ,..., tm (0 t1 t 2 ... tm ) ,任意正数 s 以及任意 i1 , i2 ,..., im , j E ,满足
第4章 马尔科夫过程
第11页
P{ X (tm s) j | X (t1 ) i1 , X (t2 ) i2 ,...., X (tm ) im } P{ X (tm s) j | X (tm ) im }
以后我们只讨论时齐马尔科夫过程。
(3.2)
根据条件概率性质。转移概率函数具有下列两条性质:
(有限多个或无限多个) (1) 0 pij ( s) 1, i , j 0,1, 2,...
(2)
p ( s) 1, i 0,1, 2, ...
ij j
第4章 马尔科夫过程 通常,我们规定
则称 { X (t ), t [0, )}为马尔科夫过程。 (3.1)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
7
此乃有名的切普曼 − 柯尔莫哥洛夫方程,
8 9
∑ p (n ) = 1.
ij
证明:利用概率公式及马尔可夫性有: pij (n ) = P{X m+ n = a j / X m = ai } = P{X m = ai , X m+ n = a j } P{X m = ai }
= =
ar ∈I
∑ P{X
ir
m+ k
= ar / X m = ai }⋅ P{X m+ n = a j / X m+ k = ar }
rj
5、初始概率与绝对概率
用矩阵形式表示为: P(n ) = P(k ) ⋅ P(n − k )
ar ∈I
∑ p (k ) ⋅ p (n − k )
(1) 定义:设{X n , n ≥ 0}为马尔可夫链, 分别称
n −1in
对任意的n ≥ 1, 和整数i0 , i1 L, in
n
(2) P{X 1 = ai , X 2 = ai ,L, X n = ai / X 0 = ai } = pi i L pi
变量序列,且ξ 0 , ξ1 , ξ 2 ,L有相同的分布,令X n= ∑ ξ k
n −1in
P{X n = in / X 0 = i0 , X 1 = i1 L, X n−1 = in−1} P{X 0 = i0 , X 1 = i1 L, X n = in } P{X 0 = i0 , X 1 = i1 L, X n−1 = in−1}
由所有的一步转移概率pij构成的矩阵
P{ X m+ k = a j / X m = ai } = pij (m, m + k )
∑p
a j∈ I
ij
=1
为马氏链在m时刻处于ai 状态经k步,在m + k时刻 转移到a j 状态的转移概率,记为:pij (m, m + k )
i, j , m, k均为正整数, 一般pij (m, m + k )与i, j , m, k有关, 若pij (m, m + k )与m无关,则称马氏链为齐次的, 下面我们仅讨论齐次马氏链。
ij
ar
aj
ai
m
m+k m+n
整数n ≥ 0, ai , a j ∈ I, 有:pij (n ) =
或 P(n ) = P(k )P(n − k ) 简称c − k方程。
ar ∈I
∑p
ir
(k ) prj (n − k )
直观解释对照图
(1) ( 2)
pij (n ) ≥ 0.
a j ∈I
ai , a j ∈ I ai ∈ I
ai ∈I ai ∈I ai ∈I
(3)马氏链的有限维分布 定理:设{X n , n ≥ 0}为齐次马氏链, 则对任意的
a j ∈ I和n ≥ 1,绝对概率p j (n)具有性质: r r (1) p (n) = ∑ pi (0) pij (n) 或 p (n) = p (0) P(n)
j
ai , ai ,L, ai ∈ I和n ≥ 1有:
17
k =0 n
=
= P{ξ n = in − in −1}.
P{ξ 0 = i0 , ξ1 = i1 − i0 L, ξ n = in − in −1} P{ξ 0 = i0 , ξ1 = i1 − i0 ,L, ξ n −1 = in −1 − in − 2 }
同理:P{X n = in / X n −1 = in −1} = P{ξ n =i n −in −1}
2
3
2、马氏链的转移概率
ai
称条件概率
aj
m+k
3、一步转移概率及矩阵
在上面转移概率中,取 k = 1即得一步转移概率
pij具有性质: (1) ( 2) pij ≥ 0 ai , a j ∈ I ai ∈ I
m
pij = pij (m, m + 1) = P{ X m+1 = a j / X m = ai }
1 n 1 n −1
14
= ∑ pi (0) pii pi i L pi
ai ∈I
1 12
n −1in
15
推论: (1) P{X 0 = ai , X 1 = ai ,L, X n = ai } = pi (0) pi i L pi
0 1 n 0 01 1 2 n 0 01
例:随机游动,设ξ 0 , ξ1 , ξ 2 ,L是整数值独立随机
11
为马氏链的初始概率和绝对概率,并分别称
{p (0), a
j
= a ∈I
r
∑ P{X
P{X m = ai }
分布和绝对分布,简记为和{p j (0)}和{p j (n)} 。
写成向量形式: r p(0) = ( p1 (0), p2 (0),L, p j (0),L) r p(n) = ( p1 (n), p2 (n),L, p j (n),L)
而X (t )在每一时刻t (n = 0,1,2,L), 所处状态记为:
= P{X m + k = ai
P{X m+ k = ai
m+k
/ X m = ai , X m−1 = ai ,L, X 1 = ai }
m m+k
/ X m = ai
m
}
m −1
1
则称此随机过程{X n , n ≥ 0}为马尔可夫链, 简称马氏链。
则称{X n , n ≥ 0}为随机游动。
k =0
=
总结: 1)齐次马氏链多步转移概率可由一步转移概率确定; P ( n) = P n 2) 绝对概率可由初始概率及n步转移概率确定
p j (n) = ∑ pi (0) pij (n)
ai ∈I
随机游动可以解释为质点在直线上的整数格点上作 运动的质点,初始位置为X 0 = ξ 0,每隔一个单位时 间质点移动一次,第k次移动的长度为整数ξ k,于是 X n = ∑ ξ k 表示在时刻n质点的位置,则{X n , n ≥ 0}是 随机游动,随机游动{X n , n ≥ 0}是时齐马氏链。
10
在上式取k = 1,P (n) = P ⋅ P (n − 1) 则n = 2时有:P (2) = P (1) P (1) = P 2 n = 3时:P (3) = P (1) P (2) = P 3 一般当n为任意整数时有:P (n) = P n
表明一步转移概率是最基本的,它确定了 马氏链的状态转移的统计规律。
3)有限维分布可完全由初始概率及一步转移概率确定。
16
一步转移概率为: pij = P{X n = j / X n−1 = i} = P{ξ n = j − i} = p j −i
18
几种特殊的随机游动 例.无限制的随机游动:质点在直线上作随机游动,
如某一时刻质点位于i,则下一步质点以概率p向 右移动一格到达i + 1, 或以概率q = 1 − p向左移一 格到达i − 1, 若以X n 表示时刻n时质点的位置,则
j
m
12
(2) 绝对概率与初始概率的关系
定理:设{X n , n ≥ 0}为马尔可夫链,则对任意的
证: (1) p j (n) = P{X n = a j } = ∑ P{X 0 = ai , X n = a j } = ∑ P{X 0 = ai }⋅ P{X n = a j / X 0 = ai } = ∑ pi (0) pij (n)
11 12 1n 21 2:
1 i = j p (0 ) = pij (m, m ) = δ ij = 0 i ≠ j (2)、切普曼—柯尔莫哥洛夫方程 (Chapman-Kolmogorov) 定理:设{X n , n ≥ 0}为齐次马氏链, 则对任意的
pi ,i +1 = p 一步转移概率为: pi ,i −1 = q p =0 ii
(i ∈ I ,0 ≤ p ≤ 1)
j ≠ i + 1, i − 1, j ∈ I
由于m1 , m2只能取整数,所以n + ( j − i )必须是 偶数,且在 n 步中哪 m1 步向右,哪m2步向左 是任意的,选取方法为:Cnm
n
则X n = ∑ξn为随机游动,它是一个 齐次马氏链
ai ∈I ai ∈I
{
1
2
n
}
n
= ∑ P{X 0 = ai }⋅ P{ X 1 = ai / X 0 = ai } ⋅ P{ X 2 = ai /
ai ∈I
1 2
n −1 n
13
= ∑ Pi (n − 1) pij
ai ∈I
X 0 = ai , X 1 = ai }L P{ X n = ai / X 0 = ai , X 1 = ai L X n −1 = ai }
20
{X n , n ≥ 0}是一随机过程。 1 , 第k次向右移动一格 令ξk = −1,第k次向左移动一格
n k =0
m1 m1 m2 Cn p q , n + ( j − i )为偶数 ∴ pij (n) = n + ( j − i )为奇数 0 n n Cn2 p 2 q 2 , n为偶数 pii (n) = n为奇数 0
p j (0) = P{X 0 = a j }和p j (n) = P{X n = a j } , (a j ∈ I ) ∈ I }和{p j (n), a j ∈ I }为马氏链的初始
= a ∈I
r
∑ P{X
m
= ai , X m+ k = ar , X m+ n = a j } = ai , X m+ k = ar }P{X m+ n = a j / X m = ai , X m+ k = ar } P{X m = ai }
此乃有名的切普曼 − 柯尔莫哥洛夫方程,
8 9
∑ p (n ) = 1.
ij
证明:利用概率公式及马尔可夫性有: pij (n ) = P{X m+ n = a j / X m = ai } = P{X m = ai , X m+ n = a j } P{X m = ai }
= =
ar ∈I
∑ P{X
ir
m+ k
= ar / X m = ai }⋅ P{X m+ n = a j / X m+ k = ar }
rj
5、初始概率与绝对概率
用矩阵形式表示为: P(n ) = P(k ) ⋅ P(n − k )
ar ∈I
∑ p (k ) ⋅ p (n − k )
(1) 定义:设{X n , n ≥ 0}为马尔可夫链, 分别称
n −1in
对任意的n ≥ 1, 和整数i0 , i1 L, in
n
(2) P{X 1 = ai , X 2 = ai ,L, X n = ai / X 0 = ai } = pi i L pi
变量序列,且ξ 0 , ξ1 , ξ 2 ,L有相同的分布,令X n= ∑ ξ k
n −1in
P{X n = in / X 0 = i0 , X 1 = i1 L, X n−1 = in−1} P{X 0 = i0 , X 1 = i1 L, X n = in } P{X 0 = i0 , X 1 = i1 L, X n−1 = in−1}
由所有的一步转移概率pij构成的矩阵
P{ X m+ k = a j / X m = ai } = pij (m, m + k )
∑p
a j∈ I
ij
=1
为马氏链在m时刻处于ai 状态经k步,在m + k时刻 转移到a j 状态的转移概率,记为:pij (m, m + k )
i, j , m, k均为正整数, 一般pij (m, m + k )与i, j , m, k有关, 若pij (m, m + k )与m无关,则称马氏链为齐次的, 下面我们仅讨论齐次马氏链。
ij
ar
aj
ai
m
m+k m+n
整数n ≥ 0, ai , a j ∈ I, 有:pij (n ) =
或 P(n ) = P(k )P(n − k ) 简称c − k方程。
ar ∈I
∑p
ir
(k ) prj (n − k )
直观解释对照图
(1) ( 2)
pij (n ) ≥ 0.
a j ∈I
ai , a j ∈ I ai ∈ I
ai ∈I ai ∈I ai ∈I
(3)马氏链的有限维分布 定理:设{X n , n ≥ 0}为齐次马氏链, 则对任意的
a j ∈ I和n ≥ 1,绝对概率p j (n)具有性质: r r (1) p (n) = ∑ pi (0) pij (n) 或 p (n) = p (0) P(n)
j
ai , ai ,L, ai ∈ I和n ≥ 1有:
17
k =0 n
=
= P{ξ n = in − in −1}.
P{ξ 0 = i0 , ξ1 = i1 − i0 L, ξ n = in − in −1} P{ξ 0 = i0 , ξ1 = i1 − i0 ,L, ξ n −1 = in −1 − in − 2 }
同理:P{X n = in / X n −1 = in −1} = P{ξ n =i n −in −1}
2
3
2、马氏链的转移概率
ai
称条件概率
aj
m+k
3、一步转移概率及矩阵
在上面转移概率中,取 k = 1即得一步转移概率
pij具有性质: (1) ( 2) pij ≥ 0 ai , a j ∈ I ai ∈ I
m
pij = pij (m, m + 1) = P{ X m+1 = a j / X m = ai }
1 n 1 n −1
14
= ∑ pi (0) pii pi i L pi
ai ∈I
1 12
n −1in
15
推论: (1) P{X 0 = ai , X 1 = ai ,L, X n = ai } = pi (0) pi i L pi
0 1 n 0 01 1 2 n 0 01
例:随机游动,设ξ 0 , ξ1 , ξ 2 ,L是整数值独立随机
11
为马氏链的初始概率和绝对概率,并分别称
{p (0), a
j
= a ∈I
r
∑ P{X
P{X m = ai }
分布和绝对分布,简记为和{p j (0)}和{p j (n)} 。
写成向量形式: r p(0) = ( p1 (0), p2 (0),L, p j (0),L) r p(n) = ( p1 (n), p2 (n),L, p j (n),L)
而X (t )在每一时刻t (n = 0,1,2,L), 所处状态记为:
= P{X m + k = ai
P{X m+ k = ai
m+k
/ X m = ai , X m−1 = ai ,L, X 1 = ai }
m m+k
/ X m = ai
m
}
m −1
1
则称此随机过程{X n , n ≥ 0}为马尔可夫链, 简称马氏链。
则称{X n , n ≥ 0}为随机游动。
k =0
=
总结: 1)齐次马氏链多步转移概率可由一步转移概率确定; P ( n) = P n 2) 绝对概率可由初始概率及n步转移概率确定
p j (n) = ∑ pi (0) pij (n)
ai ∈I
随机游动可以解释为质点在直线上的整数格点上作 运动的质点,初始位置为X 0 = ξ 0,每隔一个单位时 间质点移动一次,第k次移动的长度为整数ξ k,于是 X n = ∑ ξ k 表示在时刻n质点的位置,则{X n , n ≥ 0}是 随机游动,随机游动{X n , n ≥ 0}是时齐马氏链。
10
在上式取k = 1,P (n) = P ⋅ P (n − 1) 则n = 2时有:P (2) = P (1) P (1) = P 2 n = 3时:P (3) = P (1) P (2) = P 3 一般当n为任意整数时有:P (n) = P n
表明一步转移概率是最基本的,它确定了 马氏链的状态转移的统计规律。
3)有限维分布可完全由初始概率及一步转移概率确定。
16
一步转移概率为: pij = P{X n = j / X n−1 = i} = P{ξ n = j − i} = p j −i
18
几种特殊的随机游动 例.无限制的随机游动:质点在直线上作随机游动,
如某一时刻质点位于i,则下一步质点以概率p向 右移动一格到达i + 1, 或以概率q = 1 − p向左移一 格到达i − 1, 若以X n 表示时刻n时质点的位置,则
j
m
12
(2) 绝对概率与初始概率的关系
定理:设{X n , n ≥ 0}为马尔可夫链,则对任意的
证: (1) p j (n) = P{X n = a j } = ∑ P{X 0 = ai , X n = a j } = ∑ P{X 0 = ai }⋅ P{X n = a j / X 0 = ai } = ∑ pi (0) pij (n)
11 12 1n 21 2:
1 i = j p (0 ) = pij (m, m ) = δ ij = 0 i ≠ j (2)、切普曼—柯尔莫哥洛夫方程 (Chapman-Kolmogorov) 定理:设{X n , n ≥ 0}为齐次马氏链, 则对任意的
pi ,i +1 = p 一步转移概率为: pi ,i −1 = q p =0 ii
(i ∈ I ,0 ≤ p ≤ 1)
j ≠ i + 1, i − 1, j ∈ I
由于m1 , m2只能取整数,所以n + ( j − i )必须是 偶数,且在 n 步中哪 m1 步向右,哪m2步向左 是任意的,选取方法为:Cnm
n
则X n = ∑ξn为随机游动,它是一个 齐次马氏链
ai ∈I ai ∈I
{
1
2
n
}
n
= ∑ P{X 0 = ai }⋅ P{ X 1 = ai / X 0 = ai } ⋅ P{ X 2 = ai /
ai ∈I
1 2
n −1 n
13
= ∑ Pi (n − 1) pij
ai ∈I
X 0 = ai , X 1 = ai }L P{ X n = ai / X 0 = ai , X 1 = ai L X n −1 = ai }
20
{X n , n ≥ 0}是一随机过程。 1 , 第k次向右移动一格 令ξk = −1,第k次向左移动一格
n k =0
m1 m1 m2 Cn p q , n + ( j − i )为偶数 ∴ pij (n) = n + ( j − i )为奇数 0 n n Cn2 p 2 q 2 , n为偶数 pii (n) = n为奇数 0
p j (0) = P{X 0 = a j }和p j (n) = P{X n = a j } , (a j ∈ I ) ∈ I }和{p j (n), a j ∈ I }为马氏链的初始
= a ∈I
r
∑ P{X
m
= ai , X m+ k = ar , X m+ n = a j } = ai , X m+ k = ar }P{X m+ n = a j / X m = ai , X m+ k = ar } P{X m = ai }