3. 量子测量噪声的表征与建模

在量子计算实验中，我们最终得到的不是量子态的连续振幅，而是一串离散的经典比特。这个从量子态到经典比特的转换过程，本身就带有噪声。理解这些噪声从哪里来、如何用数学描述、怎样从真实设备数据中估计出来，是数据质量控制的第一道关口。本节围绕三类主要噪声——读出错误、弛豫与退相干——展开，并给出可操作的建模与拟合方法。

三类噪声的物理来源与区分

读出错误（readout error）发生在测量环节本身。理想情况下，量子态 |0⟩ 应读出 0，|1⟩ 应读出 1；但实际设备中，由于读出谐振腔信号区分度不足、放大器噪声或阈值判定偏差，会出现 |0⟩ 被误判为 1、|1⟩ 被误判为 0 的情况。它描述的是“测量仪器看错了”，与量子态在测量前的演化无关。

弛豫（relaxation，T₁过程）是能量耗散：激发态 |1⟩ 会以一定概率衰减到基态 |0⟩。其特征时间是 T₁，表示布居数从 1 衰减到 1/e 所需的时间。如果测量脉冲到来之前等待时间过长，|1⟩ 可能已经衰减，导致读出结果偏向 0。

退相干（dephasing，T₂过程）是相位信息的丢失。它不改变 |0⟩ 和 |1⟩ 的布居数，但会破坏叠加态中的相位关系。在 Ramsey 干涉或回波实验中，T₂ 表现为干涉条纹幅度的衰减。退相干本身不直接翻转比特，但会与弛豫、门操作误差耦合，间接影响测量结果的统计分布。

三者的关键区别在于：读出错误是“测量映射”的错误，弛豫是“测量前态已变”，退相干是“相位信息丢失”。在建模时，读出错误通常用混淆矩阵描述，而弛豫和退相干则通过时间依赖的量子通道作用于态，再经过读出映射。

用概率模型描述测量过程

最直接的建模方式是把测量看作一个概率映射。设真实制备的态为 |0⟩ 或 |1⟩，测量结果为 0 或 1。定义读出混淆矩阵 M：

M = [[P(0|0), P(0|1)],
     [P(1|0), P(1|1)]]

其中 P(r|s) 表示真实态为 s 时读出结果为 r 的条件概率。理想情况下 M 为单位矩阵。实际设备中，P(1|0) 称为“0 态误判为 1 的概率”，P(0|1) 称为“1 态误判为 0 的概率”。这两个参数就是读出错误率。

如果只考虑读出错误，测量结果的概率分布为 p_read = M · p_true，其中 p_true 是真实态布居数向量。但真实实验中，制备的态本身可能已经受到弛豫影响。例如，制备 |1⟩ 后等待时间 t，则真实态布居数变为：

p_true(1) = exp(-t / T1)
p_true(0) = 1 - exp(-t / T1)

再经过读出矩阵，最终观测到的 P(1) 为：

P_obs(1) = P(1|1) * exp(-t/T1) + P(1|0) * (1 - exp(-t/T1))

这个表达式把弛豫和读出错误统一在一个可拟合的模型中。退相干的影响通常不直接改变布居数，但在涉及叠加态的测量中，它会降低干涉可见度，从而影响从干涉条纹提取的 T₂ 估计。

从真实设备数据拟合读出错误率

实际做法是：制备一组已知态，统计测量结果，反推混淆矩阵。以超导量子比特为例，典型流程如下。

第一步，制备 |0⟩ 态并测量 N 次，记录结果为 1 的次数 n₁₀。则 P(1|0) ≈ n₁₀ / N。第二步，制备 |1⟩ 态并测量 N 次，记录结果为 0 的次数 n₀₁。则 P(0|1) ≈ n₀₁ / N。第三步，用这两个参数构造混淆矩阵，并检查其是否满足行和为 1。

假设某设备上，制备 |0⟩ 测量 10000 次，有 320 次读出 1；制备 |1⟩ 测量 10000 次，有 780 次读出 0。则：

P(1|0) = 320 / 10000 = 0.032
P(0|1) = 780 / 10000 = 0.078
M = [[0.968, 0.078],
     [0.032, 0.922]]

注意，这里 P(0|0)=1-0.032=0.968，P(1|1)=1-0.078=0.922。该矩阵不是对称的，说明 1 态更容易被误判为 0，这通常与弛豫有关——|1⟩ 在测量前部分衰减到了 |0⟩。

为了分离读出错误和弛豫，可以改变制备 |1⟩ 后的等待时间 t，重复上述测量，得到 P_obs(1) 随 t 的变化曲线。用模型 P_obs(1) = P(1|1)·exp(-t/T₁) + P(1|0)·(1-exp(-t/T₁)) 拟合，即可同时估计 T₁ 和读出错误率。如果拟合残差较大，说明还存在其他噪声源，如退相干或门操作误差。

噪声矩阵在数据质控中的使用

得到混淆矩阵后，可以用于校正测量结果。如果观测到布居数 p_obs，真实布居数 p_true 可通过 p_true = M⁻¹ · p_obs 估计。但要注意，当 M 接近奇异或噪声较大时，直接求逆会放大误差。实践中常采用最大似然估计或约束最小二乘，保证 p_true 的每个分量在 [0,1] 内且和为 1。

此外，混淆矩阵本身也是设备健康度的指标。如果 P(1|0) 或 P(0|1) 突然增大，可能意味着读出谐振腔频率漂移、放大器工作点变化或磁通偏置异常。把读出错误率随时间的变化画成趋势图，可以作为异常诊断的早期信号。

把上述流程固化为可重复的质控步骤，可以形成一份简明的读出错误率标定记录。每次标定至少记录：制备态、测量次数、等待时间、观测到的误判次数、计算出的 P(1|0) 和 P(0|1)、拟合得到的 T₁ 以及拟合残差。当残差超过预设阈值（例如 0.01）时，应检查是否存在退相干或门误差的额外贡献，而不是简单接受读出错误率结果。

最后需要强调，读出错误率不是固定常数。它随设备温度、磁通偏置、读出脉冲功率和测量时长变化。在长时间实验或批量任务中，建议每隔一段时间重新标定一次，并把标定结果与实验数据的时间戳对齐，这样才能在后续异常诊断中区分“设备漂移”和“真实量子态变化”。