Research Article | Open Access

Volume 2016 |Article ID 1487824 | https://doi.org/10.1155/2016/1487824

Michał Niezabitowski, "On the Sequences Realizing Perron and Lyapunov Exponents of Discrete Linear Time-Varying Systems", Mathematical Problems in Engineering, vol. 2016, Article ID 1487824, 7 pages, 2016. https://doi.org/10.1155/2016/1487824

# On the Sequences Realizing Perron and Lyapunov Exponents of Discrete Linear Time-Varying Systems

Accepted12 Oct 2016
Published01 Nov 2016

#### Abstract

We investigate properties of partial exponents (in particular, the Lyapunov and Perron exponents) of discrete time-varying linear systems. In the set of all increasing sequences of natural numbers, we define an equivalence relation with the property that sequences in the same equivalence class have the same partial exponent. We also define certain subclass of all increasing sequences of natural numbers, including all arithmetic sequences, such that all partial exponents are achievable on a sequence from this class. Finally, we show that the Perron and Lyapunov exponents may be approximated by partial exponents achievable on sequences in certain sense similar to geometric sequences.

#### 1. Introduction

Consider a discrete time-varying system:where is a bounded sequence of invertible -by- real matrices such that is bounded. For the coefficient matrices, denote the transition matrices and , where is the identity matrix. For an initial condition , the solution of (1) is denoted by ; that is,If is a sequence of real numbers, then the Perron exponent and the Lyapunov exponent of are defined in the following ways:By , denote the Euclidean norm in and the induced operator norm. For an initial condition , the Perron and the Lyapunov exponents of the solution of system (1) are defined (see ) asIt means that the Perron and Lyapunov exponents of the solution are the Perron and Lyapunov exponents of the sequence , respectively.

To characterize many properties of system (1) characteristics exponents, for example the Lyapunov, Perron, Bohl, general exponents may be used. These quantities describe the different types of stability and trajectories growth rate. For interesting summary on main properties of the Lyapunov, Perron, Bohl, general exponents of the discrete time-varying linear system, and relations between these exponents and different types of stability of the considered system see .

The Lyapunov  and the Perron  exponents are one of the most commonly used numerical characteristics of dynamical systems. They describe, inter alia, such important properties like stability (exponential and Poisson). Numerical calculation of them is related to two main problems. The first one is that they are very sensitive to inaccuracies in the coefficients (they are not even continuous functions of the coefficients; see ). The second problem is that these quantities are defined by the partial limits (upper and lower one), and we do not know in advance what time sequence they are achieved on; therefore, an a priori one would need to look into all increasing sequences of natural numbers [30, 31].

This paper is linked to the second problem. In the paper, we try to describe a smaller class of all growing sequences of natural numbers with the property that the Lyapunov or Perron exponents are achieved on one of the sequences in this class.

The paper is organized in the following way: in the next paragraph, we establish certain properties of partial limits of real sequences, in particular their Perron and Lyapunov exponents. In the third section, containing the main results of the work, the theorems from the second section are applied to obtain properties of the Perron and Lyapunov exponents of the solutions of system (1). The work ends with a paragraph containing conclusions and suggestions for further research.

#### 2. Preliminaries

Denote by the set of all sequences of positive real numbers such that there exist constants ,   (in general depending on the sequence ) such that By we denote the set of all increasing sequences of natural numbers. If is any sequence of real numbers and are such that there exists a finite limit then the number we will be called a partial limit or limit point of sequence and we will say that it is achieved on the sequence

The next theorem shows that each number between the Perron and Lyapunov exponents of the sequence is a partial limit of the sequence .

Theorem 1. For each sequence and each number there exists sequence such that

Proof. If or then the conclusion follows from the properties of the upper and lower limits. Suppose that Let us define sequence in the following way: By the inequality and by the definition of the upper and lower limits, it follows that the definition of is correct, the sequence is increasing, and From the above two inequalities we get Since , therefore there exists a constant , such that By the last two inequalities, we obtain Passing to the upper limit and taking into account that , we have Inequalities (11) and (16) imply thatIt means that the sequence ,   is that one from the theorem’s thesis. The proof is completed.

It is easy to construct an example showing that the theorem is no longer true without the assumption that .

Example 2. Let us define sequence in the following way: It is clear that . Moreover, Therefore, each convergent subsequence of the sequence may have as a limit only or

Theorem 1 may be generalized as follows.

Theorem 3. If is such that then for each sequence and each number there exists sequence such that

Proof. This theorem may be obtained from the general fact from the functional analysis as it was shown in [32, Lemma  7.5]. Repeating the construction from the proof of Theorem 1, we obtain instead of inequality (15) the following one:From this inequality and by assumption (20) the thesis follows. The proof is completed.

