# 广义 QSP 与 QETU 详解：新一代相位因子框架

*前置阅读：[量子信号处理（QSP）详解](qsp-tutorial.md)、[量子奇异值变换（QSVT）详解](qsvt-tutorial.md)。本篇介绍 2022 年之后出现的三个变体——QETU、QPP、GQSP——它们分别拆掉了标准 QSP 的一根"紧箍咒"。*

**课程目标：**

1. 说清标准 QSP/QSVT 的三根"紧箍咒"：奇偶性约束、对块编码的依赖、单比特信号框架的刚性。
2. 理解 **QETU**（Dong–Lin–Tong 2022）：只用时间演化算符 $e^{-iHt}$ 这一种"黑盒"实现本征值变换，是早期容错时代的关键简化。
3. 理解 **QPP**（量子相位处理，Wang 等 2023）与 **GQSP**（Motlagh–Wiebe 2024）：摆脱奇偶性约束的一般化信号处理，及其带来的线路节省与"鲁棒相位估计"。
4. 会按"手里有什么黑盒、目标多项式什么性质"查表选择框架。
5. 了解这些变体与第 6 章（基态制备、能量估计）的接口。

## 1. 回顾：标准 QSP 的三根紧箍咒

[QSP 详解](qsp-tutorial.md)的可达集定理要求目标多项式：

**(a) 奇偶性约束**：$P$ 的奇偶性必须与层数 $d$ 相同。遇到无确定奇偶的目标（如 $e^{-i\tau x}$、截断的 $1/x$），标准做法是拆成偶、奇两支分别实现再用 LCU 合并——**代价近乎翻倍**（块编码归一化常数 $\alpha$ 变为两倍）。

**(b) 依赖块编码**：QSVT 把信号 $x$ 装进一个块编码 $U$（要求 $\|A\| \le \alpha$、构造 $(\alpha, a, 0)$ 块编码）。很多物理问题里，$H$ 的块编码构造（LCU、oracle 化）正是资源大头；而**时间演化** $e^{-iHt}$ 反而常有更便宜的直接实现（Trotter、模拟器原生指令）——信号处理却用不上它。

**(c) 单比特刚性**：信号处理固定为"单辅助比特上的 $S(\phi) = e^{i\phi Z}$ 与 $W(x)$ 交替"。这个形式优雅但约束了可实现多项式类与线路形状，也导致某些目标（一般复多项式、Laurent 多项式）需要间接绕路。

下面三个变体分别主打 (b)、(c)、(a)+(c)。

## 2. QETU：只用手里的演化算符

**Quantum Eigenvalue Transformation of Unitary matrices**（Dong, Lin & Tong, PRX Quantum 3, 040305, 2022）针对紧箍咒 (b)：

> **QETU 定理（表述一）**。设 $H$ 厄米、$U = e^{-iH/2}$ 可作受控调用。对任意次数 $\le d$ 的**偶实系数**多项式 $F$，$|F| \le 1$ 于 $[-1,1]$，存在对称相位序列 $(\phi_0, \phi_1, \ldots, \phi_d, \ldots, \phi_1, \phi_0)$，使得由 $d$ 次受控演化与相位门组成的单辅助比特线路满足
> $$ \langle 0|\, U_{\text{QETU}}\, |0\rangle \;=\; F\big(\cos(H/2)\big). $$

三个值得咀嚼的细节：

- **信号是 $\cos(H/2)$ 而非 $H$**：本征值 $\lambda$ 映到 $\cos(\lambda/2) \in [-1,1]$，信息经"余弦折叠"进入信号区间——目标函数需要先做这层变量替换再设计多项式。
- **只用偶实多项式看似很弱，实则够用**：$\cos$ 的偶次幂展开能表达 $\cos(H)$ 的任意多项式，而滤波、虚时间衰减 $e^{-\beta x}$ 的近似、以及"模糊二分"阈值函数都是实偶型目标。奇性部分确有损失——这正是后来 GQSP 补的位。
- **对称相位 + 单辅助比特**：$d$ 次受控 $e^{-iH/2}$ 调用、无块编码、无 QPE 的长寄存器，资源轮廓对早期容错硬件友好得多。

**应用（本篇看门戏）：基态制备与能量估计的"模糊二分法"**。取 $F$ 为过渡带在 $[\mu-\delta, \mu+\delta]$ 的阶跃函数近似，QETU 线路成为"本征值是否低于阈值 $\mu$"的过滤器。从 $\mu$ 取在整个谱的高处开始，二分式下降、逐层收窄 $\delta$，可以把初态投影到能量低于 $\mu$ 的子空间，最终收敛到基态——整个过程**不需要相位估计的长深线路**，且对初态重叠的依赖远比 QPE 温和。论文证明在"早期容错"资源模型下（逻辑比特贵、Clifford 便宜），基于 QETU 的基态制备+能量估计总代价优于基于 QPE 的方案。这条路线随后成为量子化学资源估计论文里的常客。

## 3. QPP 与 GQSP：拆掉奇偶性约束

**量子相位处理（QPP, Wang 等 2023）** 先把 QSP 的舞台从"块编码的信号比特"搬到"相位寄存器"上：对 $e^{-iH}$ 的本征相位 $\theta$ 直接实现 Laurent 多项式 $P(e^{i\theta})$。它给出了后来 GQSP 的核心机制——**每个信号门带两个相位自由度**（而不是一个 $S(\phi)$）——从而：

- 目标多项式不必满足奇偶性（$P$ 与 $Q$ 的约束由"单位圆上的模长条件"$|P(e^{i\theta})|^2 + |Q(e^{i\theta})|^2 = 1$ 取代）；
- 相位的数值求解（上篇的主题）也变得更稳定。

**广义 QSP（GQSP, Motlagh & Wiebe, PRX Quantum 5, 020368, 2024）** 把这套机制整理成与标准 QSP 平行的极简定理。定义含三参数的单比特旋转

$$
R(\theta, \phi, \lambda) \;=\; \begin{pmatrix} e^{i(\lambda+\phi)}\cos\theta & e^{i\phi}\sin\theta \\ e^{i\lambda}\sin\theta & -\cos\theta \end{pmatrix},
$$

> **GQSP 定理（表述一）**。设 $W$ 为任意酉算符。对任意次数 $\le d$ 的多项式 $P$（**无奇偶性要求**）满足 $|P(e^{i\vartheta})| \le 1$ 对 $W$ 的全部本征相位 $e^{i\vartheta}$ 成立，存在参数组 $\{(\theta_j, \phi_j)\}_{j=0}^d$ 与 $\lambda$，使得
> $$ \Big(\prod_{j=1}^{d} R(\theta_j, \phi_j, 0)\, W\Big)\, R(\theta_0, \phi_0, \lambda) \;=\; \begin{pmatrix} P(W) & \cdot \\ \cdot & \cdot \end{pmatrix}, $$
> 共使用 $W$ 恰好 $d$ 次。

对比标准 QSP 可达集：奇偶性条件消失了，代价是模长条件从 $|P(x)|\le1$（$x\in[-1,1]$）变为"单位圆上"的版本（对块编码型信号通过变量替换相连）。**标准 QSP 是 GQSP 的特例**（把 $R$ 退化成单参数 $Z$-相位即可）。工程收益立竿见影：

1. **免掉 LCU 合并**：$e^{-i\tau x}$ 这类无奇偶目标一步到位，块编码归一化常数不翻倍——复杂度常数直接减半量级；
2. **相位更好算**：GQSP 论文配套给出相位求解的稳定递推（结合上篇的确定性构造传统）；
3. **鲁棒相位估计**：GQSP 框架下可以构造对态制备误差、对相位寄存器长度误差都一阶不敏感的"鲁棒版 QPE"，替换第 3 章标准 QPE 后能降低容错线路的纠错开销。

## 4. 框架选择速查表

