Next Article in Journal
Hermite-Hadamard-Fejér Type Inequalities for Preinvex Functions Using Fractional Integrals
Next Article in Special Issue
(C , Ψ * , G ) Class of Contractions and Fixed Points in a Metric Space Endowed with a Graph
Previous Article in Journal
Multi-Product Production System with the Reduced Failure Rate and the Optimum Energy Consumption under Variable Demand
Previous Article in Special Issue
Reich, Jungck, and Berinde Common Fixed Point Results on ℱ-Metric Spaces and an Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Inertial Forward-Backward Mid-Point Methods for Sum of Infinitely Many Accretive Mappings, Variational Inequalities, and Applications

1
School of Mathematics and Statistics, Hebei University of Economics and Business, Shijiazhuang 050061, China
2
Department of Mathematics, Texas A & M University-Kingsville, Kingsville, TX 78363, USA
3
Florida Institute of Technology, Melbourne, FL 32901, USA
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(5), 466; https://doi.org/10.3390/math7050466
Submission received: 19 April 2019 / Accepted: 17 May 2019 / Published: 24 May 2019
(This article belongs to the Special Issue Fixed Point Theory and Related Nonlinear Problems with Applications)

Abstract

:
Some new inertial forward-backward projection iterative algorithms are designed in a real Hilbert space. Under mild assumptions, some strong convergence theorems for common zero points of the sum of two kinds of infinitely many accretive mappings are proved. New projection sets are constructed which provide multiple choices of the iterative sequences. Some already existing iterative algorithms are demonstrated to be special cases of ours. Some inequalities of metric projection and real number sequences are widely used in the proof of the main results. The iterative algorithms have also been modified and extended from pure discussion on the sum of accretive mappings or pure study on variational inequalities to that for both, which complements the previous work. Moreover, the applications of the abstract results on nonlinear capillarity systems are exemplified.

1. Introduction and Preliminaries

