1. 引言
许多反问题都可以被建模成不适定方程
(1.1)
其中,
是一个连续且Fréchet可微的非线性算子,
是Hilbert空间。在实际中,我们一般不会获得精确数据y,而是需要处理满足下式的噪音数据
:
(1.2)
其中
为噪音水平。因为方程(1.1)是不适定的,即解不连续依赖于数据,直接求解往往导致近似解严重偏离真实解,所以需要正则化方法来获得解的合理近似。而Landweber迭代法 [1] 就是一种有效的正则化方法,其经典格式为
(1.3)
其中,
为迭代变量,
是初始迭代点,包含着对方程(1.1)精确解
的先验信息。到目前为止,Landweber迭代已被广泛研究,例如 [2] - [8] 。但是这些文献所考虑的向前算子F大部分都是精确的,而实际中算子一般是带有扰动的,从而方程(1.1)就变成了
(1.4)
其中
为F的扰动算子。因此,本文针对不适定的非线性方程(1.4),在经典Landweber迭代法的基础上构建如下Landweber迭代
(1.5)
其中
是一个被合适选择的固定步长。
需要指出的是,对方程(1.4)的研究已出现一些文献:1994年,陈宏和侯宗义 [9] 提出了方程(1.4)在线性情况下的Tikhonov正则化方法;2000年,Q。Jin [10] 提出了非线性方程(1.4)在离散后求解的迭代正则化高斯–牛顿方法,该问题中扰动算子是通过有限维的子空间去逼近函数空间形成的;2002年,韩波,刘家琦等 [11] 提出了求解非线性方程(1.4)的Landweber迭代法;文献 [12] [13] [14] 给出了Banach空间中当算子带统计误差数据时变分正则化方法的收敛阶;文献 [15] 推导了在Hilbert空间和Banach空间中两种算子噪音模型的收敛速度;文献 [16] 给出了当扰动算子具有一般数据拟合项时变分正则化的误差界;文献 [17] 提供了推导带扰动算子变分正则化误差界的一般框架。
由于反问题的内在不适定性,Landweber迭代序列呈现出半收敛现象。因此,大部分文献(比如 [2] - [8] [10] [11] )采用最经典的偏差原则作为Landweber迭代的停止准则。但是偏差原则严重依赖于噪音水平
,而实际应用中噪音水平往往是不可获得或不可靠的,所以不依赖噪音水平的启发式停止准则被提出:1996年,Martin Hanke与Toomas Raus [18] 首次对多种正则化方法提出了一种选择适当正则化参数的启发式方法,即Hanke-Raus准则;2011年,Bangti Jin与Dirk A. Lorenz [19] 导出并推广了二次变分正则化的两个启发式参数选取准则,即Hanke-Raus准则和拟最优性准则;2016年,Qinian Jin [20] 将关于二次变分正则化的Hanke-Raus启发式准则推广到求解Banach空间不适定反问题的一般变分正则化;2018年,Zhengqiang Zhang与Qinian Jin [21] 针对Banach空间中非平稳迭代Tikhonov正则化方法提出了Hanke-Raus准则;之后,Qinian Jin与Wei Wang [22] ,Zhenwu Fu等人 [23] 分别针对Hilbert空间和Banach空间中的高斯–牛顿方法提出了启发式参数选取准则;2022年,Simon Hubmer与Ekaterina Sherina等人 [24] 对Landweber迭代法的若干启发式停止准则进行了比较。本文为Landweber迭代(1.5)提出以下Hanke-Raus准则:
(1.6)
其中,
为一固定的常数,
满足
,有
。
下面将分析Landweber迭代(1.5)在Hanke-Raus准则(1.6)下的收敛阶。
2. 主要假设条件
本节介绍对Landweber迭代(1.5)进行收敛阶证明的一些假设条件,部分细节可参考 [4] [25] 。
假设2-1 假设扰动算子及其导算子在方程(1.1)的精确解
处的误差估计为:
(2.1)
在下面的理论分析中,令
,
,注意此噪声水平
不用于实际计算。为使Landweber迭代法(1.5)在Hanke-Raus准则下收敛,我们还需要对
附加假设条件。
假设2-2 存在一个常数
,使得对于任意的
,有
其中,
。
下面给出对迭代正则化方法进行收敛阶证明的一些常用的假设条件。
假设2-3 存在
,使得方程(1.1)在
有一个解
,其中
。
切锥条件:存在
,使得
(2.2)
扰动算子的导算子的一致有界性:
(2.3)
存在一类有界线性算子
使得对于某个正常数
,有
(2.4)
源条件:存在
,使得
(2.5)
3. 收敛阶分析
本节将会给出迭代法(1.5)收敛阶的分析过程。为了简化符号,我们令
,
,
。
引理3-1 令假设2-3中条件1),2),3)成立。
是由偏差原则
所定义的,其中
。那么
是个有限数且对于任意的
,
。进一步若源条件(2.5)成立,且
充分小,那么
其中,
是仅依赖于
的某个正常数。
证明:采用与文献 [1] 收敛阶证明过程类似的方法,即可证得结论。
引理3-2 令假设2-3中条件1),2),3)成立。
是由Hanke-Raus准则(1.6)所确定的。进一步若源条件(2.5)成立,且
充分小,那么
(3.1)
其中
仅依赖于
。
证明:可参考文献 [25] 的引理4.3.1的证明过程,其中需要用到本文中的引理3-1。
引理3-3 ( [4] 引理2.9)假设
是正的,存在一个与k无关的正常数
使得
其中
是一个仅依赖于
的正常数。
引理3-4 ( [4] 引理2.10)令
,
,使得
,且设
,
,那么有下面的不等式成立:
引理3-5 令引理3-2的所有条件都成立。若源条件(2.5)成立,且
充分小,那么对于任意的
,有
其中,
是一个仅依赖于
。
证明:首先我们需要找到对于任意的
,
和
的表达式,根据迭代法(1.5),可得
利用(2.4),可得
(3.2)
令
(3.3)
其中
。结合(3.2)和(3.3),可得
由
,通过归纳假设,我们可以得到对于
,有
与
接下来,我们将证明对于
,有
(3.4)
对于
,结合(2.1),(2.3)与(2.5),有
,所以(3.4)成立是显然的。我们采用数学归纳法,假设对于
,(3.4)成立,只需要证明
时结论成立即可。应用(2.4)我们可以得到
(3.5)
通过三角不等式,(2.2),(3.1)与(3.4),可得
,有
将上式结合(2.4),(3.5)与(3.4),则
我们有
其中
。通过上式以及引理3-3,令
,可得
因为
充分小,所以存在一个常数
,使得
,有
(3.6)
基于以上的断言,通过引理3-2,引理3-4与(2.1),可得
与
令
,结论成立。
定理3-1 令假设2-3中条件1),2),3)成立,
是由Hanke-Raus准则(1.6)所确定的。进一步若源条件(2.5)成立,且
充分小,那么
其中的
仅依赖于
且
。
证明:通过引理3-5证明的
的表达式,可得
其中
。
由于引理3-5成立,再根据引理3-5中
的估计证明过程,可知
,
估计成立。
因此应用引理3-4以及
充分小(可认为
),可得
其中
。
基于
的表达式,应用引理3-4,(2.1),(2.2)以及(3.6)可得
其中
。
由插值不等式以及上面的断言可得
其中
。
基于
的表达式,根据以上的断言,引理3-4与(2.1),可得
根据(3.1)以及(3.6),可得
以及
所以结论得证。
4. 总结与展望
本文针对带有扰动算子的非线性反问题提出了一种Landweber迭代法,并在一定的假设条件下,证明了此迭代法在不依赖于噪音水平的Hanke-Raus准则下的收敛阶。由于该方法在每步迭代过程中需要计算
,我们可借鉴文献 [7] 简化的思想,将每步计算的导数固定为某个合适的导数来减少计算量。另一方面,我们也可引入加速收敛的策略(比如采用文献 [8] 的两点梯度格式等)提出新的Landweber迭代格式。证明这些迭代法在Hanke-Raus准则下的收敛阶将是我们后续的工作。
基金项目
国家自然科学基金(批准号:12071184)资助项目。
参考文献