4. 红外系统噪声分析：1/f噪声、约翰逊噪声、散粒噪声、背景限噪声与NETD
噪声这东西，说白了就是限制我们红外系统性能的“天花板”。你光学再好，探测器再贵，噪声一上来，图像照样糊成一片。我做了这么多年红外系统，最深的体会就是：搞懂噪声，你就搞懂了系统极限的一半。
这一章我们重点拆解红外系统中常见的几种噪声，以及那个绕不开的核心指标——NETD（噪声等效温差）。
4.1 红外系统的噪声源分类
红外系统的噪声来源很多，但归纳起来，主要就这几类：
| 噪声类型 | 物理来源 | 频率特性 | 典型量级 |
|---|---|---|---|
| 1/f 噪声 | 材料缺陷、接触不良 | 低频主导（1/f^α） | 低频段显著 |
| 约翰逊噪声 | 载流子热运动 | 白噪声（平坦） | 与电阻、温度相关 |
| 散粒噪声 | 光生载流子的随机性 | 白噪声 | 与光电流平方根成正比 |
| 背景限噪声 | 背景光子通量涨落 | 白噪声 | 系统理论极限 |
| 读出噪声 | ROIC 电路引入 | 与电路设计相关 | 通常固定值 |
嗯，这里要注意：实际系统中这些噪声是叠加在一起的，你测到的总噪声是它们的均方根和。我见过不少新手，拿着频谱分析仪看到低频有抬升，就以为是电路问题，其实很可能就是 1/f 噪声在作怪。
4.2 1/f 噪声（闪烁噪声）
1/f 噪声，也叫闪烁噪声。它的功率谱密度与频率成反比，频率越低，噪声越大。
数学表达式很简单：
S(f) = K / f^α其中 K 是与材料和工艺相关的常数，α 通常在 0.8~1.2 之间，理想情况下取 1。
关键点：1/f 噪声在低频段（通常 < 1 kHz）占主导。对于凝视型红外焦平面阵列，帧频一般在 30~60 Hz，正好落在 1/f 噪声的“重灾区”。
我个人习惯在系统设计时，会采用相关双采样（CDS）或斩波稳定技术来抑制 1/f 噪声。我曾经在一个项目中，探测器低频噪声特别大，图像上全是“闪烁”的条纹。后来加了 CDS 电路，噪声直接降了一个数量级。
避坑指南：我曾经遇到过，1/f 噪声的 α 值偏离 1 太多，结果用标准滤波器怎么都滤不干净。后来发现是探测器封装接触不良导致的。所以，先排除工艺和接触问题，再谈电路设计。
4.3 约翰逊噪声（热噪声）
约翰逊噪声，也叫热噪声。它是导体内部载流子随机热运动产生的。只要温度高于绝对零度，它就存在。
它的均方根电压为：
V_j = √(4 k T R Δf)其中：
- k：玻尔兹曼常数（1.38×10⁻²³ J/K）
- T：绝对温度（K）
- R：电阻值（Ω）
- Δf：系统带宽（Hz）
你想想看，这个噪声跟频率没关系，是典型的白噪声。它没法通过滤波完全消除，只能通过降低温度或减小电阻来抑制。
我做过一个制冷型红外系统，探测器工作在 77 K 液氮温度下。这时候约翰逊噪声已经很小了，但读出电路的输入级电阻如果太大，照样会引入可观的热噪声。所以，前级放大器的设计非常关键。
4.4 散粒噪声
散粒噪声来源于光生载流子的随机产生过程。光子到达探测器是随机的，服从泊松分布。
它的均方根电流为：
I_shot = √(2 q I_dc Δf)其中：
- q：电子电荷（1.6×10⁻¹⁹ C）
- I_dc：平均光电流（A）
- Δf：系统带宽（Hz）
注意：散粒噪声也是白噪声。它和约翰逊噪声的区别在于：散粒噪声与光电流有关，而约翰逊噪声与电阻和温度有关。
说白了，信号越强，散粒噪声越大。但信噪比 SNR 会随着信号增强而提高（因为信号是线性增长，噪声是平方根增长）。
4.5 背景限噪声（BLIP）
背景限噪声，全称 Background Limited Infrared Photodetection。这是红外系统的理论极限。
当探测器自身的噪声（1/f、约翰逊、散粒等）远小于背景光子涨落引入的噪声时，系统就进入了背景限工作状态。这时候，你再怎么优化探测器，性能也上不去了——因为天花板是背景光子决定的。
背景限噪声的表达式：
N_BLIP = √(2 η φ_B A_d Δf)其中：
- η：量子效率
- φ_B：背景光子通量密度（photon/(cm²·s)）
- A_d：像元面积（cm²）
- Δf：系统带宽（Hz）
重要概念：背景限探测率 D*_BLIP 是衡量探测器是否达到理论极限的标尺。如果实测 D* 接近 D*_BLIP，说明你的探测器已经“压榨”到极致了。
我记得有一次测试一个长波红外探测器，D* 值怎么都上不去。后来一算，背景限 D* 就在那个量级。嗯，不是探测器不行，是背景光子太“吵”了。
4.6 NETD（噪声等效温差）
NETD 是红外系统最核心的性能指标之一。它的定义是：当系统输出信噪比为 1 时，目标与背景之间的温差。
说白了，就是系统能“分辨”的最小温差。
NETD 的计算公式：
NETD = (V_noise) / (ΔV/ΔT)其中：
- V_noise：系统总噪声电压（均方根值）
- ΔV/ΔT：系统对温度的响应率（V/K）
更完整的表达式（考虑光学系统）：
NETD = (4 F² √Δf) / (τ₀ A_d D* (ΔP/ΔT))其中：
- F：光学系统 F 数
- τ₀：光学透过率
- A_d：像元面积
- D*：探测率
- ΔP/ΔT：辐射功率随温度的变化率
经验值：非制冷型红外系统的 NETD 通常在 30~100 mK，制冷型可以达到 10~20 mK，甚至更低。
4.7 NETD 的测试方法
我常用的 NETD 测试步骤是这样的：
- 准备黑体辐射源：设置两个温度点，比如 T₁ = 25°C，T₂ = 35°C。
- 采集图像数据：分别采集两个温度下的图像序列（至少 100 帧）。
- 计算响应率：ΔV/ΔT = (V₂ - V₁) / (T₂ - T₁)。
- 计算噪声：取同一像元在固定温度下的时域标准差。
- 计算 NETD：NETD = V_noise / (ΔV/ΔT)。
测试技巧：我曾经在测试时发现 NETD 值偏大，排查了半天，结果是黑体温度不稳定。所以，测试前一定要让黑体充分预热稳定，至少 30 分钟。
4.8 知识结构图
下面这张图帮你理清本章的知识脉络：
4.9 小结
这一章我们拆解了红外系统的五大噪声源，以及那个绕不开的 NETD。你记住一句话：噪声分析的本质，就是找到系统的短板。是 1/f 噪声限制了低频性能？还是背景光子限制了理论极限？找到它，你才知道往哪个方向优化。
我个人的经验是：先测 NETD，再拆噪声谱。NETD 告诉你系统“有多差”，噪声谱告诉你“差在哪里”。两者结合，问题基本就浮出水面了。