第二章

2.1 一般线性过程

一、线性平稳序列

线性平稳序列是白噪声的线性组合得到的序列。 最简单的线性平稳序列是有限运动平均。

{εt}={εt: tZ}\{\varepsilon_t\}=\{\varepsilon_t: \ t \in \mathbb Z\}WN(0,σ2)WN(0,\sigma^2).对于非负整数qq和常数a0,a1,,aq(a00,aq0)a_0,a_1,\ldots,a_q(a_0 \neq 0, a_q \neq 0),称

Xt=j=0qajεtj=a0εt+a1εt1++aqεtq,  tZX_t = \sum_{j=0}^q a_j \varepsilon_{t-j} = a_0\varepsilon_t +a_1\varepsilon_{t-1}+\cdots + a_q\varepsilon_{t-q}, \ \ t \in \mathbb Z

是白噪声{εt}\{\varepsilon_t\}的(有限)运动平均或者滑动平均,简称MA(Moving  Average)MA(Moving \ \ Average)
注:可将{Xt}\{X_t\}视为线性滤波器的输出,白噪声看成驱动系统的扰动序列(激励).

MAMA的平稳性
EXt=0E X_t = 0
E(Xt+kXt)={σ2j=0qkajaj+k,0kq,0,k>qE (X_{t+k} X_t) = \begin{cases} \sigma^2 \sum_{j=0}^{q-k} a_j a_{j+k}, & 0 \leq k \leq q, \\ 0, & k > q \end{cases}

可见{Xt}\{X_t\}平稳。γk=0,k>q\gamma_k=0,\forall k>q,称这样的序列为qq相关的。

时间序列——第二章1

线性滤波器的原始数据与滤波结果是一种算术运算,即用加减乘除等运算实现。包括 方框滤波(boxFilter)、均值滤波(blur)、高斯滤波(GaussianBlur);
非线性滤波器的原始数据与滤波结果是一种逻辑关系,即用逻辑运算实现。包括最大值滤波中值滤波(medianBlur)和双边滤波(bilateralFilter)。

参考资料:
常见线性滤波
线性滤波和非线性滤波

二、线性过程的两种等价形式

1、 传递形式

时间序列分析中建立随机模型的思想:
将顺序值之间高度依赖的时间序列{Xt}\{ X_t \}看成由一系列独立“冲击”序列所生成.
通常将“冲击”序列理想化为白噪声过程{εt, tZ}\{\varepsilon_t, \ t \in \mathbb Z\},时间序列Xt{ X_t }取为白噪声的加权和.

WoldWold分解式(正交分解):
任意零均值纯非确定的平稳过程都可表示为线性形式(无穷阶滑动平均):
Xt=j=0φjεtjtZX_t=\sum_{j=0}^{\infty}φ_j\varepsilon_{t-j},t\in \mathbb Z
称为沃尔德系数(格林函数、传递函数),其中{εt, tZ}\{\varepsilon_t, \ t \in \mathbb Z\}是白噪声序列。

注1:上式可表示为算子形式:Xt=φ(B)εtX_t=φ(\mathscr B)\varepsilon_t,其中φ(B)=1+j=1φjBj=j=0φjBjφ(\mathscr B)=1+\sum_{j=1}^{\infty}φ_j\mathscr B^j=\sum_{j=0}^{\infty}φ_j\mathscr B^j
称为传递函数或沃尔德系数的生成函数.

注2:对动态数据进行适当的预处理,可将非平稳序列 平稳化与零均值化.

注3:上式的工程解释:Xt=φ(B)εtX_t=φ(\mathscr B)\varepsilon_t (线性滤波器模型)
将时间序${X_t } $ 视为线性滤波器的输出,白噪声看成驱动系统的扰动序列(激励)
时间序列——第二章1

2、自回归形式

在适当条件下{Xt,tZ}\{X_t,t\in \mathbb Z \}可表示为线性形式
Xt=εt+j=1πjXtj,tZX_t=\varepsilon_t+\sum_{j=1}^{\infty}\pi_jX_{t-j},t\in \mathbb Z
εt=(1j=1πjBj)Xt=π(B)Xt,tZ\varepsilon_t=(1-\sum_{j=1}^{\infty}\pi_j\mathscr B^j)X_t=\pi(\mathscr B)X_t,t\in \mathbb Z
时间序列——第二章1

理解:回归这一词的意思是从已有的数据发现规律。与过去的XtX_t有联系,即相关。

3、 两者之间的关系

Xt=φ(B)εt,tZX_t=φ(\mathscr B)\varepsilon_t,t\in \mathbb Z(*)εt=π(B)Xt,tZ\varepsilon_t=\pi(\mathscr B)X_t,t\in \mathbb Z(**)

​ 可以看出两者互为逆转形式。

​ 将生成函数φ(B)φ(\mathscr B)作用于(**)

Xt=φ(B)εt=φ(B)π(B)XtX_t=φ(\mathscr B) \varepsilon_t =φ(\mathscr B)\pi(\mathscr B)X_t
​ 得φ(B)π(B)=1φ(\mathscr B)\pi(\mathscr B)=1
​ 或φ(B)=π(B)1,π(B)=φ(B)1φ(\mathscr B)=\pi(\mathscr B)^{-1},\pi(\mathscr B)=φ(\mathscr B)^{-1}

传递形式:
时间序列——第二章1

逆转形式:
时间序列——第二章1

:若平稳序列得自回归形式为Xt=πXt1+εtX_t =\pi X_{t-1}+\varepsilon_tπ(B)Xt=εt\pi(\mathscr B)X_t=\varepsilon_t
其中 π(B)=1πB,\pi(\mathscr B)=1-\pi \mathscr B,有:
φ(B)=π(B)1=11πB=1+πB+π2B2+...+πjBj+...φ(\mathscr B)=\pi(\mathscr B)^{-1}=\frac{1}{1-\pi \mathscr B}=1+\pi\mathscr B+\pi^2\mathscr B^2+...+\pi^j\mathscr B^j+...

求得传递形式为:
Xt=π1(B)εt=φ(B)=j=0πjεtjX_t=\pi^{-1}(\mathscr B)\varepsilon_t=φ(\mathscr B)=\sum_{j=0}^{\infty}\pi^j\varepsilon_{t-j}

三、线性过程的均值函数与自协方差函数

2.1.1 定理 (Levi 单调收敛定理)

如果非负随机变量序列{ξn}\{\xi_n \}单调不减: 0ξ1ξ20\leq \xi_1 \leq \xi_2 \leq \cdots
则当ξnξ  a.s.ξn→ξ\ \ a.s.时, 有 limnEξn=Eξ\lim_{n\to \infty} E\xi_n=E\xi

由于单调不减序列必有极限(极限也可以是正无穷),所以上面的随机变量ξ\xi一定存在。但是可能会有p(ξ=)>0p(\xi=\infty)>0这时Eξ=E\xi=\infty,可以取值\infty的随机变量通常被称为广义随机变量

对于任何时间序列 {Yt}\{Y_t\}, 利用单调收敛定理得到

E[t=Yt]=limnE[t=nnYt]=limnt=nnEYt=t=EYt.\begin{aligned} E\left[ \sum_{t=-\infty}^{\infty} |Y_t| \right]=& \lim_{n\to \infty} E \left[ \sum_{t=-n}^{n} |Y_t| \right] \\ =& \lim_{n\to \infty}\sum_{t=-n}^{n} E|Y_t| = \sum_{t=-\infty}^{\infty} E |Y_t|. \end{aligned}

2.1.2定理 (控制收敛定理)

如果随机变量序列{ξn}\{\xi_n\}满足ξnξ0  a.s.|\xi_n|\leq \xi_0\ \ a.s.Eξ0<E|\xi_0|< \infty,
则当ξnξ  a.s.\xi_n \to \xi\ \ a.s.时,有 Eξ<E|\xi|<\infty,并且limnEξnEξlim_{n\to \infty}E\xi_n \to E\xi.

在上面两个定理成立条件下均有,

limnEξn=Elimnξn.\begin{aligned} \lim_{n \to \infty} E \xi_n = E \lim_{n \to \infty} \xi_n. \end{aligned}

即期望与极限可以交换次序。

极限与积分交换次序的三个定理:

Lebesgue控制收敛定理是这样的描述的:
{fn}\{f_{n}\}是可测集EE上的可测函数列,FF是可测函数列的控制函数并且可积,{fn}\{f_{n}\}依测度收敛到f,则ff在E上可积并且
limnEfndx=Efdx\lim_{n \to \infty} \int_{E}f_{n}dx=\int_{E}fdx

Levi定理是这样叙述的:
{fn}\{f_{n}\}是可测集EE上的可测函数列,并且是非负的单调递增的即0f1(x)f2(x)fn(x)fn+1(x)0\leq f_{1}(x)\leq f_{2}(x) \cdots f_{n}(x) \leq f_{n+1}(x)\leq \cdots,那么{fn(x)}\{f_{n}(x)\}收敛(几乎处处收敛)于一个可积函数f(x)f(x)并且上面5楼的交换仍然成立。
limnEfndx=Efdx\lim_{n \to \infty} \int_{E}f_{n}dx=\int_{E}fdx

Fatou引理是这样说的:
{fn}\{f_{n}\}是可测集EE上的非负可测函数列, 有下面的非严格等式交换
Elimfn(x)dxlimEfn(x)dx\int_{E}\underline{\lim}f_{n}(x)dx \leq \underline{\lim}\int_{E}f_{n}(x)dx
详情可见
积分的极限定理

2.1.3定理

设时间序列{Xt}\{X_t\}满足suptEXt<+sup_tE|X_t|<+\infty,如果j=ψj<+\sum_{j=-\infty}^{\infty}|\psi_j|<+\infty,则序列:
ψ(B)Xt=j=ψjBjXt=j=ψjXtj  ()\psi(\mathscr B)X_t=\sum_{j=-\infty}^{\infty}\psi_j\mathscr B^jX_t=\sum_{j=-\infty}^{\infty}\psi_jX_{t-j}\ \ (*)
依概率1收敛(sup 上确界supremum)

进而,如果suptEXt2<+sup_tE|X_t|^2<+\infty,则序列均方收敛到同一极限。

:由单调收敛定理和suptEXt<+sup_tE|X_t|<+\infty,有
E(j=ψjXtj)=limnE(j=nnψjXtj)limnE(j=nnψj)suptEXtj<\begin{matrix}E\left(\sum_{j=-\infty}^{\infty}|\psi_j||X_{t-j} |\right)\end{matrix}=lim_{n\to \infty}E\left(\sum_{j=-n}^{n}|\psi_j||X_{t-j} |\right) \\≤lim_{n\to \infty}E\left(\sum_{j=-n}^{n}|\psi_j|\right)sup_tE|X_{t-j}|<\infty

因此j=ψjXtj\sum_{j=-\infty}^{\infty}|\psi_j||X_{t-j}|ψ(B)Xt=j=ψjBjXt=j=ψjXtj\psi(\mathscr B)X_t=\sum_{j=-\infty}^{\infty}\psi_j\mathscr B^jX_t=\sum_{j=-\infty}^{\infty}\psi_jX_{t-j}依概率1有限

如果suptEXt2<+,n>m>0,m,nsup_tE|X_t|^2<+\infty,n>m>0,则当m,n\to\infty时,有
Em<j<nψjXtj2=m<j<nm<k<nψjψ^kXtjX^tksuptEXt2(m<j<nψj)20\begin{matrix}E\left|\sum_{m<|j|<n}\psi_jX_{t-j}\right|^2&=&\sum_{m<|j|<n}\sum_{m<|k|<n}\psi_j\hat \psi_kX_{t-j}\hat X_{t-k}\\ &≤&sup_tE|X_t|^2(\sum_{m<|j|<n}|\psi_j|)^2\to0\end{matrix}

由Cauchy收敛准则,序列()(*)均方收敛。

令???? 表示序列()(*)的均方极限,由Fatou引理,有
ESψ(B)Xt2=Elimninf(sj=nnψjXtj)infE(sj=nnψjXtj)=0\begin{matrix}E\left|S-\psi(\mathscr B)X_{t}\right|^2&=&Elim_{n\to \infty}inf\left( s-\sum_{j=-n}^{n} \psi_j X_{t-j} \right)\\ &≤&infE\left( s-\sum_{j=-n}^{n} \psi_j X_{t-j} \right)=0\end{matrix}

即证明了S和ψ(B)Xt\psi(\mathscr B)X_t依概率1相等

推论2.1

如果随机变量序列{Yn}\{Y_n\}满足i=EYi\sum_{i=-\infty}^\infty E|Y_i|\leq \infty,则 i=Yi,a.s.\sum_{i=-\infty}^\infty Y_i ,a.s.收敛, Ei=Yi<E|\sum_{i=-\infty}^\infty Y_i| < \infty
En=Yn=n=EYn\begin{aligned} E \sum_{n=-\infty}^\infty Y_n = \sum_{n=-\infty}^\infty E Y_n \end{aligned}

证明 令
η=i=Yi\eta = \sum_{i=-\infty}^\infty |Y_i|
η\eta是随机变量, 由单调收敛定理知Eη<,η<,a.s.E\eta < \infty, \eta < \infty, a.s.

ξn=i=nnYi\xi_n = \sum_{i=-n}^n Y_i

ξnη,|\xi_n| \leq \eta,

记事件A={η<}A = \{ \eta < \infty \}, 则P(A)=1P(A) = 1, 对ωA\omega \in A, 有i=Yi(ω)<\sum_{i=-\infty}^\infty |Y_i(\omega)| < \infty所以i=Yi(ω)\sum_{i=-\infty}^\infty Y_i(\omega)收敛, limnξn(ω)\lim_{n\to\infty} \xi_n(\omega)收敛,记为$\xi(\omega); \omega \notin A,, 定义\xi(\omega) = 0, 则\xi是随机变量, 且\lim_{n\to\infty} \xi_n = \xi, \text{ a.s.}$

ξ=i=ξi.\xi = \sum_{i=-\infty}^\infty \xi_i .

由控制收敛定理 limnEξn=Eξ\lim_{n\to\infty} E \xi_n = E \xi,

limnEξn=limnEi=nnYi=limni=nnEYi=i=EYi,\lim_{n\to\infty} E \xi_n = \lim_{n\to\infty} E \sum_{i=-n}^n Y_i = \lim_{n\to\infty} \sum_{i=-n}^n E Y_i = \sum_{i=-\infty}^\infty E Y_i,

于是有Ei=ξi=i=Eξi.E \sum_{i=-\infty}^\infty \xi_i = \sum_{i=-\infty}^\infty E \xi_i .

四、线性平稳序列

如果实数列{aj}\{a_j\}满足 j=aj<,\sum_{j=-\infty}^{\infty} |a_j| < \infty,则称{aj}\{a_j\}绝对可和的. 记{aj}l1\{a_j\} \in l_1

注意:{aj}l1\{a_j\} \in l_1{aj}l2\{a_j\} \in l_2(即jaj2<\sum_j a_j^2 < \infty). 反之不一定成立。

对于绝对可和的实数列{aj}\{a_j\}, 定义零均值白噪声{εt}\{\varepsilon_t\}的无穷滑动和如下
Xt=j=ajεtj,  tZ.X_t=\sum_{j=-\infty}^{\infty} a_j \varepsilon_{t-j}, \ \ t \in \mathbb Z.
{Xt}\{X_t\}是平稳序列。 EXt=0E X_t = 0,
γk=σ2j=ajaj+k,kZ.\begin{aligned} \gamma_k = \sigma^2 \sum_{j=-\infty}^\infty a_j a_{j+k}, k \in \mathbb Z .\end{aligned}

ck=j=ajaj+k,kZc_k = \sum_{j=-\infty}^\infty a_j a_{j+k}, k \in \mathbb Z

线性序列的a.s.收敛性

作为无穷和{Xt}\{X_t\}有没有定义? 由SchwarzSchwarz不等式,
Eεt=E(εt1)Eεt21=σ,E|\varepsilon_t| = E(|\varepsilon_t| \cdot 1) \leq \sqrt{E \varepsilon_t^2 \cdot 1} = \sigma,
所以
j=Eajεtj=j=ajEεtjσj=aj<\sum_{j=-\infty}^\infty E |a_j \varepsilon_{t-j}| = \sum_{j=-\infty}^\infty |a_j| E|\varepsilon_{t-j}| \leq \sigma \sum_{j=-\infty}^\infty |a_j| < \infty,
由推论2.1可知j=ajεtj\sum_{j=-\infty}^\infty a_j \varepsilon_{t-j} a.s.收敛。

线性序列的L1L_1收敛性

级数的余项的L1L_1
Ej>Najεtjj>NajEεtjσj>Naj0\begin{aligned} & E \left| \sum_{|j|>N} a_j \varepsilon_{t-j} \right| \leq \sum_{|j|>N} |a_j| E |\varepsilon_{t-j}| \leq \sigma \sum_{|j|>N} |a_j| \to 0 \end{aligned}
即线性平稳列在L1L_1意义下收敛。

线性序列的平稳性

由推论2.1可知

Ej=ajεtj=j=ajEεtj=0\begin{aligned} E \sum_{j=-\infty}^\infty a_j \varepsilon_{t-j} &= \sum_{j=-\infty}^\infty a_j E \varepsilon_{t-j} = 0 \end{aligned}

SchwarzSchwarz不等式知Eϵtjϵt+klσ2E|\epsilon_{t-j} \epsilon_{t+k-l}| \leq \sigma^2, 于是
j=l=Eajalϵtjϵt+klσ2j=l=aj  al=σ2(j=aj)2<\sum_{j=-\infty}^\infty \sum_{l=-\infty}^\infty E |a_j a_l \epsilon_{t-j} \epsilon_{t+k-l}| \leq \sigma^2 \sum_{j=-\infty}^\infty \sum_{l=-\infty}^\infty |a_j|\; |a_l| = \sigma^2 \left( \sum_{j=-\infty}^\infty |a_j| \right)^2 < \infty,
由推论2.1,
EXtXt+k=Ej=ajεtjl=alεt+kl=j=l=ajalE(εtjεt+kl)=σ2j=ajaj+k\begin{aligned} E X_t X_{t+k} &= E \sum_{j=-\infty}^\infty a_j \varepsilon_{t-j} \sum_{l=-\infty}^\infty a_l \varepsilon_{t+k-l} \\ &= \sum_{j=-\infty}^\infty \sum_{l=-\infty}^\infty a_j a_l E(\varepsilon_{t-j} \varepsilon_{t+k-l}) \\ &= \sigma^2 \sum_{j=-\infty}^\infty a_j a_{j+k} \end{aligned}
{Xt}\{X_t\}是平稳序列, EXt=0γk=σ2j=ajaj+k.E X_t=0, \gamma_k = \sigma^2 \sum_{j=-\infty}^\infty a_j a_{j+k}.

线性序列的L2收敛性

{aj}l2\{a_j\} \in l_2,即jaj2<\sum_j a_j^2 < \infty,则
Xt=j=ajεtj (L2)\begin{aligned} X_t=\sum_{j=-\infty}^{\infty} a_j \varepsilon_{t-j}\ (L^2) \end{aligned}
也是平稳序列。期望为零,自协方差函数同上。

XtX_t定义的无穷级数是L2L^2收敛的。 证明需要应用Hilbert空间性质。

注意{aj}L1{aj}L2\{a_j \} \in L_1 \Longrightarrow \{a_j \} \in L_2

线性序列的自协方差函数收敛性

{aj}l1\{a_j\} \in l_1
k=γk<\sum_{k=-\infty}^{\infty} |\gamma_k| < \infty
事实上,
k=γkσ2k=j=ajaj+k=σ2j=ajk=aj+k=σ2(j=aj)2<\begin{aligned} & \sum_{k=-\infty}^{\infty} |\gamma_k| \leq \sigma^2 \sum_{k=-\infty}^{\infty} \sum_{j=-\infty}^{\infty} |a_j a_{j+k}| \\ & = \sigma^2 \sum_{j=-\infty}^{\infty} |a_j| \sum_{k=-\infty}^{\infty} |a_{j+k}| = \sigma^2 (\sum_{j=-\infty}^{\infty}|a_j|)^2 < \infty \end{aligned}

○○○○○○

线性序列的应用:

线性序列描述了自协方差函数衰减到零的时间序列。 只要样本自协方差函数衰减到零就可以用线性序列来描述。

单边线性序列:

Xt=j=0ajεtj, tZX_t = \sum_{j={\mathbf 0}}^\infty a_j \varepsilon_{t-j}, \ t \in \mathbb Z
称为单边运动平均(MA),或单边无穷滑动和。 这样的X_t有因果性:XtX_t只受sts \leq tεs\varepsilon_s影响而不受 tt时刻以后的εs\varepsilon_s影响。
γk={j=0ajaj+k, k0γk k<0\gamma_k = \begin{cases} \sum_{j=0}^\infty a_j a_{j+k},\ & k \geq 0 \\ \gamma_{-k} \ & k < 0 \end{cases}

相关文章: