Conditional stability for an inverse problem of a fully-discrete stochastic hyperbolic equation

技术摘要：全离散随机双曲方程逆问题的条件稳定性

问题陈述
本文研究了一个与一维随机双曲方程的全离散近似相关的离散逆问题。定义在时空域 $Q = (0,1) \times (0,T)$ 上的连续模型描述了一个受随机力作用的系统。主要目标是同时确定三个未知量：初始位移（$y_0$）、初始速度（$y_1$）以及随机源项（$g$）。

该逆问题利用两类观测数据构建：

1. 左边界（$x=0$）处的离散空间导数数据。
2. 最终时刻（$t=T$）的解的测量值。

作者聚焦于全离散设定，采用显式有限差分格式对随机偏微分方程（SPDE）进行空间和时间的离散化。这与主要处理连续设定或半离散近似的前人工作形成对比。

方法论
核心方法论依赖于推导并应用专为全离散随机双曲方程设计的Carleman 估计。该方法涉及若干技术步骤：

1. 离散微积分框架：作者建立了一个严格的离散微积分框架，定义了离散函数空间、平移算子、平均算子与差分算子，以及离散分部积分公式。特别关注了处理离散导数所需的双重网格。
2. 权函数构造：利用大参数 $s$ 和 $\lambda$ 构造了一个特定的权函数 $r = e^{s\phi(x,t)}$。该函数旨在处理方程的双曲性质及随机项。
3. Carleman 估计的推导：
   • 引入变量代换 $y = \rho Y$（其中 $\rho = r^{-1}$）以变换算子。
   • 证明过程涉及将变换后的算子分解为项 $A(Y)$、$B(Y,Z)$ 以及余项 $R(Y,Z)$。
   • 关键的是，作者利用了解与布朗运动未来增量（$B^{n+1} - B^n$）之间的独立性，证明了随机积分项的期望为零（$E[\dots] = 0$）。这使得消除离散 Itô 积分项成为可能，这是建立估计的必要步骤。
   • 证明过程仔细处理了由离散化引起的误差项，特别是那些涉及网格尺寸 $\Delta x$ 和 $\Delta t$ 的项。
4. 稳定性分析：将推导出的 Carleman 估计应用于对应不同未知量的两个解之差。通过结合离散能量估计，作者得出了该逆问题的条件稳定性结果。

主要贡献与结果

• 新的离散 Carleman 估计：本文建立了针对全离散随机双曲方程的新颖 Carleman 估计。与连续情形不同，该估计在右侧包含依赖于网格尺寸（$\Delta x$ 和 $\Delta t$）的额外项。具体而言，它包含一个涉及空间导数的离散时间导数（$|D_t D_x y|^2$）的项，以及一个在 $t=T$ 处局部化且涉及源项 $g$ 的边界项。
• 条件稳定性：基于 Carleman 估计，作者证明了一个条件稳定性结果。他们表明，如果未知量属于具有先验界限的容许集，则未知量（$g, y_0, y_1$）之差可由观测数据（边界通量和最终时刻解）之差界定。
• 网格依赖性：稳定性估计明确考虑了离散化参数。与最终时刻源项 $g$ 相关的误差项按 $(\Delta t)^{1/2}$ 缩放。作者指出，随着 $\Delta x, \Delta t \to 0$，这些离散项消失，从而恢复了连续稳定性结果。
• 放宽的最终时刻假设：相较于之前的连续研究（例如 [28]），在最终时刻条件的处理上取得了显著的技术改进。虽然连续唯一性通常假设 $y^{(1)}(T) = y^{(2)}(T)$，但本文推导出的稳定性估计涉及项 $\|y^{(1)}(T) - y^{(2)}(T)\|^{1/2}$。这种 Hölder 型依赖关系更适合数值分析，因为在数值分析中，最终时刻的精确相等很少实现，通常仅存在微小误差。

意义与主张
本文声称填补了关于全离散随机逆问题的文献空白。虽然针对确定性离散逆问题和连续随机逆问题已有大量工作，但本文是首个处理全离散随机双曲情形的工作。

作者强调，离散设定引入了特定的挑战，例如需要处理由随机项离散化引起的额外误差项，以及离散 Itô 公式的特定结构。推导出的稳定性结果被呈现为对数值实现中的误差分析特别有价值，因为它量化了离散化误差如何影响未知参数的重构。

文章最后指出，虽然确定性情形预期具有 Lipschitz 稳定性，但随机性质和离散设定目前仅产生一个关于最终时刻误差具有 $1/2$ 幂次的条件稳定性结果。正如作者所建议，未来的工作将涉及研究更强的稳定性性质以及这些逆问题的数值实现。