ABSTRACT
In 44. S. Elaydi. Stability of Volterra difference equations of convolution type, Dynamical Systems. Nankai Ser. Pure Appl. Math. Theoret. Phys., 4 (1993), 66-72., S. Elaydi obtained a characterization of the stability of the null solution of the Volterra difference equation
by localizing the roots of its characteristic equation
The assumption that (an ) ∈ ℓ1 was the single hypothesis considered for the validity of that characterization, which is an insufficient condition if the ratio R of convergence of the power series of the previous equation equals one. In fact, when R = 1, this characterization conflicts with a result obtained by Erdo¨ s et al. in 88. P. Erdös, W. Feller & H. Pollard. A property of power series with positive coefficients. Bull. Amer. Math. Soc., 55 (1949), 201-204.. Here, we analyze the R = 1 case and show that some parts of that characterization still hold. Furthermore, studies on stability for the R < 1 case are presented. Finally, we study some results related to stability via finite approximation.
Keywords:
difference equation; stability; convolution
RESUMO
In 44. S. Elaydi. Stability of Volterra difference equations of convolution type, Dynamical Systems. Nankai Ser. Pure Appl. Math. Theoret. Phys., 4 (1993), 66-72., S. Elaydi obteve uma caracterização da estabilidade da solução nula da equação a diferenças de Volterra
localizando as raízes de sua equação característica
A suposição de que (a n ) ∈ ℓ1 foi a única hipótese considerada para a validade daquela caracterização, que é uma condição insuficiente se o raio R de convergência da série de potência da equação anterior é igual a um. De fato, quando R = 1, esta caracterização entra em conflito com um resultado obtido por Erdös e colaboradores em 88. P. Erdös, W. Feller & H. Pollard. A property of power series with positive coefficients. Bull. Amer. Math. Soc., 55 (1949), 201-204.. Aqui, nós analisamos o caso R = 1 e mostramos que algumas partes daquela caracterização ainda se mantêm. Ainda, são apresentados estudos sobre a estabilidade para o caso R < 1. Finalmente, estudamos alguns resultados relativos a estabilidade via aproximações finitas.
Palavras-chave:
equação a diferenças; estabilidade; convolução
1 INTRODUCTION
In the present work, we analyze the stability of the null solution of Volterra difference equations of convolution type,
whose recursive process starts at x 0 ∈ ℝ. Several results related to this subject matter circulates in the specialized scientific literature. One of most well-known is the following theorem:
Theorem 1 (See1010. V.B. Kolmanovskii, E.C. Velasco & J.A.T. Muñoz. A survey: stability and boundedness of Volterra difference equations. Nonlinear Analysis, 53 (2003), 861-928.).If|an | < 1, then the null solution of(1.1)is asymptotically stable.
We now present another important characterization of the stability of the null solution of (1.1) as obtained by S. Elaydi in 44. S. Elaydi. Stability of Volterra difference equations of convolution type, Dynamical Systems. Nankai Ser. Pure Appl. Math. Theoret. Phys., 4 (1993), 66-72.. Let (x n ) be a solution of (1.1) with initial condition x 0 = 1.
Consider the two power series
Then, formally, such series satisfies
Thus the coefficients of x(z) can be found by determining the coefficients of the power series representation of the function (1 − a(z))−1. Hence the roots of the characteristic equation
play an important role in this sense. By this reasoning, S. Elaydi obtained necessary and sufficient conditions for the stability of the null solution by localizing the roots of 1 − a(1/z) = 0 with respect to the set {z ∈ ℂ : ║z║ ≥ 1}. We now enunciate the result obtained by Elaydi with a variable change by writing z in place of 1/z (as considered in 44. S. Elaydi. Stability of Volterra difference equations of convolution type, Dynamical Systems. Nankai Ser. Pure Appl. Math. Theoret. Phys., 4 (1993), 66-72.). In the following, we use the notation:
Theorem 2 (See44. S. Elaydi. Stability of Volterra difference equations of convolution type, Dynamical Systems. Nankai Ser. Pure Appl. Math. Theoret. Phys., 4 (1993), 66-72.), (66. S. Elaydi. Stability and asymptoticity of Volterra difference equations. A progress report. J. Comp. Appl. Math., 228 (2009) 504-513.).Let (a n ) ∈ ℓ1. Then:
-
(a) The null solution of(1.1)is stable if, and only if, the characteristicequation (1.3)has no roots in B 1(0) and its possible roots in |z| = 1 are of order 1.
-
(b) The null solution of(1.1)is asymptotically stable if, and only if, the characteristicequation (1.3)has no roots in .
It is worth mentioning that the previous theorem has also appeared as theorems 6.16 and 6.17 in 55. S. Elaydi. An Introduction to Difference Equations, Undergraduate Texts in Mathematics, Springer Verlag, (2005).. Furthermore, as a consequence of theorem 2, Elaydi set the following result on asymptotic instability:
Theorem 3 (See44. S. Elaydi. Stability of Volterra difference equations of convolution type, Dynamical Systems. Nankai Ser. Pure Appl. Math. Theoret. Phys., 4 (1993), 66-72.), (66. S. Elaydi. Stability and asymptoticity of Volterra difference equations. A progress report. J. Comp. Appl. Math., 228 (2009) 504-513.).If (a n ) ∈ ℓ1 is a sequence whose terms do not change signs for n ≥ 1, then the null solution of(1.1)is not asymptotically stable if one of the following conditions is satisfied:
-
(a) a n ≥ 1;
-
(b) a n ≤ −1 and a n > 0 for some n ≥ 1;
-
(c) a n ≤ −1 and a n < 0 for some n ≥ 1 anda n is sufficiently small.
At this point, if we consider the sequence
the null solution of (1.1) is not asymptotically stable by the item (b) of theorem 2 or the item (a) of theorem 3. On the other hand, (1.4) satisfies the conditions for asymptotic stability of the null solution of (1.1) as given by the following theorem due to Erdo¨ s, Feller e Pollard:
Theorem 4 (See88. P. Erdös, W. Feller & H. Pollard. A property of power series with positive coefficients. Bull. Amer. Math. Soc., 55 (1949), 201-204.).Let (a n ) be a sequence of nonnegative terms such that
Then the null solution of (1.1) is asymptotically stable.
Therefore there exists a contradiction between the previous theorem and theorems 2 and 3. An analysis of the proof of theorem 2 makes clear that the analyticity of the power series a(z) on the circumference |z| = 1 was strongly used. But this fact is not a consequence of the assumption that (a n ) ∈ ℓ1 , as we can see in the example (1.4). Hence a simple correction can be made by introducing the radius of convergence of the series a(z),
and replacing the hypothesis that (a n ) ∈ ℓ1 by R > 1. In fact, if R > 1, then the function a(z) is analytic in |z| = 1 and (a n ) ∈ ℓ1, which are conditions that assure us of the validity of theorem 2.
Furthermore, by applying this new hypothesis, there is no contradiction between theorem 4 and theorems 2 and 3 since we may easily show that the conditions of theorem 4 implies that R = 1. In this sense, in order to avoid conflicts, we rewrite theorem 2 as follows:
Theorem 2’.Replace (a n ) ∈ ℓ1 by R > 1. Then the items (a) and (b) of theorem 2 are valid.
Therefore, if R = 1, the validity of theorem 2 is an open problem since the analyticity on the unit circumference can not be applied. Besides theorem 4, some other results give us some sufficient conditions for the asymptotic stability of the null solution of (1.1) and can be found in 11. K.S. Berenhaut & N.G. Vish. Equations of convolution type with monotone coefficients. Journal of Difference Equations and Applications, 17 (2011), 555-566.), (77. S. Elaydi, E. Messina & A. Vecchio. A note on the asymptotic stability of linear Volterra difference equations of convolution type. J. Difference Equ. Appl., 13 (2007), 1079-1084.), (1111. E. Messina, Y. Muroya, E. Russo & A. Vecchio. Asymptotic Behavior of Solutions for Nonlinear Volterra Discrete Equations, Discrete Dynamics in Nature and Society (2008) Article ID 867623, 18 pp.), (99. V.B. Kolmanovskii & N.P. Kosareva. Stability of Volterra Difference Equations. Differential Equations, 37 (2001), 1773-1782.), (22. S.K. Choi, Y.H. Goo, Y. Hoe & N.J. Koo. Asymptotic behavior of nonlinear Volterra difference systems. Bull. Korean Math. Soc., 44 (2007), 177-184.), (1414. X.H. Tang & Z. Jiang. Asymptotic behavior of Volterra difference equation, J. Difference Equ. Appl., 13 (2007), 25-40.), (1212. N. Van Minh. On the asymptotic behaviour of Volterra difference equations, J. Difference Equ. Appl., 19 (2013), 1317-1330.), (1313. R. Nigmatulin. Asymptotic behavior of solutions of a nonlinear Volterra difference equation, Int. Electron. J. Pure Appl. Math., 6 (2013), 123-125..
Finally, note that the stability of the null solution of (1.1) depends on the behavior of the particular solution (x n ) with initial condition x 0 = 1. In fact, an arbitrary solution to the equation (1.1) with initial condition β is given by (x n β). Thus, in what follows, (x n ) denotes the solution of (1.1) with initial condition x 0 = 1. Therefore, just for future reference, we have the following elementar result, which is already known in a more general case (see e.g. theorem 3.1 in 33. M.A. Crisci, V.B. Kolmanovskii, E. Russo & A. Vecchio. Boundedness of Discrete Volterra Equations. Journal of Mathematical Analysis and Applications, 211 (1997), 106-130.):
Theorem 5. The null solution of (1.1) is:
-
stable if, and only if, (x n ) is bounded;
-
asymptotically stable if, and only if, xn → 0.
The rest of this paper is divided as follows: In section 2 we analyze the stability for R = 1. Furthermore, we study possible characterizations for the null solution of (1.1) to be stable if R < 1. In section 3 we analyze the stability via finite approximations.
2 RESULTS ON STABILITY WHEN R = 1 OR R < 1
In this section, we provide some results concerning the stability/instability for the null solution of (1.1) when R = 1 and R < 1. It is worth pointing out that the arguments to be used are still valid
in the case where R > 1. Initially we give two examples on stability/instability for the R < 1 case.
Example. Let p > 1. At first consider the sequence an = − pn . Then R < 1 and the solution of (1.1) with initial condition x0 = 1 is (xn ) = (1, − p, 0, 0, .. .), which converges to zero. So
the null solution is asymptotically stable. On the other hand, if we consider
we also have R < 1. However, since each term in this sequence is positive, it follows (by induction) that each term of (x n ) is positive. Hence x n ≥ a n for each n ≥ 1. Thus (x n ) is unbounded and consequently the null solution of (1.1) is unstable by theorem 5.
The previous example shows that the stability or instability of the null solution does not depend on the radius of convergence of a(z). In what follows we present a result on instability which does not depend on the size of R. (In particular, it holds if R = 1 or R < 1, which is the case where Theorem 2_ is not applicable.)
Theorem 6.Let Da = . If the characteristicequation (1.3)has a root in D a ∩ B 1(0), then the null solution of(1.1)is unstable.
Proof. Denote by ρ 0 e iθ , ρ 0 ∈]0, 1[, one of the zeros of 1 − a(z) of smallest modulus. Hence 1 − a(z) ≠ 0 for each z ∈ B ρ0 (0). Since x(z) = (1 − a(z))−1 for every z ∈ B ρ0 (0), we have that
Therefore the radius of convergence of x(z) is not greater than ρ 0. Then
Thus, if > β > 1, there exists a subsequence () for which > β. So we conclude that . Therefore the sequence (x n ) is not bounded. □
Remark. The converse of theorem 6 is not valid. In fact, for the sequence given in (2.1), the null solution of (1.1) is unstable. On the other hand, since D a = and
it follows that 1 − a(z) is not zero in D a ∪ B 1(0).
Corollary 7.Ifan > 1 converges, then the null solution of(1.1)is unstable.
Proof. Since a n converges, it follows from Abel’s theorem that the power series a(z) is continuous in [0, 1]. Specifically, one has that . So consider the function b(z) = 1 − a(z). Then b(0) = 1 and b(1) < 0. Therefore b(z) has a zero in the interval ]0, 1[⊂ B 1 (0). It follows from the previous theorem that the null solution of (1.1) is not stable. □
Remark. By item (a) of theorem 3, we have that the lack of asymptotic stability takes place when, in particular, the hypothesis of the previous corollary holds, provided that the terms of the sequence (a n ) do not change signs. So, the previous corollary states the lack of stability (and therefore the lack of asymptotic stability) without any sign-preserving condition. Furthermore, that corollary remains valid if we replace the hypothesis .
The following theorem shows that part of what was stated in the item (b) of theorem 2’ remains valid if R = 1.
Theorem 8.If 1 − a(z) is continuous and not zero in, then the null solution of(1.1)is asymptotically stable.
Proof. First we observe that the function x(z) = (1 − a(z))−1 is uniformly continuous on . Now, for each ρ ∈ [0, 1[ and n ∈ ℕ, define the following two functions on the interval [0, 2π]:
It follows from the uniform continuity of x(z) on that f ρ → f uniformly on [0, 2π] as ρ → 1−. Therefore
From the Cauchy’s Integral Formula, we have that, for every ρ ∈]0, 1[,
Applying the limit when ρ → 1−, it follows that
So, by the Riemann-Lebesgue Lemma, one has x n → 0, which shows that the null solution of is asymptotically stable. □
Remark. The conclusion of the preceding theorem is not valid if the radius of convergence of a(z), R, is less than one. In other words, even if 1 − a(z) is continuous and not zero in , the null solution may not be asymptotically stable. To illustrate this statement, it suffices to consider the sequence given in (2.1). Additionally, note that, if R = 1, the converse of the preceding theorem is not valid, as shown in example (1.4).
To finalize this section, we enunciate an auxiliary lemma for the characterization of the stability (not necessarily an asymptotic one) of the null solution when R = 1.
Lemma 9. Consider the following power series
Suppose that (y n ) is bounded. Then:
-
If p(z) = (1 − e −iθ z)y(z), then (p n ) is bounded.
-
If p(z) = (1 − z)y(z), then p(z) is bounded on the interval [0, 1[.
Proof. First note that p n = y n − e −iθ y n−1 for each n ≥ 1. Hence, if C = supn≥0 |y n |, then |p n | ≤ 2C for every n ≥ 1, which demonstrates the item 1. Consider now that the hypothesis of the item 2 is valid. It follows that z ∈ [0, 1[ implies
In what follows, suppose that the power series a(z) converges on and the possible zeros of 1 − a(z) occur at , in other words,
where q(z) is not zero in . Furthermore, consider the space
Next, it follows a characterization of the stability for R = 1:
Theorem 10.Assume that the power series a(z) converges onand q ∈ 𝕃1. The null solution of(1.1)is stable if, and only if, the possible zeros of 1 − a(z) as given in(2.2)are of order 1.
Proof. First consider that m 1 = ··· = m s = 1. Since q ∈ 𝕃1 and q is not zero in , by Wiener’s Theorem, we have that
On the other hand,
Since x(z) = (1 − a(z))−1 for z ∈ B 1(0), it follows from (2.3) and (2.4) that
As a result of equating coefficients, we have that x n = . Therefore
where C is a constant which is an upper bound for the sequence (|α n |). So (x n ) is bounded and then, by theorem 5, the null solution of (1.1) is stable. Conversely, suppose that the null solution of (1.1) is stable and m 1 ≥ 2. Replacing z by in x(z)(1 − a(z)) = 1, it follows that
with
Since, by theorem 5, the sequence (x n ) is bounded, by applying the item 1 several times and finally the item 2 of the preceding lemma, one has that p(z) is bounded on [0, 1[. Therefore
which is absurd. Hence m 1 = 1. □
Corollary 11. Let n|a n | < ∞. If 1 − a(z) is not zero in B 1(0) and has a finite number of zeros of order one in |z| = 1, then the null solution of(1.1)is stable.
Proof. Suppose that z = 1 is a zero of 1 − a(z) and consider
Then q 0 = 1 and, for n ≥ 1, one has that
It is easy to verify that, for n ≥ 1, we have
So
Thus 1 − a(z) = (1 − z)q(z) with q ∈ 𝕃1. On the other hand, if z = e iθ is a zero of 1 − a(z), one has that z = 1 is a zero of 1 − ã(z), ã(z) = a(e iθ z), which satisfies the hypotheses of this corollary. Hence 1 − ã(z) = (1 − z) (z) with ∈ 𝕃1 or, equivalently, (1 − a(z)) = (1 − e −iθ z)q(z), where q(z) = (e −iθ z) and therefore q ∈ 𝕃1. Now consider the set {e iθ j: j = 1, ... , s} consisting of all zeros of 1 − a(z). Then, by partial fractions, we have
where each q j ∈ 𝕃1. So, if q = β j q j ∈ 𝕃1, then
It follows from the preceding theorem that the null solution of (1.1) is stable. □
Example. The sequence , satisfies the hypotheses of the previous corollary. So the null solution of (1.1) is stable. Note that, in this case, theorem 2’ cannot be used for obtaining this result since R = 1.
3 STABILITY VIA APPROXIMATION
In this final section we state some conditions for stability via polynomial approximation by applying the following theorem (known as Rouche´’s Theorem):
Theorem 12. If f and f + h are analytic functions onsuch that
then f and f + h have the same number of zeros in B ρ (z 0).
Now, for each n ∈ ℕ, consider the polynomial
where a 1, .. ., a n are the first n coefficients of the power series a(z). Define
and z 1, ... , z n as the n zeros of p n (z), that is,
In what follows we enunciate some results of stability via finite approximations of the characteristic equation.
Theorem 13. If there exists an index n such that
then the null solution of (1.1) is asymptotically stable.
Proof. For |z| = 1, we have that
So (1 − r n )n ≤ | p n (z)| for |z| = 1. Consider the nth partial sum of 1 − a(z), that is,
Hence, since s n (z) = z n p n (1/z) for z ≠ 0, one has that s n (z) is not zero in and (1−r n )n ≤ |s n (z)| for |z| = 1. Then
By Rouché’s Theorem, 1 − a(z) and s n (z) have the same number of zeros in . Therefore 1 − a(z) is not zero in . It follows from theorem 8 that the null solution of (1.1) is asymptotically stable. □
Example. Let (β n ) be a sequence with β n ∈ {−1, 1}. The sequence
does not satisfy the hypothesis of theorem 1. However, by considering the polynomial p 2(z), we obtain r 2 = 3/4 and, since
the null solution of (1.1) is asymptotically stable.
Example. The sequence
does not satisfy the hypothesis of theorem 1. By a computational calculus, the values of r n and are as shown in the table that follows. Note that the hypothesis of the preceding theorem is satisfied for n = 6. So the null solution of (1.1) is asymptotically stable.
Theorem 14.If there exists an index n such that rn > 1 and
where ρ0is the point that maximizes the function δn: defined by
the null solution of (1.1) is unstable.
Proof. Since δ() = 0, one has < ρ 0 ≤ 1. If |z| = ρ 0, then |1 − zz i | ≥ |1 − ρ 0|z i ||. So the partial sum considered in the preceding theorem satisfies
Therefore, for |z| = ρ 0, we have that
Now, let j with |z j | = r n . Then ∈ B ρ0 (0) and s n () = 0. Hence, by Rouche´’s Theorem, 1 − a(z) has at least a zero in B ρ0 (0) ⊂ B 1(0). Therefore, by theorem 6, the null solution of (1.1) is unstable. □
Remark. Let us put the moduli of the zeros of p n in descending order, say |z 1| ≥ ··· ≥ |z n |. Assume that the hypotheses of the preceding theorem hold. Consider i 0 is the highest index with |z i0 | > 1. Then |z i0 +1| ≤ 1. Hence one has the following estimates:
-
If i 0 = n, then E 1: = |1 − |z n ||n ≤ δ n (1) ≤ δ n (ρ 0).
-
If |z i0 +1| < 1, then E 2: = min{|1 − |z i0 ||n , |1 − |z i0 +1||n } ≤ δ n (1) ≤ δ n (ρ 0).
-
If |z i0 +1 | = 1, then E 3: = |1 − ρ 1|z i0 ||n ≤ δ n (ρ 1) ≤ δ n (ρ 0) with ρ 1 = 2/(|z i0 |+ 1).
As a consequence of the previous remark, we may state the following corollary:
Corollary 15. With the same assumptions of the preceding theorem, if
where E = E1, E2 or E3 are given as in the previous remark, then the null solution of(1.1)is unstable.
Example. Consider the sequence (a n ) given by
(a n ) does not satisfy the assumptions of corollary 7. However the zeros of p 2(z) are z 1 = z 2 = 2. Since
the null solution of (1.1) is unstable by the previous corollary.
REFERENCES
-
1K.S. Berenhaut & N.G. Vish. Equations of convolution type with monotone coefficients. Journal of Difference Equations and Applications, 17 (2011), 555-566.
-
2S.K. Choi, Y.H. Goo, Y. Hoe & N.J. Koo. Asymptotic behavior of nonlinear Volterra difference systems. Bull. Korean Math. Soc., 44 (2007), 177-184.
-
3M.A. Crisci, V.B. Kolmanovskii, E. Russo & A. Vecchio. Boundedness of Discrete Volterra Equations. Journal of Mathematical Analysis and Applications, 211 (1997), 106-130.
-
4S. Elaydi. Stability of Volterra difference equations of convolution type, Dynamical Systems. Nankai Ser. Pure Appl. Math. Theoret. Phys., 4 (1993), 66-72.
-
5S. Elaydi. An Introduction to Difference Equations, Undergraduate Texts in Mathematics, Springer Verlag, (2005).
-
6S. Elaydi. Stability and asymptoticity of Volterra difference equations. A progress report. J. Comp. Appl. Math., 228 (2009) 504-513.
-
7S. Elaydi, E. Messina & A. Vecchio. A note on the asymptotic stability of linear Volterra difference equations of convolution type. J. Difference Equ. Appl., 13 (2007), 1079-1084.
-
8P. Erdös, W. Feller & H. Pollard. A property of power series with positive coefficients. Bull. Amer. Math. Soc., 55 (1949), 201-204.
-
9V.B. Kolmanovskii & N.P. Kosareva. Stability of Volterra Difference Equations. Differential Equations, 37 (2001), 1773-1782.
-
10V.B. Kolmanovskii, E.C. Velasco & J.A.T. Muñoz. A survey: stability and boundedness of Volterra difference equations. Nonlinear Analysis, 53 (2003), 861-928.
-
11E. Messina, Y. Muroya, E. Russo & A. Vecchio. Asymptotic Behavior of Solutions for Nonlinear Volterra Discrete Equations, Discrete Dynamics in Nature and Society (2008) Article ID 867623, 18 pp.
-
12N. Van Minh. On the asymptotic behaviour of Volterra difference equations, J. Difference Equ. Appl., 19 (2013), 1317-1330.
-
13R. Nigmatulin. Asymptotic behavior of solutions of a nonlinear Volterra difference equation, Int. Electron. J. Pure Appl. Math., 6 (2013), 123-125.
-
14X.H. Tang & Z. Jiang. Asymptotic behavior of Volterra difference equation, J. Difference Equ. Appl., 13 (2007), 25-40.
Publication Dates
-
Publication in this collection
Dec 2017
History
-
Received
05 Aug 2015 -
Accepted
15 July 2017