Skip to content

状态空间时间序列模型

State-space time-series model · State space model · 状态空间模型

用隐藏状态的 Markov 演化与给定状态的观测分布分解时间依赖的序列模型。

领域
统计学
条目类型
模型

形式陈述

状态空间时间序列模型用隐藏状态过程 (αt) 生成观测 (Yt)。给定初始分布 p(α0)、转移条件分布 p(αtαt1) 与观测分布 p(ytαt),并假设给定当前状态后,Yt 与其他状态和观测条件独立。联合分布因式分解为

p(α0:n,y1:n)=p(α0)t=1np(αtαt1)p(ytαt).

线性高斯形式写成

αt=Ttαt1+ct+Rtηt,Yt=Ztαt+dt+εt,

其中过程扰动 ηtN(0,Qt)、观测扰动 εtN(0,Ht) 跨期且彼此独立,α0N(a0,P0) 并与全部扰动独立。当 Qt,Ht 不随时间改变时,两列扰动是向量版白噪声过程;一般时变情形则是独立但异方差的高斯扰动序列,不能按恒定方差的白噪声定义称呼。一般状态空间模型可以非线性、非高斯,只有上述子类拥有标准 Kalman 递推的精确有限维高斯闭包。

直觉

状态是对预测未来而言足够的动态摘要。过去的全部影响先压入 αt,观测再以噪声方式读取它;于是长程观测依赖通过低维 Markov 状态传递。状态可以是不可直接观测的水平、趋势、季节分量或物理位置,而不必与 Yt 同维、同单位,也不必具有单独可解释性。

同一观测过程常有多种状态表示,增加不可观测坐标可保持似然不变。因此状态维数和参数需要可观测性、可控性或统计识别检查;“把所有历史塞进状态”虽形式上可行,却失去压缩和稳定计算的价值。状态空间是模型表示,不是某个滤波算法;过滤、预测和平滑是在该模型上求不同条件分布。

例子与边界

标量局部水平模型取

αt=αt1+ηt,Yt=αt+εt,

其中 α0N(10,4)Var(ηt)=1Var(εt)=4。因此先验下 Var(α1)=5Var(α2)=6,而

Var(Y1)=5+4=9,Var(Y2)=6+4=10,Cov(Y1,Y2)=Var(α1)=5.

共同的隐藏水平使观测相关,观测噪声只增加各自主对角方差。看到 Y1=12 后,状态的高斯更新会在先验均值 10 与观测 12 之间加权;具体增益和后验方差由Kalman 滤波器计算。

过程方差 Q=0 时,状态按已知递推演化,所有不确定性来自初值和测量;观测方差 H=0 时,若 Z 可逆则状态被无噪声读出。两者都为零可能导致奇异预测协方差,普通密度公式失效。短序列中较大的 Q 与较大的 H 可能产生相似粗糙度,却含义相反;模型错设或弱识别不能靠滤波递推本身修复。

推论与应用

过滤求 p(αty1:t),一步预测求 p(αt+1y1:t),平滑则利用未来观测求 p(αty1:n)。线性高斯情形由 Kalman 滤波与平滑精确递推,并可累积预测误差密度计算似然。非线性或非高斯模型可用扩展/无迹近似、网格或粒子方法,但这些算法引入各自近似条件,不改变原模型的条件独立结构。

ARMA、结构时间序列、动态回归和缺失观测可放进统一状态表示;这种表示等价需给出具体最小性、初始化与噪声假设,因此本文不在 metadata 中笼统声明 state-space 与 ARMA 等价。初始化尤其关键:稳定状态可用平稳协方差,随机游走水平需弥散或显式先验。错误地用固定小方差初始化,会让早期滤波过度自信并扭曲似然。

参考资料
  • James Durbin and Siem Jan Koopman, Time Series Analysis by State Space Methods, 2nd ed., Oxford University Press, 2012,Chs. 2–3,state-space formulation and linear Gaussian models。
  • Robert H. Shumway and David S. Stoffer, Time Series Analysis and Its Applications, 4th ed., Springer, 2017,Ch. 6,state-space models and dynamic linear models。
  • Andrew C. Harvey, Forecasting, Structural Time Series Models and the Kalman Filter, Cambridge University Press, 1989,Chs. 2–3,structural state-space modeling。
关系图谱11 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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

使用的工具