arxiv: v1 [math.ds] 21 Sep 2017

Transcription

1 UNBOUNDED AND BLOW-UP SOLUTIONS FOR A DELAY LOGISTIC EQUATION WITH POSITIVE FEEDBACK arxiv: v [math.ds] 2 Sep 207 ISTVÁN GYŐRI, YUKIHIKO NAKATA, AND GERGELY RÖST Abstract. We stuy boune, unboune an blow-up solutions of a elay logistic equation without assuming the ominance of the instantaneous feeback. It is shown that there can exist an exponential (thus unboune) solution for the nonlinear problem, an in this case the positive equilibrium is always unstable. We obtain a necessary an sufficient conition for the existence of blow-up solutions, an characterize a wie class of such solutions. There is a parameter set such that the non-trivial equilibrium is locally stable but not globally stable ue to the co-existence with blow-up solutions.. Introuction There is a vast literature on the stuy of elay logistic equations escribing population growth of a single species [6, ]. In [4], the qualitative stuies of such logistic type elay ifferential equations have been summarize, see also [, 4, 7, 2, 6] an references therein. In [2] the authors stuy the global asymptotic stability of a logistic equation with multiple elays. Their global stability result is generalize in Theorem 5.6 in Chapter 2 in [], see also the iscussion in [7]. Those conitions presente in [2] an in Theorem 5.6 in Chapter 2 in [] are elay inepenent conitions, exploiting the ominance of the instantaneous feeback. In [7], applying the oscillation theory of elay ifferential equations [9], the first author of this paper obtains a global stability conition for a logistic equation without assuming the ominance of the instantaneous feeback. See also [9] an references therein for the stuy of logistic equations without instantaneous feeback. When the instantaneous feeback term is small compare to the elaye feeback, the positive equilibrium is not always globally stable. Our motivation of this note comes from interesting examples shown in [7], where the author shows the existence of an exponential solution for a logistic equation with elay. Here we wish to investigate the properties of positive solutions of such an equation in etail. To be more specific, we consier the logistic equation (.) x(t) = rx(t) ( + αx(t) x(t )), t where r > 0, α R. The equation (.) is a special case of the equations stuie in [7]. Note that (.) is a normalize form of the following elay ifferential equation (.2) s N(s) = N(s) ( r + an(s) bn(s τ)),

2 UNBOUNDED AND BLOW-UP SOLUTIONS FOR A DELAY LOGISTIC EQUATION 2 where r > 0, a R an b > 0. If we efine y(s) := b r N(s) an α := a b, then we obtain y(s) = ry(s) ( + αy(s) y(s τ)). s Next we scale the time so that the elay is normalize to be one, by letting s := tτ an x(t) := y(s), then y(s τ) = y(τ(t )) = x(t ) an by calculating t x(t), we obtain (.) with r = τ r. For (.) we show that there exist some unboune solutions, when α is allowe to be positive. More precisely, it is shown that a blow-up solution (i.e. a solution that converges to infinity in finite time) exists if an only if α > 0 hols. We then show that an exponential solution, namely x(t) = ce rt, c > 0, exists when α = e r hols, which is an unboune but not blow-up solution. The case is further elaborate, as we also fin solutions which blows up faster than the exponential solutions. This paper is organize as follows. In Section 2 we collect previous results on bouneness an stability, which are known in the literature, with the exception about the existence of blow-up solutions. In Section 3, we focus on the exponential solution for the nonlinear ifferential equation (.) an its relation to stability. In Section 4 we characterize a large class of initial functions that generate superexponential blow-up solutions, an we also fin a class of subexponential solutions. Section 5 is evote to a summary an iscussions. 2. Bouneness, stability an blow-up Denote by C the Banach space C([, 0], R) of continuous functions mapping the interval [, 0] into R an esignate the norm of an element φ C by φ = sup θ 0 φ(θ). The initial conition for (.) is a positive continuous function given as x(θ) = φ(θ), θ [, 0]. There exists a unique positive equilibrium given by x = α if an only if α < hols. In the following theorem we characterize global an local ynamics of the solutions. The result on the existence of a blow-up solution seems to be new. Theorem. The following statements are true. () If (2.) α, then the positive equilibrium is globally asymptotically stable. (2) If < α <, then the positive equilibrium is locally asymptotically stable for α (2.2) r < + α arccos(α), an it is unstable for (2.3) r > Moreover, α + α arccos(α).