Let us now introduce certain relation in the set (Definition 4). It will appear to be an equivalence relation (Theorem 5) and if two sequences belong to the same equivalence class, then corresponding to them subsequences of have the same exponents (Theorem 6).

Definition 4. We say that the sequence is close to the sequence ifThis fact will be denoted in the following way

Theorem 5. The relation is an equivalence relation in the set .

Proof. Reflexivity of the relation is obvious. Suppose that For a , denote by any natural number satisfying the condition We will show that the set is infinite. On the contrary, suppose that it is finite and denote its elements by . Then, there exists infinite set such that for all and certain . For , denote Then, The last two equalities are in contradiction with the facts that is infinite and tends to infinity. Now we show symmetry of the relation . Suppose that but the sequence is not close to the sequence Denote The fact that the sequence is not close to the sequence implies that Let us fix ,   By the definition of upper limit we know that there exists sequence such that It means that or equivalently that for all . For the fixed , the second inequality may hold only for finite many (in the opposite case, after passing to the limit with we obtain ). Moreover, if for certain the first inequality holds, then for all ,   Therefore, for all there exists such that for all ,   By the definition of the relation , the fact that and the definition of the limit it follows that there exists such that for all ,  . The last inequality implies that for all Finally, notice that It means that the inequalities (33) and (35) are in contradiction. Therefore, the relation is in fact symmetric.
Now we show the transitivity of the relation . Suppose that we have three sequences such that and From the fact , we conclude that where is any natural number satisfying the condition Let us fix an arbitrary By the definition of the limit and the equality (36), it follows that there exists such that for all ,   The last inequality implies that From the fact that , it follows that To obtain the last equality nondecreaseness of the sequence is necessary. If this is not the case, then we can choose the nondecreasing subsequence from and the further reasoning lead for it. Denote by any natural number such that Applying the introduced notation and the definition of the limit, we conclude that there exists such that for all ,   The last inequality implies that From inequalities (39) and (43), we get for ,  ; that is,Due to arbitrariness of selection of , the last inequality means that However, since thenthat is, . The proof of transitivity of the relation is finished.

Theorem 6. If ,  ,  , and   and there exists the limit then there exists the limit and the limits are equal.

Proof. Since , then there exists a constant such that For , denote by any natural number satisfying the condition We haveTo obtain the last inequality, we used inequality (51). Sinceand there exists the limit then inequality (53) implies the thesis of the theorem. The proof is completed.

Denote by the greatest integer no greater than . The next theorem describes and by the partial limits of which correspond to time subsequences of the form , where ,  .

Theorem 7. For any sequence , the following equalities hold:

Proof. Let be such that Without loss of generality, for further consideration, we may assume that, for fixed in each interval ,  , there are no more than one element of the sequence For , denote by such a number that Additionally, denote by a function given by Since , then there exists a constant , such that In particular, taking and , we get Using this inequality, the introduced notation, and the definition of upper limit, we have Denotingwe get where From inequalities (63) and (65), we have Passing in the last inequality to upper limit with and taking into account that we get Analogically, passing to the lower limit with , we obtain Equality (56) follows from the equalities (69a) and (69b). In the same way, one can prove (57). The proof is completed.

#### 3. Main Results

Consider a solution of system (1) and denote by a common bound for the sequences and . We have The two above inequalities show that .

Applying Theorems 1, 3, and 6 to the sequence , we get the following result.

Theorem 8. The set of limit points of the sequence is the interval If the sequence satisfies assumption (20), then for any number there exists a sequence such that Moreover, if and , then

Notice that each arithmetic sequence satisfies condition (20). Then, taking in the above theorem or we conclude that the Lyapunov and Perron exponents are achieved at a certain subsequence of any arithmetic sequence. We do not know whether the analogous statement is true for geometric sequences. But, applying Theorem 7 to the sequence , we may formulate the following result.

Theorem 9. For any solution of system (1), we have

#### 4. Conclusions

In the paper for the discrete time-varying linear system, we described the limits points of the sequence . This set is equal to the interval Moreover, we proved that each partial limit of this sequence is achievable on a certain subsequence of any sequence satisfying condition (20), in particular on certain subsequence of any arithmetic sequence (Theorem 8). Finally, we showed that the Perron and Lyapunov exponents may be approximated by subsequences in certain sense similar to geometric sequences (Theorem 9). The objective of future works will be the investigation of the possibility of omitting limits with in equalities (73).

#### Competing Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The research presented here was done as part of the project funded by the National Science Centre in Poland granted according to decision DEC-2012/07/N/ST7/03236. Moreover, the Article Processing Charge was covered by funds of the Silesian University of Technology Rector’s Habilitation Grant no. 02/010/RGH16/0048.

1. A. Czornik, Perturbation Theory for Lyapunov Exponents of Discrete Linear Systems, Wydawnictwa AGH, 2012.
2. A. Babiarz, “On the Lyapunov, Perron, Bohl and general exponents of discrete linear time-varying diagonal systems,” in Proceedings of the 21th International Conference on Methods and Models in Automation and Robotics (MMAR '16), pp. 1045–1050, August 2016. View at: Google Scholar
3. A. Lyapunov, A general task about the stability of motion [Ph.D. thesis], 1892.
4. L. Arnold, H. Crauel, and J. P. Eckmann, Lyapunov Exponents: Proceedings of a Conference held in Oberwolfach, May 28–June 2, 1990, vol. 1486 of Lecture Notes in Mathematics, Springer, Berlin, Germany, 1991. View at: Publisher Site
5. L. Barreira and Y. B. Pesin, Lyapunov Exponents and Smooth Ergodic Theory, vol. 23 of University Lecture Series, American Mathematical Society, Providence, RI, USA, 2002. View at: MathSciNet
6. J. L. Daleckii and M. G. Krein, Stability of Solutions of Differential Equations in Banach Space, American Mathematical Society, Providence, RI, USA, 1974. View at: MathSciNet
7. L. Barreira and C. Valls, “Lyapunov's coefficient of regularity and nonuniform hyperbolicity,” Linear and Multilinear Algebra, vol. 64, no. 6, pp. 1124–1144, 2016. View at: Publisher Site | Google Scholar | MathSciNet
8. L. Barreira, D. Dragicevic, and C. Valls, “Tempered exponential dichotomies and Lyapunov exponents for perturbations,” Communications in Contemporary Mathematics, vol. 18, no. 5, article 1550058, 2016. View at: Publisher Site | Google Scholar | MathSciNet
9. N. D. Cong, T. S. Doan, and S. Siegmund, “On Lyapunov exponents of difference equations with random delay,” Discrete and Continuous Dynamical Systems, Series B, vol. 20, no. 3, pp. 861–874, 2015. View at: Google Scholar | MathSciNet
10. N. V. Kuznetsov, T. A. Alexeeva, and G. A. Leonov, “Invariance of lyapunov exponents and lyapunov dimension for regular and irregular linearizations,” Nonlinear Dynamics, vol. 85, no. 1, pp. 195–201, 2016. View at: Publisher Site | Google Scholar | MathSciNet
11. G. A. Leonov, N. V. Kuznetsov, N. A. Korzhemanova, and D. V. Kusakin, “Lyapunov dimension formula for the global attractor of the Lorenz system,” Communications in Nonlinear Science and Numerical Simulation, vol. 41, pp. 84–103, 2016. View at: Publisher Site | Google Scholar | MathSciNet
12. A. Czornik, A. Nawrat, and M. Niezabitowski, “On the Lyapunov exponents of a class of second-order discrete time linear systems with bounded perturbations,” Dynamical Systems, vol. 28, no. 4, pp. 473–483, 2013. View at: Publisher Site | Google Scholar | MathSciNet
13. L. H. Duc, J. P. Chávez, D. T. Son, and S. Siegmund, “Finite-time Lyapunov exponents and metabolic control coefficients for threshold detection of stimulus—response curves,” Journal of Biological Dynamics, vol. 10, no. 1, pp. 379–394, 2016. View at: Publisher Site | Google Scholar | MathSciNet
14. A. Czornik and M. Niezabitowski, “Lyapunov exponents for systems with unbounded coefficients,” Dynamical Systems, vol. 28, no. 2, pp. 140–153, 2013. View at: Publisher Site | Google Scholar | MathSciNet
15. A. Czornik and P. Jurgas, “Set of possible values of maximal Lyapunov exponents of discrete time-varying linear system,” Automatica, vol. 44, no. 2, pp. 580–583, 2008. View at: Publisher Site | Google Scholar | MathSciNet
16. O. Perron, “Uber stabilitat und asymptotisches verhalten der integrale von differentialgleichungssystemen,” Mathematische Zeitschrift, vol. 29, no. 1, pp. 129–160, 1929. View at: Publisher Site | Google Scholar | MathSciNet
17. O. Perron, “Die Stabilitätsfrage bei Differentialgleichungen,” Mathematische Zeitschrift, vol. 32, no. 1, pp. 703–728, 1930. View at: Publisher Site | Google Scholar | MathSciNet
18. E. A. Barabanov, “Attainability of the bound for the number of lower exponents of a linear differential diagonal system,” Doklady Akademii Nauk BSSR, vol. 26, no. 12, pp. 1069–1072, 1982. View at: Google Scholar | MathSciNet
19. A. G. Gargyants, “On the metric typicality of the Perron exponent on solutions of unbounded systems,” Differential Equations, vol. 50, no. 11, pp. 1557–1558, 2014. View at: Publisher Site | Google Scholar | MathSciNet
20. N. A. Izobov and A. V. Filiptsov, “Computation of the maximal lower perron exponent of a linear system,” Differential Equations, vol. 36, no. 11, pp. 1719–1721, 2000. View at: Publisher Site | Google Scholar
21. A. Czornik, “On the Perron exponents of discrete linear systems,” Linear Algebra and its Applications, vol. 432, no. 1, pp. 394–401, 2010. View at: Publisher Site | Google Scholar | MathSciNet
22. A. Czornik, J. Klamka, and M. Niezabitowski, “On the set of Perron exponents of discrete linear systems,” in Proceedings of the 19th IFAC World Congress on International Federation of Automatic Control (IFAC '14), pp. 11740–11742, Cape Town, South Africa, August 2014. View at: Google Scholar
23. A. Konyukh and M. Niezabitowski, “On the attainability of an estimate of a number of Perron exponents for discrete linear diagonal systems with bounded coefficients,” in Proceedings of the 2015 20th International Conference on Methods and Models in Automation and Robotics (MMAR '15), pp. 994–999, Miedzyzdroje, Poland, August 2015. View at: Publisher Site | Google Scholar
24. A. Czornik, P. Mokry, and A. Nawrat, “On the sigma exponent of discrete linear systems,” IEEE Transactions on Automatic Control, vol. 55, no. 6, pp. 1511–1515, 2010. View at: Publisher Site | Google Scholar | MathSciNet
25. A. Czornik, A. Nawrat, and M. Niezabitowski, “On the stability of Lyapunov exponents of discrete linear systems,” in Proceedings of the 12th European Control Conference (ECC '13), pp. 2210–2213, Zurich, Switzerland, July 2013. View at: Google Scholar
26. A. Czornik and M. Niezabitowski, “On the spectrum of discrete time-varying linear systems,” Nonlinear Analysis. Hybrid Systems. An International Multidisciplinary Journal, vol. 9, pp. 27–41, 2013. View at: Publisher Site | Google Scholar | MathSciNet
27. L. Dieci, R. D. Russell, and E. S. Van Vleck, “On the compuation of Lyapunov exponents for continuous dynamical systems,” SIAM Journal on Numerical Analysis, vol. 34, no. 1, pp. 402–423, 1997. View at: Publisher Site | Google Scholar | MathSciNet
28. E. Barabanov, A. Vaidzelevich, A. Czornik, and M. Niezabitowski, “On the influence of the exponentially decaying perturbations of matrices coefficients on the Lyapunov exponents of the discrete time-varying linear system,” in Proceedings of the 17th International Carpathian Control Conference (ICCC '16), pp. 29–34, High Tatras, Slovakia, May 2016. View at: Publisher Site | Google Scholar
29. E. Barabanov, A. Vaidzelevich, A. Czornik, and M. Niezabitowski, “Exact boundaries of mobility of the Lyapunov, Perron and Bohl exponents of linear difference systems under small perturbations of the coefficient matrices,” in Proceedings of the International Conference of Numerical Analysis and Applied Mathematics (ICNAAM '15), vol. 1738 of AIP Conference Proceedings, Rhodes, Greece, September 2015. View at: Publisher Site | Google Scholar
30. E. A. Barabanov, “Calculation of exponents of solutions of linear differential systems from time geometric progressions,” Differential Equations, vol. 33, no. 12, pp. 1596–1603, 1997, translation from Differentsial'nye Uravneniya vol. 33, no. 12, pp. 1592–1600, 1997 View at: Google Scholar
31. E. K. Makarov, “A representation of partial exponents of solutions to linear differential systems via geometric progressions,” Differential Equations, vol. 32, no. 12, pp. 1705–1706, 1996, Translation from Differentsial'nye Uravneniya, vol. 32, no. 12, pp. 1710–1711, 1996 View at: Google Scholar
32. L. Olsen and S. Winter, “Multifractal analysis of divergence points of deformed measure theoretical Birkhoff averages. II: non-linearity, divergence points and Banach space valued spectra,” Bulletin des Sciences Mathématiques, vol. 131, no. 6, pp. 518–558, 2007. View at: Publisher Site | Google Scholar

#### More related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.