Improved stability of low Fourier modes in inverse problems for potentials

技术摘要：势函数反问题中低阶傅里叶模态的稳定性改进

问题陈述
本研究探讨了椭圆方程反问题的稳定性，特别关注如何从狄利克雷-到-纽曼（DN）映射 $\Lambda_q$ 中恢复未知的势函数 $q$。此类反问题的核心挑战在于其固有的病态性，通常导致对测量噪声高度敏感的对数稳定性估计。作者研究了势函数的特定结构特征——即其低频傅里叶模态——是否可以在对两个势函数之差 $Q = q - q'$ 施加特定正则性假设的情况下，实现改进的稳定性（Lipschitz、亚 Hölder 或 Hölder 稳定性）。

研究考虑了关于 $Q$ 的三种不同情景：

1. 带限差异（Bandlimited differences）： $Q$ 是 $M$-带限的（当 $|k| > M$ 时，傅里叶模态消失）。
2. 实解析差异（Real-analytic differences）： $Q$ 是实解析且周期的。
3. 超指数衰减差异（Super-exponentially decaying differences）： $Q$ 的傅里叶系数以 $e^{-c|k|^{n/2}}$ 的速率衰减。

控制方程为亥姆霍兹方程 $(-\Delta - \lambda^2 + q)u = 0$（对于大频率 $\lambda$）和薛定谔方程 $(-\Delta + q)u = 0$（对于 $\lambda=0$）。

方法论
证明过程依赖于 Alessandrini 恒等式，该恒等式将 DN 映射之差与涉及势函数差和前向问题解的积分联系起来：
$$\langle (\Lambda_q - \Lambda_{q'})u, u' \rangle_{\partial \Omega} = \int_\Omega (q - q') u \bar{u}' \, dx$$
为了提取关于 $Q$ 的傅里叶系数的信息，作者构造了特定的测试函数（解 $u, u'$），这些函数在高频下振荡以隔离特定的傅里叶模态。

• 实几何光学解（RGO）： 对于带限情况（定理 1.1），作者采用了形式为 $u = e^{i\eta \cdot x}(1+r)$ 的实几何光学解，其中 $|\eta| = \lambda$。这些解被用于探测高频下的势函数差异。分析过程涉及估计余项 $r$，并利用类 Landau–Pollak–Slepian 型算子的正定性来反转 DN 映射差与傅里叶系数之间的关系。
• 复几何光学解（CGO）： 对于实解析和超指数衰减情况（定理 1.2 和 1.3），作者使用了具有复相位 $\rho, \rho' \in \mathbb{C}^n$ 的 CGO 解。这些解允许在周期性设定 $\Omega = (-\pi, \pi)^n$ 中隔离傅里叶模态。
• 傅里叶衰减估计： 方法论的一个关键组成部分是界定“尾部”（模态 $|k| > M$）的傅里叶级数。
  • 对于实解析函数，作者利用全纯延拓的柯西估计来导出傅里叶系数的指数衰减。
  • 对于超指数衰减，他们假设 $|Q_k| \leq C e^{-c|k|^{n/2}}$ 这一拉伸指数衰减界限。
  • 引理 3.2 提供了在高频壳层上对衰减指数尾部求和的估计，这对于控制稳定性估计中的误差项至关重要。

主要贡献与结果

1. 带限差异的 Lipschitz 稳定性（定理 1.1）：
   作者证明，如果差异 $Q = q - q'$ 是 $M$-带限的，那么只要频率 $\lambda$ 相对于带宽足够大（$\lambda \gtrsim M^{n/2}$），其所有傅里叶模态都可以以 Lipschitz 稳定性 被恢复。

   • 意义： 不同于典型的有限维 Lipschitz 结果（其常数随维度呈指数增长），这里的 Lipschitz 常数是关于 $M$ 一致的（仅取决于 $\Omega, n, E$）。
   • 假设： 仅需要 差异 $Q$ 是带限的；单个势函数 $q, q'$ 仅需在 $L^\infty$ 中有界。

2. 实解析差异的亚 Hölder 稳定性（定理 1.2）：
   对于 $Q$ 为实解析且周期的情形，作者确立了 亚 Hölder 稳定性。其稳定性模态表现为 $\exp(-c(\log(1/\delta))^{2/n})$，优于对数稳定性，但差于任何 Hölder 模态 $\delta^\alpha$。

   • 分辨率： 可恢复的傅里叶模态数量 $M(\delta)$ 随着测量误差 $\delta$ 的减小而增加。具体而言，$M(\delta) \sim (\log(1/\delta))^{2/n}$。
   • 假设： $Q$ 可全纯延拓至半径为 $\epsilon$ 的复邻域。

3. 超指数衰减的 Hölder 稳定性（定理 1.3）：
   如果 $Q$ 的傅里叶系数以超指数速率 $e^{-c|k|^{n/2}}$ 衰减，则稳定性提升至 Hölder 稳定性 ($\delta^\alpha$)。

   • 指数： Hölder 指数为 $\alpha = c/(\gamma + c)$，其中 $c$ 是衰减率，$\gamma$ 是几何常数。当 $c \to \infty$ 时，$\alpha \to 1$，趋于 Lipschitz 稳定性。
   • 分辨率： 与解析情况类似，随着 $\delta \to 0$，可恢复的模态数量 $M(\delta)$ 趋于无穷大，即 $M(\delta) \sim (\log(1/\delta))^{2/n}$。

意义与主张
本文声称推进了对非线性反问题中“增加的稳定性”和“增加的分辨率”的理解。

• 放宽假设： 一个主要的贡献是，正则性假设（带限、解析或超指数衰减）仅适用于 差异 $q - q'$，而非单个势函数。势函数本身仅要求是有界的（$L^\infty$）。这避免了未知量必须属于已知有限维子空间的严格假设。
• 一致常数： 在带限情况下，稳定性常数不会随未知量的维度（即模态数 $M$）呈指数级增长，这与标准的有限维反问题结果形成对比。
• 定量分辨率： 结果提供了测量精度（$\delta$）与分辨率（可恢复的傅里叶模态数量）之间的定量关系。随着测量变得更加精确，更多的低频模态可以被稳定地恢复，根据差异的正则性，其在对数、Hölder 或 Lipschitz 稳定性区间之间进行转换。
• 非线性： 本工作直接处理非线性反问题，超越了线性化设置（例如线性散射中的奇异值估计），在这些设置中此前已观察到类似的增加稳定性现象。

作者指出，估计中的常数在恢复模态的数量上是一致的，且结果适用于一般有界区域（定理 2.1）或特别是在傅里叶分析最为自然的环面/立方体设定中。本文并未提出实验应用，而是侧重于数学反问题理论稳定性估计的研究。