随机过程报告

马尔可夫链马尔可夫链是一种特殊的随机过程，最初由A.A.Markov所研究。它的直观背景如下:设有一随机运动的系统E(例如运动着的质点等)，它可能处的状态记为E,E 总共有可数个或者有穷个。这系统只可能在时刻t=l,2,・・・n,…上01n改变它的状态。随着工的运动进程，定义一列随机变量Xn,n=0,l,2,其中Xn=k，如在t二n时，工位于Ek。定义1.1设有随机过程&，neT}，假设对任意的整数neT和任意的ni,i,...ieI,条件概率满足0 1n+1P{X二iX二i，…,X二i}二P{X二i|X二i}n+1n+1 0 0 nn n+1n+1nn则称&，neT}为马尔可夫链，简称为马氏链。n实际中常常碰到具有以下性质的运动系统工。如果己知它在t二n时的状态，则关于它在n时以前所处的状态的补充知识，对预言工在n时以后所处的状态，不起任何作用。或者说，在己知的“现在”的条件下，“将来”与“过去”是无关的。这种性质，就是直观意义上的“马尔可夫性”，或者称为“无后效性”。假设马尔可夫过程&，neT}的参数集T是离散时间集合，即T={0,1,2,・・・},n其相应Xn可能取值的全体组成的状态空间是离散状态空间1={1,2,..}。定义1.2条件概率P(n)=p{X=j|X=i}ij n+1 n称为马尔可夫链&，neT}在时刻n的一步转移矩阵，其中i，jeI,简称为转n移概率。一般地，转移概率P(n)不仅与状态i,j有关，而且与时刻n有关。当P(n)不依赖ij ij于时刻n时，表示马尔可夫链具有平稳转移概率。假设对任意的i,jeI，马尔可夫链Xn,nGT｝的转移概率P（n）与n无关，则称马尔可夫链是齐次的。定义1.3设p表示一步转移概率p,所组成的矩阵，且状态空间l=｛l,2,・・・n｝，则称为马尔可夫链的一步转移概率矩阵。它具有如下性质：⑵丫P二1,ieI⑴P>0,i,jeI;. 丿刃jij定理1.1设&,neT｝为马尔可夫链，则对任意的ii，….ieI和n>1,有n 1,2 n巩船二爪…也7｝二工口附…口皿O这说明马尔可夫链的有限维分布完全由它的初始概率和一部转移概率所决定。因此，只要知道初始概率和一部转移概率，就可以知道马尔可夫链的统计特性。定义1.4假设｛Xn，n>0｝是齐次马尔可夫链，其状态空间为I,转移概率为Pij,称概率分布rj，je1 ｝为马尔可夫链的平稳分布，假设它满足◎二》兀舄工巧二込>0对于不可约马尔可夫链，假设它的状态是非周期，正常返的，则它是遍历的；对于不可约马尔可夫链，假设它的状态是有限且非周期的，则它是遍历的。值得注意的是，对于一个马尔可夫链，并不是一定存在limp（n）。例如设马尔可ns夫链的一部转移矩阵为：（01）P-U。丿易知p（2n）二I（单位矩阵），p（2n+1）二p，所以limp（n）不存在。在随机过程理论中，马尔可夫链是一类占有重要地位，具有普遍意义的随机

