1引言

正如我们所知,最初的分裂可行性问题(SFP)是由Censor和Elfving首先引入的[1]自1994年成立以来,备受关注。这是由于它在信号处理和图像重建方面的应用,特别是在调强放射治疗方面的进展;请看[26].

由于SFP是凸可行性问题(CFP)的一个特例,即在有限多个闭集和凸集的非空交点上寻找一个点,因此我们接下来简要回顾了一些与CFP相关的历史方法。CFP是一个重要的问题,因为工程和数学中的许多实际反演或估计问题都可以纳入这个框架;请参见,例如,组合框[7]博斯克和博文[8]和Kiwiel[9]. 传统上,求解CFP的迭代投影方法使用凸集上的正交投影(,关于欧氏距离函数的最近点投影);请参见,例如, [1014]. Bregman提出的广义距离函数及其相关的广义投影已经做了很多工作[15].

1994年,Censor和Elfving[1]研究了在求解SFP的单个迭代过程中使用不同类型的广义投影。他们的建议是一种迭代算法,其中涉及矩阵逆的计算,这是一项众所周知的困难任务。这就是为什么拜恩[16,17]提出了所谓的CQ公司该算法通过具有适当步长的递归过程生成序列。这个CQ公司算法只涉及到对集合的投影的计算C类因此,在这些投影具有封闭形式表达式的情况下是可以实现的(例如,C类是闭合的球或半场)。关于CQ公司文献方法;例如,请参见[1834]. 然而,我们必须指出,步长的确定取决于算子(矩阵)范数(或矩阵乘积的主要特征值)。这意味着为了实现CQ公司算法,首先必须计算(或至少估计)算子的矩阵范数,这在实践中通常不是一件容易的工作。

为了克服上述困难,开发了一种允许自适应选择步长的所谓自适应方法。注意,该方法是Goldstein投影方法的应用[35]莱维汀和波利亚克[36]一个合适的变分不等式问题,这是求解变分不等式的最简单的数值方法之一。然而,这种投影方法的效率在很大程度上取决于步长参数的选择。如果选择一个小参数以保证迭代序列的收敛,那么递归会导致收敛速度较慢。另一方面,如果选择较大的步长来提高收敛速度,则生成的序列可能不会收敛。在解决变分不等式问题的实际应用中,Lipschitz常数可能很难估计,即使底层映射是线性的,例如SFP。根据原有的Goldstein-Levitin-Polyak方法,发展了一些求解变分不等式问题的自适应方法[35,36]. 请参阅,例如, [3745].

受自适应策略的激励,张等。[45]提出了使用可变步长代替Censor中固定步长的方法等。[46]. 此外,Zhao和Yang还介绍了一种自适应投影方法[29]通过使用类似Armijo的搜索,它被采用。这些算法的优点在于,既没有矩阵范数的先验信息A类也没有任何其他条件A类是必需的,并且仍然保证收敛。

本文进一步发展和改进了求解SFP的自适应方法。提出了一种改进的自适应方法来求解SFP。作为一种特殊情况,SFP的最小范数解可以迭代求解。

2框架和初步结果

H(H) 1 H(H) 2 是两个希尔伯特空间,让C类是的两个闭凸子集 H(H) 1 H(H) 2 分别是。A类: H(H) 1 H(H) 2 是一个有界线性算子。拆分可行性问题(SFP)是要找到一个点 x个 这样的话

x个 C类A类 x个 .
(1)

接下来,我们用Γ表示SFP的解集,,Γ={x个C类:A类x个}.

1994年,Censor和Elfving[1]研究了在求解SFP的单个迭代过程中使用不同类型的广义投影。他们是第一个提出以下涉及逆运算的算法的人 A类 1 :

x个 k个 + 1 = A类 1 P(P) ( P(P) A类 ( C类 ) ( A类 x个 k个 ) ) ,k个0,

哪里C类是封闭的,并且凸集在 n个 ,同时A类是一个完整的等级n个×n个矩阵和A类(C类)={ n个 =A类x个,x个C类}。请注意 A类 1 不容易执行。因此,拜恩[16,17]提出了所谓的问题生成序列的算法{ x个 n个 }通过递归过程

x个 n个 + 1 = P(P) C类 ( x个 n个 τ n个 A类 ( P(P) ) A类 x个 n个 ) ,
(2)

