国内成品油调价与定价机制模型的建立与求解

1、数学建模竞赛承诺书我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮 件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问 题。我们知道，抄袭别人的成果是违反竞赛规则的，如果引用别人的成果或其他 公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正 文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反 竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B中选择一项填写）：A我们的队号为：53参赛队员：1.钱燕兰吴文倩吴晗指导教师或指导教师组负责人：数模组日期：年月日评阅编号（由评阅老师

2、评阅前进行编号）数学建模竞赛编号专用页评阅编号:评阅记录:评阅人评分备注国内成品油调价与定价机制模型的建立与求解摘要随着我国石油消费量不断增长，石油对外依存度不断上升，国内的成品油价 格成了广大群众十分关心的话题。本文主要针对我国的成品油调价和定价机制问 题，依次建立了模型，分析了现行的成品油定价机制，指出了其中的不合理之处， 并对制定新的定价机制给出了我们的建议。对于问题一，首先我们通过EXCEL简单绘制出了三地原油期货价格与杭州 93号汽油价格随时间变动的走势图，简要分析了三者与杭州93号汽油的关系， 然后通过SPSS作出了杭州93号汽油价格分别与三地原油期货价格的一元线性关 系，再利用S

3、PSS的多元回归分析求出杭州93号汽油价格与三地原油价格总的函 数关系式。为了分析原有定价机制的不合理之处，我们引入了 LOGIT模型，验证 了 22天为调价周期的滞后性，进而深入分析了成品油现行定价机制的弊端。对于问题二，我们首先选取了几个可能对国内油价变动有影响的几个因素， 通过主成分因子分析，剔除了影响较小的因子，再通过MATLAB进行一元线性回 归分析得出各影响因素对目标函数成品油成品油价的影响权重，最终得出成品油 价与各主要影响因素的总关系。再用MATLAB拟合出主要的影响因素随时间变动 的函数关系，最后利用穷举法得到最佳的调价周期11天，表示出了价格的调节 幅度。基于上述模型，我们

4、提出的建议是将原本22天的调价周期缩短为11天，根 据主要影响因素的变化量确定调价幅度。最后，我们也对模型进行了检验与评价。关键词：SPSS LOGIT模型主成分因子分析多元线性回归MATLAB问题背景、成品油（汽油、柴油、煤油等）与国家的经济发展，人民的生产生活直接相 关，成品油的合理定价对国家经济的发展以及社会和谐稳定具有重要的意义。中 国成品油市场运行机制先后经历了完全计划经济阶段、双轨价格过渡阶段、与国 际油价间接接轨阶段等多个主要阶段。目前我国实行的是2009年出台的成品 油价格管理办法（试行）。根据有关统计数据，自2009年以来，国内成品油价格共调整了 17次，其中 12次上调，5

5、次下调，但是，国内成品油价格一直居高不下，甚至在国际原油价 格大幅下跌的时候，国内成品油价格却出现了不降反涨的反常局面，引发人们对 现行油价机制议论纷纷，议论的焦点在于国内成品油价涨多跌少，与国际油价变 化不合拍。现行的油价体制遭到了广泛质疑，人们争相呼吁新的定价体制的出台。但是，要提出更为合理的成品油定价机制，首先必须明确现有定价机制中存 在哪些不足，通过改进这些不足才能达到优化定价机制的目的，但是定价机制是 由许许多多的因素决定的，将所有对油价变动有影响的因素都考虑进去是不可能 的。因此，我们在这里主要通过数学建模的思想，把复杂问题具体化，从中选取 主要的影响因素。通过分析这些主要的影响因

6、素分别对成品油定价机制的影响， 提出更能让大众接受的定价措施。二、问题重述汽柴油价格一直广大群众关心的问题，而我国关于油价的定价机制其实一直 并不是很明确，国务院公布的成品油定价机制简单的来说就是：1、当国际市场三地原油价格连续 22个工作日移动平均价格变化超过 4%时，可相应调整国内成品油价格。2、当国际市场原油价格低于每桶80美元时，按正常加工利润率计算成 品油价格。高于每桶80美元时，开始扣减加工利润率，直至按加工零利润 计算成品油价格。高于每桶130美元时，按照兼顾生产者、消费者利益， 保持国民经济平稳运行的原则，采取适当财税政策保证成品油生产和供应，汽、柴油价格原则上不提或少提。附件

