时间序列分解法文档格式.docx

1、这45个数据组成的序列我们称之为移动平均数序列，用MA来表示，MATC。2季节性 由于（1） 因此将观察值除以移动平均数得到的比率值就只包含季节性和随机性，从而这些比率包括了确定季节性因素所需要的信息。如果某个比率的值100，意味着实际值X比移动平均数（TC）要大。由于X中包含季节性和随机性，因而当比率值大于100时，就意味着这个季度的季节性和随机性高于平均数。反之，如果比率小于100，则表示季节性和随机性低于平均数。表.2 某产品48个季度的销售数据及数据分解 季度观察值Xt移动平均值TCSI比率长期趋势T循环变动C 13017.602774.81 23043.542813.7732094.

2、352741.33476.3392852.7396.10 42809.842805.632100.1502891.6997.02 53274.802835.569115.4902930.6596.76 63163.282840.558111.3612969.6195.65 72114.312894.24073.0523008.5796.20 83024.572907.411104.0303047.5395.40 93327.482989.961111.2883086.4996.87 103493.483071.367113.7443125.4598.27 113439.933187.92176

3、.5373164.41100.74 123490.793277.322106.5143203.37102.31 133685.083319.258111.0213242.33102.37 143661.233303.883110.8163281.29100.69 152378.433296.07372.1593320.2599.27 163459.553337.209103.6663359.2199.34 173849.633347.198115.0103398.1798.50 183701.183413.185108.4383437.1399.30 192642.383444.67876.7

4、063476.0999.10 203585.523501.936102.3873515.0599.63 214078.663553.405114.7823554.0199.98 223907.063597.425108.6073592.97100.12 232818.463723.42175.6953631.93102.52 244089.503788.657107.9413670.89103.21 254339.613849.043112.7453709.85103.75 264148.603874.540107.1013748.81103.35 272976.453872.32575.31

5、53787.77102.23 284084.643848.029106.1493826.73100.56 294242.423810.274111.3423865.6998.57 303997.583801.414105.1603904.6597.36 312881.013789.31176.0303943.6196.09 324036.233818.788105.6943982.5795.89 334360.333909.526111.5314021.5397.21 344360.533982.320109.4974060.4998.07 353172.184029.20378.730409

6、9.4598.29 364223.764111.740102.7244138.4199.36 374690.484195.228111.8054177.37100.43 384694.484237.770110.7774216.33100.51 393342.354326.23777.2584255.29101.67 404577.634394.982104.1564294.25102.35 414965.464477.872110.8894333.21103.34 425026.054509.818111.4474372.17103.15 433470.144496.89577.167441

7、1.13101.94 444525.944570.21099.0314450.09102.70 455258.714611.094114.0454489.05102.72 465489.584642.750111.7784528.01102.53 473596.764481.66780.2554566.9798.13 483881.604605.93由式（1）可知，如果能将SI中的随机性部分去掉，则就得到了季节性指数。要做到这一点，只需注意到随机性指的是偶然性、没有一定模式、围绕中间值0上下波动。因此通过平均就能去掉随机性的影响。将表4.3中“SI比率”这一栏列成表4.6的形式，将各年同一季度

8、的数据放在同一列之中，求相同各季度的平均值，得第一至第四季度的平均数分别为112.72，109.88，76.28，103.86。由于从1991年至2002年各年中相同季度的数值加以平均消除了大部分随机性，因此这四个平均数仅仅代表了季节性。用代数式表示即为 （2） 其中中上面的横线表示季节平均。表3 产品销售额的季节性指数 年份各季度季节指数 第一季度第二季度第三季度第四季度 199176.40100.15 1992115.49111.3673.05104.03 1993111.29113.7476.54106.51 1994111.02110.8272.16103.67 1995115.011

9、08.4476.71102.39 1996114.78108.6175.70107.94 1997112.75107.1075.32106.15 1998111.34105.1676.03105.69 1999111.53109.5078.732000111.81110.7877.26104.16 2001110.89111.4577.1799.03 2002111.84111.7880.26 平均数112.72109.8876.28103.86 修正平均数111.95109.1375.76103.16 表3中的四个平均值相加的和为402.74，它不等于400。为了使各季节指数的平均数等于10

10、0，必须进行简单的调整。如果400被合计数402.74来除，结果是0.9932。以0.9932乘以各季节的平均数得到111.95，109.13，75.76，103.16等（见表中最后一行）。现在这四个季节指数的和为400，它们的含义就更加清楚了，例如第二季度的109.13就表示第二季度比全年平均数高出9.13，第三季度的75.76表示第三季度比全年低24.24。3长期趋势和循环变动 前面介绍的公式MA=TC表示了一组循环变动长期趋势数值。在多数情况下这样已能满足要求，但有时仍需要把循环变动和长期趋势分离开来。为了做到这一点，我们只需确定一种能最好的描述数据长期趋势的类型。例如长期趋势可以是线性