过程。它广泛应用于现代社会的各个领域，尤其在预测领域有着广泛的应用。马尔可夫链的预测方法分为很多种。根据指标值序列分组有3种。1〕数据序列约定俗成的分组方法：根据人们长久的经验进行分组：由于人们在现实生活中积累了生活经验，人们对认识的事物有了感性的了解，就可以对现象进行分组。2〕样本均值一均方差分组法：对于数据序列x,xx，可看作是一个时间序列的前n个观测值，算出样本均值X和1 2 n样本均方差S,根据具体情况以样本均值为中心，S为标准进行分组。3〕有序聚类分组法：有序聚类是对有序样品进行分类的一种方法，更加充分地考虑序列的数据结构，使划分的区间更加合理。有序聚类实现的经典算法是Fisher算法，其基本原理为:设时间序列x,xx的某一归类是1 2 n(兀］X什人J>i定义其均值向量为将公式D亿沪士区-和⑱—茅i=f定义为｛x,x x｝的直径，其含义表示该变量段内部各变量之间的差异情况。1 2 n其值越小，表示该段内变量之间差异越小，或说相互间越接近;反之，表示该段内变量之间差异越大，或说相互间越分散。三种马氏链预测方法：1)基于绝对分布的马尔可夫链预测步骤1对历史数据进行分组；步骤2确定观测值的状态，写出频数矩阵(n)，和一步转移矩阵(f),iji，jeE iji，jeE其中f= ，其中n为样本容量，当时nTa，可用频数估计概率p=f，从ijn-1 ijij而得到一步转移概率矩阵0=B。1lj步骤3“马氏性”检验步骤4已知时刻l时系统取各个状态的概率可视为马尔可夫链的初始分布，比方xl取状态2,m=5，则始分布P(0)=〔0,1,0,0,0〕，于是l+l时的绝对分布P=PP=(Pd),P(2),P(3),P(4),P(5))，可认为时刻1+1时系统所取的状态j满足(1) (0) 1 2 3 4 5Pj二max{P(i)}，从而预测1+t时刻的状态。1 1士51步骤5还可以用马氏链的平稳性，遍历性对系统分析。2)叠加马氏链预测步骤1对历史数据进行分组；步骤2计算各阶的一步转移矩阵P，P，…P，I={l,2,...k}，其中P=(f2)，1 2k 2 iji,jeEn(2)f2= ij，其他类推。ij n-2步骤3“马氏性”检验步骤4如果要预测时刻1+1的状态，可分别利用1,1T,,1-k+1作为初始态，，1+1所处的状态j满足P二max{P}。列表分析⑴1士5⑴初始甘段状怂1734 \来源I4%1耳1-23l-~i21-43件计%、%)A5)图1叠加马氏链预测分析表步骤5重复步骤1-4递推预测；步骤6进行平稳性，遍历性及其他分析。3〕加权马氏链预测步骤1对历史数据进行分组；步骤2计算各阶的一步转移矩阵P,P，…P,I={l,2,...k}，其中P=(f2)1 2k 2 iji,jeEf2亠,其他类推。an-2步骤3“马氏性”检验；步骤4计算各阶相关系数："上 _ —工（兀-―刘匕戟一兀）%=——--——-_— eE、£也-xfr=l计算标准的相关系数：渔 ―"刀口I粽£步骤5预测n+1时刻的状态步骤6重复1-5,预测n+2时刻的状态，其余类推步骤7讨论其他性质。马尔可夫预测方法是马尔可夫链在预测领域的一种应用方法。最初这种方法在水文，气象，地震等方面有广泛的应用，之后经济学家将其应用于研究市场占有率，预测经营利润等方面。在马尔可夫预测方法中，一个非常重要的问题就是对一步状态转移概率矩阵的估算。下面以实例分析马尔可夫链在现实生活中的应用。下面给出长江水域6类水质所占的比例。年符时间11类IU类IV覽'类劣V类014.659.5血92.71.72.5映112.216.254.49.53.83,91W7110.513.643.226.(J3.22g1%925.64fi.28.22.73.4199945,fl-)1.730.312.44,95.757.935.2Si'i14.95.96価］654.3Bfi.1156.37.420027o.s35.73L1J2.9134S1.3卫1,54L7g.&3.4145200490.926.740.215.J5.211.9现在要对长江未来10年的水质污染的发展趋势做一个总体的预测。为此可建立长江水质污染的马尔可夫链趋势预测的一步转移概率矩阵估计的最优化模型。设枯水期长江全流域水质在第t年属于I类、II类、III类、W类、V类、劣V类这6类状态的比例向量分别为a(t)=(P(l),P(2),...P(6)),t=0,1,2,....9•设P=(p)为6类状态矩阵的一步转t t t ij6x6移概率，根据误差平方和到达最小的准则，建立如下最优化模型：min/'(P)=工(£十I)—q(f)尸)(◎(［十I)-a(t)P)1t=06J,5"=1I2….6y=用matlab软件求解得r0.46070.00000.53930.00000.00000.0000-10,03450.27780.50490.06120.06030.0614P=0.00000.27050.33080.32300,00000.0142().()(X)0I.()0000.00000JW00().00()00.0000L().00()0(L11070.17470.00000.00000.1746」由下式q(9+k)二q(9)P\A=I,2,…JO可以对长江未来10年的水质污染属于I类、II类、III类、W类、V类、劣V类这6类状态的比例向量作出预测，预测结果见下表年恃时间11类II类in类IV类Y类劣*类20)15Lh3.76fig30.338336268】14.a.K25011.]S572U06]|5. 7655L58.005U]3»5713J.S27a]也SIUE2007J2553962R.235137.9222J4.O2933.R220JO.460421HJK135.784928.219L牺一刪3IS.^76*3.79(2]10-189I2期145歸928.I7&O38.1442I4.0J693.70aI9.W6020L01562S.18133K.19J514.OU63.80049.85762011165.9701?28.19433S.229814.M>6?.血72012175.W7228.206738.255114.07373.80669.69342013185.997328,217338.2729I4.QS273.8087鼻碌72014196.005828.225238-2855H.卿0S.SIO39.6144从预测计算结果可以看出：枯水期长江全流域水质属于W类、V类、劣V类这3类状态的比例并没有发生根本性的减少，水质污染程度依然十分严重。因此我们要采取积极措施，例如要严加控制企业废水和城市生活垃圾乱排乱放，政