设H 是具有内积· , · 和范数· 的实Hilbert空间,C 是H 上的一个非空闭凸子集,映射A : H → H 是一个映射,研究的变分不等式问题是指寻找x * ∈ C ,使得A ( x * ) , x - x * ≥ 0 , x ∈ C ,记其解集为VI (C ,A )。
近年来,求解Hilbert空间中变分不等式问题和不动点问题受到了众多学者的关注
[1 -12 ] 。杨静等
[13 ] 提出了一种投影算法:
y n = P C ( x n - λ n A x n ) , z n = y n + λ n ( A x n - A y n ) , q n = ( 1 - β n ) z n + β n T z n , x n + 1 = α n f ( x n ) + 1 - α n q n ,
在该算法中,
λ n = γ l m n ,
m n 是满足
γ l m A x n - A y n ≤ μ x n - y n 的最小非负整数,并且
T : H → H 是
α - 半压缩映射,其中,
A 是一个伪单调满足一致连续的映射,该算法证明了产成的迭代序列的强收敛性。2023年,叶明露等
[14 ] 提出了求解准单调且Lipschitz连续的变分不等式问题,投影算法如下:
x 0 ∈ R n , α > 0,0 < η < 1,0 < σ < 1 , k = 0 , α η m k F ( x k ) - F ( P C ( x k - α η m k F ( x k ) ) ) ≤ σ x k - P C ( x k - α η m k F ( x k ) ) , r ( x k , λ k ) = x k - y k , y k = P C ( x k - λ k F ( x k ) ) , h k ( v ) = x k - y k - λ k ( F ( x k ) - F ( y k ) ) , v - y k , x k + 1 = P H k ( y k - λ k ( F ( y k ) - F ( x k ) ) ) ,
其中,λ k = α η m k ,H k = v ∈ R n : h k ( v ) ≤ 0 ,映射F : R n → R n 是连续的,映射F 在R n 是准单调映射,其在准单调且对偶变分不等式解集非空的条件下,得到了准单调变分不等式解集的压缩投影算法所生成的聚点是解的结论。
最近,高兴慧等
[15 ] 提出了求解伪单调变分不等式解集和半压缩映射有限族公共不动点集的公共元的新投影算法:
ω n = x n + σ n ( x n - x n - 1 ) , t n = P C ( ω n - λ n A ω n ) , u n = t n - λ n ( A t n - A ω n ) , y n 1 = ( 1 - β n 1 ) u n + β n 1 T 1 u n , y n 2 = ( 1 - β n 2 ) y n 1 + β n 2 T 2 y n 1 , … y n N = ( 1 - β n N ) y n N - 1 + β n N T N y n N - 1 , x n + 1 = α n f ( x n ) + ( 1 - α n ) y n N ,
其中,A 是伪单调、L -Lipschitz的,T i : H → H 是半压缩映射族,其自适应满足
λ n + 1 = m i n μ ω n - t n A ω n - A t n , λ n , A ω n - A t n ≠ 0 ; λ n , A ω n - A t n = 0 。
受到上述研究的启发,本文提出了一种新的惯性自适应迭代算法,并证明了该算法的强收敛性,最后给出数值实验验证本文算法收敛速度更快,且得到的结论从伪单调映射推广至准单调映射,进而改进和推广了文献[
15 ]的相关结果。
1 预备知识
假设H 是一个实的Hilbert空间,T : H → H 是一个非线性映射,如果条件T x * = x * 能够满足,那么x * 就叫做映射T 的一个不动点,并用F i x ( T ) 表示映射T 的不动点集,即满足:F i x ( T ) = x * ∈ H , T x * = x * 。假设C 是H 上的一个非空闭凸子集,满足序列x n ⊂ C ,用x n ⇀ x 表示x n 弱收敛于x ,用x n → x 表示x n 强收敛于x 。如果对任意的x ∈ H ,可以在C 中找到唯一的逼近点,记作P C x ,满足
其中,P C 叫做H 到C 上的映射。
定义1[16 ] 假设
T : H → H 是一个非线性映射,且满足
F i x ( T ) ≠ ∅ ,有
1)若对任意的x ∈ H , z ∈ F i x ( T ) ,满足
T x - z 2 ≤ x - z 2 + β ( I - T ) x 2 ,0 ≤ β < 1 ,
则称T 是β - 半压缩映射。
2)若对任意的x , y ∈ H 以及存在常数L > 0 ,满足
则称T 是L -Lipschitz连续的。
3)若对任意的x , y ∈ H ,满足
则称T 是伪单调的。
4)若对任意的x , y ∈ H ,满足
则称T 是准单调的。
5)若对每一个序列x n ⊂ H , x n 弱收敛到某个点x ∈ H ,有T x n 弱收敛到T x ,则称T 在H 上是序列弱连续的。
注1[14 ] 映射是伪单调的一定是准单调的,反之则不成立。
定义2[16 ] 假设
T : H → H 是一个非线性映射,
F i x ( T ) ≠ ∅ ,如果对任意的序列
x n ⊂ H , 当满足
x n ⇀ x 和
I - T x n → 0 的条件时,有
x ∈ F i x T 成立,那么称
I - T 在零点是半闭的。
引理1[16 ] 假设
C 是
H 上的非空闭凸子集,对于给定的
x ∈ H ,
p ∈ C ,有
引理2[17 ] 假设
T : H → H 是一个
β -半压缩映射且满足映射
T 的不动点集非空,那么
T 的不动点集是
H 中的闭凸集。
引理3[18 ] 假设
c n 是(0,1)上的一个实数序列,
∑ n = 1 ∞ θ n = ∞ ,
d n 是一个实数序列,假设
当序列c n 的每一个子序列c n k 均满足l i m k → ∞ i n f ( c n k + 1 - c n k ) ≥ 0 时,有l i m k → ∞ s u p d n k ≤ 0 ,则l i m n → ∞ c n = 0 。
2 主要结果
在给出主要结果之前,先给出以下条件:
(C1) A : H → H 为准单调、满足L -Lipschitz连续以及序列弱连续,f : H → H 是一个具有常数σ σ ∈ 0,1 的压缩映射;
(C2) T i : H → H 是满足∩ i = 1 N F ( T i ) ≠ ∅ 的ε i -半压缩映射,I - T i 在零点是半闭的,其中,i = 1,2 , … , N ;
(C3) V I ( C , A ) ⋂ ⋂ F i = 1 N ( T i ) ≠ ∅ ;
(C4)序列α n ,β n i ,η n ⊆ 0,1 且l i m n → ∞ α n = 0 ,∑ n = 1 ∞ α n = ∞ ,0 < α i ≤ β n i ≤ 1 - ε i ≤ 1 ,i = 1,2 , … , N ,a i > 0 。
算法1 λ 1 > 0 , μ ∈ ( 0,1 ) ,取x 0 , x 1 ∈ H , n : = 1 。
第1步 令ω n = x n + η n ( x n - x n - 1 ) 。
第2步 计算
u n = P C ( ω n - k λ n A ω n ) ,
y n = u n - k λ n ( A u n - A ω n ) ,其中,
λ n + 1 = m a x A ω n - A u n μ ω n - u n , λ n , A ω n - A u n ≠ 0 ; λ n , A ω n - A u n = 0 。
第3步 计算
t n 1 = ( 1 - β n 1 ) y n + β n 1 T 1 y n ,
t n 2 = ( 1 - β n 2 ) t n 1 + β n 2 T 2 t n 1 ,
…
t n N - 1 = ( 1 - β n N - 1 ) t n N - 2 + β n N - 1 T N - 1 t n N - 2 ,
t n N = ( 1 - β n N ) t n N - 1 + β n N T N t n N - 1 。
第4步 计算
令n : = n + 1 ,回到第1步。
引理4[19 ] 如果满足条件(C1),那么由自适应步长产成的序列
λ n 非递减且
引理5[4 ] 如果满足条件(C1)~(C4),
ω n 是算法1产成的序列,假设存在子列
ω n k 弱收敛到
z ∈ H 和
l i m n → ∞ u n k - ω n k = 0 ,则
A z =0或者
z ∈ VI (
C ,
A )。
引理6 如果满足条件(C1)~(C4),y n 是由算法产成的序列,则
y n - q 2 ≤ ω n - q 2 - ( 1 - μ 2 k 2 λ n + 1 2 λ n 2 ) u n - ω n 2 , [ q ∈ V I ( C , A ) ] 。 (1)
证明 由算法1可知
根据y n 的定义和式(2) 知
y n - q 2 = u n - k λ n ( A u n - A ω n ) - q 2 = u n - q 2 + k λ n 2 A u n - A ω n 2 - 2 k λ n A u n - A ω n , u n - q = u n - ω n + ω n - q 2 + k λ n 2 A u n - A ω n 2 - 2 k λ n A u n - A ω n , u n - q = u n - ω n 2 + ω n - q 2 + 2 u n - ω n , ω n - q +
k λ n 2 A u n - A ω n 2 - 2 k λ n A u n - A ω n , u n - q = u n - ω n 2 + ω n - q 2 - 2 u n - ω n , u n - ω n + 2 u n - ω n , u n - q + k λ n 2 A u n - A ω n 2 - 2 k λ n A u n - A ω n , u n - q = u n - ω n 2 + ω n - q 2 - 2 u n - ω n 2 + 2 u n - ω n , u n - q + k λ n 2 A u n - A ω n 2 - 2 k λ n A u n - A ω n , u n - q = ω n - q 2 - u n - ω n 2 + 2 u n - ω n , u n - q + k λ n 2 A u n - A ω n 2 - 2 k λ n A u n - A ω n , u n - q 。 (3)
由引理1以及u n = P C ( ω n - k λ n A ω n ) ,可知
ω n - k λ n A ω n - u n , u n - q ≥ 0 ,即
u n - ω n , u n - q ≤ - k λ n A ω n , u n - q ,(4)
由q ∈ VI (C ,A ),得A q , u n - q ≥ 0 。
当A q , u n - q > 0 时,由A 是准单调映射,有
将式(4) 代入式(3) 并结合式(5) ,可得
y n - q 2 ≤ ω n - q 2 - u n - ω n 2 - 2 k λ n A ω n , u n - q + k λ n 2 A u n - A ω n 2 - 2 k λ n A u n - A ω n , u n - q ≤ ω n - q 2 - u n - ω n 2 - 2 k λ n A u n , u n - q + k λ n 2 A u n - A ω n 2 ≤ ω n - q 2 - u n - ω n 2 - 2 k λ n A u n , u n - q + μ 2 k 2 λ n + 1 2 λ n 2 u n - ω n 2 ≤ ω n - q 2 - 1 - μ 2 k 2 λ n + 1 2 λ n 2 u n - ω n 2 - 2 k λ n A u n , u n - q ≤ ω n - q 2 - 1 - μ 2 k 2 λ n + 1 2 λ n 2 u n - ω n 2 ,
式(1) 得证。
当A q , u n - q = 0 时,设对任意的k ∈ ( 0 , + ∞ ) ,其子列记为k i ,有1 k i > 0 ,则有
0 < A q , u n - q + 1 k i = A q , u n + A q A q 2 k i - q ,
结合A 是准单调得出
A u n + A q A q 2 k i , u n + A q A q 2 k i - q ≥ 0 ,
令i → ∞ ,则有A u n , u n - q ≥ 0 ,与上述证明一致,式(1) 得证。
定理1 如果满足条件(C1)~(C4),l i m n → ∞ η n α n · x n - x n - 1 = 0 ,那么由算法1产成的迭代序列x n 强收敛到q ∈ V I ( C , A ) ⋂ ∩ i = 1 N F ( T i ) ,q ∈ P V I ( C , A ) ⋂ ∩ i = 1 N F ( T i ) f ( q ) 。
证明 由引理2,得到
∩ i = 1 N F ( T i ) 是闭凸集,由文献[
20 ]得到
VI (
C ,
A )也是闭凸集,再由条件(C3),注意到
P V I ( C , A ) ⋂ ∩ i = 1 N F ( T i ) f ( · ) 是压缩映射,根据Banach压缩原理,可知存在
q ∈ V I ( C , A ) ⋂ ∩ i = 1 N F ( T i ) ,使得
q = P V I ( C , A ) ⋂ ∩ i = 1 N F ( T i ) f ( q ) 。
步骤1 证明序列x n 有界。由引理4和0 < μ < 1 ,可得l i m n → ∞ 1 - μ 2 k 2 λ n + 1 2 λ n 2 = 1 - μ 2 k 2 > 0 。
由t n 1 , t n 2 , … , t n N 的构造,可得
t n N - q 2 = ( 1 - β n N ) ( t n N - 1 - q ) + β n N ( T N t n N - 1 - q ) 2 = ( 1 - β n N ) t n N - 1 - q 2 + β n N T N t n N - 1 - q 2 - β n N ( 1 - β n N ) t n N - 1 - T N t n N - 1 2 ≤ ( 1 - β n N ) t n N - 1 - q 2 + β n N t n N - 1 - q 2 + ε N ( I - T N ) t n N - 1 2 - β n N ( 1 - β n N ) t n N - 1 - T N t n N - 1 2 = t n N - 1 - q 2 - β n N ( 1 - ε N - β n N ) t n N - 1 - T N t n N - 1 2 , (7)
t n N - 1 - q 2 = ( 1 - β n N - 1 ) ( t n N - 2 - q ) + β n N - 1 ( T N - 1 t n N - 2 - q ) 2 = ( 1 - β n N - 1 ) t n N - 2 - q 2 + β n N - 1 T N - 1 t n N - 2 - q 2 - β n N - 1 ( 1 - β n N - 1 ) · t n N - 2 - T N - 1 t n N - 2 2 ≤ ( 1 - β n N - 1 ) t n N - 2 - q 2 + β n N - 1 t n N - 2 - q 2 + ε N - 1 ( I - T N - 1 ) t n N - 2 2 - β n N - 1 ( 1 - β n N - 1 ) t n N - 2 - T N - 1 t n N - 2 2 = t n N - 2 - q 2 - β n N - 1 ( 1 - ε N - 1 - β n N - 1 ) t n N - 2 - T N - 1 t n N - 2 2 , (8)
…
t n 1 - q 2 = ( 1 - β n 1 ) ( y n - q ) + β n 1 ( T 1 y n - q ) 2 = ( 1 - β n 1 ) y n - q 2 + β n 1 T 1 y n - q 2 - β n 1 ( 1 - β n 1 ) y n - T 1 y n 2 ≤ ( 1 - β n 1 ) y n - q 2 + β n 1 y n - q 2 + ε 1 ( I - T 1 ) y n 2 - β n 1 ( 1 - β n 1 ) y n - T 1 y n 2 = y n - q 2 - β n 1 ( 1 - ε 1 - β n 1 ) y n - T 1 y n 2 。 (9)
联立式(7)~式(9) 以及条件(C4),可得
t n N - q ≤ t n N - 1 - q ≤ ⋅ ⋅ ⋅ ≤ t n 1 - q 2 ≤ y n - q 。 (10)
又ω n - q = x n - q + η n ( x n - x n - 1 ) ≤
由于l i m n → ∞ η n α n x n - x n - 1 = 0 ,则存在M 1 > 0 ,使得对任何n ≥ 1 ,均有η n α n x n - x n - 1 ≤ M 1 ,则有
结合式(6) 、式(10) 和式(11) ,可得
t n N - q ≤ t n N - 1 - q ≤ … ≤ t n 1 - q 2 ≤ y n - q ≤ ω n - q ≤ x n - q + α n M 1 ,
由x n + 1 的构造,可得
x n + 1 - q = α n f ( x n ) - q + ( 1 - α n ) ( t n N - q ) ≤ α n f ( x n ) - q + ( 1 - α n ) t n N - q = α n f ( x n ) - f ( q ) + f ( q ) - q + ( 1 - α n ) t n N - q ≤ α n σ x n - q + α n f ( q ) - q + ( 1 - α n ) x n - q + α n M 1 ≤ α n σ x n - q + α n f ( q ) - q + ( 1 - α n ) x n - q + α n M 1 = 1 - α n ( 1 - σ ) x n - q + α n ( 1 - σ ) f ( q ) - q + M 1 ( 1 - σ ) ≤ m a x x n - q , f ( q ) - q + M 1 ( 1 - σ ) ≤ … ≤ m a x x 0 - q , f ( q ) - q + M 1 ( 1 - σ ) ,
故序列x n 有界。
步骤2 证明
( 1 - α n ) 1 - λ n + 1 2 λ n 2 μ 2 k 2 u n - ω n 2 + ( 1 - α n ) β n 1 ( 1 - ε 1 - β n 1 ) y n - T 1 y n 2 + ( 1 - α n ) ∑ i = 2 N β n i ( 1 - ε i - β n i ) t n i - 1 - T i t n i - 1 2 ≤ x n - q 2 - x n + 1 - q 2 + α n M 4 。 (12)
由式(11) ,可得
ω n - q 2 ≤ x n - q + α n M 1 2 = x n - q 2 + α n α n M 1 2 + 2 M 1 x n - q 2 ≤ x n - q 2 + α n M 2 ,(13)
其中,M 2 = s u p n ≥ 1 ( α n M 1 2 + 2 M 1 x n - q ) 。根据· 2 的凸性以及式(1) 、式(7)~式(11) 、式(13) 可得
x n + 1 - q 2 = α n f ( x n ) - q + ( 1 - α n ) ( t n N - q ) 2 = α n f ( x n ) - f ( q ) + f ( q ) - q 2 + ( 1 - α n ) t n N - q 2 - α n ( 1 - α n ) f ( x n ) - t n N 2 ≤ α n f ( x n ) - f ( q ) 2 + 2 f ( q ) - q , f ( x n ) - q + ( 1 - α n ) t n N - 1 - q 2 - β n N ( 1 - ε N - β n N ) t n N - 1 - T N t n N - 1 2 ≤ α n f ( x n ) - f ( q ) 2 + 2 f ( q ) - q , f ( x n ) - q + ( 1 - α n ) t n N - 2 - q 2 - β n N - 1 ( 1 - ε N - 1 - β n N - 1 ) t n N - 2 - T N - 1 t n N - 2 2 - β n N ( 1 - ε N - β n N ) t n N - 1 - T N t n N - 1 2 … ≤ α n f ( x n ) - f ( q ) 2 + 2 f ( q ) - q , f ( x n ) - q + ( 1 - α n ) y n - q 2 - β n 1 ( 1 - ε 1 - β n 1 ) y n - T 1 y n 2 - ∑ i = 2 N β n i ( 1 - ε i - β n i ) t n i - 1 - T i t n i - 1 2 ≤ α n f ( x n ) - f ( q ) 2 + 2 f ( q ) - q , f ( x n ) - q + ( 1 - α n ) x n - q 2 + α n M 2 - 1 - λ n + 1 2 λ n 2 μ 2 k 2 · u n - ω n 2 - β n 1 ( 1 - ε 1 - β n 1 ) y n - T 1 y n 2 - ∑ i = 2 N β n i ( 1 - ε i - β n i ) t n i - 1 - T i t n i - 1 2 ≤ α n f ( x n ) - f ( q ) 2 + 2 f ( q ) - q , f ( x n ) - q +
( 1 - α n ) x n - q 2 - 1 - λ n + 1 2 λ n 2 μ 2 k 2 u n - ω n 2 - β n 1 ( 1 - ε 1 - β n 1 ) y n - T 1 y n 2 - ∑ i = 2 N β n i ( 1 - ε i - β n i ) t n i - 1 - T i t n i - 1 2 + α n M 2 ≤ α n x n - q 2 + 2 α n f ( q ) - q f ( x n ) - q + ( 1 - α n ) x n - q 2 - ( 1 - α n ) 1 - λ n + 1 2 λ n 2 μ 2 k 2 · u n - ω n 2 - ( 1 - α n ) β n 1 ( 1 - ε 1 - β n 1 ) y n - T 1 y n 2 - ( 1 - α n ) ∑ i = 2 N β n i ( 1 - ε i - β n i ) t n i - 1 - T i t n i - 1 2 + α n M 2 ≤ x n - q 2 - ( 1 - α n ) 1 - λ n + 1 2 λ n 2 μ 2 k 2 ⋅ u n - ω n 2 - ( 1 - α n ) β n 1 ( 1 - ε 1 - β n 1 ) y n - T 1 y n 2 - ( 1 - α n ) ∑ i = 2 N β n i ( 1 - ε i - β n i ) t n i - 1 - T i t n i - 1 2 + α n M 2 + α n M 3 ,
其中,M 3 = s u p n ≥ 1 2 f ( q ) - q f ( x n ) - q ,令M 4 = M 2 + M 3 。故式(12) 成立。
步骤3 证明
x n + 1 - q 2 ≤ 1 - α n ( 1 - σ ) x n - q 2 + α n ( 1 - σ ) M 5 η n α n ( 1 - σ ) x n - x n - 1 + 2 1 - σ f ( q ) - q , x n + 1 - q 。(14)
由x n + 1 、 ω n 的构造可得
x n + 1 - q 2 = α n f ( x n ) - f ( q ) + α n f ( q ) - q + ( 1 - α n ) ( t n N - q ) 2 ≤ α n f ( x n ) - f ( q ) + ( 1 - α n ) ( t n N - q ) 2 + 2 α n f ( q ) - q , x n + 1 - q ≤ α n σ x n - q 2 + ( 1 - α n ) t n N - q 2 + 2 α n f ( q ) - q , x n + 1 - q , (15)
t n N - q 2 ≤ ω n - q 2 ≤ x n - q + η n x n - x n - 1 2 ≤ x n - q 2 + 2 η n x n - x n - 1 , ω n - q ≤ x n - q 2 + 2 η n x n - x n - 1 ω n - q ≤ x n - q 2 + η n x n - x n - 1 M 5 , (16)
其中,M 5 = s u p n ≥ 1 2 ω n - q 。将式(16) 代入式(15) ,可得
x n + 1 - q 2 ≤ α n σ x n - q 2 + ( 1 - α n ) x n - q 2 + η n x n - x n - 1 M 5 + 2 α n f ( q ) - q , x n + 1 - q ≤ α n σ x n - q 2 + ( 1 - α n ) x n - q 2 + η n x n - x n - 1 M 5 + 2 α n f ( q ) - q , x n + 1 - q ≤ 1 - α n ( 1 - σ ) x n - q 2 + α n ( 1 - σ ) M 5 η n α n ( 1 - σ ) x n - x n - 1 + 2 1 - σ f ( q ) - q , x n + 1 - q ,
故式(14) 成立。
步骤4 证明x n - q 2 → 0 ( n → ∞ ) 。
根据引理3可知,只要证明当序列x n - q 中的每个子列x n k - q 均能满足l i m k → ∞ i n f x n k + 1 - q - x n k - q ≥ 0 时,有l i m k → ∞ s u p f ( q ) - q , x n k + 1 - q ≤ 0 。
首先,假设序列x n k - q 是x n - q 的子列,且满足l i m k → ∞ i n f x n k + 1 - q - x n k - q ≥ 0 ,就有
l i m k → ∞ i n f x n k + 1 - q 2 - x n k - q 2 = l i m k → ∞ i n f x n k + 1 - q + x n k - q · x n k + 1 - q - x n k - q ≥ 0 。
由式(12) ,可得
l i m k → ∞ s u p ( 1 - α n k ) 1 - λ n k + 1 2 λ n k 2 μ 2 k 2 u n k - ω n k 2 + ( 1 - α n k ) β n k 1 ( 1 - ε 1 - β n k 1 ) y n k - T 1 y n k 2 + ( 1 - α n k ) ∑ i = 2 N β n k i ( 1 - ε i - β n k i ) t n k i - 1 - T i t n k i - 1 2 ≤ l i m k → ∞ s u p ( x n k - q 2 - x n k + 1 - q 2 + α n M 4 ) = - l i m k → ∞ i n f x n k + 1 - q 2 - x n k - q 2 ≤ 0 ,
于是有
l i m k → ∞ u n k - ω n k = 0 , (17)
l i m k → ∞ y n k - T 1 y n k = 0 , (18)
l i m k → ∞ t n k i - T i + 1 t n k i = 0 , i = 1,2 , ⋅ ⋅ ⋅ , N - 1 。 (19)
由l i m n → ∞ η n α n x n - x n - 1 = 0 和式(17) ,可得
u n k - x n k ≤ u n k - ω n k + ω n k - x n k = u n k - ω n k + α n k η n k α n k x n k - x n k - 1 → 0 ( k → ∞ ) 。(20)
由式(17)~式(20) 以及条件(C4),可得
x n k + 1 - x n k = x n k + 1 - t n k N + t n k N - t n k N - 1 + … + t n k 2 - t n k 1 + t n k 1 - y n k + y n k - x n k = α n k f ( x n k ) - t n k N + β n k N t n k N - 1 - T N t n k N - 1 + … + β n k 2 t n k 1 - T 2 t n k 1 + β n k 1 y n k - T 1 y n k + u n k - x n k + λ n k + 1 λ n k μ k u n k - ω n k → 0 ( k → ∞ ) 。 (21)
因为序列x n k 是有界的,所以存在子列x n k j ,使得x n k j ⇀ z ∈ H , 且满足
l i m k → ∞ s u p f ( q ) - q , x n k - q = l i m j → ∞ f ( q ) - q , x n k j - q = f ( q ) - q , z - q 。 (22)
联立式(17)~式(20) ,可得
y n k - x n k ≤ u n k - x n k + λ n k + 1 λ n k μ k u n k - ω n k → 0 ( k → ∞ ) , t n k 1 - x n k ≤ t n k 1 - y n k + y n k - x n k = β n k 1 y n k - T 1 y n k + y n k - x n k → 0 ( k → ∞ ) , t n k 2 - x n k ≤ t n k 2 - t n k 1 + t n k 1 - x n k = β n k 2 t n k 1 - T 2 t n k 1 + t n k 1 - x n k → 0 ( k → ∞ ) , … t n k N - x n k ≤ t n k N - t n k N - 1 + t n k N - 1 - x n k = β n k N t n k N - 1 - T N t n k N - 1 + t n k N - 1 - x n k → 0 ( k → ∞ ) 。
综上,可得
l i m k → ∞ t n k 1 - x n k = 0 , l i m k → ∞ t n k 2 - x n k = 0 , … , l i m k → ∞ t n k N - x n k = 0 , l i m k → ∞ y n k - x n k = 0 。 (23)
注意到x n k ⇀ z ,由x n k - ω n k → 0 可知ω n k ⇀ z ( n → ∞ ) ,根据式(17) 以及引理5,得到z ∈ VI (C ,A )。另一方面,根据式(23) ,可得y n k ⇀ z , t n k 1 ⇀ z , … , t n k N - 1 ⇀ z , 再结合式(18) 和式(19) ,根据定义2,有z ∈ F ( T i ) ( i = 1,2 , … , N ) ,于是有z ∈ ∩ i = 1 N F ( T i ) ( i = 1,2 , … , N ) 。因此z ∈ V I ( C , A ) ⋂ ∩ i = 1 N F ( T i ) 。根据q ∈ p V I ( C , A ) ⋂ ∩ i = 1 N F ( T i ) f ( q ) 和引理1,可得
l i m k → ∞ s u p f ( q ) - q , x n k - q = f ( q ) - q , z - q ≤ 0 , (24)
由式(21) 和式(24) ,可得
l i m k → ∞ s u p f ( q ) - q , x n k + 1 - q ≤ l i m k → ∞ s u p f ( q ) - q , x n k + 1 - x n k + l i m k → ∞ s u p f ( q ) - q , x n k - q ≤ 0 ,(25)
联立式(14) 、式(25) 和引理3,可得l i m n → ∞ x n - q 2 = 0 。
证毕。
注2 本文将文献[
15 ]算法中的最小值自适应步长改为最大值自适应步长,并在算法[
15 ]步骤2的自适应步长上增加了系数,从而加快了本文算法的收敛速度;本文将文献[
15 ]的伪单调映射推广为准单调映射。
3 数值实验
在数值实验中,用n 代表程序的迭代步数,t 代表运行时间,x n + 1 - x n 代表测量第n 步误差。
例1 假设
A : R 2 → R 2 是一个单调且满足
L -Lipschitz连续的映射,其中,
L -Lipschitz常数是
L = M ,满足
A ( x ) =
M x + q ,其中,
q ∈ R 2 ,
M =
RR T +
Q +
N ,定义
C = x = x 1 , x 2 ∈ R 2 x i ≤ 2 , i = 1,2 为非空必凸集。
R 是
2 × 2 阶矩阵,且
R 的元素在(-2,2)中随机生成;
Q 是
2 × 2 阶斜对称矩阵,其中,
Q 的元素同样在(-2,2)中随机生成;
N 是
2 × 2 阶对角元素非负的对角矩阵,其中,对角元素在(0,2)中随机生成。对于本次数值实验,以
x 0 = x 1 = 1,1 T 为初始值,
q = ( 0,0 ) T 。令
f ( x ) = 1 20 x ,定义映射
T i : R 2 → R 2 ( i = 1,2 , 3,4 ) ,分别为
T 1 x = - x 2 ,
T 2 x = - x 6 ,
T 3 x = x 3 ,
T 4 x = x 8 ,将迭代终止条件
eps 分别设置为
x n + 1 - x n ≤ 10 - 50 ,
x n + 1 - x n ≤ 10 - 60 ,
x n + 1 - x n ≤ 10 - 70 。将本文算法与文献[
10 ]算法进行比较,选取参数如下:
本文算法A l g : k = 1.2 ,α n = 1 n + 1 3 ,λ 1 = 0.5 L ,η n = 1 n + 1 ,μ = 0.99 ,β n 1 = n 2 n + 1 ,β n 2 = n n + 2 ,β n 3 = n 3 n + 2 ,β n 4 = n 5 n + 5 ;
文献[
10 ]算法
A l g . F :
α n = 1 n + 1 3 ,
λ 1 = 0.5 L ,
σ n = 1 n + 1 ,
μ = 0.99 ,
β n 1 = n 2 n + 1 ,
β n 2 = n n + 2 ,
β n 3 = n 3 n + 2 ,
β n 4 = n 5 n + 5 。
终止条件
e p s = x n + 1 - x n ≤ 10 - 50 时,本文算法对比文献[
15 ],其算法收敛速度更快,数值实验对比结果见
图1 。
注3 根据程序实验,随着迭代步数增加,本文算法程序的误差逐渐减小并且趋向0,验证了本文算法的可行性和有效性。同时,通过对比
图1 可以看出随着迭代次数增加,本文算法的迭代收敛速度优于文献[
15 ]算法的迭代收敛速度。
4 结束语
本文证明了所构造的惯性自适应迭代序列,强收敛到准单调变分不等式解集和半压缩映射有限族公共不动点集的公共元。在接下来的研究中,可以进一步考虑将半压缩映射不动点与均衡问题相结合,并证明其强收敛性。
国家自然科学基金项目(61866038)
陕西省大学生创新训练计划项目(202510719024)
延安大学研究生教育创新计划项目(YCX2024046)