Suppose H is a real Hilbert space with norm · and inner product · , · . Let K be the non-empty closed and convex subset of H . We use → and ⇀ to denote the strong and weak convergence in H , respectively.
We know that Hilbert space H satisfies Opial’s condition in the sense that l i m i n f n x n z < l i m i n f n x n y for { x n } H with x n z and y z (see [1]).
The inclusion problem for finding u H such that
0 S u + T u
is studied intensively, where S : H H is a mapping and T : H 2 H is a multi-valued mapping. This is mainly because many problems appear in convex programming, variational inequalities, split feasibility problems, minimization problem, inverse problem and image processing can be modeled by (1).
A mapping T : D ( T ) H 2 H is said to be an accretive mapping (see [2]) if for each x , y D ( T ) , there exist u T x and v T y such that x y , u v 0 . An accretive mapping T : D ( T ) H 2 H is said to be m-accretive if R ( I + k T ) = H , for k > 0 .
A mapping S : D ( S ) H H is said to be μ -inversely strongly accretive mapping (see [3]) if for each x , y D ( S ) and μ > 0 , x y , S x S y μ S x S y 2 .
For a mapping W : D ( W ) H H , a point x D ( W ) is called a zero point of W if W x = 0 . The set of zero points of W is denoted by W 1 0 . If x D ( W ) H satisfies that W x = x , then x is called a fixed point of W. The set of fixed points of W is denoted by F i x ( W ) .
The study of the special case of inclusion problem (1), where T is accretive and S is μ -inversely strongly accretive, has been a hot topic during the past few years. In particular, the constructions of the iterative algorithms for approximating the zero point of the sum of T and S are focused, see [3,4,5,6,7,8,9,10,11,12] and the references therein. The inertial forward-backward splitting method is one of the important iterative algorithms studied by some authors, see [7,8,9,13,14].
In 2015, Lorenz and Pock [9] proposed the following inertial forward-backward algorithm for approximating zero points of T + S , where T : H 2 H is m-accretive and S : H H is μ -inversely strongly accretive:
u 0 , u 1 H chosen arbitrarily , v n = u n + θ n ( u n u n 1 ) , u n + 1 = ( I + r n T ) 1 ( v n r n S v n ) , n N .
In addition, the result that u n p ( T + S ) 1 0 , as n , is proved under some conditions.
To get strong convergence, Dong et al. proposed the following inertial forward-backward projection algorithm in Hilbert spaces in [14]:
u 0 , u 1 H chosen arbitrarily , v n = u n + α n ( u n u n 1 ) , w n = ( I + r n T ) 1 ( v n r n S v n ) , C n = { p H : w n p 2 u n p 2 2 α n u n p , u n 1 u n + α n 2 u n 1 u n 2 } , Q n = { p H : u n p , u n u 0 0 } , u n + 1 = P C n Q n ( u 0 ) , n N ,
where T and S are the same as those in (2) and P C n Q n is the metric projection whose meaning can be seen in Definition 1. The projection sets C n and Q n play an important role in the iterative construction to ensure the strong convergence. The result that u n P ( T + S ) 1 0 ( u 0 ) , as n , is proved under some conditions.
In 2018, Khan et al. proposed the following one in which the projection set Q n is deleted (see [7]):
u 0 , u 1 H chosen arbitrarily , v n = u n + θ n ( u n u n 1 ) , w n = α n u n + ( 1 α n ) ( I + r n T ) 1 ( v n r n S v n ) , C 1 : = H , C n + 1 = { p C n : w n p 2 u n p 2 2 θ n ( 1 α n ) u n p , u n 1 u n + 2 θ n 2 u n 1 u n 2 } , u n + 1 = P C n + 1 ( u 0 ) , n N ,
where T and S are the same as those in (3). The strong convergence that u n P ( T + S ) 1 0 ( u 0 ) , as n , is also obtained under some conditions.
On the other hand, the inclusion problem (1) is extended to the system of inclusion problems:
0 S i u + T i u ,
where T i is m-accretive and S i is μ i -inversely strongly accretive for i { 1 , 2 , , m } or i N . In addition, some iterative algorithms for approximating common zero points of T i + S i are constructed in [3,15,16,17]. In particular, Wei et al. proposed the following implicit mid-point forward-backward projection algorithm in [17]:
u 1 H chosen arbitrarily , w n = α n η f ( u n ) + ( I α n F ) u n , v n = β n w n + δ n i = 1 c n , i ( I + r n , i T i ) 1 [ w n + v n 2 r n , i S i ( w n + v n 2 ) ] + λ n e n , C 1 = H , C n + 1 = { p C n : v n p 2 1 + β n λ n 2 δ n w n p 2 + 2 λ n 2 δ n e n p 2 } , u n + 1 = P C n + 1 ( u 1 ) , n N ,
where f : H H is a contraction, F : H H is strongly positive linear bounded mapping, e n is the computational error, T i is m-accretive and S i is μ i -inversely strongly accretive for i N . The result that u n P i = 1 ( T i + S i ) 1 0 ( u 1 ) , as n , is proved under some conditions.
Recall that f : H H is called a contraction (see [17]) if there exists a constant l ( 0 , 1 ) such that f ( x ) f ( y ) l x y for x , y H .
A mapping F : H H is called strongly positive (see [17]) if there exists ξ > 0 such that x , F x ξ x 2 for x H . In this case,
a I b F = s u p x 1 | ( a I b F ) x , x | ,
where I is the identity mapping, a [ 0 , 1 ] and b [ 1 , 1 ] .
A mapping U : H H is said to be non-expansive (see [17]) if for each x , y H , U x U y x y .
In 2018, Wei et al. proposed some new hybrid iterative algorithms to approximate the common element of the set of zero points of infinitely many m-accretive mappings T i : H H and the set of fixed points of infinitely many non-expansive mappings B i : H H . A special case (see Corollary 3.6 in [18]) is presented as follows:
u 1 H chosen arbitrarily , y n = α n u n + ( 1 α n ) i = 1 c n , i ( I + r n , i T i ) 1 ( u n + e n ) , z n = β n u n + ( 1 β n ) i = 1 b i B i y n , C 1 = H = Q 1 , C n + 1 = { p C n : y n p 2 α n u n p 2 + ( 1 α n ) u n + e n p 2 , z n p 2 β n u n p 2 + ( 1 β n ) y n p 2 } , Q n + 1 = { p C n + 1 : u 1 p 2 u 1 P C n + 1 ( u 1 ) 2 + λ n + 1 } , u n + 1 Q n + 1 , n N .
The result that u n P m = 1 C m ( u 1 ) ( i = 1 T i 1 0 ) ( i = 1 F i x ( B i ) ) , as n , is proved under some conditions. We may notice that infinite choices of { u n } can be made, which is totally different from traditional projection iterative algorithms, e.g., (3).
In 2016, Wei et al. proposed an implicit forward-backward mid-point iterative algorithm for approximating common zero points of T i + S i , where T i is m-accretive and S i is μ i -inversely strongly accretive, for i N . A special case in [19] in the frame of Hilbert space is presented as follows:
u 1 K H chosen arbitrarily , y n = P K [ ( 1 α n ) ( u n + e n ) ] , z n = δ n y n + β n i = 1 a i ( I + r n , i T i ) 1 [ ( y n + z n 2 r n , i S i ( y n + z n 2 ) ] + ξ n e n , u n + 1 = γ n η f ( u n ) + ( I γ n F ) z n + e n , n N ,
where f and F are the same as those in (6), { e n } , { e n } and { e n } are the error sequences. Under some conditions, { u n } is proved to be convergent strongly to u 0 i = 1 ( T i + S i ) 1 0 , which also satisfies the following variational inequality:
F u 0 η f ( u 0 ) , u 0 z 0 , z i = 1 ( T i + S i ) 1 0 .
We may notice that the connection between the common element of ( T i + S i ) 1 0 for i N and the solution of one kind variational inequality is set up in [19].
In this paper, our main purpose is formulated as follows: (1) obtain strong convergence theorems instead of weak ones; (2) construct new projection sets, which ensure that infinitely many iterative sequences can be generated compared to traditional projection iterative algorithms (3), (4) and (6); (3) inject the idea of inertial forward-backward algorithm into the iterative construction, compared to iterative algorithms (6)–(8); (4) set up the connection between the common zero point of the sum of two kinds of infinitely many accretive mappings and the solution of one kind variational inequality, which complements the corresponding work since rare studies of the projection iterative algorithms (e.g., (3)–(7)) have mentioned that; (5) provide the application of the abstract result to capillarity systems.
To begin our study, we need some preliminaries.
Definition 1.
(see [2]) For the Hilbert space H and its non-empty closed and convex subset K, there exists a unique point x 0 K such that x x 0 = i n f { x y : y K } , for each x H . In this case, the metric projection mapping P K : H K is defined by P K x = x 0 , for x H .
Definition 2.
(see [20]) Let { K n } be a sequence of non-empty closed and convex subsets of H. Then
(1) 
the strong lower limit of { K n } , s l i m i n f K n , is defined as the set of all x H such that there exists x n K n for almost all n and it tends to x as n in the norm;
(2) 
the weak upper limit of { K n } , w l i m s u p K n , is defined as the set of all x H such that there exists a subsequence { K n m } of { K n } and x n m K n m for every n m and it tends to x as n m in the weak topology;
(3) 
the limit of { K n } , l i m K n , is the common value when s l i m i n f K n = w l i m s u p K n .
Lemma 1.
(see [20]) Let { K n } be a decreasing sequence of closed and convex subsets of H, i.e., K n K m if n m . Then { K n } converges in H and l i m K n = n = 1 K n .
Lemma 2.
(see [21]) Suppose H is a real Hilbert space. If l i m K n exists and is not empty, then P K n x P l i m K n x for every x H , as n .
Lemma 3.
(see [2,19]) If B : H H is accretive, then ( I + r B ) 1 : H H is non-expansive.
Lemma 4.
(see [22]) If H is a real Hilbert space with K its non-empty closed and convex subset, T i : K K is non-expansive for i N and i = 1 a i = 1 for { a i } ( 0 , 1 ) , then i = 1 a i T i is non-expansive with F i x ( i = 1 a i T i ) = i = 1 F i x ( T i ) under the assumption that i = 1 F i x ( T i ) .
Lemma 5.
(see [19]) If H is a real Hilbert space with K its non-empty closed and convex subset, S : K H is a single-valued mapping and T : H 2 H is an m-accretive mapping, then
F i x ( ( I + r T ) 1 ( I r S ) ) = ( T + S ) 1 0 ,
for r > 0 .
Lemma 6.
(see [23]) Let H be a real Hilbert space and r ( 0 , + ) . Then there exists a continuous, strictly increasing and convex function g : [ 0 , 2 r ] [ 0 , + ) with g ( 0 ) = 0 such that k x + ( 1 k ) y 2 k x 2 + ( 1 k ) y 2 k ( 1 k ) g ( x y ) , for k [ 0 , 1 ] , x , y H with x r and y r .
Lemma 7.
(see [24]) Let K be the non-empty closed and convex subset of Hilbert space H and P K : H K be the metric projection. Then
(1) 
for x H and y K , P K x y 2 + P K x x 2 y x 2 .
(2) 
y = P K x if and only if there holds the following inequality x y , y z 0 , for z K .
Lemma 8.
(see [25]) If f : H H is a contraction, then there is a unique element x H such that f ( x ) = x .

2. Some Inertial Forward-Backward Algorithms

In this section, unless otherwise stated, we always assume that:
(1)
H is a real Hilbert space;
(2)
A i : H H is μ i -inversely strongly accretive and B i : H H is m-accretive, for each i N . In addition, i = 1 ( A i + B i ) 1 0 ;
(3)
{ e n } H is the computational error;
(4)
{ σ n } , { s n , i } and { μ i } are three real number sequences in ( 0 , + ) for i , n N ;
(5)
{ α n } , { β n } and { γ n } are three real number sequences in ( 0 , 1 ) with α n + β n + γ n 1 , for n N ;
(6)
{ ω n , i } is a real number sequence in ( 0 , 1 ) with i = 1 ω n , i = 1 , for n N ;
(7)
{ k n } is a real number sequence in [ 0 , k ] for some k [ 0 , 1 ) .

2.1. New Inertial Forward-Backward Projection Algorithms

Theorem 1.
Let { u n } be generated by the following iterative algorithm:
u 0 , u 1 H c h o s e n a r b i t r a r i l y , e 1 H , v n = u n + k n ( u n u n 1 ) , w n = α n u n + β n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) v n + γ n e n , C 1 = H = Q 1 , C n + 1 = { p C n : w n p 2 ( 1 γ n ) u n p 2 + γ n e n p 2 + k n 2 u n u n 1 2 2 β n k n u n p , u n 1 u n } , Q n + 1 = { p C n + 1 : u 1 p 2 P C n + 1 ( u 1 ) u 1 2 + σ n + 1 } , u n + 1 Q n + 1 , n N .
Under the assumptions that: ( i ) s n , i μ i for i , n N ; ( i i ) σ n 0 , γ n 0 , as n ; ( i i i ) 0 < i n f n β n < 1 ; ( i v ) there exists M > 0 such that e n M , for n N , we have: u n P m = 1 C m ( u 1 ) = P i = 1 ( A i + B i ) 1 0 ( u 1 ) i = 1 ( A i + B i ) 1 0 , as n .
Proof. 
We split the proof into nine steps.
Step 1. i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) : H H is non-expansive, for n N .
The proof of Step 1 is essentially from that of Step 1 in Theorem 2.1 of [17]. For the sake of completeness, we present it as follows.
Since s n , i 2 μ i , then for each x , y H ,
( I s n , i A i ) x ( I s n , i A i ) y 2 = ( x y ) s n , i ( A i x A i y ) 2 = x y 2 2 s n , i x y , A i x A i y + s n , i 2 A i x A i y 2 x y 2 + s n , i ( s n , i 2 μ i ) A i x A i y 2 x y 2 .
Thus, ( I s n , i A i ) : H H is non-expansive, for i , n N . It then follows from Lemmas 3 and 4 that i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) : H H is non-expansive, for n N . Combining with Lemma 5, i = 1 F i x ( ( I + s n , i B i ) 1 ( I s n , i A i ) ) = i = 1 ( A i + B i ) 1 0 .
Step 2. i = 1 ( A i + B i ) 1 0 C n , for n N .
If n = 1 , it is obvious that i = 1 ( A i + B i ) 1 0 C 1 .
Now, q i = 1 ( A i + B i ) 1 0 , suppose the result is true for n = m + 1 , then if n = m + 2 , it follows from (10) that
v m + 1 q 2 = u m + 1 q 2 + 2 k m + 1 u m + 1 q , u m + 1 u m + k m + 1 2 u m + 1 u m 2 .
Using Step 1, we have:
w m + 1 q 2 α m + 1 u m + 1 q 2 + β m + 1 v m + 1 q 2 + γ m + 1 e m + 1 q 2 .
Combining (11) and (12),
w m + 1 q 2 ( α m + 1 + β m + 1 ) u m + 1 q 2 + γ m + 1 e m + 1 q 2 + k m + 1 2 u m + 1 u m 2 2 k m + 1 β m + 1 u m + 1 q , u m u m + 1 ,
which ensures that q C m + 2 . Then by induction, q C n , for n N .
Step 3. C n is a closed and convex subset of H , for each n N .
It is not difficult to see that
w n p 2 ( 1 γ n ) u n p 2 + γ n e n p 2 + k n 2 u n u n 1 2 2 β n k n u n p , u n 1 u n w n 2 ( 1 γ n ) u n 2 γ n e n 2 k n 2 u n u n 1 2 + 2 β n k n u n , u n 1 u n 2 p , w n 2 ( 1 γ n ) p , u n 2 γ n p , e n + 2 β n k n p , u n 1 u n .
Then C n is a closed and convex subset of H , for each n N .
Step 4. Q n is non-empty for each n N , which ensures that { u n } is well-defined.
From Step 3 and the definition of metric projection, for σ n + 1 , there exists δ n + 1 C n + 1 such that u 1 δ n + 1 2 ( i n f z C n + 1 u 1 z ) 2 + σ n + 1 = P C n + 1 ( u 1 ) u 1 2 + σ n + 1 . Thus, Q n + 1 , for n N . And then { u n } is well-defined.
Step 5. P C n + 1 ( u 1 ) P m = 1 C m ( u 1 ) , as n .
The proof of Step 5 is similar to Step 2 of Theorem 3.1 in [18]. It follows from Lemma 1 that l i m C n exists and l i m C n = n = 1 C n . Then Lemma 2 implies that P C n + 1 ( u 1 ) P m = 1 C m ( u 1 ) , as n .
Step 6. u n P m = 1 C m ( u 1 ) , as n .
Since u n + 1 Q n + 1 C n + 1 and C n is a convex subset of H, then for t ( 0 , 1 ) , t P C n + 1 ( u 1 ) + ( 1 t ) u n + 1 C n + 1 , which implies that
P C n + 1 ( u 1 ) u 1 t P C n + 1 ( x 1 ) + ( 1 t ) u n + 1 u 1 .
Using Lemma 6, we have
t P C n + 1 ( u 1 ) + ( 1 t ) u n + 1 u 1 2 = t ( P C n + 1 ( u 1 ) u 1 ) + ( 1 t ) ( u n + 1 u 1 ) 2 t P C n + 1 ( u 1 ) u 1 2 + ( 1 t ) u n + 1 u 1 2 t ( 1 t ) g ( P C n + 1 ( u 1 ) u n + 1 ) .
Then (13) and (14) ensure that t g ( P C n + 1 ( u 1 ) u n + 1 ) u n + 1 u 1 2 P C n + 1 ( u 1 ) u 1 2 σ n + 1 . Letting t 1 first and then n , we know that P C n + 1 ( u 1 ) u n + 1 0 as n . From Step 5, u n P m = 1 C m ( u 1 ) , as n .
Step 7. v n P m = 1 C m ( u 1 ) and w n P m = 1 C m ( u 1 ) , as n .
Since u n + 1 u n 0 and u n + 1 Q n + 1 C n + 1 , then from (10), u n + 1 v n + 1 = k n + 1 u n + 1 u n 0 , as n . Thus, v n P m = 1 C m ( u 1 ) , as n . Since u n + 1 Q n + 1 C n + 1 , then
w n u n + 1 2 ( 1 γ n ) u n + 1 u n 2 + γ n e n u n + 1 2 + k n 2 u n u n 1 2 2 β n k n u n u n + 1 , u n 1 u n ( 1 γ n ) u n + 1 u n 2 + γ n e n u n + 1 2 + k n 2 u n u n 1 2 + 2 u n + 1 u n u n 1 u n 0 ,
as n . Thus, w n P m = 1 C m ( u 1 ) , as n .
Step 8. P m = 1 C m ( u 1 ) i = 1 ( A i + B i ) 1 0 .
In fact, if, otherwise, P m = 1 C m ( u 1 ) ¯ i = 1 ( A i + B i ) 1 0 . Then P m = 1 C m ( u 1 ) i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P m = 1 C m ( u 1 ) .
Since w n = α n u n + β n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) v n + γ n e n , then β n [ i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) v n w n ] = α n ( w n u n ) + γ n ( w n e n ) 0 , as n .
Since i n f n β n > 0 , then there exists a subsequence of { n } , which is still denoted by { n } such that i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) v n w n 0 , as n .
Thus, i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) v n P m = 1 C m ( u 1 ) , as n .
Since H satisfies Opial’s condition, then
l i m i n f n v n P m = 1 C m ( u 1 ) < l i m i n f n v n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P m = 1 C m ( u 1 ) l i m i n f n v n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) v n + l i m i n f n v n P m = 1 C m ( u 1 ) l i m i n f n v n P m = 1 C m ( u 1 ) ,
which makes a contradiction! Thus, P m = 1 C m ( u 1 ) i = 1 ( A i + B i ) 1 0 .
Step 9. P m = 1 C m ( u 1 ) = P i = 1 ( A i + B i ) 1 0 ( u 1 ) .
From Step 8, P i = 1 ( A i + B i ) 1 0 ( u 1 ) u 1 P m = 1 C m ( u 1 ) u 1 .
On the other hand, since i = 1 ( A i + B i ) 1 0 m = 1 C m , then P m = 1 C m ( u 1 ) u 1 P i = 1 ( A i + B i ) 1 0 ( u 1 ) u 1 . Thus,
P i = 1 ( A i + B i ) 1 0 ( u 1 ) u 1 = P m = 1 C m ( u 1 ) u 1 .
Using Lemma 7, we have
P i = 1 ( A i + B i ) 1 0 ( u 1 ) P m = 1 C m ( u 1 ) 2 + P m = 1 C m ( u 1 ) u 1 2 P i = 1 ( A i + B i ) 1 0 ( u 1 ) u 1 2 = P m = 1 C m ( u 1 ) u 1 2 ,
which implies that P i = 1 ( A i + B i ) 1 0 ( u 1 ) = P m = 1 C m ( u 1 ) .  □
Remark 1.
Compared to (3) and (4), infinitely many m-accretive mappings and infinitely many μ i -inversely strongly accretive mappings are considered in (10). Compared to (6), the idea of inertial forward-backward algorithm is embodied in (10). Compared to (3), (4) and (6), infinite choices of the iterative sequences { u n } are defined.
Remark 2.
The traditional idea for choosing the unique iterative element u n + 1 as the metric projection of the initial element in iterative algorithms (3), (4) and (6) is contained in the ideas of (10) in our paper.
In fact, if take u n + 1 = P C n + 1 ( u 1 ) , we can easily see that
u n + 1 u 1 2 P C n + 1 ( u 1 ) u 1 2 + σ n + 1 .
Thus, u n + 1 = P C n + 1 ( u 1 ) Q n + 1 , which means that this u n + 1 is a kind of choice of (10).
Corollary 1.
If k n 0 , then (10) in Theorem 1 becomes to the traditional forward-backward iterative algorithm:
u 0 , u 1 H c h o s e n a r b i t r a r i l y , e 1 H , w n = α n u n + β n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) u n + γ n e n , C 1 = H = Q 1 , C n + 1 = { p C n : w n p 2 ( 1 γ n ) u n p 2 + γ n e n p 2 } , Q n + 1 = { p C n + 1 : u 1 p 2 P C n + 1 ( u 1 ) u 1 2 + σ n + 1 } , u n + 1 Q n + 1 , n N .
Corollary 2.
If i 1 , then (10) in Theorem 1 becomes to the following iterative algorithm:
u 0 , u 1 H c h o s e n a r b i t r a r i l y , e 1 H , v n = u n + k n ( u n u n 1 ) , w n = α n u n + β n ( I + s n B ) 1 ( I s n A ) v n + γ n e n , C 1 = H = Q 1 , C n + 1 = { p C n : w n p 2 ( 1 γ n ) u n p 2 + γ n e n p 2 + k n 2 u n u n 1 2 2 β n k n u n p , u n 1 u n } , Q n + 1 = { p C n + 1 : u 1 p 2 P C n + 1 ( u 1 ) u 1 2 + σ n + 1 } , u n + 1 Q n + 1 , n N .
Remark 3.
Let e n 0 , α n + β n 1 , for n N . After taking u n + 1 = P C n + 1 ( u 1 ) in (15), we may see that Q n + 1 can be deleted which implies that (15) reduces to (4). However, the strong assumption that n = 1 k n u n u n 1 < + in [7] is no longer needed in our paper.

