物流需求预测

1、第1 1章 物流需求预测 内容提要 物流需求预测的含义，预测是指从已知事件测定未知事件。物流需求预测的含义，预测是指从已知事件测定未知事件。 时间序列预测法（算术平均法、加权平均法、移动平均法、时间序列预测法（算术平均法、加权平均法、移动平均法、指数平滑法）。指数平滑法）。 相关因素预测法（一元线性回归法、多元线性回归法等）。相关因素预测法（一元线性回归法、多元线性回归法等）。 预测方法在预测方法在Excel 2003Excel 2003中的具体实现。中的具体实现。CHAPTERq 预测理论作为通用的方法论，既可以应用于研究自预测理论作为通用的方法论，既可以应用于研究自然现象，又可以应用于研究

2、社会现象。将预测理论和然现象，又可以应用于研究社会现象。将预测理论和方法与具体领域的实际相结合，就形成了预测学科的方法与具体领域的实际相结合，就形成了预测学科的各个分支，如社会预测、人口预测、经济预测、政治各个分支，如社会预测、人口预测、经济预测、政治预测、科技预测、军事预测、气象预测等。预测、科技预测、军事预测、气象预测等。 q 本章内容仅讨论定量分析预测方法。本章内容仅讨论定量分析预测方法。q 在定量分析预测中应用比较广泛的有时间序列预测在定量分析预测中应用比较广泛的有时间序列预测法（包括算术平均法、加权平均法、移动平均法、指法（包括算术平均法、加权平均法、移动平均法、指数平滑法、最小二乘

3、法等）和相关因素预测法（包括数平滑法、最小二乘法等）和相关因素预测法（包括一元线性回归法、多元线性回归法等）。一元线性回归法、多元线性回归法等）。1.1 时间序列预测法 序列又称时间数列，是指观察或记录到的一组按时间顺序列又称时间数列，是指观察或记录到的一组按时间顺序排列的数据。例如，某段时间内，某类产品产量的统序排列的数据。例如，某段时间内，某类产品产量的统计数据；某企业产品销售量、利润、成本的历史统计数计数据；某企业产品销售量、利润、成本的历史统计数据；某地区人均收入的历史统计数据等。据；某地区人均收入的历史统计数据等。 时间序列预测方法，假设预测对象的变化仅与时间有关。时间序列预测方法，

4、假设预测对象的变化仅与时间有关。根据它的变化特征，以惯性原理推测其未来状态。事实根据它的变化特征，以惯性原理推测其未来状态。事实上，预测对象与外部因素有着密切而复杂的联系。时间上，预测对象与外部因素有着密切而复杂的联系。时间序列中的每一个数据都反映了当时许多因素综合作用的序列中的每一个数据都反映了当时许多因素综合作用的结果。整个时间序列则反映了外部因素综合作用下预测结果。整个时间序列则反映了外部因素综合作用下预测对象的变化过程。因此，预测对象仅与时间有关的假设，对象的变化过程。因此，预测对象仅与时间有关的假设，是对外部因素复杂作用的简化，从而使预测的研究更为是对外部因素复杂作用的简化，从而使预

5、测的研究更为直接和简便。直接和简便。 时间序列预测方法分为确定性方法和随机性方时间序列预测方法分为确定性方法和随机性方法两类。随机性时间序列预测方法是研究由随机过法两类。随机性时间序列预测方法是研究由随机过程产生的时间序列的预测问题。本教材仅介绍确定程产生的时间序列的预测问题。本教材仅介绍确定性时间序列预测方法。现实中的时间序列的变化受性时间序列预测方法。现实中的时间序列的变化受许多因素的影响，有些起着长期的、决定性的作用，许多因素的影响，有些起着长期的、决定性的作用，使时间序列的变化呈现出某种趋势和一定的规律性；使时间序列的变化呈现出某种趋势和一定的规律性；有些则起着短期的、非决定性的作用，

6、使时间序列有些则起着短期的、非决定性的作用，使时间序列的变化呈现出某种不规则性。的变化呈现出某种不规则性。时间序列的变动情况时间序列的变动情况大体可分解为以下四种。大体可分解为以下四种。（1 1）趋势变化。趋势变化指现象随时间变化朝着一定方）趋势变化。趋势变化指现象随时间变化朝着一定方向呈现出持续稳定地上升、下降或平稳的趋势。向呈现出持续稳定地上升、下降或平稳的趋势。（2 2）周期变化。周期变化（季节变化）指现象受季节性）周期变化。周期变化（季节变化）指现象受季节性影响，按某一固定周期呈现波动变化。影响，按某一固定周期呈现波动变化。 （3 3）循环变动。循环变动指现象按不固定的周期呈现波）循环

7、变动。循环变动指现象按不固定的周期呈现波动变化。动变化。（4 4）随机变动。随机变动指现象受偶然因素的影响而呈）随机变动。随机变动指现象受偶然因素的影响而呈现出不规则波动现出不规则波动 1.2 平均数预测法 平均数预测法是最简单的定量预测平均数预测法是最简单的定量预测方法，它只在时间序列主要表现为随机方法，它只在时间序列主要表现为随机变动时采用。下面我们具体介绍算术平变动时采用。下面我们具体介绍算术平均数和加权平均数两种预测方法。均数和加权平均数两种预测方法。 1.2.1 算术平均数预测法时间（t）123n观察值（yt） y1y2y3 yn 设有时间序列如表设有时间序列如表1.1所示，对应于时