3 UNBOUNDED AND BLOW-UP SOLUTIONS FOR A DELAY LOGISTIC EQUATION 3 (a) If < α 0, then every solution is boune. (b) If 0 < α, then there exists a blow-up solution in a finite time. (3) If α, then for every solution, which exists globally, one has lim sup t x(t) =. Proof. ) For the global stability of the equilibrium we refer to the proof of Theorem 5.6 in Chapter 2 in []. 2) The result for local asymptotic stability is well-known, see for example Theorems 2 an 3 in [4]. 2-a) Notice that a solution of (.) satisfies x(t) = x(0)e r t +αx(s) x(s )s 0, hence positive solutions remain positive. When < α < 0, the result follows from a simple comparison principle applie for the inequality x (t) x(t)r( αx(t)). In the case α = 0 we obtain the Wright s equation for which bouneness is known, see [5]. 2-b) Let us assume that α > 0. We show that there exists a solution that blows up in a finite time. Consier a positive continuous initial function satisfying {, θ [, φ(θ) = 2], q, θ = 0, where q is a positive constant to be etermine later. Since one has it hols x(t ) = 0, 0 t 2, (2.4) x (t) = rαx(t) 2, 0 t 2. Then equation (2.4) is easily integrate as for t < qrα then we have x(t) = q rαt, where the solution exists. Let us set lim x(t) = lim t h t h q = h rα, h 2, rα( =, h t) so the solution blows up at t = h < 2. 3) The result follows from Theorem 5. in [7]. 3. Instability an exponential solutions A remarkable feature of the logistic equation with positive feeback (.) is the possible existence of exponential solutions, espite the equation being nonlinear. As in [7], we fin the following.

4 UNBOUNDED AND BLOW-UP SOLUTIONS FOR A DELAY LOGISTIC EQUATION 4 r 5 Unstable 0 x x xxxx xxxx GAS xxxxx xx xx xxx xxx 5 xxxx xxxxxx š xxxx xx e rt Stable xxxx xxx _ Figure 2.. Stability region for the positive equilibrium in the (α, r)-parameter plane. The shae region is the stability region given by (2.) an (2.2). The positive equilibrium is globally stable for α an is unstable above the stability bounary. Exponential solutions exist on the enote curve. Blow-up solutions exist for α > 0, hence we can observe a region where the positive equilibrium is locally stable yet blow-up solutions also exist. Proposition 2. There exists an exponential solution (3.) x c (t) := ce rt, t of the generalize logistic equation ( (3.2) x(t) = rx(t) + t if an only if (3.3) ) k a i x(t τ i ), i= k a i e rτi = 0. i= The proof is straightforwar, an for (.) it means that an exponential solution exists if an only if (3.4) α = e r hols. It can be shown that the existence of the exponential solution implies the instability of the positive equilibrium. Proposition 3. Let us assume that (3.4) hols. Then the positive equilibrium of (.) is unstable.

