第二章
2.1 一般线性过程
一、线性平稳序列
线性平稳序列是白噪声的线性组合得到的序列。 最简单的线性平稳序列是有限运动平均。
设{εt}={εt: t∈Z}是WN(0,σ2).对于非负整数q和常数a0,a1,…,aq(a0=0,aq=0),称
Xt=j=0∑qajεt−j=a0εt+a1εt−1+⋯+aqεt−q, t∈Z
是白噪声{εt}的(有限)运动平均或者滑动平均,简称MA(Moving Average)
注:可将{Xt}视为线性滤波器的输出,白噪声看成驱动系统的扰动序列(激励).
MA的平稳性:
EXt=0
E(Xt+kXt)={σ2∑j=0q−kajaj+k,0,0≤k≤q,k>q
可见{Xt}平稳。γk=0,∀k>q,称这样的序列为q相关的。

线性滤波器的原始数据与滤波结果是一种算术运算,即用加减乘除等运算实现。包括 方框滤波(boxFilter)、均值滤波(blur)、高斯滤波(GaussianBlur);
非线性滤波器的原始数据与滤波结果是一种逻辑关系,即用逻辑运算实现。包括最大值滤波、中值滤波(medianBlur)和双边滤波(bilateralFilter)。
参考资料:
常见线性滤波
线性滤波和非线性滤波
二、线性过程的两种等价形式
1、 传递形式
时间序列分析中建立随机模型的思想:
将顺序值之间高度依赖的时间序列{Xt}看成由一系列独立“冲击”序列所生成.
通常将“冲击”序列理想化为白噪声过程{εt, t∈Z},时间序列Xt取为白噪声的加权和.
Wold分解式(正交分解):
任意零均值纯非确定的平稳过程都可表示为线性形式(无穷阶滑动平均):
Xt=j=0∑∞φjεt−j,t∈Z
称为沃尔德系数(格林函数、传递函数),其中{εt, t∈Z}是白噪声序列。
注1:上式可表示为算子形式:Xt=φ(B)εt,其中φ(B)=1+∑j=1∞φjBj=∑j=0∞φjBj
称为传递函数或沃尔德系数的生成函数.
注2:对动态数据进行适当的预处理,可将非平稳序列 平稳化与零均值化.
注3:上式的工程解释:Xt=φ(B)εt (线性滤波器模型)
将时间序${X_t } $ 视为线性滤波器的输出,白噪声看成驱动系统的扰动序列(激励)

2、自回归形式
在适当条件下{Xt,t∈Z}可表示为线性形式
Xt=εt+j=1∑∞πjXt−j,t∈Z
或 εt=(1−∑j=1∞πjBj)Xt=π(B)Xt,t∈Z

理解:回归这一词的意思是从已有的数据发现规律。与过去的Xt有联系,即相关。
3、 两者之间的关系
Xt=φ(B)εt,t∈Z(∗) 与 εt=π(B)Xt,t∈Z(∗∗)
可以看出两者互为逆转形式。
将生成函数φ(B)作用于(∗∗)式
Xt=φ(B)εt=φ(B)π(B)Xt
得φ(B)π(B)=1
或φ(B)=π(B)−1,π(B)=φ(B)−1
传递形式:

逆转形式:

例:若平稳序列得自回归形式为Xt=πXt−1+εt或π(B)Xt=εt
其中 π(B)=1−πB,有:
φ(B)=π(B)−1=1−πB1=1+πB+π2B2+...+πjBj+...
求得传递形式为:
Xt=π−1(B)εt=φ(B)=∑j=0∞πjεt−j
三、线性过程的均值函数与自协方差函数
2.1.1 定理 (Levi 单调收敛定理)
如果非负随机变量序列{ξn}单调不减: 0≤ξ1≤ξ2≤⋯
则当ξn→ξ a.s.时, 有 limn→∞Eξn=Eξ
由于单调不减序列必有极限(极限也可以是正无穷),所以上面的随机变量ξ一定存在。但是可能会有p(ξ=∞)>0这时Eξ=∞,可以取值∞的随机变量通常被称为广义随机变量。
对于任何时间序列 {Yt}, 利用单调收敛定理得到
E[t=−∞∑∞∣Yt∣]==n→∞limE[t=−n∑n∣Yt∣]n→∞limt=−n∑nE∣Yt∣=t=−∞∑∞E∣Yt∣.
2.1.2定理 (控制收敛定理)
如果随机变量序列{ξn}满足∣ξn∣≤ξ0 a.s.和 E∣ξ0∣<∞,
则当ξn→ξ a.s.时,有 E∣ξ∣<∞,并且limn→∞Eξn→Eξ.
在上面两个定理成立条件下均有,
n→∞limEξn=En→∞limξn.
即期望与极限可以交换次序。
极限与积分交换次序的三个定理:
Lebesgue控制收敛定理是这样的描述的:
{fn}是可测集E上的可测函数列,F是可测函数列的控制函数并且可积,{fn}依测度收敛到f,则f在E上可积并且
limn→∞∫Efndx=∫Efdx
Levi定理是这样叙述的:
{fn}是可测集E上的可测函数列,并且是非负的单调递增的即0≤f1(x)≤f2(x)⋯fn(x)≤fn+1(x)≤⋯,那么{fn(x)}收敛(几乎处处收敛)于一个可积函数f(x)并且上面5楼的交换仍然成立。
limn→∞∫Efndx=∫Efdx
Fatou引理是这样说的:
{fn}是可测集E上的非负可测函数列, 有下面的非严格等式交换
∫Elimfn(x)dx≤lim∫Efn(x)dx
详情可见
积分的极限定理
2.1.3定理
设时间序列{Xt}满足suptE∣Xt∣<+∞,如果∑j=−∞∞∣ψj∣<+∞,则序列:
ψ(B)Xt=j=−∞∑∞ψjBjXt=j=−∞∑∞ψjXt−j (∗)
依概率1收敛(sup 上确界supremum)
进而,如果suptE∣Xt∣2<+∞,则序列均方收敛到同一极限。
证:由单调收敛定理和suptE∣Xt∣<+∞,有
E(∑j=−∞∞∣ψj∣∣Xt−j∣)=limn→∞E(j=−n∑n∣ψj∣∣Xt−j∣)≤limn→∞E(j=−n∑n∣ψj∣)suptE∣Xt−j∣<∞
因此∑j=−∞∞∣ψj∣∣Xt−j∣和ψ(B)Xt=∑j=−∞∞ψjBjXt=∑j=−∞∞ψjXt−j依概率1有限
如果suptE∣Xt∣2<+∞,n>m>0,则当m,n→∞时,有
E∣∣∣∑m<∣j∣<nψjXt−j∣∣∣2=≤∑m<∣j∣<n∑m<∣k∣<nψjψ^kXt−jX^t−ksuptE∣Xt∣2(∑m<∣j∣<n∣ψj∣)2→0
由Cauchy收敛准则,序列(∗)均方收敛。
令???? 表示序列(∗)的均方极限,由Fatou引理,有
E∣S−ψ(B)Xt∣2=≤Elimn→∞inf(s−∑j=−nnψjXt−j)infE(s−∑j=−nnψjXt−j)=0
即证明了S和ψ(B)Xt依概率1相等
推论2.1
如果随机变量序列{Yn}满足∑i=−∞∞E∣Yi∣≤∞,则 ∑i=−∞∞Yi,a.s.收敛, E∣∑i=−∞∞Yi∣<∞且
En=−∞∑∞Yn=n=−∞∑∞EYn
证明 令
η=∑i=−∞∞∣Yi∣
则η是随机变量, 由单调收敛定理知Eη<∞,η<∞,a.s.。
令 ξn=∑i=−nnYi
则 ∣ξn∣≤η,
记事件A={η<∞}, 则P(A)=1, 对ω∈A, 有∑i=−∞∞∣Yi(ω)∣<∞所以∑i=−∞∞Yi(ω)收敛, limn→∞ξn(ω)收敛,记为$\xi(\omega); 对\omega \notin A,定义\xi(\omega) = 0,则\xi是随机变量,且\lim_{n\to\infty} \xi_n = \xi, \text{ a.s.}$
记 ξ=∑i=−∞∞ξi.
由控制收敛定理 limn→∞Eξn=Eξ,
而 limn→∞Eξn=limn→∞E∑i=−nnYi=limn→∞∑i=−nnEYi=∑i=−∞∞EYi,
于是有E∑i=−∞∞ξi=∑i=−∞∞Eξi.
四、线性平稳序列
如果实数列{aj}满足 ∑j=−∞∞∣aj∣<∞,则称{aj}是绝对可和的. 记{aj}∈l1
注意:{aj}∈l1则{aj}∈l2(即∑jaj2<∞). 反之不一定成立。
对于绝对可和的实数列{aj}, 定义零均值白噪声{εt}的无穷滑动和如下
Xt=∑j=−∞∞ajεt−j, t∈Z.
则{Xt}是平稳序列。 EXt=0,
γk=σ2j=−∞∑∞ajaj+k,k∈Z.
令ck=∑j=−∞∞ajaj+k,k∈Z。
线性序列的a.s.收敛性
作为无穷和{Xt}有没有定义? 由Schwarz不等式,
E∣εt∣=E(∣εt∣⋅1)≤Eεt2⋅1=σ,
所以
∑j=−∞∞E∣ajεt−j∣=∑j=−∞∞∣aj∣E∣εt−j∣≤σ∑j=−∞∞∣aj∣<∞,
由推论2.1可知∑j=−∞∞ajεt−j a.s.收敛。
线性序列的L1收敛性
级数的余项的L1模
E∣∣∣∣∣∣∣j∣>N∑ajεt−j∣∣∣∣∣∣≤∣j∣>N∑∣aj∣E∣εt−j∣≤σ∣j∣>N∑∣aj∣→0
即线性平稳列在L1意义下收敛。
线性序列的平稳性
由推论2.1可知
Ej=−∞∑∞ajεt−j=j=−∞∑∞ajEεt−j=0
由Schwarz不等式知E∣ϵt−jϵt+k−l∣≤σ2, 于是
∑j=−∞∞∑l=−∞∞E∣ajalϵt−jϵt+k−l∣≤σ2∑j=−∞∞∑l=−∞∞∣aj∣∣al∣=σ2(∑j=−∞∞∣aj∣)2<∞,
由推论2.1,
EXtXt+k=Ej=−∞∑∞ajεt−jl=−∞∑∞alεt+k−l=j=−∞∑∞l=−∞∑∞ajalE(εt−jεt+k−l)=σ2j=−∞∑∞ajaj+k
即{Xt}是平稳序列, EXt=0,γk=σ2∑j=−∞∞ajaj+k.
线性序列的L2收敛性
设{aj}∈l2,即∑jaj2<∞,则
Xt=j=−∞∑∞ajεt−j (L2)
也是平稳序列。期望为零,自协方差函数同上。
Xt定义的无穷级数是L2收敛的。 证明需要应用Hilbert空间性质。
注意{aj}∈L1⟹{aj}∈L2。
线性序列的自协方差函数收敛性
当{aj}∈l1时
∑k=−∞∞∣γk∣<∞
事实上,
k=−∞∑∞∣γk∣≤σ2k=−∞∑∞j=−∞∑∞∣ajaj+k∣=σ2j=−∞∑∞∣aj∣k=−∞∑∞∣aj+k∣=σ2(j=−∞∑∞∣aj∣)2<∞
○○○○○○
线性序列的应用:
线性序列描述了自协方差函数衰减到零的时间序列。 只要样本自协方差函数衰减到零就可以用线性序列来描述。
单边线性序列:
Xt=j=0∑∞ajεt−j, t∈Z
称为单边运动平均(MA),或单边无穷滑动和。 这样的X_t有因果性:Xt只受s≤t的εs影响而不受 t时刻以后的εs影响。
γk={∑j=0∞ajaj+k, γ−k k≥0k<0