线性平稳时间序列分析

中心化的AR(2)模型： X t1X t 12X t 2t
非中心化的AR(2)模型：X t c1 X t 12 X t 2t
其中εt为随机扰动，一般为零均值的白噪声序列。
AR(2)模型的平稳性条件
X t c1 X t 12 X t 2t
平稳条件：对应齐次差分方程的特征根在单位圆内
对应齐次差分方程的特征多项式
其根互 为倒数
() p 1p 1 2p 2 L p
AR模型平稳性判别方法
特征根判别
AR(p)模型平稳的充要条件是它的p个特征根都在单位 圆内
根据特征根和自回归系数多项式的根成倒数的性质， 等价判别条件是该模型的自回归系数多项式的根都在 单位圆外
AR(1)模型：一阶自回归模型
描述这种关系的数学模型就是一阶自回归模型，简 记为AR(1)，即
Xt 1Xt1t
其中Xt为零均值(即中心化处理 后的)平稳序列。φ1为Xt对Xt-1的 依赖程度,εt为随机扰动，一般 为零均值的白噪声序列。
AR(1)的中心化变换
一般情形： Xt c1Xt1t
z t a 1 z t 1 a 2 z t 2 a p z t p h ( t )
非齐次线性差分方程的通解 zt
齐次线性差分方程的通解 z 程的特解 z t 之和
t
和非齐次线性差分方
zt ztzt
一阶差分方程
yt yt1t
P33
用递归替代法解差分方程：假设已知y-1和ω的各期
复根场合 z t r t( c 1 e i t c 2 e i t ) c 3 t 3 c p t p
z t a 1 z t 1 a 2 z t 2 a p z t p h ( t )
非齐次线性差分方程的解
非齐次线性差分方程的特解
使得非齐次线性差分方程成立的任意一个解
非中心化的MA(q)模型：
X t t 1 t 1 2 t 2 L q t q ,t: W N 0 ,2
引进延迟算子， MA(q)模型又可以简记为：
Xt (B)t
q阶移动平均系数多项式：
(B ) 1 1 B 2 B 2 L q B q
MA(q)模型的统计性质
特征方程： 212 0
特征根:
1 12 42
2
AR(2)模型平稳
1
AR(2)模型的平稳性条件
平稳域 AR(2)平稳性判别：
特征根 平稳域
1 ,22 1 ， 且 2 1 1
考察下列模型的平稳性：
(3 )X tX t 1 0 .5 X t 2t
(4 )X t X t 1 0 .5 X t 2t
ω将对系统产生持久性影响。
线性过程
定义：{Xt}称为线性过程，若 Xt Gjt j ，其中
j
{εt}是白噪声序列，系数序列{Gj}满足
G
2 j
。
j
系统是因果性的：若系数序列Gj满足Gj=0, j<0，即
Xt Gj tj G0 t G1 t1L j0
◦ 定理3.1：线性过程肯定是平稳过程，且是均方收敛的。
线性差分方程
z t a 1 z t 1 a 2 z t 2 a p z t p h ( t )
齐次线性差分方程
z t a 1 z t 1 a 2 z t 2 a p z t p 0
z t a 1 z t 1 a 2 z t 2 L a p z t p 0
X t t 1 t 1 2 t 2 L q t q ,t: W N 0 ,2
常数均值：模型两边求期望可得
EX t
常数方差：
varXt vart1t1Lqtq
1122 2Lq 2
2 a
【注】MA(q)模型一定为平稳模型。
MA(q)模型的可逆性
可逆MA模型定义
序列的期望和方差如何求？
AR(p) 模型：一般自回归模型
中心化的AR(p)模型：
Xt 1Xt12Xt2 L pXtp t
t : WN 0,2
EXs t 0, st
说明当前期的随 机扰动与过去的
序列值无关
非中心化的AR(p)模型：
X t c 1 X t 1 2 X t 2 L p X t p t
用ARMA模型去近似地描述动态数据在实际应用 中有许多优点，例如它是线性模型，只要给出少 量参数就可完全确定模型形式；另外，便于分析 数据的结构和内在性质，也便于在最小方差意义 下进行最佳预测和控制。
线性过程
方法性工具
这些工具会使得时间序列模型表达和分析更为简 洁和方便。 延迟算子 线性差分方程
引进延迟算子，ARMA(p,q)模型又可以简记为 ：
(B )Xt c (B )t
p阶自回归系数多项式：
( B ) 1 1 B 2 B 2 L p B p
q阶移动平均系数多项式：
( B ) 1 1 B 2 B 2 q B q
AR、MA和ARMA之间的关系
ARMA(p,q)模型： X t c 1 X t 1 L p X t p t 1 t 1 L q t q
yt yt1t P33
动态乘子(动态乘子为输入ω对输出yt的影响)
yt t 或ytj j
0
t
当0<φ<1，动态乘子按几何方式衰减到零；当-1<φ<0，动 态乘子振荡衰减到零；
当φ>1，动态乘子指数增加；当φ<-1，动态乘子发散性振 荡；
当︱φ︱<1，动态系统稳定，即给定的ω的影响将逐渐消 失；当︱φ︱>1，动态系统发散；当︱φ︱=1，输入变量
若一个MA模型能够表示成无穷阶的自回归模型，则称 该MA模型称为可逆的。
例：（1）Xt t 2t1 （2）Xt t 0.5t1
（1）Xt 12Bt （2）Xt 10.5Bt
（ 1 ） t 1/12BXt
（ 2） t 1/10.5BXt 0.5BnXt 0.5nXtn
n0
n0
ARMA模型
自回归移动平均模型
Autoregressive-Moving Average Model
ARMA模型的背景
一个系统，如果它在t时刻的响应 Xt 不仅与其以前 时刻的响应有关，而且还与其以前时刻进入系统的 扰动存在着一定的相关关系，那么这个系统就是自 回归移动平均系统，相应的模型记作ARMA模型。
在此模型下，一个影响系统的扰动εt 被“牢记”一 定时期，从而影响系统的后继行为。正是系统的这 种动态性，引起了时间序列中的依存关系，从而决 定了序列中的依存关系不能用普通静态回归模型来 描述，而只能用ARMA模型。
此时
EXt
c
11
0
中心化：令Yt=Xt- ，Yt即为Xt的中心化序列，
此时有
E Yt 0
AR模型平稳性的判别
判别原因
AR模型是常用的平稳序列的拟合模型之一， 但并非所有的AR模型都是平稳的。
判别方法
特征根判别法
AR(1)模型的平稳性条件
Xt c1Xt1t
平稳条件：对应齐次差分方程的特征根在单位圆内
X t 1 X t 1 2 X t 2 L p X t p t
AR(p) 的自回归系数多项式
引进延迟算子，中心化的AR(p)模型又可以简记为
(B)Xt t
其 中 (B)11B2B2LpBp
自回归系数多项式
( u ) 1 1 u 2 u 2 L p u p
Average Model)
AR(p)模型：p 阶自回归模型
AR(1)模型的背景
如果时间序列是独立的，没有任何依赖关系，这样 的资料所揭示的统计规律就是事物独立的随机变动， 系统无记忆能力。如果情况不是这样，资料之间有 一定的依存性，那么最简单的关系就是后一时刻的 行为主要与其前一时刻的行为有关，而与其前一时 刻以前的行为无直接关系，即已知Xt-1，Xt主要与 Xt-1相关。用记忆性来说，就是最短的记忆，即一 期记忆，也就是一阶动态性。
值
y t t 1 y 1 t 0 t 11 L t
动态乘子
yt t 或ytj j
0
t
动态乘子为输入ω对输出yt的影响，依赖于j，即输入ωt
和输出yt+j观察值之间的时间间隔。
当参数φ取不同的值，系统最后的状态也不同。
一阶差分方程 y t t 1 y 1 t 0 t 11 L t
延迟算子
定义：设B为一步延迟算子，如果当前序列乘以
一个延迟算子，就表示把当前序列值的时间向过
去拨一个时刻，即 BXt=Xt-1。
性质： B0 1
B(c
X
t
)
c
B(
X
t
)
c
X
t
1,
c为任意常数
B(
Xt
Yt
)
X t 1
Yt1
(1
B)n
n
(1)i Cni Bi
Bn Xt
i0
Xtn
线性差分方程
线性过程的因果性
在应用时间序列分析去解决实际问题时，所使用 的线性过程是因果性的，即：
X t G j tjtG 1t 1 L且G 2 j
j 0
j 0
用延迟算子表示：
Xt
GjBj
t
GBt
j0
条件： G j j0
线性过程逆转形式
用t时刻及其以前时刻的Xt-j(j=0,1, …)来表示白噪声
当p=0时，ARMA(p,q) 模型就退化为MA(q)模型； 当q=0时，ARMA(p,q) 模型就退化为AR(p)模型； AR(p)模型和MA(q) 模型实际上是ARMA(p,q)模型
MA(1)模型：一阶移动平均模型
如果一个系统在t时刻的响应Xt仅与其前一时刻进入系 统的扰动εt-1存在着一定的相关关系，描述这种关系的 数学模型就是一阶移动平均模型，记作MA(1)，即