8、间所示，对应于时间t t 1, 2, , 1, 2, , n n ；时间序列的数据为：时间序列的数据为：y y1 1, , y y2 2, , y y3 3, , , , y yn n。 表表1.1 时间序列数据时间序列数据 设算术平均数为，其计算公式为设算术平均数为，其计算公式为 （1-11-1）12111nnnttttyyyyyynnn式中：式中： 算术平均数；算术平均数； y yt t 第第t t周期的实际值；周期的实际值； t t 时间下标变量，表示周期序号；时间下标变量，表示周期序号； n n 时间序列的周期个数，也即数据个数。时间序列的周期个数，也即数据个数。 连加号，连加号， 中

9、的中的t t 1 1表示首项的时间表示首项的时间变量取值为变量取值为1 1，n n表示末项的时间下标变量取值为表示末项的时间下标变量取值为n.n.y1nt1.2.1 算术平均数预测法例例1.1 某公司的产品在某城市最近某公司的产品在某城市最近6个月的销售量如表个月的销售量如表1.2所示。所示。 表表1.2 最近最近6个月的销售量个月的销售量月月 份份 123456销售量销售量/台台 105010801030107010501060 观察实际数据序列，其变动特征主要为随机变动。观察实际数据序列，其变动特征主要为随机变动。因此，可采用算术平均数方法预测下个月的销售量，因此，可采用算术平均数方法预测

10、下个月的销售量，即即 （台）（台） 算术平均数算术平均数1 0571 057台反映了最近台反映了最近6 6个月产品销量的个月产品销量的平均水平。如果判断影响产品销量的外界因素无重大平均水平。如果判断影响产品销量的外界因素无重大变化，即可预测下个月的销售量为变化，即可预测下个月的销售量为1 0571 057台。台。1212661050 1080 1030 1070 1050 106010576nyyyyyyyn 算术平均数算术平均数10571057台反映了最近台反映了最近6 6个月产品销量的平均水平。个月产品销量的平均水平。如果判断影响产品销量的外界因素无重大变化，即可预测下个月如果判断影响产品

11、销量的外界因素无重大变化，即可预测下个月的销售量为的销售量为10571057台。台。 此外，还可用极差和标准偏差研究数据的变动情况。此外，还可用极差和标准偏差研究数据的变动情况。 设极差为设极差为R R，其计算公式，其计算公式 R R x xmaxmax x xmin min （1-21-2） 式中：式中：R R极差；极差； x xmaxmax所有观察值中的最大值；所有观察值中的最大值； x xminmin所有观察值中的最小值。所有观察值中的最小值。 极差极差R R表示数据的最大变动幅度。本例中表示数据的最大变动幅度。本例中，x xmax max 1 0801 080，x xmin min 1

12、0301030，因此，因此R R 5050（台），即销售量的最大变动幅度为（台），即销售量的最大变动幅度为5050。 设标准偏差为设标准偏差为S S，其计算公式为，其计算公式为 （1-31-3）211()nttSyyn 式中式中S 表示数据的离散程度。由公式可知，表示数据的离散程度。由公式可知，S 越小，数据的越小，数据的离散程度越小。离散程度越小。 本例标准偏差计算如下：本例标准偏差计算如下： （台）（台） S 16台反映了数据的离散程度。台反映了数据的离散程度。 算术平均数、极差算术平均数、极差R、标准偏差、标准偏差S 是统计学中描述一组数据变动是统计学中描述一组数据变动状况的三个数字特征

13、状况的三个数字特征。62111534()1666ttSyy1.2.2 加权平均数预测法 如果认为近期的数据反映的预测信息比较早的数据更重要，也就是说各个如果认为近期的数据反映的预测信息比较早的数据更重要，也就是说各个时间数据的重要程度不相等，则可以采用加权的方法对数据进行处理。时间数据的重要程度不相等，则可以采用加权的方法对数据进行处理。 设设y1 , y2 , y3 , , yn为各个时间的观察值，对每一个数据估计相应的权数，分为各个时间的观察值，对每一个数据估计相应的权数，分别表示为别表示为W1, W2 , W3 , , Wn 。则加权平均数。则加权平均数 的计算公式为的计算公式为 （1-

14、4）式中：式中： 加权平均数；加权平均数； yt第第t周期的观察值；周期的观察值； Wt第第t周期观察值的权数，通常设权数之和等于周期观察值的权数，通常设权数之和等于1，即，即 。y11221121nttnntnnttW yW yWyWyWWWW y yy11nttW 例1.2 沿用例1.1的数据，并设各时期销售量的权数依次为0.1，0.1，0.15，0.15，0.25，0.25。其加权平均数计算如下： （台） 同理，如果判断影响预测对象变化的外部因素无重大变化，即可将1 056台作为下个月产品销量的预测值。y 611611nttttttnttttW yW yWW105625. 025. 01