7、中是三地原油期货价格的日数据（2007年1月2011年12月），总 共五年的数据，根据这五年的数据，请回答以下问题：解答下面问题：查询杭州市93号汽油在这五年中的数据，最终给出杭州市93号汽油价格和三地原油期货价格之间的函数关系。并分析此定价机制的合理性。由于国际原油价格在22个工作日内经常会变动很大，现在大家都觉得22个工作日太长，应该缩短，那么考虑到原油价格波动性以及 实际成品油的滞后性（现在所卖成品油其实是若十天以前购入的原油）， 你觉得定价和调价机制如何制定是比较合理的。三、模型假设1、搜集到的数据全部有效。2、无投机抛售等利用市场漏洞的现象。3、政府干预没有对成品油价产生影响。4、假

8、设主要因素对成品油价格的影响占足够大的比例，且这些因素与成品油价格之间均为线性关系。5、假设全国成品油油价受主要因素的影响程度与各地区相同，根据全国成 品油价格所建立的新的定价机制适用于全国各地区。四、符号说明为了便于描述问题，我们用一些符号来代替问题中涉及的一些基本变量，如 表1所示。其他一些变量将在文中陆续说明。符号意义y杭州93号汽油价格pWWTI的原油期货价格pBBrent的原油期货价格POOPEC的原油期货价格Pi成品油价调整概率x1国际原油价格d虚拟变量X1国际每季度原油价格向量X2中国每季度DGP向量X3中国每季度CPI向量X4美元汇率构成的向量Y线性变换形成的新综合变量R相关系

9、数矩阵P1国内成品油价格x2国内生产总值x3国内消费价格指数（表1：符号说明表）五、模型建立与求解5.1问题一的模型5.1.1问题一模型的分析问题一首先让我们根据题中所给2007年1月到2011年12月的三地原油 期货价格以及自行搜集到的杭州市这些年内93号汽油的价格作出函数关系，但 是由于题中所给的三地原油期货价格五年的数据数量较多，所以我们事先对数据 进行了初步的处理，我们通过季度平均求出每年内四个季度三地原油的平均油 价，把原本以日为单位的数据转化成了以季度为单位的数据，2007年至2011年 共20个季度，然后通过查阅数据得到了杭州2007年至2011年的每年每个季度 的平均值，然而由

10、于国内油价与国际油价的定价单位不同，国内是以元/升计量， 而国外是以美元/桶计算，所以我们首先将国内价格按美元汇率以及桶与升之间 的转换关系进行了单位替换，使杭州市的93号汽油价格计价单位与国际相同， 具体数据见表2：年份.季 度WTK美元/桶）Brent （美元/桶）OPEC （美元/桶）杭州93汽油换算后 价格（美元/桶）杭州93汽油换算 前价格（元/升）美元汇率2007.158.0258.5154.3294.814.437.782007.264.9568.6064.6896.114.527.732007.375.1574.4371.1298.664.677.682007.490.7688

11、.3784.63101.384.837.632008.197.7596.2292.5299.575.027.212008.2123.90122.20116.35109.815.637.092008.3118.67117.79115.09105.525.4872008.458.2257.9754.3094.014.916.962009.142.9645.6942.4492.574.926.842009.259.5259.5057.93105.035.596.832009.368.1168.8867.74108.415.776.832009.476.0375.4273.92117.996.286.

12、832010.178.7377.3177.31122.696.536.832010.277.8679.6076.74119.316.356.832010.375.9876.8573.64122.336.536.812010.485.2287.2683.50126.456.786.782011.193.88105.16100.59127.997.066.592011.2102.15117.17112.37134.627.466.562011.389.53112.24108.53129.657.246.512011.494.12108.94107.80128.507.226.47（表2：三地原油期

13、货价格与杭州93汽油价格）【1】5.1.2问题一模型的建立与求解5.1.2.1三地油价与杭州93汽油价总体趋势对比：将三地原油期货价格与杭州市93号汽油的价格通过EXCEL绘制出油价关于 时间的走势图如图一所示：杭州93号汽油价与国际三地油价的对比桶4兀美价油150.00100.0050.000.00t-WTI BrentOPEC*杭州93汽油135791113 15 17 19季度（图1：三地油价与杭州93号汽油价格走势图）从图中我们可以发现三地原油期货价格随时间的变动幅度以及方向几乎是 一致的，但是杭州93号汽油的价格的走势却与三地相差很大。5.1.2.2杭州93号汽油价分别与三地原油期货

14、价格的关系为了分析它与三地油价之间的具体关系，我们下面分别分析三地油价与杭州 93号汽油之间的关系，首先通过SPSS【2】分别做出杭州93号汽油价格随着三 地原油价格的变化关系图如下所示：90.00-140.00=130.00-120.00-110.00100.00-IIIIHI40.0060.0080.00100.00120.00140.00WTI原油期货价格（图2：杭州93号汽油价格与WTI油价关系图）模型汇总植刑RR方调整R方标淮警的误更改统计里R方更改F更改df1df2Sig. F更改1.445 日.198.1531279037.1984.445118.049a.预删寰里:禅里），州1

15、原油侑格。（表3：杭州93号汽油价格与WTI油价的相关性分析）系折#苜刑非标准化系数标淮系数tSig.相关性B标准误茅试用场零盼偏部分1mi原油价格37.150.30212.023.143.4457.2492.108.000.049.445.445.445a.因变重:杭州犯号汽油侑格（表4：杭州93号汽油价格与WTI油价的函数关系）由图2我们可以看出杭州93号汽油价格与WTI的具有线性关系，相关系数R2 = 0.198，相关性系数相对来讲还是比较小的，说明WTI的原油期货价格对杭 州93号汽油价格的影响不是很大。从表3中我们看出sig是小于0.05的，说明 拟合出的线性回归方程是有效的，同时我

16、们根据表4中的数据得到杭州93号汽 油价格与WTI原油期货价格的线性方程为：j = 0.302七+ 87.15（ Pw表示WTI的 原油期货价格）。（图3：杭州93号汽油价格与Brent油价关系图）模型汇总RR方调整R方标准警的误审改瑜计单R方更改F更改dfldf2Sig.F更改1.593.351.31 511.50399.3519.745110.006a.预删安里:慌幽,Brent原油侑格。（表5：杭州93号汽油价格与Brent油价的相关性分析）系数日非标淮佛菜数标淮菜数tSig.相关伎B标淮误差试用版零盼偏部分1常里）Brent原油价格81.083.3611 0.1 61.116.5937

17、.9803.122.000.006.593.593.593a.因变量:杭州弱号汽油侑格（表6：杭州93号汽油价格与Brent油价的函数关系）从图3拟合出杭州93号汽油的价格与Brent的原油期货价格函数关系为 j = 0.361P + 81.083（七代表Brent的原油价格），相关性系数R2 = 0.351，这表 明Brent的原油期货价格对杭州93汽油价格的影响比WTI大。（图4：杭州93号汽油价格与OPEC油价关系图）模型汇总植刑RR方调整R方标淮警的误更改统计里R方更改F更改dt1df2Sig. F更改1.61 V.373.33811.31083.37310.701118.004a.预

18、孤变重:禅里）,QFEC原油价格。（表7：杭州93号汽油价格与OPEC油价的相关性分析）系部非标淮化至数标淮嘉数tSig.相关伎B标准误差试用版零盼偏部分1禅里） FEC原油价格80.753.3799.814.116.6118.2293.271.000.004.611.611.611a.因变里:杭州珀号汽油侑格（表8：杭州93号汽油价格与OPEC油价的函数关系）由图4拟合得到的函数关系为j = 0.379p + 80.753（ P代表OPEC原油期货 价格），相关性系数R2 = 0.373，这个数据相比Brent又大了一点，由三幅图拟 合的结果可以看出，对杭州市93号汽油价格影响最大的是OPE

19、C，其次是Brent， 最小的是WTI。5.1.2.3杭州93号汽油价与三地原油期货价格的多元线性关系在前面我们通过SPSS拟合了分别以三地原油期货价格作为自变量，杭州93 号汽油价格作为因变量的简单线性关系，但是这样拟合出的函数关系并不能很好 的表示杭州93号汽油价格与国际原油期货价格之间的关系，三地的原油期货价 格不能代表国际的原油价格，所以，我们继续利用SPSS对所有数据进行多元线 性回归，以WTI、Brent、OPEC的原油期货价格为自变量，杭州93号汽油价格作 为因变量，作出多元回归线性方程，具体如下图所示：模型汇总RR方调整R方标准管的误更改统计里R方更改F更改df1df2Sig.

