利用SPSS和Matlab进行时间序列预测.pdf

《利用SPSS和Matlab进行时间序列预测.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实习指导——《计量地理学》（徐建华，华东师范大学）§7.利用SPSS和Matlab进行时间序列预测1.移动平均和滑动平均计算例1：下表给出了某地区1990～2004年粮食产量数据（表1）。试分别用Matlab和SPSS软件，对该地区的粮食产量进行移动平均和和滑动平均计算。表1某地区1990～2004年粮食产量及其平滑结果粮食产量y（单移动平均滑动平均年份自然序号4位：10t）三点移动五点移动三点滑动五点滑动199013149.44199123303.663154.47199233010.303010.303141.193242.44199343109.613154.473253.0432

2、63.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.613733.693618.814052.514019.782000114253.253914.723692.894121.454075.782001124101.504052

3、.513892.784158.214148.582002134119.884121.454019.784160.014227.012003144258.654158.214075.784260.112004154401.794160.014148.58利用spss进行移动平均计算主要有以下步骤：（1）在菜单中依次选择transform->Createtimeseries…,在弹出的对话框中的单击function下面的下拉条，选择Priormovingaverage方法，span框中输入数值为3（表示进行三点滑动平均）。32实习指导——《计量地理学》（徐建华，华东师范大学）（2）在左边变量

4、框中选择“粮食产量”变量，点击“添加”按钮，则在右边的新变量框中自动生成一个新变量名，在下面的name框中可以对此变量名进行改变，本例采用默认变量名。（3）在span框中输入数值5，重复步骤（2）。（4）在function中选择Centeredmovingaverage方法，span框中分别输入3和5，分别重复步骤（2）。，此时新变量框有四个变量，前两个分别储存三点和五点移动平均计算后的结果，后两个分别储存三点和五点滑动平均计算后的结果。（5）点击ok按钮进行计算。在原始数据表后新加的四列变量极为所求的移动和滑动平均结果。2.季节性预测4例2：某旅游景点2002～2004年各季度客流量y

5、i（10人次）如表2所示，下面我们用季节性预测方法，预测该旅游景点2005年各季度的客流量。4表2某市2002～2004年各季度客流量及其三点滑动平均值（10人次）年份季度t游客人数三点滑动平均11260.0022375.00325.00200233340.00312.6744223.00279.3333实习指导——《计量地理学》（徐建华，华东师范大学）25275.00303.3326412.00346.33200337352.00331.6748231.00290.0019287.00315.33210428.00359.672004311364.00345.00412243.00使用

6、季节变动预测法对该旅游景点2005年各季度的客流量的步骤为：（1）对原始数据求三次滑动平均值。该过程可以使用spss来实现，方法如例1。（2）用二次指数平滑法求预测模型系数。二次指数平滑法的预测公式为yˆ=a+bkt+ktt式中：k代表从基期t到预测期的期数。(1)(2)a=2S−Stttα()(1)(2)b=S−Sttt1−α式中：(1)(1)S=αy+(1−α)Sttt−1(2)(1)(2)S=αS+(1−α)Sttt−1指数平滑和方程系数可以使用Matlab计算，代码如下（已经输入原始数据34实习指导——《计量地理学》（徐建华，华东师范大学）矩阵a，平滑系数为alf＝0.2，m为原

7、始数据矩阵的行数）：[m,n]=size(a)alf=0.2;%设定平滑系数的数值forj=1:2s(1,j)=a(1,1);endfori=2:mforj=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);％将循环式中每次运