(完整word版)马尔可夫链应用于天气预报

1、马尔可夫链应用于天气预报摘要：在概率论与随机过程课中学习了马尔可夫链，马尔可夫过程因其无后效性、遍历性和时齐性，在科学研究、天气预测、农业预测、市场预测等方面应用非常广泛。本文通过对马尔可夫链理论和切普曼柯尔莫哥洛夫方程的探讨，结合天气因素、 降水情况的不确定性和无后效性等诸多特点，构建了基于天气预报的马尔可夫链预测模型，文中给出了马尔可夫链的一步转移概率矩阵和多重转移概率的计算方法，根据此算法可以预报短期天气情况，达到预测天气的目的。关键字：马尔可夫链天气预报转移概率切普曼柯尔莫哥洛夫方程1 引言天气变化情况与人们的生产、生活息息相关， 是人们普遍关注的重点问题之一。所以天气预报的准确性与时

2、效性就显得尤为重要，否则将对人们带来不便，甚至有可能带来重大经济和人员损失。 本文借助随机过程中著名的马尔可夫链模型，以某日天气的状态转移数据为例，建立了天气情况预测模型，并借助该模型应用马尔可夫链的遍历性，对未来天气的变化趋势作出了预测分析。由于马尔可夫过程应用广泛，它的重要特征是无后效性和遍历性。因此，运用马尔可夫链，只需要最近或现在的动态资料则可按转移概率可预测将来，这样就可以很方便地达到预测天气变化的目的。2 马尔可夫链预测模型2.1马尔可夫链的概念和特性马尔可夫过程是指具有以下特性的过程：过程X(t) （或系统）在时刻t 0 所处的状态为已知的条件下, 过程在时刻t t 0 所处状态

3、的条件分布与过程在时刻t 0 之前所处的状态无关，只与时刻t 0 所处的状态有关，这种特性称为马尔可夫性或无后效性。则称X(t) 为马尔可夫过程。马尔可夫链实际上就是状态和时间都是离散的马尔可夫过程。这一特性可用分布函数来确切地表出： 设随机过程 X(t),tT ，状态空间为，若对于 t的任意 n 个值 t 1t 20,i , j =1, ,N ，则此链具有遍历性， 且有极限分布=(1,2,N) ，它是方程组NP或ji pijj1,2,., Ni1的满足条件j 0,j1的唯一解。j在定理条件下，马氏链的极限分布是平稳分布，即若用作为链的初始分布，p(0)= ，则链在任一时刻n 的分布 p(n)

4、 永远与一致.因为p( n )p( 0)P( n)P nP n 1P2.2Chapman-Kolmogorov （切普曼柯尔莫哥洛夫）方程设 Xn， n=0,1, 为齐次马氏链，则对于任意的正整数k， m，有Pij ( mk )Pir ( m)Prj ( k ) ，此方程称为Chapman-kolmogorov( 切普曼柯尔莫哥r洛夫 ) 方程，简称C-K 方程。如果把转移概率写成矩阵的形式，那么CK 方程具有以下简单的形式nP(m+k)=P(m)P(k)，m, k 0，特别地， P(n)=P , n步转移概率由一步转移概率完全决定。这样，如果知道了马尔可夫链的初始概率， 即初始时刻各个状态的

5、概率， 并且知道它的一步转移概率矩阵，进而求得所有有限维概率分布，由此便可进行未来天气情况的预测，进行预报 。3 马尔可夫链预测天气如果明天是否有雨仅于今天的天气（是否有雨）有关，而与过去的天气无关，并设今天下雨的情况下，明天有雨的概率为；今天无雨的情况下，而明天有雨的概率为 ；又假定把有雨称为0状态天气， 把无雨称为 1状态天气， 则本例是一个两状态的马尔可夫链，其一步转移概率矩阵为Pp00p011p10p1113设 =0.7 ， =0.4 ，则一步转移概率矩阵为0. 70. 3P0. 60. 4于是两步转移概率矩阵为P 2P2PP0. 70. 30. 70. 30. 610. 390. 4

6、0. 60. 40. 60. 520. 48由此，可预报后天的天气情况，今天有雨，后天仍有雨的概率为p （ 2）=0.61 ，后天无雨概00率为 p01（ 2）=0.39 ；今天无雨，后天有雨的概率是p10（ 2）=0.52 ，后天无雨的概率是p11（ 2）=0.48 。同理四步转移概率矩阵为P 4(P( 2) 20. 57490. 42510. 56680. 4332据此可预报四日后的天气状况，今日有雨，第五日仍有雨的概率是p00 （4） =0.5749 ，今天无雨，则第五日的降雨概率p（4） =0.5668 。10六步转移概率矩阵为P 6P(4) * P(2)0. 57490. 42510

7、. 610. 390. 57170. 42830. 56680. 43320. 520. 480. 57100. 4290据此可预报六天后的天气状况，今日有雨，第七日仍有雨的概率是p00 （6） =0.5717 ，今天无雨，则第七日的降雨概率p（6） =0.5710 。10因此，根据今日的天气状况，可由马尔可夫链的多步转移概率矩阵求得短期内的天气情况的概率，从而实现对天气的预报。这是马尔科夫链在天气预报方面的应用。如果我们一直重复作下去，推得每一天的状况，一直到最后数字趋于稳定，我们称之“ 终极状态 ”。 这个状态由马尔可夫链的遍历性可以求得。设其极限分布为112 112(1 ,2)，则由P

8、得到，0. 720. 40. 57140. 30.6 ，解得120. 428612由此可预测将来有雨的概率比较大，为0.5714 。4 结论4马尔可夫链是一种很重要的离散随机过程， 简单说来就是用来求出一个事件的后续发展可能，当已知“现在”时刻过程的状态，那么过程“将来”的情况与“过去”的情况是无关的。或者说，这种随机过程的将来只通过“现在”与“过去”发生联系，如果一旦知道“现在”，那么“将来”就和“过去”无关了。通过研究马尔可夫链，可观察后面的发展状况，利用初始分布和一步转移概率就可以完全确定马尔可夫链的统计规律。天气预报是建立在对天气数据的分析统计之上的。因此数据越多， 越精确， 预测也就越可靠。利用马尔可夫模型对天气情况做出合理科学的预测分析， 在各种天气变化预测中有着实际的重要意义，为气象管理部门预报天气和决策提供科学的理论依据。当然，马尔可夫链的理论有趣，所用到的数学工具也较简单，广泛应用于