[论文解读] Continuous regularized Gauss-Newton-type algorithm for nonlinear ill-posed equations with simultaneous updates of inverse derivative

А. Г. Рамм, Alexandra Smirnova|ArXiv.org|Mar 28, 2001

Numerical methods in inverse problems参考文献 9被引用 10

一句话总结

本文提出了一种用于在希尔伯特空间中求解非线性不适定方程的连续正则化高斯-牛顿型算法，通过微分系统同时更新算子 $(F'^*(x(t))F'(x(t)) + \varepsilon(t)I)$ 的逆。该方法在不依赖 $F'(x)$ 的谱假设条件下确保收敛，转而依赖于一种源型条件，并证明了解轨迹的全局存在性与收敛性，满足对所有 $t \geq 0$ 有 $\|x(t) - \hat{x}\| < R\varepsilon(t)$。其主要贡献在于为非线性反问题构建了一个稳健、无谱依赖的正则化框架。

ABSTRACT

A new continuous regularized Gauss-Newton-type method with simultaneous updates of the operator $(F^{\pr*}(x(t))F'(x(t))+\ep(t) I)^{-1}$ for solving nonlinear ill-posed equations in a Hilbert space is proposed. A convergence theorem is proved. An attractive and novel feature of the proposed method is the absence of the assumptions about the location of the spectrum of the operator $F'(x)$. The absence of such assumptions is made possible by a source-type condition.

研究动机与目标

• 开发一种用于求解希尔伯特空间中非线性不适定算子方程的连续正则化方法，其中弗雷chet导数 $F'(x)$ 不具备有界逆。
• 消除对 $F'(x)$ 谱的限制性假设，如非负性或位置约束，这些假设在先前的正则化方法中常见。
• 设计一种新算法，通过耦合微分系统同时更新解 $x(t)$ 和正则化逆算子 $(F'^*F' + \varepsilon(t)I)^{-1}$。
• 在对 $F$ 的结构假设最小的前提下，建立解轨迹的全局存在性与收敛性。
• 为一种连续、稳定且收敛的方法提供理论基础，适用于光滑性有限或导数病态的反问题。

提出的方法

• 该方法被表述为微分方程组：$\dot{x}(t) = -B(t)[F'^*(x(t))F(x(t)) + \varepsilon(t)(x(t) - x_0)]$ 和 $\dot{B}(t) = -[(F'^*F' + \varepsilon(t)I)B(t) - I]$，其中 $B(0) \in L(H)$ 且 $\varepsilon(t) > 0$ 单调递减至零。
• 逆算子 $B(t) := (F'^*F' + \varepsilon(t)I)^{-1}$ 与 $x(t)$ 一同连续更新，避免重复矩阵求逆，提升数值稳定性。
• 收敛性分析依赖于一种源型条件，并利用积分不等式将解轨迹控制在真解 $\hat{x}$ 周围半径为 $R\varepsilon(t)$ 的球内。
• 对系统的积分形式应用逐次逼近法，以证明在 $[0, \infty)$ 上解的存在性与唯一性。
• 证明利用了非线性算子的利普希茨连续性，并在连续函数的巴拿赫空间中使用压缩映射原理，确保迭代格式的一致收敛。
• 该方法通过依赖源条件和 $\varepsilon(t)$ 的衰减，避免了对 $F'(x)$ 谱的假设，如非负性。

实验结果

研究问题

• RQ1 能否设计一种连续正则化高斯-牛顿方法，实现解与法方程算子逆的同步更新，且无需对 $F'(x)$ 的谱施加假设？
• RQ2 当 $F'(x)$ 不具备有界逆时，如何保证非线性不适定方程的收敛性？
• RQ3 在正则化与动态逆更新存在的情况下，何种条件可确保解轨迹的全局存在性与有界性？
• RQ4 缺乏谱假设（如 $F'(x) \geq 0$）对方法的收敛行为与稳定性有何影响？
• RQ5 是否能以正则化参数 $\varepsilon(t)$ 和源条件参数的形式量化收敛速率？

主要发现

• 当 $t \to \infty$ 时，解轨迹 $x(t)$ 收敛于方程 $F(x) = 0$ 的真解 $\hat{x}$，且对所有 $t \geq 0$ 满足 $\|x(t) - \hat{x}\| < R\varepsilon(t)$，其中 $R$ 为正常数。
• 该方法确保了耦合微分系统 (2.3)–(2.4) 在 $[0, \infty)$ 上的全局存在性与唯一性，通过逐次逼近与压缩映射原理建立。
• 收敛性不依赖于假设 $F'(x)$ 的谱位于特定区域（如非负实数），相较于先前工作具有显著改进。
• 证明依赖于源型条件与积分不等式，以控制解与逆算子的增长，确保有界性与收敛性。
• 在 $\varepsilon(t)$ 为正且单调递减至零的假设下，推导出界 $\|x(t) - \hat{x}\| < R\varepsilon(t)$，其中 $\varepsilon(0)$ 与 $\|B(0)\/$ 影响半径 $R$。
• 该方法对 $F'$ 缺乏结构假设（如单调性或可逆性）具有鲁棒性，并通过连续正则化与动态逆更新保持稳定性。