每日精学 157｜相关样本均值的有效样本数

2026-09-25

雷达 · 通信 · 电子战　/　基础精学 157

相关数据的算术均值不再具有独立样本的方差缩减速度。有效样本数将全部样本间协方差折算为等效独立观测数量。

实际记录只包含有限个采样值。由记录计算得到的数值与总体参数并不相同：前者受记录长度、噪声、边界约定及运算精度影响，后者属于给定模型的属性。建立可复核结论需要先固定输入含义、统计假设和归一化方式，再比较可计算量与理论量。

01
模型与符号

实值宽平稳过程具有有限方差；样本均值估计总体常均值；实验使用平稳AR(1)过程。

N：实际记录长度（样本）

ρ[k]：归一化自协方差（无量纲）

N_eff：有效样本数（样本）

02
从模型到公式

均值定义

对全部样本等权平均。

$\bar x=N^{-1}\sum_{n=0}^{N-1}x[n]$

方差展开

保留不同样本之间的协方差，不预先假设独立。

$\operatorname{var}(\bar x)=N^{-2}\sum_{i=0}^{N-1}\sum_{j=0}^{N-1}R_x[i-j]$

按滞后计数

滞后正负成对，各有相同有效样本对数。

$\operatorname{var}(\bar x)=\frac{\sigma^2}{N}\left[1+2\sum_{k=1}^{N-1}(1-k/N)\rho[k]\right]$

等效样本数

以独立观测的均值方差定义等效数量。

$N_{\mathrm{eff}}=\frac{N}{1+2\sum_{k=1}^{N-1}(1-k/N)\rho[k]}$

AR相关代入

一阶正相关模型形成几何序列。

$\rho[k]=a^{|k|},\quad|a|<1$

长记录近似

忽略边界项后由几何级数求和得到近似。

$N_{\mathrm{eff}}\simeq N\frac{1-a}{1+a}\quad(N\text{ 足够大})$

负相关也会改变有效样本数

对方差为 σ² 的平稳序列，样本均值方差由所有协方差项决定：

$\operatorname{var}(\bar x)=\frac{\sigma^2}{N} \left[1+2\sum_{k=1}^{N-1}\left(1-\frac{k}{N}\right)\rho[k]\right].$

括号内不必大于一。

若相邻样本倾向于相互抵消，有效样本数可以大于记录长度；它是与均值方差相匹配的等效量，不是实际获得了更多独立记录。对稳定 AR 一阶相关 ρ[k]=a^k，长记录近似为

$N_{\mathrm{eff}}\approx N\frac{1-a}{1+a},\qquad |a|<1.$

该近似需要记录足够覆盖相关衰减。

当 a 极接近正负一时，有限长度权重和边界项不可忽略；直接代入渐近公式可能严重误估均值的不确定性。

03
物理含义

强正相关意味着相邻样本携带相近的信息，实际记录长度不能直接视作独立证据数量。负相关可能使均值噪声抵消，从而有效样本数大于原始数量；因此该量不是物理存储容量，而是针对均值方差定义的统计量。

04
固定参数算例

取长度128、相关系数0、0.5、0.9的平稳AR序列，各重复4000次，以精确有限和预测均值方差并比较有效样本数。数值脚本使用固定随机种子；随机实验的判据采用有限容差，代数关系的判据采用浮点残差。图中同时展示参数变化趋势与对应误差，关键数值写入独立结果文件。

固定参数数值验证结果

量	数值
记录长度	128
repetitions	4000
independent variance	0.00807919715923
a09 mean variance	0.137487661328
a09 theory variance	0.137451187147
a09 effective count	7.27531002647
a05 effective count	43.1157894737

图 1　相关样本均值的有效样本数的固定参数数值验证

表中数值和图中曲线由 MATLAB 实际运行得到，随机种子为 20261076；运行版本与完整结果另存为独立文件。

05
数字实现的边界

SystemVerilog实现输入样本及其外部延迟副本的相关乘积、平方与样本计数累加。端口为16位有符号码，内部64位，同步复位，valid空拍保持状态；

按外部首尾标志进行块累加并在块尾输出。自检覆盖多组数据、空拍、边界和复位，仅验证所述整数子核，不代表完整算法的硬件验证。

独立自检程序已在 Icarus Verilog 13.0 中实际运行并通过。该结果不包含器件综合、静态时序分析或板级测试。

06
适用范围

用于低通噪声平均、缓慢漂移测量及置信区间中独立性假设的审查。实际估计相关长度时应限制高滞后噪声，以免有效样本数被不稳定尾部主导。

07
工程上需要注意

非平稳漂移不适合用固定自协方差描述。

不同统计量拥有不同等效自由度，均值的有效样本数不能直接用于峰值概率。

有限精度实现必须区分输入量化、内部舍入和最终结果格式，不能仅依据浮点曲线估计硬件误差。

均值方差由完整协方差和决定。独立采样只是相关结构为零的特殊情形。

内容依据：每日精学第 157 课《相关样本均值的有效样本数》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。

阅读10

写评论...

每日精学 157｜相关样本均值的有效样本数

01模型与符号

02从模型到公式