[论文解读] Increasing stability in an inverse problem for the acoustic equation

Sei Nagayasu, Günther Uhlmann+1 more|arXiv (Cornell University)|Oct 24, 2011

Numerical methods in inverse problems参考文献 9被引用 4

一句话总结

本文通过证明稳定性估计随频率 $k$ 增加而从对数型过渡到Lipschitz型，建立了声学方程逆问题的递增稳定性。利用复几何光学（CGO）解和精细的算子范数估计，作者推导出一个界，其中折射率 $q$ 的误差以 $k^{-2}\exp(Ck^2)\|\Lambda_1 - \Lambda_2\|_*$（Lipschitz部分）和 $\left(k^2 + \log \frac{1}{\|\Lambda_1 - \Lambda_2\|_*}\right)^{-(2s-n)}$（对数部分）的形式衰减，表明在高频下稳定性得到改善。

ABSTRACT

In this work we study the inverse boundary value problem of determining the refractive index in the acoustic equation. It is known that this inverse problem is ill-posed. Nonetheless, we show that the ill-posedness decreases when we increase the frequency and the stability estimate changes from logarithmic type for low frequencies to a Lipschitz estimate for large frequencies.

研究动机与目标

• 分析声学方程关于折射率 $q(x)$ 的反边界值问题的稳定性。
• 研究当频率 $k$ 增加时，反问题的病态性如何变化。
• 推导一个随 $k$ 增加而从对数型过渡到Lipschitz型行为的稳定性估计。
• 建立 $\widetilde{q} = q_1 - q_2$ 差异的定量界，以 $\|\Lambda_1 - \Lambda_2\|_*$ 表示。

提出的方法

• 作者在有界区域 $\Omega \subset \mathbb{R}^n$，$n \geq 3$ 中，使用声学方程 $(\Delta + k^2 q(x))u = 0$ 的复几何光学（CGO）解。
• 应用Alessandrini的稳定性估计方法，追踪常数对频率 $k$ 的显式依赖关系。
• 该方法涉及构造具有显式 $k$-依赖性的CGO解，并通过加权 $L^2$-型估计来估计其 $H^s$-范数。
• 关键步骤是在CGO构造中引入参数 $T$，以平衡最终估计中Lipschitz项与对数项的权衡。
• 证明分为两种情况，基于 $k^2$ 相对于 $\log(1/\|\Lambda_1 - \Lambda_2\|_*)$ 的大小，分别推导出各自的界。
• 最终估计通过优化 $T$ 并选择参数以最小化误差界获得，从而得到一个结合Lipschitz与对数衰减的混合估计。

实验结果

研究问题

• RQ1 声学方程反问题的稳定性是否随频率 $k$ 增加而改善？
• RQ2 能否推导出一个单一的稳定性估计，以捕捉高频区域中从对数型到Lipschitz型行为的过渡？
• RQ3 在存在俘获射线或一般 $q(x)$ 的情况下，$k$ 的依赖性在稳定性估计的常数中如何体现？
• RQ4 复几何光学（CGO）解在实现高频下更高稳定性中起到什么作用？

主要发现

• 反问题的稳定性估计随 $k$ 增加而从对数型过渡到Lipschitz型，其中对数部分衰减，Lipschitz部分在高频下占主导地位。
• $\|\widetilde{q}\|_{H^{-s}(\mathbb{R}^n)}$ 的界为 $\frac{C}{k^2}\exp(Ck^2)\|\Lambda_1 - \Lambda_2\|_* + C\left(k^2 + \log\frac{1}{\|\Lambda_1 - \Lambda_2\|_*}\right)^{-(2s-n)}$，其中 $C$ 仅依赖于 $n$、$s$、$\Omega$、$M$ 和 $\text{supp}(q_1 - q_2)$。
• Lipschitz项 $\frac{C}{k^2}\exp(Ck^2)\|\Lambda_1 - \Lambda_2\|_*$ 源于CGO解的使用以及对 $q(x)$ 无几何限制，尽管在更简单条件下该指数因子预期可被移除。
• 该估计对所有 $k^2 \geq 1/(C_1 M)$ 一致成立，无需对不同频率范围分别分析，这与以往结果不同。
• 对数部分 $\left(k^2 + \log\frac{1}{\|\Lambda_1 - \Lambda_2\|_*}\right)^{-(2s-n)}$ 随 $k$ 增加而衰减，表明病态性降低。
• 该结果确认了数值观测到的反散射中稳定性递增的现象，并将严格的稳定性分析扩展至具有折射率的声波方程。