Loading [MathJax]/jax/output/SVG/jax.js
Research article

Asymptotic for a second order evolution equation with damping and regularizing terms

  • Received: 05 December 2020 Accepted: 18 February 2021 Published: 01 March 2021
  • MSC : 34Exx, 46N10

  • Let H be a real Hilbert space. We investigate the long time behavior of the trajectories x(.) of the vanishing damped nonlinear dynamical system with regularizing term

    x(t)+γ(t)x(t)+Φ(x(t))+ε(t)U(x(t))=0,(GAVDγ,ε)

    where Φ,U:HR are two convex continuously differentiable functions, ε(.) is a decreasing function satisfying limt+ε(t)=0, and γ(.) is a nonnegative function which behaves, for t large enough, like Ktθ where K>0 and 0θ1. The main contribution of this paper is the following control result: If +0ε(t)γ(t)dt=+, U is strongly convex and its unique minimizer x is also a minimizer of Φ then every trajectory x(.) of (GAVDγ,ε) converges strongly to x and the rate of convergence to 0 of its energy function

    W(t)=12x(t)2+Φ(x(t))Φ+ε(t)(U(x(t))U)

    is of order to (1/t1+θ). Moreover, we prove a new result concerning the weak convergence of the trajectories of (GAVDγ,ε) to a common minimizer of Φ and U (if one exists) under a simple condition on the speed of decay of the regularizing factor ε(t) to 0.

    Citation: Ramzi May, Chokri Mnasri, Mounir Elloumi. Asymptotic for a second order evolution equation with damping and regularizing terms[J]. AIMS Mathematics, 2021, 6(5): 4901-4914. doi: 10.3934/math.2021287

    Related Papers:

    [1] Ziran Yin, Chongyang Liu, Xiaoyu Chen, Jihong Zhang, Jinlong Yuan . A comprehensive characterization of the robust isolated calmness of Ky Fan k-norm regularized convex matrix optimization problems. AIMS Mathematics, 2025, 10(3): 4955-4969. doi: 10.3934/math.2025227
    [2] Aníbal Coronel, Fernando Huancas . New results for the non-oscillatory asymptotic behavior of high order differential equations of Poincaré type. AIMS Mathematics, 2022, 7(4): 6420-6444. doi: 10.3934/math.2022358
    [3] Woocheol Choi, Young-Pil Choi . A sharp error analysis for the DG method of optimal control problems. AIMS Mathematics, 2022, 7(5): 9117-9155. doi: 10.3934/math.2022506
    [4] Adel M. Al-Mahdi, Maher Noor, Mohammed M. Al-Gharabli, Baowei Feng, Abdelaziz Soufyane . Stability analysis for a Rao-Nakra sandwich beam equation with time-varying weights and frictional dampings. AIMS Mathematics, 2024, 9(5): 12570-12587. doi: 10.3934/math.2024615
    [5] M. Ali Akbar, Norhashidah Hj. Mohd. Ali, M. Tarikul Islam . Multiple closed form solutions to some fractional order nonlinear evolution equations in physics and plasma physics. AIMS Mathematics, 2019, 4(3): 397-411. doi: 10.3934/math.2019.3.397
    [6] Jianbo Yuan, Shixuan Zhang, Yongqin Xie, Jiangwei Zhang . Exponential attractors for nonclassical diffusion equations with arbitrary polynomial growth nonlinearity. AIMS Mathematics, 2021, 6(11): 11778-11795. doi: 10.3934/math.2021684
    [7] M. TarikulIslam, M. AliAkbar, M. Abul Kalam Azad . Traveling wave solutions in closed form for some nonlinear fractional evolution equations related to conformable fractional derivative. AIMS Mathematics, 2018, 3(4): 625-646. doi: 10.3934/Math.2018.4.625
    [8] Tingting Liu, Tasneem Mustafa Hussain Sharfi, Qiaozhen Ma . Time-dependent asymptotic behavior of the solution for evolution equation with linear memory. AIMS Mathematics, 2023, 8(7): 16208-16227. doi: 10.3934/math.2023829
    [9] Huiyan Peng, Xuemei Wei . Qualitative analysis of a time-delayed free boundary problem for tumor growth with Gibbs-Thomson relation in the presence of inhibitors. AIMS Mathematics, 2023, 8(9): 22354-22370. doi: 10.3934/math.20231140
    [10] Zayd Hajjej . On the exponential decay of a Balakrishnan-Taylor plate with strong damping. AIMS Mathematics, 2024, 9(6): 14026-14042. doi: 10.3934/math.2024682
  • Let H be a real Hilbert space. We investigate the long time behavior of the trajectories x(.) of the vanishing damped nonlinear dynamical system with regularizing term

    x(t)+γ(t)x(t)+Φ(x(t))+ε(t)U(x(t))=0,(GAVDγ,ε)

    where Φ,U:HR are two convex continuously differentiable functions, ε(.) is a decreasing function satisfying limt+ε(t)=0, and γ(.) is a nonnegative function which behaves, for t large enough, like Ktθ where K>0 and 0θ1. The main contribution of this paper is the following control result: If +0ε(t)γ(t)dt=+, U is strongly convex and its unique minimizer x is also a minimizer of Φ then every trajectory x(.) of (GAVDγ,ε) converges strongly to x and the rate of convergence to 0 of its energy function

    W(t)=12x(t)2+Φ(x(t))Φ+ε(t)(U(x(t))U)

    is of order to (1/t1+θ). Moreover, we prove a new result concerning the weak convergence of the trajectories of (GAVDγ,ε) to a common minimizer of Φ and U (if one exists) under a simple condition on the speed of decay of the regularizing factor ε(t) to 0.



    Let H be a real Hilbert space endowed with the inner product .,. and the associated norm .. Let Φ,U:HR be two convex continuously differentiable functions and γ,ε two real positive functions defined on a fixed time interval [t0,+) for some t0>0. Motivated by the work [3] of Attouch, Chbani, and Riahi on the asymptotic behavior of the trajectories of the asymptotic vanishing damping dynamical system with regularizing regularizing term

    x(t)+αtx(t)+Φ(x(t))+ε(t)x(t)=0,(AVDγ,ε)

    we investigate in this paper the long time behavior, as t+, of the trajectories of the following generalized version of the (AVDα,ε) dynamical system

    x(t)+γ(t)x(t)+Φ(x(t))+ε(t)U(x(t))=0.(GAVDγ,ε)

    For the importance and the applications of these two dynamical systems and many other related dynamical systems in Mechanics and Optimization, we refer the reader to [2,5,6,14] and references therein.

    Throughout this paper, we assume the following general hypothesis:

    (H1) The functions Φ,U:HR are convex, differentiable, and bounded from below. We set Φ=infxHΦ(x) and U=infxHU(x).

    (H2) The set SΦ:=argminΦ={zH:Φ(z)=Φ} is nonempty.

    (H3) The gradient functions Φ and U of Φ and U are Lipschitz on bounded subsets of H.

    (H4) The function γ:[t0,+)(0,+) is absolutely continuous and satisfies the following property: there exist t1t0 and two real constants K1,K2>0 such that

    γ(t)K1t and γ(t)K2t2

    for almost every tt1.

    (H5) The function ε:[t0,+)(0,+) is absolutely continuous, nonincreasing and satisfies

    limt+ε(t)=0.

    Proceeding as in the proof of [Theorem 3.1, [6]] and using the classical Cauchy-Lipschitz theorem and the energy function

    W(t)=12x(t)2+Φ(x(t))Φ+ε(t)(U(x(t))U), (1.1)

    one can easily prove that for every initial data (x0,v0)H×H, the dynamical system (GAVDγ,ε) has a unique solution x(.)C2(t0,+;H) which satisfies x(t0)=x0 and x(t0)=v0. Therefore, we assume in what follows that x(.) is a global solution of (GAVDγ,ε) and focus our attention on the study of the long time behavior of x(t) as t goes to infinity. Before starting the presentation of the main contributions of this work in this direction, let us first recall some well known results on the asymptotic behavior of solutions of a variant dynamical systems related to (GAVDγ,ε). In the pioneer work [1], Alvarez considered the case where γ(.) is constant and ε=0. He established that the trajectory x(t) converges weakly to some element ˉx of SΦ. He also proved that the rate of convergence of Φ(x(t)) to Φ is of order (1/t) (see [2]). To overcome the drawback of the weak convergence to a non identified minimizer of Φ, Attouch and Cazerniki [5] proved that adding a regularizing term ε(t)x(t) forces any trajectory x(t) of the system

    x(t)+γx(t)+Φ(x(t))+ε(t)x(t)=0 (1.2)

    to converge strongly to the element x of minimum norm of the set SΦ provided that +t0ε(t)dt=+. Using a different approach, Jendoubi and May [10] proved that this strong convergence result remains true even if a perturbation integrable term g(t) is added to the equation (1.2). In an other direction, in order to improve the rate of convergence of Φ(x(t)) to Φ, Su, Boyd, and Candes [14] introduced the following dynamical system which is the continuous version of the Nestrov's accelerated minimization method [12]:

    x(t)+αtx(t)+Φ(x(t))=0. (1.3)

    They proved that if α3 then

    Φ(x(t))Φ=O(1/t2).

    This result was later improved in [4] and [11]. In fact it was proved that if α>3 then x(t) converges weakly to some element ˉx of SΦ and that

    Φ(x(t))Φ=(1/t2).

    In order to benefit at the same time of the quick minimization property Φ(x(t))Φ=(1/t2) due to the presence of the vanishing damping term γ(t)=αt in (1.3) and the strong convergence of the trajectories of (1.2) to a particular minimizer of Φ which is a consequence of the regularizing term ε(t)x(t), Attouch, Chbani, and Riahi [3] have considered the dynamical system (AVDα,ε) and have established some properties of the asymptotic behavior of its trajectories which are summarized in the following theorem.

    Theorem 1.1 (Attouch, Chbani and Riahi). Let xC2(t0,+;H) be a solution of (AVDα,ε). The following assertions hold:

    (A) If α>1 and +t0ε(t)tdt<+, then +t0x(t)2tdt<+, limt+x(t)=0 and limt+Φ(x(t))=Φ.

    (B) If α>3 and +t0tε(t)dt<+, then x(t) converges weakly to some element of SΦ.

    Furthermore, the associated energy function E(t)=12x(t)2+Φ(x(t))Φ satisfies E(t)=(1/t2) and +t0tE(t)dt<+.

    (C) If the function ε satisfies moreover one of the following hypothesis

    (H5a) limt+t2ε(t)=+ if α=3

    (H5b) t2ε(t)c>49α(α3) if α>3

    (H5c) +t0ε(t)tdt=+

    then lim inft+x(t)x=0 where x is the element of minimal norm of the set SΦ.

    In our present paper, we improve and extend these results to the general dynamical system (GAVDγ,ε). Moreover, we highlight some new asymptotic properties of the trajectories of (GAVDγ,ε).

    Our first main result is a general minimization property of (GAVDγ,ε) which improves the assertion (A) in the previous theorem.

    Theorem 1.2 (A general minimization property of (GAVDγ,ε)). Let x(.) be a solution of (GAVDγ,ε). Then +t0γ(t)x(t)2dt<+, and the energy function W(t), defined by (1.1) decreases and converges to 0 as t+. In particular limt+x(t)=0 and limt+Φ(x(t))=Φ.

    The second result concerns the weak convergence properties of the trajectories of (GAVDγ,ε). The first part of this result is similar to the assertion (B) in Theorem 1.1. Our proof, which is different from the arguments given by Attouch, Chbani, and Riahi [Theorem 3.1, [3]], provides an other confirmation of the fact, noticed recently in many works as [2,3,11] and [14], that the value α=3 in the the system (1.3) is critical and somehow mysterious. The second part of the theorem is a simple result on the weak convergence to a common minimizer of the two convex functions Φ and U. At our knowledge, this result was not known even in the case where the damping term γ is constant. A comparable result has been proved by Cabot(see [Proposition 2.5, [7]]) for the first order system x(t)+Φ(x(t))+ε(t)U(x(t))=0.

    Theorem 1.3 (Weak convergence properties of (GAVDγ,ε)). Assume that there exist t1t0,0θ1,α>0 with α>3 if θ=1 such that

    γ(t)αtθtt1and+t0[(tθγ(t))]+dt<+ (1.4)

    where [(tθγ(t))]+=max{0,(tθγ(t))}. Let x(.) be a solution of (GAVDγ,ε). Then the two following properties hold:

    (P1) If +t0tθε(t)dt<+ then x(t) converges weakly to some element of SΦ.

    (P2) If SΦSU and lim inft+t1+θε(t)>0then x(t) converges weakly to some element of SΦSU.

    Moreover, in both case, the energy function W satisfies the following asymptotical behavior

    W(t)=(1/t1+θ)and+t0tθW(t)dt<+. (1.5)

    The last result deals with the strong convergence of the trajectories of (GAVDγ,ε) to a minimizer of the function U on the set of minimizers of Φ.

    Theorem 1.4 (Strong convergence properties of (GAVDγ,ε)). Assume that U is strongly convex and γ(t)=αtθ such that α>0 and 0θ<1 or α>3 if θ=1. Suppose in addition that +t0tθε(t)dt=+. Let x(.) be a solution of (GAVDγ,ε). Then the two following assertions hold true:

    (Q1) If x(t)=(1/tθ) and +t0tθx(t)2dt<+ then x(t) converges strongly to the unique minimizer p of U on SΦ.

    (Q2) If the unique minimizer x of U on H belongs to SΦ then x(t) converges strongly to x and the energy function W satisfies the asymptotic properties (1.5).

    Combining Theorem 1.2 and Theorem 1.4 provides a new proof of following important result due to Attouch and Cazernicki (see [Theorem 2.3, [5]]).

    Theorem 1.5 (Attouch and Cazernicki). Let α>0. If +t0ε(t)dt=+, then any trajectory x(.) to the dynamical system

    x(t)+αx(t)+Φ(x(t))+ε(t)x(t)=0, (1.6)

    converges strongly to the projection of zero on the closed and convex subset SΦ.

    This section is devoted to the proof of Theorem 1.2. The main idea of the proof is inspired by the paper [9].

    Proof. Differentiating the energy function W defined by (1.1) and using the equation (GAVDγ,ε), we obtain

    W(t)=γ(t)x(t)2+ε(t)(U(x(t))U)γ(t)x(t)2. (2.1)

    Hence

    +t0γ(t)x(t)2dt<, (2.2)

    and the function W(t) decreases and converges to some positive real number W as t+. Therefore, to conclude, we just have to show that W0. Let v be an arbitrarily element of H. Consider the function

    hv(t)12x(t)v2.

    Using the equation (GAVDγ,ε) and the convexity of Φ and U, one can easily check that

    hv(t)+γ(t)hv(t)=x(t)2+Φ(x(t)),vx(t)+ε(t)U(x(t)),vx(t)x(t)2+Φ(v)Φ(x(t))+ε(t)(U(v)U(x(t)))=32x(t)2W(t)+Φ(v)Φ+ε(t)(U(v)U). (2.3)

    Recalling that W(t)W, we get

    Ahv(t)γ(t)hv(t)+32x(t)2+ε(t)(U(v)U)

    where A=W+ΦΦ(v).

    Integrating the last inequality from t0 to t>t0 and using the fact that γhv0 and the assumption γ(t)K2t2, we find

    (tt0)Ahv(t0)+γ(t0)hv(t0)hv(t)+32tt0x(s)2ds+tt0fv(s)ds, (2.4)

    where

    fv(s)=ε(s)(U(v)U)+K2s2hv(s). (2.5)

    Since γ(t)K1t, we deduce from (2.2) that +t0x(s)2sds< which, thanks to [Lemma 3.2, [9]], implies that

    tt0x(s)2ds=(t). (2.6)

    Using now the Cauchy-Schwarz inequality, we infer

    x(t)x(t0)+tt0(tt0x(s)2ds)12=(t), (2.7)

    which implies that limt+fv(t)=0. Hence, we deduce that

    tt0fv(s)ds=(t). (2.8)

    Recalling that since W is bounded, x is also bounded. Thus, from (2.7), we infer that

    hv(t)=2x(t),x(t)v=(t). (2.9)

    Finally, dividing the inequality (2.4) by t, using the estimates (2.6), (2.8), (2.9) and letting t+, we obtain A0, which implies that WΦ(v)Φ. Since this holds for every vH, the required result W0 follows.

    Remark 2.1. Let us notice that in the proof of Theorem 1.2, we did not use the hypothesis (H2). Moreover, we can prove that if SΦ is empty, then any solution x(.) of the (GAVDγ,ε) system satisfies x(t)+ as t+. Indeed, otherwise there exists a sequence (tn)n tending to + so that (x(tn))n converges weakly to an element ˉxH. From the lower semi-continuity property it follows that

    Φ(ˉx)lim infn+Φ(x(tn)),

    which, thanks to the fact limt+Φ(x(t))=Φ, implies that Φ(ˉx)Φ and contradicts the assumption SΦ=.

    In this section, we prove Theorem 1.3. The proof relies on the classical Opial's lemma and the following technical lemma which will be also useful in the study of the strong convergence properties of the trajectories of (GAVDγ,ε) in the next section.

    Lemma 3.1. Assume that the function γ(.) satisfies the assumption (1.4) in Theorem 1.3. Let x(.) be a solution of (GAVDγ,ε) and let vSΦ such that the positive function tθrv(t) belongs to L1(t0,+;R) where rv(t)=ε(t)(U(v)U). Then the function hv(t)=12x(t)v2 converges as t+ and the energy function W satisfies the asymptotic property (1.5).

    Proof. First, we notice that up to take t1 large enough we can assume that

    γ(t)Kt for every tt1

    with K>3 and K=α if θ=1.

    Let λ(t)=t1+θ. Using (2.1) and the above inequality, we find

    (λW)λWλγx2λWK1+θλx2λWK2λx2. (3.1)

    Therefore,

    32λx23KλW3K(λW).

    Multiplying (2.3) by λ(t) (we recall that, since vSΦ, Φ(v)=Φ) and using the above inequality, we obtain

    (13K)λW+3K(λW)λhvλγhv+λrv.

    Integrating this inequality from t1 to t>t1 leads after some simple computations to the following inequality

    (13K)tt1λ(s)W(s)ds+3Kλ(t)W(t)C0λ(t)hv(t)+fθ(t)hv(t)+tt1gθ(s)hv(s)ds, (3.2)

    where

    fθ(t)=λ(t)(λγ)(t), (3.3)
    gθ(t)=[(λγ)]+(t)λ(t), (3.4)

    and

    C0=λ(t1)hv(t1)λ(t1)hv(t1)+3Kλ(t1)W(t1)+λ(t1)hv(t1)++t1λ(s)rv(s)ds

    which is a finite real constant thanks to the assumption on the function rv. Let A(θ) and μ(θ) be two strictly positive constants such that A(θ)+μ(θ)<α(θ+1) if θ<1 and A(θ)+μ(θ)=2(α1) if θ=1. Since γ(t)αtθ, we have

    fθ(t)(1+θ)(θtθ1α).

    Therefore, up to take t1 large enough in the case θ<1, we can assume that

    fθ(t)A(θ)μ(θ) tt1. (3.5)

    Using now the fact that

    |hv(t)|x(t)x(t)v2W(t)hv(t),

    it follows, from the estimate (3.5), that

    λ(t)hv(t)+fθ(t)hv(t)A(θ)hv(t)+2λ(t)W(t)hv(t)μ(θ)hv(t).

    Applying now the elementary inequality

    bxax2b24a   a>0, (x,b)R2

    with x=hv(t), we deduce that for every tt1

    λ(t)hv(t)+fθ(t)hv(t)(λ(t))2W(t)A(θ)μ(θ)hv(t)=B(θ,t)λ(t)W(t)μ(θ)hv(t), (3.6)

    where

    B(θ,t)=(θ+1)2tθ1A(θ).

    Inserting (3.6) in the inequality (3.2), we obtain

    (13K)tt1λ(s)W(s)ds+(3KB(θ,t))λ(t)W(t)+μ(θ)hv(t)C0+tt1gθ(s)hv(s)ds (3.7)

    Let us notice that if 0θ<1 then limt+B(θ,t)=0 and in the case where θ=1, since α>3, one can choose 0<μ(1)<23(α3) to get

    3KB(1,t)=3α4A(1)>0.

    Hence, up to take t1 large enough, we can assume that, for every 0θ1, there exists a constant ν(θ)>0 such that

    3KB(θ,t)ν(θ), for all tt1.

    In particular, the inequality (3.7) implies

    μ(θ)hv(t)C0+tt1gθ(s)hv(s)ds.

    Recalling that the function gθ is integrable over [t1,+) and applying the Gronwall's lemma, we deduce that the function hv is bounded. Hence, by going back to the inequality (3.7), we infer that

    suptt1λ(t)W(t)<+

    and

    +t1λ(s)W(s)ds<+. (3.8)

    Now, using the fact that the energy function W is decreasing, we deduce from (3.8) that t1+θW(t)0 as t+; in fact, for every t2t1, we have

    (1+θ)(t2)1+θW(t)tt2λ(s)W(s)ds.

    To conclude, it remains to prove that limt+hv(t) exists. From (2.3), the function hv satisfies the differential inequality

    hv(t)+γ(t)hv(t)ζ(t)

    where

    ζ(t)=32x(t)2+rv(t).

    The assumption on the function rv and the estimate (3.8) imply that tθζ(t)L1(t0,+;R+), thus the existence of limt+hv(t) follows from the following lemma.

    Lemma 3.2. Let a>0 and w:[a,+)R+ be a continuous function satisfying

    w(t)αtθta

    where α and θ are nonnegative constants with 0θ1 and α>1 if θ=1. Let φC2(a,+;R+) satisfy a differential inequality

    φ(t)+w(t)φ(t)ψ(t) (3.9)

    with tθψ(t)L1(a,+;R+). Then limt+φ(t) exists.

    Proof. Multiplying each member of the differential inequality (3.9) by eΓ(t,a), where

    Γ(t,s)=tsw(τ)dτ,

    and integrating on [t0,t], we get

    φ(t)eΓ(t,a)φ(a)+taeΓ(t,s)ψ(s)ds. (3.10)

    Applying Fubini's Theorem and using [Lemma 3.14, [8]], we deduce that there is a real constant M>0 such that

    +ataeΓ(t,s)ψ(s)dsdtM+asθψ(s)ds.

    We therefore infer from (2.2) that the positive part [φ]+ of φ belongs to L1(a,+;R+) which implies that limt+φ(t) exists.

    Before starting the proof of Theorem 1.3, let us recall the classical Opial's lemma.

    Lemma 3.3 (Opial's lemma).Let x:[t0,+)H. Assume that there exists a nonempty subset S of H such that:

    i) if tn+ and x(tn)x weakly in H, then xS,

    ii) for every zS, limt+x(t)z exists.

    Then there exists zS such that x(t)z weakly in H as t+.

    For a simple proof of Opial's lemma, we refer the reader to [13].

    Proof of Theorem 1.3. :

    Step 1: Proof of the property (P1).

    According to Lemma 3.1, limt+hv(t) exists for every vSΦ and the energy function W satisfies (1.5). Let tn+ such that x(tn) converges weakly in H to some ˉx. Since Φ(x(t))Φ as t+, the weak lower semi-continuity of the convex function Φ implies that Φ(ˉx)Φ which means that ˉxSΦ. By Opial's lemma, we deduce that x(t) converges weakly in H as t+ to some element of SΦ.

    Step2: Proof of the property (P2).

    Let vS=SΦSU. Since rv=0, Lemma 3.1 implies that limt+hv(t) exists and W satisfies (1.5). Thus, in view of the assumption lim inft+tθ+1ε(t)>0, we have U(x(t))U as t+. Therefore, as in the above step, the weak lower semi-continuity of the convex functions Φ and U yields that every sequential weak cluster point of x(t),as t+, belongs to the subset S. This completes the proof of the property (P2) thanks to Opial's lemma.

    This section is devoted to the proof of Theorem 1.4. Before proving separately the two properties (Q1) and (Q2), let us first recall some general facts about strongly convex functions and a Tikhonov approximation method [15]. Since the function U is strongly convex, there exists a positive real m such that U(x)m2x2 is convex (we say that U is mstrongly convex). Moreover, for every nonempty, convex and closed subset C of H, the function U has a unique minimizer xC on C. Let x be the minimizer of U on H and p its minimizer on SΦ. For every tt0, we consider the function Φt defined on H by

    Φt(x)=Φ(x)+ε(t)U(x).

    Clearly, Φt is ε(t)m-strongly convex. Therefore, Φt satisfies the convex inequality

    Φt(z)Φt(y)+Φt(y),zy+m2ε(t)zy2, (4.1)

    and has a unique global minimizer which we denote by xε(t). Adopting the Tikhonov method, we can prove that xε(t) converges strongly to p as t+. Indeed, since

    Φt(xε(t))Φt(p) (4.2)

    and

    Φ(p)Φ(xε(t)),

    then

    U(xε(t))U(p). (4.3)

    Furthermore, since U is strongly convex, U(x)+ as x+; hence from the inequality 4.3 we deduce that (xε(t))tt0 is bounded. So, let ˜xH be a weak limit of a sequence (xε(tn)) where tn+. Using the weak lower semi-continuity of the two convex functions Φ and U and letting t=tn+ in the inequalities (4.2) and (4.3), we deduce that Φ(˜x)Φ(p) and U(˜x)U(p) which is, from the definition of p, is equivalent to ˜x=p. Consequently, we infer that xε(t) converges weakly to p as tn+. Now, since U is mstrongly convex, we have

    U(xε(t))U(p)+U(p),xε(t)p+m2xε(t)p2.

    Hence, by (4.3), we conclude that limt+xε(t)p=0 which completes the proof of the claim.

    Proof of Theorem 1.4. Let us first prove the assertion (Q1). We consider the function h(t)=hp(t)=12x(t)p2. Using the equation (GAVDγ,ε) and the convex inequality (4.1), we obtain

    h(t)+γ(t)h(t)=x(t)2+Φt(x(t)),px(t)x(t)2+Φt(p)Φt(x(t))mε(t)h(t)x(t)2+Φt(p)Φt(xε(t))mε(t)h(t)x(t)2+ε(t)(U(p)U(xε(t)))mε(t)h(t). (4.4)

    In the last inequality we have used the fact that p is also a minimizer of Φ. Set

    σ(t)U(xε(t))U(p)+mh(t).

    The inequality (4.4) becomes

    h(t)+γ(t)h(t)+ε(t)σ(t)x(t)2. (4.5)

    Let us prove that lim inft+h(t)=0. We argue by contradiction. As consequence of

    limt+U(xε(t))U(p)=0,

    there exists t2t0 large enough and μ>0 such that σ(t)μ for every tt2. Therefore the differential inequality (4.5) implies that, for every tt2, we have

    h(t)+μtt2τt2eΓ(τ,s)ε(s)dsdτh(t2)+tt2eΓ(τ,t2)dτh(t2)+tt2τt2eΓ(τ,s)x(s)2dsdτ,

    where

    Γ(t,s)=tsγ(τ)dτ.

    Applying Fubini's theorem, we then infer that

    μ+t2ε(s)+seΓ(τ,s)dτdsh(t2)+|h(t2)|+t2eΓ(τ,t2)dτ++t2x(s)2+seΓ(τ,s)dτds. (4.6)

    Since γ(t)=αtθ, a simple integration by parts ensures the existence of two real constants Bθ>Aθ>0 so that

    Aθsθ+seΓ(τ,s)dτBθsθ, for every st0.

    Hence, by combining the inequality (4.6) and the assumption

    +t0sθx(s)2ds<+,

    we get

    +t0sθε(s)ds<+,

    which contradicts our assumption on the function ε(.). We therefore deduce that

    lim inft+h(t)=0. (4.7)

    Now let us suppose that

    lim supt+h(t)>0. (4.8)

    The continuity of the function h combined with (4.7) and (4.8) ensures the existence of two real numbers λ<δ and two positive real sequences (sn)n and (tn)n such that for every nN we have

    max{t,n}<sn<tn,h(tn)=δ,h(sn)=λ,h(s)[λ,δ] on [sn,tn],

    where t>t2 is a fixed positive number such that for every tt

    U(xε(t))U(p)mλ.

    We deduce from (4.5) that for every nN and for all t[sn,tn]

    h(t)+αtθh(t)x(t)2.

    Multiplying the last differential inequality by tθ and integrating over [sn,tn], we obtain

    tθnh(tn)sθnh(sn)+θsθ1nλθtθ1nδ+α(δλ)+θ(θ1)tnsntθ2h(t)dttnsntθx(t)2. (4.9)

    Using now the three facts

    |h(tn)|x(tn)2h(tn)=x(tn)2δ,|h(sn)|x(sn)2λ,tnsntθ2h(t)dtδsθ1n1θ if 0θ<1,

    and letting n goes to + in the the inequality (4.9), we get

    (α1)(δλ)0 if θ=1,α(δλ)0 if 0θ<1.

    This contradicts the assumption δ>λ. We therefore conclude that limt+h(t)=0, which completes the proof of the assertion (Q1).

    In order to prove the assertion (Q2), we first apply 3.1 with v=x to deduce that the energy function W satisfies the asymptotic property (1.5) which implies in particular the solution x(.) fulfils the assumption in the previous assertion (Q1). Therefore, we conclude that x(t) converges strongly to p which is, in the present case, equal to x. This ends the proof of our main theorem.

    Let us now prove the theorem 1.5 of Attouch and Cazernicki.

    Proof. The dynamical system (1.6) is the particular case of the general system (GAVDγ,ε) corresponding to θ=0 and U(x)=12x2. It follows from Theorem 1.2 that any solution x(.) of the dynamical system (1.6) satisfies the assumptions of the assertion (Q1) of Theorem 1.4. Therefore, we deduce that x(t) converges strongly to the unique minimizer p of the function of U. To conclude we notice that, in this case, p is the projection of zero on SΦ.

    Remark 4.1. Without the regularizing term ϵ(t)U(x(t)) and under appropriate assumptions on the damping term, the trajectories of the dynamical system (GAVDγ,ε) weakly converge to a non-specified minimizer of the objective functional Φ (See Theorem 1.3 with ϵ(t)=0). As it is shown in the main theorem (Theorem 1.4), if ϵ(t) vanishes slowly at infinity, then te regularizing term ϵ(t)U(x(t)) forces the trajectories of the differential system (GAVDγ,ε) to converges strongly to a particular minimizer of Φ. In this sense, The added term ϵ(t)U(x(t)) may be considered as a stabilizer factor for our dynamical system.

    The authors acknowledge the Deanship of Scientific Research at King Faisal University for the financial support under the annual research project (Grant No. 170065). The authors are also indebted to the referee for his careful reading and insightful comments.

    The authors declare that they have no conflict of interest.



    [1] F. Alvarez, On the minimizing properties of a second order dissipative system in Hilbert spaces, SIAM J. Cont. Optim., 38 (2000), 1102-1119.
    [2] H. Attouch, A. Cabot, Asymptotic stabilization of inertial gradient dynamics wuth time-dependent viscosity, J. Differ. Equations, 263 (2017), 5412-5458.
    [3] H. Attouch, Z. Chbani, H. Riahi, Combining fast inertial dynamics for convex optimization with Tikhonov regularization, J. Math. Anal. Appl., 457 (2018), 1065-1094.
    [4] H. Attouch, Z. Chbani, J. Peypouquet, P. Redont, Fast convergence of inertial dynamics and algorithms with asymptotic vanishing viscosity, Math. Program., 168 (2018), 123-175.
    [5] H. Attouch, M. O. Czarnecki, Asymptotic control and stabilization of nonlinear oscillators with non-isolated equilibria, J. Differ. Equations, 179 (2002), 278-310.
    [6] H. Attouch, X. Goudou, P. Redont, The heavy ball with friction method, Ⅰ: The continuous dynamical system: Global exploration of the local minima of a real valued function by asymptotic analysis of a dissipative dynamical system, Commun. Contemp. Math., 2 (2000), 1-34.
    [7] A. Cabot, The steepest descent dynamical system with control. Applications to constrained minimization, Control Optimization and Calculus of Variations, 10 (2004), 243-258.
    [8] A. Cabot, P. Frankel, Asymptotics for some semilinear hyperbolic equations with non-autonomous damping, J. Differ. Equations, 252 (2012), 294-322.
    [9] A. Haraux, M. A. Jendoubi, On a second order dissipative ODE in Hilbert space with an integrable source term, Acta Math. Sci., 32 (2012), 155-163.
    [10] M. A. Jendoubi, R. May, On an asymptotically autonomous system with Tikhonov type regularizing term, Arch. Math., 95 (2010), 389-399.
    [11] R. May, Asymptotic for a second order evolution equation with convex potential and vanishing damping, Turk. J. Math., 41 (2017), 681-685.
    [12] Y. Nestrov, A method of solving a convex programming problem with convergence rate O(1/$k^2$), Sov. Math. Dokl., 27 (1983), 372-376.
    [13] Z. Opial, Weak convergence of the sequence of successive approximation for nonexpansive mappings, Bull. Amer. Math. Soc., 73 (1967), 591-597.
    [14] W. Su, S. Boyd, E. J. Candès, A differential equation for modeling Nestrov's accelerated gradient method: Theory and Insights, J. Mach. Learn. Res., 17 (2016), 1-43.
    [15] A. N. Tikhonov, V. Y. Arsenin, Solutions of Ill-Posed Problems, New York: Halsted Press, 1977.
  • This article has been cited by:

    1. Haixin Ren, Bin Ge, Xiangwu Zhuge, Fast Convergence of Inertial Gradient Dynamics with Multiscale Aspects, 2023, 196, 0022-3239, 461, 10.1007/s10957-022-02124-w
    2. Hedy Attouch, Aïcha Balhag, Zaki Chbani, Hassan Riahi, Damped inertial dynamics with vanishing Tikhonov regularization: Strong asymptotic convergence towards the minimum norm solution, 2022, 311, 00220396, 29, 10.1016/j.jde.2021.12.005
    3. Gangfan Zhong, Xiaozhe Hu, Ming Tang, Liuqiang Zhong, Fast Convex Optimization via Differential Equation with Hessian-Driven Damping and Tikhonov Regularization, 2024, 203, 0022-3239, 42, 10.1007/s10957-024-02462-x
    4. Hedy Attouch, Aïcha Balhag, Zaki Chbani, Hassan Riahi, Accelerated Gradient Methods Combining Tikhonov Regularization with Geometric Damping Driven by the Hessian, 2023, 88, 0095-4616, 10.1007/s00245-023-09997-x
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2915) PDF downloads(197) Cited by(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog