货运量预测是制定合理物流计划的重要内容。本文以A物流公司的历史货运量为例,主要基于趋势外推法和季节分解法来进行研究,建立趋势预测单项模型和季节分解-趋势外推的组合模型。探讨对时间序列的短期趋势预测是直接曲线拟合还是“剔除”季节性后再拟合,对于趋势模型的选择,对比多项式拟合和Logistic曲线拟合的效果,同时假设季节因素对每个月的影响是不变的,最终构建四个预测模型,通过比较四种模型的拟合优度来评判出合适的预测模型。预测结果证明组合模型的拟合效果要比单项模型的要好。 Freight volume forecasting is an important part of the reasonable logistics planning. In this paper, a historical freight volume of A company is regarded as research data, mainly basing on the trend extrapolation method and the seasonal decomposition method, establishing a single model and a combination model. Discussion on the time series of short-term forecasting is a direct curve or “remove” seasonal before fitting. For the trend model, there are two choices: One is polynomial model; another is Logistic model. While assuming the influence of monthly seasonal factors is same. We could build four predictive models ultimately, compare the goodness of fitting, then chose the best one. The result demonstrates that the combination model is better than the single model.
罗维,方晓平
中南大学交通运输工程学院,湖南 长沙
收稿日期:2016年2月19日;录用日期:2016年3月5日;发布日期:2016年3月9日
货运量预测是制定合理物流计划的重要内容。本文以A物流公司的历史货运量为例,主要基于趋势外推法和季节分解法来进行研究,建立趋势预测单项模型和季节分解-趋势外推的组合模型。探讨对时间序列的短期趋势预测是直接曲线拟合还是“剔除”季节性后再拟合,对于趋势模型的选择,对比多项式拟合和Logistic曲线拟合的效果,同时假设季节因素对每个月的影响是不变的,最终构建四个预测模型,通过比较四种模型的拟合优度来评判出合适的预测模型。预测结果证明组合模型的拟合效果要比单项模型的要好。
关键词 :趋势外推法,季节分解法,Logistic曲线,货运量预测
货运量是指在一定时期内运输企业实际产生的货物的数量。对物流企业来说,一定时期内的货运量的大小能反映企业的运营的好坏。对货运量进行预测是物流企业经营决策的重要前提,通过合理预测,物流企业还可以有效利用各种资源来降低成本,在满足顾客需求的前提下获得最大利润。所以,货运量预测对物流企业是十分重要的。
据统计,到目前为止,国内外应用于货运量预测的方法多达300多种。主要有时间序列预测法、灰色预测法、回归分析法、神经网络法和组合预测法。从相关性研究方面来看,Russo F在考虑了供需关系、经济政策等多种影响因素的情况下,结合宏观经济以及运输模式建立了意大利综合货运需求模型[
研究货运量预测的文献较多的是从国家、区域和交通通道等宏观层面,针对物流企业进行货运量预测的较少[
本文主要关注时间序列预测法中的趋势外推模型和季节性波动分解模型的组合模型。时间序列是指将某一事物发展的数值按其发生的时间的先后顺序排列形成的数列。时间序列预测法即根据历史数据对未来进行预测。任何时间序列从要素上可以被划分为四个部分:趋势(T)、周期(C)、季节性(S)和随机元素(R)。这四种成分的随意组合,可以将季节性时间序列分为加法模型和乘法模型[
本文不考虑经济长周期,而只考虑年为周期的正常经济活动,故季节性也体现在该周期内。公式(1)和(2)里面的C与S就合二为一了。
在对社会统计资料进行分析时可以发现,经济现象总是随着时间的变化呈现出某种上升或是下降的趋势,并且没有明显的季节波动,那么我们可以通过一条曲线来描述这种变化。再对这条曲线建立相应的数学模型,即趋势预测模型。在本文中主要讨论多项式模型和Logistic模型。
(1) 多项式模型:
(2) Logistic模型:
季节分解预测法又称为季节因素预测法或是季节指数预测法。系统内部的周期性运动、变化规律和系统外部的环境的季节性作用共同导致了季节性的波动。季节分解预测法能够既考虑到季节指数的趋势变化,还能充分发挥已知数据的信息价值。季节分解预测法的思路分为三个步骤:第一步,“剔除”季节性并预测长期趋势值;第二步,计算季节指数;第三步,通过方法让季节性在趋势上反应出来,达到预测的效果。
首先“剔除”趋势,在选择究竟是加法还是乘法模型时,高孝伟认为加法模型比乘法模型更加合理[
考虑到货运量与经济之间的密切关系,而且依据数据的周期特性,确定周期
得到季节性与随机元素的组合
假设不同年度的相同月份的季节影响是一样的,即
即当
按照以上方法,我们求出了趋势模型和季节因子,故我们建立预测模型:
A物流公司主要经营国内的公路零担货物运输业务。目前运输服务的线路基本上覆盖了全国各个省、市、自治区,服务网络遍及全国,已开设直营网点5000多家,自有的营运车辆近9000台,分布在全国的转运中心的总面积超过100万平方米。A公司在2011年的货物吞吐量达到近500万吨。本文数据是A公司自2009年1月至2011年12月共36个月的货运量记录。时间序列记为:
样本的时间序列图如图1所示。从图可以看出货运量随着时间的推移有明显的不断上升的趋势,而且每年的1~2月和6~7月的货运量都是全年最低的,4~5月和10~11月的货运量是一年中最多的,故序列存在着明显的季节波动现象。
本节内容分为两部分,第一部分是建立趋势模型,第二部分是计算季节因子。其中,趋势模型的建立运用趋势外推单项模型和外推模型与季节分解模型的组合模型。单项模型基于多项式和Logistic曲线的直接趋势外推模型,组合模型先运用季节分解模型得到“剔除”了季节性的趋势序列,再进行趋势拟合。
(1) 直接趋势外推建模
A) 用多项式模型分别进行拟合,拟合效果如图2所示。对拟合模型进行参数检验,其中二次、三次模型的拟合度相等,而且拟合优度要高于线性模型。如表1所示。
图1. 时序图
进一步发现二次模型和三次模型的系数t值和相应p值均大于临界值,说明当使用二次和三次模型时系数是不显著的,即统计没有意义。故选用线性方程来建立趋势模型。
模型一:
B) Logistic曲线拟合要输入上限值,故令K取不同值。输出结果如表2所示。
当K = 3700000时R2是最大的,建立Logistic曲线模型为。
模型二:
(2) 季节分解模型与趋势外推模型组合建模
运用Eviews8.0的X-12法进行季节分解。SER01_SA表示季节调整后的序列,SER01_SF表示季节因子序列,SER01_TC表示趋势与循环波动叠加序列,SER01_IR表示不规则的波动。如图3所示。
A) 对SER01_TC进行多项式趋势拟合,结果如表3。
其中二次和三次模型的系数t值和相应p值均大于临界值,故选用线性模型来建立趋势模型。
模型三:
图2. 多项式模型拟合图
方程 | 模型汇总 | 参数估计值 | |||||||
---|---|---|---|---|---|---|---|---|---|
R方 | F | df1 | df2 | Sig. | 常数 | b1 | b2 | b3 | |
线性 | 0.676 | 70.780 | 1 | 34 | 0.000 | 902,948.933 | 62,328.423 | ||
二次 | 0.678 | 34.775 | 2 | 33 | 0.000 | 803,570.406 | 78,019.769 | −424.090 | |
三次 | 0.678 | 22.483 | 3 | 32 | 0.000 | 791,283.313 | 81,751.564 | −672.817 | 4.482 |
表1. 多项式模型汇总和参数检验
K | 3,500,000 | 3,600,000 | 3,700,000 | 3,750,000 | 3,800,000 | 3,900,000 |
---|---|---|---|---|---|---|
R2 | 0.66 | 0.669 | 0.671 | 0.67 | 0.67 | 0.669 |
表2. 不同K值的模型拟合度
图3. 货运量序列季节分解后得到的各个序列
方程 | 模型汇总 | 参数估计值 | |||||||
---|---|---|---|---|---|---|---|---|---|
R 方 | F | df1 | df2 | Sig. | 常数 | b1 | b2 | b3 | |
线性 | 0.923 | 408.095 | 1 | 34 | 0.000 | 1,018,216.812 | 54,055.968 | ||
二次 | 0.935 | 236.900 | 2 | 33 | 0.000 | 863,953.195 | 78,413.382 | −658.308 | |
三次 | 0.976 | 434.120 | 3 | 32 | 0.000 | 1,234,551.298 | −34,143.441 | 6843.677 | −135.171 |
表3. SER01_TC序列多项式模型汇总和参数估计
B) 对SER01_TC序列进行Logistic模型拟合,令K取如下值,输出不同K值时的R2。见表4。
可以发现当K = 2790000,2800000,2810000时R2是最大的,且R2都为0.960。在这里综合实际情况选择建立Logistic曲线模型。
模型四:
根据求季节因子的公式(8),分别求得各个模型的季节因子如表5所示。
根据所求得趋势模型和季节因子,由公式(10)得到预测模型分别为:
模型五:
K | 2730000 | 2750000 | 2770000 | 2780000 | 2790000 | 2800000 | 2810000 |
---|---|---|---|---|---|---|---|
R方 | 0.938 | 0.953 | 0.958 | 0.959 | 0.960 | 0.960 | 0.960 |
表4. SER01_TC序列不同K值模型的拟合度
一月 | 二月 | 三月 | 四月 | 五月 | 六月 | 七月 | 八月 | 九月 | 十月 | 十一月 | 十二月 | |
---|---|---|---|---|---|---|---|---|---|---|---|---|
模型一 | −0.268 | −0.412 | 0.075 | 0.061 | 0.064 | 0.001 | 0.062 | 0.055 | 0.222 | 0.014 | 0.049 | 0.048 |
模型二 | −0.266 | −0.41 | 0.073 | 0.061 | 0.062 | −0.003 | 0.057 | 0.05 | 0.216 | 0.007 | 0.044 | 0.043 |
模型三 | −0.29 | −0.43 | 0.064 | 0.054 | 0.065 | 0.009 | 0.073 | 0.072 | 0.246 | 0.039 | 0.077 | 0.079 |
模型四 | −0.258 | −0.393 | 0.086 | 0.069 | 0.065 | 0.001 | 0.061 | 0.059 | 0.229 | 0.026 | 0.067 | 0.069 |
表5. 各模型的季节因子
直接趋势外推法 | 组合模型 | 文献[ | ||||
---|---|---|---|---|---|---|
多项式模型 | Logistic模型 | 多项式模型 | Logistic模型 | GM(1,1)周期外延组合模型 | SARIMA 模型 | |
(模型五) | (模型六) | (模型七) | (模型八) | |||
R2 | 0.892 | 0.899 | 0.900 | 0.904 | 0.960 | 0.445 |
表6. 预测模型的R2
模型六:
模型七:
模型八:
本文采用
对不同预测模型运用拟合模型利用式进行拟合优度对比分析,计算结果如下:[
由表6知,组合模型的拟合优度要高于直接趋势外推模型和SARIMA模型,但不如GM(1,1)周期外延组合模型。且Logistic趋势模型的拟合度要比多项式模型高。
本文基于趋势外推法和季节分解法建立了货运量单项预测模型和组合预测模型,并对四种预测模型拟合优度进行了横向和纵向比较分析。从纵向看,Logistic曲线对序列趋势的拟合效果要好于多项式拟合,说明本案例中货运量时间序列的发展变化与Logistic曲线变化更相似;从横向看,组合模型的拟合优度要高于单项预测模型,表明在对序列趋势进行拟合之前,“剔除”序列的季节性是有必要的。
湖南省科技计划项目,项目编号2015ZK3001。
罗维,方晓平. 基于季节性时间序列的物流企业货运需求预测研究与应用 Research and Application of Logistics Enterprise Freight Demand Forecasting Based on Seasonal Time Series[J]. 管理科学与工程, 2016, 05(01): 7-14. http://dx.doi.org/10.12677/MSE.2016.51002