【时间序列分析】04. 自回归模型

自回归模型

推移算子

为了方便对时间序列的讨论,需要对时间序列的时间指标 t t t 引入向后推移算子 B \mathscr{B} B

直观上看,向后推移算子的作用方式可以如下表示:
B X t = X t − 1 ​ \mathscr{B}X_t=X_{t-1}​ BXt=Xt1

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(Xt1)=Xt2

推移算子可以和多项式或无穷级数相结合,下面给出定义:

对任何时间序列 { 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=bjXtj 在某种意义下收敛,就定义
ψ ( 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=bjXtj

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[a1X1+a2Xt2+...+apXtp]=0 ,    tZ
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)=1j=1pajzj ,可以把齐次差分方程等价的写为 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,...,Xp1 逐步递推得到,因此 { 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=1kr(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=1k(1zj

  • 2
    点赞
  • 12
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值