时间序列分析试验报告.doc

时间序列分析实验报告（一）实验目 的时间序列是按时间次序排列的随机变量序列，任何时间序列经过合理的函数变换后都可以被认为是由三个部分叠加而成，这三个部分别是趋势项部分、周期项部分和随机噪声项部分，从时间序列中把这三个部分分解出来是时间序列分析的首要任务。本实验目的学习时间序列的分解方法。实验内容下面的表7.1.1中的数据是某城市19911996年中每个季度的民用煤消耗量（单位：吨）。数据图形由图1.1.2给出。表7.1.1某城市居民季度用煤消耗量（单位：吨）年份1季度2季度3季度4季度年平均19916878.45343.74847.96421.95873.0 19926815.45532.64745.66406.25875.0 19936634.45658.54674.86445.55853.319947130.25532.64989.66642.36073.719957413.55863.14997.46776.16262.619967476.55965.55202.16894.16384.5季平均7058.15649.34909.66597.7（1）由表8.1.1中每年每季的数据计算年平均值与季平均值，并绘出19911996年中每个季度的数据的散点图。（2）用回归直线趋势法对序列进行分解。（3）若1997年四季的数据分别为：7720.5 5973.3 5304.4 7075.1，运用（2）对1997年数据作预测并分析误差。实验步骤(1)首先做理论准备。每个时间序列，或经过适当的函数变换的时间序列，都可以分解成三个部分的叠加(1)其中是趋势项，是季节项，是随机项。然后，利用MATLAB软件，编写程序进行计算，即A=6878.4 5343.7 4847.9 6421.96815.4 5532.6 4745.6 6406.26634.4 5658.5 4674.8 6445.57130.2 5532.6 4989.6 6642.37413.5 5863.1 4997.4 6776.17476.5 5965.5 5202.1 6894.1;A1=mean(A) %计算季平均值B=A;A2=mean(B) %计算年平均值plot(B(:),+-) %画出每个季度的散点图则得:季平均值为：7058.1 5649.3 4909.6 6597.7年平均值为:5873.0 5875.0 5853.3 6073.7 6262.5 6384.5每个季度的数据的散点图：图1 城市居民季度用煤消耗量散点图（2）分解回归直线趋势。由于数据有缓慢的上升趋势，可以试用回归直线表示趋势项，这时认为（满足一元线性回归模型在Matlab命令窗口中继续输入下列命令：polyfit(1:24,B(:),1)输出： ans =1.0e+003 *0.0219 5.7801所以得： ，这时，趋势项的估计值是回归直线：，利用原始数据减去趋势项的估计后得到的数据基本只含有季节项和随机项了。分解季节项：用第k季度的平均值作为季节项的估计。如果用分别表示第j年第k个季度的数据和趋势项，则时刻（j,k）的时间次序指标为.在Matlab命令窗口中继续输入下列命令：dx=B(:)-(5780.1+21.9*(1:24)C=dx(:,1:4);dx(:,5:8);dx(:,9:12);dx(:,13:16);dx(:,17:20);dx(:,21:24);s=mean(C) %季节项估计则得 s = 1.0e+003 * 1.0371 -0.3936 -1.1552 0.5110即季节项估计为分解随机项：利用原始数据减去趋势项的估计和季节项的估计后得到的数据就是随机项的估计.在Matlab命令窗口中继续输入下列命令：for j=1:6for k=1:4St(k+4*(j-1)=s(k); %求季节项值StendendRt=dx-St; %求随机项估计plot(1:24,St,*-,1:24,Rt,-) %画出季节项和随机项图形图2季节项和随机项散点图预测:为得到1997年的预报值，可以利用公式， 这里，是用例中的24个观测数据对第个数据的预测值，利用MATLAB编写命令：for i=25:28 m=5780.1+21.9*(i)+s(i-24) %计算1997年四个季度的预测值end得到结果：m=7364.7 5955.9 5216.2 6904.3将1997年的预测值与真实值比较