20、F更改1782a.611.5389.44755.611S.379316.001a.预删变里:岸里）,OPEC原油伯格,WIT原油价格,Brent原油侑格。（表9：杭州93号汽油价格与三地原油价格的相关性分析）系粉但刑非标准化菜数标淮系数tSiq.相关性B标淮误差试用场偏部升1岸里）93.5749.1781 0.1 95.000如IT原油侑格-.911.354-1.343-2.573.020.445-.541-.401Brent原油侑格-1.9851.602-3.255-1.239.233.593-.296-.193PEC原油价格3.1921.6005.1391.996.063.611.446.

21、311a.因变里:杭州93号汽油侑格（表10：杭州93号汽油价格与三地原油价格的函数关系）从表10中所得数据我们可以列出杭州93号汽油价格与三地原油价格的多元 线性函数关系，j = 0.911。1.985。+3.192%+ 93.574，sig 的数值是 0.001， 是小于0.05的，说明拟合得到的方程式是有效地，并且R2 = 0.611，相关系数还 是很大的。5.1.3合理性分析通过以上的图表以及函数关系分析，我们可以概括出国内的成品油定价仍然 存在以下不合理之处：一、数学模型决定易涨难跌自2007年新调价体制实施后，发改委一共17次调整油价，其中上调次数 12次，单从数学模型来看，假设原

22、油价格从100美元上涨到104美元，涨幅4%， 我们将上调价格，但如果从104美元跌到100美元，跌幅不足4%，不具备下调 条件，再从100美元涨到104美元，涨幅又达4%，再度涨价。这充分说明，单 纯用4%公式来决定涨跌，必然会造成涨多跌少，引起非议。据卓创资讯能源行 业分析师朱春凯统计，2010年8月6日和2011年8月8日的WTI原油期货价 格分别为每桶80. 70美元和81.31美元，二者间相差不到1美元。同时间国内 93号汽油均价分别为每吨7717元和9306元;柴油均价分别为每吨6736元和8307 元。在国际原油价格相差无几的情况下，国内却出现了汽油每吨1589元、柴油 每吨15

23、71元的差价。布伦特原油期货价格与国内成品油价格对比也呈类似情况。据朱春凯统计， 2008年9月24日、2011年2月17日、2011年8月9日，布伦特原油期货价格 分别为每桶102.45美元、102.59美元、102.57美元。同时间国内93号汽油均 价分别为每吨7750元、8681元和9306元，柴油均价分别为每吨7300元、8069 元和8307元。不难看出，今年2月17日与8月9日国际原油价格仅差0.02美 元，但不到半年的时间国内汽油价格每吨上涨625元，柴油价格每吨上涨238 元。相同国际原油水平下，成品油价格涨幅依然明显。二、“三地”参考标准存在质疑我们通过对比2012年2月与2

24、008年7月的国际原油价格与国内成品油价格 发现，在这两个时间点，WTI与Brent原油价格分别下跌25.25%和12.8%，但是 同期的国内汽柴油价格却分别上涨22.9%和8.39%，其中柴油价格甚至涨至国内 历史最高。从价格看来，原料成本表现为下降，而产出产品价格却表现为上升。我国现行成品油定价机制采用的是迪拜、辛塔、布伦特三地原油现货价格为 基准，但是每年我国进口的上述三种原油数量在总进口额的占比都比较小。尽管 迪拜、辛塔、布伦特都代表着世界各地油价，但是对我国成品油价格影响较小， 通过2011年国内部分原油进口份额来看，从沙特进口 了 33.1亿元的原油，位列 第一，其次是安哥拉，24

25、.4亿元，第三是伊朗为19.4亿元。相比于这些进口大国，去年我国从迪拜进口的原油只有67万吨，仅占2.67%， 辛塔和布伦特占比更是微乎其微，辛塔只有0.26%，布伦特dtd甚至是零进口。三、限定的22天调价周期太长不少人指出，正是由于22个工作日的长周期外加4%的大幅度，导致了自今 年6月1日成品油调整价格之后，虽然国际原油价格上下波动，但是国内油价却 4个多月未进行任何调整，其间，部分地区甚至出现了油荒问题。我们从图1可 以看出，三地的油价还是比较能代表国际原油价格的，波动幅度几乎一致，而杭 州近五年的93号汽油的价格波动幅度与国际油价相差差距还是很大的，从图中 我们可以看出当国际原油降价

