基于灰色马尔科夫链模型的交通事故伤亡人数预测

1、基于灰色马尔科夫链模型的交通事故伤亡人数预测摘要：道路交通系统是一个基于人、车、路的动态系统，影响交通安全的因素很多，作用机理复杂，因此道路交通事故的发生具有很大的随机性和偶然性。传统的GM(1,1)模型和马尔科夫模型都能单独解决有关时间序列的预测问题，但各有优缺点：GM(1,1)模型能预测出事物发展的总体趋势和大体方向，对预期远、波动大的数据的预测误差较大；而马尔科夫模型对于波动性大的数据序列的预测精度较高，但其主要是对具有平稳随机过程的问题进行的预测，对现实问题中占绝大多数的非平稳过程问题的预测存在局限性。本文以灰色GM(1,1)模型为基础，利用马尔科夫链模型对灰色GM(1,1)模型的预测

2、结果进行误差修正，并利用某市交通事故伤亡人数的数据对之后几年的伤亡人数进行预测。通过对比，证明基于灰色马尔科夫链模型的交通事故伤亡人数的预测更加准确。关键词：交通事故预测；马尔科夫链；灰色GM(1,1)模型；误差修正1、引言交通安全是国民经济发展和社会安定的重要方面，也是道路交通管理的两项基本任务之一。道路交通事故预测是道路交通安全研究的一项重要内容，它的目的是为了掌握交通事故的未来状况，以便及时采取相应的对策，有效地控制各影响因素，避免工作中的盲目性和被动性，减少交通事故的发生。因此，准确地对交通事故进行预测具有重要的现实意义。道路交通系统的非线性、随机性、动态性以及不确定性等特点，决定了作

3、为道路交通系统行为特征量的道路交通事故预测的复杂性。本文根据现实生活中交通系统非线性、随机性和动态性的特点，将灰色GM(1,1)模型和马尔科夫模型的结合起来，使其优势互补，提高对交通事故预测的准确性。2、GM(1,1)模型客观世界的很多实际问题，其内部的结构、参数以及特征并未全部被人们了解，人们不可能象研究白箱问题那样将其内部机理研究清楚，只能依据某种思维逻辑与推断来构造模型。对这类部分信息已知而部分信息未知的系统，我们称之为灰色系统。灰色系统的研究对象是“部分信息已知、部分信息未知”的“小样本”、“贫信息”不确定性系统，它通过对“部分”已知信息的生成、开发实现对现实世界的确切描述和认识。信息

4、不完全是“灰”的基本含义。灰色系统理论建模的主要任务是根据具体灰色系统的行为特征数据，充分开发并利用不多的数据中的显信息和隐信息，寻找因素间或因素本身的数学关系。通常的办法是采用离散模型，建立一个按时间作逐段分析的模型。但是，离散模型只能对客观系统的发展做短期分析，适应不了从现在起做较长远的分析、规划、决策的要求。事实上，微分方程的系统描述了我们所希望辨识的系统内部的物理或化学过程的本质。由于灰系统对一切随机量都可看作是在一定范围内变化的灰色量，因此，为适应灰系统建模需要，提出“生成”的概念，“生成”即指对原始数据做累加（或累减）处理。累加生成一般可写成AGO。若计为原始数列，为次累加生成后数

