人大王燕时间序列课后习题答案25含上机的.docx
- 文档编号:568887
- 上传时间:2022-10-11
- 格式:DOCX
- 页数:31
- 大小:797.69KB
人大王燕时间序列课后习题答案25含上机的.docx
《人大王燕时间序列课后习题答案25含上机的.docx》由会员分享,可在线阅读,更多相关《人大王燕时间序列课后习题答案25含上机的.docx(31页珍藏版)》请在冰豆网上搜索。
人大王燕时间序列课后习题答案25含上机的
第二章P34
1、
(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)样本自相关图:
Autocorrelation
PartialCorrelation
AC
PAC
Q-Stat
Prob
.|*******|
.|*******|
1
0.850
0.850
16.732
0.000
.|*****|
.*|.|
2
0.702
-0.076
28.761
0.000
.|****|
.*|.|
3
0.556
-0.076
36.762
0.000
.|***|
.*|.|
4
0.415
-0.077
41.500
0.000
.|**.|
.*|.|
5
0.280
-0.077
43.800
0.000
.|*.|
.*|.|
6
0.153
-0.078
44.533
0.000
.|.|
.*|.|
7
0.034
-0.077
44.572
0.000
.*|.|
.*|.|
8
-0.074
-0.077
44.771
0.000
.*|.|
.*|.|
9
-0.170
-0.075
45.921
0.000
.**|.|
.*|.|
10
-0.252
-0.072
48.713
0.000
.**|.|
.*|.|
11
-0.319
-0.067
53.693
0.000
***|.|
.*|.|
12
-0.370
-0.060
61.220
0.000
该图的自相关系数衰减为0的速度缓慢,可认为非平稳。
4、
LB(6)=1.6747LB(12)=4.9895
(6)=12.59
(12)=21.0
显然,LB统计量小于对应的临界值,该序列为纯随机序列。
第三章P97
1、解:
2、解:
对于AR
(2)模型:
解得:
3、解:
根据该AR
(2)模型的形式,易得:
原模型可变为:
=1.9823
4、解:
原模型可变形为:
由其平稳域判别条件知:
当
,
且
时,模型平稳。
由此可知c应满足:
,
且
即当-1 (2)模型平稳。 5、证明: 已知原模型可变形为: 其特征方程为: 不论c取何值,都会有一特征根等于1,因此模型非平稳。 6、解: (1)错, 。 (2)错, 。 (3)错, 。 (4)错, (5)错, 。 7、解: MA (1)模型的表达式为: 。 8、解: 原模型可变为: 显然,当 能够整除1-0.5B时,模型为MA (2)模型,由此得B=2是 =0的根,故C=0.275。 9、解: : 10、解: (1) 即 显然模型的AR部分的特征根是1,模型非平稳。 (2) 为MA (1)模型,平稳。 11、解: (1) ,模型非平稳; 1.3738 -0.8736 (2) , , ,模型平稳。 0.6 0.5 (3) , , ,模型可逆。 0.45+0.2693i 0.45-0.2693i (4) , , ,模型不可逆。 0.2569 -1.5569 (5) ,模型平稳; 0.7 ,模型可逆; 0.6 (6) , , ,模型非平稳。 0.4124 -1.2124 ,模型不可逆; 1.1 12、解: , 13、解: 14、证明: ; 15、解: (1)错; (2)对;(3)对;(4)错。 16、解: (1) , 已知AR (1)模型的Green函数为: , [9.9892-1.96* ,9.9892+1.96* ] 即[3.8275,16.1509] (2) [10.045-1.96× ,10.045+1.96* ] 即[3.9061,16.1839] 习题4 1、 所以,在 中 与 前面的系数均为 。 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: length(income) t2<-t^2 z<-lm(income~t+t2) summary(z) lines(z$fitted.values,col=2) 方法二: 移动平滑法拟合 选取N=5 income.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,start=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<-t^2 t3<-t^3 r1<-lm(z.2~t) r2<-lm(z.2~t+t2) r3<-lm(z.2~t+t2+t3) summary(r1) summary(r2) summary(r3)##发现3次拟合效果最佳,故选用三次拟合 ts.plot(z.2) lines(r3$fitt,col=4) pt<-(length(z.2)+1): (length(z.2)+12) pt1<-pt##预测下一年序列 pt2<-pt^2 pt3<-pt^3 pt<-matrix(c(pt1,pt2,pt3),byrow=T,nrow=3)/*为预测时间的矩阵。 */ p<-r3$coef[2: 4]%*%pt+r3$coef[1]/*矩阵的乘法为%*%;coef【1】为其截距项,coef【2: 4】为其系数*/ p1<-z.s$sea[1: 12]+p/*加回原有季节系数,因为原来是加法模型*/ ts.plot(ts(z),xlim=c(1,123),ylim=c(550,950)) lines(pt1,p1,col=2) ##包含季节效应的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模型乘以季节因素 result<-arima(z,order=c(0,1,1),seasonal=list(order=c(0,1,1),period=12))/*季节因素里的order为阶数的意思,与前面的airma模型的阶数含义同*/ tsdiag(result)//诊断 ##下图为预测后的图 4.8 z<-scan('4.8.txt') adf.test(z)##单位根检验。 比较科学的定量的方法 ##其原假设: 具有单位根,即不平稳。 此题中接受备则假设: 平稳。 指数平滑预测 ffe<-function(z,a)##定义指数平滑预测。 其中a为平滑项 { y<-c() y<-z[1] for(iin1: length(z)) y<-c(y,a*z[i]+(1-a)*y[i]) return(y) } y<-ffe(z,0.6)##执行上述定义的function ts.plot(z) lines(y,col=3) y[length(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.1[length(z.1)] z.2<-z[(length(z)-10): length(z)] for(iin1: n) { m<-sum(rep(1/12,12-i)*z.2[i: 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) ##SARIMA par(mfrow=c(2,1)) ds<-diff(z) acf(ds,40) pacf(ds,40) ##可以看出有一些不明显的周期性,故采用sarima拟合 result<-arima(z,order=c(2,1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 人大 时间 序列 课后 习题 答案 25 上机