其中台阶尺寸 τ n个 在间隔中选择(0,2/ A类 2 )值得注意的是CQ公司算法只涉及投影的计算 P(P) C类 P(P) 到集合上C类,因此在以下情况下是可实现的 P(P) C类 P(P) 具有封闭式表达式(例如,C类是闭合的球或半场)。然而,我们观察到,台阶尺寸的确定 τ n个 取决于算子(矩阵)范数A类(或最大特征值 A类 A类). 这意味着为了实际实施问题算法,必须首先计算(或至少估计)A类这在实践中通常不是一件容易的任务。

为了克服上述困难,所谓的自适应方法允许步长 τ n个 开发了自适应选择。如果我们设置

(f)(x个):= 1 2 A类 x个 P(P) A类 x个 2 ,

然后是凸目标(f)是可微的,并且具有由下式给出的Lipschitz梯度

(f)(x个)= A类 ( P(P) )A类.

因此CQ公司算法(2)可以通过最小化以下凸极小化问题得到

最小值 x个 C类 (f)(x个).
(3)

我们知道这一点 x个 C类是问题(3)的驻点,如果它满足

(f) ( x个 ) , x个 x个 0,x个C类.
(4)

因此,我们可以使用下面的梯度投影算法来求解(SFP)

x个 n个 + 1 = P(P) C类 ( x个 n个 τ n个 (f) ( x个 n个 ) ) ,
(5)

哪里 τ n个 ,迭代时的步长n个,在间隔中选择(0,2/L(左)),其中L(左)是的Lipschitz常数(f).

必须将上述方法(5)视为Goldstein投影方法的应用[35]莱维汀和波利亚克[36]变分不等式问题(4)是求解变分不等式的最简单的数值方法之一。然而,这种投影方法的效率在很大程度上取决于参数的选择 τ n个 .一个小的 τ n个 保证了迭代序列的收敛性,但递归导致收敛速度较慢。另一方面,较大的步长将提高收敛速度,但生成的序列可能不会收敛。在解决变分不等式问题的实际应用中,Lipschitz常数可能很难估计,即使底层映射是线性的,例如SFP。

张中的方法等。[45]和审查员等。[46]提出了解决多集分裂可行性问题的方法。

算法2.1S1.给定一个非负序列 τ n个 这样的话 n个 = 0 τ n个 <,δ(0,1),μ(0,1),ρ(0,1),ϵ>0, β 0 >0,和任意初始点 x个 0 .设置 γ 0 = β 0 n个=0.

S2.求最小非负整数 n个 这样的话 β n个 + 1 = μ k个 γ k个

x个 n个 + 1 = P(P) C类 ( x个 n个 β n个 + 1 (f) ( x个 n个 ) ) ,

这满足了

β n个 + 1 (f) ( x个 n个 ) (f) ( x个 n个 + 1 ) 2 (2δ) x个 n个 x个 n个 + 1 , (f) ( x个 n个 ) (f) ( x个 n个 + 1 ) .

S3.如果

β n个 + 1 (f) ( x个 n个 ) (f) ( x个 n个 + 1 ) 2 ρ x个 n个 x个 n个 + 1 , (f) ( x个 n个 ) (f) ( x个 n个 + 1 ) ,

然后设置 γ n个 + 1 =(1+ τ n个 + 1 ) β n个 + 1 ; 否则,设置 γ n个 + 1 = β n个 + 1 .

S4.如果e(电子)( x个 n个 , β n个 )ϵ,停止;否则,设置n个:=n个+1然后转到S2。

赵和杨介绍了以下自适应投影方法[29]通过使用类似Armijo的搜索,它被采用。

算法2.2给定常数β>0,σ(0,1)γ(0,1).让 x个 0 随心所欲。对于n个=0,1, , 计算

x个 n个 + 1 = P(P) C类 ( x个 n个 τ n个 (f) ( x个 n个 ) ) ,

哪里 τ n个 =β γ n个 n个 是最小的非负整数这样的话

(f) ( P(P) C类 ( x个 n个 β γ (f) ( x个 n个 ) ) ) (f)( x个 n个 )σ (f) ( x个 n个 ) , x个 n个 P(P) C类 ( x个 n个 β γ (f) ( x个 n个 ) ) .

算法2.1和算法2.2的优点在于,关于矩阵范数的先验信息A类也没有任何其他条件A类是必需的,并且仍然保证收敛。

我们将介绍我们改进的自适应方法来求解SFP。在这方面,我们需要立即引入配料。

