利用spss和matlab进行时间序列预测

《利用spss和matlab进行时间序列预测》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、《计量地理学》（徐建华，高等教育出版社，2005）配套实习指导§7.利用SPSS和Matlab进行时间序列预测1.移动平均和滑动平均计算例1：表1给出了某地区1990一2004年粮食产量数据。试分别用Matlab和SPSS软件，对该地区的粮食产量进行移动平均和和滑动平均计算。表1某地区1990一2004年粮食产量及其平滑结果年份自然序号粮食产量y/104t移动平均滑动平均三点移动五点移动三点滑动五点滑动199013149.44199123303.663154.47199233010.303010.303

2、141.193242.44199343109.613154.473253.043263.32199453639.213141.193334.213295.88199563253.803253.043242.443453.173461.80199673466.503334.213263.323520.073618.81199783839.903453.173295.883733.693692.89199893894.663520.073461.803914.723892.781999104009.6137

3、33.693618.814052.514019.782000114253.253914.723692.894121.454075.782001124101.504052.513892.784158.214148.582002134119.884121.454019.784160.014227.012003144258.654158.214075.784260.112004154401.794160.014148.58利用SPSS进行移动平均计算主要有以下步骤：（1）在菜单中依次选择Transform->

4、Createtimeseries…,在弹出的对话框中的单击Function下面的下拉条，选择Priormovingaverage方法，span框中输入数值为3（表示进行三点滑动平均）。36《计量地理学》（徐建华，高等教育出版社，2005）配套实习指导（2）在左边变量框中选择“粮食产量”变量，点击“添加”按钮，则在右边的新变量框中自动生成一个新变量名，在下面的Name框中可以对此变量名进行改变，本例采用默认变量名。（3）在Span框中输入数值5，重复步骤（2）。（4）在Function中选择Centere

5、dmovingaverage方法，Span框中分别输入3和5，分别重复步骤（2）。，此时新变量框有四个变量，前两个分别储存三点和五点移动平均计算后的结果，后两个分别储存三点和五点滑动平均计算后的结果。（5）点击OK按钮进行计算。在原始数据表后新加的四列变量即为所求的移动和滑动平均结果。2.季节性预测例2：某旅游景点2002一2004年各季度客流量yi（104人次）如表2所示，下面我们用季节性预测方法，预测该旅游景点2005年各季度的客流量。表2某市2002一2004年各季度客流量及其三点滑动平均值（10

6、4人次）年份季度t游客人数三点滑动平均200211260.0022375.00325.0036《计量地理学》（徐建华，高等教育出版社，2005）配套实习指导33340.00312.6744223.00279.33200325275.00303.3326412.00346.3337352.00331.6748231.00290.00200419287.00315.33210428.00359.67311364.00345.00412243.00使用季节变动预测法对该旅游景点2005年各季度的客流量的步骤为

7、：（1）对原始数据求三次滑动平均值。该过程可以使用SPSS来实现，方法如例1。（2）用二次指数平滑法求预测模型系数。二次指数平滑法的预测公式为式中：k代表从基期t到预测期的期数。式中：36《计量地理学》（徐建华，高等教育出版社，2005）配套实习指导指数平滑和方程系数可以使用Matlab计算，代码如下（已经输入原始数据矩阵a，平滑系数为alf＝0.2，m为原始数据矩阵的行数）：[m,n]=size(a)alf=0.2;%设定平滑系数的数值forj=1:2s(1,j)=a(1,1);endfori=2:m

8、forj=1:2ifj==1s(i,j)=alf*a(i,1)+(1-alf)*s(i-1,j);elses(i,j)=alf*s(i,j-1)+(1-alf)*s(i-1,j);endendend%计算预测方程的系数a(t),b(t)temp=alf/(1-alf);％将循环式中每次运算的常数提取出来计算，以减少计算量fori=1:mat(i,1)=2*s(i,1)-s(i,2);bt(i,1)=temp*(s(i,1)-s(i,2));