具有季节性特点的时间序列的预测

3.2 具有季节性特点的时间序列的预测

这里提到的季节，可以是自然季节，也可以是某种产品的销售季节等。显然，在现实的经济活动中，表现为季节性的时间序列是非常多的。比如，空调、取暖设备、季节性服装的生产与销售所产生的数据等。对于季节性时间序列的预测，要从数学上完全拟合其变化曲线是非常困难的。但预测的目的是为了找到时间序列的变化趋势，尽可能地做到精确。从这个意义上来讲，可以有多种方法，下面介绍其中一种，即所谓季节系数法。季节系数法的具体计算步骤如下：

1．收集m年的每年各季度或者各月份（每年n个季度）的时间序列样本数据xij。 2．计算每年所有的季度或所有月份的算术平均值，即：

1mn

xij,

ki 1j 1

m

k mn

3．计算同季度或同月份数据的算术平均值.j

x

i 1

ij

,j 1,2, ,n

4．计算季节系数或月份系数 j .j/。其中j 1,2, ,n为季度或者月份的序号。 5．预测计算。当时间序列是按季度列出时，先求出预测年份（下一年）的年加权平均：

ym 1

n

w1y1 w2y2 wmym

w1 w2 wm

式中，yi

x

j 1

ij

为i年份的年合计数：wi为i年份权数，按自然数列取值。再计算预

测年份的季度平均值m 1:m 1 ym 14。最后，预测年份第i季度的预测值为：

ym 1,i m 1 i

季节系数法的Matlab程序如下。 funjie.m

%简单季节系数法，文件名funjie.m

function JiJie=funjie(x) %输入m年，每年n个季节的历史数据 [m,n]=size(x);

BarX=mean(mean(x)) %计算所有数据的算术平均值 BarXj=mean(x) %计算同季节的算术平均值 Betaj=BarXj./BarX %计算季节系数 y1=[1:m];

y=y1*sum(x,2)/sum(y1) %计算预测下一年的年加权平均值 y2=y/n %计算预测年份的季节平均值 y3=y2*Betaj %预测年份的季节预测值 end

【例3-11】某商店某类商品1999-2003年各季度的销售额如表3-6所示。试预测2004

年第3、4季度的销售额。

手工计算如下： 1．

3967761

198338.05； 20

2．季节系数 i：

3．y2004

143799.8198713.4

0.7248， 2 1.0018

198338.05198338.05269374.2181664.8 3 1.3578， 4 0.9517

198338.05198338.05

1 774656 2 790177 3 834871 5 825665 803333

1 2 3 4 5

1

4．2004年季度平均数：2004 803333/4 200833.25 于是，2004年第1、2、3、4季度预测值分别为：

y2004,1 200833.25 0.7248 145563.94，y2004,2 200833.25 1.0018 201194.75 y2004,3 200833.25 1.3578 272691.39，y2004,4 200833.25 0.9157 183903.01

利用Matlab程序funjie.m计算如下： >>

x=[137920,186742,274561,175433;142814,198423,265419,183521;131002,193987,247556,169847;157436,200144,283002,194319;149827,214301,276333,185204]; >> funjie(x)

BarX = 1.9839e+005 BarXj = 1.0e+005 *

1.4380 1.9872 2.6937 1.8166 Betaj = 0.7248 1.0017 1.3578 0.9157 y = 803341 y2 = 2.0084e+005 y3 = 1.0e+005 *

1.4557 2.0117 2.7269 1.8390

如果在求出季节系数后，还已知预测年份某季节的实际值，那么，可以利用季节系数进行调节，以使得在随后的预测中能更准确地反映其变化趋势。假如2004年第二季度的实际销售额是223045元，二季度季节系数为1.0017，则2004年第三季度预测值修正为：

y2004,3

223045

1.3578 3.0234e 005

1.0017

2004年第四季度的预测值修正为：

y2004,4

223045

0.9157 2.0389e 005

1.0017

【例3-12】某服装企业2001-2003年每月针织夏衣销售额如表3-7所示。预测2004年5、6、7月销售额。

表3-7 某服装企业2001-2003年每月针织夏衣销售额 （单位：万元）

可利用Matlab程序funjie.m计算如下： >>

x=[240,360,640,1080,2720,4480,4040,2360,800,520,360,280;240,460,860,1240,2800,5060,4980,2440,760,480,260,300;340,420,1080,1760,2580,5620,5600,2100,1340,840,420,220];

>> funjie(x)

BarX = 1.6689e+003 BarXj = 1.0e+003 * Columns 1 through 7

0.2733 0.4133 0.8600 1.3600 2.7000 5.0533 4.8733 Columns 8 through 12

2.3000 0.9667 0.6133 0.3467 0.2667 Betaj =

Columns 1 through 7

0.1638 0.2477 0.5153 0.8149 1.6178 3.0280 2.9201 Columns 8 through 12

1.3782 0.5792 0.3675 0.2077 0.1598 y = 2.0767e+004 y2 = 1.7306e+003 y3 = 1.0e+003 * Columns 1 through 7

0.2834 0.4286 0.8918 1.4103 2.7998 5.2401 5.0534 Columns 8 through 12

2.3850 1.0024 0.6360 0.3595 0.2765

本文来源：https://www.bwwdw.com/article/niv4.html