《FINTS第四章线性ARMA模型.ppt》由会员分享，可在线阅读，更多相关《FINTS第四章线性ARMA模型.ppt（28页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、模型定阶或识别n假设数据已经平稳化，下一步是确定模型的阶数。有两种方法，一种是根据随机过程的参数特征，一种是根据信息准则。n下面是几类随机过程的参数特征：三种随机过程偏自相关函数的特点三类过程的偏自相关函数和自相关函数MA（q）AR（p）ARMA（p，q）自相关函数q步截尾拖尾拖尾偏自相关函数拖尾p步截尾拖尾定阶样本自相关函数的计算和判断定阶H0：i=i+1=0用前面介绍的方法计算出样本自相关系数，零假设成立时近似服从正态分布N(0,1/T)所以近似5%显著水平下，每个 在两倍标准差之间，则不能拒绝零假设。具体地说，如果-2/p那么估计的偏相关系数近似服从正态分布N（0，1/T）所以近似5%显
2、著水平下，如果-2/T1/2*kp成立n评价模型的优劣准则根据信息准则进行模型识别(定阶)残差平方和最小化AIC(Akaikes information criterion)BIC(Schwartz Bayesian information criterion)准则n对自由度进行调整nk是模型中未知参数的个数,et是估计出的误差定阶：AIC准则和准则和BIC准则准则不同的书对AIC和BIC使用不同的变形。经常使用的有两种AIC（p，q）=ln()+2(p+q)/TBIC（p，q）=ln()+(p+q)ln(T)/TT样本长度，如果有常数项p+q被p+q+1代替，ln表示自然对数。在ARMA模型
3、中需要选择p和q，所以用p+q代替k。是对噪声项方差的估计定阶：AIC准则和准则和BIC准则准则AIC（p，q）=2lnL/T+2(p+q)/TBIC（p，q）=-2lnL/T+(p+q)ln(T)/TLnL是模型的对数似然函数值Q是与参数无关的量。因为我们只关心使得AIC或BIC最小的值，所以忽略Q.带入对数似然函数表达式中，可以发现与前面的AIC和BIC的表达是一致的。AIC和BIC判断步骤（1）给定滞后长度的上限P和Q，一般取为T/10，Ln(T),（2）修改样本区间使得滞后长度不出现负值。（3）对任意一对滞后长度p=0，1，P，q=0，1，Q，分别估计模型ARMA（p，q）（4）代入上
4、面的公式，计算出AIC(p，q)和BIC(p，q)（5）最小值对应的p，q值作为ARMA模型的阶数。用AIC和BIC准则确定阶数AIC准则-MA(1)q0123P0-7.4152-7.433-7.383-7.174-7.221用AIC和BIC准则确定阶数BIC-白噪声q0123P0-7.4152-7.345-7.251-6.998-7.001练习：nP17915(9)极大似然估计:以AR(1)为例t=c+t-1+t假设 i.i.d.N(0,2)估计:=(c,2)已知:y1,y2,yTE(1)=c/(1-)E(1-)2=2/(1-2)极大似然估计当1的观测已知时，2的条件分布2=c+1+2（2|
5、1=y1）N(c+y1,2)极大似然估计Y1,Y2的联合分布密度函数，是条件密度和边际密度相乘f2，Y1(y2，y1;)=f2|Y1(y2|y1;)f1(y1;)类似的，已知y1,y2，3的条件分布极大似然估计三者的联合分布f3，2，Y1(y3，y2，y1;)=f3|Y2，Y1(y3|y2，y1;)f2|Y1(y2|y1;)f1(y1;)一般给定y1,y2,yt-1，t的条件分布只和yt-1有关极大似然估计ft,Yt-1，,Y1(yt,yt-1，,y1;)=f1(y1;)ft|Yt-1(yt|yt-1;)（）估计：满足下面的条件的解求解未知参数的方程是非线性的，如果只关心（2，T）的条件联合分
6、布，得到条件极大似然函数。极大似然估计同样通过解方程来得到未知参数的估计：这时得到的是线性方程组最小二乘估计法：计算例子-produceastationaryAR(2)process:yii-1i-2+xiandfindtheestimateofparametersMatlabcode:nx=randn(1000,1);ny=filter(1,1-0.60.08,x);%产生上述AR(2)过程nm=ar(y,2)模型的检验检验残差是否是白噪声过程1）画出残差的折线图2）画出残差的ACF,PACF3)计算统计量QBox-PierceQ-检验其中，T为样本容量Ljung and Box检验Q检验1
7、）m主观给定，一般在15到30之间，可令m=T1/22）H0：t是白噪声过程3)H0成立时，统计量Q渐进服从2（m-p-q)，如果 模型中包括常数项，那么Q渐进服从2(m-1-p-q)4）Q检验会给出相应的P-值（P-值拒绝H0 0）Q检验图示真实临界值计算值卡方分布临界Q检验存在缺陷:经常不能拒绝零假设。把不是白噪声时，也误认为是白噪声检验练习例m=6，模型中有常数项，考虑下面的几个模型，哪个模型是合格的模型？给出其它几个模型Q检验统计量的自由度。(p+q)Q自由度P-value模型选择一个好模型满足的条件n每个解释变量都显著不等于0.n残差是白噪声过程n具有最小的AIC或BIC值练习：从下面的几个模型中选择一个最优模型AR(1)AR(2)AR(3)ARMA(1,1)MA(2)1 0.17 0.21 0.3 0.19 (0.0000)(0.0004)(0.002)(0.0024)2 0.06 0.04 (0.0005)(0.003)3 (0.44)1 (0.0007)(0.0034)2 (0.009)Q(8)P-值值 Q(16)P-值值 1.对AR(2)模型其中求的条件极大似然函数并写出应当满足的方程。作业的极大似然估计2.