# 量子经典混合神经网络详解：原理、架构与前沿

量子经典混合神经网络（Quantum-Classical Hybrid Neural Network）是把变分量子电路（Variational Quantum Circuit, VQC）嵌入经典深度学习架构的机器学习范式。它利用量子电路的高维希尔伯特空间作为特征映射，结合经典神经网络的表达力和成熟的优化技术，是 NISQ（Noisy Intermediate-Scale Quantum，近期含噪中等规模量子）时代量子机器学习最具前景的研究方向之一。

:::{admonition} 本课知识点
:class: tip

1. **[混合架构的动机与量子的可能性](#motivation-quantum-possibilities)**——能列出经典深度学习面临的维度、量子数据与结构先验三类挑战，写出量子电路提供的三种可能性，并解释"指数维状态空间"为何尚不等于已被证明的优势。
2. **[表达力的傅里叶刻画](#fourier-expressivity)**——能推导输出 $f(x)$ 是频率集合有限的三角多项式，由编码门生成元的谱计算可用频率，并比较扩大频带的两种办法。
3. **[参数化电路与三种数据编码](#data-encoding)**——能写出角度、振幅与 IQP 编码的公式，比较各自的量子比特数与制备门开销，并说明 IQP 核含两两乘积特征的来源。
4. **[经典-量子前馈网络与非线性来源](#hybrid-feedforward-network)**——能依次写出编码层、量子层、测量层与输出层的公式链，并解释非线性来自读出对参数的正弦型依赖而非测量本身。
5. **[QCNN 与序列、生成变体](#qcnn-variants)**——能描述 QCNN 的卷积-池化流程并解释参数共享为何使参数量随层数线性增长，说明 QRNN 隐藏态与 QGAN 训练的要点。
6. **[混合反向传播与参数平移](#backprop-chain)**——能逐环推导损失对经典输出层、量子参数与编码层参数的梯度，写出 $G^2=I$ 与一般谱隙两种情形的参数平移公式，并核算电路执行的规模。
7. **[最小示例的端到端训练](#pennylane-example)**——能解释示例代码中经典层与量子层的衔接方式，说明两类参数如何在同一优化步中分别由参数平移与自动微分更新，并通过修改缩放试验解释收敛变化。
8. **[Barren Plateaus 与泛化上界](#barren-plateaus)**——能陈述梯度方差随量子比特数指数衰减的结论并列举缓解策略，写出泛化差距的有效维度上界并解释其含义。

:::

## 问题背景

(motivation-quantum-possibilities)=
### 为什么需要量子神经网络

经典深度学习在高维数据处理上已经取得巨大成功，但仍然面临以下挑战：

- 维度问题：若想显式使用指数维的特征映射（例如高阶多项式特征），经典模型的参数量或计算量随特征维数快速增长；
- 量子数据：当训练对象本身是量子态或量子过程时，经典网络必须先做层析（tomography，代价随量子比特数指数增长）才能拿到经典描述，而量子模型可以直接消费量子态；
- 结构先验：某些物理问题的对称性与纠缠结构，用量子电路作为模型先验可能比手工设计的经典结构更自然。

作为对比，量子电路提供的可能性包括：

- 指数维状态空间：$n$ 个量子比特携带 $2^n$ 维复向量，用 $O(n)$ 个量子比特与多项式门数即可制备其中多项式复杂的态；
- 纠缠：电路天然产生非局域关联，这是经典概率模型难以直接表示的；
- 干涉：输出的期望值是各路径振幅的相干叠加，这既是表达力的来源，也是分析（与攻击）的难点。

需要强调的是，这些"可能性"目前多数没有在可证明的意义上转化为对经典模型的优势；本文末尾"当前局限性"一节会回到这一点。

(fourier-expressivity)=
### 表达力的傅里叶视角：一个可推导的刻画

在正式介绍架构之前，我们先给出一个可以完整推导的 expressivity（表达力）刻画，它回答"一个固定的含数据电路到底能表示什么样的函数"。设模型输出为

$$f(x) = \langle 0|\, U^\dagger(x)\, M\, U(x)\, |0\rangle,$$

其中数据 $x$ 通过若干门 $e^{-ix G^{(k)}/2}$（$k = 1, \ldots, K$）进入电路，$M$ 是被测的可观测量，电路中其余门与 $x$ 无关（为叙述简单，我们把所有含 $x$ 的门集中写出并在其间插入与 $x$ 无关的酉阵 $W_j$；一般布局只是多几个 $W$ 因子，结论不变）：

$$U(x) = W_K\, e^{-ixG^{(K)}/2} \cdots W_1\, e^{-ixG^{(1)}/2}\, W_0.$$

对每个生成元做谱分解 $G^{(k)} = \sum_\lambda \lambda\,|\lambda^{(k)}\rangle\langle\lambda^{(k)}|$，则每个含参门的指数也是可分解的：

$$e^{-ixG^{(k)}/2} = \sum_{\lambda \in \mathrm{spec}(G^{(k)})} e^{-ix\lambda/2}\, |\lambda^{(k)}\rangle\langle\lambda^{(k)}|.$$

把 $U(x)$ 与 $U^\dagger(x)$ 的乘积按这些投影算符逐项展开：每一项是一个与 $x$ 无关的矩阵元乘上一个标量相位 $e^{ix(\sum_k \lambda'_k - \lambda_k)/2}$（$U^\dagger$ 中的门贡献共轭相位 $e^{+ix\lambda'/2}$）。于是

$$f(x) = \sum_{\vec\lambda, \vec\lambda'} C_{\vec\lambda, \vec\lambda'}\; e^{\,i\omega(\vec\lambda,\vec\lambda')\,x}, \qquad \omega(\vec\lambda, \vec\lambda') = \frac{1}{2}\sum_{k=1}^{K}\big(\lambda'_k - \lambda_k\big),$$

其中系数 $C_{\vec\lambda,\vec\lambda'}$ 完全由 $W_j$、初态与 $M$ 决定。结论：$f(x)$ 是有限个频率 $\omega \in \Omega = \{\frac12\sum_k(\lambda'_k - \lambda_k)\}$ 的三角多项式，**可用的频率集合完全由编码门的生成元谱决定，线性读出层只能组合这些频率成分，不能创造新频率**。

我们用一个最小例子核对：单个编码门 $R_y(x) = e^{-ixY/2}$，$\mathrm{spec}(Y) = \{+1, -1\}$，频率取值 $(\lambda' - \lambda)/2 \in \{0, \pm 1\}$，因此 $f(x) = a\sin x + b\cos x + C$——这与第 2 章《参数平移法则》附录中对单参数门导出的正弦响应完全一致。要提高频率分辨率或扩大频带，有两个办法：使用谱隙更大的生成元（例如两比特生成元 $Z_1 + Z_2$，谱为 $\{-2, 0, +2\}$，频率可达 $\pm 2$），或做数据重上传（data re-uploading）——把同一个 $x$ 编码进多个门，频率集合变为多次差组合，频带按重上传次数倍增，代价是更深的电路。

## 核心架构：参数化量子电路（PQC）

### 变分量子电路

变分量子电路是量子神经网络的基本构建块：

$$|\psi(\vec{\theta})\rangle = U_L(\theta_L)\cdot U_{L-1}(\theta_{L-1}) \cdots U_1(\theta_1)\, |0\rangle^{\otimes n},$$

其中每一层由单比特旋转和两比特纠缠门组成，例如

$$U_l(\theta_l) = \prod_{(i,j)} \mathrm{CNOT}_{ij} \cdot \prod_i R_y(\theta_{i,l})\, R_z(\theta'_{i,l}).$$

注意参数以 $R_y(\theta) = e^{-i\theta Y/2}$ 的方式进入门（本教程统一约定），这对后文梯度公式的形式至关重要。

(data-encoding)=
### 数据编码

将经典数据 $\mathbf{x} \in \mathbb{R}^d$ 编码为量子态有三种常用方式。

角度编码（Angle Encoding）把每个特征当作一个旋转角：

$$|\mathbf{x}\rangle = \bigotimes_{i=1}^{d} R_y(x_i)\, |0\rangle,$$

需要 $n \ge d$ 个量子比特，每比特一个旋转门，编码开销 $O(d)$。

振幅编码（Amplitude Encoding）把数据写进幅度向量：

$$|\mathbf{x}\rangle = \frac{1}{\|\mathbf{x}\|} \sum_{i=0}^{N-1} x_i |i\rangle, \qquad N = 2^n \ge d,$$

只需 $n = \lceil \log_2 d \rceil$ 个量子比特，但制备一般态需要 $O(N)$ 个门，指数的制备成本往往抵消比特数上的节省。

IQP 编码（Instantaneous Quantum Polynomial）用对角的相位门编码，是量子核方法中常用的高维特征映射：

$$|\mathbf{x}\rangle = \Big(\prod_{i<j} e^{i x_i x_j Z_i Z_j/2}\Big)\Big(\prod_i e^{i x_i Z_i /2}\Big) H^{\otimes n}|0\rangle^{\otimes n}.$$

所有因子都是对角的（因此互相内对易，"instantaneous"由此得名）。其中的两比特因子与受控相位门只差单比特相位：记 $\mathrm{CP}(\theta) = \mathrm{diag}(1,1,1,e^{i\theta}) = e^{i\theta |11\rangle\langle 11|}$，利用 $|11\rangle\langle 11| = \frac14(I - Z_i)(I - Z_j)$ 与指数的逐因子拆开（各因子互相对易），可得

$$\mathrm{CP}(\theta) = e^{i\theta/4}\, e^{-i\theta Z_i/4}\, e^{-i\theta Z_j/4}\, e^{i\theta Z_iZ_j/4},$$

因此 $e^{ix_ix_jZ_iZ_j/2} = e^{-i\theta/4}\,e^{i\theta Z_i/4}e^{i\theta Z_j/4}\,\mathrm{CP}(\theta)$（取 $\theta = 2x_ix_j$）：IQP 的两比特因子确实可以用一个受控相位门加上可吸收的单比特 $R_z$ 实现。由于二次相位 $e^{ix_ix_j}$ 的存在，IQP 编码的核 $|\langle\mathbf{x}|\mathbf{x}'\rangle|^2$ 含有特征的所有两两乘积 $x_ix_j$，这就是它作为高维特征映射的来源。

## 混合架构类型

(hybrid-feedforward-network)=
### 架构一：经典-量子前馈网络

```
输入 x → 经典编码层 → 量子电路(θ) → 经典测量层 → 输出 y
```

我们把每一段的公式明确写出来，这组公式就是后文反向传播推导的对象。

1. 经典编码层：线性变换加缩放，把 $d$ 维输入映射为 $n$ 个旋转角
   $$\vec{\phi} = \sigma(W_1 \mathbf{x} + \mathbf{b}_1), \qquad W_1 \in \mathbb{R}^{n\times d},$$
   其中 $\sigma$ 是逐元素缩放函数（例如 $\tanh$ 后乘 $\pi$，把角度限制在 $(-\pi, \pi)$）。
2. 量子电路层：角度编码加变分层
   $$\rho = U(\vec{\theta})\, \Big(\bigotimes_i R_y(\phi_i)\Big)\, |0\rangle\langle 0|^{\otimes n}\, \Big(\bigotimes_i R_y(\phi_i)\Big)^\dagger U^\dagger(\vec{\theta}).$$
3. 测量层：对每个量子比特测 $Z$ 期望
   $$z_i = \mathrm{Tr}\big[\rho\, Z_i\big] = \langle Z_i\rangle \in [-1, 1].$$
4. 经典输出层：线性读出
   $$\hat{y} = \mathbf{w}_2^\top \vec{z} + b_2 = \sum_i w_{2,i}\, z_i + b_2.$$

对回归任务，损失取均方误差 $\mathcal{L} = \frac12(\hat y - t)^2$；对 $C$ 分类任务，把读出改为 $C$ 个量子比特的期望并接 softmax 与交叉熵，推导完全平行。

**期望值读出在什么意义上是"激活函数"？** 我们需要澄清一个容易误解的点：测量映射 $\rho \mapsto \mathrm{Tr}[Z_i\rho]$ 对 $\rho$ 是线性的，量子层并不在"对态做非线性变换"的意义上模仿经典激活函数。非线性来自另一个地方——态对参数的依赖。以单个参数门为例（这正是第 2 章附录的正弦响应引理）：若 $U(\theta) = e^{-i\theta G/2}$ 且 $G^2 = I$，则对任意初态与观测量

$$z(\theta) = a\sin\theta + b\cos\theta + C,$$

系数与 $\theta$ 无关。也就是说，读出值是参数的正弦型非线性函数。混合模型中真正的分工是：量子电路实现一个（对参数、对输入都是非线性的）特征映射 $\mathbf{x} \mapsto \vec z(\mathbf{x}, \vec\theta) \in [-1,1]^n$，经典层承担可微的线性组合与损失计算。这与经典网络中"固定的非线性激活 + 可学习的线性变换"的组合在功能上是对应的。

(qcnn-variants)=
### 架构二：量子卷积神经网络（QCNN）

量子卷积神经网络模仿经典 CNN 的层级结构：

1. 量子卷积层：在相邻量子比特对上施加同一个参数化两比特门（参数跨位置共享，这是"卷积"的含义）；
2. 量子池化层：测量部分量子比特，根据测量结果对其余量子比特做受控旋转（或做后选择），有效自由度减半；
3. 重复卷积-池化直到只剩 $O(1)$ 个量子比特；
4. 最终测量得到输出。

优势：参数随层数线性增长（每层 $O(1)$ 个共享参数），电路深度 $O(\log n)$ 级别的纠缠结构由池化自然产生；由于每层只涉及局域自由度且参数量小，QCNN 的训练对 Barren Plateaus 相对不敏感（见理论分析一节）。

### 架构三：量子循环神经网络（QRNN）

把 PQC 用于序列数据处理，隐藏态是量子态：

$$|\psi_t\rangle = \mathrm{PQC}(\vec{\theta};\, |\psi_{t-1}\rangle, |x_t\rangle),$$

其中 $|x_t\rangle$ 是 $t$ 时刻的输入态，$|\psi_{t-1}\rangle$ 是上一时刻的隐藏态，PQC 在每一步共享同一组参数 $\vec\theta$。量子 LSTM 用量子电路实现门控机制，让隐藏态保持相干性；主要工程困难在于隐藏态不能被测量破坏（否则退化为经典隐藏态），因此读出要用弱测量或辅助比特。

### 架构四：量子生成对抗网络（QGAN）

生成器是量子电路 $G(\vec\theta)$，从 $|0\rangle^{\otimes n}$（加上噪声输入）生成样本；判别器可以是经典网络或量子电路，负责区分真实数据与生成样本。训练时交替优化生成器与判别器的参数（minimax 博弈），与经典 GAN 的框架相同，区别在于生成器的样本空间是 $2^n$ 个计算基态上的一个可调分布。

## 训练方法：前向与反向传播的完整推导

我们以架构一（单输出回归）为例，把训练一轮所需的全部公式链逐步写出。记可学习参数为经典层参数 $\{W_1, \mathbf{b}_1, \mathbf{w}_2, b_2\}$ 与量子层参数 $\vec{\theta}$。

### 前向传播链

前向传播就是把架构一的四个步骤依次执行一遍：

$$\mathbf{x} \;\xrightarrow{\;W_1, \mathbf{b}_1\;}\; \vec{\phi} = \sigma(W_1\mathbf{x} + \mathbf{b}_1) \;\xrightarrow{\;R_y\;\circ\; U(\vec\theta)\;}\; \rho \;\xrightarrow{\;\langle Z_i\rangle\;}\; \vec{z} \;\xrightarrow{\;\mathbf{w}_2, b_2\;}\; \hat{y} \;\xrightarrow{\;\mathcal{L}\;}\; \tfrac12(\hat y - t)^2.$$

其中 $\vec z$ 的每个分量 $z_i = \langle 0|V^\dagger Z_i V|0\rangle$，$V$ 是整个电路（编码门加变分门）的乘积。

(backprop-chain)=
### 反向传播：从损失到各层参数

反向传播就是链式法则的逐层应用，我们把每一环都显式写出。

第一环：损失对输出的导数。

$$\frac{\partial \mathcal{L}}{\partial \hat{y}} = \hat{y} - t.$$

第二环：经典输出层的参数梯度。由 $\hat y = \sum_i w_{2,i}z_i + b_2$ 对各参数求偏导（$\partial \hat y/\partial w_{2,i} = z_i$、$\partial \hat y/\partial b_2 = 1$），代入链式法则：

$$\frac{\partial \mathcal{L}}{\partial w_{2,i}} = (\hat y - t)\, z_i, \qquad \frac{\partial \mathcal{L}}{\partial b_2} = \hat y - t.$$

第三环：把误差信号传到读出值。定义 $\delta_i := \frac{\partial\mathcal{L}}{\partial z_i}$。由 $\hat y = \sum_j w_{2,j}z_j + b_2$ 得 $\partial \hat y/\partial z_i = w_{2,i}$，于是

$$\delta_i = (\hat y - t)\, w_{2,i}.$$

第四环：损失对量子参数的梯度。$z_i$ 通过态 $\rho$ 依赖于每个 $\theta_k$，链式法则按求和展开：

$$\frac{\partial \mathcal{L}}{\partial \theta_k} = \sum_{i=1}^{n} \delta_i\, \frac{\partial z_i}{\partial \theta_k}.$$

其中每个 $\partial z_i/\partial\theta_k$ 是单个参数门的响应导数，由参数平移规则（Parameter-Shift Rule, PSR）精确计算。按本教程的约定 $U_k(\theta_k) = e^{-i\theta_k G_k/2}$：当生成元满足 $G_k^2 = I$（涵盖 $R_x, R_y, R_z$、$e^{-i\theta Z_iZ_j/2}$ 等常用门）时，

$$\frac{\partial z_i}{\partial \theta_k} = \frac{1}{2}\Big[z_i\Big(\theta_k + \frac{\pi}{2}\Big) - z_i\Big(\theta_k - \frac{\pi}{2}\Big)\Big];$$

若 $G_k$ 恰有两个不同本征值、谱隙 $\Delta = \lambda_1 - \lambda_2$，则推广为

$$\frac{\partial z_i}{\partial \theta_k} = \frac{\Delta}{4}\Big[z_i\Big(\theta_k + \frac{\pi}{\Delta}\Big) - z_i\Big(\theta_k - \frac{\pi}{\Delta}\Big)\Big].$$

两者的推导（期望值的正弦响应 $z_i = a\sin\theta_k + b\cos\theta_k + C$，以及用两次平移点的差分重建导数）见[第 2 章附录：参数平移法则](../ch02-quantum-nn/parameter-shift-appendix.md)。

第五环：梯度继续穿过编码层回到经典参数。角度 $\phi_j$ 也是量子门的参数（$R_y(\phi_j)$ 的生成元为 $Y$，谱 $\{\pm1\}$），同样用 PSR 得到 $\partial z_i/\partial \phi_j = \frac12[z_i(\phi_j + \frac\pi2) - z_i(\phi_j - \frac\pi2)]$，于是

$$\frac{\partial \mathcal{L}}{\partial \phi_j} = \sum_{i=1}^{n} \delta_i\, \frac{\partial z_i}{\partial \phi_j},$$

再对 $\vec\phi = \sigma(W_1\mathbf{x} + \mathbf{b}_1)$ 用一次链式法则（$\sigma$ 逐元素作用，导数 $\sigma'$ 逐元素相乘）：

$$\frac{\partial \mathcal{L}}{\partial W_{1,j\mu}} = \frac{\partial \mathcal{L}}{\partial \phi_j}\, \sigma'(\cdot)\, x_\mu, \qquad \frac{\partial \mathcal{L}}{\partial b_{1,j}} = \frac{\partial \mathcal{L}}{\partial \phi_j}\, \sigma'(\cdot),$$

其中 $(\cdot) = (W_1\mathbf{x} + \mathbf{b}_1)_j$。至此整条链闭合：损失对**所有**经典参数与量子参数的梯度都有了显式表达式。

执行成本可以这样核算：$\frac{\partial\mathcal{L}}{\partial\theta_k}$ 的每一次平移评估需要对整个批量数据求和，即 $P$ 个量子参数与 $n + 1$ 个读出相关角度共 $O(P + n)$ 个待估导数，每个需要 2 次电路执行（每次执行含统计测量次数），单样本单参数的电路执行数是常数，总执行次数 $O((P + n)\cdot N_{\text{batch}})$；经典参数的梯度在框架内自动微分无需额外量子执行。这与经典反向传播"每参数常数次前向"的代价结构同阶，只是这里的"前向"是真实量子电路的重复执行。

### 参数平移与自动微分的分工

在 PennyLane、TensorFlow Quantum、JAX 等框架中，上述链式法则由自动微分系统自动完成：量子参数的局部导数由 PSR 给出（解析、无偏差，且天然兼容硬件上的有限次数测量估计），经典参数的导数由标准反向传播给出，两者在"角度 $\vec\phi$"这一接口处拼接。用户只需要把量子电路写成 QNode（或 Keras 层），即可端到端训练。

常用的经典优化器：

- Adam：自适应学习率，混合模型默认首选；
- 自然梯度（Natural Gradient）：用量子 Fisher 信息矩阵预条件梯度，考虑参数空间的几何，小模型上收敛更快；
- SPSA（同时扰动随机近似）：每次迭代只需 2 次电路执行即可得到一个随机梯度估计，对测量噪声与门噪声鲁棒。

## 量子经典混合训练范式

### 范式一：前向-反向混合训练

即上文完整推导的模式：每轮前向执行电路得到 $\vec z$ 与损失，反向时量子参数走参数平移、经典参数走反向传播。这是 PennyLane/TensorFlow Quantum 的默认范式。

优势：梯度解析精确；代价：每个参数每批需要 2 次电路执行，总电路执行数随参数数线性增长。

### 范式二：基于样本复用的估计（classical shadows）

同一批测量样本不止能估计一个期望值：经典影子（classical shadow）技术用随机基测量 + 经典后处理，从 $O(\log M/\varepsilon^2)$ 次测量同时估计 $M$ 个可观测量。由于混合模型的读出与梯度都只涉及泡利字符串期望，影子方法可以显著减少总测量次数，把"每参数 2 次执行"的常数摊薄。

### 范式三：核方法（Kernel Methods）

把电路固定（不再训练 $\vec\theta$），只把它当作特征映射，定义量子核

$$K(\mathbf{x}, \mathbf{x}') = \big|\langle \psi(\mathbf{x}) | \psi(\mathbf{x}')\rangle\big|^2,$$

再用经典的 SVM 或高斯过程在核上训练。核函数的估计可以通过交换测试（SWAP test）或对角元素的重构在硬件上完成。

优势：完全绕开变分优化，因此不受 Barren Plateaus 影响，理论性质（如核的泛化界）可以直接搬用；代价：核矩阵需要 $O(N_{\text{train}}^2)$ 次电路评估，且模型表达力受限于固定电路的核。

(pennylane-example)=
## 一个最小可运行示例

下面用 PennyLane 实现架构一：2 维经典输入、4 个量子比特、期望值读出、线性输出层，端到端训练一个小回归任务。量子参数的梯度由 PSR 自动计算，经典层参数的梯度由 autograd 反向传播给出，两者在前文推导的"角度"接口处自动衔接。

```python
import pennylane as qml
from pennylane import numpy as np

n_qubits = 4
dev = qml.device("default.qubit", wires=n_qubits)

def circuit(theta, angles):
    """角度编码层 + 两层硬件高效拟设"""
    for i in range(n_qubits):
        qml.RY(angles[i], wires=i)          # 数据编码：phi_i 来自经典层
    for l in range(2):                       # 变分层：单比特旋转 + 环式 CNOT
        for i in range(n_qubits):
            qml.RY(theta[l, i], wires=i)
            qml.RZ(theta[l, i + n_qubits], wires=i)
        for i in range(n_qubits):
            qml.CNOT(wires=[i, (i + 1) % n_qubits])

@qml.qnode(dev)
def quantum_layer(theta, angles):
    circuit(theta, angles)
    return [qml.expval(qml.PauliZ(i)) for i in range(n_qubits)]

def model(classical_params, theta, x):
    """经典层 -> 量子层 -> 经典读出，即前向传播链"""
    W1 = classical_params[: n_qubits * 2].reshape((n_qubits, 2))
    b1 = classical_params[n_qubits * 2 : n_qubits * 3]
    w2 = classical_params[n_qubits * 3 : n_qubits * 4]
    b2 = classical_params[-1]
    angles = np.tanh(W1 @ x + b1) * np.pi    # 编码层：线性 + 缩放到 (-pi, pi)
    z = quantum_layer(theta, angles)          # 量子层：z_i = <Z_i>
    return np.dot(w2, z) + b2                 # 输出层：线性读出

def loss(classical_params, theta, x, y):
    return (model(classical_params, theta, x) - y) ** 2

# 训练数据：4 个样本的小回归任务
X = np.array([[0.5, -1.0], [-0.3, 0.8], [1.2, 0.4], [-1.5, -0.6]], requires_grad=False)
Y = np.array([0.7, -0.2, 1.0, -0.8], requires_grad=False)

# 参数初始化：经典层（W1/b1/w2/b2 展平）与量子层
np.random.seed(0)
classical_params = np.random.uniform(-1, 1, n_qubits * 4 + 1, requires_grad=True)
theta = np.random.uniform(0, 2 * np.pi, (2, 2 * n_qubits), requires_grad=True)

def cost(classical_params, theta):
    return np.mean([loss(classical_params, theta, X[i], Y[i]) for i in range(len(X))])

opt = qml.AdamOptimizer(stepsize=0.05)
for step in range(200):
    classical_params, theta = opt.step(cost, classical_params, theta)
    if (step + 1) % 40 == 0:
        print(f"step {step+1:3d}: loss = {cost(classical_params, theta):.6f}")
```

运行这段代码，损失会随训练下降。值得对照前文的推导检查两点：其一，`opt.step(cost, classical_params, theta)` 同时更新两类参数——量子参数走 PSR、经典参数走反向传播，这正是"混合"二字的实现；其二，把 `np.tanh(...) * np.pi` 中的缩放去掉后收敛通常变慢，因为角度会漂出正弦响应的一个周期，梯度的有效幅度变小。

## 理论分析

### 表达力（Expressibility）

表达力刻画一个电路族在参数扫动时能覆盖多大的态空间。常用的度量是：随机采样参数得到的态分布与 Haar 均匀分布之间的偏差（例如用态的酉不变函数的分布距离衡量）。经验规律是：表达力随层数增加而快速饱和，深层硬件高效拟设的参数化行为接近 Haar 随机酉；太浅的电路表达力不足（如上文傅里叶刻画所示，可用频率太少），太深的电路则落入下述 Barren Plateaus 的陷阱。表达力与可训练性之间存在张力，需要按任务折中。

(barren-plateaus)=
### Barren Plateaus

**定理（McClean et al., 2018，大意）** 对深度为 $\mathrm{poly}(n)$ 的全局随机式硬件高效拟设，损失函数关于任一参数的偏导数满足

$$\mathbb{E}\left[\frac{\partial \mathcal{L}}{\partial \theta_k}\right] = 0, \qquad \mathrm{Var}\left[\frac{\partial \mathcal{L}}{\partial \theta_k}\right] = O(2^{-n}),$$

即梯度幅度随量子比特数指数衰减，优化景观指数平坦，从随机初始化出发几乎不可能训练。

常用的缓解策略：

1. 局域代价函数：让损失只由少数量子比特上的测量组成（如 QCNN 的逐层读出），局域观测的方差不随 $n$ 指数衰减；
2. 结构化初始化：从零角度或小的随机扰动出发（正弦响应 $a\sin\theta + b\cos\theta + C$ 在 $\theta = 0$ 附近的梯度由系数 $a$ 决定，物理动机的电路常有多项式级的 $a$），或用问题对称性约束初值；
3. 限制深度与作用范围：逐层训练、块.identity 初始化；
4. 换训练范式：如核方法完全绕开参数优化（见范式三）。

### 泛化能力

**定理（Abbas et al., 2021，大意）** 量子神经网络在训练集与测试集上的期望损失差距可以上界为

$$\mathbb{E}\big[\mathcal{L}_{\text{test}} - \mathcal{L}_{\text{train}}\big] = O\Big(\sqrt{d_{\text{eff}} / N_{\text{train}}}\Big),$$

其中 $d_{\text{eff}}$ 是模型的有效维度（由 Fisher 信息谱决定，不超过参数数量 $P$），$N_{\text{train}}$ 是训练样本数。这与经典统计学习理论的形状一致：参数（更准确地说，有效维度）越多、训练数据越少，泛化差距越大。它同时提示：量子模型若真能用少量参数达到经典大模型的效果，其泛化上界也会相应更好——但这只是上界论证，不构成优势的证明。

## 具体应用

### 应用一：量子态分类

把量子态 $|\psi_i\rangle$（或其混态版本）分类到标签 $y_i \in \{0, 1, \ldots, C-1\}$。

QCNN 方法：把 $|\psi_i\rangle$ 直接作为输入态，通过多层量子卷积和池化提取特征，测量最后剩下的量子比特得到分类概率。

优势：处理量子数据时无需先做层析（其代价随 $n$ 指数增长），量子网络直接在量子态上操作；这类任务被认为是最有可能率先体现量子优势的场景，因为输入本身就是量子的。

### 应用二：量子化学性质预测

预测分子的性质（如 HOMO-LUMO 能隙、偶极矩）。

架构：角度编码（分子轨道特征）→ PQC → 测量 → 经典回归层。规模上 $n$ 随轨道数多项式增长、深度 $D = O(n^2)$ 量级、参数 $P = O(D\cdot n)$。与 VQE 的区别在于：这里学习的是"性质对结构的映射"，训练好之后对新分子只需一次前向，而 VQE 对每个分子都要做完整的变分优化。

### 应用三：量子纠错解码

用量子神经网络学习量子纠错码的解码器：输入是综合征（syndrome）测量结果，输出是最可能的错误算符。量子解码器的潜在优势主要在于处理量子数据（如软信息、多次往返的综合征历史）时可以保留相干信息；对纯经典输入的解码问题，它要与成熟的经典解码器竞争。

### 应用四：量子强化学习

用 VQC 作为策略网络，在量子环境中执行强化学习：策略 $\pi(a|s) = |\langle a|U(\vec\theta)|s\rangle|^2$，即对环境态 $|s\rangle$ 施加参数化电路后在动作基上测量。对量子环境（如量子控制问题），策略网络直接在量子态上操作，避免了"测量-经典化-再制备"的中间步骤。

## 当前局限性

### 1. Barren Plateaus

对深层全局电路，梯度随 $n$ 指数衰减，这是量子神经网络训练的核心障碍。缓解手段（局域损失、结构化初始化、浅层化、核方法）各有代价，尚无普适解法。

### 2. 量子优势不确定

目前没有理论证明量子神经网络在标准的经典机器学习任务上相对最佳经典方法具有多项式加速。反方向的结果确实存在：Tang（2019）等人给出了一系列"去量子化"（dequantization）算法，把若干量子推荐/线性代数算法的关键增益用经典的采样技术再现，这提醒我们：在被证明之前，"指数维状态空间"并不自动等于"指数级表达力优势"。

### 3. 噪声敏感性

NISQ 设备的门错误率 $\sim 10^{-3}$ 限制了电路深度，且噪声会让梯度估计的方差增大、有效学习率变小，训练往往在噪声地板上提前停滞。

### 4. 数据加载瓶颈

振幅编码等高效率编码的态制备本身需要 $O(2^n)$ 级别的门数；如果模型的优势依赖大规模经典数据的量子编码，制备开销可能完全抵消量子层的收益。这也是"量子数据进、量子模型学"被认为比"经典数据进"更可行的原因。

### 5. 测量开销

每个期望值的估计需要 $O(1/\varepsilon^2)$ 次测量（推导见 VQE 一篇的测量开销一节）；混合模型每步训练要估计 $O(P)$ 个导数，总测量次数容易超过同规模经典网络的一次前向+反向的成本。

## 当前进展

| 年份 | 贡献 |
|------|------|
| 2018 | Farhi & Neven 提出用于分类的量子神经网络 |
| 2019 | Havlíček et al. 提出量子核方法（Nature） |
| 2020 | TensorFlow Quantum 发布，经典-量子端到端可微训练框架落地 |
| 2021 | Abbas et al. 给出量子神经网络的有效维度与泛化分析 |
| 2022–2023 | 多个实验组在 10+ 量子比特设备上演示混合模型训练；QCNN 的可训练性分析陆续出现 |

## 总结

量子经典混合神经网络把变分量子电路作为可微分层嵌入经典架构：经典层把输入加工成旋转角，量子电路以正弦响应的非线性方式把角度变成期望值读出，经典输出层再做线性组合。前向传播是"经典线性-量子非线性-经典线性"的清晰分工，反向传播由参数平移规则（量子参数）与自动微分（经典参数）在角度接口处拼接而成，整条梯度链都可以解析写出。表达力可以用编码门的生成元谱严格刻画（傅里叶视角），可训练性受 Barren Plateaus 约束，泛化有与经典学习理论同形的上界。它的最终实用价值取决于能否在具体任务——尤其是量子数据任务——上证明相对于经典方法的明确优势。

## 练习题

**练习 1【混合架构的动机与量子的可能性】**（→ [为什么需要量子神经网络](#motivation-quantum-possibilities)）

1. 列出正文给出的经典深度学习面临的三类挑战（维度问题、量子数据、结构先验），各用一句话说明；并写出 $n$ 个量子比特的状态空间维数。
2. 解释：当训练对象本身是量子态时，为什么经典网络必须先做层析、而量子模型可以直接消费量子态？说明对 $n$ 比特系统做完整层析需要确定的实参数个数随 $n$ 如何增长。

> 提示：$n$ 比特的密度矩阵是 $2^n\times 2^n$ 的厄米矩阵，由约 $4^n$ 个实参数确定。

**练习 2【表达力的傅里叶刻画】**（→ [表达力的傅里叶视角](#fourier-expressivity)）

1. 写出 $f(x)=\langle 0|U^\dagger(x)MU(x)|0\rangle$ 的三角多项式展开与频率集合 $\Omega=\{\frac12\sum_k(\lambda'_k-\lambda_k)\}$；对单个编码门 $R_y(x)=e^{-ixY/2}$（$\mathrm{spec}(Y)=\{+1,-1\}$）计算可用频率，并写出 $f(x)$ 的一般形式。
2. 设编码门生成元取两比特算符 $Z_1+Z_2$，由其谱 $\{-2,0,+2\}$ 计算可用频率集合；再解释数据重上传为什么能把频带按重上传次数倍增、代价是什么。

> 提示：频率是 $(\lambda'-\lambda)/2$ 对所有谱值对的取值；重上传让多个门的相位差相加组合出新的频率。

**练习 3【参数化电路与三种数据编码】**（→ [数据编码](#data-encoding)）

1. 写出角度编码与振幅编码把 $\mathbf{x}\in\mathbb{R}^d$ 变成量子态的公式；对 $d=8$ 分别计算两种编码所需的量子比特数，并写出振幅编码制备一般态的门开销量级。
2. 解释 IQP 编码的核 $|\langle\mathbf{x}|\mathbf{x}'\rangle|^2$ 为什么包含特征的所有两两乘积 $x_ix_j$，并说明其两比特因子 $e^{ix_ix_jZ_iZ_j/2}$ 如何用一个受控相位门加可吸收的单比特旋转实现。

> 提示：从二次相位 $e^{ix_ix_j}$ 与恒等式 $|11\rangle\langle 11|=\frac14(I-Z_i)(I-Z_j)$ 出发，仿照正文对 $\mathrm{CP}(\theta)$ 的分解。

**练习 4【经典-量子前馈网络与非线性来源】**（→ [架构一：经典-量子前馈网络](#hybrid-feedforward-network)）

1. 依次写出架构一四层的公式：$\vec\phi=\sigma(W_1\mathbf{x}+\mathbf{b}_1)$、密度矩阵 $\rho$、读出 $z_i=\mathrm{Tr}[\rho Z_i]$ 与输出 $\hat y=\mathbf{w}_2^\top\vec z+b_2$，并指出 $z_i$ 的取值范围。
2. 说明测量映射 $\rho\mapsto\mathrm{Tr}[Z_i\rho]$ 对 $\rho$ 是线性的、而混合模型仍是非线性的：对 $U(\theta)=e^{-i\theta G/2}$、$G^2=I$，推导读出具有 $z(\theta)=a\sin\theta+b\cos\theta+C$ 的正弦响应，并指出非线性来自哪里。

> 提示：展开 $e^{-i\theta G/2}=\cos\frac{\theta}{2}\,I-i\sin\frac{\theta}{2}\,G$，代入 $z=\langle0|U^\dagger M U|0\rangle$ 后用半角恒等式归并。

**练习 5【QCNN 与序列、生成变体】**（→ [架构二：QCNN](#qcnn-variants)）

1. 按顺序描述 QCNN 中一层"量子卷积"与一层"量子池化"分别做什么，并说明"参数跨位置共享"对应经典 CNN 中的什么机制。
2. 解释：为什么 QRNN 的隐藏态不能用强测量直接读出（后果是什么）？为什么 QCNN 的参数量随层数只线性增长、且池化能把有效自由度减半？

> 提示：测量坍缩会破坏叠加态，隐藏态因此退化为经典隐藏态；QCNN 每层只有 $O(1)$ 个共享参数，池化测量掉约一半量子比特。

**练习 6【混合反向传播与参数平移】**（→ [反向传播：从损失到各层参数](#backprop-chain)）

1. 对 $\mathcal{L}=\frac12(\hat y-t)^2$、$\hat y=\sum_i w_{2,i}z_i+b_2$，推导 $\frac{\partial\mathcal{L}}{\partial w_{2,i}}=(\hat y-t)z_i$、$\frac{\partial\mathcal{L}}{\partial b_2}=\hat y-t$，并写出误差信号 $\delta_i$ 与 $\frac{\partial\mathcal{L}}{\partial\theta_k}=\sum_i\delta_i\frac{\partial z_i}{\partial\theta_k}$。
2. 写出 $G_k^2=I$ 时的参数平移公式；再把 $R_z(\theta)=e^{-i\theta Z/2}$（生成元谱 $\{+1,-1\}$、谱隙 $\Delta=2$）代入一般谱隙公式 $\frac{\Delta}{4}[z(\theta+\frac{\pi}{\Delta})-z(\theta-\frac{\pi}{\Delta})]$，验证两种形式一致。

> 提示：$\Delta=\lambda_1-\lambda_2=2$ 时 $\frac{\Delta}{4}=\frac12$、平移角 $\frac{\pi}{\Delta}=\frac{\pi}{2}$，恰好回到第一个公式。

**练习 7【最小示例的端到端训练】**（→ [一个最小可运行示例](#pennylane-example)）

1. 指出示例代码中 `angles = np.tanh(W1 @ x + b1) * np.pi` 对应公式链的哪一段、`quantum_layer` 返回的是什么物理量，并说明 `opt.step(cost, classical_params, theta)` 一次更新中两类参数的梯度分别来自哪里。
2. 修改代码：把编码层改为不做缩放的 `W1 @ x + b1`（角度不再被限制在 $(-\pi,\pi)$），重新训练并解释收敛速度的变化；再解释为什么角度漂出正弦响应的一个周期后，梯度的有效幅度会变小。

> 提示：正弦响应的梯度含 $\cos\theta$ 因子；无界角度使参数更新在不同周期间来回，有效梯度被平均掉。

**练习 8【Barren Plateaus 与泛化上界】**（→ [Barren Plateaus](#barren-plateaus)）

1. 陈述 McClean 等人定理中梯度期望与方差的结论；写出 Abbas 等人的泛化差距上界 $O(\sqrt{d_{\text{eff}}/N_{\text{train}}})$，并指出 $d_{\text{eff}}$ 与 $N_{\text{train}}$ 各由什么决定。
2. 解释局域代价函数与结构化初始化为什么能缓解 Barren Plateaus（各给出一两句机理），并说明核方法为什么完全绕开参数优化、其 $O(N_{\text{train}}^2)$ 代价从何而来。

> 提示：局域观测的方差不随 $n$ 指数衰减；核矩阵的每个元素 $K(\mathbf{x},\mathbf{x}')$ 都要一次电路评估，训练集内共 $N_{\text{train}}^2$ 对。

---

**参考文献：**

1. Farhi, E., & Neven, H. (2018). *Classification with quantum neural networks on near term processors.* arXiv:1802.06002.
2. Havlíček, V., et al. (2019). *Supervised learning with quantum-enhanced feature spaces.* Nature, 567(7747), 209-212.
3. McClean, J. R., et al. (2018). *Barren plateaus in quantum neural network training landscapes.* Nature Communications, 9, 4812.
4. Abbas, A., et al. (2021). *The power of quantum neural networks.* Nature Computational Science, 1, 403-409.
---

> 返回目录：[量子计算算法教程系列](https://chenzhaoyun.com/index.php/archives/54/)