时间序列ARMA模型及分析

1、ARMA模型及分析本次试验主要是通过等时间间隔，连续读取70个某次化学反应的过程数据，构成一个时间序列。试对该时间序列进行模型拟合以及模型的优化，最后进行预测。以下本次试验的数据：表1连续读取70个化学反应数据资料来源：下面的分析及检验、预测均是基于上述数据进行的，本次试验是在6.上0完成的。一、序列预处理由于只有对平稳的时间序列才能建立ARMA模型，因此在建立模型之前有必要对序列进行预处理，主要包括了平稳性检验和纯随机检验。图化学反应过程时序图序列时序图显示此化学反应过程无明显趋势或周期，波动稳定。见图1。AutocorrelationPartialCorrelationACPACCl-St

2、atProb11111-0.390-0.39011.1030.0011二1120.3040.18017.9700.000匚1113-0.1660.00220.0320.0001111140.071-0.0U20.4140.000|E1115-0.097-0.06921.1440.001|111匚16-0.04-7-0.121213190.0021111170.0350.02021.4190.003111110-0.0430.00521.5720.00611119-0.005-0.05621.5740.0101111100.0U0.00421.5920.0171J111110.1100.1432

3、2.6240.020111112-0.069-0.00923.0350.027111J1130.1400.09224.9720.02311111140.0360.16725.0880.034111115-0.007-0.00125.0920.049111Z3160.1730.22127.8850.0331匚111I17-0.1110.05329.0640.034111匚1180.020-0.10529.1030.0471111119-0.04-70.042293240.061111200.0160.050293500.08111111210.0220.05629.4020.1051111122

4、-0.079-0.04230.0520.117111匚123-0.010-0.13730.0620.143111匚124-0.073-0.16330.6480.164111c125-0.020-0.00530.6900.19911111260.041-0.00430.8870233111匚127-0.022-0.12730.9450.2731111280.0090.02531.0930.2791111290.0160.01731.925032311|L1300.00斗-0.09331.9270.3711111310.0050.00431.9300.4201111132-0.025-0.0573

5、2.0110.4-66图化学反应过程相关图和统计量从图2的序列的相关分析结果:1.可以看出自相关系数始终在0周围波动，判定该序列为平稳时间序列看统计量的值：该统计量的原假设为的期，期期的自相关系数均等于0备择假设为自相关系数中至少有一个不等于，因此如图知，该值在滞后23期是都为所以拒接原假设即序列是非纯随机序列,即非白噪声序列（因为序列值之间彼此之间存在关联,所以说过去的行为对将来的发展有一定的影响,因此为非纯随机序列,即非白噪声序列。二、模型识别由于检验出时间序列是平稳的，且是非白噪声序列，因此可以建立模型，在建立模型之前需要识别模型阶数即确定阶数。阶数确定要借助于时间序列的相关图，即序列的

6、自相关函数和偏自相关函数，并根据他们之间的理论模式进行阶数最后的确定。下面给出自相关函数和偏自相关函数之间的理论模式：表时间序列的与理论模式自相关系数偏相关系数模型定阶拖尾P阶截尾AR(p模型q阶截尾拖尾MA(q)模型拖尾拖尾ARMA(p,q)模型根据时间序列的相关图图2显示的自相关系数的2阶截尾的性质以及偏自相关系数1阶截尾性，我们尝试拟合ARMA(1,2)模型。三、模型参数估计在识别了模型的形式后，我们就可以使用Eviews估计方程参数。下面就对ARMA(1,2)模型其参数估计的结果。R-squaredAdjustedR-squaredS.E.ofregressionSumsquaredr

7、esidLoglikelihoodF-statisticProb(F-statistic)0.1940750.157715R-squaredAdjustedR-squaredS.E.ofregressionSumsquaredresidLoglikelihoodF-statisticProb(F-statistic)0.1940750.15771510999947864.908-261.30085.2442580.002651Meandependsnt炖S.D.dependentvarAkaikeinfocriterionSchwarzcriterionH日nnan-Quinncriter.D

8、urbin-Watsonstat511884111.905627.6893797.8193927.7412611.877349VariableCoefficientStd.Errort-StatisticProL.C51.194061.21312742.200070.0000AR(1)-0.3049250359240-0.3430050.3991MA(1)-0.0506990350557-0.144R230.8355MA(20.2506640.1656961.51279301352InvertedARRoots-.30InvertedMARoots,03-.50i,03+.50i图(，)模型估

9、计结果以上就是拟合ARMA(1,2)的结果，我们用来表示时间序列，于是我们基于上述结果写出ARMA(1,2)的估计结果：(1-0.304925B)yield(=51.19406,(1-0.050699B,0.250664B2)w(-22-R20.1948R0.1577AIC准则7.6898SC准则7.8193对于ARMA(1,2)模型估计，其命令形式为：Isyieldcar(1)ma(1)ma(2)。四、模型诊断检验ARMA模型参数估计后，应该检验模型的确认是否正确，通常是对模型的残差序列进行白噪声检验。图4模型图4模型的残差相关图和统计量对残差序列进行白噪声检验，可以看出和都没有显著异于零，

10、统计量的值都远远大于，因此可以认为残差序列为白噪声序列，模型信息提取比较充分。常数和滞后一阶参数的值都很小，参数显著；因此整个模型比较精简，模型较优。在模型检验之后，我们还可以对模型优化，模型优化的主要判断标准就是AIC准则和SC准则。在几个模型都符合要求，且也都有效参数显著，这个时候我就要通过比较AIC准则和SC准则，从而来确定最终的模型，当然是AIC准则和SC准则越小越优。五、模型预测通过上述的分析我们知道，模型ARMA(1,2)是合适的，因此，我们就基于它来进行预测。XF2SF在这我们利用模型对65到70的这几个数据进行预测，预测结果如下图：XF2SFForecast:XFActual:XForecastsample:6570Ineludedobservations:6RootMeanSquaredError12.85548MeanAbsoluteError9.475397MeanAbs.PercentError30.97303TheilInequalityCoefficient0.131232BiasProportion0.181260VarianeeProportion0.586456CovarianeeProportion0.232284图5