3. 量子测量与估计理论

量子传感的精度上限，本质上由“测量怎么读、参数怎么估”共同决定。前两章我们建立了量子态与演化的语言，这一章要把语言接到工程指标上：给定一个待估参数（如磁场 \(B\)、相位 \(\phi\)、频率 \(\omega\)），探测器最终给出一个随机结果，我们如何从统计上判断“这套方案最好能准到什么程度”。核心工具链是：投影测量与 POVM 描述“能读到什么”，Fisher 信息刻画“读数对参数有多敏感”，Cramér-Rao 界给出“估计方差的下限”，量子 Fisher 信息与量子 Cramér-Rao 界则回答“在量子力学允许范围内，极限在哪里”。

投影测量：最基础的一类量子读数

投影测量由一组正交投影算符 \(\{P_m\}\) 描述，满足 \(P_m^\dagger=P_m\)、\(P_mP_{m'}=\delta_{mm'}P_m\)、\(\sum_m P_m=I\)。若系统态为 \(\rho\)，得到结果 \(m\) 的概率为 \(p(m)=\mathrm{Tr}(P_m\rho)\)，测量后态坍缩为 \(P_m\rho P_m/p(m)\)。在量子传感中，最典型的投影测量是：先施加参数化演化 \(U(\phi)\)，再用一组正交基（如 \(|0\rangle,|1\rangle\)）做投影，统计两端口计数。它的优点是物理实现直接、噪声模型清晰；缺点是只能读取一组互斥结果，若参数信息分散在非正交的多个通道中，投影测量会丢失部分可提取信息。

POVM：把“测量”推广到非正交结果

POVM（正算符值测量）用一组厄米正定算符 \(\{E_m\}\) 描述，满足 \(E_m\ge 0\)、\(\sum_m E_m=I\)，但不再要求 \(E_m\) 是投影。概率仍为 \(p(m)=\mathrm{Tr}(E_m\rho)\)。POVM 的价值在于：它允许测量结果对应的“效应”彼此不正交，从而在一次测量中同时获取多个非对易可观测量的部分信息。例如在相位估计中，若只用投影测量，相位信息可能只出现在一个 quadrature 上；采用合适的 POVM（如外差或自适应测量）可以更均衡地提取信息。工程上，POVM 常作为“理想测量”的上界模型，用来判断当前投影方案离最优测量还差多少。

Fisher 信息：敏感度的量化

设测量结果为 \(x\)，其概率分布 \(p(x|\theta)\) 依赖待估参数 \(\theta\)。经典 Fisher 信息定义为

\[ F(\theta)=\sum_x p(x|\theta)\left(\frac{\partial \ln p(x|\theta)}{\partial \theta}\right)^2 \]

连续变量时求和改为积分。它衡量的是：分布随 \(\theta\) 变化时，对数似然函数的斜率平方的期望。\(F(\theta)\) 越大，说明结果分布对参数越敏感，单次测量能携带的信息越多。对 \(N\) 次独立重复测量，总 Fisher 信息为 \(N F(\theta)\)。在量子情形中，若固定测量策略，则 \(F(\theta)\) 是经典 Fisher 信息；若允许在所有合法 POVM 中取最大，则得到量子 Fisher 信息 \(F_Q(\theta)\)，满足 \(F(\theta)\le F_Q(\theta)\)。

Cramér-Rao 界：估计方差的下限

对任意无偏估计量 \(\hat{\theta}\)，其方差满足

\[ \mathrm{Var}(\hat{\theta})\ge \frac{1}{F(\theta)} \]

这就是 Cramér-Rao 界（CRB）。若考虑 \(N\) 次独立测量，则 \(\mathrm{Var}(\hat{\theta})\ge 1/(N F(\theta))\)。它告诉我们：无论后处理算法多聪明，只要测量方案固定，估计精度不可能突破这个下限。达到 CRB 的估计量称为有效估计量，最大似然估计在样本量较大时通常渐近达到该界。量子版本为量子 Cramér-Rao 界：\(\mathrm{Var}(\hat{\theta})\ge 1/(N F_Q(\theta))\)，其中 \(F_Q(\theta)\) 由对称对数导数决定。对纯态 \(|\psi(\theta)\rangle\)，有