2.2. New Mid-Point Inertial Forward-Backward Projection Algorithms

Theorem 2.
Suppose f : H H is a contraction with k ( 0 , 1 ) and F : H H is a strongly positive linear bounded operator with coefficient ξ > 0 . Let { u n } be generated by the following iterative algorithm:
u 0 , u 1 H c h o s e n a r b i t r a r i l y , e 1 H , z 0 = u 0 , z n = δ n λ f ( u n ) + ( I δ n F ) u n , v n = z n + k n ( z n z n 1 ) , w n = α n v n + β n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) ( v n + w n 2 ) + γ n e n , C 1 = H = Q 1 , C n + 1 = { p C n : w n p 2 2 α n + β n 2 β n z n p 2 + 2 γ n 2 β n e n p 2 + 2 α n + β n 2 β n k n 2 z n z n 1 2 2 k n 2 α n + β n 2 β n z n p , z n 1 z n } , Q n + 1 = { p C n + 1 : u 1 p 2 P C n + 1 ( u 1 ) u 1 2 + σ n + 1 } , u n + 1 Q n + 1 , n N .
Under the assumptions of ( i ) ( i v ) in Theorem 1 and ( v ) λ > 0 and ( v i ) δ n 0 , as n , we have: u n P m = 1 C m ( u 1 ) = P i = 1 ( A i + B i ) 1 0 ( u 1 ) i = 1 ( A i + B i ) 1 0 , as n .
Proof. 
We split the proof into ten steps.
Step 1. i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) : H H is non-expansive, for n N .
Copy Step 1 in Theorem 1.
Step 2. { w n } is well-defined.
Define U : H H by
U x = t y + s T ( x + y 2 ) + ( 1 t s ) v ,
where T : H H is non-expansive and t , s ( 0 , 1 ) .
It is easy to check that U is a contraction since
U x U z s x + y 2 z + y 2 = s 2 x z < x z ,
for x , z H .
In view of Lemma 8, there exists a unique element x H such that x = U x . Combining with the result of Step 1, { w n } is well-defined.
Step 3. i = 1 ( A i + B i ) 1 0 C n , for n N .
If n = 1 , it is obvious that i = 1 ( A i + B i ) 1 0 C 1 .
Now, q i = 1 ( A i + B i ) 1 0 , suppose the result is true for n = m + 1 , then if n = m + 2 , in view of Lemma 4, we have:
w m + 1 q 2 α m + 1 v m + 1 q 2 + β m + 1 v m + 1 + w m + 1 2 q 2 + γ m + 1 e m + 1 q 2 ,
which ensures that
w m + 1 q 2 2 α m + 1 + β m + 1 2 β m + 1 v m + 1 q 2 + 2 γ m + 1 2 β m + 1 e m + 1 q 2 .
It follows from (16) that
v m + 1 q 2 = z m + 1 q 2 + 2 k m + 1 z m + 1 q , z m + 1 z m + k m + 1 2 z m + 1 z m 2 .
Combining (17) and (18),
w m + 1 q 2 2 α m + 1 + β m + 1 2 β m + 1 z m + 1 q 2 + 2 γ m + 1 2 β m + 1 e m + 1 q 2 + 2 α m + 1 + β m + 1 2 β m + 1 k m + 1 2 z m + 1 z m 2 2 k m + 1 2 α m + 1 + β m + 1 2 β m + 1 z m + 1 q , z m z m + 1 .
Thus, q C m + 2 . Then by induction, q C n , for n N .
Step 4. C n is a closed and convex subset of H , for each n N . It is not difficult to see that
w n p 2 2 α n + β n 2 β n z n p 2 + 2 γ n 2 β n e n p 2 + k n 2 2 α n + β n 2 β n z n z n 1 2 2 k n 2 α n + β n 2 β n z n p , z n 1 z n w n 2 2 α n + β n 2 β n z n 2 2 γ n 2 β n e n 2 k n 2 2 α n + β n 2 β n z n z n 1 2 + 2 k n 2 α n + β n 2 β n z n , z n 1 z n 2 p , w n 2 2 α n + β n 2 β n p , z n 4 γ n 2 β n p , e n + 2 k n 2 α n + β n 2 β n p , z n 1 z n .
Then C n is a closed and convex subset of H , for each n N .
Step 5. Q n is non-empty for each n N , which ensures that { u n } is well-defined.
Copy Step 4 in Theorem 1.
Step 6. P C n + 1 ( u 1 ) P m = 1 C m ( u 1 ) , as n .
Copy Step 5 in Theorem 1.
Step 7. u n P m = 1 C m ( u 1 ) , as n .
Copy Step 6 in Theorem 1.
Step 8. z n P m = 1 C m ( u 1 ) , v n P m = 1 C m ( u 1 ) and w n P m = 1 C m ( u 1 ) , as n .
Since z n u n = δ n ( λ f ( u n ) F u n ) and δ n 0 , then it is easy to see that z n P m = 1 C m ( u 1 ) , as n .
Since v n = z n + k n ( z n z n 1 ) , then v n P m = 1 C m ( u 1 ) , as n .
Since u n + 1 u n 0 and u n + 1 Q n + 1 C n + 1 , then
w n u n + 1 2 2 α n + β n 2 β n z n u n + 1 2 + 2 γ n 2 β n e n u n + 1 2 + k n 2 2 α n + β n 2 β n z n z n 1 2 2 k n 2 α n + β n 2 β n z n u n + 1 , z n 1 z n z n u n + 1 2 + 2 γ n e n z n + 1 2 + k n 2 z n z n 1 2 + 2 u n + 1 z n z n 1 z n 0 .
Thus, w n P m = 1 C m ( u 1 ) , as n .
Step 9. P m = 1 C m ( u 1 ) i = 1 ( A i + B i ) 1 0 .
In fact, if, otherwise, P m = 1 C m ( u 1 ) ¯ i = 1 ( A i + B i ) 1 0 . Then P m = 1 C m ( u 1 ) i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P m = 1 C m ( u 1 ) .
Since w n = α n v n + β n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) ( v n + w n 2 ) + γ n e n , then β n [ i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) ( v n + w n 2 ) w n ] = α n ( w n v n ) + γ n ( w n e n ) 0 , as n .
Since i n f n β n > 0 , then there exists a subsequence of { n } , which is still denoted by { n } such that i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) ( v n + w n 2 ) w n 0 , as n .
Thus, i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) ( v n + w n 2 ) P m = 1 C m ( u 1 ) , as n .
Since H satisfies Opial’s condition, then
l i m i n f n v n + w n 2 P m = 1 C m ( u 1 ) < l i m i n f n v n + w n 2 i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P m = 1 C m ( u 1 ) l i m i n f n v n + w n 2 i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) v n + w n 2 + l i m i n f n v n + w n 2 P m = 1 C m ( u 1 ) l i m i n f n v n + w n 2 P m = 1 C m ( u 1 ) ,
which makes a contradiction. Thus, P m = 1 C m ( u 1 ) i = 1 ( A i + B i ) 1 0 .
Step 10. P m = 1 C m ( u 1 ) = P i = 1 ( A i + B i ) 1 0 ( u 1 ) .
From Step 9, P i = 1 ( A i + B i ) 1 0 ( u 1 ) u 1 P m = 1 C m ( u 1 ) u 1 .
On the other hand, since i = 1 ( A i + B i ) 1 0 m = 1 C m , then P m = 1 C m ( u 1 ) u 1 P i = 1 ( A i + B i ) 1 0 ( u 1 ) u 1 . Thus,
P i = 1 ( A i + B i ) 1 0 ( u 1 ) u 1 = P m = 1 C m ( u 1 ) u 1 .
Using Lemma 7, we have
P i = 1 ( A i + B i ) 1 0 ( u 1 ) P m = 1 C m ( u 1 ) 2 + P m = 1 C m ( u 1 ) u 1 2 P i = 1 ( A i + B i ) 1 0 ( u 1 ) u 1 2 = P m = 1 C m ( u 1 ) u 1 2 ,
which implies that P i = 1 ( A i + B i ) 1 0 ( u 1 ) = P m = 1 C m ( u 1 ) .  □
Remark 4.
Similar to Remark 2, u n + 1 = P C n + 1 ( u 1 ) is also a possible choice of u n + 1 Q n + 1 in Theorem 2.
Corollary 3.
If δ n 0 , then (16) becomes to the following traditional mid-point inertial forward-backward projection iterative algorithm:
u 0 , u 1 H c h o s e n a r b i t r a r i l y , e 1 H , v n = u n + k n ( u n u n 1 ) , w n = α n v n + β n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) ( v n + w n 2 ) + γ n e n , C 1 = H = Q 1 , C n + 1 = { p C n : w n p 2 2 α n + β n 2 β n u n p 2 + 2 γ n 2 β n e n p 2 + 2 α n + β n 2 β n k n 2 u n u n 1 2 2 k n 2 α n + β n 2 β n u n p , u n 1 u n } , Q n + 1 = { p C n + 1 : u 1 p 2 P C n + 1 ( u 1 ) u 1 2 + σ n + 1 } , u n + 1 Q n + 1 , n N .
If, moreover, k n 0 in (19), then it becomes to the following traditional forward-backward mid-point iterative algorithm:
u 0 , u 1 H c h o s e n a r b i t r a r i l y , e 1 H , w n = α n u n + β n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) ( u n + w n 2 ) + γ n e n , C 1 = H = Q 1 , C n + 1 = { p C n : w n p 2 2 α n + β n 2 β n u n p 2 + 2 γ n 2 β n e n p 2 } , Q n + 1 = { p C n + 1 : u 1 p 2 P C n + 1 ( u 1 ) u 1 2 + σ n + 1 } , u n + 1 Q n + 1 , n N .
Corollary 4.
If k n 0 in (15), then it becomes to the following one:
u 0 , u 1 H c h o s e n a r b i t r a r i l y , e 1 H , z n = δ n λ f ( u n ) + ( I δ n F ) u n , w n = α n z n + β n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) ( w n + z n 2 ) + γ n e n , C 1 = H = Q 1 , C n + 1 = { p C n : w n p 2 2 α n + β n 2 β n z n p 2 + 2 γ n 2 β n e n p 2 } , Q n + 1 = { p C n + 1 : u 1 p 2 P C n + 1 ( u 1 ) u 1 2 + σ n + 1 } , u n + 1 Q n + 1 , n N .
If, moreover, take u n + 1 = P C n + 1 ( u 1 ) in (21), then it becomes to (6).

