3. 量子相位估计的标准线路与工作流程

量子相位估计（Quantum Phase Estimation, QPE）是量子算法中最重要的基础子程序之一，Shor 算法、量子模拟、HHL 线性方程求解等都依赖它。本节聚焦 QPE 的标准线路结构，把整条线路拆成四个关键环节：相位 kickback 机制、受控酉算子应用、逆量子傅里叶变换、测量与相位读出。理解这条线路，关键不在于记住门序列，而在于理解"相位信息如何从本征态转移到计数寄存器"这一核心物理过程。

一、问题设定与寄存器划分

QPE 要解决的问题是：给定一个酉算子 U 及其本征态 |ψ⟩，满足 U|ψ⟩ = e^{2πiφ}|ψ⟩，估计相位 φ ∈ [0,1)。线路使用两个寄存器：

• 计数寄存器：t 个量子比特，初始化为 |0⟩^{⊗t}，最终用于存放 φ 的二进制近似。
• 目标寄存器：存放本征态 |ψ⟩，位数取决于 U 的作用空间，全程保持在本征态上（理想情况下不被改变）。

估计精度与 t 直接相关：若 φ 可精确表示为 t 位二进制小数，则测量以概率 1 得到精确结果；否则以高概率得到最接近的 t 位近似，误差量级约为 2^{-t}。这是选择 t 的基本依据。

二、第一步：Hadamard 层与叠加制备

对计数寄存器的每个比特施加 Hadamard 门，得到均匀叠加态：

|0⟩^{⊗t}  --H^{⊗t}-->  (1/√(2^t)) Σ_{k=0}^{2^t-1} |k⟩

这一步本身没有"计算"任何东西，它的作用是把计数寄存器变成一把"相位标尺"——每个基态 |k⟩ 后续会以不同次数调用 U，从而累积不同倍数的相位。可以把它理解为：叠加态让线路一次性并行地"询问"了所有 2^t 个相位倍数。

三、第二步：受控酉算子与相位 kickback

这是整条线路的物理核心。对计数寄存器的第 j 个比特（从 0 开始计，对应权重 2^j），以它为控制位、目标寄存器为受控目标，施加受控 U^{2^j}：

控制位为 |1⟩ 时：|1⟩|ψ⟩ --> |1⟩(U^{2^j}|ψ⟩) = e^{2πi·2^j·φ}|1⟩|ψ⟩
控制位为 |0⟩ 时：|0⟩|ψ⟩ --> |0⟩|ψ⟩（不变）

由于 |ψ⟩ 是 U 的本征态，U^{2^j} 作用后只产生一个全局相位因子 e^{2πi·2^j·φ}。这个相位因子"回踢"（kickback）到控制位上，等价于对控制位施加了一个绕 Z 轴的相位旋转：

|0⟩ --> |0⟩
|1⟩ --> e^{2πi·2^j·φ}|1⟩

目标寄存器始终停留在 |ψ⟩，相位信息被完整地"搬"到了计数寄存器上。全部 t 个受控操作完成后，计数寄存器的联合态为：

(1/√(2^t)) Σ_{k=0}^{2^t-1} e^{2πi·k·φ} |k⟩ ⊗ |ψ⟩

注意相位因子 e^{2πi·k·φ} 恰好是 k 的线性函数——这正是逆量子傅里叶变换能够"解码"它的原因。工程实现上，U^{2^j} 通常不是直接构造大矩阵，而是通过重复施加 U 共 2^j 次实现；对稀疏或结构化酉算子，也可用快速幂或专用线路优化。

四、第三步：逆量子傅里叶变换

此时计数寄存器的态可以写成 QFT|φ̃⟩ 的形式，其中 |φ̃⟩ 是 φ 的 t 位二进制近似态。因此施加逆量子傅里叶变换 QFT† 即可把它"还原"为计算基态：

QFT† · (1/√(2^t)) Σ_k e^{2πi·k·φ} |k⟩  ≈  |φ̃⟩ = |φ_1 φ_2 ... φ_t⟩

若 φ = 0.φ_1φ_2...φ_t（二进制小数），则 QFT† 后计数寄存器精确处于 |φ_1φ_2...φ_t⟩。若 φ 不能精确表示，则得到的是集中在最近二进制小数附近的概率分布。逆 QFT 的标准实现由 Hadamard 门与受控相位门级联构成，门数约为 O(t^2)，是整条线路中门数增长最快的部分，也是精度与资源权衡的关键点。

五、第四步：测量与相位读出

对计数寄存器做计算基测量，得到 t 位二进制串 m = m_1m_2...m_t，读出相位估计值：

φ_est = 0.m_1m_2...m_t = Σ_{i=1}^{t} m_i · 2^{-i}

例如 t = 4，测得 m = 1011，则 φ_est = 1/2 + 1/8 + 1/16 = 0.6875。若真实 φ = 0.7，误差为 0.0125，在 2^{-4} = 0.0625 的量化步长之内，属于合理近似。实际运行时单次测量结果具有随机性，通常需要多次重复实验统计频率分布，取概率峰值对应的比特串作为估计值，并观察分布是否集中在单一峰值附近来判断线路是否工作正常。

六、标准线路全景

下图把上述四个环节按信号流方向串起来，并标出每个阶段计数寄存器与目标寄存器的状态变化，便于对照线路逐段排查。

把线路按信号流读一遍：Hadamard 层把计数寄存器铺成"相位标尺"，受控 U^{2^j} 把目标寄存器的本征相位按不同倍数回踢到各控制位，逆 QFT 把线性相位编码解码为计算基比特串，测量则把比特串翻译成十进制相位估计。四个环节中，受控酉算子的幂次对应关系和逆 QFT 的相位门符号是最容易出错的两处，调试时应优先核对。