第三章时间序列分析.ppt

第一节线性差分方程 一 后移算子B定义为 从而 前面的MA m 模型 AR n 模型和ARMA n m 模型可分别表示为 其中 二 线性差分方程 差分方程的通解为 可将写成 这里 这里 C t 是齐次方程通解解 I t 是特解 三 齐次方程解的计算 假定G1 G2 Gn是互不相同 则在时刻t的通解 其中Ai为常数 可由初始条件确定 无重根考虑齐次差分方程 重根设 有d个相等的根 可验证通解为 对一般情形 因此 齐次方程解是由衰减指数项Gt 多项式tj 衰减正弦项Dt ksin 2 f0t F 以及这些函数的组合混合生成的 齐次方程解便是 定义 设零均值平稳序列 第二节格林函数 Green sfunction 和平稳性 Stationarity 一 格林函数 Green sfunction 能够表示为 则称上式为平稳序列 的传递形式 式中的加权系数 称为格林 Green 函数 其中 格林函数的含义 格林函数是描述系统记忆扰动程度的函数 1 式可以记为 其中 式 1 表明具有传递形式的平稳序列可以由现在时刻以前的白噪声通过系统 的作用而生成 是j个单位时间以前加入系统的干扰项对现实响应的权 亦即系统对的 记忆 一 AR 1 系统的格林函数 由AR 1 模型 即 则AR 1 模型的格林函数 例 下面是参数分别为0 9 0 1和 0 9的AR 1 系统对 扰动的记忆情况 演示试验 AR n 模型 即 其中 的平稳性条件为 的根在单位圆外 或 的根在单位圆内 AR n 系统的平稳性条件 请同学们观察平稳性AR n 与非平稳性AR n 的区别 格林函数与AR n 系统的平稳性 平稳性的涵义就是干扰项对系统的影响逐渐减弱 直到消失 对于一个AR n 系统 将其写成格林函数的表示形式 如果系统是平稳的 则预示随着j 扰动的权数 上面结论也可以用来求AR n 系统的系数平稳性条件 请同学们思考MA m 系统的平稳性条件 ARMA模型格林函数的通用解法 ARMA n m 模型 且 则 令 则 化为 比较等式两边B的同次幂的系数 可得 由上式 格林函数可从 开始依次递推算出 例 求AR 2 1 系统的格林函数 是零均值平稳序列 如果白噪声序列 第三节逆函数和可逆性 Invertibility 能够表示为 一 逆函数的定义 设 则称上式为平稳序列 式中的加权系数 称为逆函数 ARMA n m 模型逆函数通用解法对于ARMA n m 模型的逆函数求解模型格林函数求解方法相同 令 二 ARMA模型的逆函数 的逆转形式 则平稳序列 可表示为 由ARMA n m 模型 可得 仍由先前定义的 和 则上式可化为 比较上式两边B的同次幂的系数 得到 即 可从 由此 开始推算出 对于MA m 模型的可逆性讨论与AR n 模型平稳性的讨论是类似的 即 MA m 模型的可逆性条件为其特征方程 的特征根 满足 ARMA n m 系统格林函数与逆函数的关系 在格林函数的表达式中 用 代替 代替 代替 即可得到相对应的逆函数 理论自协方差函数和自相关函数对于ARMA系统来说 设序列的均值为零 则自协方差函数 第四节自相关函数与偏自相关函数 自相关函数 样本自相关函数的计算在拟合模型之前 我们所有的只是序列的一个有限样本数据 无法求得理论自相关函数 只能求样本的自协方差函数和自相关函数 样本自协方差有两种形式 一 自相关函数 则相应的自相关函数为 在通常情况下 我们采用第一种算法 1 AR p 过程自相关函数ACF 1阶自回归模型AR 1 Xt Xt 1 t的k阶滞后自协方差为 1 2 因此 AR 1 模型的自相关函数为 1 2 由AR 1 的稳定性知 1 因此 k 时 呈指数形衰减 直到零 这种现象称为拖尾或称AR 1 有无穷记忆 infinitememory 注意 0时 呈振荡衰减状 Xt 1Xt 1 2Xt 2 t该模型的方差 0以及滞后1期与2期的自协方差 1 2分别为 阶自回归模型AR 2 类似地 可写出一般的k期滞后自协方差 K 2 3 于是 AR 2 的k阶自相关函数为 K 2 3 其中 1 1 1 2 0 1 如果AR 2 稳定 则由 1 2 1知 k 衰减趋于零 呈拖尾状 至于衰减的形式 要看AR 2 特征根的实虚性 若为实根 则呈单调或振荡型衰减 若为虚根 则呈正弦波型衰减 一般地 p阶自回归模型AR p Xt 1Xt 1 2Xt 2 pXt p t k期滞后协方差为 从而有自相关函数 可见 无论k有多大 k的计算均与其 到p阶滞后的自相关函数有关 因此呈拖尾状 如果AR p 是稳定的 则 k 递减且趋于零 其中 1 zi是AR p 特征方程 z 0的特征根 由AR p 平稳的条件知 zi 1 因此 当1 zi均为实数根时 k呈几何型衰减 单调或振荡 当存在虚数根时 则一对共扼复根构成通解中的一个阻尼正弦波项 k呈正弦波衰减 事实上 自相关函数 是一p阶差分方程 其通解为 对MA 1 过程 2 MA q 过程 可容易地写出它的自协方差系数 于是 MA 1 过程的自相关函数为 可见 当k 1时 k 0 即Xt与Xt k不相关 MA 1 自相关函数是截尾的 其自协方差系数为 一般地 q阶移动平均过程MA q 相应的自相关函数为 可见 当k q时 Xt与Xt k不相关 即存在截尾现象 因此 当k q时 k 0是MA q 的一个特征 于是 可以根据自相关系数是否从某一点开始一直为0来判断MA q 模型的阶 二 偏自相关函数 自相关函数ACF k 给出了Xt与Xt 1的总体相关性 但总体相关性可能掩盖了变量间完全不同的隐含关系 例如 在AR 1 随机过程中 Xt与Xt 2间有相关性可能主要是由于它们各自与Xt 1间的相关性带来的 即自相关函数中包含了这种所有的 间接 相关 与之相反 Xt与Xt k间的偏自相关函数 partialautocorrelation 简记为PACF 则是消除了中间变量Xt 1 Xt k 1带来的间接相关后的直接相关性 它是在已知序列值Xt 1 Xt k 1的条件下 Xt与Xt k间关系的度量 从Xt中去掉Xt 1的影响 则只剩下随机扰动项 t 显然它与Xt 2无关 因此我们说Xt与Xt 2的偏自相关系数为零 记为 在AR 1 中 同样地 在AR p 过程中 对所有的k p Xt与Xt k间的偏自相关系数为零 AR p 的一个主要特征是 k p时 k Corr Xt Xt k 0即 k 在p以后是截尾的 一随机时间序列的识别原则 若Xt的偏自相关函数在p以后截尾 即k p时 k 0 而它的自相关函数 k是拖尾的 则此序列是自回归AR p 序列 对于一个k阶AR模型 有 由此得到Yule Walker方程 记为 已知时 由该方程组可以解出 遗憾的是 用该方程组求解时 需要知道自回归过程的阶数 因此 我们可以对连续的k值求解Yule Walker方程 对k 1 2 3 依次求解方程 得 上述 序列为AR模型的偏自相关函数 偏自相关性是条件相关 是在给定 的条件下 和 的条件相关 换名话说 偏自相关 函数是对 和 所解释的相关的度量 之间未被 由最小二乘原理易得 是作为 关于 线性回归的回归系数 如果自回归过程的阶数为n 则对于k n应该有 kk 0 MA 1 过程可以等价地写成 t关于无穷序列Xt Xt 1 的线性组合的形式 或 是一个AR 过程 它的偏自相关函数非截尾但却趋于零 因此MA 1 的偏自相关函数是非截尾但却趋于零的 注意 式只有当 1时才有意义 否则意味着距Xt越远的X值 对Xt的影响越大 显然不符合常理 因此 我们把 1称为MA 1 的可逆性条件 invertibilitycondition 或可逆域 与MA 1 相仿 可以验证MA q 过程的偏自相关函数是非截尾但趋于零的 MA q 模型的识别规则 若随机序列的自相关函数截尾 即自q以后 k 0 k q 而它的偏自相关函数是拖尾的 则此序列是滑动平均MA q 序列 同样需要注意的是 在实际识别时 由于样本自相关函数rk是总体自相关函数 k的一个估计 由于样本的随机性 当k q时 rk不会全为0 而是在0的上下波动 但可以证明 当k q时 rk服从如下渐近正态分布 rk N 0 1 n 式中n表示样本容量 因此 如果计算的rk满足 我们就有95 5 的把握判断原时间序列在q之后截尾 ARMA p q 的自相关函数 可以看作MA q 的自相关函数和AR p 的