5、列，即(2-1)则次累加生成算式为(2-2)一般常用的是一次累加生成，即(2-3)建立GM模型，实际就是将原始数列经过累加生成后，建立具有微分、差分近似指数规律兼容的方程，称为灰色建模。GM(m,n)称为m阶n个变量的灰色模型，其中GM(1,1)模型是GM(1,n)模型的特例，是灰色系统最基本的模型，也是常用的预测模型，其简单的微分方程形式（白化形式的微分方程）是(2-4)利用常数变易法解得，通解为(2-5)若初始条件为，则可得到微分方程的特解为(2-6)或时间响应函数(2-7)其中白化微分方程中的项中的为的背景值，也称为初始值;为常数（有时也将写成）。按白化导数定义有差分形式的微分方程，即(

6、2-8)显然，当时间密化值定义为1，即当时，上式可记为(2-9)记为离散形式(2-10)这显然表明是一次累计生成，因此上述方程可改写为(2-11)这实际也表明，模型是以生成数（是以的一次累加）为基础的。 当足够小时，到不会发生突变，因此可取与的平均值作为时的背景值，因此，背景值便可记为(2-12)或(2-13)于是白化的微分方程可改写为(2-14)或(2-15)即(2-16)因此，上述方程可以改写为矩阵方程形式，即(2-17)引入下列符号，设(2-18)于是便有(2-19)令(2-20)则(2-21)解得(2-22)将求解得到的代入微分方程的解式（也称时间响应函数），则(2-23)由于，因此求

7、导还原得(2-24)上述两式便为GM(1,1)的时间响应式，及灰色系统预测模型的基本算式，当然上述两式计算结果只是近似计算值。3、马尔科夫模型马尔科夫过程是由俄国著名数学家马尔科夫提出的，马尔科夫过程既适用于区间序列同样也适用于时间序列，是一个典型的随机过程，该理论主要是研究序列的状态和转移规律。假设一个序列有几种状态，该序列目前处于某种状态，而下一时间段可能会转移到另一个状态，通过研究各状态的初始概率及各状态之间的转移概率，来确定各状态的变化趋势，从而进行预测，这样离散时间之下的随机过程就是马尔科夫过程。马尔科夫链是最简单的马尔科夫过程。马尔科夫过程的数值是连续的，任意两值之间都可以无限分割

8、，状态也有无限多个。而马尔科夫链模型的时间以及状态参数都是离散数值，状态也是有限可列的。马尔科夫过程的特点是将来的状态只与现在有关，而与历史无关。也就是说，系统在时刻的状态，仅仅与时刻所处的状态有关，而不受时刻之前所处状态的影响。这种特点就是马尔科夫过程的无后效性也称为马氏性，当然马尔科夫链模型是特殊的马尔科夫过程同样也具有无后效性。3.1马尔科夫过程的无后效性假设马尔科夫过程为离散的时间集合，即，其状态空间为可数的。考虑有限维分布函数，对，及状态。由乘法公式可以得到：(3-1)上式最简单的情况为：(3-2)此时：(3-3) 显然上式中任意的都成立，那么相互独立。但在实际中这种都相互独立的情况

9、很少见，下式为不独立的情况：(3-4)此时：(3-5) 将上式中的性质称之为无后效性，即表示在当前的情况下，系统未来的变化不受过去的影响，只依赖于目前系统所处的状态。 设 是定义在概率空间上的随机过程，状态空间为，若对任意的，以及状态均有：(3-6)就称为马尔科夫链，此过程为马尔科夫过程。马尔科夫过程具有无后效性。上式的等价简化形式为：(3-7)3.2状态转移概率及其转移概率矩阵 设随机过程为一个马尔科夫链，状态空间为，每一个时间点，只能处于其中一种状态，因此每个状态都有个转向（包括转向自身），即：(3-8)则一步转移概率(3-9)称为时刻从状态经过一步转移到状态的概率。系统所有状态一步转移概

10、率集合所组成的矩阵称为一步状态转移概率矩阵。其形式如下：(3-10)此矩阵具有以下两个性质： 1）非负性： 2）行元素和为1，即那么，两步转移概率：(3-11)称为时刻从状态经过两步转移到状态概率。那么随机过程的两步转移概率矩阵为：(3-12)同理，步转移概率：(3-13)称为时刻从状态经过个时刻到状态的概率。因此，系统的步转移概率矩阵就是由所有状态的步转移概率集合所组成的矩阵。其形式如下：(3-14)那么，式(3-14)则为的步转移概率矩阵。此矩阵同样具有以下两个性质： 1）非负性： 2）行元素和为1，即 对于式(3-11)，一般情况下，不仅和状态有关，而且和时刻有关。当和时刻无关时，表示马