15、5. 015. 01 . 01 . 0106025. 0108011. 010501 . 0 1.3 移动平均预测法 移动平均预测法是一种简单平滑预测技术，它是在算术平均数的基础上发展起来的一种预测方法。它的基本思想是：根据时间序列资料，逐项推移，依次计算包含一定项数的时间序列平均值，以反映长期趋势的方法。因此，当时间序列的数值由于受周期变动和随机波动的影响起伏较大，不易显示出事件的发展趋势时，使用移动平均预测法可以消除这些因素的影响，显示出事件的发展方向与趋势（即趋势线），然后依趋势线分析预测序列的长期趋势。1.3.1 一次移动平均法一次移动平均数的计算公式为 （1-5）式中：t 周期序号；

16、 第t周期的平均数； 第t周期的一次移动平均数； 第t周期的实际值； N 计算移动平均数所选定的数据个数。以第t周期的一次移动平均数作为第t1周期的预测值，即 。 11(1)ttt NtyyyMN tM(1)tMty(1)1ttyM例1.3 已知某种产品最近15个月的销售量如表1.3所示。 表1.3 最近15个月的销售量 (1)tM(1)tM时间序号（时间序号（t） 123456789101112131415销售量（销售量（yt ） 1015 8 20 1016 18 20 22 24 20 26 27 29 29 (N 3) 11.0 14.3 12.7 15.3 14.7 18.0 20.

17、0 22.0 22.0 23.3 24.4 27.3 28.3 (N 5) 12.6 13.8 14.4 16.8 17.2 20.0 20.8 22.4 23.8 25.2 26.2 分别取N 3，N 5，计算一次移动平均数，结果见表1.3。 解解 取N 3，N 5，按式（1-5）计算，并把计算结果填入表中。表中一次移动平均数的计算过程举例如下（N 3）： 323 3 1321(1)3432(1)48151011.03332081514.333yyyyyyMyyyM 其余类推。 由上述计算过程可看出， 比 增添了一个新数据y4，同时去掉了一个旧数据y1，这就是移动平均的意思。 移动平均数也构

18、成了时间序列。为了研究移动平均法的特点，将实际数据序列和N 3，N 5的移动平均序列均绘于图1.1中。(1)4M(1)3M分析图1.1，可看出：（1）移动平均法可以削弱随机变动的影响，具有平滑数据的作用。移动平均数序列比实际数据序列平滑，能在一定程度上描述时间序列的变化趋势。（2）合理地选择模型参数N值，是用好移动平均法的关键。N越大，平滑作用越强，对新数据的反应越不灵敏；N越小，则效果相反。（3）在实际序列的线性增长部分，移动平均数的变化总是落后于实际数据的变化，存在着滞后偏差。当N越大时，滞后偏差也越大。 应注意到：当N 1时，移动平均数序列即为实际序列；当N等于全部数据的个数n时，移动平

19、均数即为算术平均数。通常根据实际序列的特征和经验选择模型参数N。N的取值范围可在320之间。 如果实际时间序列没有明显的周期变动和倾向变动，即可用最近时间的一次移动平均数作为下一周期的预测值。 若不考虑例1.3中线性增长趋势的影响，则可求得上一个月的销售量的预测值 。 取N 3，得出 （万台） 显然，从实际数据序列的变动情况来看， 的数值偏低。这是因为滞后偏差引起的。改进的方法是采用二次移动平均等方法。16 y1514131629 29 2728.333yyyy16 y1.3.2 二次移动平均法 如果时间序列具有明显的线性变化趋势，则不宜采用一次移动平均法进行预测。原因是滞后偏差将使预测值偏低

20、，不能合理地进行趋势外推；因此，需要进行修正。修正的方法是在一次移动平均的基础上再进行二次移动平均，利用移动平均滞后偏差的规律找出曲线的发展方向和发展趋势，然后才建立直线趋势的预测模型。故也称为趋势移动平均法。 二次移动平均数是在一次移动平均数的基础上经过计算得到的。其计算公式为 （1-6）式中： 第t周期的一次移动平均数； 第t周期的二次移动平均数； N计算移动平均数所选定的数据个数。(1)(1)(1)(2)11ttt NtMMMMN(1)tM(2)tM例1.4 沿用例1.3的数据，若取N 3，其二次移动平均数的计算结果如表1.4所示。表1.4 二次移动平均数的计算结果 实际序列和一次、二次

21、移动平均数序列的图形如图1.2所示。试建立预测模型，并求第16、17周期的预测值。 时间序号（时间序号（t） 123456789101112131415销售量（销售量（yt ）101582010 16 18 20 22 24 20 26 27 29 29 (N 3) 11.0 14.312.7 15.3 14.7 18.0 20.0 22.0 22.0 23.3 24.4 27.3 28.3 (N 3) 12.7 14.1 14.2 16.0 17.6 20.0 21.3 22.4 23.3 25.0 26.6 (1)tM(2)tM 观察图1.2，一次平均数序列总是落后于实际数据序列，出现了滞