\[ F_Q(\theta)=4\left(\langle\partial_\theta\psi|\partial_\theta\psi\rangle-|\langle\psi|\partial_\theta\psi\rangle|^2\right) \]

这一表达式是量子传感中判断“是否接近海森堡极限”的常用入口。

从标准量子极限到海森堡极限

对 \(N\) 个独立粒子（或 \(N\) 次独立测量）的经典策略，Fisher 信息通常正比于 \(N\)，于是 CRB 给出 \(\Delta\theta\propto 1/\sqrt{N}\)，称为标准量子极限（SQL）。若使用纠缠态（如 GHZ 态）或压缩态，量子 Fisher 信息可正比于 \(N^2\)，CRB 变为 \(\Delta\theta\propto 1/N\)，即海森堡极限。需要强调：海森堡极限是量子力学允许的渐近上界，实际系统受退相干、损耗、探测效率限制，往往只能部分逼近。工程上判断“是否接近量子极限”，常用方法是：先由理论或标定得到 \(F_Q\)，再实测估计方差，比较 \(\mathrm{Var}(\hat{\theta})\) 与 \(1/(N F_Q)\) 的差距；若差距主要来自测量策略而非探针态，则说明还有 POVM 优化空间。

一个可复现的相位估计示例

考虑 Mach-Zehnder 干涉仪，两臂相位差 \(\phi\) 为待估参数。输入 \(|1,0\rangle\)（单光子）时，输出端口概率为 \(p_0=(1+\cos\phi)/2\)、\(p_1=(1-\cos\phi)/2\)。计算经典 Fisher 信息：

p0 = (1 + cos φ)/2,  p1 = (1 - cos φ)/2
∂ln p0/∂φ = -sin φ / (1 + cos φ)
∂ln p1/∂φ =  sin φ / (1 - cos φ)
F(φ) = p0·(∂ln p0)² + p1·(∂ln p1)²
     = sin²φ / (1 - cos²φ) = 1

即单光子投影测量下 \(F=1\)，CRB 为 \(\mathrm{Var}(\hat{\phi})\ge 1/N\)。若改用压缩真空注入，理论 \(F_Q\) 可随压缩因子指数提升；若用 \(N\) 光子 NOON 态，理想无损耗下 \(F_Q=N^2\)，CRB 为 \(1/N^2\)。实际实验中，探测器暗计数、干涉可见度 \(V<1\) 会引入额外噪声，使有效 Fisher 信息下降。一个实用的标定流程是：固定 \(\phi\) 扫描若干点，记录两端口计数 \(n_0,n_1\)，用最大似然估计反推 \(\hat{\phi}\)，再对多次重复求方差，与 \(1/(N F_{\text{eff}})\) 对比，其中 \(F_{\text{eff}}\) 可由可见度修正 \(F_{\text{eff}}=V^2 F_{\text{ideal}}\) 近似估计。若实测方差明显高于该界，优先检查测量基是否最优、是否存在经典串扰或相位漂移。

常见误区与检查清单

• 把 CRB 当成“一定能达到”：CRB 是无偏估计的方差下限，只有有效估计量或大样本 MLE 才渐近达到；小样本、有偏估计或非正则分布下需谨慎。
• 忽略测量策略对 F 的影响：同一量子态，不同 POVM 给出的经典 Fisher 信息不同；优化测量基往往比换探针态更易落地。
• 把 \(F_Q\) 直接当实验指标：\(F_Q\) 是理想上界，实际需扣除损耗、暗计数、有限探测效率，得到有效 Fisher 信息再比较。
• 混淆 \(N\) 的含义：\(N\) 是独立重复次数还是粒子数，直接影响 SQL 与海森堡极限的标度判断。