LINEARIZED STABILITY PRINCIPLE FOR DIFFERENTIAL EQUATIONS WITH DELAYS

Transcription

1 DOI: /ctcr LINEARIZED STABILITY PRINCIPLE FOR DIFFERENTIAL EQUATIONS WITH DELAYS Leonid Berezansky 1, *, 2, ** Elena Braverman 1 Ben-Gurion University of the Negev, Beer-Sheva, Israel, 2 University of Calgary, Calgary, Canada In this article а linearized global stability principle is announced for nonlinear delay differential equations which is illustrated by several models of Population Dynamics. Is given а review of some mathematical models with possible applications of the linearized principle is presented. Keywords: delay differential equations, linearized global stability principle. 1. Introduction Many applied scientists ask the following question: why do we need to study differential equations with delays when so much is known about equations without delays and these equations are much easier? The answer is: because so many processes, both natural and man-made, in biology, medicine, chemistry, engineering, economics etc. involve time delays. Whether you like it or not, time delays occur so often in almost every situation that to ignore them means to ignore the reality. A simple example in nature is reforestation. After cutting a forest and replanting, it will take at least 20 years before reaching any kind of maturity. Hence any mathematical model of forest harvesting should involve time delays. Similarly any model of population dynamics without delays is an approximation at best. Another motivation for studying DDE comes from their applications in feedback control theory. In his study of ship stabilization Minorsky (1942) pointed out very clearly the importance of incorporating delays in the feedback mechanism. We would like to present here several applied models with time delays. 1. The model of population growth (the delay logistic equation) () = rn(t) 1. involves a constant concentrated delay. This equation was first introduced by Hutchinson in Another equation is the logistic equation with a distributed delay () = rn(t) 1 N(s)G(t s)ds. This equation was introduced in 1978 by MacDonalds. 2. Lotka-Volterra predator-prey system with distributed delays x (t) = x(t)(a bx(t) y (t) = y(t)( c + dy(t) was introduced in F (s)y(t + s)ds, F (s)x(t + s)ds, 3. Delay models in medicine allowed to describe dynamic diseases which correspond to sustainable oscillations of the following equations: * Partially supported by the Israeli Ministry of Absorption. ** Partially supported by the NSERC Research Grant. Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника» Т. 17, 3. С

2 Краткие сообщения x (t) = a () (), () x (t) = cx(t), () x (t) = () cx(t). () These equations were introduced by Mackey and Glass in Ship stabilization model x (t) + a(t)x (t h) + b(t)x(t g) = f(t) was introduced by Minorsky in Comparison Between Ordinary and Delay Differential Equations Consider a simple linear equation without delay x (t) = x(t); x(0) = x > 0. This equation has exact solution x(t) = x e, this solution is positive, monotone and tends to zero as t. Consider now a simple linear equation with delay x (t) = x(t 2). To find a solution of this equation we need to assign an initial function instead of an initial value. Suppose that we have the following initial function x(t) = 1; t [ 2; 0]. There is no exact formula for the solution of this problem, but we can find this solution by the method of steps. Suppose first that t ϵ [0; 2], then x (t) = 1; x(0) = 1, hence x(t) = 1 t. If we continue this process we obtain the solution which oscillates about zero and is unstable. So the behavior of solutions for equations with delay and without delay is quite different. Actually the behavior of solutions also depends on the value of delay. In particular, consider the equation with any constant delay x (t) = x(t τ). If τ, then the solution of this equation with the same initial function as above is positive, monotone and tends to zero. But there are also oscillating solutions of this equation. If τ, then all solutions for any initial function oscillate and tend to zero. If τ, then all solutions oscillate and the equation is unstable. There are several important questions on the asymptotic behavior of solutions of delay differential equations. Here we consider only one property which is a global stability of the equation. What does it mean? Consider the nondelay logistic equation N (t) = r(t)n(t) 1 (). The carrying capacity K is a positive equilibrium. It is a simple fact that for r(t) > r > 0 this equilibrium is locally asymptotically stable. But what is more interesting, this solution attracts all positive solutions. In this case we say that x(t) = K is a global attractor or that the equation is globally asymptotically stable. 3. Nonlinear and Linear Equations Consider the nonlinear differential equation with several delays x (t) + f t, xh (t),, xh (t) = 0, t 0, (1) where f (t, u,, u ) is a Carathéodory function, f (t, 0,,0) = 0, h (t) t, lim h (t) =. We will assume that these conditions hold for all nonlinear functions and all delays. 130 Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics. 2017, vol. 17, no. 3, pp

3 Березанский Л., Браверман Е. Линеаризованный принцип устойчивости для дифференциальных уравнений с задержками Together with equation (1) we consider for each t 0 0 an initial value problem x (t) + f t, xh (t),, xh (t) = 0, t t, (2) x(t) = φ(t), t < t, x(t ) = x. (3) We also assume that the following hypothesis holds: φ (, t ) R is a Borel measurable bounded function. Definition. A locally absolutely continuous function x R R is called a global solution of problem (2), (3), if it satisfies equation (2) for almost all t (t ; ) and equalities (3) for t t. We assume that a global solution of (2) (3) exists and is unique Linear equations Consider the linear equation with several delays x (t) + a (t)x(g (t)) = 0, (4) which is a partial case of equation (1). We will assume that a k (t) and delays g k (t) t are Lebesgue measurable essentially bounded on [0; ) functions. The following stability definition is for nonlinear equation (2); it also is applies to linear equation (4). Definition [1]. Equation (2) is stable if for any initial point t and number ε > 0 there exists δ > 0 such that the inequality sup φ(t) + x(t ) < δ implies x(t) < ε, t t, for the solution (2) (3). Equation (2) is asymptotically stable if it is stable and all solutions of (2) (3) for any initial point t tend to zero as t. The following lemmas are modifications of the results of paper [2]. Lemma 1. Suppose a (s) ds <. Then all solutions of equation (4) are bounded. { ()} (5) Lemma 2. Suppose a (t) 0, a (s)ds =, lim sup a (s)ds < 1 +. Then equation (4) is asymptotically stable. Consider now a linear equation in the following form x (t) + a (t)xh (t) + b (t)xg (t) = 0. Lemma 3. Suppose a (t) 0, (6) () a (t)dt =, lim sup a (s)ds <, (7) there exists λ, 0 < λ < 1, such that b (t) < λ a (t). (8) If there exists an eventually positive solution of equation (4) then equation (6) is asymptotically stable. Consider a partial case of equation (6) with a nondelay term x (t) + a(t)x(t) + b (t)xg (t) = 0. (9) Lemma 4. Suppose a(t) 0, () a(t)dt =, lim sup a(s)ds <, (10) there exists λ, 0 < λ < 1, such that b (t) < λa(t). Then equation (9) is asymptotically stable Main result Now we can formulate our main result. To this end consider the following condition. There exists a partition {1,2,, m} = I I I, and indices i k, such that (11) Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника» Т. 17, 3. С

4 Краткие сообщения (,,, ) β (t), u 0, k I, (12) α (t) (,,, ), u 0, k I, (13) α (t) (,,, ) β (t), u 0, k I. (14) Theorem Suppose there exist essentially bounded functions α (t), β (t) such that conditions (12) (14) hold. Suppose also that all solutions of (4) are bounded for any a k (t) such that a (t) β (t), k I ; a (t) α (t), k I ; α (t) a (t) β (t), k I. (15) Then all solutions of equation (1) are bounded. 2. Suppose there exist essentially bounded functions α (t), β (t) and numbers A, B, A < B such that conditions (12) (14) hold for all u k satisfying the inequality A u k B. Suppose also that equation (4) is asymptotically stable for any a k (t) for which condition (15) holds. Then lim x(t) = 0 for all solutions x(t) of equation (1) such that A x(t) B. 4. Examples In this section we will apply Theorem 1 to obtain global stability results for some nonlinear delay equations. It means that we will prove that all solutions or all positive solutions tend to an equilibrium as t. Such applications consist of two steps. First by Theorem 1, Part 1 we prove that all solutions of the equation are bounded. Then by the second part of Theorem 1 we prove that every bounded solution tends to zero or to another equilibrium. Example 1. Consider the equation x (t) = (()) (()) ax(t), (16) where 0 b < a, lim supt g(t) <. We have only one nonlinear term. Denote f(t, u, u ) =. We have b (,, ) b. Consider now the linear equation y (t) = b(t)yg(t) ay(t), (17) where b b(t) b < a. There exists λ, 0 < λ < 1, such that b < λa. Then b(t) < λa. By Lemma 4 equation (17) is asymptotically stable. In particular all solutions of this equation are bounded. If we take α(t) = b, β(t) = b, then condition (14) is satisfied. By the first part of Theorem 1 all solution of equation (16) are bounded. If we take again α(t) = b, β(t) = b and two arbitrary numbers A, B, A < B as a lower and an upper bounds, respectively, then condition (14) is satisfied for A u i B; i = 1, 2. Hence by the second part of Theorem 1 for any solution of equation (16) such that A x(t) B we have lim x(t) = 0. Since A and B are arbitrary numbers, then for any solution of equation (16) we have lim x(t) = 0. If equation (16) involves an arbitrary positive degree γ > 0 rather than the square in the denominator x (t) = (()) (()) ax(t), (18) then it is a generalization of the Mackey-Glass equation [4, 5], where g(t) h(t) t τ. Since all solutions of (18) with nonnegative initial conditions and a positive initial value are positive (see, for example, Lemma 1 in [6]), then the same argument as above leads to the conclusion that if 0 b < a, then all solutions (with nonnegative initial conditions) tend to zero. Example 2. Consider the equation x (t) = b1 + x g(t)x(t) + axh(t), (19) where 0 a < b, lim supt h(t) <. Denote f(t, u, u ) = b(1 + u )u. 132 Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics. 2017, vol. 17, no. 3, pp

5 Березанский Л., Браверман Е. We have (,, ) b > 0. Consider now the linear equation Линеаризованный принцип устойчивости для дифференциальных уравнений с задержками y (t) = b(t)y(t) + ayh(t), (20) where b(t) b > 0. There exists λ, 0 < λ < 1, such that a < λb λb(t). By Lemma 4 equation (20) is asymptotically stable. Hence all solutions of this equation are bounded. If we take α(t) = b, then condition (13) is satisfied. By the first part of Theorem 1 all solutions of equation (19) are bounded. If we take again α(t) = b and two arbitrary numbers A, B, A < B as a lower and an upper bounds, respectively, then condition (13) is satisfied for any A u i B, i = 1, 2. Hence by the second part of Theorem 1 for any solution of equation (19) such that A x(t) B we have lim x(t) = 0. Since A and B are arbitrary numbers, then for any solution of equation (19) we have lim x(t) = 0. Example 3. Consider the equation where x (t) = a(t) arctan xh(t), (21) a(t) 0, a(s)ds =, Denote f(t, u) = a(t) arctan u. We have 0 (,) a(t), u 0. Consider now the linear equation lim sup a(s)ds < 1 + (). (22) y (t) = b(t)yh(t), (23) where 0 b(t) a(t). By Lemmas 1 and 2 all solutions of equation (23) are bounded. If we take α(t) = 0, β(t) = a(t), then condition (14) is satisfied. By the first part of Theorem 1 all solutions of equation (21) are bounded. Suppose x(t) is a solution of (21). Then x(t) H for some H > 0. There exists α > 0 such that α, u H. Hence 0 αa(t) (,) a(t), u H. Consider again linear equation (23), where αa(t) b(t) a(t). By Lemma 2 equation (23) is asymptotically stable. If we take α(t) = αa(t), β(t) = a(t), and as lower and an upper bounds two numbers A = H; B = H, then condition (14) is satisfied for any -H u i H, i = 1, 2. Hence by the second part of Theorem 1 for any solution of equation (21) such that -H x(t) H we have lim x(t) = 0. Since x(t) is an arbitrary solution, then for any solution of equation (21) we have lim x(t) = 0. References 1. Hale J.K., Verduyn Lunel S.M. Introduction to Functional Differential equations. Applied Mathematical Sciences, vol. 99. Springer-Verlag, New York, 1993, pp DOI: / Berezansky L., Braverman E. New Stability Conditions for Linear Differential Equations with Several Delays, arxiv: v1 [math.ds], June 20, p. DOI: /2011/ So J.W.H., Yu J.S., Chen M.P. Asymptotic Stability for Scalar Delay Differential Equations, Funkcial. Ekvac, 1996, vol. 39, pp Mackey M.C., Glass L. Oscillation and Chaos in Physiological Control Systems, Science, 1977, vol. 197, pp DOI: /science Losson J., Mackey M.C., Longtin A. Solution Multistability in First order Nonlinear Differential Delay Equations, Chaos, 1993, vol. 3, no. 2, pp DOI: / Berezansky L., Braverman E. Mackey-Glass Equation with Variable Coefficients, Comput. Math. Appl., 2006, vol. 51, iss. 1, pp Received 5 May 2017 Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника» Т. 17, 3. С

6 Краткие сообщения УДК DOI: /ctcr ЛИНЕАРИЗОВАННЫЙ ПРИНЦИП УСТОЙЧИВОСТИ ДЛЯ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ С ЗАДЕРЖКАМИ Л. Березанский 1, Е. Браверман 2 1 Университет им. Бен-Гуриона в Негеве, Беэр-Шева, Израиль, 2 Университет Калгари, Калгари, Канада Предложен линеаризованный принцип глобальной устойчивости для нелинейной задержки в дифференциальных уравнениях, которые иллюстрируются несколькими моделями динамики населения. Приводится обзор некоторых математических моделей с возможным применением линеаризованного принципа. Ключевые слова: дифференциальные уравнения с задержками, линеаризованный принцип глобальной устойчивости Березанский Леонид, отделение математики, Университет им. Бен-Гуриона в Негеве, Беэр- Шева, Израиль; brznsky@cs.bgu.ac.il. Браверман Елена, кафедра математики и статистики, Университет Калгари, Калгари, Канада. Поступила в редакцию 5 мая 2017 г. ОБРАЗЕЦ ЦИТИРОВАНИЯ Berezansky, L. Linearized Stability Principle for Differential Equations with Delays / L. Berezansky, E. Braverman // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника» Т. 17, 3. С DOI: /ctcr FOR CITATION Berezansky L., Braverman E. Linearized Stability Principle for Differential Equations with Delays. Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics, 2017, vol. 17, no. 3, pp DOI: /ctcr Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics. 2017, vol. 17, no. 3, pp