📚 课程目录第4章振幅估计的标准量子电路

4. 振幅估计的标准量子电路

振幅估计（Amplitude Estimation, AE）的核心任务，是估计一个量子态在某个好子空间上的振幅 a，从而得到概率 a²。它之所以能相对经典蒙特卡洛获得二次加速，关键不在于某个单独的门，而在于把「算符 A 的制备」「受控 A 与受控 A†」「反射算符」以及「量子相位估计（QPE）」组合成一个结构严谨的标准电路。本节围绕这条电路主线，讲清每个部件的构造方式、它们之间的接口关系，以及组合时必须满足的条件。

一、算符 A 的构造：把问题编码进振幅

振幅估计的起点是定义一个算符 A，使得

A|0⟩_n = √(1−a)|ψ₀⟩|0⟩ + √a|ψ₁⟩|1⟩

其中最后一个辅助量子比特是「目标比特」，它的 |1⟩ 分量振幅恰为 √a。也就是说，A 必须把待估计的概率「装」进一个辅助比特的振幅里。构造 A 的常见做法是：先用一组态制备门把数据寄存器制备到叠加态，再用一个多控旋转门（如 R_y(2θ) 或多控 R_y）把目标比特旋转到所需角度，使 sin²θ = a。例如在金融衍生品定价中，A 负责把资产价格路径的分布制备到数据寄存器，并把「期权收益是否为正」编码到目标比特；在 Grover 型搜索中，A 就是均匀叠加制备加一个标记预言机。

构造 A 时要特别注意两点：其一，A 必须是可逆的酉算符，任何非酉的经典采样步骤都要改写成酉形式（例如用比较器 + 受控旋转代替直接赋值）；其二，A 的代价直接决定整个算法的门数，因此实际实现中应尽量复用已有制备电路，避免为每个振幅单独构造。

二、受控 A 与受控 A†：QPE 的调用接口

QPE 需要把 A 的幂次以受控方式作用到工作寄存器上。标准做法是：在 QPE 的计数寄存器（t 个比特）与工作寄存器之间，施加受控-A^{2^j}（j = 0, …, t−1），即第 j 个计数比特控制 A 作用 2^j 次。由于 A 是酉算符，A^{2^j} 可通过重复平方得到，但更常见的实现是直接构造受控版本，或利用 A 的结构做受控分解。

受控 A† 则出现在 QPE 之后的「逆运算」环节：在完成相位估计、提取出相位信息后，需要把工作寄存器恢复到初始状态，以便后续的反射与迭代不破坏已积累的相位。若 A 由若干基本门复合而成，受控 A† 就是各基本门受控版本的逆序共轭。实现时要注意控制比特的「开」条件与目标比特的相位约定一致，否则会在多次迭代中累积系统性偏差。

工程提示：受控 A 与受控 A† 的控制线应尽量复用同一组计数比特，并在电路图上明确标注每个受控门的幂次，避免把 A^{2^j} 误写成 A^j。这是初学者最常见的接口错误。

三、反射算符：把相位差放大成可测信号

仅有 QPE 只能得到 A 的某个本征相位，要把它转化为对 a 的估计，还需要反射算符。标准 AE 使用两个反射：

S₀：对「全零态」反射，即 S₀ = I − 2|0⟩⟨0|，作用在计数寄存器与工作寄存器的联合空间上；
S_χ：对目标比特为 |1⟩ 的子空间反射，即 S_χ = I − 2|1⟩⟨1|_目标。

两者组合成 Grover 型算符 Q = −A S₀ A† S_χ。Q 在由 |ψ₀⟩|0⟩ 与 |ψ₁⟩|1⟩ 张成的二维子空间内是一个旋转，旋转角为 2θ，其中 sin²θ = a。于是 Q 的本征相位为 e^±2iθ，QPE 对 Q 做相位估计即可反解出 θ，进而得到 a = sin²θ。这正是「反射算符 + QPE」组合的数学本质：反射把振幅信息转成相位，QPE 把相位转成可读的比特串。

四、标准电路的整体组合

把上述部件按顺序拼接，就得到标准振幅估计电路：计数寄存器初始化为 |0⟩^⊗t，工作寄存器由 A 制备；随后施加受控 Q^{2^j}（j = 0…t−1），其中每个 Q 内部包含 A、S₀、A†、S_χ；接着对计数寄存器做逆量子傅里叶变换并测量，得到相位估计值；最后根据 θ 与 a 的关系还原振幅。下图给出这条电路的信号与算符流向，便于对照各部件的位置与接口。

五、组合时的关键约束与常见错误

把四个部件拼起来时，有几条约束必须同时满足，否则估计结果会系统性偏离：

相位约定一致：S₀ 与 S_χ 的负号约定必须与 Q 的本征相位推导一致。若把 Q 写成 A S₀ A† S_χ 而漏掉整体负号，QPE 读出的相位会整体平移 π，反解出的 θ 出现偏差。
受控幂次正确：第 j 个计数比特控制的是 Q^{2^j}，不是 Q^j。若误用 Q^j，QPE 的频率刻度会被压缩，估计值随 t 增大反而收敛到错误值。
A† 的逆序共轭：受控 A† 必须是 A 各基本门受控版本的逆序共轭，不能简单地把 A 的受控版本取逆后原序拼接。
工作寄存器复位：每轮 Q 迭代后工作寄存器应回到 A|0⟩ 附近，若复位不彻底，残余纠缠会污染后续迭代的相位累积。

排障要点：若 QPE 读出的相位分布出现双峰且间距异常，优先检查 S₀ 是否作用在「计数寄存器 + 工作寄存器」的联合全零态上，而不是只作用在计数寄存器上。这是标准 AE 电路最常见的实现错误之一。

六、一个可对照的示例参数

设待估计振幅 a = 0.25，则 θ = arcsin(0.5) = π/6，Q 的旋转角为 2θ = π/3。取计数寄存器 t = 4 比特，QPE 的相位分辨率为 1/16。理论上 QPE 会以较高概率读出 φ ≈ (π/3)/(2π) = 1/6，对应比特串约 0010 或 0011（取决于舍入与相位约定）。由 φ 反解 θ = πφ，得 θ ≈ π/6，a = sin²θ ≈ 0.25。若实际读出集中在 0000 附近，说明受控 Q 的幂次或 S₀ 的作用范围有误；若读出均匀散布，则多半是 A 的制备电路未正确把振幅编码进目标比特。这个示例仅用于说明参数关系与判读逻辑，实际数值会随具体实现的门分解与噪声而浮动。

掌握这条标准电路后，后续的迭代振幅估计（IAE）与最大似然振幅估计（MLAE）都可以看作在「Q 的调用次数」与「经典后处理」两个维度上对它的改造：前者用逐次比特判定的方式减少 Q 的调用，后者用似然函数替代逆 QFT 的硬判决。因此，把本节这条「A → 受控 Q → QFT† → 测量」的骨架吃透，是理解整个振幅估计家族的前提。