在看 LT1021 的资料的时候，我的爬虫给我带回来一篇论文，真的是太棒了！它第一次把离我们非常遥远的天文和精密信号链耦合在一起，甚至里面用到的一些实验方式都是让我心潮澎湃。
题目：Long term characterization of voltage references
作者：Hubert Halloin et al.
年份：2013
来源：arXiv:1312.5101
研究领域：超低频（mHz–Hz）电压基准噪声与温度特性
我先介绍一下背景：
激光干涉空间天线开路者号（LISA Pathfinder）是一个欧洲空间局的空间探测器，于2015年12月3日04:04 UTC时间发射升空。该任务将会对于演化激光干涉空间天线所需的科技进行检验。
演化激光干涉空间天线是欧洲空间局计划于2034年射入空间的一座引力波天文台；该任务预算费用为4亿欧元。
简单来说，就是一个引力波探测器，但是这个波的范围是：探测频率范围为0.1毫赫兹至10赫兹；这个东西可以追溯到宇宙诞生之初。
当然这个东西就是放大版的Virgo，放大了 60 倍，在深空中，更大，更灵敏。
人类第一次探测到了拂过地球的引力波。激光干涉引力波天文台（LIGO）的两台具有极高灵敏度的探测器（分别位于美国华盛顿州和路易斯安那州）捕获到了两个黑洞合并所引起的时空扭曲。5个月后，科学家宣布了这一消息。世界为之轰动，这一发现成为了2016年最重要的一则物理学新闻。自1916年爱因斯坦首次预言引力波以来，物理学家已经努力了近百年的时间，希望能找到引力波的直接证据。
引力波探测器的外形就好比一个双臂张开数千米的大写字母“L”，它们能探测到引力波引发的比质子直径还小的臂长变化。但如果只有一台超灵敏探测器，科学家就无法区分引力波导致的时空扭曲和环境中其他振动导致的臂长变化。此外，每台探测器还要负责一片不小的星空范围：其视野需要覆盖地球周围40%的星空——类似于你站在一片沙漠里，抬头、原地转圈所看到的星空范围，然后试图从这片星空中找到一颗暗弱的恒星。这就好比一个人身处广袤的沙漠，一边打转一边还要找到一颗晦暗不明的星星。
LIGO探测器需要两根探测臂还有另一个原因。引力波以光速传播，但除非引力波同时正对着两个探测器，否则两个探测器检测到的扭曲信号总会有毫秒级的时差。科学家就可以利用这个时间延迟计算出碰撞所在的方向，这便可以缩小寻找引力波源的范围。根据2015年的观测数据，该范围已经缩减到星空的2%，但对于目标搜寻而言这还是太大了。
里面的核心是一个干涉仪：
利用干涉仪，研究人员可以测量由引力波引起的自由质量的运动。10E 量级的变化-18 这里测量的是米——这对应的距离是 0.000 000 000 000 000 001 米的距离变化。因此，这些测量对外部干扰因素(如振动、温度)或内部因素(如光损失(吸收))极其敏感。因此，优化所有组件和光学器件是成功的重要因素。
开始建模
磨磨唧唧的，终于到这里了，先算这个引力波的情况，最后给出信号链的设计情况，那引力波咋了，照样也得咱们精密信号链的工程师算！
先通用的建模这个引力波的模型，然后再收集公开资料对卫星建模。（真是太鸡儿严谨了，呜呼～）
第一部分是按照 wiki 搞得，没有人比这个民间机构更加权威吧！
建模链条是：从 （应变）一路推到 （等效位移）→ （相位）→ 光电探测电流/电压→ ADC 端等效噪声/带宽指标。
三航天器等边三角形、臂长约 500 万 km、单程光程约 16 s、目标能在该臂长上分辨到 10 pm 量级长度变化，并给出探测频段与应变灵敏度量级。
从引力波 到干涉仪“等效位移”
把引力波看作 TT 规范下的平面波（对于写“信号链”，这个模型足够用，毕竟没有教授叼我）：
在低频近似（引力波波长 ）下，单臂的长度变化可以写成：
：臂长。LISA 量级为 （500 万 km）。
：天线方向图（入射方向、偏振、臂的几何投影）；其量级通常 。对“信号链指标估算”，可以先用 给出上界、再用 RMS 给出“平均响应”的量级。
应变直接变成位移，比例系数就是臂长
位移 到相位
激光干涉测的是相位。单程传播的相位变化：
如果等效为“往返”（类似迈克耳孙直觉），相位因子乘 2：
其中 是激光波长（如 1064 nm；文章里不写具体值，留作符号，别忘了我们是信号链，不是天文学家）。
相位到光电探测电流/电压（把“光学”接到“电学”）
在干涉仪工作点附近（锁定在正交点/暗条纹附近），光电探测输出可线性化：
把光电二极管 responsivity 记为 ，跨阻增益 ，则：
于是应变到电压的链路增益可以写成：
可以把每个乘子对应到一个子系统（光功率、读出、TIA、ADC）。
数字量级算： 到底意味着多少 pm、多少 rad？
应变 → 位移
取 ，：
这和材料里提到的“在 500 万 km 上探测 10 pm 量级变化”的目标量级是同一数量级。
反过来，如果系统的等效位移噪声是 ，对应的应变噪声约为
位移 → 相位
若取 （常见空间干涉用激光），往返等效相位摆幅：
也就是说： 级应变对应的是几十微弧度量级的相位摆动（极其小），这就是为什么后端信号链要把“相位噪声、频率噪声、时延噪声、读出噪声”压到极限。
频段与“为什么信号链长得像低频计量系统”
材料给出 LISA 探测频段大约在 Hz 到 Hz，且在 Hz 附近的可探测应变灵敏度量级给到 ****。
对于头秃的工程师面对的是 mHz 量级的超低频信号（周期：几分钟到几小时），另外单程光程约 16 s（500 万 km 的量级）使系统天然带有“时延/多普勒/相位跟踪”的复杂度。 （这个可是光哎，还能能 16S）
大一统的信号链指标
因为肯定这么复杂的东西，子单元很多，所以要把引力波探测写成一个“等效位移计”的指标链：（方便甩锅啊，可能机械在这里是最微不足道的东西了，电子系统那可太损了）
物理输入：（应变）
机械等效：
光学读出：
光电转换：
噪声指标统一到应变域：
二进宫
看到这里都是非常爱爱爱爱学习的宝宝，让我们继续建模。这次按 LPF 的实际架构：异频（heterodyne）干涉 + 数字相位计（phasemeter）来建模（这是 LPF OMS 的核心路线）。
LPF 关键背景参数
OMS（Optical Metrology System）用异频干涉测两块自由漂浮测试质量（TM）的相对位移与姿态，达到 pm / nrad量级精度，测量带宽约 1–30 mHz；LPF 的异频频率较低，文献中提到约 ~2 kHz量级（相对 LISA 正式任务 MHz 级更低）。
ESA 对 OMS/phasemeter 的公开描述：phasemeter 输出速率约 100 Hz，并用“Single Bin DFT”测异频信号相位（更像数字锁相/锁频的单频点解调）。
LPF 目标位移读出噪声：在 mHz 带宽内达到 10 pm/√Hz量级（常在 3–30 mHz 带内引用）。
物理量到相位：
LPF 的 OMS 读出本质是“长度变化等效相位变化”。对等效光程差 ：
这里用“单程等效”，如果文章采用“等效迈克耳孙往返”直觉，也可写 。两者差一个常数因子，所以在整篇文章一致即可。
把 LPF 的 10 pm/√Hz 位移噪声换成相位噪声取 （LPF/LISA 系列 Nd:YAG 1064 nm 是公开一致的描述）。
换成“cycle/√Hz”（phasemeter领域常用）：
也就是说10 pm/√Hz 对应的相位噪声只有几十微弧度/√Hz级别。
光电探测的异频信号模型（干涉强度 → 正弦载波）
异频干涉的 PD 输出可写成：
：异频频率（LPF ~ kHz 级）。
：我们想测的“相位”，其变化由 产生（上一节）。
：包含 shot noise、TIA 噪声、RIN 耦合项、拍频相关噪声等。
LPF/ LISA 不是“直接读干涉条纹强度”，而是把信息搬到一个窄带载波上，再用相位计提取 。
Phasemeter（数字相位计）的等效算法模型：单频点锁相/单 Bin DFT
ESA 的公开描述是：phasemeter 用 Single Bin DFT 在异频频点提取相位，输出速率约 100 Hz。 把它写成一个工程可读的“离散锁相/锁频解调”：对每个输出时刻 （更新周期 ），在一个长度为 的窗口内做单频点相关（等效数字锁相）：
这就是“single-bin DFT/lock-in”的标准形态；其中 是 ADC 采样周期； 是每次更新跨越的采样点数（决定 100 Hz 输出）。
这组式子把“相位估计噪声”清晰地和以下因素挂钩：
输入 SNR（载波幅度 与噪声谱）
窗口长度 （等效积分时间）
采样率/更新率（带宽与数据率），到这里的时候终于就到了信号链工程师熟悉的内容了。
把“电压噪声/ADC噪声”折算到“位移噪声”的统一公式
我想最终要在文章里给一个“可对比”的指标：等效位移噪声 或 **相位噪声 **。
链路最干净的写法是：
而相位估计误差的一个常用小信号近似（对高 SNR 正弦相位估计）：
：载波附近的等效输入电压噪声密度（包括 TIA + ADC + 1f-RIN 等）。
：相位计等效噪声带宽，约与积分时间成反比（窗口越长，相位噪声越低），这样俺就就能把“前端电子”与“最终 10 pm/√Hz 指标”直接连起来。
可视化
时域：真实位移 vs 相位计估计位移
这里故意设置了一个比较现实的情形：注入信号只有 5 pm peak @ 20 mHz，但相位计输入端（载波）有一定的加性噪声 + 低频 1/f 相位漂移
所以时域上会出现：“真实信号被噪声淹没”的效果，呈现一种时域看不见，但是不代表频域/统计工具看不见。
Welch 的平均把随机噪声“压平”，而线谱仍然突出，所以在 ASD 上能看到注入频点。
我在上面的文章这样写的：
异频干涉把位移信息搬到载波相位上，数字相位计再把相位恢复成位移；最终在 mHz 带内用 PSD/ASD 看信号与噪声地板。
总结
信号链的物理—数学起点（决定一切指标）
被测物理量
目标量：测试质量相对位移
关注频段：
所以这是一个极端超超超超超低频测量系统
位移 → 相位（第一关键放大）
激光干涉的本质优势在这里：
数值上：，不是在直接测 pm，而是在测 μrad 级相位；后面所有指标，本质都是“相位噪声预算”
相位计（Phasemeter）对信号链的硬性要求
相位计 ≈ 在载波 处做 单频点锁相 / 单 bin DFT：
本质是一个极窄带数字锁相放大器
相位噪声 → 位移噪声指标反推
目标位移噪声（LPF）：
对应相位噪声要求：
这是整个电子信号链必须满足的“终极指标”
对 ADC或者TIA 的直接要求
ADC 不是“位数问题”，而是：
相位抖动（aperture jitter），近载波相位噪声，低频 flicker 调制
在异频下，ADC 抖动 → 相位噪声：
所以：kHz 载波 + mHz 目标，这个时基与时钟质量比 ENOB 更关键
后记
LISA Pathfinder 的信号链不是在“提高瞬时分辨率”，而是在一个横跨 6–7 个数量级频率范围的系统中，保证相位信息在极低频处仍然是可统计、可建模、可验证的。
而我们只是研究其中一个微不足道的 REF 而已，嘻嘻嘻（不嘻嘻）
看回主角
这个信号链的难度我们也看到了，复杂无比，那我们开头的论文呢，其实是聚焦在里面的 REF，就是文章里面经常出现的东西：
研究背景
论文服务的直接背景是 eLISA（空间引力波干涉仪） 项目，其关键测量频段为：0.1 mHz – 10 Hz
在这个频段内电压基准被用于：ADC 参考，光电二极管偏置，电容/惯性质量读出，系统对电压稳定性的要求是频谱意义上的稳定，而不是传统的“短期噪声”或“多年漂移”。
作者明确指出一个现实问题：
器件厂商几乎不在 mHz 频段给出噪声或稳定性指标。
因此现有 datasheet 无法直接判断器件是否适合 eLISA必须做专门的超低频实验表征。
与以往工作的区别
论文强调与既有研究的差异：
以往研究：0.1 Hz（短期）或 数月/数年（长期标定、计量）
本文：10⁻⁵ Hz – 10 Hz，覆盖“分钟—小时—天”尺度
作者将其定义为“低频连续谱稳定性”研究，而不是传统的“漂移”或“噪声”单一概念。
被测试器件（Section II + Table I/II）
被测电压基准一览（共 8 种）
论文测试了 8 款电压基准，覆盖三类技术：
Datasheet 参数对比（Table II）
作者整理了各器件 datasheet 中的关键指标（统一换算到 ppm）：温度漂移（ppm/K），线性调整率（ppm/V），0.1–10 Hz 峰峰值噪声，长期漂移（ppm/kHr）。
但作者明确说明：
这些参数仅用于“事前筛选”，并不代表 mHz 频段下的真实性能。
实验方法（Section II：Experimental Setup）
核心测量思想：差分法
论文没有使用“更高级的基准”作为参考，而是采用：两颗同型号电压基准，测量它们的差分输出
假设：
这是在无法获得更优参考源情况下的标准噪声测量方法。
为什么不使用 trim？
作者明确说明外部 trim 网络会：增加温漂，其实是会增加低频噪声，因此：所有 DUT 均未使用 trim，即使输出存在 mV 级偏差也接受，这是一个重要的实验设计选择，而非疏忽。
放大与采集链路（别急我会大写特写）
前端放大器：ADA4528，零漂移，超低噪声，使用了非反相放大，增益 401。
数据采集：Agilent 34411A，采样率 50 Hz（整周期工频平均）作者对放大器与采集链路噪声进行了单独标定，确保其噪声低于 DUT 噪声至少 2 个数量级。
温度控制策略
为消除非本征噪声：
DUT 与放大器：浸没在 Galden HT110 恒温液体中
温控：恒温 23°C或叠加 0.2 mHz 正弦温度调制
温度测量：使用双 Pt100mK/√Hz 级噪声;这样可以：抑制空气流动，抑制湿度变化，抑制热电势不均匀，反正就是把所有的不可靠元素都控制。
实验噪声源分析
作者系统分解了测量系统的噪声来源，并给出 Fig.2：
噪声源分类
采集链路噪声，电源线性调整噪声，温度引入噪声（恒温 vs 空气）
关键结论
采集链路噪声：10⁻³ – 2×10⁻² ppm/√Hz，远低于目标
电源影响：可忽略（这里是怎么说呢，大概意思就是用了很低噪音的电源）
温度噪声：在空气中 → 超过目标，在恒温液体中 → 低于目标
作者结论：
在该实验中，主导噪声源是温度波动。
温度系数测量
为什么不采用 datasheet 的 TC？
作者指出Datasheet 使用的是 box method：全温区最大变化 / 温区宽度；而实际工作点：电压–温度曲线通常是 S 形，局部斜率（dV/dT）与 box TC 可能差别很大
实际测量方法
对 DUT 施加：18–28°C，使用周期 5000 s 的正弦温度调制，对电压输出做：正弦 + 二次多项式拟合
得到：
单颗 TC
差分 TC
结果解读（以 LT1021 为例）
单颗 TC：≈ −0.18 ppm/K
差分 TC：≈ −0.67 ppm/K
Datasheet 最大 TC：5 ppm/K
作者强调：实测结果与 datasheet 数量级一致，差异主要来自测量方法不同以及样本数量有限（每种仅 2 颗）。
长期稳定性与低频噪声
短时间筛选（20 ks，Fig.3）
XFET（ADR435/445）：噪声显著偏高，低频斜率陡
其他器件：在 ~1 mHz 附近接近目标
长时间测量（55 h，Fig.4）
对筛选后的器件进行更长时间测量，得到关键结论：
MAX6126：全频段最低噪声，温度影响几乎可忽略
LTC6655（Bandgap）：噪声接近最佳 buried zener，温度影响小
LT1021：与 AD587UQ 相当，在 0.5 mHz 处超过目标，温度贡献在 0.1 mHz 附近显著
MAX6350：低频呈现较陡 1/f 行为
论文结论（Conclusion）
作者的结论非常克制、明确：
没有任何测试器件完全满足 eLISA 要求，这个MAX6126 最接近目标，LT1021 与 AD587 仍是可考虑的 buried zener 方案，LTC6655 是一种有前景的 bandgap 替代路线，XFET 在 <0.1 Hz 频段噪声偏大
作者也明确声明：
本文不讨论辐射效应、空间可靠性，仅限本征噪声与温度特性。
真没了
一篇方法论扎实、结论保守的实验论文，系统性地证明了：在 10⁻⁵–10 Hz 频段，电压基准的真实性能必须通过专门实验获得，而不能从 datasheet 推断。
可写的东西还有很多，先看看这个结果咋样。（来看看真真正正的信号链，天天抱着个三极管水文真没有意思。）
https:///abs/1312.5101
https://www./watch?v=cS7s_nYnA0E
https://www./media/en/technical-documentation/data-sheets/adr430_431_433_434_435.pdf
https:///products/galden-ht-110?srsltid=AfmBOopcp3hoUlOmJnIiN3BOfrOa2TZ99RYqIJhP4kEn7HVNCZtvytcf
https://www./past-job-offers/?_paged=2