22、后偏差；二次移动平均数序列也与一次平均数序列形成了滞后偏差。二次移动平均法正是利用这种滞后偏差的演变规律建立线性预测模型的。线性预测模型为 （1-7）式中：t目前的周期序号； T由目前周期t到预测周期的周期间隔个数，即预测超前周期数； 第t T周期的预测值； at线性模型的截距； bt线性模型的斜率，即单位周期的变化量。ttt TyabTt Ty图1.2 例1.4图 其中at、bt的计算公式为 （1-8） （1-9） (1)(2)2tttaMM(1)(2)2()1tttbMMN 下面介绍对例1.4的线性预测模型的求解过程。 已知目前周期序号t 15，将第15周期的一次、二次、移动平均数代入式（

23、1-8）和式（1-9）得(1)(2)151515(1)(2)1515152228.326.630.022()(28.326.6)1.712aMMbMMN得到线性预测模型为1530.01.7ttTyabTT 试求下一个月的销售量预测值。下一个月的周期序号t 16，即周期间隔数T16151，故 （万台）求第17周期的销售量预测值。此时，周期间隔数T17152，故 （万台）15 1 y161515131.7yab 15 2 y171515233.4yab于是利用线性预测模型，求得第16周期和第17周期的产品销售量预测值分别为31.7万台和33.4万台。 预测模型，即式（1-7）是一个直线方程，可以预

24、测最近几个周期的销售量（模型的推导过程从略）。如果实际数据序列具有曲线变化趋势时，就需要利用三次移动平均法。但在这种场合下，下一节介绍的指数平滑预测法更为常用。 1.4 指数平滑预测法 短期预测中最有效的方法可能就是指数平滑预测法，也可以简称为指数平滑法。该方法很简单，只需要得到很小的数据量就可以连续使用。 1.4.1 指数平滑法的基本理论 指数平滑法是移动平均法中的一种，其特点在于给过去的观测值不一样的权重，即较近期观测值的权数比较远期观测值的权数要大。根据平滑次数不同，指数平滑法分为一次指数平滑法、二次指数平滑法和三次指数平滑法等。但它们的基本思想都是：预测值是以前观测值的加权和，且对不同

25、的数据给予不同的权数，新数据给予较大的权数，旧数据给予较小的权数。1.4.2 一次指数平滑法 设时间序列为 ,则一次指数平滑公式为 （1-10）式中： 第t周期的一次指数平滑值； 加权系数（也称为平滑系数），01。12,tyyy(1)(1)1(1)tttSyS(1)tS用上述平滑值进行预测，就是一次指数平滑法。其预测模型为 （1-12）即以第t周期的一次指数平滑值作为第t1期的预测值。1ty(1)(1)1(1)(1)tttttSyyyS例1.5 用一次指数平滑值预测例1.3中下个月的销售量 。 为了分析加权系数的不同取值的特点，分别取 0.1, 0.3, 0.5计算一次指数平滑值，并设初始值为

26、最早的三个数据的平均值，即： 。以 0.5的一次指数平滑值计算为例，有16 y123(1)011.03yyyS 0.5100.511.010.5 0.5150.510.512.8(1)(1)110(1)SyS(1)(1)221(1)SyS依次类推，求得 0.5、 0.1和 0.3时的一次指数平滑值数列，计算结果如表1.5所示。 表1.5 一次指数平滑值的计算表 万台 指数平滑法的基本特点如下。（1）指数平滑法对实际序列具有平滑作用，权系数（平滑系数） 越小，平滑作用越强，但对实际数据的变动反应较迟缓。（2）在实际序列的线性变动部分，指数平滑值序列出现一定的滞后偏差的程度随着权系数（平滑系数）

27、的增大而减少。指数平滑法的主要优点如下 （1）对不同时间的数据的非等权处理较符合实际情况。 （2）实用中仅需选择一个模型参数 即可进行预测，简便易行。 （3）具有适应性，也就是说预测模型能自动识别数据模式的变化而加以调整。 指数平滑法的缺点 一是对数据的转折点缺乏鉴别能力，但这一点可通过调查预测法或专家预测法加以弥补。 二是长期预测的效果较差，故多用于短期预测。1.4.3 二次指数平滑法 如果实际数据具有明显的趋势变动时，使用第t周期一次指数平滑就能直接预测第t1周期之值。但当时间序列的变动出现直线趋势时，用一次指数平滑法来进行预测仍将存在着明显的滞后偏差。因此，也需要进行修正。修正的方法也是

28、在一次指数平滑的基础上再进行二次指数平滑，利用滞后偏差的规律找出曲线的发展方向和发展趋势，然后建立直线趋势预测模型，故称为二次指数平滑法。 设一次指数平滑为 ，则二次指数平滑 的计算公式为 （1-13）式中： 第t周期的二次指数平滑值； 第t周期的一次指数平滑值； 第t1周期的二次指数平滑值； 加权系数（也称为平滑系数）。(1)tS(2)tS(2)(1)(2)1(1)tttSSS(2)tS(1)tS(2)1tS 若时间序列 从某时期开始具有直线趋势，且认为未来时期亦按此直线趋势变化，则与趋势移动平均类似，可用如下的直线趋势模型来进行预测。 （1-14） 式中：t当前周期数； T由当前周期数t到

