马尔可夫链在天气预测中的应用.docxVIP

马尔可夫链在天气预测中的应用 龚海涛（数学系,093班25号） 摘要：马尔可夫链是一种预测方法，模式先假设某一时间各种状态之间的转移概率是基于 当询状态的而为其他因索无关，然后利用这一转移概率来推测未来状态的分布情况。木文将 利用马尔可夫链对鞍山市区天气状态进行探究，通过对鞍山市区从2010年2月7号到2012 年2月6号共730天的天气历史经验数据进行马尔可夫链分析,得到鞍山市天气状况的稳定 分布。 关键字：马尔可夫链；转移概率矩阵 马尔可夫链模型（Markov Chain Model）是一种常用的概率模型也叫马尔可夫分析 （Markov Chain Analysis）,其原理为利用概率转移矩阵所进行的模拟分析。此模型为一动态模 型，参数可随时间而变，故可以用来预测未来事物变化状态的趋势。 马尔可夫链的基本概念是在1907年由俄国数学家马尔可夫（Markov）从布朗运动 （Brown motion）的研究中提出的，后经由 Wiener、Kolmogorve Feller Doeblin 及 Lery 等人的研究整理而于1930到1940年代建立此模型（杨超然，1977）。 二、马尔可夫链的基本介绍 定义2.1 （Markov过程）随机过程｛Xn ,n=0,l,2,3,-）若它只取有限或可列个值E0,EbE2,-俄 们用｛0, 1, 2,???｝来标记Eo,Ei,E2,…,并称它们是过程的状态。｛（）, 1, 2, ???｝或其子 集记为S,称为过程的状态空间）对任意的2 0及状态… 诂有 P｛Xn+】=jlX（）=i）,X]=i|, -Xn.1=in.1,Xn=i｝=P｛ X^^jlXn^｝ （2.1） 式（2.1）刻画的Markov链的特性称为Markov性川。 Markov链表示一个随机序列的条件概率只与最近的系统状态有关，而与先询系统状态 无关，所以Markov性也被称为无后效性向。Markov性也可以用一句通俗的话来概括一一已 知现在，将来与过去无关。 定义2.2（转移概率）称式（2.1）中的条件概率P｛ 乂”冃风弓｝为Markov链｛X“ ,n=0,l,2,3,…｝ 的一步转移概率，简称转移概率山。 定义2.3 （时齐马尔可夫链）当Markov链的转移概率P｛ X.+冃IXR｝只与状态i,j有关，而 与n无关时，称Markov链为时齐的，并记二P｛ X.+冃IXR｝（nn 0）。 不管Markov链的状态是否有限，我们都可以将P? （i,jw S）排成一个矩阵的形式，令 ^PooP01P02 ^Poo P01 P02 P03 Pio Pl! P12 P13 p = (Pij)= P20 P21 P22 P23 P30 Psi P32 P33 P40 P41 P42 P43 (2.2) 定义2.4（转移概率矩阵）称式（2.2）为转移概率矩阵，容易看出Pij （i,jeS）有性质 （1） P?二 0, i,jeS （2） 三 Pij=l, Vi e S （2.3） i 定理 2.1 （Chapman-kolmogorov,C-K 方程） （…、 m ?或p的叫^（叫阿 （2.4） x U J r ik r kj k?S 其中P； =P{ Xn+m=jlXn=i}为m步转移概率，p m = 为m步转移概率矩阵。 三、数据分析 因为今天的天气状况很显然与昨天有一定关系而与前天及更早前的关系不大，即天气具 有无后效性，所以我们就可以用马尔可夫模型来对未来天气进行模拟预报。而且这种预测也 是很有意义的，因为有一句老话说“天有不测之风云”，所以如果我们能将未来的天气状况 预测出来，那对我们的生产生活都很有帮助。当然了天气预报更科学而且也更准确，我做此 文是想仅从历史经验数据出发来预测，因为两种预测方法迥然不同，所以我的预测与天气预 报没有可比性。 状态空间的类 天气有很多种状态，比如说晴、晴转多云、多云、小雨、中雨、大雨等等。为了简 化研究我按降水与否以及FI照或降水强度将天气状态简单分为以下四类，具体分类标准 见表1： 衣1:天气分类标准农 状态编号 天气简称 具体天气 0 晴 晴、晴转多云等 1 云 多云、阴、多云转晴、多云转阴、多云转小雨（雪）、雾等 2 小 小雨（雪）、小雨（雪）转晴、小雨（雪）转多云、小雨（雪）到中雨（雪）、阵雨等 3 大 中雨（雪）、中到大雨（雪）、大雨（雪）、大到暴雨（雪）、暴雨（雪）、雷阵雨等 原始历史天气数据来口 天气风雨录”网站（/anshan/tianqi）。 一共录得从2010年2月7日到2012年2月6日共计730日的历史天气状况，根据上表 的分类标准我们可以将原始数据转换成如表2所示的天气状况数据。 表2 : 2010年2月7日到2012年2月6日的天气状况衣 月份 天气状况（按时