一、时间序列分析---滞后算子（lag operator）

§1.基本概念
时间序列是以观测值发生的时期作为标记的数据集合。一般情况下，我们是从某个特定的时间开始采集数据，直到另一个固定的时间为止，我们可以将获得的数据表示为:
$$(y_1,y_2,...,y_T)$$

如果能够从更早的时间开始观测，或者观测到更晚的时期,那么上面的数据区间可以进一·步扩充。相对而言，上述数据只是一个数据的片段，整个数据序列可以表示为:
( . . . , y 1 , y 2 , . . . , y T , . . . ) = { y t } t = − ∞ t = + ∞ (...,y_1,y_2,...,y_T,...)=\{y_t\}_{t=-\infty}^{t=+\infty} (...,y1​,y2​,...,yT​,...)={yt​}t=−∞t=+∞​

例1.1:几种代表性的时间序列
(1)时间趋势本身也可以构成一个时间序列，此时: $y_t=t$ ;
(2) 另一种特殊的时间序列是常数时间序列，即: $y_t=c$，$c$是常数，这种时间的取值不受时间的影响;
(3)在随机分析中常用的一种时间序列是高斯白噪声过程，表示为: $y_t={\epsilon }_t$， { ε t } t = − ∞ t = + ∞ \{ε_t\}_{t=-\infty}^{t=+\infty} {εt​}t=−∞t=+∞​是一个独立随机变量序列，每个随机变量都服从$N(0,{\sigma }^2)$分布。
时间序列之间也可以进行转换,类似于使用函数关系进行转换。它是将输入时间序列转换为输出时间序列。
例1.2:几种代表性的时间序列转换
(1) 假设$x_t$是一个时间序列，假设转换关系为: $y_t=\beta x_t$，这种算子是将一个时间序列的每一个时期的值乘以常数转换为一个新的时间序列。
(2) 假设$x_t$和$w_t$是两个时间序列，算子转换方式为: $y_t=x_t+w_t$，此算子是将两个时间序列求和。
定义1.1: 如果算子运算是将一个时间序列的前一期值转化为当期值，则称此算子为滞后算子，记做$L$。即对任意时间序列$x_t$，滞后算子满足:
$$L(x_t)\equiv x_{t-1}(1)$$

类似地，可以定义高阶滞后算子，例如二阶滞后算子记为$L^2$，对任意时间序列$x_t$，二阶滞后算子满足:
$$L^2(x_t)\equiv L[L(x_t)]=x_{t-2}(2)$$

一般地，对于任意正整数$k$，有：
$$L^k(x_t)=x_{t-k}$$

命题1.1 滞后算子满足线性性质
$$\begin{array}{l}(1)L(\beta x_t)=\beta L(x_t)\\ (2)L(x_t+w_t)=L(x_t)+L(w_t)\end{array}$$

$proof:$
$$\begin{array}{l}(1)L(\beta x_t)=\beta x_{t-1}=\beta L(x_t)\\ (2)L(x_t+w_t)=x_{t-1}+w_{t-1}=L(x_t)+L(w_t)\end{array}$$

由于滞后算子具有上述运算性质和乘法的交换性质，因此可以定义滞后算子多项式，它的作用是通过它对时间序列的作用获得一个新的时间序列，并且揭示这两个时间序列之间的关系。
显然，滞后算子作用到常数时间序列上，时间序列仍然保持常数，即: $L(c)=c$。

§2.一阶差分方程
利用滞后算子，可以将前面的一阶差分方程表示成为滞后算子形式:
$$y_t=\varphi y_{t-1}+w_t=\varphi L{y}_t+w_t(4)$$

也可以表示为：
$$(1-\varphi L{y}_t)=w_t(5)$$

在上述等式两边同时作用算子：$(1+\varphi L+{\varphi }^2{L}^2+\cdots +{\varphi }^t{L}^t)$，可以得到：
$$(1+\varphi L+{\varphi }^2{L}^2+\cdots +{\varphi }^t{L}^t)(1-\varphi L{y}_t)=(1+\varphi L+{\varphi }^2{L}^2+\cdots +{\varphi }^t{L}^t)w_t$$

