《非平稳序列随机分析PPT课件.ppt》由会员分享，可在线阅读，更多相关《非平稳序列随机分析PPT课件.ppt（118页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、关于非平稳序列的随机分析第一张，PPT共一百一十八页，创作于2022年6月本章结构n差分运算nARIMA模型nAuto-Regressive模型n异方差的性质n方差齐性变化n条件异方差模型第二张，PPT共一百一十八页，创作于2022年6月5.1 差分运算n差分运算的实质n差分方式的选择n过差分第三张，PPT共一百一十八页，创作于2022年6月差分运算的实质n差分方法是一种非常简便、有效的确定性信息提取方法nCramer分解定理在理论上保证了适当阶数的差分一定可以充分提取确定性信息n差分运算的实质是使用自回归的方式提取确定性信息 第四张，PPT共一百一十八页，创作于2022年6月差分方式的选择n
2、序列蕴含着显著的线性趋势，一阶差分就可以实现趋势平稳 n序列蕴含着曲线趋势，通常低阶（二阶或三阶）差分就可以提取出曲线趋势的影响 n对于蕴含着固定周期的序列进行步长为周期长度的差分运算，通常可以较好地提取周期信息 第五张，PPT共一百一十八页，创作于2022年6月例5.1【例1.1】1964年1999年中国纱年产量序列蕴含着一个近似线性的递增趋势。对该序列进行一阶差分运算 考察差分运算对该序列线性趋势信息的提取作用 第六张，PPT共一百一十八页，创作于2022年6月差分前后时序图n原序列时序图n差分后序列时序图第七张，PPT共一百一十八页，创作于2022年6月例5.2n尝试提取1950年199
3、9年北京市民用车辆拥有量序列的确定性信息第八张，PPT共一百一十八页，创作于2022年6月差分后序列时序图n一阶差分n二阶差分第九张，PPT共一百一十八页，创作于2022年6月例5.3n差分运算提取1962年1月1975年12月平均每头奶牛的月产奶量序列中的确定性信息 第十张，PPT共一百一十八页，创作于2022年6月差分后序列时序图n一阶差分n1阶12步差分第十一张，PPT共一百一十八页，创作于2022年6月过差分 n足够多次的差分运算可以充分地提取原序列中的非平稳确定性信息n但过度的差分会造成有用信息的浪费 第十二张，PPT共一百一十八页，创作于2022年6月例5.4n假设序列如下 n考察
4、一阶差分后序列和二阶差分序列 的平稳性与方差 第十三张，PPT共一百一十八页，创作于2022年6月比较n一阶差分n平稳n方差小n二阶差分（过差分）n平稳n方差大第十四张，PPT共一百一十八页，创作于2022年6月5.2 ARIMA模型nARIMA模型结构nARIMA模型性质nARIMA模型建模nARIMA模型预测n疏系数模型n季节模型第十五张，PPT共一百一十八页，创作于2022年6月ARIMA模型结构n使用场合n差分平稳序列拟合n模型结构第十六张，PPT共一百一十八页，创作于2022年6月ARIMA 模型族nd=0ARIMA(p,d,q)=ARMA(p,q)nP=0ARIMA(P,d,q)=
5、IMA(d,q)nq=0ARIMA(P,d,q)=ARI(p,d)nd=1,P=q=0ARIMA(P,d,q)=random walk model第十七张，PPT共一百一十八页，创作于2022年6月随机游走模型(random walk)n模型结构n模型产生典故nKarl Pearson(1905)在自然杂志上提问：假如有个醉汉醉得非常严重，完全丧失方向感，把他放在荒郊野外，一段时间之后再去找他，在什么地方找到他的概率最大呢？第十八张，PPT共一百一十八页，创作于2022年6月ARIMA模型的平稳性nARIMA(p,d,q)模 型共有p+d个特征根，其中p个在单位圆内，d个在单位圆上。所以当 时
6、ARIMA(p,d,q)模 型非平稳。n例5.5ARIMA(0,1,0)时序图第十九张，PPT共一百一十八页，创作于2022年6月ARIMA模型的方差齐性n 时，原序列方差非齐性nd阶差分后，差分后序列方差齐性第二十张，PPT共一百一十八页，创作于2022年6月ARIMA模型建模步骤获获得得观观察察值值序序列列平稳性平稳性检验检验差分差分运算运算YN白噪声白噪声检验检验Y分分析析结结束束N拟合拟合ARMA模型模型第二十一张，PPT共一百一十八页，创作于2022年6月例5.6n对1952年1988年中国农业实际国民收入指数序列建模 第二十二张，PPT共一百一十八页，创作于2022年6月一阶差分序
7、列时序图第二十三张，PPT共一百一十八页，创作于2022年6月一阶差分序列自相关图第二十四张，PPT共一百一十八页，创作于2022年6月一阶差分后序列白噪声检验延迟阶数 统计量P值615.330.01781218.330.10601824.660.1344第二十五张，PPT共一百一十八页，创作于2022年6月拟合ARMA模型n偏自相关图第二十六张，PPT共一百一十八页，创作于2022年6月建模n定阶nARIMA(0,1,1)n参数估计n模型检验n模型显著n参数显著第二十七张，PPT共一百一十八页，创作于2022年6月ARIMA模型预测n原则n最小均方误差预测原理 nGreen函数递推公式第二十
8、八张，PPT共一百一十八页，创作于2022年6月预测值第二十九张，PPT共一百一十八页，创作于2022年6月例5.7n已知ARIMA(1,1,1)模型为 且n求 的95的置信区间 第三十张，PPT共一百一十八页，创作于2022年6月预测值n等价形式n计算预测值第三十一张，PPT共一百一十八页，创作于2022年6月计算置信区间nGreen函数值n方差n95置信区间第三十二张，PPT共一百一十八页，创作于2022年6月例5.6续：对中国农业实际国民收入指数序列做为期10年的预测 第三十三张，PPT共一百一十八页，创作于2022年6月疏系数模型nARIMA(p,d,q)模型是指d阶差分后自相关最高阶
9、数为p，移动平均最高阶数为q的模型，通常它包含p+q个独立的未知系数：n如果该模型中有部分自相关系数 或部分移动平滑系数 为零，即原模型中有部分系数省缺了，那么该模型称为疏系数模型。第三十四张，PPT共一百一十八页，创作于2022年6月疏系数模型类型n如果只是自相关部分有省缺系数，那么该疏系数模型可以简记为n 为非零自相关系数的阶数n如果只是移动平滑部分有省缺系数，那么该疏系数模型可以简记为n 为非零移动平均系数的阶数n如果自相关和移动平滑部分都有省缺，可以简记为第三十五张，PPT共一百一十八页，创作于2022年6月例5.8n对1917年1975年美国23岁妇女每万人生育率序列建模 第三十六张
10、，PPT共一百一十八页，创作于2022年6月一阶差分第三十七张，PPT共一百一十八页，创作于2022年6月自相关图第三十八张，PPT共一百一十八页，创作于2022年6月偏自相关图第三十九张，PPT共一百一十八页，创作于2022年6月建模n定阶nARIMA(1,4),1,0)n参数估计n模型检验n模型显著n参数显著第四十张，PPT共一百一十八页，创作于2022年6月季节模型n简单季节模型n乘积季节模型 第四十一张，PPT共一百一十八页，创作于2022年6月简单季节模型n简单季节模型是指序列中的季节效应和其它效应之间是加法关系n简单季节模型通过简单的趋势差分、季节差分之后序列即可转化为平稳，它的模
11、型结构通常如下 第四十二张，PPT共一百一十八页，创作于2022年6月例5.9n拟合19621991年德国工人季度失业率序列 第四十三张，PPT共一百一十八页，创作于2022年6月差分平稳n对原序列作一阶差分消除趋势，再作4步差分消除季节效应的影响，差分后序列的时序图如下 第四十四张，PPT共一百一十八页，创作于2022年6月白噪声检验延迟阶数 统计量P值643.840.00011251.710.00011854.480.0001第四十五张，PPT共一百一十八页，创作于2022年6月差分后序列自相关图第四十六张，PPT共一百一十八页，创作于2022年6月差分后序列偏自相关图第四十七张，PPT共
12、一百一十八页，创作于2022年6月模型拟合n定阶nARIMA(1,4),(1,4),0)n参数估计第四十八张，PPT共一百一十八页，创作于2022年6月模型检验残差白噪声检验参数显著性检验延迟阶数 统计量P值待估参数 统计量P值62.090.71915.480.00011210.990.3584-3.410.0001第四十九张，PPT共一百一十八页，创作于2022年6月拟合效果图第五十张，PPT共一百一十八页，创作于2022年6月乘积季节模型n使用场合n序列的季节效应、长期趋势效应和随机波动之间有着复杂地相互关联性，简单的季节模型不能充分地提取其中的相关关系 n构造原理n短期相关性用低阶ARM
13、A(p,q)模型提取n季节相关性用以周期步长S为单位的ARMA(P,Q)模型提取n假设短期相关和季节效应之间具有乘积关系，模型结构如下 第五十一张，PPT共一百一十八页，创作于2022年6月例5.10:拟合19481981年美国女性月度失业率序列 第五十二张，PPT共一百一十八页，创作于2022年6月差分平稳n一阶、12步差分第五十三张，PPT共一百一十八页，创作于2022年6月差分后序列自相关图第五十四张，PPT共一百一十八页，创作于2022年6月差分后序列偏自相关图第五十五张，PPT共一百一十八页，创作于2022年6月简单季节模型拟合结果延迟阶数拟合模型残差白噪声检验AR(1,12)MA(
14、1,2,12)ARMA(1,12),(1,12)值P值 值P值 值P值614.580.00579.50.023315.770.00041216.420.088314.190.115817.990.0213结果拟合模型均不显著第五十六张，PPT共一百一十八页，创作于2022年6月乘积季节模型拟合n模型定阶nARIMA(1,1,1)(0,1,1)12n参数估计第五十七张，PPT共一百一十八页，创作于2022年6月模型检验残差白噪声检验参数显著性检验延迟阶数 统计量P值待估参数 统计量P值64.500.2120-4.660.0001129.420.400223.030.00011820.580.15
15、07-6.810.0001结果模型显著参数均显著第五十八张，PPT共一百一十八页，创作于2022年6月乘积季节模型拟合效果图第五十九张，PPT共一百一十八页，创作于2022年6月5.3 Auto-Regressive模型n构造思想n首先通过确定性因素分解方法提取序列中主要的确定性信息n然后对残差序列拟合自回归模型，以便充分提取相关信息 第六十张，PPT共一百一十八页，创作于2022年6月Auto-Regressive模型结构第六十一张，PPT共一百一十八页，创作于2022年6月对趋势效应的常用拟合方法n自变量为时间t的幂函数n自变量为历史观察值第六十二张，PPT共一百一十八页，创作于2022年
16、6月对季节效应的常用拟合方法n给定季节指数n建立季节自回归模型第六十三张，PPT共一百一十八页，创作于2022年6月例5.6续n使用Auto-Regressive模型分析1952年1988年中国农业实际国民收入指数序列。n时序图显示该序列有显著的线性递增趋势，但没有季节效应，所以考虑建立如下结构的Auto-Regressive模型 第六十四张，PPT共一百一十八页，创作于2022年6月趋势拟合n方法一：变量为时间t的幂函数n方法二：变量为一阶延迟序列值 第六十五张，PPT共一百一十八页，创作于2022年6月趋势拟合效果图第六十六张，PPT共一百一十八页，创作于2022年6月残差自相关检验n检验
17、原理n回归模型拟合充分，残差的性质n回归模型拟合得不充分，残差的性质第六十七张，PPT共一百一十八页，创作于2022年6月Durbin-Waston检验（DW检验）n假设条件n原假设：残差序列不存在一阶自相关性 n备择假设：残差序列存在一阶自相关性 第六十八张，PPT共一百一十八页，创作于2022年6月DW统计量n构造统计量nDW统计量和自相关系数的关系第六十九张，PPT共一百一十八页，创作于2022年6月DW统计量的判定结果正相关相关性待定不相关相关性待定负相关042第七十张，PPT共一百一十八页，创作于2022年6月例5.6续 n检验第一个确定性趋势模型 残差序列的自相关性。第七十一张，P
18、PT共一百一十八页，创作于2022年6月DW检验结果n检验结果n检验结论n检验结果显示残差序列高度正自相关。DW统计量的值P值0.13781.421.530.0001第七十二张，PPT共一百一十八页，创作于2022年6月Durbin h检验 nDW统计量的缺陷n当回归因子包含延迟因变量时，残差序列的DW统计量是一个有偏统计量。在这种场合下使用DW统计量容易产生残差序列正自相关性不显著的误判 nDurbin h检验第七十三张，PPT共一百一十八页，创作于2022年6月例5.6续n检验第二个确定性趋势模型 残差序列的自相关性。第七十四张，PPT共一百一十八页，创作于2022年6月Dh检验结果n检验
19、结果n检验结论n检验结果显示残差序列高度正自相关。Dh统计量的值P值2.80380.0025第七十五张，PPT共一百一十八页，创作于2022年6月残差序列拟合n确定自回归模型的阶数n参数估计n模型检验第七十六张，PPT共一百一十八页，创作于2022年6月例5.6续n对第一个确定性趋势模型的残差序列 进行拟合第七十七张，PPT共一百一十八页，创作于2022年6月残差序列自相关图第七十八张，PPT共一百一十八页，创作于2022年6月残差序列偏自相关图第七十九张，PPT共一百一十八页，创作于2022年6月模型拟合n定阶nAR(2)n参数估计方法n极大似然估计n最终拟合模型口径第八十张，PPT共一百一
20、十八页，创作于2022年6月例5.6n第二个AutoRegressive模型的拟合结果第八十一张，PPT共一百一十八页，创作于2022年6月三个拟合模型的比较模型AICSBCARIMA(0,1,1)模型：249.3305252.4976AutoRegressive模型一：260.8454267.2891AutoRegressive模型二：250.6317253.7987第八十二张，PPT共一百一十八页，创作于2022年6月5.4 异方差的性质n异方差的定义n如果随机误差序列的方差会随着时间的变化而变化，这种情况被称作为异方差n异方差的影响n忽视异方差的存在会导致残差的方差会被严重低估，继而参数
21、显著性检验容易犯纳伪错误，这使得参数的显著性检验失去意义，最终导致模型的拟合精度受影响。第八十三张，PPT共一百一十八页，创作于2022年6月异方差直观诊断n残差图n残差平方图第八十四张，PPT共一百一十八页，创作于2022年6月残差图n方差齐性残差图n递增型异方差残差图第八十五张，PPT共一百一十八页，创作于2022年6月残差平方图n原理n残差序列的方差实际上就是它平方的期望。n所以考察残差序列是否方差齐性，主要是考察残差平方序列是否平稳 第八十六张，PPT共一百一十八页，创作于2022年6月例5.11n直观考察美国1963年4月1971年7月短期国库券的月度收益率序列的方差齐性。第八十七张
22、，PPT共一百一十八页，创作于2022年6月一阶差分后残差图第八十八张，PPT共一百一十八页，创作于2022年6月一阶差分后残差平方图第八十九张，PPT共一百一十八页，创作于2022年6月异方差处理方法n假如已知异方差函数具体形式，进行方差齐性变化n假如不知异方差函数的具体形式，拟合条件异方差模型 第九十张，PPT共一百一十八页，创作于2022年6月5.5 方差齐性变换n使用场合n序列显示出显著的异方差性，且方差与均值之间具有某种函数关系 其中：是某个已知函数n处理思路n尝试寻找一个转换函数 ，使得经转换后的变量满足方差齐性第九十一张，PPT共一百一十八页，创作于2022年6月转换函数的确定原
23、理n转换函数 在 附近作一阶泰勒展开n求转换函数的方差n转换函数的确定第九十二张，PPT共一百一十八页，创作于2022年6月常用转换函数的确定n假定n转换函数的确定第九十三张，PPT共一百一十八页，创作于2022年6月例5.11续n对美国1963年4月1971年7月短期国库券的月度收益率序列使用方差齐性变换方法进行分析 n假定n函数变换第九十四张，PPT共一百一十八页，创作于2022年6月对数序列时序图第九十五张，PPT共一百一十八页，创作于2022年6月一阶差分后序列图第九十六张，PPT共一百一十八页，创作于2022年6月白噪声检验延迟阶数LB统计量P值63.580.73371210.820
24、.54411821.710.2452第九十七张，PPT共一百一十八页，创作于2022年6月拟合模型口径及拟合效果图第九十八张，PPT共一百一十八页，创作于2022年6月5.6 条件异方差模型nARCH模型nGARCH模型nGARCH模型的变体nEGARCH模型nIGARCH模型nGARCH-M模型nAR-GARCH模型第九十九张，PPT共一百一十八页，创作于2022年6月ARCH模型n假定n原理n通过构造残差平方序列的自回归模型来拟合异方差函数 nARCH(q)模型结构第一百张，PPT共一百一十八页，创作于2022年6月GARCH 模型结构n使用场合nARCH模型实际上适用于异方差函数短期自相
25、关过程 nGARCH模型实际上适用于异方差函数长期自相关过程 n模型结构第一百零一张，PPT共一百一十八页，创作于2022年6月GARCH模型的约束条件n参数非负 n参数有界 第一百零二张，PPT共一百一十八页，创作于2022年6月EGARCH模型第一百零三张，PPT共一百一十八页，创作于2022年6月IGARCH模型第一百零四张，PPT共一百一十八页，创作于2022年6月GARCH-M模型第一百零五张，PPT共一百一十八页，创作于2022年6月AR-GARCH模型第一百零六张，PPT共一百一十八页，创作于2022年6月GARCH模型拟合步骤n回归拟合n残差自相关性检验n异方差自相关性检验nA
26、RCH模型定阶n参数估计n正态性检验第一百零七张，PPT共一百一十八页，创作于2022年6月例5.12n使用条件异方差模型拟合某金融时间序列。第一百零八张，PPT共一百一十八页，创作于2022年6月回归拟合n拟合模型n参数估计n参数显著性检验nP值0.0001,参数高度显著 第一百零九张，PPT共一百一十八页，创作于2022年6月残差自相关性检验n残差序列DW检验结果nDurbin h=-2.6011n n拟合残差自回归模型n方法：逐步回归n模型口径第一百一十张，PPT共一百一十八页，创作于2022年6月异方差自相关检验nPortmantea Q检验n拉格朗日乘子（LM）检验 第一百一十一张，
27、PPT共一百一十八页，创作于2022年6月Portmantea Q检验n假设条件n检验统计量n检验结果n拒绝原假设n接受原假设第一百一十二张，PPT共一百一十八页，创作于2022年6月LM检验n假设条件n检验统计量n检验结果n拒绝原假设n接受原假设第一百一十三张，PPT共一百一十八页，创作于2022年6月例5.12残差序列异方差检验第一百一十四张，PPT共一百一十八页，创作于2022年6月ARCH模型拟合n定阶：GARCH(1,1)n参数估计：极大似然估计n拟合模型口径：AR(2)-GARCH(1,1)第一百一十五张，PPT共一百一十八页，创作于2022年6月模型检验n检验方法：正态性检验n假设条件：n检验统计量n检验结果n拒绝原假设n接受原假设第一百一十六张，PPT共一百一十八页，创作于2022年6月例5.13正态性检验结果 P值0.5603 nAR(2)-GARCH(1,1)模型显著成立第一百一十七张，PPT共一百一十八页，创作于2022年6月感谢大家观看第一百一十八张，PPT共一百一十八页，创作于2022年6月