5 UNBOUNDED AND BLOW-UP SOLUTIONS FOR A DELAY LOGISTIC EQUATION 5 Proof. We compare the two conitions (3.4) an (2.3). We set (3.5) ω = arccos(α) for α (0, ). Note that ω ( ) 0, π 2. Using the parameter transformation (3.5) we get α arccos(α) + α = ω cos ω sin ω an the conition (3.4) is written as r = ln(cos ω). Define g (ω) := ω cos ω sin ω, g 2 (ω) := ln(cos ω) for ω ( ) 0, π 2. We claim that ( (3.6) g 2 (ω) > g (ω), ω 0, π ). 2 It is easy to see that lim ω 0 g (ω) = lim ω 0 g 2 (ω) = 0. Straightforwar calculations show g cos ω ( (ω) = + ω ), sin ω sin ω g 2(ω) = sin ω cos ω. Then we see g 2(ω) g (ω) { ( = sin 2 ω cos ω ( cos ω) + ω )} cos ω sin ω sin ω = cos ω { ( ( + cos ω) cos ω + ω )} cos ω sin ω sin ω = cos ω ( ω cos ω ) cos ω sin ω sin ω > 0 for ω ( 0, π 2 ). Thus we get (3.6) an obtain the conclusion. In Figure 2. we visualize the conition (3.4) in (α, r) parameter plane. Figure 2. clearly shows that if (3.4) hols, then the positive equilibrium is unstable. 4. A new class of blow-up solutions We investigate other solutions of (.) when (3.4) hols. Theorem 4. Let the conition (3.4) hol. Consier a solution with the initial function satisfying (4.) 0 < φ(s) ce rs, s [, 0] with (4.2) (4.3) φ(0) = c, φ( ) < ce r for c > 0. Then one has x(t) > x c (t) = ce rt for t > 0 an the solution blows up at a finite time.

6 UNBOUNDED AND BLOW-UP SOLUTIONS FOR A DELAY LOGISTIC EQUATION 6 Proof. Looking for a contraiction, we assume that x(t) exists on [0, ). Then x(t) > 0, t. Define ( ) ( ) x(t) x(t) z(t) := ln = ln x c (t) ce rt, t. Let ψ(s) := ln Then z(s) = ψ(s) for s [, 0], an we have ( ) φ(s) ce rs, s [, 0]. 0 = ψ(0) ψ(s), 0 = ψ(0) > ψ( ), s [, 0] from (4.), (4.2) an (4.3). Now we obtain the relation cert x (t)ce rt x(t)rce rt z(t) = t x(t) c 2 e 2rt which can be rewritten by using (.) an α = e r as t z(t) = r(e r x(t) x(t )), = x (t) x(t) r, an by x(t) = e z(t) ce rt we obtain a nonautonomous ifferential equation for z: (4.4) z (t) = rce r(t ) ( e z(t) e z(t )), t > 0 using (3.4). First we show that z (t) > 0 for any t 0. Since ( ) ( ) φ(0) φ( ) ψ(0) = ln = 0 > ψ( ) = ln c ce r follows from (4.2) an (4.3), one fins z (0) = rce r ( e ψ(0) e ψ( )) > 0. Assume that there exists t > 0 such that z (t) > 0 for 0 t < t an z (t ) = 0 hol. If t (0, ) then, since t (, 0), while z (t ) = rce r(t ) ( e z(t) e z(t )), z(t ) z(t ) = z(t ) z(0) + ψ(0) ψ(t ) > 0, thus we obtain a contraiction. If t, then z(t ) z(t ) = t t z (s)s > 0, which leas a contraiction again. Therefore, we obtain z (t) > 0 for t 0. We can fix a T > 2 such that < ( α) re r ce rt. Since z(t) exists on [0, ) an z (t) > 0 for t 0, z (t) > 0 for 0 t T. Thus m := min 0 t T z (t) > 0.

7 UNBOUNDED AND BLOW-UP SOLUTIONS FOR A DELAY LOGISTIC EQUATION 7 By the intermeiate value theorem, for each 0 t T, there exists ξ(t) [z(t ), z(t)] such that e z(t) e z(t ) = e ξ(t) (z(t) z(t )). Since ξ(t) z(t ) > 0, we have This yiels (4.5) z (t) > z (t) = rce r(t ) e ξ(t) (z(t) z(t )) rce r(t ) (z(t) z(t )) = rce r(t ) t t t t z (s)s. z (s)s, t T an hence z (T ) > T T z (s)s m. This implies that z (t) > m for t T. Otherwise there is a t > T > such that z (t) > m for 0 < t < t an z (t ) = m. But from (4.5) which is a contraiction. Thus for any t we have Therefore, or equivalently z (t ) > t t z (s)s > m, e z(t) e z(t ) = e z(t) ( e (z(t) z(t ))) = e z(t) ( e t e z(t) ( e m). t z (s)s z (t) rce r(t ) e z(t) ( e m), t, z (t)e z(t) ( rce r ( e m ) ) e rt, t. Integrating both sies of the above equation, This yiels t z (s)e z(s) s = [ e z(s)] s=t ) s= = e z() e z(t) ( rce r ( e m ) ) t e rs s = ( ce r ( e m ) ) (e rt e r ), t >. e z() ( ce r ( e m ) ) (e rt ) + e z(t) > ( ce r ( e m ) ) (e rt ), t, which is a contraiction. Therefore, z oes not exist on [0, ), moreover z (t) 0, t 0. Consequently we shoul have a T (0, ) such that lim t T z(t) = +, that is it is a blow-up solution. Corresponing to this blow-up solution, we also have that is x(t) is also a blow-up solution. x(t) = ce z(t)+rt, t T,

