第 20 课 · 模块整合实操——本地模拟与含噪模拟 · 模块三:现代算法体系

模拟后端:态矢量与密度矩阵

  • 命题框:混合态 $\rho = \sum_i p_i|\psi_i\rangle\langle\psi_i|$;纯态是特例 $\rho = |\psi\rangle\langle\psi|$
  • 推导步骤(为何态矢量装不下噪声):
    1. 比特翻转通道作用于 $|0\rangle$:$\rho = (1-p)|0\rangle\langle 0| + p|1\rangle\langle 1| = \mathrm{diag}(1-p, p)$
    2. 反证:若 $\rho = |\phi\rangle\langle\phi| = (a,b)^T(a^*,b^*)$,两个对角元 $|a|^2, |b|^2$ 均非零 $\Rightarrow$ 非对角元 $ab^* \ne 0$;但这里 $\rho_{12} = 0$ $\Rightarrow$ 矛盾 $\Rightarrow$ 非纯态,只能用 $\rho$
    3. 纯度检验:$\mathrm{Tr}(\rho^2) = (1-p)^2 + p^2 = 1 - 2p(1-p)$;代入 $p = 0.02$ → $\mathrm{Tr}(\rho^2) = 0.9608 < 1$(纯态恒为 1)
  • 开销演算:态矢量 $2^n$ 个幅值 vs 密度矩阵 $4^n$ 个矩阵元;$n=3$:8 vs 64;每加 1 比特 $\times 2$ vs $\times 4$;$n=20$:$\approx 1.0\times 10^6$ vs $\approx 1.1\times 10^{12}$
  • 验证行:同一 Bell 线路双后端 pmeasure 完全一致:['0.5000','0.0000','0.0000','0.5000'] ✓
  • 报错卡片:NoisySimulator 配 statevector 直接抛 ValueError(原文:simulate_pmeasure is only available for density_operator type OpcodeSimulator backend)
  • 结果框:噪声通道 $\Rightarrow$ 混合态 $\Rightarrow$ 必须 density_matrix;后端表压缩为一行(statevector:无噪/$O(2^n)$/默认;density_matrix:含噪/$O(4^n)$/NoisySimulator 标配;qutip 变体:复杂噪声)
类型definition-derivation