第 22 课 · 量子 ODE/PDE 求解器——从哈密顿量模拟到 Schrödingerization · 模块四:应用与前沿专题

PDE II:Poisson 方程与条件数

  • 命题框:中心差分 $A = (1/h^2)\cdot\mathrm{tridiag}(-1, 2, -1)$(精度 $O(h^2)$);特征值 $\lambda_k = (4/h^2)\cdot\sin^2(k\pi/(2(N+1)))$,$k = 1, \ldots, N$
  • 推导/演算步骤:
    1. Taylor 展开相加消奇次项:$u''(x_i) = (u_{i+1} - 2u_i + u_{i-1})/h^2 - (h^2/12)\cdot u'''' + O(h^4)$
    2. 特征向量 $v_i = \sin(i\theta)$:和差恒等式给 $\lambda(\theta) = (2-2\cos\theta)/h^2 = 4\sin^2(\theta/2)/h^2$;边界条件 $\sin((N+1)\theta) = 0$ → $\theta = k\pi/(N+1)$
    3. 两端估计:$\lambda_{\min} \approx (4/h^2)\cdot(\pi/(2(N+1)))^2 = \pi^2$(与连续算符 $-d^2/dx^2$ 的第一特征值自洽);$\lambda_{\max} \approx 4/h^2 = 4(N+1)^2$ → $\kappa(A) = 4(N+1)^2/\pi^2 = O(N^2)$
    4. 代入演算(PDE 教材练习 7 第 2 题):$N = 9$、$h = 0.1$ → $\lambda_1 = 400\cdot\sin^2(\pi/20) = 400\times(0.156)^2 \approx 9.79$ vs $\pi^2 \approx 9.87$(偏差 $\approx 0.8\%$);$\lambda_9 = 400\cdot\cos^2(\pi/20) \approx 390$ → $\kappa \approx 390/9.79 \approx 40$
    5. 量子 vs 经典:$O(\kappa\cdot\mathrm{polylog}(N^d)\cdot\log(1/\varepsilon))$ vs $O(N^d\cdot\sqrt{\kappa}\cdot\log(1/\varepsilon))$ → 优势要求 $N^d \gg N^2\cdot\mathrm{polylog}$,即 $d \ge 3$;$d = 1, 2$ 时条件数抵消对数编码收益
  • 结果框:$\kappa = O(N^2)$ 是命门;缓解 = 预条件化 / 谱方法(下页);输出必须是泛函而非全场
类型definition-derivation