2.3. Relationship with Variational Inequalities

A lot work has been done on designing iterative algorithms to approximate solution of variational inequalities due to their wide applications (e.g., [26,27]). A classical variational inequality is to find u K such that for v K ,
v u , T u 0 ,
where T : K H is a nonlinear mapping. The symbol V I ( K , T ) denotes the solution of the above variational inequality.

2.3.1. The First Kind Iteration Theorems

Definition 3.
Let H be a real Hilbert space with K being its non-empty closed and convex subset. T : K H is called a τ-Lipschitz continuous mapping if T x T y τ x y , for x , y K .
Theorem 3.
Let H be a real Hilbert space with K being its non-empty closed and convex subset. Suppose A i : K H is μ i -inversely strongly accretive and B i : K H is m-accretive, for each i N . Suppose T : K H is an accretive and τ-Lipschitz continuous mapping. Let { u n } be generated by the following iterative algorithm:
u 0 K , u 1 K c h o s e n a r b i t r a r i l y , e 1 H , y 0 = P K ( u 0 λ 0 T u 0 ) , y n = P K ( u n λ n T u n ) , v n = y n + k n ( y n y n 1 ) , w n = α n v n + β n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P K ( u n λ n T y n ) + γ n e n , C 1 = K = Q 1 , C n + 1 = { p C n : w n p 2 α n y n p 2 + β n u n p 2 + γ n e n p 2 + k n 2 y n y n 1 2 2 α n k n y n p , y n 1 y n } , Q n + 1 = { p C n + 1 : u 1 p 2 P C n + 1 ( u 1 ) u 1 2 + σ n + 1 } , u n + 1 Q n + 1 , n N .
Under the assumptions of (i)–(iv) in Theorem 1 and the additional assumptions ( v ) λ n 0 , ( v i ) i = 1 ( A i + B i ) 1 0 V I ( K , T ) , we have: u n P m = 1 C m ( u 1 ) = P i = 1 ( A i + B i ) 1 0 V I ( K , T ) ( u 1 ) , as n .
Proof. 
We split the proof into nine steps.
Step 1. i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) : K K is non-expansive, for n N .
Copy the proof of Step 1 in Theorem 1.
Step 2. i = 1 ( A i + B i ) 1 0 V I ( K , T ) C n , for n N .
We first show that q i = 1 ( A i + B i ) 1 0 V I ( K , T ) ,
P K ( u n λ n T y n ) q u n q .
In fact, in view of Lemma 7, we have:
P K ( u n λ n T y n ) q 2 u n q λ n T y n 2 u n P K ( u n λ n T y n ) λ n T y n 2 = u n q 2 u n P K ( u n λ n T y n ) 2 + 2 λ n [ T y n T q , q y n + T q , q y n + T y n , y n P K ( u n λ n T y n ) ] u n q 2 u n P K ( u n λ n T y n ) 2 + 2 λ n T y n , y n P K ( u n λ n T y n ) = u n q 2 ( u n y n 2 + y n P K ( u n λ n T y n ) 2 ) + 2 λ n T u n T y n , P K ( u n λ n T y n ) y n + 2 u n y n λ n T u n , P K ( u n λ n T y n ) y n u n q 2 ( u n y n 2 + y n P K ( u n λ n T y n ) 2 ) + 2 λ n τ u n y n y n P K ( u n λ n T y n ) u n q 2 + ( λ n 2 τ 2 1 ) u n y n 2 u n q 2 ,
which implies that (24) is true.
Next, we can easily check the following by noticing the result of Step 1 and (24):
w n q 2 α n v n q 2 + β n u n q 2 + γ n e n q 2 α n y n q 2 + β n u n q 2 + γ n e n q 2 + k n 2 y n y n 1 2 2 α n k n y n q , y n 1 y n .
Thus, by induction as that in Theorem 1, q C n , for n N .
Step 3. C n is a closed and convex subset of H , for each n N .
It is not difficult to see that
w n p 2 α n y n p 2 + β n u n p 2 + γ n e n p 2 + k n 2 y n y n 1 2 2 α n k n y n p , y n 1 y n w n 2 α n y n 2 β n u n 2 γ n e n 2 k n 2 y n y n 1 2 + 2 α n k n y n , y n 1 y n 2 p , w n 2 α n p , y n 2 β n p , u n 2 γ n p , e n + 2 α n k n p , y n 1 y n .
Then C n is a closed and convex subset of H , for each n N .
Copy the results of Steps 4–6 in Theorem 1, we have:
Step 4. Q n is non-empty for each n N , which ensures that { u n } is well-defined.
Step 5. P C n + 1 ( u 1 ) P m = 1 C m ( u 1 ) , as n .
Step 6. u n P m = 1 C m ( u 1 ) , as n .
Step 7. y n P m = 1 C m ( u 1 ) , v n P m = 1 C m ( u 1 ) and w n P m = 1 C m ( u 1 ) , as n .
It is easy to see that Q n is a closed subset of K. Then u n + 1 Q n + 1 K and u n P m = 1 C m ( u 1 ) imply that P m = 1 C m ( u 1 ) K . Therefore,
y n P m = 1 C m ( u 1 ) u n λ n T u n P m = 1 C m ( u 1 ) u n P m = 1 C m ( u 1 ) + λ n T u n 0 ,
as n . Thus, y n P m = 1 C m ( u 1 ) , as n .
Since y n + 1 y n 0 and v n = y n + k n ( y n y n 1 ) , then v n P m = 1 C m ( u 1 ) , as n .
Since u n + 1 Q n + 1 C n + 1 , then from (23), w n u n + 1 2 α n y n u n + 1 2 + β n u n + 1 u n 2 + γ n e n u n + 1 2 + k n 2 y n y n 1 2 2 α n k n y n u n + 1 , y n 1 y n 0 , as n . Thus, w n P m = 1 C m ( u 1 ) , as n .
Step 8. P m = 1 C m ( u 1 ) i = 1 ( A i + B i ) 1 0 V I ( K , T ) .
We shall first show that P m = 1 C m ( u 1 ) V I ( K , T ) .
For this, define
B v = T v + N K v , i f v K , , i f v ¯ K ,
where N K v = { w H : v u , w 0 , u K } is the normal cone to K at v K . It is well-known that B : H H is m-accretive and 0 B v if and only if v V I ( K , T ) [28].
Let z B v = T v + N K v , then z T v N K v . From the definition of the normal cone, we have
v y n , z T v 0 .
From Lemma 7, we have:
u n λ n T u n P K ( u n λ n T u n ) , P K ( u n λ n T u n ) v 0 , v K ,
which implies that
v P K ( u n λ n T u n ) , P K ( u n λ n T u n ) u n λ n + T u n 0 , v K .
In view of (25) and (26), we know that
v y n , z v y n , T v = v P K ( u n λ n T u n ) , T v v P K ( u n λ n T u n ) , T v v P K ( u n λ n T u n ) , P K ( u n λ n T u n ) u n λ n + T u n = v P K ( u n λ n T u n ) , T v T P K ( u n λ n T u n ) + v P K ( u n λ n T u n ) , T P K ( u n λ n T u n ) T u n v P K ( u n λ n T u n ) , P K ( u n λ n T u n ) u n λ n v P K ( u n λ n T u n ) , T P K ( u n λ n T u n ) T u n v P K ( u n λ n T u n ) , P K ( u n λ n T u n ) u n λ n = v y n , T y n T u n v y n , y n u n λ n .
Taking limit on both sides of the above inequality, we have: v P m = 1 C m ( u 1 ) , z 0 , which implies from the fact B is m-accretive that P m = 1 C m ( u 1 ) B 1 0 , and then P m = 1 C m ( u 1 ) V I ( K , T ) .
Next, we shall show that P m = 1 C m ( u 1 ) i = 1 ( A i + B i ) 1 0 .
In fact, if, otherwise, P m = 1 C m ( u 1 ) ¯ i = 1 ( A i + B i ) 1 0 . Then P m = 1 C m ( u 1 ) i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P m = 1 C m ( u 1 ) .
Since w n = α n v n + β n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P K ( u n λ n T y n ) + γ n e n , then β n [ i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P K ( u n λ n T y n ) w n ] = α n ( w n v n ) + γ n ( w n e n ) 0 , as n .
Since i n f n β n > 0 , then there exists a subsequence of { n } , which is still denoted by { n } such that i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P K ( u n λ n T y n ) w n 0 , as n .
Thus, i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P K ( u n λ n T y n ) P m = 1 C m ( u 1 ) , as n . From Step 1 and u n y n 0 , we can also know that i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P K ( u n λ n T u n ) P m = 1 C m ( u 1 ) , as n .
Since H satisfies Opial’s condition, λ n 0 and P m = 1 C m ( u 1 ) K , then
l i m i n f n u n λ n T u n P m = 1 C m ( u 1 ) < l i m i n f n u n λ n T u n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P m = 1 C m ( u 1 ) l i m i n f n u n λ n T u n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P K ( u n λ n T u n ) + l i m i n f n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P K ( u n λ n T u n ) i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) P m = 1 C m ( u 1 ) l i m i n f n u n λ n T u n P m = 1 C m ( u 1 ) ,
which makes a contradiction! Thus, P m = 1 C m ( u 1 ) i = 1 ( A i + B i ) 1 0 .
Step 9. P m = 1 C m ( u 1 ) = P i = 1 ( A i + B i ) 1 0 V I ( K , T ) ( u 1 ) .
From Step 8, P i = 1 ( A i + B i ) 1 0 V I ( K , T ) ( u 1 ) u 1 P m = 1 C m ( u 1 ) u 1 .
On the other hand, since i = 1 ( A i + B i ) 1 0 V I ( K , T ) m = 1 C m , then
P m = 1 C m ( u 1 ) u 1 P i = 1 ( A i + B i ) 1 0 V I ( K , T ) ( u 1 ) u 1 .
Thus
P i = 1 ( A i + B i ) 1 0 V I ( K , T ) ( u 1 ) u 1 = P m = 1 C m ( u 1 ) u 1 .
Using Lemma 7, we have
P i = 1 ( A i + B i ) 1 0 V I ( K , T ) ( u 1 ) P m = 1 C m ( u 1 ) 2 + P m = 1 C m ( u 1 ) u 1 2 P i = 1 ( A i + B i ) 1 0 V I ( K , T ) ( u 1 ) u 1 2 = P m = 1 C m ( u 1 ) u 1 2 ,
which implies that P i = 1 ( A i + B i ) 1 0 V I ( K , T ) ( u 1 ) = P m = 1 C m ( u 1 ) .
This completes the proof. □

