第三讲：随机抖动的数学模型——高斯分布、热噪声与散粒噪声
说到随机抖动（Random Jitter, RJ），好多新手第一反应就是“随机嘛，就是碰运气”。但做信号完整性可不能靠猜。RJ 的根源在物理层——热噪声和散粒噪声，这两个都是高斯分布的嫡系。今天我们就把它扒开看看。
我刚开始接触抖动分解时，也觉得 RJ 不就一滩噪声吗？直到有一次调试 10Gbps SerDes，眼图死活关不拢，用频谱仪一看，RJ 的尾巴把眼图压得快闭上了。那时候我才意识到：不搞懂 RJ 的数学底子，你连分离 RJ/DJ 的门都摸不着。
3.1 RJ 的数学本质：高斯分布
随机抖动的瞬时值服从高斯（正态）分布。说白了，它的概率密度函数就是那个铃铛曲线：
p(x) = (1 / (σ√(2π))) * exp(-x² / (2σ²))这里的 σ 就是标准偏差，单位是秒或 UI。RJ 的均值为 0，因为它是零均值随机过程。
为什么 RJ 对系统这么致命？关键是它的尾巴——高斯分布的尾巴延伸到无穷远。你想想看，哪怕 σ 只有 1ps，理论上也可能出现 10σ 的抖动事件，虽然概率极小，但数字系统里几亿次跳变，总有一次会撞上。
我个人的习惯是：在做链路预算时，RJ 的峰峰值不是直接取 ±3σ，而是根据误码率要求（比如 BER=10⁻¹²）反推对应的 σ 倍数。一般 BER=10⁻¹² 对应约 ±7σ，而不是 ±3σ。很多人忽略这点，结果设计余量不足。
3.2 热噪声：电阻里的随机骚动
热噪声（Johnson-Nyquist noise）是导体内部载流子的无规则热运动产生的。温度越高，振动越剧烈，噪声功率越大。它的功率谱密度是平坦的（白噪声）：
S_v(f) = 4k_B T R (V²/Hz)k_B 是玻尔兹曼常数（1.38e-23 J/K），T 是温度（K），R 是电阻（Ω）。
把热噪声耦合到时钟或数据通道上，就变成了时间域的随机抖动。怎么算？有个常用公式：
σ_RJ (热噪声) = (V_n_rms) / (2πf_c · A)其中 V_n_rms 是噪声电压均方根，f_c 是信号频率，A 是信号摆幅。
3.3 散粒噪声：电流中的离散颗粒感
散粒噪声（Shot noise）源于电流的不连续性——电子是一颗一颗通过的。当电流流过势垒（比如 PN 结、肖特基二极管）时，到达时间的随机性就产生散粒噪声。它的功率谱密度也是平坦的：
S_i(f) = 2qI_dc (A²/Hz)q 是电子电荷（1.6e-19 C），I_dc 是直流电流（A）。
散粒噪声转换成时间抖动的公式与热噪声类似，但噪声源是电流噪声：
σ_RJ (散粒噪声) = (i_n_rms) / (2πf_c · A)注意：散粒噪声只在有直流偏置的器件中才显著。无源电阻只产生热噪声，不产生散粒噪声。
3.4 两者的区别与联系
| 特性 | 热噪声 | 散粒噪声 |
|---|---|---|
| 根本原因 | 载流子热运动 | 电流的离散电荷效应 |
| 存在条件 | 任何带电阻的导体 | 有势垒的器件（PN结、肖特基） |
| 功率谱密度 | 4k_B T R | 2qI_dc |
| 频率特性 | 白噪声（上限 GHz） | 白噪声（上限极高） |
| 温度依赖性 | 强烈（正比于 T） | 弱（仅通过 I_dc 间接影响） |
| 时间域分布 | 高斯 | 高斯（大电流下近似） |
你可能会问：两种噪声叠加后还是高斯吗？答案是肯定的。两个独立高斯随机变量之和仍是高斯，方差相加。所以整体 RJ 仍可用单一 σ 描述。
3.5 知识体系总览
下面这张图把 RJ 的两个物理源、数学描述和它们在抖动分解中的位置串了起来。我建议你把它存下来，后面讲 DJ 分离时还会反复引用。
3.6 从数学到工程：如何测量 RJ
实际测量 RJ 时，示波器的直方图功能会把 TIE（时间间隔误差）统计出来。理想的高斯直方图应该是完全对称的钟形曲线。但现实总有偏差，因为 DJ 的尾巴会叠加上去。
我习惯用双狄拉克方法（Dual-Dirac）分离 RJ 和 DJ。简单说：从直方图的左右两边分别拟合高斯尾巴，得到左右两个 σ。如果左右 σ 相等，那就是纯高斯 RJ。一旦不等，就说明有 DJ 成分。具体操作后面章节会细讲。
RJ_pp(BER) = 2 · Q_factor(BER) · σ_RJ
其中 Q(10⁻¹²) ≈ 7.03，Q(10⁻¹⁵) ≈ 7.94。
别死记，查表或者用 math.erfinv 算一下就好。
3.7 高频 I/O 中的 RJ 实战案例
几年前我做 28Gbps NRZ 收发器。实验室测试发现，当温度从 25°C 升到 85°C 时，RJ 几乎翻倍。一开始我怀疑是电源噪声，排查无果。后来算了一下热噪声随温度变化：
σ_RJ (85°C) / σ_RJ (25°C) = √((273+85)/(273+25)) ≈ √(358/298) ≈ 1.097也就 10% 的增幅，不可能翻倍。最后发现是 VCO 的变容管在高温下漏电流剧增，产生了额外的散粒噪声。那个变容管的偏置电路设计有缺陷，高温下 PN 结反向漏电从 nA 级升到了 μA 级，散粒噪声变成原来的几十倍。换了低漏电工艺后问题解决。
所以 RJ 排查不能只看温度系数，要同时考虑热噪声和散粒噪声，尤其注意有源器件和势垒结构。
3.8 小结一下
RJ 的两个物理来源——热噪声和散粒噪声——都是高斯分布的白噪声。它们在时间域叠加后仍是高斯，方差相加。工程上通过测量 TIE 直方图的尾巴来提取 σ，再用 BER 对应的 Q 因子换算峰峰值。记住：RJ 不可消除，只能通过低噪声设计尽量降低。真正麻烦的是 DJ，因为 DJ 有确定性的结构，分离起来更讲究。
嗯，这一讲就到这里。下回我们聊聊确定性抖动 DJ 的种类和建模方法。