《时间序列分析第四章均值和自协方差函数的估计.ppt》由会员分享,可在线阅读,更多相关《时间序列分析第四章均值和自协方差函数的估计.ppt(73页珍藏版)》请在第壹文秘上搜索。
1、均值和自协方差函数的估计均值和自协方差函数的估计+均值的估计均值的估计+自协方差函数的估计自协方差函数的估计+白噪声检验白噪声检验+相合性相合性+中心极限定理中心极限定理+收敛速度收敛速度+的模拟计算的模拟计算X+AR,MA,ARMA模型的参数可以由自协方差函数唯一确定。+有了样本之后,可以先估计均值和自协方差函数。+然后由均值和自协方差函数解出模型参数。+均值和自协方差可以用矩估计法求。+还要考虑相合性,渐进分布,收敛速度等问题。+设 是平稳列 的观测。+的点估计为+把观测样本看成随机样本时记作大写的12,Nx xxtXtEX11NNkkxxN12,NXXX+设统计量 是 的估计,在统计学中
2、有如下的定义+1 如果 ,则称 是 的无偏估计。+2 如果当 则称 是 的渐 进无偏估计。+3 如果 依概率收敛到 ,则称 是 的相合估计。+4如果 收敛到 ,则称 是 的强相合估计。NNENE,.NNE NNN.Na sN+一般情况下,无偏估计比有偏估计来得好,对于由(1.1)定义的 。有所以 是均值 的无偏估计。_NX1111.NNNkkkE XE XNNNX+好的估计量起码应是相合的。否则,估计量不收敛到要估计的参数,它无助于实际问题的解决。+对于平稳序列 ,如果它的自协方差函数 收敛到零,则:tXk+利用切比雪夫不等式+得到 依概率收敛到 。于是 是 的相合估计。22()Pr(|)0.
3、(0)NNEXXNXNX+定理1.1 设平稳序列 有均值 和自协方差函数 。则 1 是 的无偏估计。2 如果 则 是 的相合估计。3 如果 还是严平稳遍历序列,则 是 的强相合估计。tXkNX0,kNXtXNX+第三条结论利用1.5的遍历定理5.1可得。一般地,任何强相合估计一定是相合估计。线性平稳列的均值估计是相合估计。ARMA模型的均值估计是相合估计。+若 。则+可以据此计算 的 置信区间。(1.3)其中的1.96也经常用2近似代替。212,(,)NXXX idd 2()(0,)NNXd N 95%1.96/,1.96/.NNXNXN+定理1.2 设 是独立同分布的 ,线性平稳序列 由 (
4、1.5)定义。其中 平方可和。如果 的谱密度 (1.6)在 连续,并且 则当 时,t2(0,)WNtX,tktkkXtZ ktX2()|,2ktkkftZ 0(0)0.fN ()(0,2(0)dNN XNf+当 绝对可和时,连续。+推论1.3 如果 和 成立,则当 时 并且 (1.7)k()f|kk 0kkN ()(0,2(0)dNN XNf 012(0)2.jjf+相合的估计量渐进性质除了是否服从中心极限定理外,还包括这个估计量的收敛速度。+收敛速度的描述方法之一是所谓的重对数律。+重对数律成立时,得到的收敛速度的阶数一般是+除了个别情况,这个阶数一般不能再被改进。2ln ln().NoN+
5、定理1.4 设 是独立同分布的 。线性平稳序列 由(1.5)定义。谱密度 。当以下的条件之一成立时:1 当 以负指数阶收敛于0.2 谱密度 在 连续。并且 对某个 成立。t2(0,)WNtX(0)0f|,kk()f0|rtE 2r+则有重对数律 (1.8)(1.9)易见重对数律满足时 不收敛。()2(0),.2 ln lnlim supNNNXfa sN()2(0),.2 ln lnlim infNNNXfa sN ln ln()(1)0(),()/(1)2ln lnNnNNXoXoNN+令 考虑AR(2)模型为模拟方便设 。()(1)(1)iiA zezez()ttAXB2122costtt
6、tXXX2(0,)tiidN1111,NNNNttttXXNN+为了观察 时 的收敛可以模拟L个值然后观察 的变化。+为了研究固定N情况下 的精度以至于抽样分布。可以进行M次独立的随机模拟,得到M个 的观察值。这种方法对于难以得到估计量的理论分布的情况是很有用的。N NX0,0,1,NXNn nLNXNX+自协方差估计公式及正定性+的相合性+的渐进分布+模拟计算kk+(2.2)样本自相关系数(ACF)估计为 (2.3)11()(),01,NkNNkjjkjxxxxkNNkk0,|1kkkN+估计 一般不使用除了 的估计形式:(2.4)因为:我们不对大的k值计算 更重要的是只有除以N的估计式才是
7、正定的。kNk11()()NkNNjjkjxxxxNkk+只要观测 不全相同则 正定。+令 记 (2.5)只要 不全是零则A满秩。12,Nx xx,1,2,()Nkjk jN.Njjyxx121123110000000NNNNNyyyyyyyyAyyy1TNA ANiy+事实上,设 则A矩阵左面会出现一个以 值开始非零的斜面。显然是满秩的。+故 不全相同时 正定。+作为 的主子式也是正定的。110,0.kkyyyky1,NxxN(1)nnNN+定理2.1 设平稳序列的样本自协方差函数 由式(2.2)或(2.4)定义。1 如果当 时,则对每个确定的k,是 的渐进无偏估计:kkk 0.kkk.li
8、mkkNE+2如果 是严平稳遍历序列。则对每个确定的k,和 分别是 和 的强相合估计:tXkkkk,.,.limlimkkkkNNa sa s+下面只对由(2.2)定义的样本自协方差函数证明定理2.1。对由(2.4)定义的 的证明是一样的。设 则 是零均值的平稳序列。利用 (2.7)k1.EX ttYX112111()()1().NN kNNNNkjjj kjjN kNNjj kj kjjYYXYYYYNNYYY YYYN+只考虑线性序列。设 是4阶矩有限的独立同分布的 实数列 平方可和。线性平稳序列 (2.8)t22(0,)(0).WNk,.tjtjjXtZ+有自协方差函数 (2.9)有谱密
9、度 (2.10)tX2kjjkj tX22()|.2ijjjfe+设自协方差函数列 平方可和。+设 为独立同分布的 。+令+定义正态时间序列 (2.11)(2.12)ktW(0,1)N4441/210421,()0EM001()().0jjtjtjttMWWj1(2),1jtjtjtjttRWj+定理2.2 设 是独立同分布的 。满足 。如果线性平稳序列(2.8)的谱密度(2.10)平方可积:t2(0,)WN441E 2()fd+则对任何正整数h,当 时,有以下结果 1 依分布收敛到 2 依分布收敛到 N 0101(,)hhN01(,).h1212(,)hhN12(,).hR RR+例2.1(
10、接第三章例1.1)对MA(q)序列 。利用定理2.2得到,只要当 依分布收敛到 的分布。注意 时,中的 应属于 ,所以令 有 tX:mmqNmR1(2),1mtmt mtmttRW mq 1mq0,0,mtmt mtm,q qltmqmllmlqRW+为期望为0,方差为 的正态分布。+在假设 是MA(q)下,对mq有mR221122q0:tHX221|P r(1.96)0.05122mqN+现在用 表示第三章例1.1中差分后的化学浓度数据。在 是MA(q)下。用 代替真值 后分别对 计算出tX0:tHXkk0,1q 22212(),1,6.1222mqqqNTmm 在q=0的假设下,所以应当否
11、定q=0.12345605.7780.2810.9510.1211.0710.11610.2430.8210.1040.9250.1001.631mqq(1)5.7781.96.aT+实际工作中人们还计算概率 并且把p称为检验的p值。明显p值越小,数据提供的否定原假设的依据越充分。现在在 下 ,近似服从标准正态分布。所以p值几乎是零,因而必须拒绝 是MA(0)的假设。取q=1时,所以不能拒绝 是MA(1)的假设。1(|5.778|).pPN 0H1NtX1|()|1.96(16).T mmtX+对于实际工作者来讲谱密度平方可积的条件通常很难验证。于是希望能把定理2.2中谱密度平方可积的条件改加
12、在自协方差函数 的收敛速度上。+定理2.3 对于一平稳序列 它的自协方差函数平方可积的充分必要条件是它的谱密度平方可积。k.tX+这个结论主要是利用实变函数论中Fourier级数的理论。只有证明 时用了周期图(如P.67定理3.1的证明,那里 绝对可和)。证明略。+推论2.4 设 是独立同分布的白噪声 满足 如果线性平稳序列(2.8)的自协方差函数平方可和:则定理2.2中的结论成立。()0fk t2(0,).WN44.tE 2.kk+定理2.2 要求白噪声的方差有4阶矩。下面关于线性平稳序列的样本自相关系数的中心极限定理不要求噪声项的4阶矩有限。+定理2.5 设 是独立同分布的 线性平稳序列
13、由(2.8)定义。如果自协方差函数 平方可和,并且对某个常数 (2.13)k t2(0,),WNtXk0.5,2|0,.kkmmm+则对任何正数h.当 时,依分布收敛到 ARMA序列的 满足(2.13).ARMA序列的白噪声列是独立同分布序列时定理2.5结论成立。N 1212(,)hhN12(,).hR RRj+推论2.6 如果 是独立同分布的白噪声,是样本自相关系数,则对任何正整数h:1:依分布收敛到多元标准正态分布 这里 是 的单位矩阵。tX121()()()NkNNttktkNNttxxxxxx12(,)hN(0,).hNIhIhh+2:如果 则 依分布收敛到44.tE 201(,)hN
14、2001(,).hM W WW+对白噪声,定理2.5的条件满足。第二条满足推论2.4的条件。20120001200000000001(2),1()()()()jtjtjtjtjtjjtjtjtjjtjtttRWWjMWWWWMWWMWM W+首先用图形表示N不同时 的误差。+然后重复M=1000次计算1000个 的标准差(称为标准误差)。发现N增大时标准误差减小。+误差随N减小的速度为 。+根离单位圆近的模型其估计标准误差大。kk12N+白噪声的白噪声的 检验检验+样本自相关置信区间检样本自相关置信区间检验法验法2+若 是独立同分布的白噪声,根据推论2.6,N足够大时 服从iid标准正态分布。
15、于是近似服从 分布。2tX12(,)mN222212()()mmN2()m+对于AR(2)模型取不同根离单位圆距离实验。根离单位圆越近与白噪声差别越大。+对AR(1)模型用不同的b模拟。B接近于1时与白噪声差别明显。+关于 中项数m的选取:m=5比m=20有效。注意以ARMA模型为例,当k较大时 已经很小,所以 贡献不大,取太大的m容易使检验不敏感。2()mk2k+白噪声的 检验法:是独立白噪声;是相关序列。下,拒绝域为 其中 20:tHX1:tHX0.05,5m22()().mm2222221250.050.05(5)(),(5)11.071,(20)31.41N+当 为独立同分布白噪声时 近似m维标准正态分布。+如果超过5%的 可否定 为独立同分布白噪声。+与 检验理由类似,m不应取太大。tX12(,)mNPr(|1.96)0.05.kN|1.96/jn0:tHX2+正态分布检验法:1()#|1.96,1.JQ mjNjmm+对 的检验可以比较成功。+对MA(1)的检验如果取m=20则很可能不成功。因为一般只有 超过界限。+对AR的检验一般成功。因为其相关系数不截尾。+演示(0,1)WN1