11、尔科夫链具有平稳转移概率。如果，马尔科夫链的转移概率与无关，就称马尔科夫链是齐次的。应用上主要研究的是齐次马尔科夫链。马尔科夫链的一步转移概率矩阵和步转移概率矩阵之间的关系：设是马尔科夫链的步转移概率，则，且有： (3-15)当马尔科夫链具有齐次性时，则转移概率就具备了平稳性，也就是仅与时间、间距、 和有关，即。此时，C-K方程可简写为，可得步转移概率矩阵为一步转移概率矩阵的次方。即： (3-16)3.3马尔科夫链预测模型 一个时间序列，其可能的观测数据可以取个离散的值，即序列可以处于个状态 。设序列处于状态的状态概率为。序列从转到一下状态的转移概率为，则由马尔科夫链的无后效性可知。引入状态概

12、率向量和一步转移概率矩阵式： (3-17) (3-18) 由状态转移的马氏性以及式(3-16)可以写出马尔科夫链的基本方程：(3-19)4、灰色马尔科夫模型4.1灰色GM(1,1)模型灰色GM(1,1)模型首先通过对原始数据进行累加，建立均值生成序列和矩阵B与Y，然后通过最小二乘回归和微分等数学方法建立模型，最后通过模型得到的值经过还原数据，得到预测结果。它的建模过程为：1）根据模型在各个时刻的值，建立如式(4-1)所示的原始数据序列(4-1)2）对原始数据序列进行累加，得(4-2)3）对(4-2)式序列作均值，生成序列(4-3)4）利用式(4-1)与式(4-3)，建立矩阵Y与B，得(4-4)

13、5）对参数进行最小二乘估计，得出a与b的值(4-5)6）确定模型形式，并还原得到的灰色预测值，如式(4-6)、式(4-7)所示(4-6)(4-7)4.2马尔科夫模型马尔科夫链是根据所观察的离散状态，以经验为主的估计转移概率参数 化的随机过程。它是对原始数据进行状态划分，求出转移概率矩阵，得出未来的预测值。以灰色马尔科夫链模型为例，其一般步骤如下：4.2.1状态划分根据灰色模型预测值与实际值间的相对误差，把相对误差分成r类状态。状态划分数量并无严格规定，是综合考量样本数量、拟合的误差范围等相关因素而确定，一般分成35类比较合适。4.2.2建立状态转移概率矩阵 假设是状态到的步转移概率，是状态到的

14、步转移次数，属于个状态的数量，状态转移概率矩阵如式(4-8)所示(4-8)式中：。4.2.3计算预测值 假设时间序列在时刻处于状态，根据状态的残差区间的中值，与灰色预测值，可以得出灰色马尔科夫链模型的预测值为，如式(4-9)所示:(4-9)4.3对模型精度的检验灰色预测模型建立以后，对模型的实用性以及模型的精度进行验证。GM(1,1)模型通过计算残差、平均相对误差、均方差比值、小误差概率等指标后，查找灰色预测模型精度检验等级表（见表1），从而可以判断模型的精度等级。表1 灰色模型的精度检测表 精度等级指标范围相对误差（）均方差比值（）小误差概率（ ）一级（好） 二级（合格） 三级（勉强合格）

15、四级（不合格）计算过程和算式如下：1）分别计算出原始数据序列的残差，相对误差与平均相对误差(4-10)(4-11)(4-12)2）分别算出原始数据与残差的标准差，。根据，分别算出均方差比值和小误差概率(4-13)(4-14)(4-15)(4-16)5、案例分析利用马尔科夫链对灰色GM(1,1)模型的预测误差进行修正，以某市20072013年的伤亡人数为基础，对某市20142016年的交通事故伤亡人数进行预测。5.1建立GM(1,1)预测模型灰色GM(1,1)模型的建立过程如下：1）原始数据序列为：；2）对数据进行累加得：；3）建立均值生成序列，矩阵与为(5-1)4）对参数进行最小二乘估计，得出