26、时，杭州市93号汽油的相对于国际原油价格的滞 后时间约为33天，可见，中国现行的定价体制对国际原油价格变化的灵敏度较 低。可见，22天不能及时反映国际市场石油供求情况，由上面的图形进行分析 我们也得到从原油到成品油存在一个月左右的滞后期，所以原先设定的22天会 导致国内油价与国际原油动态不匹配。为了进一步分析国内成品油定价机制存在的缺陷，我们通过LOGIT模型【3】 进行条件概率估计，找出了现有成品油定价机制存在的调价周期滞后，调价难以 反映国内市场供求关系及消费结构变化等缺陷，从而为建立合理的成品油定价机 制模型提供参考。以2007年3月至2008年5月的国内成品油价格（这里暂时将杭州市93

27、号 汽油当做国内成品油价格）是否调整为因变量，设当成品油价格未调整时为 0，当成品油价格调整时d为1，以国际原油（三地原油期货市场平均价格）为 自变量构造线性回归模型（LRM）。通过查阅历年我国成品油价格调整表发现在 2007年3月至2008年5月间成品油价格仅有一次调整，因此借用LOGIT模型进 行条件概率估计，具体分析如下：1、P = E（Y = 11 X ） = 1 （1）户代表概率，B、B为常数。i12上式也可以改写为：i 1 + e -Z1 + eZi对式（2）做线性变换：P _ 1 + ezt _ ”，T-P 1- ez eei 对式（3）取对数，得到：P 、L =ln（ ） =

28、Z = B + B X i1 - Pi 12 ii利用LOGIT模型的上述形式，即：- -，d 、 一 一-L = In（厂j）= Z. = BOY七表示国际原油价格。（2）（3）（4）（5）利用式（5）进行回归分析，经eviews6.0分析结果如表3所示:变项系数标准差P值0Y1530.26-5.860.0000Akaik系敬0.0559Schwarz 系数0.0678（表11：国内成品油价调整概率与国际原油价格LOGIT分析结果）表中是2007年3月至2008年5月国内成品油价调整概率与国际原油价格 LOGIT的分析结果，下面的曲线图则反映了在这期间国内成品油价格调整的概率 的预测值。20

29、0 ()22m? fl32DO? 0421X)7 0!（图5：国内成品油价格调整概率预测值）从模型分析结果来看，变量的Z检验结果显著，匕与国际成品油价格变化 呈负相关，在这一时期，虽然国际原油价格不断上涨,但我国国内成品油价格上调 的概率却在不断下降，对d做概率预测发现，2007年11月的调整概率仅有 0.002。由此可以看出，按照现行的定价机制，国家发改委以22天为一个定价周期 是不合理的，国内价格的调整与国际市场变化相比至少滞后22天，这样导致的 油价迟滞主要有以下两方面的弊端：一、调价迟滞难以反映国内市场供求关系及 消费结构变化。二、由于价格调整与国际市场变化滞后过多，调价不能及时反映

30、国际原油价格影响。四、盲目跟随国际定价某种程度上，国际定价是具有一定代表性的，但是，在市场经济条件下，价 格随供求关系的变化而变化，我国当前的成品油定价机制参考的是国际石油价 格，而没有反映国内的成品油供求关系，没有考虑国内的消费现状，所以不能反 映我国市场上的石油供需，价格失真成分较大。5.2问题二模型的分析、建立与求解5.2.1模型一主成分分析5.2.1.1模型的分析影响油价的因素有很多，不同专家研究的结果也不尽相同。本文主要采用了 普遍认可的几个影响油价的指标来进行分析。主要研究国内成品油价与国际原油 价格、本国GDP、CPI以及美元汇率之间的关系。根据相关文献制得2007年至2011年

31、每季度各指标数据表如下【4】：时间（年.季度）国内（杭州） 油价（美元/ 桶）国际原油价 格（美元/桶）GDP （亿美元/季度）CPI美元汇率2007.194.8156.717038.05102.437.782007.296.1166.237922.77102.977.732007.398.6673.248346.61103.837.682007.4101.3887.3711233.16104.67.632008.199.5795.059193.34107.677.212008.2109.81120.2010464.60108.077.092008.3105.52116.8910935.431

32、07.3372008.494.0156.5413939.51106.26.962009.192.5743.799611.84100.036.842009.2105.0358.8710851.6899.076.832009.3108.4168.2511413.6298.836.832009.4117.9974.9418021.3899.16.832010.1122.6977.6111950.56101.936.832010.2119.3178.0713355.45102.56.832010.3122.3375.3714070.48102.86.812010.4126.4585.3019074.1