11、的、二次的、S曲线或其它。对于本例，如果将数据在图上画出来，可以看出线性的长期趋势是比较合适的：Tt = a + bt（.3） t = 1,2,348。用最小二乘法可求得模型的最佳拟合参数为：a = 2735.85, b = 38.96 因此趋势直线方程为 T_t=2735.85+38.96t 如图4所示。用此方程即可求得每个季度的趋势值。如第20季度（2000年的第四季度）趋势值为 T_20=a+bt=3515.05 由于MA=TC，因此 MA/T=（Ttimes C）/T=C （4） 应用上式即可求得循环变动值C。如第45季度的循环变动值C_45等于表3中的移动平均数除以T_45，即 如同

12、季节指数，循环指数也采取百分比率。其值大于100的表明该季度经济活动水平高于所有季度的平均值，而小于100的循环指数所表明的情况则刚好相反。循环因子比较复杂，且其变动周期较长，因而在短期预测中可以忽略不计，或将其归入到趋势变化之中（称为趋势循环因子）。人们更关心的是趋势和季节的识别。至此我们完成了对原始数据Xt的分解工作，其步骤总结如下：1） 用MA=TC分析长期趋势和循环变动；2） 用分析季节性和随机性；3） 用分析季节性；4） 用趋势外推法中介绍的方法来分析长期趋势；5） 用MA / T = C分析循环变动。总之，分解法提供了分析时间序列各种因素的手段，它使用简单，只需用加法、乘法和除法等

13、简单代数运算即可，而且分解法非常直观，能给企业提供一定时期内的大量信息。根据时间序列分解法进行预测 用分解法确定了季节指数、趋势值和循环指数之后，就可以根据上面总结的步骤进行预测了。我们对2003年第一季度（第49季度）进行预测。数据的基本关系式为 X=TCI 由于随机性无法直接进行预测，进行预测的关系式为：S 于是，计算出第49季度的T49,C49,S49值即可求得第49季度的预测值。表3中已得到第一季度的季节指数为111.95，由趋势方程求得 最后循环指数通常要根据判断估算出来，或者用某种方法预测得到。这里我们假定通过判断为：C49 = 98，于是 同样可以对第50、51季度进行预测。时间

14、序列分解法的进一步说明 1. 居中移动平均数 为了求得移动平均数MA，上面我们是将相邻的4个原始数据相加取平均得到一个数，这样在表4.5的第三列中就少了三个数据。于是产生了这样一个问题：最初的四个数据被平均时，它们的平均数应该置于何处？严格讲应该放在第二季度和第三季度的中间（14）/22.5，第2.5个季度）。其余数据取平均时也有类似的问题。但实际数据是表示各个季度而不是半个季度的，这里我们只好将平均数放在靠后半个季度的地方。假如对平均数再取平均的话就不会产生这样的问题了，因为如第一季度至第四季度的平均数2741.34是指第2.5季度，而第二季度至第五季度的平均数是指第 3.5季度，则它们的平

15、均数就是指第3个季度（2.53.5）/2=3）。称如此的平均数为居中移动平均数，于是居中移动平均数比原始数据少四个（首尾各两个）。现在，实际值除以居中移动平均值所得的比率（还是SI）也可以用来计算季度指数，具体的与上面所述完全一样。这样求得的四个季度的季节指数分别为112.20，109.44，75.37，103.17，其和为400.18，非常接近于400，这是因为移动平均数居中的缘故。2分解法的改进 在上面所叙述的分解法基础上，我们也可作一些改进，如：1） 修正原始数据中工作日或营业日的差额。由于各个月度（或季度）的工作日是不尽相同的，这就会影响到销售额或别的所要预测的变量。因此首先必须对数据

16、进行校正。如对月度数据的校正可通过原始数据乘以30对工作日的比率来进行，即将各月度的原始数据折算到工作日均为30天的统一情况。2） 利用统计方法来淘汰极值（即修改或舍去超出标准差的三倍范围的数值），在分解法实施之前先对数据进行预处理。3） 按上一节求得的季节性指数还可进一步改进，并进行动态的调整，因为实际上季节指数并不一定是一成不变的，它本身亦是一个变化的时间序列。还应注意到用分解法进行预测时，循环因素的确定是最为困难的。如有什么秘诀的话，那就是应具备足够数量的历史数据，以使管理人员了解循环模式是从哪里开始重复的，必要时可用图表方法来帮助确定。由于循环模式可能会发生变化，按照管理人员的判断对循