29、预测期的时期数； 第tT期的预测值；12,tyyy,1, 2,tt Ttyab TTt Ty 截距， 斜率，其计算公式为 （1-15） （1-16） 模型公式（1-14）与二次移动平均法建立的线性模型的形式是一样的。 tatb(1)(2)2tttaSS(1)(2)1tttbSS1.4.4 三次指数平滑法 若时间序列的变动呈现出二次曲线趋势，则需要采用三次指数平滑法进行预测。三次指数平滑是在二次指数平滑的基础上再进行一次平滑，其计算公式为 （1-17）(3)(2)(3)1(1)tttSSS三次指数平滑法的预测模型为 （1-18） 式中： （1-19） （1-20） （1-21） 2tttt Ty

30、ab Tc T(1)(2)(3)33ttttaSSS(1)(2)(3)2(65 )2(54 )(43 )2(1)ttttbSSS2(1)(2)(3)222(1)ttttcSSS例1.6 我国某种耐用消费品1996年至2006年的销售量如表1.6所示。试预测2007年、2008年的销售量。表1.6 三次指数平滑值的计算表 万台 时间序号（时间序号（t） 01234567891011年份年份 1996 1997199819992000200120022003200420052006销售量（销售量（yt ） 225.2 249.9 263.2 293.6 318.9 356.7 363.3 424.

31、2 466.5 582.0 750.0 ( 0.3) 246.1 237.8 242.8 248.9 262.3 279.3 302.5 320.8 351.9 386.3 445.0 536.5 246.1 244.2 243.7 245.3 250.4 259.1 272.1 286.7 306.3 330.3 364.7 416.2 244.5 244.4 244.3 244.5 246.3 250.1 256.7 265.7 277.9 293.6 314.9 345.3 (1)tS(2)tS(3)tS解解 通过初步分析，实际数据序列呈非线性递增趋势，故采用三次指数平滑预测方法。解题步

32、骤如下。（1）确定指数平滑的初始值和权系数（平滑系数）。设一次、二次指数平滑的初始值为最早三个数据的平均值，即 ， 取 。 实际数据序列的倾向性变动较明显，权系数（平滑系数） 不宜取太小，故取 0.3。 123(1)(2)00225.2249.9263.2246.133yyySS(3)0244.5S（2）根据指数平滑值计算公式依次计算一次、二次、三次指数平滑值 结果见表1.6。 （3）计算非线性预测模型的系数at, bt, ct。目前周期数t 11，将表1.6中的有关数据代入式（1-19）、式（1-20）、式（1-21）后分别得 (1)(2)(3)11111111(1)(2)(3)111111

33、11222(1)(2)1111112333 536.5 3 416.2 345.3 706.2(6 5 )2(5 4 )(4 3 )2(1)0.3(6 5 0.3)536.5 2(5 4 0.3)416.2 (4 3 0.3)345.3 98.42(1 0.3)(22(1)aSSSbSSScSSS (3)11) 4.4（4）建立非线性预测模型。 将各系数代入式（1-18）得22706.2 98.44.4( 11 )tttt Tya bT cTTTt （5）预测2007年和2008年的产品销售量。2007年，其预测超前周期为T 1；2005年，其预测超前周期为T 2。代入模型，得 706.298

34、.44.412 809（万台） 706.298.424.422 920（万台） 于是得到2007年的产品销售量的预测值为809万台，2008年的产品销售量的预测值为920万台。预测人员可以根据市场需求因素的变动情况，对上述预测结果进行评价和修正。 2004 y11 1 y2005 y11 2 y1.4.5 加权系数的选择 在指数平滑法中，预测成功的关键是 的选择。 的大小规定了在新预测值中新数据和原预测值所占的比例。 值愈大，新数据所占的比重就愈大，原预测值所占比重就愈小，反之亦然。 若把一次指数平滑法的预测公式改写为 则从上式可以看出，新预测值是根据预测误差对原预测值进行修正后得到的。 的大

35、小表明了修正的幅度。 值愈大，修正的幅度愈大； 值愈小，修正的幅度愈小。因此， 值既代表了预测模型对时间序列数据变化的反应速度，又体现了预测模型修正误差的能力。 1ty(ttyy)ty 在实际应用中， 值是根据时间序列的变化特性来选取的。若时间序列的波动不大，比较平稳，则 应取小一些，如0.10.3；若时间序列具有迅速且明显的变动倾向，则 应取大一些，如0.60.9。实质上， 是一个经验数据，通过多个 值进行试算比较而定，哪个 值引起的预测误差小，就采用哪个。 例1.7 根据表1.7给出的1月份至11月份餐刀需求量的观察值，分别取 0.1，0.5，0.9，预测12月份的餐刀需求量，并对不同的

36、值进行误差比较。表1.7 用指数平滑法预测12月份的餐刀需求量 月份 时期 需求量的 观察值 指数平滑值 0.1 0.5 0.9 1月12 000 2月21 3502 0002 0002 0003月31 9501 9351 6751 4154月41 9751 9371 8131 8975月53 1001 9401 8941 9676月61 7502 0562 4972 9877月71 5502 0262 1231 8748月81 3001 9781 8371 5829月92 2001 9101 5681 32810月102 7751 9391 8842 11311月112 3502 0232

37、3302 70912月12 2 0562 3402 386 表1.7给出了 值为0.1, 0.5和0.9时所计算出来的预测值。表中最后三栏的数据可用指数平滑法计算预测值的通式 或 来进行计算。这时应注意，第一个时期没有前期预测值可以利用，最好能利用观察值。因此，1 935这个数，在 0.1这栏是用2 000（前期预测值）再加上0.1（1 3502 000）来得到的。然后把这个数用做第三时期的预测值。 值对前期观察值进行平滑修正所产生的效果在表1.7中可以看出，较大的 值（0.1）对预测值的平滑作用很小，而较小的 值则有着相当大的平滑作用。 1(1)tttSyS1()ttttSSyS 表1.8表

38、明，在此例中较小的 值比较大的 值能给出更好的预测值。 表1.8 指数平滑法预测误差比较 误差 绝对误差 误差平方 0.1 总数 551 4 771 3 431 255 平均数 55 477 343 126 0.5 总数674 5 688 4 338 332 平均数 67 569 433 833 0.9 总数423 6 127 5 034 081 平均数 42 613 503 408 1.5 回归分析预测方法1.5.1 回归分析模型概述回归分析模型概述 回归分析预测法是从各种经济现象之间的相互关系出发，通过对与预测对象有联系的现象的变动趋势的分析，推算预测对象未来状态数量表现的一种预测法。所谓

39、回归分析就是研究某一个随机变量（因变量）与其他一个或几个变量（自变量）之间的数量变动关系，由回归分析求出的关系式通常称为回归模型。 回归模型的一般分为如下几类。 （1）根据自变量个数的多少，回归模型可以分为一元回归模型和多元回归模型。 （2）根据回归模型是否线性，回归模型可以分为线性回归模型和非线性回归模型。所谓线性回归模型就是指因变量与自变量之间是直线关系。（3）根据回归模型是否带虚拟变量，回归模型可以分为普通回归模型和虚拟变量回归模型。普通回归模型的自变量都是数量变量，而虚拟变量回归模型的自变量既有数量变量也有品质变量。 本节先讨论一元线性回归问题，然后再讨论多元线性回归问题。 1.5.2

40、 一元线性回归模型 一元线性回归预测是回归预测的基础。若预测对象只受一个主要因素影响，并且它们之间存在着明显的线性相关关系时，通常采用一元线性回归预测法。1预测模型 设变量x与变量y之间有相关关系，且当x确定之后，y有某种不确定性，如果在散点图上可以看出x与y之间有线性相关关系，其相关方程为 （1-22） yabx采用最小二乘法得到a, b的计算公式为 （1-23） （1-24）式中，xi和yi（i 1, 2, n）均是已有的历史数据。 其中：a称为截距，b称为回归直线的斜率，也称为回归系数；是变量y的估计值。 22()iiiiiinx yxybnxxaybx y 求直线回归方程 ，实际上是用

41、回归直线拟合散点图中的各观测点。常用的方法是最小二乘法，也就是使该直线与各点的垂直距离最小，即求使观察值y与回归直线 之差的平方和 达到最小时的a和b的问题。 yabx y2()yy2显著性检验 判定一个线性回归方程的拟合程度的优劣称为模型的显著性检验，即判断所建立的一元线性回归模型是否符合实际，所选的变量之间是否具有显著的线性相关关系。这就需要对建立的回归模型进行显著性检验，通常用的检验法是相关系数检验法。相关系数是一元回归模型中用来衡量两个变量之间相关程度的一个指标，其计算公式为 （1-25） 相关系数r是一个重要的判定指标。从式（1-25）中可以看出，相关系数等于回归平方和在总平方和中所

42、占的比率，即回归方程所能解释的因变量变异性的百分比。如果r 0.775，说明变量y的变异性中有77.5%是由自变量x引起的；如果r 1，表示所有的观测值全部落在回归直线上；如果r 0，则表示自变量与因变量无线性关系。 下面通过例子说明一元线性回归分析预测法的使用。22()()()()iiiixxyyrxxyy2222()()iiiiiiiinx yxynxxnyy 例1.8 根据经验，运输企业的业务收入同广告费支出之间具有相关关系。某运输企业1997年至2006年的业务收入和广告费支出的资料如表1.9所示。 表1.9 某企业业务收入与广告费支出表 年份年份 广告费支出（广告费支出（xi）/万元

43、万元 业务收入（业务收入（yi）/百万元百万元 xi yi xi2 yi2 1997年年4728 16 49 1998年年 71284 49 144 1999年年 917153 81 289 2000年年 1220240 144 400 2001年年 1423 322 196 529 2002年年 1726 442 289 676 2003年年 2029 580 400 841 2004年年 2232 704 484 1 024 2005年年 2535 875 625 1 225 2006年年 2740 1 080 729 1 600 合计合计 157241 4 508 3 013 6 777

