根据马氏链的油价预测模型

1、数学建模基于马氏链的油价预测模型油价预测模型摘要近年来中国及国际油价变化多端、震荡无常，这种“杂乱无章”的变化为油价的走势推断带来了困难。 而石油价格又是政府部门、原油生产部门和用油单位及投资者关注的焦点， 因此，准确并正确的预 测未来较长时期的油价变动趋势，将具有重要的意义。本文将使用马 尔科夫链模型来预测油价未来的变化情况。问题重述1.基于0512年的原油价格，预测下半年及13年的油价情况，并与国际油价相对比。2.对国家石油机制提出自己的建议。模型的假设1 .政府干预没有对成品油油价造成影响2 .国内进口原油占国内总原油产量的比例不变3 .国内开采原油的成本与进口原油油价相同4 .收集的数

2、据真实有效模型的建立数据的收集表一 0512年原油价格（单位：美元/桶）20052006200720082009201020112012年年年年年年年年1月40.2658.1756.9193.8143.16576.1291.29100.992月41.4856.6359.4095.4142.7577.0693.3296.915553月49.0757.8663.95101.944.8581.08102.9107.0555164月49.4664.3865.90107.349.4085.3110.5103.25565月46.6265.1764.2212059.2480.09107.8101.35555

3、6月51.5164.5967.7814069.2074.1497.885.18557月5302768.9774.65132.669.2975.995.245558月57.8388.8175.23120.370.7876.5991.89559月57.9959.3477.81108.169.2876.9683.935510月54.3954.9787.2383.1673.8581.5185.26555511月51.155.4291.0659.5777.6383.5396.2955512月52.4857.9592.6246.9778.8989.0699.9255油价转移状态划分由表可知，油价起伏变化还

4、是很大的，从2005年1月至2008年6月，原油价格从40.26美元/桶提高到140 美元/桶。此间经历了若干次跳跃式增长，每次跳跃都使 油价跃升至一个新的台阶，因而可将这些价格“台阶”看 作是油价的转移状态。由分析可得出，可将油价分为五个转 移 状 态 区 间 （0,50）? 50,70) ? 70,90) ? 90,110) ? 110, + ,并分别因而可 将70 , 90）区间看作是中高态，50, 70）区间看作是中低态，这二者可并称中态；（0 , 50）则为低态；90, 110）则为高态，由于110美元后出现超常飙升，可将110 , + ） 看作是超咼态。油价转移状态的构成如果我们用

5、E来表示原油状态（事件），用E表示油价低状态区间（0,50），用Eml表示50,70），用Emh表示70 ,90），用 Eh 表示90 , 110），用 Eeh 表示110 , + 呵，则2005年1月到2012年6月的油价转移过程可归结为:E|7 E|7 E|7 E|7 E|7 Eml Eml Eml Eml EmlEml Eml Eml Eml Eml Eml Eml Eml Eml Eml T Eml t Eml tEml t Eml t Eml t Eml tEml tEml tEml T Eml t Emh t Emh t Emh t Emh t Eht Eht EhtEhT EhT

6、 EhT Eeh t Eeh t Eeh tEeh t EhtEmh t EmlElT ElT ElT ElT ElT Eml t Eml t Eml t Emh t EmlEmh TEmh tEmh t Emh t Emh tEmh t Emh t EmhEmh T Emh t Emh t Emh t Emh t Emh t Emh t EhEhT EhT Eeh t Eh t Eht Eht Eht Emh t Emh t Eh tEhT Eh T EhT Eh t EhT Eh t Emh（此四.判断状态转移链是马尔科夫链ElEmlEmhEhEehni.El320005Eml216300

7、21Emh03182023Eh00326130Eeh0002810n.j5212430989表二 2005年1月至2012年6月油价各状态之间的转移频数马尔柯夫链的基础概念一一转移概率是条件概率，是以事件处为状态）的非独立性为前提的，因而在讨论原油价格跳动转移过程 是否是马尔柯夫链时，首先要对事件的非独立性进行检验， 这是通过 以下统计量来完成的.?=(n ij -n i. n jn) 2"£=1£=1 (n i. n .j/n)统计学上已证明，它服从自由度为(m 一 1)2的？?分布。nij ni. n.j式中m为状态数，n为状态转移总次数，ni.为第i个状态（