2.3.2. The Second Kind Iteration Theorems

The following result is a special result of Lemma 10 in [19]:
Theorem 4.
Let H be a real Hilbert space and K be a non-empty closed and convex subset of H. Suppose f : H H is a contraction with k ( 0 , 1 ) , F : H H is a strongly positive linear bounded operator with coefficient ξ and U : H H is non-expansive mapping. If 0 < λ < ξ 2 k , then there exists x t which satisfies x t = t λ f ( x t ) + ( I t F ) U x t , for 0 < t F 1 . Moreover, x t p 0 as t 0 , and p 0 satisfies the following variational inequality: for z F i x ( U ) ,
( F λ f ) p 0 , p 0 z 0 .
In Lemma 10 of [19], we can also know that the solution of the variational inequality (27) is unique.
Theorem 5.
Under the assumptions of Theorem 2, { u n } generated by (16) converges strongly to P i = 1 ( A i + B i ) 1 0 ( u 1 ) . Set x ˜ = P i = 1 ( A i + B i ) 1 0 ( u 1 ) . If x ˜ = P i = 1 ( A i + B i ) 1 0 [ λ f ( x ˜ ) F ( x ˜ ) + x ˜ ] , then x ˜ satisfies the following variational inequality: z i = 1 ( A i + B i ) 1 0 ,
( F λ f ) x ˜ , x ˜ z 0 .
Proof. 
It follows from Lemma 7 that ( F λ f ) x ˜ , x ˜ z 0 , z i = 1 ( A i + B i ) 1 0 . Since Theorem 4 tells us that (28) has a unique solution, then we know that { u n } generated by (16) converges strongly to the unique solution of variational inequality (28). □
Remark 5.
The assumption that x ˜ = P i = 1 ( A i + B i ) 1 0 [ λ f ( x ˜ ) F ( x ˜ ) + x ˜ ] is reasonable. For example, we may take f ( x ) = x 2 and F ( x ) = λ x 2 , for x H .
Remark 6.
For projection iterative algorithms such as (16), rare work can be found to show that the limit of the iterative sequences is also the solution of a kind of variational inequalities.