C类是实Hilbert空间的非空闭凸子集H(H).A映射T型:C类C类被称为非扩张如果

T型x个T型x个,x个,C类.

地图ψ:C类C类据说是δ-如果存在常数,则为收缩δ[0,1)这样的话

ψ ( x个 ) ψ ( ) δx个,x个,C类.

回忆一下(最近点或公制)投影H(H)到上面C类,表示为 P(P) C类 ,分配给每个x个H(H)独特的点 P(P) C类 (x个)C类使用该属性

x个 P(P) C类 ( x个 ) =inf公司 { x个 : C类 } .

众所周知,度量投影 P(P) C类 属于H(H)到上面C类具有以下基本属性:

  1. (a)

    P(P) C类 (x个) P(P) C类 ()x个为所有人x个,H(H);

  2. (b)

    x个, P(P) C类 (x个) P(P) C类 () P(P) C类 ( x个 ) P(P) C类 ( ) 2 对于每个x个,H(H);

  3. (c)

    x个 P(P) C类 (x个), P(P) C类 (x个)0为所有人x个H(H)C类.

接下来,我们采用以下符号:

  • x个 n个 x个意味着 x个 n个 强烈收敛于x个;

  • x个 n个 x个意味着 x个 n个 弱收敛到x个;

  • ω w个 ( x个 n个 ):={x个: x个 n个 j个 x个}是弱者ω-序列的极限集{ x个 n个 }.

回想一下一个函数(f):H(H)称为凸,如果

(f) ( λ x个 + ( 1 λ ) ) λ(f)(x个)+(1λ)(f)(),λ(0,1),x个,H(H).

众所周知,一个可微函数(f)是凸的当且仅当以下关系成立:

(f)(z(z))(f)(x个)+ (f) ( x个 ) , z(z) x个 ,z(z)H(H).

回忆一下一个元素H(H)据说是的次梯度(f):H(H)x个如果

(f)(z(z))(f)(x个)+,z(z)x个,z(z)H(H).

如果函数(f):H(H)在处至少有一个次梯度x个,据说在x个.的次梯度集(f)在这一点上x个称为的次微分(f)x个,表示为(f)(x个).A函数(f)如果它是次可微分的,则称为次可微分x个H(H).如果(f)是凸的且可微的,则其梯度和次梯度重合。A函数(f):H(H)称为弱下半连续(w-lsc)x个如果 x个 n个 x个暗示

(f)(x个) lim信息 n个 (f)( x个 n个 ).

(f)据说是w-lsc onH(H)如果每个点都是w-lscx个H(H).

第一个引理很容易证明。

引理2.1[14]

(f)(x个):= 1 2 A类 x个 P(P) A类 x个 2 .然后

  1. (i)

    (f) 是凸的和可微的;

  2. (ii)

    (f) 是w-lsc打开 C类.

引理2.2[47]

鉴于 x个 H(H) 1 .然后 x个 解决了SFP协议当且仅当 x个 求解不动点方程

x个 = P(P) C类 ( x个 γ A类 ( P(P) ) A类 x个 ) ,

哪里 γ>0.

引理2.3[48]

假设 { n个 } 是一个非负实数序列,使得

n个 + 1 (1 γ n个 ) n个 + δ n个 ,

哪里 { γ n个 } 是中的序列 (0,1) { δ n个 } 是这样一个序列

  1. (1)

    n个 = 1 γ n个 =;

  2. (2)

    酸橙酱 n个 δ n个 γ n个 0 n个 = 1 | δ n个 |<.

然后 n个 n个 =0.

引理2.4[49]

{ n个 } 是一个不在无穷远处减少的实数序列,在存在子序列的意义上 { n个 } 属于 { n个 } 这样的话 n个 n个 + 1 为所有人 0.对于每个 n个 n个 0 ,定义整数序列 {τ(n个)} 作为

τ(n个)=最大值{k个n个: n个 < n个 + 1 }.

然后 τ(n个) 作为 n个 以及所有人 n个 n个 0

最大值{ τ ( n个 ) , n个 } τ ( n个 ) + 1 .

3主要成果

在本节中,我们陈述并证明了我们的主要结果。

C类是实Hilbert空间的非空闭凸子集 H(H) 1 H(H) 2 分别是。ψ:C类 H(H) 1 成为δ-收缩δ[0, 2 2 ).让A类: H(H) 1 H(H) 2 是一个有界线性算子。

