人大(王燕)时间序列课后习题答案)2-5(含上机的)

1、第二章P341、（1）因为序列具有明显的趋势，所以序列非平稳。 （2）样本自相关系数： 35 29.75 25.9167 21.75 (4)=17.25 (5)=12.4167 (6)=7.25 =0.85（0.85） =0.7405（0.702） =0.6214（0.556） =0.4929（0.415） =0.3548（0.280） =0.2071（0.153）注：括号内的结果为近似公式所计算。 （3）样本自相关图：AutocorrelationPartial CorrelationAC PAC Q-Stat Prob . |*| . |*|10.8500.85016.7320.000 .

2、 |* | . *| . |20.702-0.07628.7610.000 . |* | . *| . |30.556-0.07636.7620.000 . |* | . *| . |40.415-0.07741.5000.000 . |*. | . *| . |50.280-0.07743.8000.000 . |* . | . *| . |60.153-0.07844.5330.000 . | . | . *| . |70.034-0.07744.5720.000 . *| . | . *| . |8-0.074-0.07744.7710.000 . *| . | . *| . |9-0.

3、170-0.07545.9210.000 .*| . | . *| . |10-0.252-0.07248.7130.000 .*| . | . *| . |11-0.319-0.06753.6930.000 *| . | . *| . |12-0.370-0.06061.2200.000该图的自相关系数衰减为0的速度缓慢，可认为非平稳。4、 LB(6)=1.6747 LB(12)=4.9895 (6)=12.59 (12)=21.0显然，LB统计量小于对应的临界值，该序列为纯随机序列。第三章P971、解： 2、解：对于AR（2）模型：解得：3、解：根据该AR(2)模型的形式，易得： 原模型可

4、变为： =1.9823 4、解：原模型可变形为： 由其平稳域判别条件知：当，且时，模型平稳。 由此可知c应满足：，且 即当1c0时，该AR(2)模型平稳。 5、证明：已知原模型可变形为： 其特征方程为： 不论c取何值，都会有一特征根等于1，因此模型非平稳。6、解：（1）错，。 （2）错，。 （3）错，。 （4）错， （5）错，。7、解： MA(1)模型的表达式为：。8、解： 原模型可变为： 显然，当能够整除10.5B时，模型为MA(2)模型，由此得B2是0的根，故C0.275。9、解： 10、解：（1） 即 显然模型的AR部分的特征根是1，模型非平稳。（2） 为MA(1)模型，平稳。 11、解

5、：（1），模型非平稳； 1.3738 -0.8736 （2），模型平稳。 0.6 0.5 （3），模型可逆。 0.450.2693i 0.450.2693i （4），模型不可逆。 0.2569 -1.5569 （5），模型平稳；0.7 ，模型可逆；0.6 （6），模型非平稳。 0.4124 -1.2124 ，模型不可逆；1.112、解： ，13、解： 14、证明：； 15、解：（1）错；（2）对；（3）对；（4）错。16、解：（1）， 已知AR(1)模型的Green函数为：， 9.9892-1.96*，9.98921.96* 即3.8275,16.1509 （2） 10.045-1.96，10

6、.0451.96* 即3.9061,16.1839习题41、所以，在中与前面的系数均为。2、由 代入数据得 解得 3、（1） （2）利用且初始值进行迭代计算即可。另外， 该题详见Excel。11.79277 （3）在移动平均法下: 在指数平滑法中: 5、由 代入数据得 解得 z-c(10,11,12,10,11,14,12,13,11,15,12,14,13,12,14,12,10,10,11,13)6、方法一：趋势拟合法income-scan(习题4.6数据.txt)ts.plot(income)由时序图可以看出，该序列呈现二次曲线的形状。于是，我们对该序列进行二次曲线拟合：t-1:leng