8、i=l , 2 ,? , m）转向其他各状态的转移次数之和，其数学表示为n i二zmi nij ,n.j为第j状态（j=l , 2, ? m）由其他各状态转来的次数之和，其数学表示为 n.j二 £5=1 nj。这里原假设Ho为状态E与E相互独立（i,j=1，2，?，5）。此原 假设下，频数矩阵应该是均匀的，其各元素值应相差不大，据这个公 式计算出的2值应该小于以自由度（m-1） 2和选定的显著性水平a （本文选定a =0 . 05）所限定的2表查值。这时在肯定原假设的同时，也肯定了状态的相互独立性，否定了马尔柯夫链的存在；如果相反，这 个计算值大于表查值，那么就要否定原假设，即是在所

9、选定的显著性 水平（a上肯定了状态转移的非独立性，也就是肯定了所考察的过程 是一个马尔柯夫链。表二原油价格状态转移链X2检验表(nij - ni.n.j/n) 2j=1( El )j=2(E ml )j=3(E ml)j=4(E h)j=5(E eh)(n i. n.j/n)i=1( E i )26.3200.570251.348311.685390.5056191538i=2 (E ml )0.570224.61941.252207.078652.12359516726i=3(Emh )0.387720.59881.653337.752802.32584871493i=4(Eh)0.0587

10、11.24373.2023810.11233.03370273868i=5(E:eh)3.681778.85470.034123.370781.01123989976m刀30.00279.42501.382445.056181.51685j=17444m 刀 i=1m 刀 j=1(nij -(nni.n.j/n) 2 i.n.j/n)117.3832这里自由度(m-1) 2=(5-1) 2=16，如取显著性水平a =0 . 05，则在 X2 (16 , a=.05 )统计表中查得 X2 (16 , a=.05 ) =26.3，而由样本计算得出的X2值为117.3832，远远大于X2(16 ,

11、a=.05 )=26.3 ,明显推翻了状态相互独立的原假设。这就是肯定了2005年1月至2012年6月原油价格的状态转移链是马尔柯夫链，因而可用马尔柯夫链的性质对其进行分析和讨论。五.马尔科夫链的极限频率马尔柯夫链的遍历性定理为震荡油价走势推断提供了理论依据。尽管马尔柯夫链的无后效性与回归思路有不 相容之处，但不等于说油价的走势是不可推断的，恰恰是 对于马尔柯夫链具有根本意义的遍历性定理为震荡期油价走势推断提供了理论依据和实现方法，从而打开了一个 新局面。遍历性的含义是：对于马尔柯夫链来说，从任意 状态E出发，当转移步骤次数逐渐增大时，转移概率都将趋于收敛，即达到它的极限值一一极限概率P,这就

