[论文解读] A backward problem for the time-fractional pseudo-parabolic equation with a variable coefficient

Arshyn Altybay|arXiv (Cornell University)|Mar 15, 2026

Numerical methods in inverse problems被引用 0

一句话总结

该论文证明了在带时变系数的时分数伪抛物方程中重建初始状态的存在性与唯一性，提出了稳定的有限差分直接求解器，并利用 Tikhonov 正则化从最终时间观测中识别初始数据，数值验证包括含噪声情形。

ABSTRACT

This work addresses an inverse reconstruction task for a time-fractional pseudo-parabolic model with a temporally varying coefficient. By imposing Dirichlet boundary conditions, we aim to recover the unknown initial state from observations collected at the final time. From a theoretical perspective, we derive existence and uniqueness results by proving that, under suitable hypotheses, the problem admits a unique solution. Computationally, we introduce a finite-difference discretisation based on a time-stepping strategy and provide a detailed stability and convergence analysis. Leveraging the resulting forward solver, we then formulate an initial-data identification procedure using Tikhonov regularisation. The proposed approach is validated with numerical simulations, and its resilience is assessed via experiments that incorporate perturbations in the final-time measurements.

研究动机与目标

• 为带变系数的时分数伪抛物方程的向后时间重构问题提供动机与形式化描述。
• 证明经典解的存在性与唯一性，并给出解与初始数据的显式重构公式。
• 为直接问题开发并分析稳定的有限差分方案，建立稳定性与收敛性结果。
• 提出基于 Tikhonov 正则化的从最终时间测量中识别未知初始态的重构程序，并进行数值验证。
• 通过数值实验评估对含噪最终时间数据的鲁棒性。

提出的方法

• 将解分解为正弦本征函数以获得模态方程。
• 使用 Caputo 导数将每一模态转化为一个Volterra积分方程。
• 给出以模态系数（A_k, B_k）和最终数据（ψ_k）表示的显式表示式 u(x,t) 与 u0(x)。
• 在分级时间网格上采用 L1 型离散方式、在空间使用中心差分进行离散化来离散化直接模型。
• 利用分级网格 L1 运算符的强化性来证明全离散方案的稳定性和收敛性。
• 基于正向求解器提出初始数据的 Tikhonov 正则化重构。

Figure 1. Reconstruction results for different fractional orders. Left: exact initial state $u_{0}(x)$ and reconstructed state $\widehat{u_{0}}(x)$ . Right: exact final-time measurement $\psi(x)$ and the final-time state $u(x,T;\widehat{u_{0}})$ generated by evolving the reconstructed initial state

实验结果

研究问题

• RQ1 在所给假设下，带时变系数的时分数伪抛物方程的向后时间问题是否存在唯一的经典解？
• RQ2 是否能够从最终时间数据导出解与初始态的显式表示，以及对数据的稳定性性质如何？
• RQ3 是否存在稳定且收敛的直接问题数值方案以实现初始数据的可靠重构？
• RQ4 是否能够利用所提正向求解器在带噪的最终时间测量下通过 Tikhonov 正则化鲁棒地恢复初始状态？

主要发现

• 在给定假设下向后问题存在唯一的经典解，并给出按模态分量的显式重构公式。
• 初始态系数通过最终时间数据的 A_k(T) > 0 与 B_k(T) 获得，从而实现稳定重构。
• 存在一个稳定性估计，将重构得到的初始态在最终数据和源项的量级下进行界定，显示数据的连续依赖性。
• 为直接问题开发了稳定的有限差分离散化，并证明其稳定性与收敛性。
• 数值实验验证了重构方法并展示了对最终时间测量扰动的鲁棒性。

Figure 2. Exact solution (left), reconstructed solution (centre), and absolute error (right) for $u(x,t)$ for several values of $\alpha$ .