Mathematical Problems in Engineering

Volume 2016 (2016), Article ID 1487824, 7 pages

http://dx.doi.org/10.1155/2016/1487824

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

Faculty of Automatic Control, Electronics and Computer Science, Silesian University of Technology, Akademicka 16 Street, 44-100 Gliwice, Poland

Received 29 July 2016; Accepted 12 October 2016

Academic Editor: Sotiris K. Ntouyas

Copyright © 2016 Michał Niezabitowski. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### 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 [1]) 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 [2].

The Lyapunov [3–15] and the Perron [16–23] 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 [24–29]). 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.

#### References

- A. Czornik,
*Perturbation Theory for Lyapunov Exponents of Discrete Linear Systems*, Wydawnictwa AGH, 2012. - 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. - A. Lyapunov,
*A general task about the stability of motion [Ph.D. thesis]*, 1892. - 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 · View at Google Scholar - 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 - 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 - 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 · View at Google Scholar · View at MathSciNet · View at Scopus - 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 · View at Google Scholar · View at MathSciNet · View at Scopus - 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 · View at MathSciNet - 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 · View at Google Scholar · View at MathSciNet - 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 · View at Google Scholar · View at MathSciNet - 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 · View at Google Scholar · View at MathSciNet · View at Scopus - 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 · View at Google Scholar · View at MathSciNet - 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 · View at Google Scholar · View at MathSciNet · View at Scopus - 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 · View at Google Scholar · View at MathSciNet · View at Scopus - O. Perron, “Uber stabilitat und asymptotisches verhalten der integrale von differentialgleichungssystemen,”
*Mathematische Zeitschrift*, vol. 29, no. 1, pp. 129–160, 1929. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - O. Perron, “Die Stabilitätsfrage bei Differentialgleichungen,”
*Mathematische Zeitschrift*, vol. 32, no. 1, pp. 703–728, 1930. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - 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 · View at MathSciNet - 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 · View at Google Scholar · View at MathSciNet · View at Scopus - 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 · View at Google Scholar - 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 · View at Google Scholar · View at MathSciNet · View at Scopus - 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 Scopus - 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 · View at Google Scholar - 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 · View at Google Scholar · View at MathSciNet · View at Scopus - 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 Scopus - 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 · View at Google Scholar · View at MathSciNet · View at Scopus - 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 · View at Google Scholar · View at MathSciNet · View at Scopus - 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 · View at Google Scholar - 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 · View at Google Scholar - 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 - 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 - 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 · View at Google Scholar