算法大全马氏链模型

《算法大全马氏链模型》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、算法大全马氏链模型第十七章马氏链模型§1随机过程的概念一个随机试验的结果有多种可能性，在数学上用一个随机变量（或随机向量）来描述。在许多情况下，人们不仅需要对随机现象进行一次观测，而且要进行多次，甚至接连不断地观测它的变化过程。这就要研究无限多个，即一族随机变量。随机过程理论就是研究随机现象变化过程的概率规律性的。定义1设{ξt,t∈T}是一族随机变量，T是一个实数集合，若对任意实数t∈T,ξt是一个随机变量，则称{ξt,t∈T}为随机过程。T称为参数集合，参数t可以看作时间。ξ的每一个可能取值称为随机过

2、程的一个状态。其全体可能取值所构成的集合称为状态空(t)间，记作E。当参数集合T为非负整数集时，随机过程又称随机序列。本章要介绍的马尔可夫链就是一类特殊的随机序列。例1在一条自动生产线上检验产品质量，每次取一个，“废品”记为1，“合格品”记为0。以ξn表示第n次检验结果，则ξn是一个随机变量。不断检验，得到一列随机变量ξ1,ξ2,L，记为{ξn,n=1,2,L}。它是一个随机序列，其状态空间E={0,1}。例2在m个商店联营出租照相机的业务中（顾客从其中一个商店租出，可以到m个商店中的任意一个归还），规定

3、一天为一个时间单位，“ξt=j”表示“第t天开始营业时照相机在第j个商店”，j=1,2,L,m。则{ξn,n=1,2,L}是一个随机序列，其状态空间E={1,2,L,m}。例3统计某种商品在t时刻的库存量，对于不同的t，得到一族随机变量，{ξt,t∈[0,+∞)}是一个随机过程，状态空间E=[0,R]，其中R为最大库存量。我们用一族分布函数来描述随机过程的统计规律。一般地，一个随机过程{ξt,t∈T}，对于任意正整数n及T中任意n个元素t1,L,tn相应的随机变量ξt,L,ξt1n的联合分布函数记为Ft1

4、Ltn(x1,L,xn)=P{ξt1≤x1,L,ξtn≤xn}（1）由于n及ti(i=1,L,n)的任意性，（1）式给出了一族分布函数。记为{F(x,L,x),t∈T,i=1,L,n;n=1,2,L}tLt1ni1n称它为随机过程{ξt,t∈T}的有穷维分布函数族。它完整地描述了这一随机过程的统计规律性。§2马尔可夫链2.1马尔可夫链的定义现实世界中有很多这样的现象：某一系统在已知现在情况的条件下，系统未来时刻的情况只与现在有关，而与过去的历史无直接关系。比如，研究一个商店的累计销售额，如果现在时刻的累计

5、销售额已知，则未来某一时刻的累计销售额与现在时刻以前的任一时刻累计销售额无关。上节中的几个例子也均属此类。描述这类随机现象的数学模型称为马氏模型。定义2设{ξn,n=1,2,L}是一个随机序列，状态空间E为有限或可列集，对于任意的正整数m,n，若i,j,ik∈E(k=1,L,n.1)，有-207-P{ξξn+m=j

6、ξn=i,ξn.1=in.1,L,ξ1=i1}=P{ξn+m=j

7、ξn=i}（2）则称{ξn,n=1,2,L}为一个马尔可夫链（简称马氏链），（2）式称为马氏性。事实上，可以证明若等式（2）对

8、于m=1成立，则它对于任意的正整数m也成立。因此，只要当m=1时（2）式成立，就可以称随机序列{ξn,n=1,2,L}具有马氏性，即{ξn,n=1,2,L}是一个马尔可夫链。定义3设{ξn,n=1,2,L}是一个马氏链。如果等式（2）右边的条件概率与n无关，即P{ξn+m=j

9、ξn=i}=pij(m)（3）则称{ξn,n=1,2,L}为时齐的马氏链。称pij(m)为系统由状态i经过m个时间间隔（或m步）转移到状态j的转移概率。（3）称为时齐性。它的含义是：系统由状态i到状态j的转移概率只依赖于时间间隔的长

10、短，与起始的时刻无关。本章介绍的马氏链假定都是时齐的，因此省略“时齐”二字。2.2转移概率矩阵及柯尔莫哥洛夫定理对于一个马尔可夫链{ξn,n=1,2,L}，称以m步转移概率pij(m)为元素的矩阵P(m)=(pij(m))为马尔可夫链的m步转移矩阵。当m=1时，记P(1)=P称为马尔可夫链的一步转移矩阵，或简称转移矩阵。它们具有下列三个基本性质：（i）对一切i,j∈E，0≤pij(m)≤1；（ii）对一切i∈E，∑pij(m)=1；j∈E（iii）对一切i,j∈E，pij(0)=δij=...10,，当当

11、ii=≠jj时时。当实际问题可以用马尔可夫链来描述时，首先要确定它的状态空间及参数集合，然后确定它的一步转移概率。关于这一概率的确定，可以由问题的内在规律得到，也可以由过去经验给出，还可以根据观测数据来估计。例4某计算机机房的一台计算机经常出故障，研究者每隔15分钟观察一次计算机的运行状态，收集了24小时的数据（共作97次观察）。用1表示正常状态，用0表示不正常状态，所得的数据序列如下：11100100111111100111