44、 预测该企业2007年的广告费支出为35万元，要求在95%的概率下预测该年的企业业务收入。解解 （1）相关分析。 在坐标系上将企业的广告费支出和业务收入的数据标出，形成散点图，可以发现图形呈直线趋势，从而判定二者呈一元回归性。（2）建立相关方程 ，关键是求a、b的值。根据表1.9中的有关资料，利用式（1-23）和式（1-24）求得 所以，所求的相关方程为 。3.361.321yx yabx22210 4508 157 2411.32110 3013 (157)()2411571.3213.361010iiiiiiiinx yxybnxxyxaybxbnn （3）检验。 计算相关系数：取显著性水

45、平 0.05，df n 2 8。查相关系数临界值表得r0.05(8) 0.632。因为r r，说明广告费支出与业务收入存在很强的正相关关系。222222()()104 5081572410.999 4103013(157)106 777(241)iiiiiiiinx yxyrnxxnyy 相关系数r2检验和F检验。相关系数r2检验和F检验均可用于检验回归方程线性关系的显著性，二者在检验原理上大体相同，均借助了方差分析，即 （1-26）式中： 总变差； 剩余变差； 回归变差。计算结果如表1.10所示。2()iyy2()iyy2()yy2()iyy2()iyy2 ()yy表1.10 检 验 结 果

46、时间 xi yi 1478.644 2.703 238.888 292.41 136.89 27 12 12.607 0.368 132.089 146.41 75.69 39 17 15.249 3.066 78 .340 50.41 44.89 412 2019.212 0.621 23. 893 16.81 13.69 514 2321.854 1.313 5.045 1.2 2.89 617 2625.817 0.033 2.948 13.61 1.69 7202929.78 0.608 32.262 24.01 18.49 8223232.422 0.178 69.256 62.41

47、 39.69 9253536.385 1.918 150.921 118.81 86.49 10274039.027 0.947 222.815 252.81 127.69 合计 157241240.997 11.755 956.457 968.92 546.10 y2()yy2()yy2()yy2()xx 相关系数r2利用回归变差、点变差、总变差的比重说明回归直线的代表性，若这个比例越大，则说明x与y之间关系越密切，回归直线代表性越好。一般地，r2的取值在01之间。222()956.54798.7%968.9()yyryy 用F检验法将自变量作为一个整体来检验与因变量之间的线性关系是否显著，

48、其计算如下： 取显著性水平 0.05，df1 1，df2 n 2 8。查F分布表得r0.05(1.8) 5.32。 因为F F，说明广告费支出与业务收入线性关系显著。这与相关系数r2检验所得结论一致。22()956.547651.20411.7558() /(2)yyFyyn3预测分析 （1）点预测。2007年的广告费支出预计为35万元。将x035（万元）代入回归方程，得 3.36 1.321 3549.595（百万元） 即2007年的业务收入可望达49.595百万元。y （2）区间预测 计算估计标准误差，即因为 0.05，df n 2 8，查t分布表得 2()11.7551.21228iyy

49、ySn,22nt0.025,82.036t 当广告费支出达到x0 35（万元）时，业务收入的预测区间为 即若以95%的把握程度预测，当广告费支出达到35万元时，企业的业务收入在45.86453.326百万元之间。20,22221()1()1(3515.7)49.5952.3061.212149.5953.73110546.1ynixxytSnxx1.5.3 多元线性回归 多元线性回归分析预测法是对自变量和因变量的n组统计数据（x1i, x2 i, , x m i ; yi ）（i1,2,n），在明确因变量y与各个自变量间存在线性相关关系的基础上，给出适宜的回归方程，并据此做出关于因变量y的发展

50、变化趋势的预测。因此，多元线性回归分析预测法的关键是找到适宜的回归方程。 类似于一元线性回归分析，可以用线性方程y a b1x1 b2x2 bmxm （1-27） 来近似描述y与x1, x2, , xm之间的线性相关关系，其中 为根据所有自变量计算出来的估计值，a为常数项， 称为y对应于x1, x2, , xm的偏回归系数。偏回归系数是假设在其他所有自变量保持不变的情况下，某一个自变量的变化引起因变量变化的比重。它的参数也可以用最小二乘法进行估计。建立一个多元回归模型需要用到复杂的统计方法，现在可以使用计算机软件包来根据统计数据建立合适的多元回归方程，这样就会方便许多。自变量个数为2个的多元线

51、性回归方程称为二元回归方程，它是多元线性回归方程中的特例。多元（以二元为例）线性回归分析的步骤如下。 y12,mb bb,（1）建立线性方程y a b1x1 b2x2式中参数a、b1、b2仍使用最小二乘法推算，即 （1-28）将相关数据代入上述方程组，得到系数a、b1、b2。所以，二元线性回归方程为 a b1x1 b2x2 （1-29 ）11222111212122211222iiiiiiiiiiiiiiiynabxbxx yaxbxbx xx yaxbx xbxy （2）显著性检验。利用复相关系数检验回归方程整体显著性。 （1-30） 取一个特定的 ，并计算出dfnk1（k为自变量个数），查

