随机过程 第4章 马尔可夫链
- 格式:pdf
- 大小:290.11 KB
- 文档页数:46
马尔可夫链马尔可夫链是一种特殊的随机过程,最初由A.A .M arkov 所研究。
它的直观背景如下:设有一随机运动的系统E (例如运动着的质点等),它可能处的状态记为,....E ,...,E ,E n 10总共有可数个或者有穷个。
这系统只可能在时刻t=1,2,…n,…上改变它的状态。
随着∑的运动进程,定义一列随机变量Xn,n=0,1, 2, ⋯其中Xn=k ,如在t=n 时,∑位于Ek 。
定义1.1 设有随机过程}{T n X n ∈,,若对任意的整数T n ∈和任意的,,...,110I i i i n ∈+条件概率满足}i {},...,i X i {1n 10001n 1n n n n n n i X X P i X X P ======++++ 则称}{T n X n ∈,为马尔可夫链,简称为马氏链。
实际中常常碰到具有下列性质的运动系统∑。
如果己知它在t=n 时的状态,则关于它在n 时以前所处的状态的补充知识,对预言∑在n 时以后所处的状态,不起任何作用。
或者说,在己知的“现在”的条件下, “将来”与“过去”是无关的。
这种性质,就是直观意义上的“马尔可夫性”,或者称为“无后效性”。
假设马尔可夫过程}{T n X n ∈,的参数集T 是离散时间集合,即T={0,1,2,…},其相应Xn 可能取值的全体组成的状态空间是离散状态空间I={1,2,..}。
定义1.2 条件概率}{P 1)(i X j X p n n n ij ===+称为马尔可夫链}{T n X n ∈,在时刻n 的一步转移矩阵,其中i ,j ∈I ,简称为转移概率。
一般地,转移概率)(P n ij 不仅与状态i,j 有关,而且与时刻n 有关。
当)(P n ij 不依赖于时刻n 时,表示马尔可夫链具有平稳转移概率。
若对任意的i ,j ∈I ,马尔可夫链Xn,n ∈T}的转移概率)(P n ij 与n 无关,则称马尔可夫链是齐次的。
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 } 为此马氏链的平稳分布。
0.5丿 当初始分布为P{ X 0 = 1} =P{X 0 =2} = 0, P{ X 0 = 3} = 1时经三步转移后处于状态 3的概率。
7 .已知本月销售状态的初始分布和转移概率矩阵如下:1•设质点在区间[0 , 4]的整数点作随机游动,到达 0点或4点后以概率1停留在原处, 1 —向左、右移动一格或停留在原处。
求质点随机游动的一 3在其它整数点分别以概率 步和二步转移的概率矩阵。
2.独立地重复抛掷一枚硬币, 1, 2或3,这些值分别对应于第 n -1次和第n 次抛掷的结果为(正,正),(正,反), (反,正)或(反,反)。
求马尔可夫链{X n ,n 0,1,2,…}的一步和二步转移的概 率矩阵。
设{X n , n _0}为马尔可夫链,试证: (1 ) P{X n.1=i n1,X n.2=i n.2, ,X n^ ~lnm |X 0 - i 0,X ^i 1, ,X n=i n }= P{X n ・1 =in1,X n 2 - i n 2 , , X n m - i n m | X n - i n }(2) P{X 0 =i°,X 1 , X n - i n , Xn 2 ~ i n 2 , , X n ~ i n m | Xn ~ i n 1}= P{X ° = i°,X 1 二「…,X n -i n |X n^^i n-1} P{X n-2 ~ i n 2 / , Xn m i n m | Xn 1 _ i n 1}设{X n , n _1}为有限齐次马尔可夫链,其初始分布和转移概率矩阵为 每次抛掷出现正面的概率为 p ,对于n 一 2求,令X n =0, 3. 4. P i 二 P{X 。
5. P{X 2=4|X 设{X(t),r T}为随机过程 立同分布随机变量序列,令 {Y n , n _0}是马尔可夫链。
1/4 1/4 1/4 1/4"1/4 1/4 1/4 1/4 1/4 1/8 1/4 3/8J/4 1/4 1/4 1/4』0=1, 1 <X 1<4^ P{X ,且 X 1 =X(t 1),X 2,试证 1 「4"3,4,八 2 = 4 |1 :: X r :: 4}= X(t 2),…,X n = X(tJ …为独 Y 0 -0,Y ^-Y(t 1W X 1,Y ncY n 4^X n, n 一2,试证0.5 0.56.已知随机游动的转移概率矩阵为0.5 0.5 ,求三步转移概率矩阵 P (3)及0.5(1) P T(O) =(0.4, 02 0.4), P 二0.80.80.1 0.10.70.2 020.20.60.7 0.1 0.1 0.1?0.1 0.6 0.2 0.1(2) P T(0)=(02 02 0.3, 0.3) , p =0.1 0.1 0.6 0.230.1 0.2 0.5」求下一、二个月的销售状态分布。
第四章 马尔可夫链随机过程在不同时刻下的状态之间一般具有某种关系,马尔可夫(Markov )过程就是描述一类状态之间具有某种特殊统计联系的随机过程.Markov 过程在近代物理学、生物学、管理科学、信息处理与数字计算方法等领域都有重要的应用.按其状态和时间参数是连续的或离散的,它可分为三类:(1)时间、状态都是离散的Markov 过程,称为Markov 链;(2)时间连续、状态离散的Markov 过程,称为连续时间的Markov 链;(3)时间、状态都连续的Markov 过程.本章主要讨论Markov 链,有关连续时间的Markov 链的相关理论将在下章讨论.4.1 马尔可夫链的概念和例子独立随机试验模型最直接的推广就是Markov 链模型,早在1906年俄国数学家Markov 对它进行研究而得名,以后Kolmogorov 、Feller 、Doob 等数学家发展了这一理论.4.1 .1 Markov 链的定义假设Markov 过程{,}n X n T ∈的参数集T 是离散时间集合,即{0,1,2,}T =,相应n X 可能取值的全体组成的状态空间是离散状态集012{,,,}I i i i =.定义 4.1 设有一随机过程{,}n X n T ∈,若对于任意整数n T ∈和任意011,,,n i i i I +∈,条件概率满足11001111{|,,,}{|}n n n n n n n n P X i X i X i X i P X i X i ++++=======则称{,}n X n T ∈为离散时间的Markov 链,简称Markov 链(Markov chains )或马氏链.从定义可以看出:Markov 链具有Markov 性(即无后效性),如果把时刻n 看作现在,那么,1n +是将来的时刻,而0,1,2,,1n -是过去的时刻.Markov 性表示在确切知道系统现在状态的条件下,系统将来的状况与过去的状况无关,而且Markov 链的统计特征完全由条件概率11{|}n n n n P X i X i ++==所决定. 因此,如何确定这个条件概率,是研究Markov 链理论和应用中十分重要的问题之一. 4.1.2 转移概率定义 4.2 称条件概率1(){|}ij n n p n P X j X i +=== (4.1)为Markov 链{,}n X n T ∈在时刻n 的一步转移概率,其中,i j I ∈,简称转移概率(transition probability ).一般地,转移概率()ij p n 不仅仅与状态,i j 有关,而且与时刻n 有关,如果()ij p n 不依赖时刻n 时,则称Markov 链具有平稳转移概率.定义 4.3 若对任意,i j I ∈,Markov 链{,}n X n T ∈的转移概率()ij p n 与n 无关,则称Markov 链是齐次的(或称时齐的)(time homogeneous -),并记()ij p n 为ij p . 下面只讨论齐次Markov 链,并且通常将“齐次”两字省去.定义 4.4 设P 表示一步转移概率ij p 所组成的矩阵,且状态空间{1,2,}I =,则1112121222...........................n n p p p P p p p ⎛⎫ ⎪= ⎪ ⎪⎝⎭称为系统状态的一步转移概率矩阵(transition probability matrix ),它具有性质: (1)0,,ij p i j I ≥∈; (2)1,ijj Ipi I ∈=∈∑.(2)式说明一步转移概率矩阵中任一行元素之和为1,通常称满足性质(1)(2)的矩阵为随机矩阵.定义 4.5 称条件概率(){|},n ij m n m p P X j X i +=== ,,0,1i j I m n ∈≥≥ (4.2)为Markov 链{,}n X n T ∈的n 步转移概率,并称()()()n n ij P p =为Markov 链{,}n X n T ∈的n 步转移矩阵.其中()()0,1n n ij ij j Ip p ∈≥=∑,即()n P 也是一个随机矩阵.特别地,当1n =时,(1)ij ij p p =,此时,一步转移矩阵(1)P P =.我们还规定(0)0,1,iji jpi j ≠⎧=⎨=⎩Markov 链n 步转移概率满足重要的Chapman Kolmogorov -方程(简称C K -方程)。
第四章习题解答4.1Y1,Y2,···是来自总体Y的随机变量,与X0独立,h(x,y)是实函数.对于n 1,取X n=h(X n−1,Y n).设{X n}的状态空间为I,验证{X n}是马氏链,给出转移概率p ij.解:由题知,Y k与X1,···,X k−1独立,k 1,∀n,i,j,i1,...,i n−1∈I有,P(X n+1=j|X n=i,X n−1=i n−1, (X0)i0)=P(h(i,Y n+1)=j|X n=i,X n−1=i n−1,···,X0=i0)=P(h(i,Y n+1)=j|X n=i)=P(h(i,Y)=j)=P(h(i,Y1)=j|X0=i)=P(X1=j|X0=i).∴X n是马氏链,P ij=P(h(i,Y)=j).4.2设{X i,i 0}是取非负整数值的独立同分布的随机变量序列,V ar(X0)>0.验证以下随机序列是马氏链:(a){X n,n 0};(b){S n,n 0},其中S n=∑ni=0X i;(c){ξn,n 0},其中ξn=∑ni=0(1+X i).解:∀n,i,j,i0,···,i n−1∈N+,(a).P(X n+1=j|X n=i,X n−1=i n−1,···,X0=i0)=P(X n+1=j)= P(X n+1=j|X n=i)=P(X1=j)=P(X1=j|X0=i).1第四章离散时间马尔可夫链第四章离散时间马尔可夫链(b).P(S n+1=j|S n=i,S n−1=i n−1,···,X0=i0)=P(X n+1=j−i|X n=i−i n−1,···,X0=i0)=P(X n+1=j−i)=P(X n+1=j−i,S n=i|S n=i)=P(S n+1=j|S n=i)=P(X1=j−i)=P(X1=j−i|X0=i)=P(S1=j|S0=i).(c).P(ξn+1=j|ξn=i,ξn−1=i n−1,···,ξ0=i0)=P(X n+1=ji −1)=P(X n+1=ji−1|ξn=i)=P(ξn+1=j|ξn=i)=P(X1=ji −1)=P(X1=ji−1|X0=i)=P(ξ1=j|ξ0=i).4.3马氏链的状态空间是I=(1,2,3,4,5),转移概率矩阵P=0.20.80000.50.5000000.50.500.20.3000.500001界定马氏链的状态。