平稳时间序列分析王燕.ppt
- 文档编号:2453604
- 上传时间:2022-10-29
- 格式:PPT
- 页数:153
- 大小:1.95MB
平稳时间序列分析王燕.ppt
《平稳时间序列分析王燕.ppt》由会员分享,可在线阅读,更多相关《平稳时间序列分析王燕.ppt(153页珍藏版)》请在冰豆网上搜索。
平稳时间序列分析平稳时间序列分析:
ARMA模型模型西北大学经济管理学院田洪志本章结构方法性工具ARMA模型平稳序列建模序列预测3.1方法性工具差分运算延迟算子线性差分方程差分运算一阶差分阶差分步差分延迟算子延迟算子类似于一个时间指针,当前序列值乘以一个延迟算子,就相当于把当前序列值的时间向过去拨了一个时刻记B为延迟算子,有延迟算子的性质,其中用延迟算子表示差分运算阶差分步差分线性差分方程线性差分方程齐次线性差分方程齐次线性差分方程的解特征方程特征方程的根称为特征根,记作齐次线性差分方程的通解不相等实数根场合有相等实根场合复根场合非齐次线性差分方程的解非齐次线性差分方程的特解使得非齐次线性差分方程成立的任意一个解非齐次线性差分方程的通解齐次线性差分方程的通解和非齐次线性差分方程的特解之和3.2ARMA模型的性质AR模型(AutoRegressionModel)MA模型(MovingAverageModel)ARMA模型(AutoRegressionMovingAveragemodel)AR模型的定义具有如下结构的模型称为阶自回归模型,简记为特别当时,称为中心化模型AR(P)序列中心化变换称为的中心化序列,令自回归系数多项式引进延迟算子,中心化模型又可以简记为自回归系数多项式AR模型平稳性判别判别原因AR模型是常用的平稳序列的拟合模型之一,但并非所有的AR模型都是平稳的判别方法单位根判别法平稳域判别法例3.1:
考察如下四个模型的平稳性例3.1平稳序列时序图例3.1非平稳序列时序图AR模型平稳性判别方法特征根判别AR(p)模型平稳的充要条件是它的p个特征根都在单位圆内根据特征根和自回归系数多项式的根成倒数的性质,等价判别条件是该模型的自回归系数多项式的根都在单位圆外平稳域判别平稳域AR
(1)模型平稳条件特征根平稳域AR
(2)模型平稳条件特征根平稳域例3.1平稳性判别模型特征根判别平稳域判别结论
(1)平稳
(2)非平稳(3)平稳(4)非平稳平稳AR模型的统计性质均值方差协方差自相关系数偏自相关系数均值如果AR(p)模型满足平稳性条件,则有根据平稳序列均值为常数,且为白噪声序列,有推导出Green函数定义AR模型的传递形式其中系数称为Green函数Green函数递推公式原理方法待定系数法递推公式方差平稳AR模型的传递形式两边求方差得例3.2:
求平稳AR
(1)模型的方差平稳AR
(1)模型的传递形式为Green函数为平稳AR
(1)模型的方差协方差函数在平稳AR(p)模型两边同乘,再求期望根据得协方差函数的递推公式例3.3:
求平稳AR
(1)模型的协方差递推公式平稳AR
(1)模型的方差为协方差函数的递推公式为例3.4:
求平稳AR
(2)模型的协方差平稳AR
(2)模型的协方差函数递推公式为自相关系数自相关系数的定义平稳AR(P)模型的自相关系数递推公式常用AR模型自相关系数递推公式AR
(1)模型AR
(2)模型AR模型自相关系数的性质拖尾性呈复指数衰减例3.5:
考察如下AR模型的自相关图例3.5自相关系数按复指数单调收敛到零例3.5:
例3.5:
自相关系数呈现出“伪周期”性例3.5:
自相关系数不规则衰减偏自相关系数定义对于平稳AR(p)序列,所谓滞后k偏自相关系数就是指在给定中间k-1个随机变量的条件下,或者说,在剔除了中间k-1个随机变量的干扰之后,对影响的相关度量。
用数学语言描述就是偏自相关系数的计算滞后k偏自相关系数实际上就等于k阶自回归模型第个k回归系数的值。
偏自相关系数的截尾性AR(p)模型偏自相关系数P阶截尾例3.5续:
考察如下AR模型的偏自相关图例3.5理论偏自相关系数样本偏自相关图例3.5:
理论偏自相关系数样本偏自相关图例3.5:
理论偏自相关系数样本偏自相关图例3.5:
理论偏自相关系数样本偏自相关系数图MA模型的定义具有如下结构的模型称为阶自回归模型,简记为特别当时,称为中心化模型移动平均系数多项式引进延迟算子,中心化模型又可以简记为阶移动平均系数多项式MA模型的统计性质常数均值常数方差MA模型的统计性质自协方差函数P阶截尾自相关系数P阶截尾常用MA模型的自相关系数MA
(1)模型MA
(2)模型MA模型的统计性质偏自相关系数拖尾例3.6:
考察如下MA模型的相关性质MA模型的自相关系数截尾MA模型的自相关系数截尾MA模型的偏自相关系数拖尾MA模型的偏自相关系数拖尾MA模型的可逆性MA模型自相关系数的不唯一性例3.6中不同的MA模型具有完全相同的自相关系数和偏自相关系数可逆的定义可逆MA模型定义若一个MA模型能够表示称为收敛的AR模型形式,那么该MA模型称为可逆MA模型可逆概念的重要性一个自相关系数列唯一对应一个可逆MA模型。
可逆MA
(1)模型MA模型的可逆条件MA(q)模型的可逆条件是:
MA(q)模型的特征根都在单位圆内等价条件是移动平滑系数多项式的根都在单位圆外逆函数的递推公式原理方法待定系数法递推公式例3.6续:
考察如下MA模型的可逆性
(1)
(2)逆函数逆转形式(3)(4)逆函数逆转形式ARMA模型的定义具有如下结构的模型称为自回归移动平均模型,简记为特别当时,称为中心化模型系数多项式引进延迟算子,中心化模型又可以简记为阶自回归系数多项式阶移动平均系数多项式平稳条件与可逆条件ARMA(p,q)模型的平稳条件P阶自回归系数多项式的根都在单位圆外即ARMA(p,q)模型的平稳性完全由其自回归部分的平稳性决定ARMA(p,q)模型的可逆条件q阶移动平均系数多项式的根都在单位圆外即ARMA(p,q)模型的可逆性完全由其移动平滑部分的可逆性决定传递形式与逆转形式传递形式逆转形式ARMA(p,q)模型的统计性质均值协方差自相关系数ARMA模型的相关性自相关系数拖尾偏自相关系数拖尾例3.7:
考察ARMA模型的相关性拟合模型ARMA(1,1):
并直观地考察该模型自相关系数和偏自相关系数的性质。
自相关系数和偏自相关系数拖尾性样本自相关图样本偏自相关图ARMA模型相关性特征模型自相关系数偏自相关系数AR(P)拖尾P阶截尾MA(q)q阶截尾拖尾ARMA(p,q)拖尾拖尾3.3平稳序列建模建模步骤模型识别参数估计模型检验模型优化序列预测建模步骤平平稳稳非非白白噪噪声声序序列列计计算算样样本本相相关关系系数数模型模型识别识别参数参数估计估计模型模型检验检验模模型型优优化化序序列列预预测测YN计算样本相关系数样本自相关系数样本偏自相关系数模型识别基本原则选择模型拖尾P阶截尾AR(P)q阶截尾拖尾MA(q)拖尾拖尾ARMA(p,q)模型定阶的困难因为由于样本的随机性,样本的相关系数不会呈现出理论截尾的完美情况,本应截尾的或仍会呈现出小值振荡的情况由于平稳时间序列通常都具有短期相关性,随着延迟阶数,与都会衰减至零值附近作小值波动?
当或在延迟若干阶之后衰减为小值波动时,什么情况下该看作为相关系数截尾,什么情况下该看作为相关系数在延迟若干阶之后正常衰减到零值附近作拖尾波动呢?
样本相关系数的近似分布BarlettQuenouille模型定阶经验方法95的置信区间模型定阶的经验方法如果样本(偏)自相关系数在最初的d阶明显大于两倍标准差范围,而后几乎95的自相关系数都落在2倍标准差的范围以内,而且通常由非零自相关系数衰减为小值波动的过程非常突然。
这时,通常视为(偏)自相关系数截尾。
截尾阶数为d。
例2.5续选择合适的模型ARMA拟合1950年1998年北京市城乡居民定期储蓄比例序列。
序列自相关图序列偏自相关图拟合模型识别自相关图显示延迟3阶之后,自相关系数全部衰减到2倍标准差范围内波动,这表明序列明显地短期相关。
但序列由显著非零的相关系数衰减为小值波动的过程相当连续,相当缓慢,该自相关系数可视为不截尾偏自相关图显示除了延迟1阶的偏自相关系数显著大于2倍标准差之外,其它的偏自相关系数都在2倍标准差范围内作小值随机波动,而且由非零相关系数衰减为小值波动的过程非常突然,所以该偏自相关系数可视为一阶截尾所以可以考虑拟合模型为AR
(1)例3.8美国科罗拉多州某一加油站连续57天的OVERSHORT序列序列自相关图序列偏自相关图拟合模型识别自相关图显示除了延迟1阶的自相关系数在2倍标准差范围之外,其它阶数的自相关系数都在2倍标准差范围内波动。
根据这个特点可以判断该序列具有短期相关性,进一步确定序列平稳。
同时,可以认为该序列自相关系数1阶截尾偏自相关系数显示出典型非截尾的性质。
综合该序列自相关系数和偏自相关系数的性质,为拟合模型定阶为MA
(1)例3.91880-1985全球气表平均温度改变值差分序列序列自相关图序列偏自相关图拟合模型识别自相关系数显示出不截尾的性质偏自相关系数也显示出不截尾的性质综合该序列自相关系数和偏自相关系数的性质,可以尝试使用ARMA(1,1?
)模型拟合该序列参数估计待估参数个未知参数常用估计方法矩估计极大似然估计最小二乘估计矩估计原理样本自相关系数估计总体自相关系数样本一阶均值估计总体均值,样本方差估计总体方差例3.10:
求AR
(2)模型系数的矩估计AR
(2)模型Yule-Walker方程矩估计(Yule-Walker方程的解)例3.11:
求MA
(1)模型系数的矩估计MA
(1)模型方程矩估计例3.12:
求ARMA(1,1)模型系数的矩估计ARMA(1,1)模型方程矩估计对矩估计的评价优点估计思想简单直观不需要假设总体分布计算量小(低阶模型场合)缺点信息浪费严重只用到了p+q个样本自相关系数信息,其他信息都被忽略估计精度差通常矩估计方法被用作极大似然估计和最小二乘估计迭代计算的初始值极大似然估计原理在极大似然准则下,认为样本来自使该样本出现概率最大的总体。
因此未知参数的极大似然估计就是使得似然函数(即联合密度函数)达到最大的参数值似然方程由于和都不是的显式表达式。
因而似然方程组实际上是由p+q+1个超越方程构成,通常需要经过复杂的迭代算法才能求出未知参数的极大似然估计值对极大似然估计的评价优点极大似然估计充分应用了每一个观察值所提供的信息,因而它的估计精度高同时还具有估计的一致性、渐近正态性和渐近有效性等许多优良的统计性质缺点需要假定总体分布最小二乘估计原理使残差平方和达到最小的那组参数值即为最小二乘估计值条件最小二乘估计实际中最常用的参数估计方法假设条件残差平方和方程解法迭代法对最小二乘估计的评价优点最小二乘估计充分应用了每一个观察值所提供的信息,因而它的估计精度高条件最小二乘估计方法使用率最高缺点需要假定总体分布例2.5续确定1950年1998年北京市城乡居民定期储蓄比例序列拟合模型的口径拟合模型:
AR
(1)估计方法:
极大似然估计模型口径例3.8续确定美国科罗拉多州某一加油站连续57天的OVERSHORTS序列拟合模型的口径拟合模型:
MA
(1)估计方法:
条件最小二乘估计模型口径例3.9续确定1880-1985全球气表平均温度改变值差分序列拟合模型的口径拟合模型:
ARMA(1,1)估计方法:
条件最小二乘估计模型口径模型检验模型的显著性检验整个模型对信息的提取是否充分参数的显著性检验模型结构是否最简模型的显著性检验目的检验模型的有效性(对信息的提取是否充分)检验对象残差序列判定原则一个好的拟合模型应该能够提取观察值序列中几乎所有的样本相关信息,即残差序列应该为白噪声序列反之,如果残差序列为非白噪声序列,那就意味着残差序列中还残留着相关信息未被提取,这就说明拟合模型不够有效假设条件原假设:
残差序列为白噪声序列备择假设:
残差序列为非白噪声序列检验统计量LB统计量例2.5续检验1950年1998年北京市城乡居民定期储蓄比例序列拟合模型的显著性残差白噪声序列检验结果延迟阶数LB统计量P值检验结论65.830.3229拟合模型显著有效1210.280.50501811.380.8361参数显著性检验目的检验每一个未知参数是否显著非零。
删除不显著参数使模型结构最精简假设条件检验统计量例2.5续检验195
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 平稳 时间 序列 分析
![提示](https://static.bdocx.com/images/bang_tan.gif)