算法3.1对于给定的 x个 0 C类,假设{ x个 n个 }已建造。如果(f)( x个 n个 )=0,然后停止并 x个 n个 是SFP(1)的解决方案。否则,继续并计算 x个 n个 + 1 通过递归

x个 n个 + 1 = P(P) C类 [ α n个 ψ ( x个 n个 ) + ( 1 α n个 ) ( x个 n个 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) ) ] ,n个0,
(6)

哪里{ α n个 }(0,1){ ρ n个 }(0,2).

定理3.1 假设SFP协议是一致的,那就是,Γ.假设以下条件成立:

  1. (a)

    n个 α n个 =0 n个 = 1 α n个 =;

  2. (b)

    inf公司 n个 ρ n个 (2 ρ n个 )>0.

然后 { x个 n个 } 由定义(6)强烈收敛于 z(z),它解决了以下变分不等式:

z(z)Γ这样的话 z(z) ψ ( z(z) ) , z(z) x个 0 为所有人 x个Γ.
(7)

证明首先,很明显,变分不等式(7)的解是唯一的(由于ψ根据变分不等式中的相关结果),表示为z(z).然后z(z)= P(P) Γ (ψ(z(z)))。我们可以假设{ x个 n个 }是无限的,即算法3.1不会在有限次数的迭代中终止。因此,(f)( x个 n个 )0为所有人n个。从(6)开始,我们有

x个 n个 + 1 z(z) 2 = P(P) C类 [ α n个 ψ ( x个 n个 ) + ( 1 α n个 ) ( x个 n个 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) ) ] z(z) 2 α n个 ( ψ ( x个 n个 ) z(z) ) + ( 1 α n个 ) ( x个 n个 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) z(z) ) 2 α n个 ψ ( x个 n个 ) z(z) 2 + ( 1 α n个 ) x个 n个 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) z(z) 2 ( 1 α n个 ) [ x个 n个 z(z) 2 + ρ n个 2 (f) 2 ( x个 n个 ) (f) ( x个 n个 ) 2 2 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) , x个 n个 z(z) ] + α n个 ( ψ ( x个 n个 ) ψ ( z(z) ) + ψ ( z(z) ) z(z) ) 2 .
(8)

通过凸性(f)(引理2.1)和事实(f)(z(z))=0对于z(z)Γ,我们推断

(f)( x个 n个 )=(f)( x个 n个 )(f)(z(z)) (f) ( x个 n个 ) , x个 n个 z(z) .
(9)

使用不等式 ( + b条 ) 2 2( 2 + b条 2 )为所有人,b条,我们有

( ψ ( x个 n个 ) ψ ( z(z) ) + ψ ( z(z) ) z(z) ) 2 2 ψ ( x个 n个 ) ψ ( z(z) ) 2 + 2 ψ ( z(z) ) z(z) 2 2 δ 2 x个 n个 z(z) 2 + 2 ψ ( z(z) ) z(z) 2 .
(10)

从(8)到(10),我们得到

x个 n个 + 1 z(z) 2 ( 1 α n个 ) [ x个 n个 z(z) 2 ρ n个 ( 2 ρ n个 ) (f) 2 ( x个 n个 ) (f) ( x个 n个 ) 2 ] + 2 δ 2 α n个 x个 n个 z(z) 2 + 2 α n个 ψ ( z(z) ) z(z) 2 [ 1 ( 1 2 δ 2 ) α n个 ] x个 n个 z(z) 2 + ( 1 2 δ 2 ) α n个 2 ψ ( z(z) ) z(z) 2 1 2 δ 2 最大值 { x个 n个 z(z) 2 , 2 ψ ( z(z) ) z(z) 2 1 2 δ 2 } .

通过归纳,我们推断

x个 n个 + 1 z(z) 2 最大值 { x个 0 z(z) 2 , 2 ψ ( z(z) ) z(z) 2 1 2 δ 2 } .

因此,{ x个 n个 }有界。

通过使用的坚定非扩张性 P(P) C类 ,我们得出