计算得到：
$$(1-{\varphi }^{t+1}{L}^{t+1})y_t=(1+\varphi L+{\varphi }^2{L}^2+\cdots +{\varphi }^t{L}^t)w_t$$

利用滞后算子性质得到：
$$\begin{gathered} y_t-{\varphi }^{t+1}{L}^{t+1}{y}_t=w_t+\varphi L{w}_t+{\varphi }^2{L}^2{w}_t+\cdots +{\varphi }^t{L}^t{w}_t \\ \Rightarrow y_t-{\varphi }^{t+1}{y}_{-1}=w_t+\varphi w_{t-1}+{\varphi }^2{w}_{t-2}+\cdots +{\varphi }^t{w}_0 \end{gathered}$$

$\Rightarrow y_t={\varphi }^{t+1}{y}_{-1}+w_t+\varphi w_{t-1}+{\varphi }^2{w}_{t-2}+\cdots +{\varphi }^t{w}_0(6)$
上述差分方程的解同利用叠代算法得到的解是一致的。 
注意算子作用后的等式：
$(1+\varphi L+{\varphi }^2{L}^2+\cdots +{\varphi }^t{L}^t)(1-\varphi L{y}_t)=(1-{\varphi }^{t+1}{L}^{t+1})y_t$
即：
$(1+\varphi L+{\varphi }^2{L}^2+\cdots +{\varphi }^t{L}^t)(1-\varphi L{y}_t)=y_t-{\varphi }^{t+1}{y}_{-1}$
如果时间序列$y_t$是有界的，即存在有限的常数$M$，使得任意时间均有: $|y_t|\le M$， 并且$|\varphi |<1$，则上式当中的尾项随着时间增加趋于零，从而有:
$$\underset{t\to \infty }{\lim }(1+\varphi L+{\varphi }^2{L}^2+\cdots +{\varphi }^t{L}^t)(1-\varphi L{y}_t)=y_t(7)$$

如果利用“1”表示恒等算子，则有：
$$\underset{t\to \infty }{\lim }(1+\varphi L+{\varphi }^2{L}^2+\cdots +{\varphi }^t{L}^t)(1-\varphi L)=1(8)$$

记：
$$(1-\varphi L{)}^{-1}=\underset{t\to \infty }{\lim }(1+\varphi L+{\varphi }^2{L}^2+\cdots +{\varphi }^t{L}^t)(9)$$

因此得到了“逆算子”的表达式，这类似于以滞后算子为变量的函数展开式。
定义2.1： 当$|\varphi |<1$时，定义算子$(1-\varphi L)$的逆算子为$(1-\varphi L{)}^{-1}$，它满足：
$(1)(1-\varphi L)(1-\varphi L{)}^{-1}=(1-\varphi L{)}^{-1}(1-\varphi L)=I(10)$
其中$I$表示单位算子，即对任意时间序列$y_t$有：$I(y_t)=y_t$。
(2)在形式上逆算子可以表示为：
$(1-\varphi L{)}^{-1}={\sum }_{j=0}^{\infty }{\varphi }^j{L}^j(11)$
这表示逆算子作为算子运算规则是:对于任意时间序列$y_t$，有:
$(1-\varphi L{)}^{-1}{y}_t={\sum }_{j=0}^{\infty }{\varphi }^j{L}^j{y}_t={\sum }_{j=0}^{\infty }{\varphi }^j{y}_{t-j}$
如果时间序列$y_t$是有界的，则一阶差分方程的解可以表示为:
$y_t=w_t+\varphi w_{t-1}+{\varphi }^2{w}_{t-2}+\cdots +{\varphi }^t{w}_0={\sum }_{j=0}^{\infty }{\varphi }^j{w}_{t-j}$
可以验算上述表达式确实满足一阶线性差分方程。但是解并不唯一，例如对于任意实数$a_0$，下述形式的表达式均是方程的解。
$y_t=a_0{\varphi }^t+{\sum }_{j=0}^{\infty }{\varphi }^j{w}_{t-j}$

