📚 课程目录第2章量子振幅估计的数学基础

2. 量子振幅估计的数学基础

量子振幅估计（Quantum Amplitude Estimation, QAE）之所以能在某些问题上取得相对经典蒙特卡洛的二次加速，根源在于它把“概率”这个统计量，巧妙地转化成了“相位”这个可由量子相位估计精确读取的量。理解这一转化链条，是掌握 QAE 全部后续内容的前提。本节围绕三个核心对象展开：振幅 a、概率 p 与相位 θ，以及承载它们的算符 A 与好态子空间。

一、从算符 A 到好态子空间

设我们有一个 n 量子比特的寄存器，以及一个辅助量子比特（ancilla）。算符 A 是一个作用在 n+1 个量子比特上的酉算符，它把初态 |0⟩_n+1 变换为：

A|0⟩_{n+1} = √(1−a²)|ψ₀⟩|0⟩ + a|ψ₁⟩|1⟩

这里 |ψ₀⟩ 与 |ψ₁⟩ 是 n 量子比特寄存器上的两个归一化量子态，分别对应辅助比特为 |0⟩ 和 |1⟩ 的分支。系数 a ∈ [0,1] 称为振幅，它正是我们想要估计的目标量。这个式子的结构非常关键：算符 A 把整个希尔伯特空间分解为两个正交子空间的直和——由 {|ψ₀⟩|0⟩} 张成的“坏态子空间”，和由 {|ψ₁⟩|1⟩} 张成的“好态子空间”。所谓“好态”，就是辅助比特处于 |1⟩、寄存器处于 |ψ₁⟩ 的那一支，它对应我们关心的“成功事件”。

由此，好态子空间可形式化定义为 𝒮_good = span{|ψ₁⟩|1⟩}，其正交补 𝒮_bad = span{|ψ₀⟩|0⟩}。二者共同张成 A 作用所涉及的一个二维子空间，记作 𝒮 = span{A|0⟩, |ψ₁⟩|1⟩}。QAE 的全部动力学都被限制在这个二维子空间内，这使得原本高维的问题可以降维成平面上的旋转来分析。

二、振幅、概率与相位的三角关系

定义相位 θ 满足：

a = sin(θ),   θ ∈ [0, π/2]

于是测量辅助比特得到 |1⟩ 的概率为 p = a² = sin²(θ)。这就是振幅与概率的第一层关系：概率是振幅的平方，而振幅是概率的算术平方根。QAE 的目标是估计 a（或等价地 p），但它实际通过量子相位估计去估计的是 θ，再由 a = sin(θ) 反推。这条“概率 → 振幅 → 相位 → 估计值 → 反三角 → 概率”的链条，是理解 QAE 精度与误差传播的主线。

为什么绕道相位反而更优？因为经典方法要靠反复采样统计频率来逼近 p，误差按 1/√N 衰减（N 为采样次数）；而相位估计利用量子干涉，让 θ 的信息以旋转角的形式累积，误差按 1/N 衰减，从而实现二次加速。代价是需要构造受控的 Q 算符并做量子傅里叶逆变换。

三、算符 Q 与相位累积

仅有 A 还不够，相位估计需要一个能反复施加、使相位线性累积的酉算符。定义：

Q = A·S₀·A⁻¹·S_χ

其中 S_χ 对好态子空间施加 −1 相位（即 S_χ = I − 2|ψ₁⟩⟨ψ₁|⊗|1⟩⟨1|），S₀ 对 |0⟩_n+1 施加 −1 相位。可以证明，Q 在二维子空间 𝒮 上等价于一个旋转矩阵：

Q = [[cos 2θ, −sin 2θ],
     [sin 2θ,  cos 2θ]]

也就是说，Q 把 A|0⟩ 绕子空间原点旋转 2θ。施加 k 次后，状态相对好态轴的夹角变为 (2k+1)θ。量子相位估计正是通过读取这个旋转的本征相位 2θ，再除以 2 得到 θ。这里有一个容易混淆的点：Q 的本征相位是 e^±i2θ，所以 QPE 直接给出的是 2θ 的二进制近似，必须记得除以 2 才能还原 θ，进而 a = sin θ。若忘记这一步，振幅会被系统性高估。

四、一个可手算的教学示例

假设某算符 A 使得真实振幅 a = 0.5，则 p = 0.25，θ = arcsin(0.5) = π/6 ≈ 0.5236 rad。若使用 m = 4 个相位估计比特，QPE 把 2θ/2π = (π/3)/(2π) = 1/6 量化为 4 位二进制小数。1/6 ≈ 0.1667，最接近的 4 位小数为 3/16 = 0.1875（对应二进制 0.0011）。于是估计的 2θ̃ = 2π·(3/16) = 3π/8，θ̃ = 3π/16 ≈ 0.5890，ã = sin(0.5890) ≈ 0.5556，p̃ ≈ 0.3086。与真值 0.25 相比，绝对误差约 0.0586。这个例子清楚展示了：相位量化误差会经 sin 函数非线性地传递到概率上，在 θ 接近 π/2（即 p 接近 1）时，同样的相位误差引起的概率误差会被压缩；而在 θ 接近 0 时则被放大。设计实验时应据此评估所需比特数。

判断精度是否够用的实用方法：先由目标概率误差 ε_p 反推相位误差容限。由 dp/dθ = 2 sinθ cosθ = sin 2θ，得 ε_θ ≈ ε_p / |sin 2θ|。当 p 很小或接近 1 时，|sin 2θ| 小，需要更多相位比特；当 p ≈ 0.5 时最宽松。

五、常见误区与检查要点

混淆 a 与 p：QAE 的原始输出是振幅 a 的估计，若业务指标是概率，必须再平方；反之若已知概率要构造 A，则需取平方根。
忘记 Q 的 2θ 因子：QPE 读的是 2θ，还原 θ 要除以 2，这是最高频的实现错误。
好态定义与测量算符不一致：好态子空间必须与最终测量算符的本征子空间严格对应，否则估计的是另一个量。
θ 的取值范围：约定 θ ∈ [0, π/2] 可保证 a = sinθ 单调、反解唯一；若允许 θ 超出该区间，arcsin 会出现多值歧义。

把上述关系串起来看：算符 A 负责把目标概率编码进振幅，好态子空间界定了“成功”的判据，算符 Q 把振幅转化为可累积的旋转相位 2θ，量子相位估计读出 θ，最后由 a = sinθ、p = a² 还原业务量。这条链路中任何一环的定义偏差，都会以系统误差的形式出现在最终估计里，因此在实现 QAE 之前，先把 A、好态、Q 与 θ 的约定逐一对齐，是最值得投入的准备工作。