33、9103.176.782011.1127.99100.8914614.72104.936.592011.2134.62112.0016485.98105.236.562011.3129.65105.8617935.50105.66.512011.4128.50105.2923237.19105.56.47（表12： 2007年至2011年每季度各指标数据）（注：国际原油价格是按照模型一中三地原油价格对杭州汽油价格的影响程 度计算相应的权数，从而三者加权平均得到的。）在考虑到变量个数太多会使得问题复杂化，无法更好地分析出各变量之间的 关系，所以我们尽可能使用较少的变量但是又能保全信息的较好完整性

34、，反映它 原有的信息。此模型通过主成分分析法来找出具有代表性的因素。主成分分析法【5】是对数学降维的一种方法。其基本思想是设法将原来众 多的具有一定相关性的p个指标，重新组合成一组新的互不相关的综合指标来代 替原指标。将原来p个指标做线性组合，作为新的综合指标。在所有的线性组合 中选取的F应是方差最大的，作为第一主成分。如果第一主成分不足以代表原1来p个指标的信息，再考虑选取F作为第二个线性组合。此时Ca（ FF） 0，称 F为第二主成分，以此类推。5.2.1.2模型的建立本模型涉及p = 4个指标，分别用X1X X X 4表示。4个指标构成的向量.，n = 20, i = 1, 2,3, 4

35、 p = 4ni ，p.X=(X , X,X,X) X = (x ,x ,.,x ) 为1 2 3 4 i1i 2i* 1 国际每季度原油价格向量X 2中国每季度DGP向量X 3中国每季度CPI向量X 4 美元汇率构成的向量Y对X进行线性变换，形成的新综合变量(1)新的综合变量可以由原来的变量线性表示，即满足下式:Y=u X+ u X1111122Y=u X+u X2211222Y=u X+u Xpp11p22 +/,组合系数满足 u 2 + u 2 +. + u 2 = 1组合系数满足：11i 2ip通过计算协方差矩阵、特征值以及相应的特征向量，从而选择主成分。5.2.1.3模型的求解将没一

36、个变量单独作为线性组合，观察每个变量的方差贡献率，由此来得到 每个变量所包含的信息量的多少。利用SPSS软件进行操作得到相关系数矩阵以及各变量的方差贡献率：(具体 输出结果见附录1)r 1.0000.393-.3270.718、0.3931.000-0.7550.075-0.327-0.7551.0000.043 0.7180.0750.0431.000/变量名称国际原油价格GDPCPI美元汇率方差贡献率0.53740.35290.06170.0480累计贡献率0.53740.89020.95201（表13:标量的方差贡献率表）从以上表格中可以看出，前三个变量的累计贡献率已经达到了 95.2%

37、,美元 汇率的贡献率很低，因此前三个变量已经包含了较多信息，这里选取前三个变量 国际原油价格、GDP、CPI作为影响国内成品油价格的主要因素。这既简化了所 求问题同时也能为模型二成品油的定价奠定数据基础。5.2.2模型二成品油定价机制模型5.2.2.1模型的分析22个工作日太长，为了能适当缩短工作日，需要制定更优的国内成品油定 价机制，而更优的定价模型需要考虑从多方面的因素保证其合理性。本模型在模 型一主成分分析结果的基础上首先将国内成品油价格与国际原油价格、国内成产 总值、国内消费价格指数分别进行一元拟合，观察其趋势与相关程度。再对其进 行多元拟合，建立合理的函数，并用最小二乘法确定函数中的

38、系数，最后再与实 际比较验证其是否优化。5.2.2.2模型的建立为了更能直观地观察国际原油价格、GDP、CPI对国内成品油价格的影响大 小，本文利用SPSS软件画出了各变量与国内成品油价格的散点图，并用线性函 数进行一元线性回归分析。见图6、图7、图8：（图6:国内成油品价格与国际原油价格的一次回归图） 从上图中看出国内成油品价格与国际原油价格的存在着正相关关系，R 2 0.331数字较大，存在的较密切的相关性。140.00130.00-120.00=110.00-100.00=90.005000.0010000.0015000.0020000.0025000.00GDP（图7:国内成油品价格

39、与GDP的一次回归图）从上图中看出国内成油品价格与 GDP的存在着明显的正相关关系，R 2 0.623说明两者有很强的相关性。140.00-130.00-120.00110.00-100.00-90.00-I 98.00100.00102.00104.00106.00108.00110.00CPI（图8:国内成油品价格与CPI的一次回归图）从上图的趋势以及R2 = 0.007看出国内成油品价格与CPI关系非常弱，但是由于在主成分分析时并未剔除，我们还是将它作为其中一个因素进行接下来的研 究。总之，从以上三幅图中我们可以清晰地看到，国内成品油价格与国际原油价 格、国内生产总值、国内消费价格指数之