8 UNBOUNDED AND BLOW-UP SOLUTIONS FOR A DELAY LOGISTIC EQUATION 8 Proposition 5. The following estimate is vali: ) ( r ln + er T, c where T is the blow-up time for a blow-up solution x in Theorem 4. Proof. We can fin the lower boun for the blow-up time T by the stanar comparison principle. Consier the following orinary ifferential equation y (t) = y(t)r ( + e r y(t) ) with y(0) = c = x(0). By the comparison theorem, we have x(t) y(t), t 0. Integrating the equation, we get y(t) = ce rt + ( e rt ) e r c for sufficiently small t. From this expression, we fin the finite blow-up time for y an then we obtain the require estimation. Similar to the proof of Theorem 4, we obtain the following theorem. Theorem 6. Let the conition (3.4) hol. Consier a solution with the initial function satisfying (4.6) ce rs φ(s), s [, 0] with (4.7) (4.8) φ(0) = c, φ( ) > ce r for c > 0. Then one has x(t) < x c (t) = ce rt for t > 0, an consequently the solution exists on [, ). x(t) x c(t) For the initial functions consiere in Theorems 4 an 6, is a monotone function for t > 0, thus the orer of the solution with respect to the exponential solution, x c (t) = ce rt, is preserve. We o not analyze the qualitative behavior of the solution with the initial conition that oscillates about the exponential solution. Numerical simulations suggest that, for many solutions, monotone function. 5. Discussion x(t) x c(t) eventually becomes a In this paper we stuy the logistic equation (.). In the stability analysis of elaye logistic equations, negative instantaneous feeback is usually assume, see [4, 7, 4, 6] an references therein. Only a few stability results are available in the literature for the case of positive instantaneous feeback e.g., [2, 3]. However, the blow-up solutions, which are present ue to the positive instantaneous feeback, have not been analyze in etail, since the publication of the paper [7]. This manuscript has been inspire by the work one in [7], especially, paying attention to the examples an open questions given in Section 5 of the paper [7]. Our primary goal was to clarify an unerstan a relation between the stability conition of the positive equilibrium an the existence conition for the exponential solution. For the logistic equation (.), we show that the existence of the exponential solution implies instability of the positive equilibrium in Proposition 2, see also Fig. 2..

9 UNBOUNDED AND BLOW-UP SOLUTIONS FOR A DELAY LOGISTIC EQUATION 9 Since stability analysis becomes extremely har for the ifferential equation with multiple elays, the comparison of the existence conition of the exponential solution to the stability conition is not straightforwar in general, thus it remains an open problem whether the positive equilibrium of (3.2) is always unstable whenever exists an (3.3) hols. Fining a global stability conition for (.) in the case of < α < 0 is still an open problem. For α = 0 the global stability problem is known as the famous Wright conjecture [2]. On the other han, for 0 < α <, ue to the existence of the blow-up solution, it is shown that the stable equilibrium can not attract every solution, thus there is no hope to obtain global stability conition for 0 < α <. Numerical simulations also suggest that there are many boune an oscillatory solutions. In Theorems 4 an 6 we fix the parameters as in (3.4) in Proposition 2 so that the exponential solution exists for the logistic equation (.). We consier some solutions that preserve the orer with respect to the exponential solution, an show that some blow up, while others exist for all positive time. The qualitative behaviour of the solution with the initial conition that oscillates about the exponential solution is not stuie. For such an initial function, careful estimation of the solution seems to be necessary to unerstan the long term solution behaviour. x(t) x c(t) Numerical simulations suggest that, for many solutions, eventually becomes a monotone function. The etaile unerstaning of the evolution of such solutions is also left for future work. Finally, one might be intereste in the equation y(t) = y(t)r ( + αy(t) + y(t )), t which has an opposite sign for the elaye feeback term. For this equation the qualitative ynamics is stuie in the literature. The positive equilibrium y = + α exists if an only if α <. Accoring to Theorem 5.6 in Chapter 2 in [], the positive equilibrium is globally asymptotically stable. When α the solutions are unboune, see again Theorem 5. in [7]. Acknowlegement. The first author s research has been supporte by the Hungarian National Research Fun Grant OTKA K2086 an the Széchenyi 2020 project EFOP The secon author was supporte by JSPS Grant-in-Ai for Young Scientists (B) 6K The thir author was supporte by Hungarian National Research Fun Grant NKFI FK 2406 an MSCA-IF The meeting of the authors have been supporte by JSPS an NKFI Hungary-Japan bilateral cooperation project. References [] J.A.D. Appleby, I. Győri, D.W. Reynols, History-epenent ecay rates for a logistic equation with infinite elay. Proc. Roy. Soc. Einburgh Sect. A 4 (20), no., [2] B. Bánhelyi, T. Csenes, T. Krisztin, an A. Neumaier,,Global attractivity of the zero solution for Wright s equation, SIAM Journal on Applie Dynamical Systems, 3(204), [3] O. Diekmann, S.A. van Gils, S.M.V. Lunel, H.O. Walther, Delay Equations Functional, Complex an Nonlinear Analysis, Springer Verlag (99).

10 UNBOUNDED AND BLOW-UP SOLUTIONS FOR A DELAY LOGISTIC EQUATION 0 [4] T. Faria, E. Liz, Bouneness an asymptotic stability for elaye equations of logistic type. Proc. Roy. Soc. Einburgh Sect. A 33 (2003), no. 5, [5] E. Liz, G. Röst, Dichotomy results for elay ifferential equations with negative Schwarzian, Nonlinear Analysis: Real Worl Applications.3 (200): [6] K. Gopalsamy, Stability an oscillation in elay ifferential equations of population ynamics. Kluwer Acaemic Publishers (992). [7] I. Győri, A new approach to the global asymptotic stability problem in a elay Lotka-Volterra ifferential equation. Mathematical an Computer Moelling 3 6 (2000) pp [8] I. Győri, F. Hartung, Funamental solution an asymptotic stability of linear elay ifferential equations, Dynamics of Continuous, Discrete an Impulsive Systems, 3:2 (2006) [9] I. Győri, G. Laas, Oscillation theory of elay ifferential equations with applications, Clarenon Press, Oxfor (99). [0] X. He, Global stability in nonautonomous Lotka-Volterra systems of pure-elay type, Differential an Integral Equations.2 (998) [] Y. Kuang, Delay ifferential equations with applications in population ynamics, Acaemic Press, San Diego (993). [2] S.M. Lenhart, C.C. Travis, Global stability of a biological moel with time elay, Proc. Amer. Math. Sot. 96 (986) [3] H. Li, R. Yuan, An affirmative answer to the extene Gopalsamy an Liu s conjecture on the global asymptotic stability in a population moel. Nonlinear Anal. Real Worl Appl. (200), no. 5, [4] S. Ruan, Delay ifferential equations in single species ynamics. In: Delay Differential Equations an Applications, Springer (2006), [5] G. Stépán, Retare Dynamical Systems: Stability an Characteristic Function, Wiley, New York, 989. [6] Z. Teng, Permanence an stability in non-autonomous logistic systems with infinite elay. Dyn. Syst. 7 (2002), no. 3, (István Győri) University of Pannonia, Veszprem, Hungary (Yukihiko Nakata) Shimane University, Japan (Gergely Röst) University of Szege, Hungary an University of Oxfor, UK aress: rost@math.u-szege.hu