每日精学 015｜白噪声周期图的分布与估计方差

2026-09-25

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

前一课程给出了白噪声频点功率的期望值，但期望正确并不意味着一次观测稳定。单帧周期图常呈现大幅锯齿，原因不是 FFT 数值误差，而是功率估计量本身具有与均值同量级的随机标准差。

谱估计需要同时讨论偏差和方差。延长记录可改善频率分辨率并减小某些偏差，却不会自动降低单个白噪声频点的归一化方差。降低方差需要独立记录平均、分段平均或其他平滑。

01
模型与符号

X[k]：未归一化 DFT 频点（样本单位）

I[k]：单帧周期图（PSD 单位）

S_x：真实白噪声功率谱密度（PSD 单位）

Y：归一化周期图 I/S_x（无量纲）

L：独立平均记录数（正整数）

ν：功率估计的等效卡方自由度（自由度）

c_v：变异系数，标准差除以均值（无量纲）

02
单频点指数分布

对复高斯白噪声，DFT 频点可写为两个独立高斯分量。适当归一化后

$Y=\frac{I[k]}{S_x}=\frac{U^2+V^2}{2}, \qquad U,V\overset{\mathrm{iid}}{\sim}\mathcal{N}(0,1).$

因此

$2Y\sim\chi_2^2$

并有

$p_Y(y)=e^{-y},\quad y\geq0,\qquad \operatorname{E}\{Y\}=1,\quad \operatorname{var}(Y)=1.$

单频点标准差等于真实谱值。

03
记录长度与非一致性

若白噪声记录长度从 N 增加到更大值，DFT 方差和归一化分母按相同比例增长，I[k]/S_x 的指数分布保持不变。于是

$\lim_{N\rightarrow\infty}\operatorname{var}\!\left(\frac{I_N[k]}{S_x}\right)=1,$

单帧周期图在逐点意义下不是一致估计量。

记录变长的主要收益是频率网格变细和谱窗主瓣变窄。

04
独立记录平均

对 L 个独立周期图求平均

$\overline I_L[k]=\frac{1}{L}\sum_{\ell=1}^{L}I_\ell[k],$

则归一化变量满足

$\frac{\overline I_L[k]}{S_x}\sim \operatorname{Gamma}\!\left(L,\frac{1}{L}\right), \quad \operatorname{var}\!\left(\frac{\overline I_L}{S_x}\right)=\frac{1}{L}.$

等效自由度为 ν=2L，变异系数为

$c_v=1/\sqrt{L}$

05
物理含义

周期图的尖峰和深谷是随机功率估计的自然表现。以分贝显示时，指数分布的非对称性会使视觉平均与线性功率平均不同，因此应在线性功率域完成平均后再取对数。独立平均将自由度从 2 提高到 2L，但需要更多时间记录或牺牲同一记录内的频率分辨率。

06
固定参数算例

数值验证生成归一化指数功率样本，并用复高斯白噪声 FFT 频点作独立交叉检查。比较 L=1,4,16 时的经验均值、方差、变异系数和解析 Gamma 密度；另以不同 N 的单帧试验确认归一化方差不会随记录长度消失。

独立双精度参考计算结果

量	数值
normalized single bin mean	1.00052467351
normalized single bin variance	1.00122532938
cv L1	1.00008775757
cv L4	0.500233516114
cv L16	0.248830217708
theoretical cv L4	0.5
theoretical cv L16	0.25
样本数	240000

图 1：单频点指数分布独立平均 Gamma 分布变异系数与记录长度比较（子图 1）

图 1：单频点指数分布独立平均 Gamma 分布变异系数与记录长度比较（子图 2）

图 1　单频点指数分布独立平均 Gamma 分布变异系数与记录长度比较

图表保留原报告的独立双精度参考结果。本次补充运行同参数 MATLAB 脚本并通过数值断言；随机算例比较统计关系，参考程序与 MATLAB 的随机样本不作逐点相同的假设。

07
数字实现的边界

模块 power_average8.sv 连续接收八个无符号已标定功率值，使用扩位累加器求和，并在第八个样本到达时加 4 后右移 3 位，实现最近整数平均。

模块输出一个周期有效脉冲。testbench 检查全零、常量、递增序列、最大输入和两组连续帧，确认累加器清零与输出时序。

RTL 只实现本课中可明确量化的确定性数据通路或控制功能。概率分布拟合、浮点函数求值和图形生成留在 MATLAB 参考程序中，防止把软件统计功能误认为可综合硬件。

08
适用范围

独立周期图平均适用于重复脉冲、稳定测量时段和多通道独立快拍。平均前必须保证各记录采用相同窗函数、FFT 缩放和 PSD 单位。若记录之间相关，L 不能直接解释为独立平均次数。

09
工程上需要注意

实信号的直流与 Nyquist 频点只有一个实自由度，其分布不同于一般复频点。
有色谱在谱窗主瓣内变化时，周期图均值还包含频域平滑偏差。
重叠或滤波会使相邻记录相关，实际自由度小于 2L。
分贝域直接平均会产生偏差，必须先在线性功率域平均。

白噪声单频点周期图在归一化后服从指数分布，均值正确而相对方差等于 1。增加单帧长度不会消除这种逐点起伏。平均 L 个独立记录可把方差降为 1/L，等效自由度增至 2L，这构成后续 Welch 分段平均的统计基础。

内容依据：每日精学第 15 课《白噪声周期图的分布与估计方差》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。

阅读11

写评论...

每日精学 015｜白噪声周期图的分布与估计方差

01模型与符号

02单频点指数分布

03记录长度与非一致性

04独立记录平均

05物理含义

06固定参数算例

07数字实现的边界

08适用范围

09工程上需要注意