40、间具有一定的相关性，其中前两者很强。 鉴于此，我们对此进行多元回归分析。Pi 国内成品油价格X2国内生产总值国内消费价格指数P0P = P +P x +P11，2，七均为参数（2）x + P x01 12 23 35.2.2.2模型的求解在MATLAB中输入相应程序程序，MATLAB输出结果整理后得到如下表：（具体程序以及输出结果见附录1）：参数P0P1P2g3估计值279.72170.42940.0018-2.1945（表14:参数估计值表）所以得到表达式为：P1 = 279.7217 + 0.4294七 + 0.0018%-2.1945%5.2.2.3模型的检验计算值与实际值对比检验运用M

41、ATLAB我们得到了国内成品油价格的表达式，为了确定此表达式的和 合理度，我们将每季度的各因素数据一一代入总而得到计算后的成品油价格，将其与实际国内成品油价格进行比较，得到如下表格:时间（年.季度）实际国内石油季度平均 价格（美元/桶）计算国内石油季度平均 价格（美元/桶）2007.194.8191.958832007.296.1196.454182007.398.6698.339922007.4101.38107.91342008.199.57100.80242008.2109.81113.01222008.3105.52114.06242008.494.0196.035192009.192

42、.5796.31062009.2105.03107.12442009.3108.41112.69032009.4117.99126.86452010.1122.69110.87312010.2119.31112.34852010.3122.33111.81782010.4126.45124.27652011.1127.99119.08152011.2134.6291.958832011.3129.6596.454182011.4128.5098.33992为更直观地找出两者的契合程度，利用MATLAB软件，由实际国内石油季度 平均价格和计算国内石油季度平均价格画出对比图如下（具体MATLA程序

43、输入见 附录2）：（图9:实际国内石油季度平均价格和计算国内石油季度平均价格对比图）（注：红色代表实际值，另一个颜色代表估计值，2007.1在横坐标以1表示， 2007.2在横坐标以2表示，以此类推。）从以上图中可知：红绿色曲线虽不是完全相同，但是趋势基本一致，这说明 上述定价机制的表达式存在一定的合理性。对残差和置信区间的指标进行残差分析检验进一步对国内成品油价格表达式进行检验，利用MATLAB得到得到一些相关数值表格以及残差分析图（输入MATLAB的程序以及具体输出结果见附录2）.表格如下：回归系数回归系数估计值回归字数置信区间P0279.7217(105.1256 454.3177)P1

44、0.42940.17790.6809)P20.00180.00090.0027)g3-2.1945(-3.9885 -0.4005)R2 =0.7929 F = 20.4212 p = 0.00005 2 =47.5233（表15:相关数值整理表）从表中可得，四个置信区间全是符合要求的不存在零点，且R 2 =0.7929 这表明了多元回归的可行性很好。用MATLAB得到残差分析结果如下图:Residual Ca&e Order Plot2468101214161820Casg Number（图10：残差分析结果）从此图我们看到：20个数据符合要求间，因此模型具有合理性，此模型是 可以应用的。综

45、上所述，经过以上计算值与实际值对比检验和对残差和置信区间的指标进 行残差分析检验，我们可以确定此模型的可行性，最终建立我国国内成品油价格 与国际原油价格、国内成产总值以及国内消费价格指数的多元线性回归模型为：P = 55.5479 + 0.1241x - 0.0018x + 0.6693x原先模型一中只按照国际原油价格来定国内成品油价格误差极大，模型二除 国际原油价格之外我们考虑了其他较为重要的因素这里主要是我国的 GDP和 CPI。加入这两个因素拟合出国内成品油价格与实际较为接近也表明了此定价机 制较于之前的确有优化。因此作为改进的定价机制。5.2.2模型三逐个试验法确定调价机制为了确定调价

46、机制，制定合理的调整周期，我们首先得确立三因素与时间之 间存在的函数关系，然后利用调价机制的模型来确定最终合适的调整周期。【6】 5.2.2.1确定三因素与时间的函数关系由上述模型二可知国内成品油价格主要是由国际原油价格、GDP、CPI作用 的，先为了寻找合适的国内成品油价格的调节周期，我们需先从三个因素着手研 究，先用MATLAB分别对国际原油价格、GDP、CPI关于时间进行曲线模拟。首先曲线模拟国际原油价格关于时间t（ t是以每季度为单位，从2007年 至2011年共20个季度）的变化曲线，（SPSS的输出结果具体见附录3），得到图 如下：国际原油价格时mi季度（图ii：国际原油价格关于时

