Linear Algebra and its Applications
|
|
- Adam Wood
- 5 years ago
- Views:
Transcription
1 Linear Algebra and its Applications 433 (200) Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: On the exponential exponents of discrete linear systems Adam Czornik, Piotr Mokry, Aleksander Nawrat Department of Automatic Control, Silesian Technical University, ul. Akademicka 6, 44-0 Gliwice, Poland A R T I C L E I N F O A B S T R A C T Article history: Received 0 September 2009 Accepted 6 April 200 Available online 0 May 200 Submitted by H. Schneider Keywords: Time varying discrete linear systems Lyapunov exponents Perturbation theory Characteristic exponents In this paper we introduce the concepts of exponential exponents of discrete linear time varying systems. It is shown that these exponents describe the possible changes in the Lyapunov exponents under perturbation decreasing at infinity at exponential rate. Finally we present formulas for the exponential exponents in terms of the transition matrix of the system. 200 Elsevier Inc. All rights reserved.. Introduction Consider the linear discrete time-varying system x(n + ) = A(n)x(n), n 0 () where (A(n)) n N is a bounded sequence of invertible s-by-s real matrices. By we denote the Euclidean norm in R s and the induced operator norm. The transition matrix is defined as A(m, k) = A(m ) A(k) for m > k and A(m, m) = I, where I is the identity matrix. For an initial condition x 0, the solution of () is denoted by x(n, x 0 ) so x(n, x 0 ) = A(n, 0)x 0. Let a = (a(n)) n N be a sequence of real numbers. The number (or the symbol ± ) defined as λ(a) = lim sup n ln a(n) Corresponding author. address: Adam.Czornik@polsl.pl /$ - see front matter 200 Elsevier Inc. All rights reserved. doi:0.06/j.laa
2 868 A. Czornik et al. / Linear Algebra and its Applications 433 (200) is called the characteristic exponent of the sequence (a(n)) n N.Forx 0 R s with x 0 /= 0, the Lyapunov exponent λ(x 0 ) of () is defined as the characteristic exponent of ( x(n, x 0 ) ) n N, that is λ(x 0 ) = lim sup n ln x(n, x 0 ). In the case of a sequence of matrices = ( (n)) n N we define the characteristic exponent λ( ) as the characteristic exponent of ( (n) ) n N. It is well known 3] that the set of all Lyapunov exponents of system () contains at most s elements, say λ (A) <λ 2 (A) < <λ r (A) < and the set {λ, λ 2,..., λ r } is called the spectrum of (). The greatest and the smallest exponent of () we will denote λ g (A) and λ s (A), respectively. Together with () we consider the following perturbed system z(n + ) = (A(n) + (n)) z(n), (2) where ( (n)) n N is a sequence of s-by-s real matrices from a certain class M. Under the influence of the perturbation ( (n)) n N, the characteristic exponents of () vary, in general, discontinuously. It is possible that a finite shift of the characteristic exponents of the original system () corresponds to an arbitrarily small sup (n). In particular, it is possible for an exponentially stable system to be perturbed by an exponentially decreasing perturbation and the resulting system is not stable. An example which illustrates this situation is given in 6]. The quantity Λ (M) = sup { λ g (A + ) : M } is referred to as the maximal upper movability boundary of the higher exponent of (2) with perturbation in the class M. In the paper we will consider exponentially small perturbations. The determination of the movability boundary of the higher exponent under various perturbations is one of important problems of the modern theory of characteristic exponents. This problem has been solved for continuous-time systems for many classes M. For example, upper bound for the higher exponent of (2) under small perturbations, the so-called central exponent (A), was constructed in 5, p. 4]. The attainability of this estimate was proved in ] with the use of the classical rotation method. This problem was solved in 7,8] for linear systems with perturbations decreasing at infinity at various rates and in 2] for linear systems with perturbations determined by integral conditions. Later in 9,0] perturbations infinitesimal in mean with a weight function have been investigated. Similar problems for discrete-time systems have been investigated in 6]. In this paper we introduce the notion and properties of exponential exponents. We will study the properties of the spectrum and the relations with the exponential exponents, and the influence of these objects on how the exponents vary when passing from ()to(2). According to the authors knowledge, exponential exponents of difference equations has not yet been investigated. 2. Exponential exponents In this paper, M is defined to be the class of perturbations M = {( (n)) n N : λ( ) <0}. Definition. Bounded sequences (r(n)) n N and (R(n)) n N are said to be lower and upper sequences for (), respectively, if for every ε>0 there are constants d r,ε and D R,ε such that for all natural m, k, m k we have m m d r,ε exp r(i) εk A(m, k) D R,ε exp R(i) + εk. (3) The numbers (A) = inf R i=k lim sup R(i), n (A) = sup r i=k lim inf r(i) n
3 A. Czornik et al. / Linear Algebra and its Applications 433 (200) are called the upper exponential exponent and the lower exponential exponent of system (). Here the infimum is taken over the set of all upper sequences whereas the supremum is taken over the set of all lower sequences. We will denote by L(A) and U(A) the set of all lower and upper sequences of (), respectively. Because system () is uniquely defined by the sequence (A(n)) n N we will also talk about lower and upper sequences and lower exponential exponent and upper exponential exponent of sequence (A(n)) n N. The following relations between the exponential and Lyapunov exponents immediately follow from the definition: (A) λ (A) λ r (A) (A). The dual or adjoint equation to () is defined to be y(n + ) = B(n)y(n), n 0, where B(n) = (A T (n)) and A T denotes the transpose of A. The transition matrix of the dual system is given by B(m, k) = B(m ) B(k) for m > k and B(m, m) = I. It appears that the lower exponential exponent does not require special consideration since the problem can be reduced to the investigation of upper exponential exponent for the adjoint system. This is the content of the next theorem. Theorem. The lower exponential exponent of () is equal to the upper exponential exponent of the adjoint system (4) taken with the opposite sign. Proof. By the definition of adjoint system we have A(m, k) = ( B(m, k) ) T. Therefore the estimate m d R,ε exp r(i) εk A(m, k) i=k implies B(m, k) m exp r(i) + εk d R,ε i=k and the latter implies our statement. Now we show that under the condition that λ( ) <0, the exponential exponents of () and (2) coincides. Theorem 2. If λ( ) <0 then L(A) = L(A + ), U(A) = U(A + ) and in particular (A) = (A + ) and (A) = (A + ). In the proof of this theorem we will use the following version of the discrete Gronwall s inequality ]. Theorem 3. Suppose that for two sequences u(n) and f (n),n= m, m +,...of nonnegative real numbers the following inequality u(n) p + q u(i)f (i) (4)
4 870 A. Czornik et al. / Linear Algebra and its Applications 433 (200) holds for certain p, q R and all n = m, m +,..., then u(n) p ( + qf (i)) for all n = m, m +,... (5) Proof. Denote by C(n, m) the transition matrix of (2) and consider an upper sequence R(n) for (). We show that it is also an upper sequence for (2). For the transition matrix C(n, m) we have C(n, m) = A(n, m) + A(n, i + ) (i)c(i, m) and therefore C(n, m) A(n, m) + A(n, i + ) (i) C(i, m). Fix an ε>0 such that ε< λ( ) and estimate A(n, i) according to (3). Then C(n, m) D R,ε exp R(i) + mε + D R,ε exp R(j) + (i + ) ε (i) C(i, m) j=i+ (6) or C(n, m) exp R(i) mε i D R,ε + D R,ε C(i, m) exp R(j) mε (i) exp ( R(i) + (i + ) ε). j=m By the Gronwall s inequality (5) with i u(i) = C(i, m) exp (R(j) + mε), j=m p = D R,ε, q = D R,ε, f(i) = (i) exp ( R(i) + (i + ) ε) we get ( C(n, m) exp R(i) mε D R,ε + DR,ε (i) exp ( R(i) + (i + ) ε) ) and C(n, m) D R,ε exp (mε) ( exp (R(i)) + DR,ε (i) exp ((i + ) ε) ). (7)
5 A. Czornik et al. / Linear Algebra and its Applications 433 (200) By the boundedness of R(n) there exist ε > 0 and C > 0 such that ( ) exp (R(i)) + DR,ε/2 ε C exp R(i) + δ, for all natural numbers n > m and all δ>0. By (6) there exists a positive integer N such that for all n > N (i) exp ((i + ) ε) <ε. By the last two inequalities we obtain from (7) that for all n > m > N the following holds C(n, m) D R,ε C exp R(i) + mε + δ. Increasing if necessary the constant C, we conclude that R(n) is the upper function for (2) because δ>0is arbitrary. Changing the roles of systems () and (2) in the above arguments, we come to the conclusion that an upper function for (2) is also an upper function for (). Thus, the sets of upper functions of these systems coincide; hence the upper exponential exponents also coincide. For the lower functions and lower exponential exponents, the arguments are carried out by passing to the adjoin system. For a sequence of perturbation matrices = ( (n)) n N, we define λ g (A + ) and λ s (A + ) as the greatest and smallest Lyapunov exponent of (2), respectively. Moreover for δ>0define Λ δ = sup λ g (A + ) and λ δ = inf λ s (A + ). λ( )<δ λ( )<δ Because the functions Λ δ and λ δ (considered as functions of δ>0) are nondecreasing and nonincreasing, respectively, the following limits lim Λ δ 0 δ and lim λ δ 0 δ exist and are equal to sup λ g (A + ) and inf λ s (A + ), λ( )<0 λ( )<0 respectively. Theorem 4. The following inequalities hold sup λ g (A + ) (A) and λ( )<0 inf λ s (A + ) (A). λ( )<0 Proof. By Theorem it is enough to show the case of the upper exponential exponent. In that end, we show that for every ε>0, there exists δ>0such that λ g (A + ) (A) + ε for all perturbation matrices = ( (n)) n N satisfying λ( ) <δ. Fix an upper sequence R(n) such that lim sup R(i) (A) + ε n 2, (8) and denote C(n, m) the transition matrix of (2). Repeating the arguments from the proof of Theorem 2 we obtain the following bounds for z(n) = C(n, 0)z(0): ( z(n) D R,ε/4 z(0) exp (R(i)) + DR,ε/4 (i) exp ((i + ) ε/4) ).
6 872 A. Czornik et al. / Linear Algebra and its Applications 433 (200) By (8) wehaveforλ( ) <δ< λ(z 0 ) = lim sup n The proof is completed. ε the inequality 4D R,ε/2 ln z(n, z(0)) (A) + ε + ε + ε = (A) + ε Formula for the exponential exponents In this paragraph we present a formula for the exponential exponents in terms of the transition matrix. We start with a result from 4] which states that we may reduce an arbitrary system to upper triangular form by an unitary transformation. Theorem 5. For each sequence (A(n)) n N there exists a sequence (U(n)) n N of orthogonal matrices such that C(n) = U T (n + )A(n)U(n) is upper triangular. Clearly the orthogonality of U(n) implies that U(A) = U(C), L(A) = L(C). In particular, the exponential exponents of (A(n)) n N and (C(n)) n N are equal. Lemma 6. If matrices (A(n)) n N are upper triangular with diagonal elements a ii (n), then U(A) = U(A d ) and L(A) = L(A d ), where A d (n) = diaga ii (n)] i=,...,n. Proof. We will proof only the part about upper sequences. The proof of the second part is similar. Denote by a ij (n) and z ij (n, k), the elements of A(n) and A(n +,k), respectively. By a straightforward induction we have that for n > m n k j p=0 l=i+ a il(n p)z lj (n p,k) n q=n p+ a ii(q) for i < j, z ij (n, k) = n q=k a ii(q) for i = j, (9) 0 for i > j. Let R U(A). By the definition of upper function, for every ε>0there is a constant D R,ε such that for all n k we have n A(n +,k) D R,ε exp R(l) + kε. (0) Because all norms in finite dimensional space are equivalent there exist positive constants D and D 2 such that D X max xij D2 X, () i,j=,...,s for every matrix X =x ij ] i,j=,...,s. Therefore from (0)wehave n z ii (n, k) D 2 D R,ε exp R(l) + kε ; max i=,...,s hence, by ()wehaver U(A d ). Let P U(A d ). By the definition of upper function, for every ε>0there is a constant D R,ε such that for all n k we have zjj (n, k) n DR,ε exp P(i) + kε (2) i=k
7 A. Czornik et al. / Linear Algebra and its Applications 433 (200) for all j =,...,s. We will show by induction that for every ε>0there is constant D R,ε such that for all n k we have zij (n, k) n D P(l) + kε (3) R,ε exp for all j = i,...,s.fori = j it follows from (2). Suppose that (3) holds for i = j, j,..., α +. By (9) and boundednes of P we obtain zαj (n, k) n k j a αl (n p) zlj (n p,k) n a αα (q) p=0 l=α+ q=n p+ n k j n p n D R,ε a αl (n p) exp P(l) + kε exp P(i) + ε (n p + ) p=0 l=α+ i=n p+ D R,ε n n k exp j (P(l) + 2εk) a αl (n p). p=0 l=α+ From the last inequality, (3) follows by the boundedness of (A(n)) n N. Applying the equivalence of all norms in finite dimensional space we obtain from (3) that P U(A). For each real number >, define α() = lim sup lim inf n ln A ( i+ ], i ]), where x] is the integer part of x. Lemma 7. If (a(n)) n N is a bounded sequence, > and A() = lim sup M M lim sup M M In particular, lim A() exists. + = lim + A(). n ln A( i+ ], i ]) and α() = M ] M, then Proof. Suppose the sequence (a(n)) n N is bounded by D > 0. For each natural number n such that M < n M+, for certain natural number M, wehave M ] n M n M + M M n M n + M ] M i=n M n M D + n M] M D 2D ( ). M ] M
8 874 A. Czornik et al. / Linear Algebra and its Applications 433 (200) Theorem 8. The limits lim α() and lim α() exist and the following equalities hold + + lim α() = (A) and lim α() = (A), + + where (A) and (A) are the upper exponential exponent and lower exponential exponent of (). Proof. By Theorem it is enough to show only the first part; moreover, according to Theorem 5, we may assume that A(n) are upper triangular with elements a ij (n). For a fixed real number >, consider sequence (R (n)) n N. defined in the following way:. If n is such that j n < j, for certain natural number j, then R (n) is equal to ln a ii (n) such that j ] k= j ] a ii(k) is maximal with respect to i =,...,s. Then we have ] ] j j exp R (k) = max a ii (k) k= j ] i=,...,s k= j ] and consequently R U(A d ), where A d (n) = diaga ii (n)] i=,...,n. According to Lemma 6 R U(A) and therefore lim sup R (i) (A). n By Lemma 7 lim lim sup + N N N ] If we define B(m, k) to be the transition matrix of (A d (n)) n N, then B ( j+ ], j ]) = max i=,...,s and consequently lim sup N N N ln R (i) (A). (4) j+ ] k= j ] By the inequality () it is also clear that lim sup N N N ln a ii (k) B ( j+ ], j ]) = lim sup N A ( j+ ], j ]) = lim sup N N N N ] N R (j). Combining (4) and (5) we obtain that the limit lim α() exists and that + lim α() (A). + Fix ε>0 and a function R ε U(A) such that ( ] ]) ln B j+, j. (5) lim sup R ε (i) < (A) + ε. (7) n By Lemma 6, R ε U(A d ) and therefore ( ] ]) B j+, j D R,ε exp j+ ] i= j ] R ε (i) + ε j] (6)
9 A. Czornik et al. / Linear Algebra and its Applications 433 (200) for certain positive D R,ε. Hence, N N ( ] ln B j+, j ]) ln D R,ε N + N ] j+ R ε (i) + ε j] i= j ] and by (5) and (7)wehave lim α() (A). + Inequalities (6) and (8) lead to the conclusion of the theorem. (8) 4. Conclusions In this paper we have considered the influence of exponentially decreasing perturbations on the Lyapunov spectrum of a linear discrete-time system with time-varying coefficients. We have shown that the changes can be described by the exponential exponents. Finally we have presented formulas for the exponential exponents in terms of the transition matrix of the original system. Acknowledgments The research presented here by the first author was done as a part of research project No. N and research of the second and third author was done as a part of research and development project O R and have been supported by Ministry of Science and Higher Education funds in the years References ] R.P. Agarwal, Difference Equations and Inequalities. Theory, Methods, and Applications, Marcel Dekker, New York, ] E.A. Barabanov, O.G. Vishnevskaya, Dokl. Akad. Nauk Belarusi 4 (5) (997) ] L. Barreira, Y.B. Pesin, Lyapunov Exponents and Smooth Ergodic Theory, University Lectures Series, vol. 23, AMS Bookstore, ] L. Barreira, C. Valls, Stability theory and Lyapunov regularity, J. Differential Equations 232 (2) (2007) ] B.F. Bylov, R.E. Vinograd, D.M. Grobman, V.V. Nemytskii, Theory of Lyapunov Exponents and its Applications to Stability Theory, Nauka, Moscow, 966 (Russian). 6] A. Czornik, P. Mokry, A. Nawrat, On the sigma exponent of discrete linear systems, IEEE Trans. Automat. Control, accepted for publication. 7] N.A. Izobov, The higher exponent of a system with perturbations of order higher than one, Vestn. Bel. Gos. Un-ta. Ser. I, 969, No. 3, pp. 6 9 (Russian). 8] N.A. Izobov, The higher exponent of a linear system with exponential perturbations, Differential Equations 5 (7) (969) ] E.K. Makarov, I.V. Marchenko, N.V. Semerikova, On an upper bound for the higher exponent of a linear differential system with integrable perturbations on the half-line, Differ. Uravn. 4 (2) (2005) ] E.K. Makarov, I.V. Marchenko, On an algorithm for constructing an attainable upper boundary for the higher exponent of perturbed systems, Differential Equations 4 (2) (2005) ] V.M. Millionshchikov, A proof of the attainability of the central exponents of linear systems, Siberian Math. J. 0 () (969) ] F. Wirth, On the calculation of time-varying stability radii, Internat. J. Robust Nonlinear Control 8 (998)
Linear Algebra and its Applications
Linear Algebra its Applications 432 21) 394 41 Contents lists available at ScienceDirect Linear Algebra its Applications journal homepage: wwwelseviercom/locate/laa On the Perron exponents of discrete
More informationAN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES
AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES JOEL A. TROPP Abstract. We present an elementary proof that the spectral radius of a matrix A may be obtained using the formula ρ(a) lim
More informationIntrinsic products and factorizations of matrices
Available online at www.sciencedirect.com Linear Algebra and its Applications 428 (2008) 5 3 www.elsevier.com/locate/laa Intrinsic products and factorizations of matrices Miroslav Fiedler Academy of Sciences
More informationPersistence and global stability in discrete models of Lotka Volterra type
J. Math. Anal. Appl. 330 2007 24 33 www.elsevier.com/locate/jmaa Persistence global stability in discrete models of Lotka Volterra type Yoshiaki Muroya 1 Department of Mathematical Sciences, Waseda University,
More informationCONTROLLABILITY OF NONLINEAR DISCRETE SYSTEMS
Int. J. Appl. Math. Comput. Sci., 2002, Vol.2, No.2, 73 80 CONTROLLABILITY OF NONLINEAR DISCRETE SYSTEMS JERZY KLAMKA Institute of Automatic Control, Silesian University of Technology ul. Akademicka 6,
More informationAn example of a convex body without symmetric projections.
An example of a convex body without symmetric projections. E. D. Gluskin A. E. Litvak N. Tomczak-Jaegermann Abstract Many crucial results of the asymptotic theory of symmetric convex bodies were extended
More informationMATH 117 LECTURE NOTES
MATH 117 LECTURE NOTES XIN ZHOU Abstract. This is the set of lecture notes for Math 117 during Fall quarter of 2017 at UC Santa Barbara. The lectures follow closely the textbook [1]. Contents 1. The set
More informationLecture 7. Econ August 18
Lecture 7 Econ 2001 2015 August 18 Lecture 7 Outline First, the theorem of the maximum, an amazing result about continuity in optimization problems. Then, we start linear algebra, mostly looking at familiar
More informationCONTINUITY OF SOLUTIONS OF RICCATI EQUATIONS FOR THE DISCRETE TIME JLQP
Int. J. Appl. Math. Comput. Sci., 2002, Vol.12, No.4, 539 543 CONTINUITY O SOLUTIONS O RICCATI EQUATIONS OR THE DISCRETE TIME JLQP ADAM CZORNIK, ANDRZEJ ŚWIERNIAK Department of Automatic Control Silesian
More informationThe Free Central Limit Theorem: A Combinatorial Approach
The Free Central Limit Theorem: A Combinatorial Approach by Dennis Stauffer A project submitted to the Department of Mathematical Sciences in conformity with the requirements for Math 4301 (Honour s Seminar)
More informationInternational Competition in Mathematics for Universtiy Students in Plovdiv, Bulgaria 1994
International Competition in Mathematics for Universtiy Students in Plovdiv, Bulgaria 1994 1 PROBLEMS AND SOLUTIONS First day July 29, 1994 Problem 1. 13 points a Let A be a n n, n 2, symmetric, invertible
More informationASYMPTOTIC BEHAVIOUR OF SECOND-ORDER DIFFERENCE EQUATIONS
ANZIAM J. 462004, 57 70 ASYMPTOTIC BEHAVIOUR OF SECOND-ORDER DIFFERENCE EQUATIONS STEVO STEVIĆ Received 9 December, 200; revised 9 September, 2003 Abstract In this paper we prove several growth theorems
More informationLinear Algebra and its Applications
Linear Algebra and its Applications 436 2012 41 52 Contents lists available at SciVerse ScienceDirect Linear Algebra and its Applications journal homepage: www.elsevier.com/locate/laa Compactness of matrix
More informationBounds for Levinger s function of nonnegative almost skew-symmetric matrices
Linear Algebra and its Applications 416 006 759 77 www.elsevier.com/locate/laa Bounds for Levinger s function of nonnegative almost skew-symmetric matrices Panayiotis J. Psarrakos a, Michael J. Tsatsomeros
More information1. Supremum and Infimum Remark: In this sections, all the subsets of R are assumed to be nonempty.
1. Supremum and Infimum Remark: In this sections, all the subsets of R are assumed to be nonempty. Let E be a subset of R. We say that E is bounded above if there exists a real number U such that x U for
More informationLinear Algebra and its Applications
Linear Algebra and its Applications 430 (2009) 579 586 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: www.elsevier.com/locate/laa Low rank perturbation
More informationBanach Algebras of Matrix Transformations Between Spaces of Strongly Bounded and Summable Sequences
Advances in Dynamical Systems and Applications ISSN 0973-532, Volume 6, Number, pp. 9 09 20 http://campus.mst.edu/adsa Banach Algebras of Matrix Transformations Between Spaces of Strongly Bounded and Summable
More informationPacific Journal of Mathematics
Pacific Journal of Mathematics OPTIMAL OSCILLATION CRITERIA FOR FIRST ORDER DIFFERENCE EQUATIONS WITH DELAY ARGUMENT GEORGE E. CHATZARAKIS, ROMAN KOPLATADZE AND IOANNIS P. STAVROULAKIS Volume 235 No. 1
More informationBoolean Inner-Product Spaces and Boolean Matrices
Boolean Inner-Product Spaces and Boolean Matrices Stan Gudder Department of Mathematics, University of Denver, Denver CO 80208 Frédéric Latrémolière Department of Mathematics, University of Denver, Denver
More informationA property of orthogonal projectors
Linear Algebra and its Applications 354 (2002) 35 39 www.elsevier.com/locate/laa A property of orthogonal projectors Jerzy K. Baksalary a,, Oskar Maria Baksalary b,tomaszszulc c a Department of Mathematics,
More informationASYMPTOTIC BEHAVIOR OF SOLUTIONS OF DISCRETE VOLTERRA EQUATIONS. Janusz Migda and Małgorzata Migda
Opuscula Math. 36, no. 2 (2016), 265 278 http://dx.doi.org/10.7494/opmath.2016.36.2.265 Opuscula Mathematica ASYMPTOTIC BEHAVIOR OF SOLUTIONS OF DISCRETE VOLTERRA EQUATIONS Janusz Migda and Małgorzata
More informationOrthogonal similarity of a real matrix and its transpose
Available online at www.sciencedirect.com Linear Algebra and its Applications 428 (2008) 382 392 www.elsevier.com/locate/laa Orthogonal similarity of a real matrix and its transpose J. Vermeer Delft University
More informationA ZERO ENTROPY T SUCH THAT THE [T,ID] ENDOMORPHISM IS NONSTANDARD
A ZERO ENTROPY T SUCH THAT THE [T,ID] ENDOMORPHISM IS NONSTANDARD CHRISTOPHER HOFFMAN Abstract. We present an example of an ergodic transformation T, a variant of a zero entropy non loosely Bernoulli map
More informationOn some linear combinations of hypergeneralized projectors
Linear Algebra and its Applications 413 (2006) 264 273 www.elsevier.com/locate/laa On some linear combinations of hypergeneralized projectors Jerzy K. Baksalary a, Oskar Maria Baksalary b,, Jürgen Groß
More informationOptimization Theory. A Concise Introduction. Jiongmin Yong
October 11, 017 16:5 ws-book9x6 Book Title Optimization Theory 017-08-Lecture Notes page 1 1 Optimization Theory A Concise Introduction Jiongmin Yong Optimization Theory 017-08-Lecture Notes page Optimization
More informationConvergence Rate of Nonlinear Switched Systems
Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the
More informationLecture 8 : Eigenvalues and Eigenvectors
CPS290: Algorithmic Foundations of Data Science February 24, 2017 Lecture 8 : Eigenvalues and Eigenvectors Lecturer: Kamesh Munagala Scribe: Kamesh Munagala Hermitian Matrices It is simpler to begin with
More informationExistence of Almost Periodic Solutions of Discrete Ricker Delay Models
International Journal of Difference Equations ISSN 0973-6069, Volume 9, Number 2, pp. 187 205 (2014) http://campus.mst.edu/ijde Existence of Almost Periodic Solutions of Discrete Ricker Delay Models Yoshihiro
More informationIntertibility and spectrum of the multiplication operator on the space of square-summable sequences
Intertibility and spectrum of the multiplication operator on the space of square-summable sequences Objectives Establish an invertibility criterion and calculate the spectrum of the multiplication operator
More informationACCURATE SOLUTION ESTIMATE AND ASYMPTOTIC BEHAVIOR OF NONLINEAR DISCRETE SYSTEM
Sutra: International Journal of Mathematical Science Education c Technomathematics Research Foundation Vol. 1, No. 1,9-15, 2008 ACCURATE SOLUTION ESTIMATE AND ASYMPTOTIC BEHAVIOR OF NONLINEAR DISCRETE
More informationLecture 5. Ch. 5, Norms for vectors and matrices. Norms for vectors and matrices Why?
KTH ROYAL INSTITUTE OF TECHNOLOGY Norms for vectors and matrices Why? Lecture 5 Ch. 5, Norms for vectors and matrices Emil Björnson/Magnus Jansson/Mats Bengtsson April 27, 2016 Problem: Measure size of
More informationON ADAPTIVE CONTROL FOR THE CONTINUOUS TIME-VARYING JLQG PROBLEM
Int. J. Appl. Math. Comput. Sci., 25, Vol. 15, No. 1, 53 62 B ON ADAPTIVE CONTROL FOR THE CONTINUOUS TIME-VARYING JLQG PROBLEM ADAM CZORNIK, ANDRZEJ ŚWIERNIAK Department of Automatic Control Silesian University
More informationA LITTLE REAL ANALYSIS AND TOPOLOGY
A LITTLE REAL ANALYSIS AND TOPOLOGY 1. NOTATION Before we begin some notational definitions are useful. (1) Z = {, 3, 2, 1, 0, 1, 2, 3, }is the set of integers. (2) Q = { a b : aεz, bεz {0}} is the set
More informationC.7. Numerical series. Pag. 147 Proof of the converging criteria for series. Theorem 5.29 (Comparison test) Let a k and b k be positive-term series
C.7 Numerical series Pag. 147 Proof of the converging criteria for series Theorem 5.29 (Comparison test) Let and be positive-term series such that 0, for any k 0. i) If the series converges, then also
More informationCovering an ellipsoid with equal balls
Journal of Combinatorial Theory, Series A 113 (2006) 1667 1676 www.elsevier.com/locate/jcta Covering an ellipsoid with equal balls Ilya Dumer College of Engineering, University of California, Riverside,
More informationFALSENESS OF THE FINITENESS PROPERTY OF THE SPECTRAL SUBRADIUS
Int J Appl Math Comput Sci, 2007, Vol 17, No 2, 173 178 DOI: 102478/v10006-007-0016-1 FALSENESS OF THE FINITENESS PROPERTY OF THE SPECTRAL SUBRADIUS ADAM CZORNIK, PIOTR JURGAŚ Department of Automatic Control,
More information5 Compact linear operators
5 Compact linear operators One of the most important results of Linear Algebra is that for every selfadjoint linear map A on a finite-dimensional space, there exists a basis consisting of eigenvectors.
More informationMAT 570 REAL ANALYSIS LECTURE NOTES. Contents. 1. Sets Functions Countability Axiom of choice Equivalence relations 9
MAT 570 REAL ANALYSIS LECTURE NOTES PROFESSOR: JOHN QUIGG SEMESTER: FALL 204 Contents. Sets 2 2. Functions 5 3. Countability 7 4. Axiom of choice 8 5. Equivalence relations 9 6. Real numbers 9 7. Extended
More informationSpectrum for compact operators on Banach spaces
Submitted to Journal of the Mathematical Society of Japan Spectrum for compact operators on Banach spaces By Luis Barreira, Davor Dragičević Claudia Valls (Received Nov. 22, 203) (Revised Jan. 23, 204)
More informationAbsolute value equations
Linear Algebra and its Applications 419 (2006) 359 367 www.elsevier.com/locate/laa Absolute value equations O.L. Mangasarian, R.R. Meyer Computer Sciences Department, University of Wisconsin, 1210 West
More informationDiscrete Applied Mathematics
Discrete Applied Mathematics 194 (015) 37 59 Contents lists available at ScienceDirect Discrete Applied Mathematics journal homepage: wwwelseviercom/locate/dam Loopy, Hankel, and combinatorially skew-hankel
More informationTHEOREM OF OSELEDETS. We recall some basic facts and terminology relative to linear cocycles and the multiplicative ergodic theorem of Oseledets [1].
THEOREM OF OSELEDETS We recall some basic facts and terminology relative to linear cocycles and the multiplicative ergodic theorem of Oseledets []. 0.. Cocycles over maps. Let µ be a probability measure
More informationMath 108b: Notes on the Spectral Theorem
Math 108b: Notes on the Spectral Theorem From section 6.3, we know that every linear operator T on a finite dimensional inner product space V has an adjoint. (T is defined as the unique linear operator
More informationYimin Wei a,b,,1, Xiezhang Li c,2, Fanbin Bu d, Fuzhen Zhang e. Abstract
Linear Algebra and its Applications 49 (006) 765 77 wwwelseviercom/locate/laa Relative perturbation bounds for the eigenvalues of diagonalizable and singular matrices Application of perturbation theory
More informationLinear Algebra and its Applications
Linear Algebra and its Applications 431 (2009) 274 296 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: www.elsevier.com/locate/laa Boolean inner-product
More informationNew Lyapunov Krasovskii functionals for stability of linear retarded and neutral type systems
Systems & Control Letters 43 (21 39 319 www.elsevier.com/locate/sysconle New Lyapunov Krasovskii functionals for stability of linear retarded and neutral type systems E. Fridman Department of Electrical
More informationConvex Geometry. Otto-von-Guericke Universität Magdeburg. Applications of the Brascamp-Lieb and Barthe inequalities. Exercise 12.
Applications of the Brascamp-Lieb and Barthe inequalities Exercise 12.1 Show that if m Ker M i {0} then both BL-I) and B-I) hold trivially. Exercise 12.2 Let λ 0, 1) and let f, g, h : R 0 R 0 be measurable
More informationShih-sen Chang, Yeol Je Cho, and Haiyun Zhou
J. Korean Math. Soc. 38 (2001), No. 6, pp. 1245 1260 DEMI-CLOSED PRINCIPLE AND WEAK CONVERGENCE PROBLEMS FOR ASYMPTOTICALLY NONEXPANSIVE MAPPINGS Shih-sen Chang, Yeol Je Cho, and Haiyun Zhou Abstract.
More informationMath Introduction to Numerical Analysis - Class Notes. Fernando Guevara Vasquez. Version Date: January 17, 2012.
Math 5620 - Introduction to Numerical Analysis - Class Notes Fernando Guevara Vasquez Version 1990. Date: January 17, 2012. 3 Contents 1. Disclaimer 4 Chapter 1. Iterative methods for solving linear systems
More informationFunctions of several variables of finite variation and their differentiability
ANNALES POLONICI MATHEMATICI LX.1 (1994) Functions of several variables of finite variation and their differentiability by Dariusz Idczak ( Lódź) Abstract. Some differentiability properties of functions
More informationEquivalence constants for certain matrix norms II
Linear Algebra and its Applications 420 (2007) 388 399 www.elsevier.com/locate/laa Equivalence constants for certain matrix norms II Bao Qi Feng a,, Andrew Tonge b a Department of Mathematical Sciences,
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More informationA property concerning the Hadamard powers of inverse M-matrices
Linear Algebra and its Applications 381 (2004 53 60 www.elsevier.com/locate/laa A property concerning the Hadamard powers of inverse M-matrices Shencan Chen Department of Mathematics, Fuzhou University,
More informationSTRONG NONNEGATIVE LINEARIZATION OF ORTHOGONAL POLYNOMIALS
STRONG NONNEGATIVE LINEARIZATION OF ORTHOGONAL POLYNOMIALS Ryszard Szwarc * Institute of Mathematics Wroclaw University pl. Grunwaldzki 2/4 50-384 Wroclaw POLAND and Institute of Mathematics Polish Academy
More informationExercises to Applied Functional Analysis
Exercises to Applied Functional Analysis Exercises to Lecture 1 Here are some exercises about metric spaces. Some of the solutions can be found in my own additional lecture notes on Blackboard, as the
More informationUniform convergence of N-dimensional Walsh Fourier series
STUDIA MATHEMATICA 68 2005 Uniform convergence of N-dimensional Walsh Fourier series by U. Goginava Tbilisi Abstract. We establish conditions on the partial moduli of continuity which guarantee uniform
More informationHilbert spaces. 1. Cauchy-Schwarz-Bunyakowsky inequality
(October 29, 2016) Hilbert spaces Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/fun/notes 2016-17/03 hsp.pdf] Hilbert spaces are
More informationLinear and Multilinear Algebra. Linear maps preserving rank of tensor products of matrices
Linear maps preserving rank of tensor products of matrices Journal: Manuscript ID: GLMA-0-0 Manuscript Type: Original Article Date Submitted by the Author: -Aug-0 Complete List of Authors: Zheng, Baodong;
More informationWEIGHTED LIMITS FOR POINCARÉ DIFFERENCE EQUATIONS. Rotchana Chieocan
WEIGHTED LIMITS FOR POINCARÉ DIFFERENCE EQUATIONS Rotchana Chieocan Department of Mathematics, Faculty of Science, Khon Kaen University, 40002 Khon Kaen, Thail rotchana@kku.ac.th Mihály Pituk* Department
More informationInterpolating the arithmetic geometric mean inequality and its operator version
Linear Algebra and its Applications 413 (006) 355 363 www.elsevier.com/locate/laa Interpolating the arithmetic geometric mean inequality and its operator version Rajendra Bhatia Indian Statistical Institute,
More informationLaplacian spectral radius of trees with given maximum degree
Available online at www.sciencedirect.com Linear Algebra and its Applications 429 (2008) 1962 1969 www.elsevier.com/locate/laa Laplacian spectral radius of trees with given maximum degree Aimei Yu a,,1,
More informationOSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix
Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra
More informationWhere is matrix multiplication locally open?
Linear Algebra and its Applications 517 (2017) 167 176 Contents lists available at ScienceDirect Linear Algebra and its Applications www.elsevier.com/locate/laa Where is matrix multiplication locally open?
More informationLinear Algebra and its Applications
Linear Algebra and its Applications 432 2010 661 669 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: wwwelseviercom/locate/laa On the characteristic and
More informationThe Drazin inverses of products and differences of orthogonal projections
J Math Anal Appl 335 7 64 71 wwwelseviercom/locate/jmaa The Drazin inverses of products and differences of orthogonal projections Chun Yuan Deng School of Mathematics Science, South China Normal University,
More informationMetric Spaces and Topology
Chapter 2 Metric Spaces and Topology From an engineering perspective, the most important way to construct a topology on a set is to define the topology in terms of a metric on the set. This approach underlies
More informationPOSITIVE PARTIAL REALIZATION PROBLEM FOR LINEAR DISCRETE TIME SYSTEMS
Int J Appl Math Comput Sci 7 Vol 7 No 65 7 DOI: 478/v6-7-5- POSITIVE PARTIAL REALIZATION PROBLEM FOR LINEAR DISCRETE TIME SYSTEMS TADEUSZ KACZOREK Faculty of Electrical Engineering Białystok Technical
More informationProjection-valued measures and spectral integrals
Projection-valued measures and spectral integrals Jordan Bell jordan.bell@gmail.com Department of Mathematics, University of Toronto April 16, 2014 Abstract The purpose of these notes is to precisely define
More informationInequalities for the numerical radius, the norm and the maximum of the real part of bounded linear operators in Hilbert spaces
Available online at www.sciencedirect.com Linear Algebra its Applications 48 008 980 994 www.elsevier.com/locate/laa Inequalities for the numerical radius, the norm the maximum of the real part of bounded
More informationOn the digital representation of smooth numbers
On the digital representation of smooth numbers Yann BUGEAUD and Haime KANEKO Abstract. Let b 2 be an integer. Among other results, we establish, in a quantitative form, that any sufficiently large integer
More informationAn introduction to some aspects of functional analysis
An introduction to some aspects of functional analysis Stephen Semmes Rice University Abstract These informal notes deal with some very basic objects in functional analysis, including norms and seminorms
More informationThis article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author s institution, sharing
More informationIowa State University. Instructor: Alex Roitershtein Summer Homework #1. Solutions
Math 501 Iowa State University Introduction to Real Analysis Department of Mathematics Instructor: Alex Roitershtein Summer 015 EXERCISES FROM CHAPTER 1 Homework #1 Solutions The following version of the
More informationFinal Project # 5 The Cartan matrix of a Root System
8.099 Final Project # 5 The Cartan matrix of a Root System Thomas R. Covert July 30, 00 A root system in a Euclidean space V with a symmetric positive definite inner product, is a finite set of elements
More informationLinear Algebra and its Applications
Linear Algebra and its Applications 433 (2010) 476 482 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: www.elsevier.com/locate/laa Nonsingularity of the
More informationRANK AND PERIMETER PRESERVER OF RANK-1 MATRICES OVER MAX ALGEBRA
Discussiones Mathematicae General Algebra and Applications 23 (2003 ) 125 137 RANK AND PERIMETER PRESERVER OF RANK-1 MATRICES OVER MAX ALGEBRA Seok-Zun Song and Kyung-Tae Kang Department of Mathematics,
More informationh(x) lim H(x) = lim Since h is nondecreasing then h(x) 0 for all x, and if h is discontinuous at a point x then H(x) > 0. Denote
Real Variables, Fall 4 Problem set 4 Solution suggestions Exercise. Let f be of bounded variation on [a, b]. Show that for each c (a, b), lim x c f(x) and lim x c f(x) exist. Prove that a monotone function
More informationCOMMON COMPLEMENTS OF TWO SUBSPACES OF A HILBERT SPACE
COMMON COMPLEMENTS OF TWO SUBSPACES OF A HILBERT SPACE MICHAEL LAUZON AND SERGEI TREIL Abstract. In this paper we find a necessary and sufficient condition for two closed subspaces, X and Y, of a Hilbert
More informationCHEBYSHEV INEQUALITIES AND SELF-DUAL CONES
CHEBYSHEV INEQUALITIES AND SELF-DUAL CONES ZDZISŁAW OTACHEL Dept. of Applied Mathematics and Computer Science University of Life Sciences in Lublin Akademicka 13, 20-950 Lublin, Poland EMail: zdzislaw.otachel@up.lublin.pl
More informationSome Results Concerning Uniqueness of Triangle Sequences
Some Results Concerning Uniqueness of Triangle Sequences T. Cheslack-Postava A. Diesl M. Lepinski A. Schuyler August 12 1999 Abstract In this paper we will begin by reviewing the triangle iteration. We
More informationSpectral Theorem for Self-adjoint Linear Operators
Notes for the undergraduate lecture by David Adams. (These are the notes I would write if I was teaching a course on this topic. I have included more material than I will cover in the 45 minute lecture;
More informationLecture notes: Applied linear algebra Part 1. Version 2
Lecture notes: Applied linear algebra Part 1. Version 2 Michael Karow Berlin University of Technology karow@math.tu-berlin.de October 2, 2008 1 Notation, basic notions and facts 1.1 Subspaces, range and
More informationGenerating Functions of Switched Linear Systems: Analysis, Computation, and Stability Applications
IEEE-TAC FP--5 Generating Functions of Switched Linear Systems: Analysis, Computation, and Stability Applications Jianghai Hu, Member, IEEE, Jinglai Shen, Member, IEEE, and Wei Zhang, Member, IEEE Abstract
More informationDISCRETE GRONWALL LEMMA AND APPLICATIONS
DISCRETE GRONWALL LEMMA AND APPLICATIONS JOHN M. HOLTE Variations of Gronwall s Lemma Gronwall s lemma, which solves a certain kind of inequality for a function, is useful in the theory of differential
More informationCONTINUITY OF SUBADDITIVE PRESSURE FOR SELF-AFFINE SETS
RESEARCH Real Analysis Exchange Vol. 34(2), 2008/2009, pp. 1 16 Kenneth Falconer, Mathematical Institute, University of St Andrews, North Haugh, St Andrews KY16 9SS, Scotland. email: kjf@st-and.ac.uk Arron
More informationLecture 7: Semidefinite programming
CS 766/QIC 820 Theory of Quantum Information (Fall 2011) Lecture 7: Semidefinite programming This lecture is on semidefinite programming, which is a powerful technique from both an analytic and computational
More informationPart 2 Continuous functions and their properties
Part 2 Continuous functions and their properties 2.1 Definition Definition A function f is continuous at a R if, and only if, that is lim f (x) = f (a), x a ε > 0, δ > 0, x, x a < δ f (x) f (a) < ε. Notice
More informationReceived 8 June 2003 Submitted by Z.-J. Ruan
J. Math. Anal. Appl. 289 2004) 266 278 www.elsevier.com/locate/jmaa The equivalence between the convergences of Ishikawa and Mann iterations for an asymptotically nonexpansive in the intermediate sense
More informationReceived: 7 September 2017; Accepted: 22 October 2017; Published: 25 October 2017
entropy Communication A New Definition of t-entropy for Transfer Operators Victor I. Bakhtin 1,2, * and Andrei V. Lebedev 2,3 1 Faculty of Mathematics, Informatics and Landscape Architecture, John Paul
More informationStability analysis of differential algebraic systems
Stability analysis of differential algebraic systems Volker Mehrmann TU Berlin, Institut für Mathematik with Vu Hoang Linh and Erik Van Vleck DFG Research Center MATHEON Mathematics for key technologies
More informationMATH5011 Real Analysis I. Exercise 1 Suggested Solution
MATH5011 Real Analysis I Exercise 1 Suggested Solution Notations in the notes are used. (1) Show that every open set in R can be written as a countable union of mutually disjoint open intervals. Hint:
More informationSupremum and Infimum
Supremum and Infimum UBC M0 Lecture Notes by Philip D. Loewen The Real Number System. Work hard to construct from the axioms a set R with special elements O and I, and a subset P R, and mappings A: R R
More informationModeling and Stability Analysis of a Communication Network System
Modeling and Stability Analysis of a Communication Network System Zvi Retchkiman Königsberg Instituto Politecnico Nacional e-mail: mzvi@cic.ipn.mx Abstract In this work, the modeling and stability problem
More informationAn estimation of the spectral radius of a product of block matrices
Linear Algebra and its Applications 379 (2004) 267 275 wwwelseviercom/locate/laa An estimation of the spectral radius of a product of block matrices Mei-Qin Chen a Xiezhang Li b a Department of Mathematics
More informationEXISTENCE OF CONJUGACIES AND STABLE MANIFOLDS VIA SUSPENSIONS
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 172, pp. 1 11. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF CONJUGACIES AND STABLE MANIFOLDS
More informationImpulsive stabilization of two kinds of second-order linear delay differential equations
J. Math. Anal. Appl. 91 (004) 70 81 www.elsevier.com/locate/jmaa Impulsive stabilization of two kinds of second-order linear delay differential equations Xiang Li a, and Peixuan Weng b,1 a Department of
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationUndergraduate Notes in Mathematics. Arkansas Tech University Department of Mathematics
Undergraduate Notes in Mathematics Arkansas Tech University Department of Mathematics An Introductory Single Variable Real Analysis: A Learning Approach through Problem Solving Marcel B. Finan c All Rights
More information(K + L)(c x) = K(c x) + L(c x) (def of K + L) = K( x) + K( y) + L( x) + L( y) (K, L are linear) = (K L)( x) + (K L)( y).
Exercise 71 We have L( x) = x 1 L( v 1 ) + x 2 L( v 2 ) + + x n L( v n ) n = x i (a 1i w 1 + a 2i w 2 + + a mi w m ) i=1 ( n ) ( n ) ( n ) = x i a 1i w 1 + x i a 2i w 2 + + x i a mi w m i=1 Therefore y
More information