x个 n个 + 1 z(z) 2 = P(P) C类 [ α n个 ψ ( x个 n个 ) + ( 1 α n个 ) ( x个 n个 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) ) ] P(P) C类 z(z) 2 α n个 ψ ( x个 n个 ) z(z) , x个 n个 + 1 z(z) + ( 1 α n个 ) x个 n个 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) z(z) , x个 n个 + 1 z(z) = α n个 ψ ( x个 n个 ) ψ ( z(z) ) , x个 n个 + 1 z(z) + α n个 ψ ( z(z) ) z(z) , x个 n个 + 1 z(z) + ( 1 α n个 ) x个 n个 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) z(z) , x个 n个 + 1 z(z) α n个 δ x个 n个 z(z) x个 n个 + 1 z(z) + α n个 ψ ( z(z) ) z(z) , x个 n个 + 1 z(z) + ( 1 α n个 ) x个 n个 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) z(z) x个 n个 + 1 z(z) = ( α n个 δ x个 n个 z(z) + ( 1 α n个 ) x个 n个 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) z(z) ) x个 n个 + 1 z(z) + α n个 ψ ( z(z) ) z(z) , x个 n个 + 1 z(z) 1 2 ( α n个 δ x个 n个 z(z) + ( 1 α n个 ) x个 n个 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) z(z) ) 2 + 1 2 x个 n个 + 1 z(z) 2 + α n个 ψ ( z(z) ) z(z) , x个 n个 + 1 z(z) .

由此可见

x个 n个 + 1 z(z) 2 ( α n个 δ x个 n个 z(z) + ( 1 α n个 ) x个 n个 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) z(z) ) 2 + 2 α n个 ψ ( z(z) ) z(z) , x个 n个 + 1 z(z) α n个 δ 2 x个 n个 z(z) 2 + ( 1 α n个 ) x个 n个 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) z(z) 2 + 2 α n个 ψ ( z(z) ) z(z) , x个 n个 + 1 z(z) α n个 δ 2 x个 n个 z(z) 2 + ( 1 α n个 ) [ x个 n个 z(z) 2 ρ n个 ( 2 ρ n个 ) (f) 2 ( x个 n个 ) (f) ( x个 n个 ) 2 ] + 2 α n个 ψ ( z(z) ) z(z) , x个 n个 + 1 z(z) = [ 1 ( 1 δ 2 ) α n个 ] x个 n个 z(z) 2 + 2 α n个 ψ ( z(z) ) z(z) , x个 n个 + 1 z(z) ( 1 α n个 ) ρ n个 ( 2 ρ n个 ) (f) 2 ( x个 n个 ) (f) ( x个 n个 ) 2 .
(11)

接下来,我们将证明这一点 x个 n个 z(z)遵循中的想法[49]. 设置 n个 = x个 n个 z(z) 2 为所有人n个0.自 α n个 0 inf公司 n个 ρ n个 (2 ρ n个 )>0,在不失一般性的情况下,我们可以假设(1 α n个 ) ρ n个 (2 ρ n个 )σ对一些人来说σ>0因此,我们可以将(11)重写为

n个 + 1 n个 + ( 1 δ 2 ) α n个 n个 + σ (f) 2 ( x个 n个 ) (f) ( x个 n个 ) 2 2 α n个 ψ ( z(z) ) z(z) , x个 n个 + 1 z(z) .
(12)

现在,我们考虑两种可能的情况。

案例1。假设{ n个 }最终会减少,,存在N个>0这样的话{ n个 }正在减少n个N个在这种情况下,{ n个 }必须收敛,从(12)可以看出

0 σ (f) 2 ( x个 n个 ) (f) ( x个 n个 ) 2 n个 n个 + 1 ( 1 δ 2 ) α n个 n个 + 2 α n个 ψ ( z(z) ) z(z) x个 n个 + 1 z(z) n个 n个 + 1 + M(M) α n个 ,
(13)

哪里M(M)>0是一个常数,因此 啜饮 n个 {2ψ(z(z))z(z) x个 n个 + 1 z(z)}M(M).出租n个在(13)中,我们得到

n个 (f)( x个 n个 )=0.

{ x个 n个 }有界,则存在子序列{ x个 n个 k个 }属于{ x个 n个 }弱收敛到 x个 ˜ C类.

从弱下半连续性(f),我们有

0(f)( x个 ˜ ) lim信息 k个 (f)( x个 n个 k个 )= n个 (f)( x个 n个 )=0.

因此,(f)( x个 ˜ )=0,,A类 x个 ˜ 。这表明

ω w个 ( x个 n个 )Γ.

此外,由于投影(c)的性质,