3. Applications

3.1. Preparation for Discussion of Capillarity Systems

To present some examples in this section, we need some basic definitions in Banach spaces.
Let E be a real Banach space with E * being its dual space and let · , · denote the generalized duality pairing between E and E * .
Definition 4.
(see [29]) Recall that J : E 2 E * is called the normalized duality mapping if x E ,
J x = { f E * : x , f = x 2 = f 2 } .
Definition 5.
(see [29]) Recall that A : E E * is said to be a monotone mapping if x i D ( A ) , i = 1 , 2 , one has
x 1 x 2 , A x 1 A x 2 0 .
A monotone mapping A : E E * is said to be maximal monotone if R ( J + r A ) = E * , r > 0 .
Definition 6.
(see [29]) Recall that a mapping B : E E * is said to be coercive if { x n } D ( B ) with lim n x n = + , then lim n x n , B x n x n = + .
Definition 7.
(see [29]) Recall that a mapping B : D ( B ) = E E * is said to be a hemi-continuous mapping if B ( x + t y ) B x , as t 0 , for x , y E .
Definition 8.
(see [29]) ψ : E ( , + ] is said to be a proper convex functional if there exists u 0 E such that ψ ( u 0 ) < + and ψ ( ( 1 λ ) u + λ v ) ( 1 λ ) ψ ( u ) + λ ψ ( v ) , for u , v E and λ [ 0 , 1 ] . ψ : E ( , + ] is said to be lower-semi-continuous: lim inf y x ψ ( y ) ψ ( x ) , for x E . The subdifferential of ψ, ψ : E E * , is defined by:
ψ ( u ) = { w E * : ψ ( u ) ψ ( v ) u v , w , v E } , u E .

