Contents
自回归模型
推移算子
为了方便对时间序列的讨论,需要对时间序列的时间指标 t t t 引入向后推移算子 B \mathscr{B} B 。
直观上看,向后推移算子的作用方式可以如下表示:
B X t = X t − 1 \mathscr{B}X_t=X_{t-1} BXt=Xt−1
B 2 X t = B ( B X t ) = B ( X t − 1 ) = X t − 2 \mathscr{B}^2X_t=\mathscr{B}(\mathscr{B}X_t)=\mathscr{B}(X_{t-1})=X_{t-2} B2Xt=B(BXt)=B(Xt−1)=Xt−2
推移算子可以和多项式或无穷级数相结合,下面给出定义:
对任何时间序列 { X t } \{X_t\} { Xt} 和无穷级数 ψ ( z ) = ∑ j = − ∞ ∞ b j z j \psi(z)=\displaystyle\sum_{j=-\infty}^\infty b_jz^j ψ(z)=j=−∞∑∞bjzj ,只要级数 ∑ j = − ∞ ∞ b j X t − j \displaystyle\sum_{j=-\infty}^\infty b_jX_{t-j} j=−∞∑∞bjXt−j 在某种意义下收敛,就定义
ψ ( B ) = ∑ j = − ∞ ∞ b j B j \psi(\mathscr{B})=\displaystyle\sum_{j=-\infty}^\infty b_j\mathscr{B}^j ψ(B)=j=−∞∑∞bjBjψ ( B ) X t = ∑ j = − ∞ ∞ b j B j X t = ∑ j = − ∞ ∞ b j X t − j \psi(\mathscr{B})X_t=\displaystyle\sum_{j=-\infty}^\infty b_j\mathscr{B}^jX_{t}=\displaystyle\sum_{j=-\infty}^\infty b_jX_{t-j} ψ(B)Xt=j=−∞∑∞bjBjXt=j=−∞∑∞bjXt−j
称 B \mathscr{B} B 为时间 t t t 的向后推移算子。
齐次常系数线性差分方程
给定 p p p 个实数 a 1 , a 2 , . . . , a p a_1,a_2,...,a_p a1,a2,...,ap , a p ≠ 0 a_p\neq0 ap=0 ,称
X t − [ a 1 X − 1 + a 2 X t − 2 + . . . + a p X t − p ] = 0 , t ∈ Z X_t-[a_1X_{-1}+a_2X_{t-2}+...+a_pX_{t-p}]=0 \ , \ \ \ \ t\in\Z Xt−[a1X−1+a2Xt−2+...+apXt−p]=0 , t∈Z
为 p p p 阶齐次常系数线性差分方程,简称为齐次差分方程。满足该方程的实值时间序列 { X t } \{X_t\} { Xt} 称为该方程的解。
利用推移算子 B \mathscr{B} B,定义多项式 A ( z ) = 1 − ∑ j = 1 p a j z j A(z)=1-\displaystyle\sum_{j=1}^pa_jz^j A(z)=1−j=1∑pajzj ,可以把齐次差分方程等价的写为 A ( B ) X t = 0 A(\mathscr{B})X_t=0 A(B)Xt=0 ,称多项式 A ( z ) A(z) A(z) 为齐次差分方程的特征多项式。
可以看出,方程 A ( B ) X t = 0 A(\mathscr{B})X_t=0 A(B)Xt=0 的解 { X t } \{X_t\} { Xt} 可以由它的 p p p 个初值 X 0 , X 1 , . . . , X p − 1 X_0,X_1,...,X_{p-1} X0,X1,...,Xp−1 逐步递推得到,因此 { X t } \{X_t\} { Xt} 可以由这 p p p 个初值唯一决定。由于初值是可以任意选择的,所以齐次差分方程的解有无穷多个,且可以看出齐次差分方程的解空间是一个线性空间。下面找出方程的所有的线性独立的解。
根据代数基本定理,设多项式 A ( z ) A(z) A(z) 有 k k k 个互不相同的根 z 1 , z 2 , . . . , z k z_1,z_2,...,z_k z1,z2,...,zk ,其中 z j z_j zj 是 r ( j ) r(j) r(j) 重根, ∑ j = 1 k r ( j ) = p \displaystyle\sum_{j=1}^kr(j)=p j=1∑kr(j)=p ,
A ( z ) = ∏ j = 1 k ( 1 − z z j ) r ( j ) , A(z)=\prod_{j=1}^k(1-\frac{z}{z_j})^{r(j)}, A(z)=j=1∏k(1−zj