酸橙酱 n个 ψ ( z(z) ) z(z) , x个 n个 + 1 z(z) = 最大值 ω ω w个 ( x个 n个 ) ψ ( z(z) ) P(P) Γ ( ψ ( z(z) ) ) , ω P(P) Γ ( ψ ( z(z) ) ) 0.

从(12)中,我们得到

n个 + 1 [ 1 ( 1 δ 2 ) α n个 ] n个 +2 α n个 ψ ( z(z) ) z(z) , x个 n个 + 1 z(z) .
(14)

将引理2.3应用于(14),我们得到 n个 0.

案例2。假设{ n个 }最终不会减少。也就是说,存在一个整数 n个 0 这样的话 n个 0 n个 0 + 1 因此,我们可以定义一个整数序列{ τ n个 }为所有人n个 n个 0 如下:

τ(n个)=最大值{k个N个 n个 0 k个n个, k个 k个 + 1 }.

显然,τ(n个)是一个非递减序列τ(n个)+作为n个

τ ( n个 ) τ ( n个 ) + 1

为所有人n个 n个 0 在这种情况下,我们从(13)中得出

σ (f) 2 ( x个 τ ( n个 ) ) (f) ( x个 τ ( n个 ) ) 2 M(M) α τ ( n个 ) 0.

由此可见

n个 (f)( x个 τ ( n个 ) )=0.

这意味着{ x个 τ ( n个 ) }在解集Γ中;, ω w个 ( x个 τ ( n个 ) )Γ.

另一方面,我们注意到

x个 τ ( n个 ) + 1 x个 τ ( n个 ) α τ ( n个 ) ψ ( x个 τ ( n个 ) ) x个 τ ( n个 ) +(1 α τ ( n个 ) ) ρ τ ( n个 ) (f) ( x个 τ ( n个 ) ) (f) ( x个 τ ( n个 ) ) 0,

从中我们可以推断出

酸橙酱 n个 ψ ( z(z) ) z(z) , x个 τ ( n个 ) + 1 z(z) = 酸橙酱 n个 ψ ( z(z) ) z(z) , x个 τ ( n个 ) z(z) = 最大值 ω ω w个 ( x个 τ ( n个 ) ) ψ ( z(z) ) P(P) Γ ( ψ ( z(z) ) ) , ω P(P) Γ ( ψ ( z(z) ) ) 0 .
(15)

τ ( n个 ) τ ( n个 ) + 1 ,我们从(12)开始

τ ( n个 ) 2 1 δ 2 ψ ( z(z) ) z(z) , x个 τ ( n个 ) + 1 z(z) .
(16)

将(15)和(16)组合得出

酸橙酱 n个 τ ( n个 ) 0,

因此

n个 τ ( n个 ) =0.

从(14)开始,我们有

酸橙酱 n个 τ ( n个 ) + 1 酸橙酱 n个 τ ( n个 ) .

因此,

n个 τ ( n个 ) + 1 =0.

根据引理2.4,我们有

0 n个 最大值{ τ ( n个 ) , τ ( n个 ) + 1 }.

因此, n个 0也就是说, x个 n个 z(z)。这就完成了证明。□

从定理3.1,我们可以很容易地推导出以下算法和推论。

算法3.2对于给定的 x个 0 C类,假设{ x个 n个 }已建造。如果(f)( x个 n个 )=0,然后停止并 x个 n个 是SFP(1)的解决方案。否则,继续并计算 x个 n个 + 1 通过递归

x个 n个 + 1 = P(P) C类 [ ( 1 α n个 ) ( x个 n个 ρ n个 (f) ( x个 n个 ) (f) ( x个 n个 ) 2 (f) ( x个 n个 ) ) ] ,n个0,
(17)

哪里{ α n个 }(0,1){ ρ n个 }(0,2).

定理3.2 假设SFP协议是一致的,那就是,Γ.假设以下条件成立:

  1. (a)

    n个 α n个 =0 n个 = 1 α n个 =;

  2. (b)

    inf公司 n个 ρ n个 (2 ρ n个 )>0.

然后 { x个 n个 } 由定义(17)强收敛于的最小范数解SFP。

4结束语

这项工作包含了我们致力于开发和改进解决分裂可行性问题的自适应方法的研究。我们引入了改进的自适应方法来解决分割可行性问题。作为一种特殊情况,分裂可行性问题的最小范数解可以迭代求解。本研究的动机是解决许多实际问题的相关应用,这些问题在变分不等式问题领域中产生了数学模型。