7、th(income)t2-t2z-lm(incomet+t2)summary(z)lines(z$fitted.values, col=2)方法二：移动平滑法拟合选取N=5income.fil-filter(income,rep(1/5,5),sides=1)lines(income.fil,col=3)7、（1）milk-scan(习题4.7数据.txt)ts.plot(milk)从该序列的时序图中，我们看到长期递增趋势和以年为固定周期的季节波动同时作用于该序列,因此我们可以采用乘积模型和加法模型。在这里以加法模型为例。z-scan(4.7.txt)ts.plot(z)z-ts(z,star

8、t=c(1962,1),frequency=12)z.s-decompose(z,type=additive) /运用加法模型进行分解z.1-z-z.s$seas /提取其中的季节系数，并在z中减去（因为是加法模/型）该季节系数ts.plot(z.1)lines(z.s$trend,col=3)z.2-ts(z.1)t-1:length(z.2)t2-t2t3-t3r1-lm(z.2t)r2-lm(z.2t+t2)r3-lm(z.2t+t2+t3)summary(r1)summary(r2)summary(r3) #发现3次拟合效果最佳，故选用三次拟合ts.plot(z.2)lines(r3$

9、fitt,col=4)pt-(length(z.2)+1) : (length(z.2)+12)pt1-pt #预测下一年序列pt2-pt2pt3-pt3pt-matrix(c(pt1,pt2,pt3),byrow=T,nrow=3)/*为预测时间的矩阵。*/p-r3$coef2:4%*%pt+r3$coef1/*矩阵的乘法为%*%；coef【1】为其截距项，coef【2：4】为其系数*/p1-z.s$sea1:12+p/*加回原有季节系数，因为原来是加法模型*/ts.plot(ts(z),xlim=c(1,123),ylim=c(550,950)lines(pt1,p1,col=2)#包含季

