形式陈述
设弱平稳过程的自协方差函数 公理库 自协方差函数 Autocovariance function · ACVF 以滞后为自变量记录弱平稳过程两个时点之间协方差的函数。 为 γ ,且 0 < γ ( 0 ) < ∞ 。自相关函数定义为
ρ ( h ) = Corr ( X t + h , X t ) = γ ( h ) γ ( 0 ) , h ∈ Z . 因此 ρ ( 0 ) = 1 、ρ ( − h ) = ρ ( h ) 且 | ρ ( h ) | ≤ 1 。它继承非负定性:任意有限矩阵 [ ρ ( t i − t j ) ] i , j 必须半正定,所以逐点落在 [ − 1 , 1 ] 只是必要而非充分条件。若过程乘以非零常数或整体加常数,ρ 不变;若 γ ( 0 ) = 0 ,过程退化为常数,相关系数因除以零而未定义。
样本 ACF 常写成
ρ ^ ( h ) = ∑ t = 1 n − h ( x t + h − x ¯ ) ( x t − x ¯ ) ∑ t = 1 n ( x t − x ¯ ) 2 , 0 ≤ h < n . 这一定义使用同一个分母,保证 ρ ^ ( 0 ) = 1 ;它是随机估计量,并非总体 ρ ( h ) ,且因估计均值、可用配对数随滞后缩短以及分母随机而在有限样本中通常有偏。对独立同分布并满足常规矩条件的白噪声,固定正滞后的渐近标准误约为 n − 1 / 2 ;但本文定义的一般白噪声只保证不相关,对拟合模型的残差、许多滞后同时查看或非白噪声过程,简单的 ± 1.96 / n 都不能原样解释。
直觉
自相关把“共同变化多少”换成“相对于序列自身波动有多大”。ρ ( 1 ) = 0.8 表示相邻时点具有强正线性关联,却不表示下一期有 80% 概率等于本期,也不表示 80% 的方差被单一滞后解释。符号给同向或反向倾向,绝对值给线性强度;非线性依赖和条件异方差可能在所有非零滞后 ACF 都近零时仍然很强。
ACF 的形状常比单个数更有信息。有限阶移动平均的总体 ACF 在某个滞后后精确截尾;因果自回归的 ACF 通常以指数或阻尼振荡方式拖尾;单位根水平序列的样本 ACF 则常在许多滞后缓慢下降。但这些是模型蕴含的诊断图样,不是仅凭一张相关图就能唯一识别模型的定理。季节峰也可能来自确定性季节均值未去除,而非随机季节动力学。
例子与边界
沿用 X t = ε t + 1 2 ε t − 1 且 Var ( ε t ) = 4 的可复算模型,其自协方差为 γ ( 0 ) = 5 、γ ( 1 ) = 2 、γ ( h ) = 0 (| h | ≥ 2 ),所以
ρ ( 0 ) = 1 , ρ ( ± 1 ) = 2 5 = 0.4 , ρ ( h ) = 0 ( | h | ≥ 2 ) . 创新方差若从 4 改成 9 ,γ ( 0 ) 与 γ ( 1 ) 同乘 9 / 4 ,ACF 仍完全相同,展示其无量纲性。三时点相关矩阵为
R 3 = ( 1 0.4 0 0.4 1 0.4 0 0.4 1 ) , 其顺序主子式为 1 , 0.84 , 0.68 ,均为正,故这段候选 ACF 与半正定约束相容。
零 ACF 不能推出独立。若 Z t 独立标准正态并令 X t = Z t Z t − 1 ,则对 h ≠ 0 有 E ( X t X t − h ) = 0 ,总体 ACF 与白噪声相同;但 X t 2 与 X t − 1 2 共享 Z t − 1 2 ,且 E ( X t 2 X t − 1 2 ) = 3 ≠ 1 = E X t 2 E X t − 1 2 ,所以相邻变量依赖。另一个边界是趋势:x t = t + e t 的高样本 ACF 主要反映共同趋势,不能解释为平稳随机记忆。
推论与应用
ACF 是线性预测、模型识别和残差诊断的共同接口。与偏自相关函数 公理库 偏自相关函数 Partial autocorrelation function · PACF 在剔除中间滞后最佳线性解释后衡量两个相隔时点剩余相关的函数。 配合时,AR(p ) 的 PACF 在 p 阶后截尾,而 MA(q ) 的 ACF 在 q 阶后截尾;有限样本只会近似呈现这一差别。对拟合后的残差,Ljung–Box 一类端口检验把多个样本自相关合并,但自由度应扣除已估参数,拒绝只说明仍有线性依赖,不自动指出应增加哪一阶。
对平稳 AR(1) X t = ϕ X t − 1 + ε t ,| ϕ | < 1 时 ρ ( h ) = ϕ | h | 。于是 ϕ = 0.8 的一阶相关为 0.8 、五阶为 0.8 5 = 0.32768 ;ϕ = − 0.8 则奇数滞后为负、偶数为正。靠近 1 的指数衰减很慢,却仍不同于单位根。做置信区间与模型比较时应把参数估计误差、数据选择和多重滞后检查纳入,而不是把越过单条参考线当成确定结构。
参考资料
Peter J. Brockwell and Richard A. Davis, Introduction to Time Series and Forecasting , 3rd ed., Springer, 2016,§§1.4.1 and 2.4.2,sample ACF and its sampling properties。
George E. P. Box, Gwilym M. Jenkins, Gregory C. Reinsel, and Greta M. Ljung, Time Series Analysis: Forecasting and Control , 5th ed., Wiley, 2015,Chs. 2, 6, and 8,correlation structure, identification, and diagnostics。
Robert H. Shumway and David S. Stoffer, Time Series Analysis and Its Applications , 4th ed., Springer, 2017,§1.3 and Ch. 3,ACF and ARMA identification。