| 框架 | 信号黑盒 | 可实现目标 | 优点 | 典型用途 |
|---|---|---|---|---|
| QSP/QSVT | 块编码 $(\alpha,a,0)$ | 奇/偶 $P$，$\|P\|_{[-1,1]}\le1$ | 体系最完整、文献最多 | 矩阵函数、HHL、求逆 |
| QETU | 受控 $e^{-iHt}$ | 偶实 $F$（作用于 $\cos(H/2)$） | 无需块编码；早期容错友好 | 基态制备、能量估计、虚时间滤波 |
| QPP | $e^{-iH}$（相位寄存器视角） | Laurent 多项式 | 直接操作相位；熵估计 | von Neumann/Rényi 熵 |
| GQSP | 任意酉 $W$ | 任意 $P$（无奇偶），单位圆 $\|P\|\le1$ | 无奇偶约束、省 LCU、相位好算、鲁棒 QPE | 模拟、求逆的现代实现 |

读表方式：先问"黑盒"——手里有的是块编码还是演化算符；再问"目标"——奇偶性是否方便拆分；最后问"资源时代"——NISQ/早期容错优先 QETU 型，完全容错且追求最优常数时 GQSP 是新默认。四个框架共享同一套数学骨架（单比特信号 + 相位处理 + 可行对表示定理），因此上篇的相位求解技术与本篇的框架选择是正交的：**任何框架落地时都绕回"算相位、验相位"**。

## 5. 与本教程其余部分的接口

- **第 6 章**：基态制备（QETU 模糊二分）、能量估计与虚时间滤波（QETU/GQSP 化的 QPE）已经进入现代资源估计的主流水线；Gibbs 态制备（[Gibbs 态详解](../ch06-scientific-computing/gibbs-state.md)）同样可用 QETU 型 $e^{-\beta H}$ 滤波实现。
- **第 8 章**：变分算法的"后变分"替代品（不做变分优化、直接做谱滤波）本质就是这些框架的浅层应用。
- **开放问题**：GQSP 与双投影子结构（QSVT 的任意矩阵奇异值变换）的完全融合——"广义 QSVT"——仍是活跃方向；非厄米、多变量（M-QSP）信号处理的实用化同样在推进中。

## 本课总结

- 标准 QSP 的三根紧箍咒：奇偶约束（LCU 翻倍）、块编码依赖、单参数相位的刚性。
- QETU 只用受控 $e^{-iHt}$、单辅助比特实现偶实 $F(\cos(H/2))$，配套"模糊二分"做基态制备，早期容错下优于 QPE 路线。
- QPP/GQSP 用双参数（三参数）相位旋转拆除奇偶约束：任意 $\|P\|\le1$ 多项式、$d$ 次黑盒调用、相位更稳，还附带鲁棒相位估计。
- 选择框架 = 看黑盒（块编码 vs 演化）+ 看目标（奇偶可拆否）+ 看资源时代；四个框架共享相位求解技术栈。

## 习题

1. 验证 $R(\theta,\phi,\lambda)$ 是酉矩阵；找出使 GQSP 退化为标准 QSP 序列的参数选取（提示：让 $R$ 只依赖一个相位并取 $\theta$ 为固定值，与 $S(\phi) = e^{i\phi Z}$ 对齐）。
2. 设 $H$ 的谱在 $[-1,1]$、谱隙 $\Delta$。设计 QETU 滤波器：$F$ 取过渡带 $[-\tfrac{\Delta}{4}, \tfrac{\Delta}{4}]$ 内由 $1$ 降到 $\approx0$ 的偶多项式（画出草图即可），说明 $\langle\psi_0|F^2(\cos(H/2))|\psi_0\rangle$（$\psi_0$ 为基态与第一激发态的叠加）如何随 $F$ 的陡度变化——由此说明"模糊二分"每一步需要多少次演化调用。
3. 用 GQSP 实现哈密顿量模拟 $e^{-i\tau H}$（块编码 $\alpha = \|H\|$）：写出目标 $P$、验证无奇偶约束带来的 LCU 节省（与 [QSP 详解](qsp-tutorial.md)第 4 节的奇偶拆分对比），估算节省的归一化常数。
4. 思考题：QETU 的"偶实多项式"限制与 GQSP 的"无奇偶限制"是否矛盾？（提示：两者黑盒不同——前者只有 $e^{-iHt}$，后者允许把 $H$ 的信息装进任意酉 $W$；把 QETU 的信号 $\cos(H/2)$ 与 GQSP 的单位圆信号 $e^{-i\vartheta}$ 通过 $x = \cos(\vartheta/2)$ 联系起来。）