Access the full text.
Sign up today, get DeepDyve free for 14 days.
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
where a, b are real numbers, was studied by Kuruklis [8] and Papanicolaou [15]. The following assertion describes the boundaries of the stability domain for (1.2) in the (a,b) plane (see also Figure 1.1). Theorem 1.1 (see [7, Theorem 2]). The zero solution of (1.2) with k > 1 is asymptotically stable if and only if the pair (a,b) is the internal point of the ï¬nite domain bounded by the Hindawi Publishing Corporation Advances in Diï¬erence Equations Volume 2006, Article ID 31409, Pages 1â9 DOI 10.1155/ADE/2006/31409 b 1 1 a b 1 1 a  1  1 (a)  1  1 (b) Figure 1.1. The stability domains of (1.2). (a) k is odd, k > 1. (b) k is even. The Cohn domain boundary |a| + |b| = 1 are shown both in (a), (b). following lines: (I) a + b = 1, (II) a = â sin kÏ/ sin(k â 1)Ï, b = (â1)k+1 sinÏ/ sin(k â 1)Ï, (III) (â1)k b â a = 1, (IV) a = sinkÏ/ sin(k â 1)Ï, b = â sinÏ/ sin(k â 1)Ï, where Ï varies between 0 and Ï/k. It follows from Cohnâs results [3] that (1.2) is asymptotically stable provided that |a| + |b| < 1. (1.3) A special case of (1.2), namely, the scalar equation xn = xnâ1 + bxnâk , was studied by Levin and May [9]. They established that the equation is asymptotically stable if and only if 2sin(Ï/2(2k â 1)) > âb > 0. Scalar equation with two delays xn = axnâm + bxnâk , more compound than (1.2), was investigated by Dannan [5], Kipnis and Nigmatulin [7], and Nikolaev [13, 14]. The matrix equation xn = xnâ1 + Bxnâk , (1.4) where B is an (m à m) real matrix, xn : N â Rm , a positive integer k is a delay, was investigated by Levitskaya [10]. She transferred Rekhlitskiiâs result [16] about the diï¬erential Ë equation x = Ax(t â Ï), where a positive Ï is a delay, A is a real matrix, to the diï¬erence equations. She established that (1.4) is asymptotically stable if and only if any eigenvalue of the matrix B lies inside the oval of the complex plane bounded by a curve Î = z â C : z = 2isin Ï iÏ Ï Ï . e , â â¤Ï⤠2k â 1 2 2 (1.5) Our ï¬rst problem is to obtain a characteristic equation for (1.1). Second problem is to obtain the necessary and suï¬cient condition in terms of the eigenvalues location of the matrix B for the asymptotic stability of the equation xn = âxnâ1 + Bxnâk , (1.6) M. Kipnis and D. Komissarova 3 where B is a real (m à m) matrix. Further, we give a suï¬cient condition for the asymptotic stability of (1.1) similar with the Cohn condition (1.3). Finally, we give explicit suï¬cient stability conditions for (1.1) in cases A I and A âI, that is, when (1.1) looks like (1.4) and (1.6), respectively. 2. Characteristic equation of (1.1) The following theorem deals with the stability of the linear system: xn = Ai xnâi , (2.1) where Ai are (m à m) real matrixes (i = 1,2,...,k), xn : N â Rm . The equations of type (2.1) in a Banach space were investigated in [1]. Theorem 2.1. Equation (2.1) is asymptotically stable if and only if all roots of the equation det Izk â k i=1 Ai zkâi = 0 (2.2) lie inside the unit disk. If at least one root of (2.2) lies outside of the unit disk, then (2.1) is unstable. Proof. Consider a generating function x(z) = â 0 xn zn of the sequence xn (n ⥠0). For n= the solution of (2.1), the following equality holds: x(z) â â Ai â i â1 j =0 â x j âi z j + x(z)zi â = 0. (2.3) Here xâ1 ,...,xâk are the initial conditions. From (2.3), we have (I â k=1 Ai zi )x(z) = i F(z), where F(z) is an (m à 1) matrix whose elements are polynomials. Hence, x(z) = 1 det I â k i i=1 Ai z Q(z), (2.4) where Q(z) is another (m à 1) matrix whose elements are polynomials. If all roots of (2.2) lie inside the unit disk, then all roots of the denominator of (2.4) lie outside the unit disk, and consequently, x(z) may be expanded in a power series with the radius of convergence greater than 1. Hence, xn tends to 0 exponentially with n â â regardless of the initial conditions. It follows that (2.1) is asymptotically stable. Let there exist a root z of (2.2), such that |z| ⥠1. We will seek for the solution of (2.1) in the form xn = Dzn , where D â Cm , D = 0. Then we obtain (Izk â k=1 Ai zkâi )D = 0. i This system is degenerate. Hence, there exists nonzero solution D, and xn = Dzn does not tend to zero as n â â. Both real sequences xn = ReDzn and xn = ImDzn are solutions of (2.1), and at least one of them does not tend to zero. Hence, (2.1) is not asymptotically stable. If, in addition, |z| > 1, then the solution xn = Dzn is unbounded. At least one of the real sequences xn = ReDzn or xn = ImDzn is unbounded, and (2.1) is unstable. Corollary 2.2. Equation (1.1) is asymptotically stable if and only if all roots of the equation det B + Azkâ1 â Izk = 0 (2.5) lie inside the unit disk. If at least one root of (2.5) lies outside the unit disk, then (1.1) is unstable. A question arises: can we formulate a necessary and suï¬cient condition for the asymptotic stability of (1.1) in terms of restrictions on the eigenvalues of matrixes A and B? The answer is no, as indicated by the following example. 0.3 Example 2.3. Consider the equation xn = Axnâ1 + Bxnâ3 , where B = ( â0.1 0.2 ). In case A = 0.1 0.1 â0.6 ( â0.7 0.2 ), the equation is asymptotically stable by the Corollary 2.2. However, if A = 0.6 ( â0.2 0.5 ), then it is unstable, although eigenvalues of the matrix A are λ1 = 0.8, λ2 = â0.5 0.5 in both cases. If k > 1 and |b| > 1, then the scalar equation (1.2) is unstable (Figure 1.1). We obtain a similar result for the matrix equations (1.1) and (2.1). Theorem 2.4. If | det Ak | > 1, then (2.1) is unstable. Proof. The characteristic equation (2.2) has the form det Izk â Ai zkâi â¡ P(z) â¡ km i =0 ai zi = 0, (2.6) where akm = 1. We put z = 0 in (2.6) and obtain |a0 | = | detAk | > 1. Hence, there exists a root of the polynomial P(z) outside of the unit disk, since |a0 | is the modulus of product of all roots of P(z). Corollary 2.5. If k > 1 and | detB | > 1, then (1.1) is unstable. 3. Stability ovals for (1.6) Theorem 3.1. The system (1.6) is asymptotically stable if and only if all eigenvalues of the matrix B lie inside the region of the complex plane bounded by the curve Î = z â C : z = (â1)k 2isin Proof. We ï¬rst consider the scalar equation xn = âxnâ1 + λxnâk , where λ â C. The characteristic equation for (3.2) is zk + zkâ1 = λ. (3.3) (3.2) Ï iÏ Ï Ï . e , â â¤Ï⤠2k â 1 2 2 (3.1) M. Kipnis and D. Komissarova 5 Writing z = âv in (3.3) gives vk â vkâ1 = (â1)k λ. But all roots of the equation z k â z k â1 = λ (3.5) (3.4) lie inside the unit disk if and only if λ lies inside the region of the complex plane bounded by the curve (1.5). A proof of this fact had been given in [10, Theorem 1]. Therefore, all roots of (3.4), as well as (3.3), lie inside the unit disk if and only if λ lies inside the oval (3.1). Return to the matrix equation (1.6). It follows from Corollary 2.2 that the characteristic equation for the system (1.6) is det B â I zk + zkâ1 = 0. (3.6) Let λ1 ,λ2 ,...,λm be the eigenvalues of B. Consider the equation (cf. (3.3)) z k + z k â1 = λ i , i = 1,2,...,m. (3.7) If each root of any equation (3.7) lies inside the unit disk, then each solution of (3.6) lies inside the unit disk, and conversely. The inclusion of all solutions of (3.7) in the unit disk is the condition for number λi to belong to the interior of oval (3.1). Thus, system (1.6) is asymptotically stable if and only if all eigenvalues of B lie inside the oval bounded by the curve (3.1). Stability ovals for (1.6) are displayed in Figure 3.1. Projection of the results of Theorem 3.1 onto the real axis gives the following addition to the classical result of Levin and May [9], mentioned in the introduction. Corollary 3.2. The scalar equation (3.2) with λ â R is asymptotically stable if and only if 2sin Ï > (â1)k+1 λ > 0. 2(2k â 1) (3.8) Remark 3.3. If k = 1, then the stability oval is the disk of radius 1 centered at (1 + 0i). If k is odd in (1.6), then stability oval is located in the right half-plane, if k is even, then it is located in the left half-plane. Therefore, it is necessary for asymptotic stability of (1.6) that the condition Re λi > 0 with k odd and the condition Re λi < 0 with k even hold. Here λi (1 ⤠i ⤠m) are the eigenvalues of B. Compare it with scalar equation (1.2) and its stability domain in Figure 1.1. We see that with a = â1 in (1.2), the condition b > 0 is necessary for asymptotic stability with k odd, while b < 0 is necessary with k even. The similarity is evident. Corollary 3.4. If there exist eigenvalues λ1 , λ2 of B, such that Re λ1 > 0, Reλ2 < 0, then (1.6) is unstable with any k. 0.4 0.2 0 k=2 k=4 k=6 k=5 k=3  0.2  0.4  1  0.8  0.6  0.4  0.2 Figure 3.1. Stability ovals for the system (1.6). The arrows represent the eigenvalues of a matrix B in Example 3.6. Now we are in a position to strengthen Corollary 2.5 for (1.4) and (1.6). Theorem 3.5. If the system (1.4) (the system (1.6)) is asymptotically stable, then all eigenvalues of B lie inside the unit disk. Proof. It is suï¬cient to consider the stability ovals (1.5) and (3.1) and to remark that |2sin(Ï/2(2k â 1))| ⤠1 for k > 1. 0.5 1.3 Example 3.6. Consider (1.6) with B = ( â0.5 â1.1 ). Eigenvalues of the matrix B are λ1,2 = â0.3 ± 0.1i. Conï¬guration of the λ1 , λ2 , and the stability ovals (see Figure 3.1) indicates that the equation is asymptotically stable for k = 2,4 and unstable for k = 1,3, and for any k > 4. 4. Explicit stability conditions for (1.1) In what follows, · is any matrix norm which satisï¬es the following conditions: (I) A ⥠0, and A = 0 if and only if A = 0, (II) for each c â R, cA = |c| · A , (III) A + B ⤠A + B , (IV) AB ⤠A · B for all (m à m) matrices A, B. In addition, matrix norm should be concordant with the vector norm · â , that is, Ax â ⤠A · x â (4.1) for all x â Rm and any (m à m) matrix A. For real (m à m) matrix A, we deï¬ne, as usual, A 1 = max1⤠j â¤m m 1 |ai j | and A â = i= max1â¤iâ¤m m 1 |ai j |. j= We will give a Cohn-type suï¬cient stability condition [3] for (2.1) (see also [2, Theorem 2.1] and [11, Theorem 2]). Theorem 4.1. If Ai < 1, then (2.1) is asymptotically stable. M. Kipnis and D. Komissarova 7 Proof. Assume that max( xâ1 ,..., xâk ) = M and n = 0,1,...,k â 1, we have xn = k i=1 Ai = b < 1. Sequentially for Ai xnâi ⤠k i=1 Ai M = bM < M. (4.2) Then for n = k, k + 1,...,2k â 1, we derive sequentially xn = Ai xnâi ⤠Ai bM = b2 M < bM. (4.3) Furthermore, for n ⥠rk, r â N, we deduce similarly that xn ⤠br+1 M. Next Corollary gives a delay-independent Cohn-type stability condition for (1.1). Corollary 4.2. If A + B < 1, then (1.1) is asymptotically stable. 0 0.3 Example 4.3. Consider (1.1) with A = ( â0.5 0.1 ), B = ( â0.2 â0.1 ). We have A 1 + B 1 = 0.2 0.4 1.4 > 1, but A â + B â = 0.9 < 1. By Corollary 4.2, the equation is asymptotically stable for any delay k. (4.4) Some additional domains to the Cohn domain for scalar variant of (2.1) were described in [2, Theorem 2.3, Corollary 2.4] employing the Halanay-type inequalities. In the spirit of [2] and earlier work [4], we will give a suï¬cient stability condition of (1.1), additional to Cohn-type condition (4.4). The following result is convenient for applica1 and B 1. tion to (1.1) when A â I Theorem 4.4. If A + B + (k â 1) B then (1.1) is asymptotically stable. Proof. Let us rewrite (1.1), xn xnâ1 â xnâk k â1 i =1 k â1 i =1 AâI + B < 1, (4.5) x n âi â x n âi â1 (A â I)xnâiâ1 + Bxnâkâi . (4.6) Combining Theorem 4.1 with (4.6), we obtain the conclusion of the theorem. 0 Example 4.5. Consider (1.1) with A = ( 1.01 â0.02 ) and B = ( â0.2 â0.2 ). It is easily seen 0 1.01 0.01 that A + B > 1 for any norm · . The Cohn-type condition (4.4) cannot be applied. However, A + B 1 + (k â 1) B 1 ( A â I 1 + B 1 ) = 0.83 + (k â 1)0.0504. Hence, Theorem 4.4 guarantees asymptotic stability of (1.1) for k = 1,2,3,4. Additional calculations, based on Corollary 2.2, give instability for k > 4. A lot of explicit stability conditions for scalar equations (see, e.g., [2, 11, 4, 12]) can be translated to the language of linear systems (1.1). However, we will indicate a new phenomenon. We introduce an additional stability condition for (1.1) depending on whether the delay k is odd or even. The condition is convenient for applications to (1.1) when A+I 1 and B 1. Theorem 4.6. If A + (â1)k+1 B + (k â 1) B then (1.1) is asymptotically stable. Proof. Assume k is odd. Let us rewrite (1.1) as xn xnâ1 â xnâk k â1 i=1 A+I + B < 1, (4.7) (â1)i+1 (A + I)xnâiâ1 + Bxnâkâi . (4.8) The conclusion of the theorem is a consequence of Theorem 4.1 and (4.8). Now suppose that k is even. We rewrite (1.1) as xn = (A â B)xnâ1 + B xnâ1 + xnâk = (A â B)xnâ1 + B k â1 i=1 (â1)i+1 (A + I)xnâiâ1 + Bxnâkâi . (4.9) The proof is completed by the same arguments as in the preceding case. 0 Example 4.7. Consider (1.1) with A = ( â1.01 â0.02 ) and B = ( â0.2 â0.2 ) (cf. Example 4.5). 0 â1.01 0.01 The Cohn-type condition (4.4) cannot be applied. For odd values of k, Theorem 4.6 is useless: A + B > 1 for any norm · , hence (4.7) fails when k is odd. For even values of k, we have A â B 1 + (k â 1) B 1 ( A + I 1 + B 1 ) = 0.83 + (k â 1)0.0504. Theorem 4.6 now implies asymptotic stability of the equation for k = 2 and k = 4. Additional calculations of the roots of characteristic equation (2.5) give instability for k = 1, k = 3, and for any k > 4. Acknowledgments This work was partially supported by the Russian Foundation for Basic Research, Grants 04-01-96069 and 07-01-96065. The authors thank V. Karachik, A. Makarov, S. Pinchuk, L. Pesin, and D. Scheglov for their helpful discussions. M. Kipnis and D. Komissarova 9 References [1] L. Berezansky and E. Braverman, On exponential dichotomy, Bohl-Perron type theorems and stability of diï¬erence equations, Journal of Mathematical Analysis and Applications 304 (2005), no. 2, 511â530. [2] L. Berezansky, E. Braverman, and E. Liz, Suï¬cient conditions for the global stability of nonautonomous higher order diï¬erence equations, Journal of Diï¬erence Equations and Applications 11 (2005), no. 9, 785â798. ¨ [3] A. Cohn, Uber die Anzahl der Wurzeln einer algebraischen Gleichung in einem Kreise, Mathematische Zeitschrift 14 (1922), no. 1, 110â148. [4] K. L. Cooke and I. Gy¨ ri, Numerical approximation of the solutions of delay diï¬erential equao tions on an inï¬nite interval using piecewise constant arguments, Computers & Mathematics with Applications 28 (1994), no. 1â3, 81â92. [5] F. M. Dannan, The asymptotic stability of x(n + k) + ax(n) + bx(n â l) = 0, Journal of Diï¬erence Equations and Applications 10 (2004), no. 6, 589â599. [6] J. Dibl´k and D. Ya. Khusainov, Representation of solutions of discrete delayed system x(k + 1) = ı Ax(k) + Bx(k â m) + f (k) with commutative matrices, Journal of Mathematical Analysis and Applications 318 (2006), no. 1, 63â76. [7] M. Kipnis and R. M. Nigmatulin, Stability of trinomial linear diï¬erence equations with two delays, Automation and Remote Control 65 (2004), no. 11, 1710â1723. [8] S. A. Kuruklis, The asymptotic stability of xn+1 â axn + bxnâk = 0, Journal of Mathematical Analysis and Applications 188 (1994), no. 3, 719â731. [9] S. A. Levin and R. M. May, A note on diï¬erence-delay equations, Theoretical Population Biology 9 (1976), no. 2, 178â187. [10] I. S. Levitskaya, A note on the stability oval for xn+1 = xn + Axnâk , Journal of Diï¬erence Equations and Applications 11 (2005), no. 8, 701â705. [11] E. Liz and J. B. Ferreiro, A note on the global stability of generalized diï¬erence equations, Applied Mathematics Letters 15 (2002), no. 6, 655â659. [12] E. Liz and M. Pituk, Asymptotic estimates and exponential stability for higher-order monotone diï¬erence equations, Advances in Diï¬erence Equations 2005 (2005), no. 1, 41â55. [13] Yu. P. Nikolaev, The set of stable polynomials of linear discrete systems: its geometry, Automation and Remote Control 63 (2002), no. 7, 1080â1088. , The geometry of D-decomposition of a two-dimensional plane of arbitrary coeï¬cients of [14] the characteristic polynomial of a discrete system, Automation and Remote Control 65 (2004), no. 12, 1904â1914. [15] V. G. Papanicolaou, On the asymptotic stability of a class of linear diï¬erence equations, Mathematics Magazine 69 (1996), no. 1, 34â43. [16] Z. I. Rekhlitskii, On the stability of solutions of certain linear diï¬erential equations with a lagging argument in the Banach space, Doklady Akademii Nauk SSSR 111 (1956), 770â773 (Russian). Mikhail Kipnis: Department of Mathematics, Chelyabinsk State Pedagogical University, 69 Lenin Avenue, Chelyabinsk 454080, Russia E-mail address: mkipnis@cspu.ru Darya Komissarova: Department of Mathematics, Southern Ural State University, 76 Lenin Avenue, Chelyabinsk 454080, Russia E-mail address: dasha@math.susu.ac.ru
Advances in Difference Equations – Hindawi Publishing Corporation
Published: Aug 2, 2006
You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!
Read and print from thousands of top scholarly journals.
Already have an account? Log in
Bookmark this article. You can see your Bookmarks on your DeepDyve Library.
To save an article, log in first, or sign up for a DeepDyve account if you don’t already have one.
Copy and paste the desired citation format or use the link below to download a file formatted for EndNote
Access the full text.
Sign up today, get DeepDyve free for 14 days.
All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.