4. 量子测量系统建模
量子测量系统的漂移补偿，第一步不是调参，而是建立一条从量子层面到经典控制层面的“可计算通路”。量子比特、探针或传感单元的真实演化由密度矩阵与量子主方程描述，但补偿器、滤波器和执行器只认识经典电压、频率或相位。因此，建模的核心任务是把量子主方程中的退相干、弛豫与测量反作用，逐步映射为经典漂移的传递函数、状态空间模型和随机微分方程，使漂移既能被解释，也能被预测和抑制。
一、量子主方程：漂移的物理源头
开放量子系统的标准描述是 Lindblad 主方程。设系统密度矩阵为 ρ，哈密顿量为 H，则
dρ/dt = -i/ħ [H, ρ] + Σ_k ( L_k ρ L_k† - 1/2 {L_k† L_k, ρ} )其中第一项是相干演化，第二项是耗散与退相干。对量子测量而言，最关键的 Lindblad 算符通常包括弛豫算符 σ₋ 和纯退相干算符 σ_z。若把测量过程视为对系统施加一个弱耦合探针，则测量强度 Γ_m 会同时带来信息增益和反作用噪声。漂移的物理源头就藏在这些算符的系数里：T₁ 弛豫导致幅值缓慢衰减，T₂ 退相干导致相位随机游走，而测量反作用则表现为额外的白噪声或 1/f 噪声。
以超导量子比特为例，假设工作频率 f₀ = 5.0 GHz，T₁ = 40 μs，T₂ = 25 μs，测量积分时间 τ_m = 300 ns。由 T₂ 与 T₁ 可估算纯退相干时间 T_φ：1/T₂ = 1/(2T₁) + 1/T_φ，得到 T_φ ≈ 33.3 μs。这些时间常数不是抽象数字，它们直接决定漂移的截止频率：T₁ 对应幅值漂移的低频分量，T₂ 对应相位漂移的带宽。若测量重复率 f_rep = 1/τ_m ≈ 3.33 MHz，则采样率远高于漂移带宽，说明经典补偿器有足够余量跟踪慢漂移。
二、从主方程到经典漂移传递函数
把量子主方程线性化，是获得传递函数的关键。设系统状态在稳态附近做小偏差展开，令 δρ 为密度矩阵偏差，δu 为控制输入偏差，δy 为测量输出偏差。对主方程在稳态处求一阶变分，可写成线性时不变形式：
d(δρ)/dt = A_q δρ + B_q δu
δy = C_q δρ + D_q δu其中 A_q 是量子 Liouvillian 矩阵，B_q、C_q、D_q 由测量算符和探针耦合决定。对上述方程做拉普拉斯变换，消去 δρ，得到从控制输入到测量输出的传递函数：
G_q(s) = C_q (sI - A_q)^(-1) B_q + D_q这个 G_q(s) 就是量子系统在经典频域中的“等效对象”。它的极点由 T₁、T₂ 和测量带宽共同决定，零点则由探针耦合方式决定。实际工程中，常把 G_q(s) 近似为低通环节与延迟环节的串联：
G_q(s) ≈ K_q / (1 + s/ω_c) · e^(-sτ_d)其中 ω_c 是漂移截止角频率，τ_d 是测量与反馈延迟。以 T₂ = 25 μs 为例，若取 ω_c ≈ 1/T₂ = 4×10⁴ rad/s，则截止频率约 6.4 kHz。若反馈延迟 τ_d = 200 ns，则相位滞后在 6.4 kHz 处约为 0.8°，在 100 kHz 处约为 7.2°。这些数值直接决定补偿器带宽上限：带宽越高，延迟引起的相位裕度损失越大。
三、状态空间模型：把漂移写成可观测状态
传递函数适合频域分析，但补偿器实现通常需要时域状态空间模型。把漂移源分为三类状态：幅值漂移 x₁、相位漂移 x₂、测量偏置 x₃。它们分别由 T₁ 弛豫、T₂ 退相干和探针零点漂移驱动。连续时间状态空间模型可写为：
dx/dt = A x + B w + B_u u
y = C x + v其中 x = [x₁, x₂, x₃]ᵀ，w 是过程噪声，v 是测量噪声，u 是补偿输入。A 矩阵的对角元由各漂移的衰减率决定，例如 A₁₁ = -1/T₁，A₂₂ = -1/T₂，A₃₃ = -1/τ_bias。B 矩阵描述噪声如何进入各状态，C 矩阵描述测量如何观测这些状态。若测量同时给出幅值和相位，则 C = [[1,0,0],[0,1,0]]；若只测幅值，则相位漂移只能通过耦合项间接观测。
离散化后，状态空间模型变为：
x[k+1] = A_d x[k] + B_d w[k] + B_ud u[k]
y[k] = C_d x[k] + v[k]以采样周期 T_s = 1 μs 为例，若 T₁ = 40 μs，则 A_d,11 = exp(-T_s/T₁) ≈ 0.9753；若 T₂ = 25 μs，则 A_d,22 ≈ 0.9608。这意味着幅值漂移每微秒衰减约 2.5%，相位漂移每微秒衰减约 3.9%。衰减越慢，漂移越接近随机游走，补偿器越需要积分作用；衰减越快，漂移越接近白噪声，补偿器越应侧重滤波而非积分。
四、随机微分方程：描述漂移的随机本质
量子测量漂移不是确定性偏差，而是随机过程。用随机微分方程描述更贴近物理本质。对幅值漂移，常用 Ornstein-Uhlenbeck 过程：
dx₁ = -(1/T₁) x₁ dt + σ₁ dW₁对相位漂移，若存在 1/f 噪声，则需用分数阶布朗运动或叠加多个 OU 过程近似：
dx₂ = -(1/T₂) x₂ dt + σ₂ dW₂ + σ_f dW_f其中 dW₁、dW₂ 是标准维纳过程增量，dW_f 是低频闪烁噪声增量。σ₁、σ₂、σ_f 分别由 T₁ 噪声、T₂ 噪声和 1/f 噪声幅度标定。若已知 Allan 方差在 τ = 1 s 处为 1×10⁻⁵，则可反推 σ_f 的近似量级。随机微分方程的价值在于：它把“漂移有多快、有多随机”变成可仿真的参数，从而让卡尔曼滤波或模型预测控制有明确的噪声协方差 Q 和 R 可用。
五、三类模型的选用与衔接
传递函数、状态空间和随机微分方程不是三选一，而是同一物理对象在不同工程环节的三种表达。频域设计补偿器时，用 G_q(s) 画 Bode 图，检查截止频率和相位裕度；实现数字滤波器或卡尔曼滤波时，用离散状态空间模型递推；做蒙特卡洛仿真或噪声预算时，用随机微分方程生成漂移轨迹。三者之间的参数必须一致：G_q(s) 的极点应等于 A 矩阵的特征值，SDE 的衰减率应等于 A 的对角元，噪声强度 σ 应通过 Q 矩阵进入状态空间。
一个常见的衔接错误是：频域设计时把截止频率设得过高，忽略了 τ_d 带来的相位滞后；时域实现时又用白噪声近似 1/f 噪声，导致低频漂移被低估。正确的做法是先用 Allan 方差或功率谱密度确认漂移的噪声类型，再决定用 OU 过程还是多阶 OU 叠加。若 Allan 方差在 0.1 s 到 10 s 区间斜率为 -1/2，说明以白噪声调频为主，可用单 OU 过程；若斜率接近 0，说明以 1/f 噪声为主，需增加闪烁噪声状态。
六、建模结果如何验证
模型建好后，不能只看公式是否漂亮，而要用可观测指标验证。第一，比较模型预测的漂移功率谱与实测功率谱，在 1 Hz 到 10 kHz 范围内误差应小于 3 dB。第二，用模型生成 100 条漂移轨迹，统计其 Allan 方差，与实测 Allan 方差在 τ = 1 ms 到 1 s 范围内对比，偏差应小于 20%。第三，把模型接入补偿器闭环，检查开环相位裕度是否大于 45°，闭环阶跃响应的超调是否小于 15%。若这三项都满足，说明从量子主方程到经典漂移模型的映射是可信的，后续的漂移补偿与稳定性控制才有可靠基础。