§3.二阶差分方程
我们考察二阶差分方程的滞后算子表达式：
$$y_t={\varphi }_1{y}_{t-1}+{\varphi }_2{y}_{t-2}+w_t$$

将其利用滞后算子表示为：
$$(1-{\varphi }_{1}L-{\varphi }_2{L}^2)=w_t(12)$$

对二阶滞后算子多项式进行因式分解，即寻求${\lambda }_1$和${\lambda }_2$使得:
$$(1-{\varphi }_{1}L-{\varphi }_2{L}^2)=(1-{\lambda }_{1}L)(1-{\lambda }_{2}L)=1-({\lambda }_1+{\lambda }_2)L+{\lambda }_1{\lambda }_2{L}^2$$

即：
$$\{\begin{array}{l}{\varphi }_1={\lambda }_1+{\lambda }_2\\ {\varphi }_2=-{\lambda }_1{\lambda }_2\end{array}$$

$\Rightarrow$ ${\lambda }_1,{\lambda }_2$是差分方程对应的特征方程的根：
$${\lambda }^2-{\varphi }_1\lambda -{\varphi }_2=0(13)$$

当特征根${\lambda }_1,{\lambda }_2$落在单位圆内的时候(这也是差分方程的稳定性条件)，滞后算子多项式分解为：
$$\begin{array}{l}(1-{\lambda }_{1}L{)}^{-1}=1+{\lambda }_{1}L+{\lambda }_1^2{L}^2+{\lambda }_1^3{L}^3+\cdots \\ (1-{\lambda }_{2}L{)}^{-1}=1+{\lambda }_{2}L+{\lambda }_2^2{L}^2+{\lambda }_2^3{L}^3+\cdots \end{array}$$

此时二阶差分方程解可表示为：
$$y_t=(1-{\lambda }_{1}L{)}^{-1}(1-{\lambda }_{2}L{)}^{-1}{w}_t$$

注意到算子分式也可以进行分项分式分解:
$$\frac{1}{(1-{\lambda }_{1}L)(1-{\lambda }_{2}L)}=\frac{1}{{\lambda }_1-{\lambda }_2}\left(\frac{{\lambda }_1}{(1-{\lambda }_{1}L)}-\frac{{\lambda }_2}{(1-{\lambda }_{2}L)}\right)$$

将上述表达式带入到二阶差分方程解中:
$$\begin{array}{l}{y}_t=\frac{1}{{\lambda }_1-{\lambda }_2}\left(\frac{{\lambda }_1}{(1-{\lambda }_{1}L)}-\frac{{\lambda }_2}{(1-{\lambda }_{2}L)}\right)w_t\\ =\frac{1}{{\lambda }_1-{\lambda }_2}[{\lambda }_1(1+{\lambda }_{1}L+{\lambda }_1^2{L}^2+{\lambda }_1^3{L}^3+\cdots )-{\lambda }_2(1+{\lambda }_{2}L+{\lambda }_2^2{L}^2+{\lambda }_2^3{L}^3+\cdots )]w_t\\ =\frac{1}{{\lambda }_1-{\lambda }_2}[({\lambda }_1-{\lambda }_2)+({\lambda }_1^2-{\lambda }_2^2)L+({\lambda }_1^3-{\lambda }_2^3)L^2+\cdots ]w_t\\ =(c_1+c_2)w_t+(c_1{\lambda }_1+c_2{\lambda }_2)w_{t-1}+(c_1{\lambda }_1^2+c_2{\lambda }_2^2)w_{t-2}+\cdots \end{array}$$

其中，$c_1=\frac{{\lambda }_1}{{\lambda }_1-{\lambda }_2},c_2=-\frac{{\lambda }_2}{{\lambda }_1-{\lambda }_2}=\frac{{\lambda }_2}{{\lambda }_2-{\lambda }_1}$

利用上述公式，可以得到外生扰动的动态反应乘子为:
$$\frac{\partial y_{t+j}}{\partial w_t}=c_1{\lambda }_1^j+c_2{\lambda }_2^j,j=0,1,2,...(14)$$