16、与的值(5-2)5）将和的值带入式(4-6)，得出模型如式(5-3)所示(5-3)根据式(5-3)，并根据式(4-6)还原数据，得出某市20072013年的伤亡人数灰色预测值，结果如表2所示。预测结果显示， 08年和09年的模型相对误差较大，分别为7.96%和-9.29%。最后可得到20142016年伤亡人数灰色预测值分别为813人、788人、763人。表2 某市交通事故实际伤亡人数与灰色模型预测值的对比年份实际伤亡人数灰色模型预测值残差误差/%20071047104700.0020081068983857.962009872953-81-9.292010902923-21-2.3320118

17、76894-18-2.052012846866-20-2.362013895839566.266）对GM(1,1)模型进行精度检验。利用式(4-10)式(4-16)可以算出，平均相对误差为4.32%，后验差比值为61.34%，小误差概率为0.7143。查找表1，可知该模型的精度为3级，说明可以用于交通事故预测，但精度较低，需要进一步优化来提高模型的精度。5.2建立马尔科夫链模型5.2.1状态划分 因为本研究样本数量较少，按照均值划分，误差可分为三个状态，分别用表示，如表3所示。表3 死亡人数状态划分表状态误差范围（-9.29%-3.54%）（-3.54%2.21%）（2.21%7.96%） 根

18、据表3中的状态划分情况，可以把20072013年交通事故伤亡人数进行状态划分，结果如表4所示。表4 某市20072013年交通事故实际伤亡人数状态划分情况年份实际伤亡人数灰色模型预测值误差/%状态2007104710470.00200810689837.962009872953-9.292010902923-2.332011876894-2.052012846866-2.3620138958396.265.2.2构建转移概率矩阵(5-4)(5-5)(5-6)5.2.3计算预测值利用式(4-9)对20072013年某市伤亡人数进行拟合。例如2008年的灰色预测值为983，处于状态，可以得出200

19、8年的灰色马尔科夫链预测值为1033人。同理，可以得出其余年份的预测值，两种模型的残差和误差情况如表5所示。表5a 灰色GM(1,1)模型预测结果年份实际伤亡人数灰色GM(1,1)模型预测值残差误差/%2007104710470020081068983857.962009872953-81-9.292010902923-21-2.332011876894-18-2.052012846866-20-2.362013895839566.26表5b 灰色马尔科夫链GM(1,1)模型预测结果年份实际伤亡人数灰色马尔科夫链GM(1,1)模型预测值残差误差/%20071047104700200810681

20、033353.282009872892-20-2.292010902917-15-1.662011876888-12-1.372012846860-14-1.652013895882131.45由上表可知，2008年和2009年的灰色GM(1,1)预测值相对误差为7.96%和-9.29%，而灰色马尔科夫链GM(1,1)预测值和相对误差降到3.28%和-2.29%。从图1可知，灰色GM(1,1)模型的预测值呈一条平滑递减曲线，而灰色马尔科夫链GM(1,1)模型的预测值具有一定的波动性，接近伤亡人数的实际值，预测结果更加可靠。图1 两种模型结果对比根据表4可知，2013年伤亡人数预测值处于状态，初

21、始行向量为。因此，说明2014年处于状态，再利用式(4-9)预测出2014的伤亡人数为761人。同理，可以预测2015年、2016年伤亡人数年所处的状态及预测值，结果如表6所示。表6 GM(1,1)模型与灰色马尔科夫链预测模型对20142016年伤亡人数预测值年份灰色GM(1,1)模型预测值状态状态区间中值灰色马尔科夫链GM(1,1)模型2014813状态1-6.37%7612015788状态2-0.665%7832016763状态3-0.665%7585.2.4对灰色马尔科夫链GM(1,1)模型进行精度检验利用马尔科夫链对灰色模型进行误差修正后，平均相对误差（越小越好）从4.32%降到1.67%，降低了2.65%；后验差比（越小越好）从61.34%降到22.04%，降低了39.3%；小误差概率（越大越好）从0.7143提高到1，提高了0.2857。查找表1，可知该模型的精度变为一级（好），说明灰色