47、间t（t是以每季度为单位）的变化曲线） 根据结果得如下表格：参数R 2b0b1b2b3数值0.36251.66613.478-1.6750.059（表16:相关参数整理表）模拟得到曲线为：尤=13.478t3 - 1.675t2 + 0.059尤 + 51.6661同理，曲线模拟GDP关于时间（t是以每季度为单位）的变化曲线，（MATLAB 中输入的程序以及输出结果具体见附录4），得到图如下：（图12:GDP关于时间t （t是以每季度为单位）的变化曲线） 根据结果得如下表格：参数R 2b0b1数值0.7426527.073615.022（表17:相关参数整理表）GDP时酣季度o U观耐我性模拟

48、得到曲线为:七=615.022t + 6527.073同理，曲线模拟CPI关于时间t（ t是以每季度为单位）的变化曲线，（EXCEL 中输入的程序以及输出结果具体见附录5），得到图如下：(图13： GDP关于时间t (t是以每季度为单位)的变化曲线)模型拟合得到函数关系为：X = 0.0019t4 + 0.0881t3 - 1.3475t2 + 7.2153t + 94.432 35.2.2.2逐个试验寻求调节周期由于国际原油价格在22个工作日内经常会变动很大，结合我国的具体 情况我们了解到很多专家都认为，22个工作日作为调节周期太长，应缩短 调节周期，大部分认为调节周期应为一周至两周之间，即

49、 7至15天左右。根据模型二中已求得的三个因素的参数 P , P , P，计算周期： 123寸(P x +P x +P X )/T= 0.041 1n2 2 n3 3 nn=1逐个试验，因为模拟曲线时以每季度为单位为，所以需将调节周期换算乘 季度来试验。T = 9/90 = 0.1，求得，连续9天移动平均价格超过0.04的个数为得到表格如下:T9101112131415应调节 次数58119968同理类推可知，符合条件的是11天，因此我们确定调节周期T为11天。调节幅度为d (P1(x )人 公 d (P1(x )人 公 d (P1(x )人AP1 =px ，Ax +P x 二Ax +px 匕

50、Ax1d(x )12 d (x )2 r3 d (X )3123六、模型评价模型的优点1、在第一个问题中用数值分析的方法客观的给出了杭州93汽油油价与国际三地 原油期货价格的直观关系，便于我们对现有定价体制进行分析。2、在建立较现行成品油定价机制更为合理的模型之前，首先借助LOGIT模型进 行了概率估计，进一步明确现有成品油定价模型存在的缺陷，为建立更合理的定 价机制奠定了基础。3、通过主成分因子分析法找出了对国内成品油价格影响较大的几个因素，抓住 主要矛盾，便于找出各因子与国内成品油价的关系和内在规律性，再通过MATLAB 求出各自所占的权重，准确表示它们之间的关系。4、由于新建立的成品油定

51、价机制中引入了一些反映国内实际情况的因素，这使 得国内成品油的定价较合理。5、通过计算得出的国内成品油价的计算公式，利用穷举法确定出了最合理的调 节周期以及调节幅度。模型的缺点1、所搜集到的数据不足量，对所建模型的准确性有一定影响。2、由于我们在定价时选取的是权重的较大的因素，而省略了其他对成品油价格 有影响的因素，使油价过分的依赖某些变量，突然性较大，可能导致模型的实用 性不强。3、各个因素的权重求取方法仍有很大的改进空间，而不是简单地运用MATLAB 自动计算，还可以配合其他分析方法，使权重更加合理。七、参考文献【1】浙江省93#汽油历次价格调整一览【2】卢纹岱，SPSS统计分析(第四版)，电子工业出版社，2010年4月【3】段龙龙，中国成品油定价机制：基于计量模型的一个解释，科学决策，2011 年4月【4】百度文库(GDP、CPI、美元汇率)【5】姜启源，谢金星，叶俊，数学建模(第三版)，高等教育出版社，2003【6】程启月，评测指标权重确定的结果熵权法，系统工程理论与实践，2010年 第07期八、附录附录1多元回归分析输入MATLAB D的程序：x1=56.71 66.23 73.24 87.37 95.05 120.20 116.89 56.