52、相关系数临界值表得到R df 。如果RR，说明x1、x2与y的线性关系显著。22()1()iiyyRyy21122221iiiiiiiyaybx ybx yyny （3）预测分析。 点预测。将x1、x2代入公式 得到预测值 。 区间预测。计算估计标准误差，即 （1-31） 取一个特定的，df n3，查t分布表得到t/2 d f 。所以，预测区间为 。 1 122 yab xb x y2()3iyySn/2df ytS例1.9 假设商品销售额除了同广告费支出相关外，还同营业网点数有一定的相关关系。资料如表1.11所示。 表1.11 例1.9资料 时间时间 广告费支出（广告费支出（x1i）/万元万

53、元 营业网点数（营业网点数（x2i）/个个 商品销售额（商品销售额（yi）/百万元百万元 1997年年 4 1 7 1998年年 7 2 121999年年 9 5172000年年 12 8202001年年 14 10232002年年 17 12 26 2003年年 20 15292004年年 22 18322005年年 25 22352006年年 27 2840合计合计 157 121 241 如果2007年该企业的广告费支出为35万元，营业网点数为34个，要求以95%的概率预测企业2007年的商品销售额。（1）建立线性方程。从表1.11中可以看出，商品销售额与广告费支出、营业网点数两个因素均

54、存在相关关系。所以拟合得到二元线性回归方程y a b1x1 b2x2 式中参数a、b1、b2仍使用最小二乘法推算，即 有关数据的计算结果如表1.12所示。 11222111212122211222iiiiiiiiiiiiiiiynabxbxx yaxbxbx xx yaxbx xbx表1.12 例1.9计算结果 时间时间 x1i x2i yi x1iyi x2iyi x1i2 x2i2 x1ix2i 1997年年 41728716141998年年 72128424494141999年年 9517153858125452000年年 1282024016014464962001年年 1410233

55、222301961001402002年年 17 12 26 442 312 289 144 204 2003年年 2015295804354002253002004年年 2218327045764843243962005年年 2522358757706254845502006年年 2726401 0801 040729676702合计合计 1571212414 5083 6393 0132 0472 451将相关数据代入上述方程组，得到解方程组，得到所以，二元线性回归方程为 。121212241 10157121450815730132451363912124512047abbabbabb12

56、6.92850.350 00.9650abb126.92850.350.965yxx（2）显著性检验。利用复相关系数检验回归方程整体显著性，即 取 0.05，并计算df n k 1 7（k为自变量个数），查相关系数临界值表得到R0.05，7 0.758。因为RR，说明x1、x2与y的线性关系显著。2211222222()11()6 7776.92852410.3545080.965 363910.990 76 77710(24.1)iiiiiiiiiyaybx ybx yyyRyyyny（3）预测分析 点预测。当广告费支出x1,0 35万元，营业网点数x2,0 34个时，有 6.928 5 0

57、.35 x1 0.965x2 6.928 5 0.35 35 0.965 34 51.988 5（百万元） y 区间预测 计算估计标准误差，即 取 0.05，df n 3 7，查t分布表得到t/2 d f t/2 7 2.365。所以，预测区间为 t/2 7 S 51.988 5 2.365 1.594 5 51.988 5 3.771 0 2()17.79651.594537iyySn y即： 当在95%的把握程度下，如果广告费支出达到35万元，营业网点有34个，则商品销售额的置信区间为（48.22，55.76）百万元。 回归分析预测法是利用变量间因果关系进行预测的重要方法之一，除了线性回归

58、分析预测法外，还有非线性回归分析预测法。 在物流需求预测中，由于企业中的物流人员经常需要的是编制短期计划所需的物流需求预测，因此常用时间序列法和回归分析法。 为了使预测结果比较切合实际，提高预测质量，为决策和计划提供可靠的依据，通常将定性预测和定量预测两种预测方法相结合。 $#YWUSPNLJHFDBzwusqomkigeb97531+)*%#YWUSQOMKHFDBzxvtrpmkigeca8631+)*%!ZXVSQOMKIGECAxvtrpnljheca86420-(%!ZXVTRPNKIGECAywuspnljhftrpnligeca86420)*%!ZXVTROMKIGECAywtr

59、pnljhfdb86420-(&$ZXVTRPNLJHECAywusqoljhfdb97530-(&$#YWUSPNLJHFDBzwusqomkigeb97531+)*$#YWUSQOMKHFDBzxvtrpmkigeca8631+)*%!ZXVSQOMKIGEChfdb97530-(&$#YWURPNLJHFDBzwusqomkigeb97531+)*$#YWUSQOMKHFDBzxvtromkigeca8631+)*%!ZXUSQOMKIGECzxvtrpnljheca86420-*%!ZXVTRPNKIGECAywurpnljhfdb96420-(&$#YVTRPNLJHFCAywusq

60、omkhfdb97531-(&$#YWUSQNLJHFDBzxvsGECAywusqoljhfdb97520-(&$#YWURPNLJHFDBzwusqomkigdb97531+)*$#YWUSQOMJHFDBzxvtromkigeca8631+)*%!ZXUSQOMKIGECzxvtrpnljgeca86420-*%!ZXVTRPNKIGECAywurpnljhfdb96420-(&$#XVTRPNLJHFCAywusqomkhfdb97531-(&$#YW0-(&!ZXVTRPNLJGECAywusqnljhfdb97520-(&$#YWURPNLJHFDBywusqomkigdb9753