设C和Q分别为Hilbert空间H1和H2的非空闭凸子集, A:H1→H2是一个有界线性算子. 分裂可行性问题(SFP)是指: 找一点x∈C满足
由于此问题在信号处理、图像重构以及其他应用领域的广泛应用而备受关注, 详见文献 [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]及其参考文献.目前有许多迭代算法求解(1.1),其中最流行的一种算法是Byrne 的CQ 算法,它以下述方式给定
2009年,Censor和Segal[16]考虑了分裂公共不动点问题(SCFPP): 找一点x∈Fix(U)满足
如果在(1.3)和(1.4)式中分别取U=PC和V=PQ则(1.3)和(1.4)式分别化为(1.1)和(1.2)式, 因此,分裂公共不动点问题(SCFPP)是分裂可行性问题(SFP)的推广(一般化), 而相应的迭代算法(1.4)是(1.2)的推广(一般化).
Li 等[14]最近考虑了带约束的分裂公共不动点问题(CSCFPP)如下: 找一点x∈C∩Fix(T)满足
注1.1 如果我们取C=H1与Q=H2,则(1.5)式化为(1.3)式. 因此带约束的分裂公共不动点问题(CSCFPP)是分裂公共不动点问题(SCFPP)的一般化, 而相应的迭代算法(1.6)是(1.4)的一般化.
我们用Γ表示问题(1.5)的解集,即 Γ={x∗∈H1:x∗∈C∩Fix(T),Ax∗∈Q∩Fix(S)}. Li 等[14]证明了,如果Γ≠∅,{αn}满足下列条件 (C1)αn→0;(C2)∞∑n=1αn=∞;(C3)αn+1αn→1, 则由(1.6)式所产生的序列 {xn}强收敛于x∗=PΓ(θ),这里PΓ表示由H到Γ上 的距离投影.
目前,我们关心的问题如下
(1) Li 等[14]的上述收敛结果能否由已知的收敛定理推得?
(2) 如何修改算法(1.6)式以便在没有假设(C3)αn+1αn→1的情况下保证有强收敛?
本文的目的是要引入三种迭代算法并分析它们的收敛性,所得结果改进并扩展了某些作者的相关结论.
设H是实Hilbert空间,用⟨⋅ ⋅⟩ 和‖⋅‖ 分别表示H中的内积和范数.
众所周知,∀x,y∈H有下列恒等式成立
设C是H中的非空闭凸子集,则∀x∈H,存在唯一的元素PCx∈C满足
称PC:H→C为从H到C上的距离投影,它享有良好的几何性质
(p1)⟨x−PCx,y−PCx⟩≤0, ∀x∈H, ∀y∈C;
(p2)‖PCx−PCy‖2≤⟨x−y,PCx−PCy⟩,∀x,y∈H,
特别地,有 ‖PCx−PCy‖≤‖x−y‖,∀x,y∈H; (p3)‖PCx−PCy‖2≤‖x−y‖2−‖(I−PC)x−(I−PC)y‖2,∀x,y∈H,
特别地,有
(p4)‖PCx−y‖2≤‖x−y‖2−‖(I−PC)x‖2,∀x∈H y∈C;
(p5)PC=12I+12S, 这里S是一个非扩张映像.
我们用dom(S)表示S的定义域,Fix(S)表示S的不动点集.
回顾映像S:dom(S)⊂H→H称为非扩张的,如果
熟知,如果S是非扩张的且Fix(S)≠∅,则有
映像S:dom(S)⊂H→H称为firmly -非扩张的,如果存在λ∈(0,1) 与另一个非扩张映像T满足
注2.1 由(p2)、(p3)和 (p5)可知,PC和I−PC 都是firmly -非扩张的与12 -平均非扩张的. 平均非扩张的重要性包含在下列三个命题之中.
命题2.1 (Byrne[4]) 如果U和V分别是 λ1 -平均非扩张的和λ2 - 平均非扩张的, 只要T=UV是适定的,则T是λ -平均非扩张的,这 里λ=λ1+λ2−λ1λ2.
命题2.2 设U是平均非扩张的,而V是非扩张的, 满足Fix(U)⋂Fix(V)≠∅. 只要UV和 VU是适定的,就有 Fix(UV)=Fix(VU)=Fix(U)⋂Fix(V).
证 先证Fix(UV)=Fix(U)⋂Fix(V).
“右边⊆左边"是明显的,只需证明”左边⊆右边".
设x=UVx,p∈Fix(U)⋂Fix(V), 因为U是平均非扩张的,故存在λ∈(0,1)及另外一个非扩张映像W满足
使用(2.1)式,x=UVx,(2.7)式和(2.8)式得 ‖x−p‖2=‖UVx−p‖2=‖[(1−λ)+λW]Vx−p‖2=‖(1−λ)Vx+λWVx−[(1−λ)Vp+λWVp]‖2=‖(1−λ)(Vx−Vp)+λ(WVx−WVp)‖2=(1−λ)‖Vx−Vp‖2+λ‖WVx−WVp‖2−λ(1−λ)‖Vx−WVx‖2, 这推得
使用(2.7)式和(2.9)式得 x=UVx=(1−λ)Vx+λWVx=(1−λ)Vx+λVx=Vx,
这推得x=Ux,从而有x∈Fix(U)⋂Fix(V).
使用一个类似的推理过程可证Fix(VU)=Fix(V)⋂Fix(U).
由于PC是12 -平均非扩张的,故使用命题2.2可以推得下述结论.
命题2.3 (Zhou等[15]) 设T:dom(T)⊇C→H是非扩张映像满足 Fix(T)≠∅,则有 Fix(PCT)=Fix(T)=Fix(TPC).
命题2.4 (Suzuki[13]) 设{xn}和{yn}是Banach空间 E中的两个有界序列,{λn} 是(0,1)中的数列满足条件 0<lim_nλn≤¯limnλn<1. 如果{xn}和{yn}满足下述关系 xn+1=λnxn+(1−λn)yn, n≥1, 则xn−yn→θ (n→∞).
命题2.5 (Xu[12]) 设{an}是一个非负数列,满足下述循环关系 an+1≤(1−tn)an+δn, n≥1, 其中tn∈[0,1],而δn∈R1满足条件 (i) ∞∑n=1tn=∞; (ii) ¯limnδntn≤0或∞∑n=1|δn|<∞, 则有an→0 (n→∞).
下面的命题在本文中起着关键的作用.
命题2.6 设C和Q分别为实Hilbert空间H1 和H2的非空闭凸子集,A:H1→H2是一个有界线性算子. 设T:dom(T)⊇C→H1和S:dom(S)⊇Q→H2分别为非扩张映像.任取固定的 δ∈(0,1‖A‖2),令U=I−δA∗(I−SPQ)A,记 A−1(Q∩Fix(S))={x∈H1:Ax∈Q∩Fix(S)}, 则下列结论成立
(1) Γ=C∩Fix(T)∩A−1(Q∩Fix(S));
(2) 若Γ≠∅,则Fix(U)=A−1(Q∩Fix(S));
(3) 若δ∈(0,1‖A‖2],则U是非扩张的; 若δ∈(0,1‖A‖2),则U是 δ‖A‖2 -平均非扩张的.
证 (1) 由Γ的定义即知.
(2) A−1(Q∩Fix(S))⊆Fix(U)是明显的, 只需证明相反的包含关系.设x=Ux,取定一点z∈Γ, 则A∗(I−SPQ)Ax=θ,Az∈Q且Az=SAz, 这推得PQAz=Az且A∗(I−SPQ)Az=θ,从而有
(3) 使用U的定义与(2.2)式,∀x,y∈H1,我们有
现在我们估计(2.16)式右边第2项.
再次使用(2.2)式和(p3)得
下面我们引入三种迭代算法求解(1.6)式.固定λ∈(0,1), 记Tλ=(1−λ)IλT.
算法A x1∈H1, u∈H1, xn+1=TPC[(1−αn)Uxn+αnu], n≥1.
算法B x1∈H1, u∈H1, xn+1=TλPC[(1−αn)Uxn+αnu], n≥1.
算法C x1∈H1, u∈H1, xn+1=λxn+(1−λ)TPC[(1−αn)Uxn+αnu], n≥1.
定理3.1 假设Γ≠∅, δ∈(0,1‖A‖2). 如果{αn} 满足下列条件
(C1) αn→0;
(C2) ∞∑n=1αn=∞;
(C3) ∞∑n=1αn+1−αn|<∞ 或αn+1αn→1, 则由算法A所产生的序列{xn}依范数收敛于x∗=PΓu.
证 由命题2.2、2.3和2.6知,U是平均非扩张的,且有 Γ=C∩Fix(T)∩Fix(U)=Fix(TPC)∩Fix(U)=Fix(UTPC).
令yn=(1−αn)Uxn+αnu, 则算法A化为xn+1=TPCyn, 从而产生出一个新的迭代序列{yn}
使用Wittmann[11]和Xu[12]的收敛定理,我们断言 yn→x∗ (n→∞), 而x∗=PFix(V)u=PΓu, 从而有xn→x∗ (n→∞).
如果在算法A中,取u=θ,则我们可得下述推论.
推论3.1 假设Γ≠∅,δ∈(0,1‖A‖2). 如果 {αn}满足定理3.1中条件(C1)--(C3). 设{xn}是由下列方式产生的序列
注3.1 推论3.1是文献[14]中的定理3.2.值得指出的 是我们所给出的收敛性分析比文献[14]中的方法要简单得多.
定理 3.2 假设Γ≠∅, δ∈(0,1/‖A‖2). 如果{αn}满足条件 (C1) αn→0 (n→∞); (C2) ∞∑n=1αn=∞, 则由算法B所产生的序列{xn}依范数收敛于x∗=PΓu. 特别地,如果我们取u=θ,则相应的迭代序列{xn}依范数收敛于 x∗=PΓθ.
证 由命题2.2、2.3和2.6知,U是平均非扩张的,而且 Γ=C∩Fix(T)∩Fix(U)=C∩Fix(Tλ)∩Fix(U)=Fix(TλPC)∩Fix(U)=Fix(TλPCU). 令yn=(1−αn)Uxn+αnu, 则算法B化为xn+1=TλPCyn, 由此产生一个新的序列{yn}
记βn=αn+1,W=UTλPC,则{βn}满足条件 (C1)与(C2),由命题2.1 知 W是平均非扩张的,而(3.4)式化为
使用Suzuki[13]的收敛定理,我们断言yn→x∗, 而x∗=PFix(W)u=PΓu, 从而xn→x∗ (n→∞). 特别地, 如果取u=θ,则迭代序列{xn}依范数收敛于x∗=PΓθ.
定理 3.3 假设Γ≠∅,δ∈(0,1/‖A‖2). 如果{αn}满足条件(C1)与(C2), 则由算法C所产生的序列{xn}依范数收敛于x∗=PΓu, 特别地,如果我们取u=θ,则相应的迭代序列{xn}依范数收敛于 x∗=PΓθ.
证 我们分四步完成证明.
第1步 证明{xn}是有界的.
使用算法C,T, PC, U的非扩张性得 ‖xn+1−x∗‖≤λ‖xn−x∗‖+(1−λ)(1−αn)‖xn−x∗‖+(1−λ)αn‖u−x∗‖=[1−(1−λ)αn]‖xn−x∗‖+(1−λ)αn‖u−x∗‖≤max{‖x1−x∗‖,‖u−x∗‖}=M, ∀n≥1. 这表明{xn}是有界的,从而{Uxn}也是有界的.
第2步 证明xn+1−xn→θ (n→∞).
令yn=TPC[(1−αn)Uxn+αnu],则算法C可写为
使用T, PC, U的非扩张性得 ‖yn+1−yn‖=‖TPC[(1−αn+1)Uxn+1+αn+1u]−TPC[(1−αn)Uxn+αnu]‖≤‖xn+1−xn‖+|αn+1−αn|(‖u‖+M), 这推得 ¯limn(‖yn+1−yn‖−‖xn+1−xn‖)≤0.
使用命题2.4,我们断言yn−xn→θ (n→∞), 从而有xn+1−xn→θ (n→∞).
由于 ‖xn−TPCUxn‖≤‖xn−yn‖+‖yn−TPCUxn‖≤‖xn−yn‖+αn(‖u‖+‖Uxn‖)→0, 故
第3步 证明¯limn⟨u−x∗xn−x∗⟩≤0,其中x∗=PΓu.
由第1步知{xn}是有界的,不失一般性,我们可以假设
第4步 证明xn→x∗ (n→∞).
由于U是δ‖A‖2 -平均非扩张的,故存在另一个非扩张映像G满足
使用算法C,‖⋅‖2的凸性,T与PC的非扩张性和(3.11)式得
使用第2步,αn→0和(3.12)式可推得xn−Gxn→θ (n→∞).
再由(3.9)式推得xn−Uxn→θ (n→∞),从而由第3步知
记tn=(1−λ)αn,δn=2tn(1−αn)⟨u−x∗,Uxn−x∗⟩+(1−λ)αn2‖u−x∗‖2,则(3.12)式可重写为
注 3.2 如果我们取T=I|H1,S=I|H2, 则算法A,B和C分别化作求解分裂可行性问题(SFP)(1.1)的相应算法.
注 3.3 由于TPC和SPQ未必是定向算子, 定理3.1--3.3不能由文献[9]推出.
本节我们将把第3节所得到的结果扩展到T和S都是λ -严格伪压缩映像的场合. T:dom(T)⊂H→H称为λ -严格伪压缩的,如果存在λ>1使得
易知,(4.1)式等价于
引理 4.1 设λ∈(−∞,1),T:dom(T)⊂H→H是λ -严格伪压缩映像. 记Tθ=(1−θ)I+θT,则下列结论成立
(1) 当θ∈(0,1−λ]时,Tθ是非扩张的;
(2) 当θ∈(0,1−λ)时,Tθ 是θ1−λ -平均非扩张的;
(3) 当θ>0时,Fix(Tθ)=Fix(T).
证 (1) 由于θ∈(0,1−λ],使用(2.1)和(4.1)式得 ‖Tθx−Tθy‖2=‖(1−θ)(x−y)+θ(Tx−Ty)‖2=(1−θ)‖x−y‖2+θ‖Tx−Ty‖2−θ(1−θ)‖(I−T)x−(I−T)y‖2≤(1−θ)‖x−y‖2+θ‖x−y‖2−θ(1−λ−θ)‖(I−T)x−(I−T)y‖2≤‖x−y‖2, ∀x,y∈dom(T), 这推得 ‖Tθx−Tθy‖≤‖x−y‖, ∀x,y∈dom(Tθ).
(2) 设θ∈(0,1−λ),则θ1−λ∈(0,1),Tθ可以重写为
由(1)知λI+λT是非扩张的,故(4.3)式表明Tθ是θ1−λ -平均非扩张的.
(3) 由于I−Tθ=θ(I−T),θ>0,故有 x∈Fix(Tθ)⇔x∈Fix(T). 证毕.
我们考虑涉及严格伪压缩映像的带约束的分裂公共不动点问题: 找一点x∈C⋂Fix(T)满足
我们用Ω表示问题(4.4)的解集,即
对θ∈(1,1−λ],ν∈(0,1−μ], 记Tθ=(1−θ)I+θT,Sν=(1−ν)I+νS.
我们引入下面的算法: u,x1∈H1,β∈(0,1),
使用定理3.1--3.3,我们可以推得下列结果.
定理4.1 假设Ω≠∅, θ∈(0,1−λ],ν∈(0,1−μ],{αn} 满足条件 (C1) αn→0; (C2) ∞∑n=1αn=∞; (C3) ∞∑n=1αn+1−αn|<∞或 αn+1αn→1, 则由(4.6)产生的序列{xn}依范数收敛于x∗=PΩu. 特别地,如果取u=θ,则{xn}依范数收敛于 x∗=PΩθ.
证 使用引理4.1(1)和(3)知,Tθ:dom(Tθ)⊇C→H1是非扩张的, Sv:dom(Sv)⊇Q→H2是非扩张的, 而且Fix(Tθ)=Fix(T),Fix(Sv)=Fix(S).记 Uv=I−δA∗(I−SvPQ)A, 则由命题2.6知,Uv是δ‖A‖2 -平均非扩张的, 且有 Fix(Uv)=A−1(Q∩Fix(Sv))=A−1(Q∩Fix(S)).
使用定理3.1知,{xn}依范数收敛于x∗=PΩu,其中 Ω=C∩Fix(T)∩A−1(Q∩Fix(S))=C∩Fix(Tθ)∩A−1(Q∩Fix(Sv))=C∩Fix(Tθ)∩A−1(Q∩Fix(Uv)).
特别地,如果取u=θ,则{xn}依范数收敛于x∗=PΩθ.
定理4.2 假设Ω≠∅,θ∈(0,1−λ), ν∈(0,1−μ],{αn}满足条件 (C1) αn→0; (C2) ∞∑n=1αn=∞, 则由(4.6)式产生的序列{xn}依范数收敛于x∗=PΩu. 特别地,如果取u=θ,则{xn}依范数收敛于x∗=PΩθ.
证 此时Tθ:dom(Tθ)⊇C→H1 是θ1−λ -平均非扩张的, Sv:dom(Sv)⊇Q→H2是非扩张的, Uv=I−δA∗(I−SvPQ)A是δ‖A‖2 -平均非扩张的,而且Fix(Tθ)=Fix(T), Fix(Sv)=Fix(S).再注意到 Ω=C∩Fix(T)∩A−1(Q∩Fix(S))=C∩Fix(Tθ)∩A−1(Q∩Fix(Sv))=C∩Fix(Tθ)∩A−1(Q∩Fix(Uv)).
使用定理3.2推知,{xn}依范数收敛于x∗=PΩu. 特别地,如果取u=θ, 则{xn}依范数收敛于x∗=PΩθ.
定理4.3 假设Ω≠∅,θ∈(0,1−λ), ν∈(0,1−μ],{αn}满足条件 (C1) αn→0; (C2) ∞∑n=1αn=∞, 则由(4.7)式产生的序列 {xn}依范数收敛于x∗=PΩu. 特别地, 如果取u=θ,则{xn}依范数收敛于x∗=PΩθ.
证 使用引理4.1(2)和(3)知,Tθ:dom(Tθ)⊇C→H1是非扩张的,Sν:dom(Sv)⊇Q→H2是非扩张的,而且Fix(Tθ)=Fix(T), Fix(Sν)=Fix(S).再注意到 Ω=C∩Fix(Tθ)∩A−1(Q∩Fix(Sv))=C∩Fix(Tθ)∩A−1(Q∩Fix(Uv)). 使用定理3.3推得,{xn}依范数收敛于x∗=PΩu, 特别地,如果取u=θ,则{xn}依范数收敛于x∗=PΩθ.