[论文解读] Iterated quasi-reversibility method applied to elliptic and parabolic data completion problems

Jérémi Dardé|arXiv (Cornell University)|Mar 30, 2015

Numerical methods in inverse problems参考文献 42被引用 29

一句话总结

本文提出了一种迭代拟逆法，用于稳定病态的椭圆与抛物型数据补全问题，如泊松方程与热方程在边界数据不完整时的情形。通过使用固定且条件良好的参数ε，迭代求解正则化变分问题，即使在高斯噪声严重的条件下，该方法仍能收敛至真实解，利用Morozov残差原则确定迭代停止时机，数值实验在腐蚀检测与一维热方程问题中均表现出高度鲁棒性。

ABSTRACT

We study the iterated quasi-reversibility method to regularize ill-posed elliptic and parabolic problems: data completion problems for Poisson's and heat equations. We define an abstract setting to treat both equations at once. We demonstrate the convergence of the regularized solution to the exact one, and propose a strategy to deal with noise on the data. We present numerical experiments for both problems: a two-dimensional corrosion detection problem and the one-dimensional heat equation with lateral data. In both cases, the method prove to be efficient even with highly corrupted data.

研究动机与目标

• 为解决泊松方程与热方程数据补全问题中因缺失边界数据而导致的严重病态性，此类问题常导致解不稳定或无解。
• 开发一种即使在测量数据含有高噪声水平时仍具有效性的正则化策略，尤其适用于腐蚀检测等反问题。
• 建立一个统一的抽象框架，适用于椭圆与抛物型反问题，实现理论与数值处理的统一。
• 证明在迭代方法中使用固定且较大的正则化参数ε，可改善矩阵条件数而不损失重建精度，优于标准拟逆法。

提出的方法

• 将数据补全问题在抽象层面上表述为Ax = y，其中A为希尔伯特空间之间连续、单射但非满射的算子，从而统一处理椭圆与抛物型情形。
• 应用标准拟逆法，引入高阶、适定的变分问题，使用小参数ε进行正则化，以稳定病态问题。
• 提出迭代扩展：在每一步中，以前一次迭代的解作为下一次正则化问题中的源项，形成一个收敛于精确解的序列（在精确数据条件下）。
• 采用Morozov残差原则在存在噪声数据时停止迭代，通过将残差与噪声水平匹配，确保解的稳定性和收敛性。
• 采用符合的有限元方法（u使用P2，p使用RT1），并采用张量积P1⊗P1对时空问题进行离散化，通过预计算矩阵分解保持计算效率。
• 将方法应用于两个测试算例：带有Robin传输条件的二维腐蚀检测问题，以及带有侧向数据的一维热方程问题，使用含噪声的Dirichlet与Neumann数据。

实验结果

研究问题

• RQ1 当提供精确数据时，即使使用固定且较大的正则化参数ε，迭代拟逆法是否仍能收敛至病态数据补全问题的真实解？
• RQ2 当输入数据被噪声污染时，该方法在数值上表现如何？是否能在高噪声水平下仍保持重建精度？
• RQ3 在存在噪声数据时，Morozov残差原则是否能有效确定最优停止迭代次数，从而同时保证稳定性和收敛性？
• RQ4 该方法在相同抽象框架下应用于椭圆（泊松）与抛物（热）方程时，是否仍保持鲁棒性与高效性？
• RQ5 该方法是否不仅能重建解本身，还能准确重建如腐蚀检测问题中不可达边界上的Robin系数等导出量？

主要发现

• 在精确数据条件下，即使使用固定且较大的正则化参数ε = 1，迭代拟逆法对椭圆与抛物型数据补全问题均能收敛至精确解，且该参数显著改善了有限元系统的条件数。
• 在噪声数据条件下，方法实现了稳定且精确的重建，即使Dirichlet数据中相对噪声达到5%，误差水平仍可接受，该结果在二维腐蚀检测与一维热方程问题中均得到验证。
• 残差量（衡量PDE与边界条件违反程度）在迭代过程中单调递减，并稳定在与噪声振幅一致的水平，证实了停止准则的有效性。
• 在高噪声条件下，对不可达边界Γc上Robin系数的重建仍具高精度，方法能成功区分腐蚀区域（μ > 0）与健康区域（μ = 0）。
• 在二维腐蚀问题中，解u在整个区域的相对误差在α = 1%与2%时低于10%，在α = 5%时低于15%，表明具有极强的鲁棒性。
• 在一维热方程情形中，解重建的相对误差在α = 2%时低于5%，在α = 5%时低于10%，证实该方法对带有噪声侧向数据的抛物型问题具有有效性。