12、是马尔柯夫链的各态遍历性定理。这时所得到的状态转移矩阵各行矢量相等，标志着过程已经稳定。如用Z(1)表示阶转移矩阵，用Z( a表示，a阶转移矩阵，则它们有如下关系。此乘幕不断进行下去就会趋向它的极限值，这为实 际应用提供了一个逼近极限值的方法。用表2中所在行的每个元素nij除以每行的行和ni.得到一阶状态转移概率矩阵:Z(1)=0.6000.4000.0000.0000.0000.0950.7610.1420.0000.0000.0000.1300.7820.0870.000 丿0.0000.0000.1000.8660.0330.0000.0000.0000.2000.8000.06910.

13、29020.31800.27670.04610.06910.29020.31800.27670.04610.06910.29020.31800.27670.04610.06910.29020.31800.27670.04610.06910.29020.31800.27670.0461此时状态转移矩阵已逼近它的收敛值一一极限概率,其二阶转移矩阵、？ ？、及收敛转移矩阵分别为:Z2=厂 0.39810.54480.05700.00000.00000.12970.63720.22070.01240.00000.01240.20140.63980.14350.00290.00000.01300.164

14、90.76650.0555Z128 =< o.。0000.00000.02000.33330.6467行概率变得完全一样，表明不论从哪个状态开始其最终结 果都是一样的。在经历了一个持续收敛的过程后，达到了 稳定，此时P,表明的是各状态Ei在整个过程中所占的分 量和比例。耐人寻味的是，马尔柯夫链从状态的非相互独 立性开始，以转移概率再不依附于它的前提条件，以状态 的相互独立为终结，在自我否定中实现了一次质的飞跃。正是这个质的飞跃使油价走势推断成为可能。这个稳定的 概率行矢量就是该马尔柯夫链的走势，在本例中也就是国 际原油价格的走势。至此我们完成了震荡油价走势推断模 型的建立，这个推断模型就

15、是状态转移收敛矩阵的行矢量, 即:马链极限概率值EEmlEmhEhEeh0.06910.29020.31800.27670.0461油价走势的解读与分析上述极限概率表明,国际油价的最终状态如下：低态油价终将占到6.91 %,中低态油价终将占到29.02 %，中高态油价终将占到31.8 %，高态油价终将占到27.67 %, 超高态油价终将占到4.61 %，其柱状图如图一所示。由图 一可见其形态为不对称的偏态分布，90美元/桶以下的油 价占到70 %,中位油价和低位油价占了大多数。构成了国 际油价的主体；而高于90美元/桶的油价占30%，处于 从属地位。图一 油价极限概率柱状图图表标题系列1少 少

16、为了对国际油价马尔柯夫链的极限概率进行解读，我们先对频数转移矩阵做如下的处理，对其中的每一个元素 nij除以状态转移总数n ,得到fij矩阵:ElEmlEmh&EehE0.03370.02250.00000.00000.0000Emi0.02250.17980.03370.00000.0000Emh0.00000.03370.20220.02250.0000Eh0.00000.00000.03370.29210.0112Eeh0.00000.00000.00000.02250.0899状态累0.05610.235950.269660.337070.10112计频率85394将状态累计频

17、率与马尔柯夫链收敛矩阵的行矢量作一个对比，即可发现它们之间的差异，见表4。表四极限概率与状态累计频率对比表日(0,50 )Eml 50.70 )Emh 70,90 )Eh90,110 )Eeh110, )马尔科夫0.0690.2900.3180.2760.046链极限频12071率值油价状态累计频率0.05620.23600.26970.33710.1011值二者之差值0.01290.05420.0483-0.0604-0.0550二者之差 值与累计 频率值之 比23.00%22.99%17.92%-17.91%-54.41%由表4可看到高于110美元/桶的超高油价状态出现的概率从0.1011

18、降低到0.0461，下降54.41 %; 90至110美元的高位油价状态出现的概率从0.3371降低到0.2767 ,下降17.91 %;而低位状0.0129、态、中位状态和中高位状态都有所上升，上升幅度分别为0.0542、0.0483，上升百分比分别为 23.00%、22.99%、17.92%。这说明油价的低位状态、中位状态和中高位状态出现的比重加大了, 这三者从56.18 %增加到67.73 %，增加了 11.55 %。油价的高位状态和超位状态出现的比重降低了，此两者从在总体中占43.82 %下降 到占32.27 %，减少了 11.55 %。这表明今后一段时间内国际原油价 格将有所下降。另

19、一方面，如果把70, 90）和90, 110）这两个状态区间合并起来看，会发现这个区间内的状态累计频率值和极限概率值 无明显变化，分别为59.47 %和60.68 %，也就是说70至110美元 区间内的油价在原油价格走高的总趋势下是相对稳定的。 在这个稳定 区间外，高于110美元/桶的油价呈下降趋势，低于70美元/桶的 油价呈上升趋势。如再将低态油价与中态油价归并于一起进行考虑, 低于90美元/桶的油价接近70 %的比重，对今后国际油价的走高趋 势起着牵制作用。如将咼位态和超咼位态也并为一组，统称咼态，可 得低、中、高三态对比柱状图，如图二所示。图二结论一目了然，今后油价的走势是，低位态油价所

20、占的百分比在总体有所上升，高位态油价所占的百分比在总体上有所下降， 中位态 油价在总体上所占百分比基本稳定。如果将EI Eml Emh Eh并在一起，其在出现的概率为 0.9539，今后一段时间，绝大多数的油价数据会落入这个区间。国0.05。即际油价高于110美元/桶的超常飙升事件出现的概率不足 所占份额不足油价总体的5%，并不能改变低、中态油价约占 70 % 的总趋势。模型中的不足和改进在上面的马尔科夫链模型中，预测未来油价的结果还是有很多问 题的，如赋予05年和12年油价以相同的权重是不合理的，明显是12年的权重应该更大些。另外，对于未来的预测是一个区间的概率,在现实生活中应用还是很困难的，如果得到的是一个具体数值，显然 会更好。而且，我发现，其实“小波分析”才是最好的预测模型，可以很 好地模拟油价走势曲线，给出相当准确的预测。但由于作者水平有限, 看了很多关于