Recovery of a Measure-valued Source in the Heat Equation from Sparse Boundary Measurements

技术摘要：从稀疏边界测量中恢复热方程的测度值源项

问题表述
本研究探讨了单位圆盘 $\Omega \subset \mathbb{R}^2$ 在时间区间 $(0, T)$ 内，热方程的反问题（ISP）。前向问题定义为如下初边值问题：
$$\begin{cases} \partial_t u - \Delta u = F(t, x) & \text{在 } (0, T) \times \Omega \text{ 内}, \\ u = 0 & \text{在 } (0, T) \times \partial\Omega \text{ 上}, \\ u(0, x) = 0 & \text{在 } \Omega \text{ 上}, \end{cases}$$
其中源项被假定具有分离变量的形式，$F(t, x) = \sigma(t)\mu(x)$。主要目标是从稀疏的边界测量中唯一确定空间分量 $\mu$，其中 $\mu$ 属于有限有符号拉东测度类（$\mu \in (C(\overline{\Omega}))'$）。具体而言，数据由在两个不同的边界点 $x_1, x_2 \in \partial\Omega$ 处观测到的法向通量 $\partial_\nu u(t, x_\ell)$ 在时间区间 $(0, T)$ 内组成。

该论文通过将可恢复源的范围从点源（狄拉克测度）或 $L^2$ 函数扩展到一般的拉东测度，从而与现有文献区别开来。这种扩展的动机在于环境建模中的应用，其中污染源可能支撑在低维子集（例如线源）上，而这些子集无法用标准的 $L^2$ 函数进行充分建模。

方法论与分析框架
理论分析依赖于泛函分析、谱理论和复分析的结合。其方法论通过以下步骤进行：

1. 适定性与正则性： 作者在“转置意义”（transposition sense）下定义解，以处理测度值源项的低正则性。他们建立了 $L^2(Q)$ 中解的唯一存在性，并证明了在远离源支撑集的区域，解具有更高的正则性（在空间靠近边界处为 $C^1$）。
2. 谱分解： 利用单位圆盘上拉普拉斯算子的特征系统 $\{(\lambda_n, \phi_n)\}$，将解 $u$ 分解为傅里叶级数。这导致边界通量 $\partial_\nu u$ 的表示涉及相关椭圆问题（$-\Delta G = \mu$）的解 $G$ 以及包含 $\sigma'$ 的时间相关卷积项。
3. 时间解析性： 一个关键步骤是假设时间系数 $\sigma(t)$ 在子区间 $(T_1, T)$ 上为常数且满足 $\sigma(0)=0$。在这些条件下，作者证明了边界通量在时间上具有解析延拓。这使得利用解析函数的唯一延拓原理，可以将通量的相等性从观测区间 $(0, T)$ 扩展到整个半直线 $(0, +\infty)$。
4. 拉普拉斯变换与唯一性： 通过对扩展后的通量相等式应用拉普拉斯变换，问题被简化为涉及 $\sigma$ 的拉普拉斯变换与 $\mu$ 的谱系数之间的代数恒等式。利用 $\sigma$ 非一致消失这一事实，作者推导出了一个要求差值测度的谱系数必须为零的基本恒等式。
5. 几何条件： 针对一般拉东测度的唯一性证明依赖于单位圆盘的具体几何结构以及特征函数（贝塞尔函数）的显式形式。证明要求观测点 $x_1 = e^{i\theta_1}$ 和 $x_2 = e^{i\theta_2}$ 满足泛型条件 $\theta_1 - \theta_2 \notin \pi\mathbb{Q}$。该条件确保了测度傅里叶系数的线性系统具有平凡解，从而意味着该测度为零。

主要贡献

• 向拉东测度的扩展： 本文提供了第一个关于在稀疏边界测量下，当源 $\mu$ 为一般有限有符号拉东测度（而非局限于点源或 $L^2$ 函数）时，反问题具有正向唯一性结果的研究。
• 放宽的时间依赖性假设： 不同于以往要求 $\sigma$ 为分段常数的工作，本研究允许 $\sigma$ 为一般的 $H^1(0, T)$ 函数，只要它在 $(T_1, T)$ 上为常数且在 $t=0$ 时为零。
• 数值可行性： 作者通过数值实验补充了理论分析，展示了对以下内容的重建：
  • 使用正则化高斯-牛顿法重建点源（狄拉克测度的和）。
  • 使用 Levenberg-Marquardt 方法重建线源（支撑在曲线上的测度）。
    数值结果表明，即使在有噪声的数据（高达 10% 的噪声）下，所提方法也能恢复源的位置和形状，尽管随着未知数数量的增加，问题的病态性也会增加。

结果与主张
主要理论结果（定理 1.1）指出：如果两个测度 $\mu_1$ 和 $\mu_2$ 在两个边界点 $x_1, x_2$ 处产生的通量在 $(0, T)$ 内相同，且时间系数 $\sigma$ 满足所述的正则性和非消失条件，那么 $\mu_1 = \mu_2$，前提是观测点的角间距不是 $\pi$ 的有理倍数。

作者明确指出了目前方法的局限性：

• $\sigma$ 在 $(T_1, T)$ 上为常数的假设是针对此类特定稀疏测量表述的所有已知工作中都存在的技术性要求，目前尚不清楚是否可以放宽。
• 分析依赖于单位圆盘上拉普拉斯特征函数的显式性质，这限制了在无需进一步开发的情况下将其直接应用于其他领域。
• 对于一般的拉东测度，要求 $\sigma(0)=0$ 是必要的，但对于更正则的源（如 $L^2$），这一要求可能是可以移除的。

意义
本文声称其意义在于弥合了理论反问题与涉及奇异源（如线状污染）的实际应用之间的鸿沟。通过建立仅凭两个边界点即可从拉东测度中恢复源的唯一性，这项工作为在无法获得或难以实现密集边界数据的场景下提供了严谨的数学基础。数值部分进一步验证了理论发现，表明尽管测量具有稀疏性，但重建此类奇异源在计算上是可行的。