3.2. Applications to Capillarity Elliptic Systems

Example 1.
Suppose Ω is bounded conical domain in R n ( n N ) with Γ C 1 , ϑ is the exterior normal derivative of Γ, ε i is a non-negative constant, λ i is a positive number, f i ( x ) L p i ( Ω ) is a given function, for i = 1 , 2 , , M . In addition, β x : R R is the subdifferential of φ x , where φ x = φ ( x , · ) : R R , for each x Γ . Suppose 2 n n + 1 < p i < + , i = 1 , 2 , , M . If p i n then 1 q i , r i < + and if p i < n then 1 q i , r i n p i n p i for i = 1 , 2 , , M .
The following capillarity system is studied in [30]:
d i v [ ( 1 + | u ( i ) | p i 1 + | u ( i ) | 2 p i ) | u ( i ) | p i 2 u ( i ) ] + λ i ( | u ( i ) | q i 2 u ( i ) + | u ( i ) | r i 2 u ( i ) ) + ε i g i ( x , u ( i ) , u ( i ) ) = f i ( x ) , x Ω < ϑ , ( 1 + | u ( i ) | p i 1 + | u ( i ) | 2 p i ) | u ( i ) | p i 2 u ( i ) > β x ( u ( i ) ( x ) ) , x Γ , i = 1 , 2 , , M ,
where | · | and < · , · > denote the norm and inner product in R n , respectively.
The study on (29) in [30] is based on the following assumptions.
(1)
x Γ , φ x = φ ( x , · ) : R R is a proper convex and lower-semi-continuous mapping with φ x ( 0 ) = 0 .
(2)
0 β x ( 0 ) , t R , x Γ ( I + λ β x ) 1 ( t ) R is measurable for λ > 0 .
(3)
For each i { 1 , 2 , , M } , g i : Ω × R n × R R satisfies Caratheodory’s conditions and satisfies that
| g i ( x , r 1 , r 2 , , r n + 1 ) g i ( x , s 1 , s 2 , , s n + 1 ) | b i | r n + 1 s n + 1 | ,
for ( r 1 , r 2 , , r n + 1 ) , ( s 1 , s 2 , , s n + 1 ) R n + 1 .
By using splitting method, the sufficient condition that (29) has non-trivial solution is obtained:
Theorem 6.
(see [30]) If u ( i ) L p i ( Ω ) satisfies that
Ω [ g i ( x , u ( i ) , u ( i ) ) f i ] | u ( i ) | p i 2 u ( i ) d x 0 ,
for i = 1 , 2 , , M , then u = ( u ( 1 ) , u ( 2 ) , , u ( M ) ) is the non-trivial solution of capillarity system (29).
Based on Example 1, we present the following example:
Example 2.
Suppose Ω, Γ and ϑ are the same as those in Example 1. Suppose λ i > 0 , 2 n n + 1 < p i < + . If p i n , then suppose 1 q i , r i < + and if p i < n , then suppose 1 q i , r i n p i n p i , for i N .
Now, we will discuss the following capillarity systems.
d i v [ ( 1 + | u ( i ) | p i 1 + | u ( i ) | 2 p i ) | u ( i ) | p i 2 u ( i ) ] + λ i ( | u ( i ) | q i 2 u ( i ) + | u ( i ) | r i 2 u ( i ) ) + u ( i ) ( x ) = f i ( x ) , x Ω < ϑ , ( 1 + | u ( i ) | p i 1 + | u ( i ) | 2 p i ) | u ( i ) | p i 2 u ( i ) > = 0 , x Γ , i N .
Please note that (30) is the extension from the finite case of (29) to that for infinite case. However, both the capillarity equations and the boundary conditions are the special case of (29) in the sense that ε i 1 and g i ( x , u ( i ) , u ( i ) ) u ( i ) for i N and β x 0 , for x Γ .
Lemma 9.
(see [30]) The mapping U i : W 1 , p i ( Ω ) ( W 1 , p i ( Ω ) ) * defined by
w , U i u = Ω < ( 1 + | u | p i 1 + | u | 2 p i ) | u | p i 2 u ,
v > d x + λ i Ω | u ( x ) | q i 2 u ( x ) v ( x ) d x + λ i Ω | u ( x ) | r i 2 u ( x ) v ( x ) d x ,
for u , w W 1 , p i ( Ω ) , is everywhere defined, hemi-continuous, monotone and coercive, for each i N .
Lemma 10.
(see [30]) Define B i : L 2 ( Ω ) L 2 ( Ω ) by
D ( B i ) = { u L 2 ( Ω ) | f L 2 ( Ω ) s u c h t h a t f U i u } .
For u D ( B i ) , B i u = { f L 2 ( Ω ) | f U i u } . Then B i : L 2 ( Ω ) L 2 ( Ω ) is m-accretive, for each i N .
Lemma 11.
(see [30]) The mapping A i : D ( A i ) L 2 ( Ω ) L 2 ( Ω ) defined by
( A i u ) ( x ) = u ( x ) f i ( x ) , u ( x ) D ( A i ) ,
is μ i -inversely strongly accretive, for μ i ( 0 , 1 ] , for i N .
Theorem 7.
If f i ( x ) λ i ( | k | q i 1 + | k | r i 1 ) s g n k + k , then { u ( i ) k : i N } is the solution of capillarity system (30). Moreover, { k } = i = 1 ( A i + B i ) 1 0 .
Proof. 
It is easy to see that { u ( i ) k : i N } is the solution of capillarity system (30) and { k } i = 1 ( A i + B i ) 1 0 . Now, we shall show that { k } i = 1 ( A i + B i ) 1 0 .
In fact, if A i u + B i u = 0 and A i v + B i v = 0 , then u + B i u = v + B i v , which implies that
0 u v , B i u B i v = u v , v u 0 .
Thus, u = v and then i = 1 ( A i + B i ) 1 0 is a singleton. Since k i = 1 ( A i + B i ) 1 0 , then the result follows. □
Theorem 8.
Let f i ( x ) λ i ( | k | q i 1 + | k | r i 1 ) s g n k + k , for i N . Suppose A i and B i are the same as those in Lemmas 10 and 11, respectively. Let F : L 2 ( Ω ) L 2 ( Ω ) be any strongly positive linear bounded operator with coefficient ξ > 0 and f : L 2 ( Ω ) L 2 ( Ω ) be a contraction with coefficient k ( 0 , 1 ) . Constructing the following iterative algorithm:
u 0 , u 1 L 2 ( Ω ) c h o s e n a r b i t r a r i l y , e 1 L 2 ( Ω ) , z 0 = u 0 , z n = δ n λ f ( u n ) + ( I δ n F ) u n , v n = z n + k n ( z n z n 1 ) , w n = α n v n + β n i = 1 ω n , i ( I + s n , i B i ) 1 ( I s n , i A i ) ( v n + w n 2 ) + γ n e n , C 1 = L 2 ( Ω ) = Q 1 , C n + 1 = { p C n : w n p 2 2 α n + β n 2 β n z n p 2 + 2 γ n 2 β n e n p 2 + 2 α n + β n 2 β n k n 2 z n z n 1 2 2 k n 2 α n + β n 2 β n z n p , z n 1 z n } , Q n + 1 = { p C n + 1 : u 1 p 2 P C n + 1 ( u 1 ) u 1 2 + σ n + 1 } , u n + 1 Q n + 1 , n N .
Under the assumptions of Theorem 2, using the result of Theorem 5, one has u n ( x ) P i = 1 ( A i + B i ) 1 0 ( u 1 ) , which is the unique solution of capillarity system (30) and satisfies the following variational inequality: For z ( x ) i = 1 ( A i + B i ) 1 0 ,
( F λ f ) P i = 1 ( A i + B i ) 1 0 ( u 1 ) , P i = 1 ( A i + B i ) 1 0 ( u 1 ) z ) 0 .
Remark 7.
From Theorem 8 we can easily see the relationship among the solution of capillarity system, the solution of variational inequality and the zero of sum of infinitely many m-accretive mappings and infinitely many μ i -inversely strongly accretive mappings.

Author Contributions

L.W. and R.P.A. contributed are equally to Section 1 and Section 2. The Y.S. contributed to Section 3.

Funding