17、环模式作一些调整无疑是必要的。在前面的两个子节中，我们是以周期为4的季度数据的一个例题来说明分解法的分解步骤和预测程序。对周期为12月度数据、周期为7的日常数据等其它情况，运用分解法的程序完全类似，在此不再举例讨论。分解法能帮助解释历史数据为什么变化，能使管理人员分别预计各局部模式的变化。这些局部模式不仅能用以预测，而且也可用于管理之中，再加上它容易被管理人员所理解，因此分解法在直观上吸引了许多管理人员的注意，从而被大量的用于实际问题的预测。经过成千上万个时间序列的反复检验，分解法被证明其效率和准确性都是较高的。当然这种证明是经验的而非理论的，这也是它的主要缺点。它不能用统计的方法来检验，也不

18、能建立置信区间。实际上，分解法仅适用于那些季节性较强的中期预测、短期预测，当预测目标受外界干扰较大时，其预测能力会明显减弱。时间序列分解法案例分析 案例一:运用时间序列分解法进行销售预测时间序列是由一个包括了4个部分的模型组合而成，即T、S、C、R。假定影响时间序列的这四个因素彼此相互作用、非独立，那么我们采用相乘的模型4,P460，即Y，=TR。该方法包括以下四个步骤：1.用4季度移动平均法确定季节性指数4,P463。其基本原理是用移动平均法来度量趋势和周期性组合（TC）。这种做法可以消除季节性和随机变动的影响，即S和R。做法如下：（1）计算时间序列中的4季度移动平均值（TC_1），例：（5

19、00+350+250+400）4=375。（2）对4季度移动平均值再求其移动平均值的中点值（TC_2）2，例：（375+362.5/2=368.75）。计算真实销售额（3）计算真实销售额（Y_1）与移动平均值（TC_2）的比率，这个比率实际上表示的季节性和随机变动综合作用的部分2,即SR=Y/TC。（4）把计算处的比率值按季度排列，例如 第一季度的SR值分别是：1.263、1.367、1.467、1.222、1.348；第二季度的SR值分别是：1.037、0.762、0.812、0.853、0.879；第三季度的SR值分别是：0.678、0.640、0.522、0.588、0.700；第四季度

20、的SR值分别是：1.003、1.067、1.206、1.275、1.116。然后后按季度分别计算平均比率以便剔除随机变动（R）的影响，而该平均比率称为季节性因子（S_1）。例：对于第一季度的计算：（1.263+1.367+1.467+1.222+1.348）5=1.333 4。依次类推计算第二、三、四季度的s1分别得：1.33 4、0.908 6、0.625 6、1.153 4。（5）对季节性因子（S1）进行调整，调整后的季节性因子（用S2表示）。例：第一季度的S2计算：1.333 40-0.005 25=1.328 15，其中0.00525=（1.333 4+0.908 6+0.625 6+

21、1.153 4-4）/5。同样，第二、三、四季度的S2分别是：0.903 35、0.620 35、1.148 15。上述（1）、（2）、（3）、（4）、（5）的计算结果如下表。2.从原始时间序列中剔除季节性变化影响，即进行非季节性处理2。从附图中观测销售额与时期是否有线性变化趋势。通过观测发现：19941995年销售额有下降的趋势，1996-1999年销售额有上升的趋势。但是，前两年下降的趋势不明显，后四年一直是上升趋势。因此，可以近似看作销售额（Y）与时间（t）有线性关系，并依据非季节性数据作线性回归方程：Y2 = a + bt其中将表3中的数据代入上述公式得：Y=149.673+20.88

22、2t（3） 3.将t=25、26、27、28代人第二步求得的回归方程（3），并乘上相应的季节性因子S2，则得出2000年每一季度的销售额预测值。即：Y25 = 892.149、Y26 = 625.665、Y_27=442.612、Y28 = 843.1664，在附图上作出上述回归方程（3）的图形。采用时间序列分解法时，必须观察时间t与非季节性数据Y_2的趋势关系。若t与Y2的趋势变化是近似于线性的，则可用线性回归方程预测;若Y2对t来说是呈几何级数增长，则宜用指数曲线回归方程;若是其他变化，必须采用相应的其他回归方程。采用线性回归方程来预测，但是由于前两年的数值是下降的，后四年的数值是上升的，因此，使得该回归方程的拟合度不是很好。但是，这种影响程度会随着t的增加而减弱，最终不会影响预测的精确度，因为总趋势是上升的。，采用时间序列分解法预测时，由于前两年的数值有下降的趋势，而拟合的是上升趋势的线性回归方程，因此，用该方法预测短期（1-2年）的销售额不准确，通过附图可以看出。