10、节效应的 SARIMA模型z-scan(4.7.txt)ts.plot(diff(z)sq-diff(diff(z),lag=12) /*12步差分*/par(mfrow=c(2,1) acf(sq,50) pacf(sq,50)#观察上图，发现ACF图12阶处明显，24阶处即变到置信区间内。#而PACF图12阶，24阶，36阶处有一个逐渐递减过程，可认为#拖尾，故可以考虑对季节效应部分采用MA(1)模型#同时，ACF图在第一阶处显著后即立刻变动到置信区间内，具有#截尾性质，PACF图在第5、6阶时变动到置信区间外，可以考虑#使用MA（1）模型，故综合可采用乘积模型#即ri1、ma1模型乘以季

11、节因素result-arima(z,order=c(0,1,1),seasonal=list(order=c(0,1,1),period=12)/*季节因素里的order为阶数的意思，与前面的airma模型的阶数含义同*/tsdiag(result)/诊断#下图为预测后的图4.8z-scan(4.8.txt)adf.test(z) #单位根检验。比较科学的定量的方法#其原假设：具有单位根，即不平稳。此题中接受备则假设：平稳。指数平滑预测ffe-function(z,a) #定义指数平滑预测。其中a为平滑项y-c()y-z1for(i in 1:length(z)y-c(y, a*zi+(1-a

12、)*yi)return(y)y-ffe(z,0.6) #执行上述定义的function ts.plot(z) lines(y,col=3) ylength(y)简单移动平均z.1-filter(z,rep(1/12,12),side=1) #side=1是指将所有算不出的序列值都空到最前面去，而在尾部没有空值。z.1-c(NA,z.1)ts.plot(z)lines(z.1,col=3)meand-function(z,z.1,n) #预测函数。以12为周期。依次为原始数据，平滑值，预测步数y-z.1length(z.1)z.2-z(length(z)-10):length(z)for(i i

13、n 1:n) m-sum(rep(1/12,12-i)*z.2i:length(z.2) n-sum(rep(1/12,i)*y) y-c(y,m+n) #一直重复：预测，原始数列取代一个，预测数列拿来一个return(y)y-meand(z,z.1,11)y-c(z.1,y) ts.plot(z,xlim=c(0,205) lines(y,col=3)#SARIMApar(mfrow=c(2,1)ds-diff(z)acf(ds,40)pacf(ds,40)#可以看出有一些不明显的周期性，故采用sarima拟合result-arima(z,order=c(2,1,0),seasonal=li

14、st(order=c(1,0,0),period=12)#在季节部分很少出现2以上的数字（指seasonal中的order部分） result-arima(z,order=c(2,1,0),seasonal=list(order=c(1,0,1),period=12) result-arima(z,order=c(4,1,0),seasonal=list(order=c(1,0,1),period=12),fixed=c(NA,NA,0,NA,NA,NA) #观察图，发现第三项在置信区间内，故认为可能为限定的sarima模型。最后两个NA指季节指数中的sar1和sma1.#第三个的aic值最小

15、，即模型拟合效果最好tsdiag(result) #检验通过1、（1）判断序列的平稳性该序列时序图如图1所示：时序图显示该序列有显著的变化趋势，为典型的非平稳序列。（2）对原序列进行差分运算：对原序列进行1阶差分运算，运算后序列时序图如图2所示：时序图显示差分后序列在均值附近比较平稳的波动。为了进一步确定平稳性，考察差分后序列的自相关图，如图三所示：自相关图显示差分后序列不存在自相关，所以可以认为1阶差分后序列平稳，从图中我们还可以判断差分后序列可以视为白噪声序列。（3）对白噪声平稳差分序列拟合AR模型原序列的自相关图和偏自相关图如图4：图中显示序列自相关系数拖尾，偏自相关系数1阶截尾，实际上

16、我们用ARIMA（1，0，0）模型拟合原序列。在最小二乘估计原理下，拟合结果为：（4）对残差序列进行检验：残差白噪声检验：参数显著性检验：图中显示：延迟6阶和12阶的P值均大于0.05，可以认为该残差序列即为白噪声序列，系数显著性检验显示两参数均显著。这说明ARIMA（1，0，0）模型对该序列建模成功。（5）模型的预测：估计下一盘的收盘价为：2、（1）绘制时序图：时序图显示该序列具有长期递增趋势和以年为周期的季节效应。（2）差分平稳化对原序列作1阶差分，希望提取原序列的趋势效应，差分后序列时序图：3、模型定阶考察差分后序列相关图和偏自相关图的性质，进一步确认平稳性判断，并估计拟合模型的阶数。自

17、相关图和偏自相关图显示延迟12阶自相关系数和偏自相关系数大于2倍标准差范围，说明差分后序列中仍有非常显著的季节效应。延迟1阶的自相关系数和偏自相关系数也大于2倍的标准差，这说明差分后序列还具有短期相关性。根据差分后序列自相关图和偏自相关图的性质，尝试拟合ARMA模型，但拟合效果均不理想，拟合残差均通不过白噪声检验。所以我们可以考虑建立乘积模型： ：（4）参数估计使用最小二乘法估计方法，得到该模型的估计方程为：（5）模型的检验对拟合模型进行检验，检验结果显示该模型顺利通过了残差白噪声检验（图21）和参数显著性检验（图22）。白噪声检验（图21）参数显著性检验（图22）（6）模型预测下一年度该城市月度婴儿出生率预测如下表：月份123456789101112预测值27.61127.60427.89527.76227.88127.80527.84827.83627.8127.84227.74827.7883、（1）展开原模型，等价形式为：即 所以 （2） 4、（1）平稳性检验：从该序列时序图中可以看到该序列为非平稳序列。（2）模型拟合：ARCH过程检验：异方差怀特检验：DW=2.05 序列中残差不存在自相关；怀特检验之后也不存在异方差；ARCH LM检验之后也不存在ARCH过程。所以确定该模型为： （3）预测：19391945年英国绵羊的数量预测如下表：年份193919401941194