Supported by the National Natural Science Foundation of China (11071053), Natural Science Foundation of Hebei Province (A2014207010), Key Project of Science and Research of Hebei Educational Department (ZD2019073), Key Project of Science and Research of Hebei University of Economics and Business (2018ZD06), Youth Project of Science and Research of Hebei University of Economics and Business (2017KYQ09) and Youth Project of Science and Research of Hebei Educational Department (QN2017328).

Acknowledgments

Thanks for the reviewers’ valuable opinions and careful work.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Opial, Z. Weak convergence of the sequence of successive approximations for nonexpansive mappings. Bull. Am. Math. Soc. 1967, 73, 591–597. [Google Scholar] [CrossRef] [Green Version]
  2. Takahashi, W. Nonlinear Functional Analysis. Fixed Point Theory and its Applications; Yokohama Publishers: Yokohama, Japan, 2000. [Google Scholar]
  3. Wei, L.; Shi, A.F. Splitting-midpoint method for zeros of the sum of accretive operator and μ-inversely strongly accretive operator in a q-uniformly smooth Banach space and its applications. J. Inequal. Appl. 2015, 2015, 183. [Google Scholar] [CrossRef]
  4. Shehu, Y. Iterative approximations for zeros of the sum of accretive operators in Banach spaces. J. Funct. Spaces 2016, 2016, 5973468. [Google Scholar] [CrossRef]
  5. Combettes, P.L.; Wajs, V.R. Signal recovery by proximal forward-backward splitting. Multiscale Model. Simul. 2005, 4, 1168–1200. [Google Scholar] [CrossRef]
  6. Tseng, P. A mordified forward-backward splitting method for maximal monotone mappings. SIAM J. Control Optim. 2000, 38, 431–446. [Google Scholar] [CrossRef]
  7. Khan, S.A.; Suantai, S.; Cholamjiak, W. Shrinking projection methods involving inertial forward-backward splitting methods for inclusion problems. RACSAM 2018, 113, 645–656. [Google Scholar] [CrossRef]
  8. Moudafi, A.; Oliny, M. Convergence of a splitting inertial proximal method for monotone operators. J. Comput. Appl. Math. 2003, 155, 447–454. [Google Scholar] [CrossRef] [Green Version]
  9. Lorenz, D.; Pock, T. An inertial forward-backward algorithm for monotone inclusions. J. Math. Imaging Vis. 2015, 51, 311–325. [Google Scholar] [CrossRef]
  10. Cholamjiak, P. A generalized forward-backward splitting method for solving quasi inclusion problems in Banach spaces. Numer. Algorithm 2016, 71, 915–932. [Google Scholar] [CrossRef]
  11. Lopez, G.; Martin-Marquez, V.; Wang, F.; Xu, H.K. Forward-backward splitting method for accretive operators in Banach spaces. Abstr. Appl. Anal. 2012, 2012, 109236. [Google Scholar] [CrossRef]
  12. Qin, X.; Cho, S.Y.; Wang, L. Convergence of splitting algorithms for the sum of two accretive operators with applications. Fixed Point Theory Appl. 2014, 2014, 166. [Google Scholar] [CrossRef] [Green Version]
  13. Bot, R.I.; Csetnek, E.R.; Hendrich, C. Inertial Douglas-Rachford splitting for monotone inclusion. Appl. Math. Comput. 2015, 256, 472–487. [Google Scholar]
  14. Dong, Q.L.; Jiang, D.; Cholamjiak, P.; Shehu, Y. A strong convergence result involving an inertial forward-backward algorithm for monotone inlusions. J. Fixed Point Theory Appl. 2017, 19, 3097–3118. [Google Scholar] [CrossRef]
  15. Wei, L.; Duan, L.L. A new iterative algorithm for the sum of two different types of finitely many accretive operators in Banach and its connection with capillarity equation. Fixed Point Theory Appl. 2015, 2015, 25. [Google Scholar] [CrossRef]
  16. Wei, L.; Duan, L.L.; Agarwal, R.P.; Chen, R.; Zheng, Y.Q. Mordified forward-backward splitting midpoint method with superposition perturbations for sum of two kind of infinite accretive mappings and its applications. J. Inequal. Appl. 2017, 2017, 227. [Google Scholar] [CrossRef] [PubMed]
  17. Wei, L.; Shen, Y.W.; Zheng, Y.Q.; Tan, R.L. A new iterative scheme for the sum of infinite m-accretive mappings and inversely strongly accretive mappings and its application. J. Nonlinear Var. Appl. 2017, 1, 345–356. [Google Scholar]
  18. Wei, L.; Agarwal, R.P. Simple form of a projection set in hybrid iterative schemes for non-linear mappings, application of inequalities and computational experiments. J. Inequal. Anal. 2018, 2018, 179. [Google Scholar] [CrossRef] [PubMed]
  19. Wei, L.; Agarwal, R.P. A new iterative algorithm for the sum of infinite m-accretive mappings and infinite μi-inversely strongly accretive mappings and its applications to integro-differentail systems. Fixed Point Theory Appl. 2016, 2016, 7. [Google Scholar] [CrossRef]
  20. Mosco, U. Convergence of convex sets and of solutions of variational inequalities. Adv. Math. 1969, 3, 510–585. [Google Scholar] [CrossRef] [Green Version]
  21. Tsukada, M. Convergence of best approximations in a smooth Banach space. J. Approx. Theory 1984, 40, 301–309. [Google Scholar] [CrossRef] [Green Version]
  22. Bruck, R.E. Properties of fixed-point sets of nonexpansive mappings in Banach spaces. Trans. Am. Math Soc. 1973, 179, 251–262. [Google Scholar] [CrossRef]
  23. Xu, H.K. Inequalities in Banach space with applications. Nonlinear Anal. 1991, 16, 1127–1138. [Google Scholar] [CrossRef]
  24. Nakajo, K.; Takahashi, W. Strong convergence theorems for nonexpansive mappings and nonexpansive semigroups. J. Math. Anal. Appl. 2003, 279, 327–379. [Google Scholar] [CrossRef]
  25. Agarwal, R.P.; O’Regan, D.; Sahu, D.R. Fixed Point Theory for Lipschtz-type Mappings with Applications; Springer: New York, NY, USA, 2009. [Google Scholar]
  26. He, S.N.; Liu, L.L.; Gibali, A. Self-adaptive iterative method for solving boundedly Lipschitz continuous and strongly monotone variational inequalities. J. Inequal. Appl. 2018, 2018, 350. [Google Scholar] [CrossRef] [PubMed]
  27. Gibali, A.; Reich, S.; Zalas, R. Ooter approximation methods for solving variational inequalities in Hilbert space. Optimization 2017, 66, 417–437. [Google Scholar] [CrossRef]
  28. Rockafellar, R.T. On the maximality of sums of nonlinear monotone operators. Trans. Am. Math. Soc. 1970, 149, 75–88. [Google Scholar] [CrossRef]
  29. Pascali, D.; Sburlan, S. Nonlinear Mappings of Monotone Type; Sijthoff Noordhoff International Publishers: Bucharest, Romania, 1978. [Google Scholar]
  30. Wei, L.; Chen, R. Study on the existence of non-trivial solution of one kind capillarity systems. J. Math. (PRC) 2017, 37, 390–400. (In Chinese) [Google Scholar]

Share and Cite

MDPI and ACS Style

Wei, L.; Shang, Y.; Agarwal, R.P. New Inertial Forward-Backward Mid-Point Methods for Sum of Infinitely Many Accretive Mappings, Variational Inequalities, and Applications. Mathematics 2019, 7, 466. https://doi.org/10.3390/math7050466

AMA Style

Wei L, Shang Y, Agarwal RP. New Inertial Forward-Backward Mid-Point Methods for Sum of Infinitely Many Accretive Mappings, Variational Inequalities, and Applications. Mathematics. 2019; 7(5):466. https://doi.org/10.3390/math7050466

Chicago/Turabian Style

Wei, Li, Yingzi Shang, and Ravi P. Agarwal. 2019. "New Inertial Forward-Backward Mid-Point Methods for Sum of Infinitely Many Accretive Mappings, Variational Inequalities, and Applications" Mathematics 7, no. 5: 466. https://doi.org/10.3390/math7050466

APA Style

Wei, L., Shang, Y., & Agarwal, R. P. (2019). New Inertial Forward-Backward Mid-Point Methods for Sum of Infinitely Many Accretive Mappings, Variational Inequalities, and Applications. Mathematics, 7(5), 466. https://doi.org/10.3390/math7050466

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop