当前位置:文档之家› 时间管理时间序列模型

时间管理时间序列模型

最新卓越管理方案您可自由编辑Wold 分解定理:任何协方差平稳过程x t ,都可以被表示为x t - μ - d t = u t + ψ1 u t -1+ ψ2 u t -2 + … + =∑∞=-0j jt j u ψ其中μ 表示x t 的期望。

d t 表示x t 的线性确定性成分,如周期性成分、时间t 的多项式和指数形式等,可以直接用x t 的滞后值预测。

ψ0 = 1,∑∞=02j j ψ< ∞。

u t 为白噪声过程。

u t 表示用x t 的滞后项预测x t 时的误差。

u t = x t - E(x t | x t -1, x t -2 , …)∑∞=-0j j t j u ψ称为x t 的线性非确定性成分。

当d t = 0时,称x t 为纯线性非确定性过程。

Wold 分解定理由Wold 在1938年提出。

Wold 分解定理只要求过程2阶平稳即可。

从原理上讲,要得到过程的Wold 分解,就必须知道无限个ψj 参数,这对于一个有限样本来说是不可能的。

实际中可以对ψj 做另一种假定,即可以把ψ (L )看作是2个有限特征多项式的比,ψ(L ) =∑∞=0j jj L ψ=)()(L L ΦΘ=pp qq L L L LL L φφφθθθ++++++++...1...1221221 注意,无论原序列中含有何种确定性成分,在前面介绍的模型种类中,还是后面介绍的自相关函数、偏自相关函数中都假设在原序列中已经剔除了所有确定性成分,是一个纯的随机过程(过程中不含有任何确定性成分)。

如果一个序列如上式, x t = μ + d t + u t + ψ1 u t -1+ ψ2 u t -2 + … + 则所有研究都是在y t = x t - μ - d t 的基础上进行。

例如前面给出的各类模型中都不含有均值项、时间趋势项就是这个道理。

2.3 自相关函数以上介绍了随机过程的几种模型。

实际中单凭对时间序列的观察很难确定其属于哪一种模型,而自相关函数和偏自相关函数是分析随机过程和识别模型的有力工具。

1. 自相关函数定义在给出自相关函数定义之前先介绍自协方差函数概念。

由第一节知随机过程{x t }中的每一个元素x t ,t = 1, 2, … 都是随机变量。

对于平稳的随机过程,其期望为常数,用 μ 表示,即E(x t ) = μ, t = 1, 2, … (2.25) 随机过程的取值将以 μ 为中心上下变动。

平稳随机过程的方差也是一个常量Var(x t ) = E [(x t - E(x t ))2 ] = E [(x t - μ)2 ] = σx 2 , t = 1, 2, … (2.26)σx 2用来度量随机过程取值对其均值 μ 的离散程度。

相隔k 期的两个随机变量x t 与x t - k 的协方差即滞后k 期的自协方差,定义为γk = Cov (x t , x t - k ) = E[(x t - μ ) (x t - k - μ ) ] (2.27) 自协方差序列γk , k = 0, 1, …, K ,称为随机过程 {x t } 的自协方差函数。

当k = 0 时γ0 = Var (x t ) = σx 2 自相关系数定义ρk =)()(),(k t t k t t x Var x ar V x x Cov -- (2.28)因为对于一个平稳过程有Var (x t ) = Var (x t - k ) = σx 2 (2.29) 所以(2.28)可以改写为 ρk =2),(x k t t x x Cov σ- =2x k σγ=γγk(2.30) 当 k = 0 时,有 ρ 0 = 1。

以滞后期k 为变量的自相关系数列ρk , k = 0, 1, …, K (2.31) 称为自相关函数。

因为ρk = ρ- k 即Cov (x t - k , x t ) = Cov (x t , x t + k ),自相关函数是零对称的,所以实际研究中只给出自相关函数的正半部分即可。

2.自回归过程的自相关函数(1) 平稳AR(1)过程的自相关函数 AR(1) 过程如下x t = φ1 x t -1 + u t , |φ1| < 1 用x t- k 同乘上式两侧x t x t- k = φ1 x t -1 x t- k + u t x t- k 两侧同取期望, γk = φ1 γk -1其中E(x t- k u t ) = 0(u t 与其t - k 期及以前各项都不相关)。

两侧同除 γ0 得, ρk = φ1 ρk -1 = φ1 φ1 ρk -2 = … = φ1k ρ0 因为 ρo = 1。

所以有 ρk = φ1k , (k ≥ 0)对于平稳序列有 | φ1| < 1。

所以当 φ1为正时,自相关函数按指数衰减至零(过阻尼情形),当 φ1为负时,自相关函数正负交错地指数衰减至零。

见图2.6。

因为对于经济时间序列,φ1一般为正,所以第一种情形常见。

指数衰减至零的表现形式说明随着时间间隔的加长,变量之间的关系变得越来越弱。

φ1 > 0 (经济问题中常见) φ1 < 0 (经济问题中少见)图2.6 AR(1) 过程的自相关函数(2)AR(p ) 过程的自相关函数用x t - k , (k > 0) 同乘平稳的 p 阶自回归过程x t = φ 1 x t -1 + φ 2 x t -2 +…+ φ p x t - p + u t (2.32) 的两侧,得x t - k x t = φ1 x t - k x t -1 + φ2 x t - k x t -2 + … + φp x t - k x t - p + x t - k u t (2.33) 对上式两侧分别求期望得γk = φ1 γk -1 + φ2 γk -2 + … + φp γk - p , k > 0 (2.34) 上式中对于 k > 0,有E(x t - k u t ) = 0。

因为当 k > 0时,x t - k 发生在u t 之前,所以 x t - k 与 u t不相关。

用 γ0分别除(2.34)式的两侧得ρk = φ1 ρk -1 + φ2 ρk -2 + … + φp ρk -p , k > 0 (2.35) 令 Φ(L ) = (1 - φ1 L - φ2 L 2 - … - φp L p )其中L 为k 的滞后算子,则上式可表达为 Φ(L ) ρk = 0 因 Φ(L ) 可因式分解为, Φ(L ) =∏=pi i L G 1) -(1,则(2.35)式的通解(证明见附录)是ρk = A 1 G 1k + A 2 G 2k + … + A p G p k . (2.36) 其中A i , i = 1, … p 为待定常数。

这里 G i -1, i = 1, 2, …, p 是特征方程 Φ(L ) = (1 - φ1 L - φ2 L 2 - … - φp L p ) = 0的根。

为保证随机过程的平稳性,要求 | G i | < 1, i = 1, 2, …, p 。

这会遇到如下两种情形。

① 当G i 为实数时,(2.36) 式中的A i G i k 将随着k 的增加而几何衰减至零,称为指数衰减(过阻尼情形)。

② 当G i 和G j 表示一对共轭复根时,设G i = a + bi , G j = a – bi,22b a += R ,则G i , G j的极座标形式是G i = R (Cos θ + i Sin θ ),G j = R (Cos θ - i Sin θ )。

若AR(p ) 过程平稳,则 |G i | < 1,所以必有R <1。

那么随着k 的增加,G i k = R k (Cosk θ + i Sink θ ),G j k = R k (Cosk θ - i Sink θ ),自相关函数(2.36)式中的相应项G i k , G j k 将按正弦振荡形式衰减(欠阻尼情形)。

实际中的平稳自回归过程的自相关函数常是由指数衰减和正弦衰减两部分混合而成。

③ 从(2.36)式可以看出,当特征方程的根取值远离单位圆时,k 不必很大,自相关函数就会衰减至零。

④ 有一个实数根接近1时,自相关函数将衰减的很慢,近似于线性衰减。

当有两个以上的根取值接近1时,自相关函数同样会衰减的很慢。

a. 两个特征根为实根b. 两个特征根为共轭复根图2.6 AR(2) 过程的自相关函数3. 移动平均过程的自相关函数 (1) MA(1) 过程的自相关函数。

对于MA(1)过程x t = u t + θ1 u t -1 有γk = E(x t x t - k ) = E [(u t + θ1 u t -1) (u t - k + θ1 u t -k -1)] 当k = 0时,γ0 = E(x t x t ) = E [(u t + θ1 u t -1) (u t + θ1 u t -1)]= E (u t 2 + θ1 u t u t -1 + θ1 u t u t -1 + θ12 u t -12 ) = (1 + θ12 ) σ 2 当k = 1时γ1 = E(x t x t - 1) = E [(u t + θ1 u t -1) (u t – 1 + θ1 u t – 2 )]= E (u t u t -1 + θ1 u t -12 + θ1 u t u t -2 + θ12 u t -1 u t -2) = θ1 E (u t -1) 2 = θ1 σ 2 当 k > 1 时,γk = E [(u t + θ1 u t -1) (u t – k + θ1 u t – k -1)] = 0 综合以上三种情形,MA(1)过程自相关函数为ρk =0γγk = 2111θθ+ , k = 1 0 , k > 1,见图2.7。

θ1 > 0 θ1 < 0图2.7 MA(1)过程的自相关函数可见MA(1) 过程的自相关函数具有截尾特征。

当k > 1时,ρk = 0。

(2) MA(q ) 过程的自相关函数 MA(q ) 过程的自相关函数是 ρk =222212211...1...q qk q k k k θθθθθθθθθθ++++++++-++, k = 1, 2, …, q ,0 k > q , 当k > q 时,ρk = 0,说明 ρk , k = 0, 1, … 具有截尾特征。

(注意:模型移动平均项的符号以及这里 ρk 的符号正好与Box-Jenkins 书中的符号相反,这样表示的好处是保持与计算机输出结果一致。

相关主题