[论文解读] Strong uniqueness principle for fractional polyharmonic operators and applications to inverse problems

Ching-Lung Lin, Hongyu Liu+1 more|arXiv (Cornell University)|Jul 3, 2023

Numerical methods in inverse problemsMathematics被引用 3

一句话总结

本文为分数阶多调和算子——由不同阶数和各向异性系数的多个分数阶拉普拉斯算子组合而成——建立了强有力的唯一性原理，使得在反问题中能够唯一恢复势函数、源项以及非各向同性系数。通过新颖的唯一延拓性质，证明了在有界区域内，Dirichlet-to-Neumann 映射唯一确定了底层的多分数阶算子及其系数。

ABSTRACT

In this work, we are concerned with inverse problems involving poly-fractional operators, where the poly-fractional operator is of the form \[P( (-Δ_g)^s)u := \sum_{i=1}^M α_i(-Δ_{g_i})^{s_i}u\] for $s=(s_1,\dots,s_M)$, $0

研究动机与目标

• 为涉及由不同阶数和各向异性结构组成的多个分数阶拉普拉斯算子的多分数阶算子的反问题建立数学框架。
• 为包含此类多分数阶算子的方程解建立强有力的唯一延拓性质（UCP），这对于证明反问题中的可识别性至关重要。
• 证明从外部测量中唯一识别势函数 $ q(x) $、半线性情况下的源项，以及与非各向同性分数阶算子相关的系数 $ \alpha_i $。
• 通过引入包含局部与非局部行为混合的通用多分数阶算子框架，扩展现有分数阶算子反问题的研究结果。

提出的方法

• 形式化多分数阶算子 $ P((- abla)^s) = \sum_{i=1}^M \alpha_i (-\Delta_{g_i})^{s_i} $，其中 $ 0 < s_1 < \cdots < s_M < \infty $，$ s_M \notin \mathbb{Z} $，且 $ g_i $ 表示不同的各向异性度量。
• 定义 Dirichlet-to-Neumann（DtN）映射 $ \mathcal{M}_{P,q} $，将外部 Dirichlet 数据 $ f \in H^{s_M}(\Omega^c) $ 映射到 $ \Omega $ 上解的非局部算子作用结果，从而实现基于测量的反演。
• 建立一种新颖的强唯一延拓原理（定理 2.11–2.12），断言若解 $ u $ 在子集 $ W \subset \Omega^c $ 上消失，则在 $ \mathbb{R}^n $ 上恒有 $ u \equiv 0 $，前提是 $ \tilde{P}((- abla_{\tilde{g}})^{\tilde{s}})u = 0 $ 在 $ \Omega $ 内成立。
• 使用线性化技术将反问题转化为分析两个不同配置下解的差值 $ \tilde{u} = u_1 - u_2 $，并利用 UCP 推导出系数相等。
• 应用 UCP 证明：若对某个非零外部数据 $ f $，DtN 映射一致，则底层算子 $ (P^1, q_1) $ 与 $ (P^2, q_2) $ 必须在 $ \Omega $ 的某个子集内完全相同。
• 利用谱理论与泛函分析工具，包括强 coercivity 假设和非特征值条件 (1.4)，以确保正向问题解的存在性与唯一性。

实验结果

研究问题

• RQ1 能否为包含不同阶数和各向异性系数的多个分数阶拉普拉斯算子的多分数阶方程解建立强唯一延拓性质？
• RQ2 此类多分数阶算子的 Dirichlet-to-Neumann 映射是否足以唯一识别势函数 $ q(x) $、半线性情况下的源函数，以及非各向同性算子的系数 $ \alpha_i $？
• RQ3 在单个方程中同时包含局部与非局部算子，如何影响唯一延拓性质与反问题的可识别性？
• RQ4 所提出的 UCP 框架能否扩展至恢复分数阶算子中的变系数或变指数？
• RQ5 唯一恢复系数 $ \alpha_i $ 所需的最少测量次数是多少？能否通过单次测量结合所提出的 UCP 实现？

主要发现

• 为多分数阶方程的解建立了一种新的强唯一延拓原理，表明若解在外部区域的非空开集上消失，且在内部满足多分数阶方程，则解在 $ \mathbb{R}^n $ 上恒为零。
• 在算子满足强 coercivity 条件且 0 不是特征值的前提下，Dirichlet-to-Neumann 映射 $ \mathcal{M}_{P,q} $ 唯一确定了 $ \Omega $ 内的势函数 $ q(x) $。
• 在正向问题解在 $ \Omega $ 的某个合适子集上非零的假设下，与非各向同性分数阶拉普拉斯算子相关的系数 $ \alpha_i(x) $ 可通过单次测量唯一恢复。
• 通过线性化方法与强 UCP，可唯一识别半线性情况下的源函数，从而消除解之间差异。
• 该框架允许在多分数阶算子同时包含非局部与局部分量时恢复系数，推广了先前关于高阶分数阶拉普拉斯算子的结果。
• 结果在最小假设下依然稳健：区域有界且为 Lipschitz 域，最高阶指数 $ s_M $ 非整数，解空间为 $ H^{s_M}(\mathbb{R}^n) $。