[论文解读] Stability result for elliptic inverse periodic coefficient problem by partial Dirichlet-to-Neumann map

Mourad Choulli, Yavar Kian+1 more|arXiv (Cornell University)|Jan 20, 2016

Numerical methods in inverse problems参考文献 26被引用 7

一句话总结

本论文为在无限圆柱形波导中，从略大于一半边界上测量的局部 Dirichlet-to-Neumann 数据，建立了恢复 1-周期势函数的对数-对数稳定性估计。通过傅里叶分析与复几何光学解法，作者证明了即使在法向量数据被限制在边界子集上时，势函数在 $ H^{-1} $-范数下仍能以对数-对数依赖于测量误差的方式稳定确定。

ABSTRACT

We study the inverse problem of identifying a periodic potential perturbation of the Dirichlet Laplacian acting in an infinite cylindrical domain, whose cross section is assumed to be bounded. We prove log-log stable determination of the potential with respect to the partial Dirichlet-to-Neumann map, where the Neumann data is taken on slightly more than half of the boundary of the domain.

研究动机与目标

• 研究在具有有界截面的波导中，从部分边界测量数据恢复 1-周期势函数的逆问题。
• 在无界区域（具体为无限圆柱形波导）中，建立逆系数问题的稳定性估计。
• 将有界区域中的部分数据稳定性结果推广至无界区域，特别是针对周期系数的情形。
• 证明当法向量数据在略大于总边界一半的子集上测量时，势函数可被稳定恢复。
• 推导出对数-对数型稳定性估计，量化重建误差对测量误差的依赖关系。

提出的方法

• 在周期方向 $ x_1 $ 上进行傅里叶分解，将三维问题转化为由傅里叶模 $ k \in \mathbb{Z} $ 参数化的二维问题族。
• 构造适用于周期结构与波导几何形状的薛定谔方程 $ (-\Delta + V)v = 0 $ 的复几何光学（CGO）解。
• 使用在子集 $ G \subset \partial\Omega $ 上定义的局部 Dirichlet-to-Neumann 映射 $ \Lambda_V $，其中 $ G $ 是边界上 $ \xi $-阴影面的邻域。
• 利用 Carleman 估计与能量估计，推导出对应于两个势函数 $ V_1 $ 与 $ V_2 $ 的解之差的 $ L^2 $-估计。
• 应用对偶论证，将势函数差的 $ H^{-1} $-范数与测量集上解差的 $ L^2 $-范数联系起来。
• 优化 CGO 假设中的参数（频率 $ \rho $、缩放参数 $ r $），以平衡指数增长与衰减项，从而导出对数-对数型稳定性估计。

实验结果

研究问题

• RQ1 在无限波导中，若在略大于一半边界的局部区域测量 Dirichlet-to-Neumann 数据，是否能稳定恢复周期势函数 $ V \in L^\infty(\Omega) $？
• RQ2 当法向量数据被限制在边界子集上时，识别周期势函数的逆问题的最优稳定性估计是什么？
• RQ3 当边界数据不完整且区域无界时，稳定性如何依赖于测量误差？
• RQ4 是否可在无限圆柱域中为周期系数建立对数-对数稳定性估计，从而将先前在有界区域的结果推广？
• RQ5 照明区与阴影区的几何结构在逆问题稳定性中起何种作用？

主要发现

• 论文建立了对数-对数型稳定性估计：$ \|V\|_{H^{-1}(\check{\Omega})} \leq C (\ln |\ln \gamma|)^{-1} $，其中 $ \gamma = \|\Lambda_{V_2,\theta} - \Lambda_{V_1,\theta}\| $，适用于小测量误差。
• 该稳定性估计在法向量数据测量区域 $ G $ 为 $ \xi $-阴影面邻域（略大于一半边界）的假设下成立。
• 结果对区域无界性具有鲁棒性，因为波导结构与周期性允许对无穷远处解的有效控制。
• 关键技术突破在于构造出在测量集上具有受控 $ L^2 $-范数的 CGO 解，并在指数权重中实现参数的最优选择。
• 最终的稳定性估计通过平衡 Carleman 估计中的指数项并优化频率参数 $ \rho $ 而导出，导致对测量误差呈现双重对数依赖。
• 稳定性估计中的常数 $ C $ 仅依赖于截面 $ \omega $ 的几何形状、势函数的 $ L^\infty $-范数上界 $ M_+ $ 以及测量区域 $ G' $，而不依赖于势函数本身。