Skip to content

定理Theorem

Wold 分解定理

施工中 · 内容逐步修订

Wold decomposition theorem · Wold representation theorem

将弱平稳过程唯一拆为由创新驱动的单边线性部分与可由无限过去预测的确定部分。

形式陈述 ​

设 (Xt) 是零均值、有限方差的弱平稳过程。这里时间取遍整数,因而可以讨论无限久以前。在平方可积随机变量的 Hilbert 空间中,内积是 E(XY);令 Ht 为 {Xs:s≤t} 的闭线性张成,即有限线性组合及其均方极限构成的空间。PHtX 是Hilbert 空间投影定理保证存在的正交投影,表示用这些变量作最佳线性预测,H−∞=⋂tHt 则保留无论把观测截止时间推得多早都能线性恢复的信息。Wold 定理断言存在唯一正交分解

Xt=Vt+∑j=0∞ψjεt−j,ψ0=1,∑j=0∞ψj2<∞,

其中 Vt∈H−∞ 是完全可预测的平稳部分;(εt) 是方差 σε2 的白噪声过程,并且 εt=Xt−PHt−1Xt 是相对于整个线性过去的一步创新。两部分在所有时点正交。

若 H−∞={0},过程称纯非确定,此时 Vt=0,得到单边无限 MA 表示。一步预测误差方差 σε2>0 时,创新由正交投影确定,归一化 ψ0=1 后系数与创新方差唯一;若该方差为零,则过程全属确定部分,零创新前的系数不再具有可识别含义。定理不声称创新独立或高斯,也不声称系数绝对可和;平方可和足以保证随机级数在均方意义收敛。

直觉

把截至 t−1 的所有观测当作一个线性信息空间,先做最佳线性一步预测。无法被这个空间解释的残差就是本期真正新增的线性信息;不同期创新彼此正交。反复把过去值也拆成“预测加创新”,就得到当前值对历次新信息的冲击响应。无限久以前仍能精确携带的部分,则被留在 Vt 中。

构造的关键是把互相正交的创新当作预测信息的坐标。当 σε2>0 时,第 j 个系数由投影给出 ψj=E(Xtεt−j)/σε2。平稳性使它只依赖滞后 j;从 Xt 中扣除所有这些投影后,剩下的就是远古信息部分。这两项极限可以由嵌套投影具体完成。因为

Ht=Ht−1⊕span{εt},

逐步正交分解给出,对每个 m≥1,

Xt=PHt−mXt+∑j=0m−1ψjεt−j.

互相正交使 σε2∑j<mψj2≤EXt2,所以创新级数在均方意义下收敛。另一方面,记 Pm=PHt−m,当 m≥l 时,嵌套子空间给出

‖PlXt−PmXt‖22=‖PlXt‖22−‖PmXt‖22.

右侧来自一个递减非负数列,故 PmXt 是 Cauchy 列。其极限属于每个闭空间 Ht−l,也就属于交集;且 Xt 减去该极限仍与交集正交。因此极限恰为 Vt=PH−∞Xt。平稳性使移位保持内积,并与这些投影相容,从而 Vt 弱平稳、系数只依赖滞后。若创新方差为零,则所有相邻过去空间已经相同,直接得到 Xt=Vt。

定理说明线性创新表示非常普遍,却不是说所有平稳序列“本质上都是某个有限 ARMA”。ψj 可无限多、衰减很慢,创新可非独立,确定部分也可能非零。Wold 是二阶投影定理;若关心非线性条件期望、尾部依赖或完整条件分布,需要更强结构。

例子与边界

取可逆 MA(1)

Xt=εt+12εt−1,Var(εt)=4.

由可逆级数

εt−1=Xt−1−12Xt−2+14Xt−3−⋯,

可知 12εt−1 属于 Ht−1,而 εt 与过去正交。因此一步预测为 X^t|t−1=12εt−1,预测误差恰为 εt、方差为 4;Wold 系数是 (1,1/2,0,…)。

同一协方差 γ(0)=5,γ(1)=2 也可写成不可逆形式 Xt=ut+2ut−1、Var(ut)=1。这里的 ut 不能由观测过去作为一步创新恢复;Wold 规范选择前一组方差 4、系数 1/2 的可逆表示。这具体区分了“方程里写的白噪声”和“由过程过去定义的创新”。

若 A,B 零均值、方差同为 s2<∞ 且互不相关,并令 Xt=Acos⁡(λt)+Bsin⁡(λt),则 γ(h)=s2cos⁡(λh),所以过程弱平稳。当 sin⁡λ≠0 时,任意两个相邻值就能线性确定 A,B 及整条序列。其一步创新方差为零,全部落入 Vt;强行要求正方差白噪声 MA 表示会失败。把这一随机正弦与独立白噪声相加,就同时具有确定项和创新项;例如谱中的离散谱线与正的常数谱密度分别对应这两部分。一般情形不能仅凭“连续谱”三个字断定纯非确定性,还需相应的谱密度条件。

推论与应用

Wold 表示为线性预测提供统一基础。在时点 t 预测 Xt+h,保留确定项和 j≥h 的已知创新项;误差为 ∑j=0h−1ψjεt+h−j,方差为 σε2∑j=0h−1ψj2。因此上面的 MA(1) 一步预测误差方差是 4,两步及更远均为 4(1+1/4)=5:两步之后已经没有能由当前过去解释的随机部分。因果且可逆 ARMA 的 MA(∞) 展开是这一创新表示的参数化特例;ARMA 还施加有理传递函数和更强衰减结构,因此可用有限参数估计。

在经验分析中,残差若要解释成创新,至少应接近与完整已用过去正交,而不只是与一两个滞后不相关。可逆性保证模型冲击与 Wold 创新对齐;非可逆参数化虽可复制同一二阶分布,却赋予冲击错误的时间方向。定理也不保证从有限样本稳定估计无限组 ψj,截断阶数仍需依据预测误差和正则性选择。

参考资料
  • Anna Mikusheva(授课)、Paul Schrimpf(记录),Lecture 10: Introduction to VARs,MIT 14.384,2007 年讲义、2012 年修订,pp. 1–2,Theorem 1:正交创新构造与平方可和性。

  • Peter J. Brockwell and Richard A. Davis, Time Series: Theory and Methods, 2nd ed., Springer, 1991,Theorem 5.7.1,Wold decomposition of stationary processes。

  • Helmut Lütkepohl, New Introduction to Multiple Time Series Analysis, Springer, 2005,§2.1.2,Wold representation and innovations。

  • James D. Hamilton, Time Series Analysis, Princeton University Press, 1994,§§4.1 and 4.8,linear projections, innovations, and Wold decomposition。

关系图谱12 个相邻概念 · 1 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组