Mathematical Problems in Engineering
Volume 2016 (2016), Article ID 1487824, 7 pages
http://dx.doi.org/10.1155/2016/1487824

Research Article

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

Michał Niezabitowski

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]) as

It 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 that

It 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 if

This 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

then

that 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 have

To obtain the last inequality, we used inequality (51). Since

and 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

Denoting

we 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.

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.

J. L. Daleckii and M. G. Krein, Stability of Solutions of Differential Equations in Banach Space, American Mathematical Society, Providence, RI, USA, 1974.

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.

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.

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.

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.

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.

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.

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.

A. Czornik and M. Niezabitowski, “Lyapunov exponents for systems with unbounded coefficients,” Dynamical Systems, vol. 28, no. 2, pp. 140–153, 2013.

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.

O. Perron, “Uber stabilitat und asymptotisches verhalten der integrale von differentialgleichungssystemen,” Mathematische Zeitschrift, vol. 29, no. 1, pp. 129–160, 1929.

O. Perron, “Die Stabilitätsfrage bei Differentialgleichungen,” Mathematische Zeitschrift, vol. 32, no. 1, pp. 703–728, 1930.

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.

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.

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.

A. Czornik, “On the Perron exponents of discrete linear systems,” Linear Algebra and its Applications, vol. 432, no. 1, pp. 394–401, 2010.

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.

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.

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.

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.

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.

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.

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.

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.

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

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

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.