A New Invariance Property of Lyapunov Characteristic Directions S. Bharadwaj and K.D. Mease Mechanical and Aerospace Engineering University of Califor
|
|
- Ruth Bryant
- 5 years ago
- Views:
Transcription
1 A New Invariance Property of Lyapunov Characteristic Directions S. Bharadwaj and K.D. Mease Mechanical and Aerospace Engineering University of California, Irvine, California, Abstract Lyapunov exponents and direction æelds are used to characterize the time-scales and geometry of general linear time-varying LTV systems of diæerential equations. Lyapunov exponents are already known to correctly characterize the time-scales present in a general LTV system; they reduce to real parts of eigenvalues when computed for linear time-invariantlti systems and real parts of Floquet exponents when computed for periodic LTV systems. Here, we bring to light new invariance properties of Lyapunov direction æelds to show that they are analogous to the Schur vectors of an LTI system and reduce to the Schur vectors when computed for LTI systems. We also show that the Lyapunov direction æeld corresponding to the smallest Lyapunov exponent when computed for an LTI system with real distinct eigenvalues reduces to the eigenvector corresponding to the smallest eigenvalue and when computed for a periodic LTV system with real distinct Floquet exponents, reduces to the Floquet direction æeld corresponding 1 Introduction In the eæort towards understanding the geometric structure and qualitative behavior of the state-space æow of smooth ænite-dimensional nonlinear dynamical systems of the form _x = fx ; x 2 R n 1 where fæ : R n! TR n is a smooth vector æeld, a commonly adopted ærst step is to study the linearized dynamics about a reference orbit ææ of Eq.1, i.e., _æ = f æt. Linear time-varying LTV systems of interest to us typically arise from such a linearization about a reference orbit which we assume to be smooth and bounded. The resulting LTV system is given æ _x where æx 2 T x R n and where At 2 R næn to have bounded entries. = Atæx 2 is assumed Recall that the transition matrix æt; associated with the LTV system is = Atæ ; æ;=i n 3 where I n is the identity matrix of order n. Unlike LTI systems, no closed form expression exists in general relating the transition matrix æ of an LTV system to At. Also note that æxt =æt; æx for all t; 2 R. The LTV system determines how tangent vectors are mapped between the diæerent tangent spaces T æt R n along the reference trajectory. Knowledge of the directions of tangent vectors that growdecay at extremal rates when mapped by the æow determined by the LTV system and their corresponding rates of expansioncontraction yields information about convergencedivergence of state trajectories neighboring ææ and their corresponding rates. A complete spectral characterization of the LTV dynamics aims at uncovering the diæerent asymptotic rates of growthdecay of solutions of the LTV system along with their corresponding directions. Geometrically speaking, we seek linearly independent time-varying direction æelds along ææ that identify the diæerent characteristic directions on each tangent space T æt R n along the reference trajectory. When the reference trajectory is an isolated equilibrium point of Eq.1, i.e., æt = æ 0 for all t 2 R, the resulting LTV system becomes linear time-invariant LTI and is determined by At = A = constant for all t 2 R. The spectral structure of an LTI system, æ _x = Aæx, is completely characterized by the eigenvalues and eigenvectors of the matrix A. When the reference trajectory is an isolated periodic orbit of Eq.1, i.e, æt + kt =æt for all t 2 R;k 2Zwhere T 2 R + is the minimal period, the spectral structure of the resulting periodic LTV system is determined by Floquet theory ë3ë. Avast body of literature ë1, 8, 10, 11, 12ë exists on eæorts to characterize the spectral structure of general LTV systems which arise when the reference trajectory is not necessarily either an equilibrium point or a periodic orbit, beginning with the pioneering work
2 of Lyapunov ë7ë in his 1892 thesis. Lyapunov introduced the notion characteristic exponents to characterize solutions of LTV dynamics with extremal evolution rates which have since received a lot of attention from theoretical and computational perspectives ë1, 2ë. The same cannot be said of the associated Lyapunov direction æelds. By investigating the invariance properties of Lyapunov direction æelds we show that in addition to identifying directions of extremal growthdecay of solutions, they span invariant distributions see ë6ë for a description of distributions analogous to the invariant subspaces spanned by Schur vectors of LTI dynamics. 1.1 Lyapunov Exponents The theory of Lyapunov exponents is described in wellknown texts ë5, 9ë. For regular LTV systems ë1, 2, 9ë, the Lyapunov exponents are well-deæned limits on tangent subbundles of ææ on which the evolution rates are extremal. They are given by ë5ë 1 i ëæx i ë = lim ln kæx i tk t!1 t 1 i ëæx i ë = lim ln kæt; æx i k t!1 t, where æx i t ; 1 i n constitute a set of normal basis solutions ë5ë of Eq.2. In the sequel, we assume that LTV systems of interest possess the property of regularity. It is important to note that Lyapunov exponents are constants over the entire trajectory ææ and are independent of the starting point æ. Also, the exponents are known to be invariant under Lyapunov transformations ë5ë. Lyapunov exponents computed for an LTI system yield the real parts of the eigenvalues and when computed for a periodic LTV system yield the real parts of the Floquet exponents see ë5ë,thm Lyapunov Directions The Lyapunov directions at any given time along the reference trajectory ææ point in the direction of extremal average rates of growthdecay. So, they can be determined at any initial time by vectors æx i which extremize kæt; æx i k as t!1. The solutions æx i to this extremization problem is given by the eigenvectors of lim t!1 æ T t; æt; 1=2t ë2ë. Greene and Kim ë4ë have shown that when the timescales given by the Lyapunov exponents are distinct, the corresponding Lyapunov directions depend only on the position in state space given by æ. The Lyapunov direction æelds are everywhere orthogonal and determine directions at each point along the reference trajectory along which the average growthdecay rate is extremal. However, Lyapunov direction æelds in general are not solutions of the LTV system. For example, when the Lyapunov directions are propagated forward from some time they do not in general remain orthogonal. In other words, not all of the Lyapunov direction æelds are invariant under the linear æow of the LTV dynamics. In the following section we investigate the invariance properties of Lyapunov direction æelds for regular LTV systems with distinct Lyapunov exponents. 2 Invariance Properties of Lyapunov Direction æelds The following lemma shows that when a particular Lyapunov direction is propagated forward by the linear æow, the time-evolved direction can develop components only along Lyapunov direction æelds associated with smaller Lyapunov exponents. Lemma 1 Consider a regular LTV system of the form 2 which has the distinct ordered Lyapunov exponents i ;i=1;::: ;n where 1 é:::é n and corresponding time-varying directions l i t;i=1;::: ;n for all t 2 R. Propagate these directions l i forward from an initial time to t; t é using the transition matrix so that m i t =æt; l i. Then 0 if iéj; hm i t;l j ti= 4 nonzero if i j where hæ; æi denotes the inner product and 1 i; j n. Proof: Recall that the Lyapunov exponents i are computed by 1 i ël i ë= lim ln kæt; l i k 5 t!1 t, where the Lyapunov directions l i t are mutually orthogonal at each t 2 R because they arise from the computation of eigenvectors of a symmetric matrix. Therefore 0 if i 6= j; hl i ;l j i =æ ij = 1 if i = j for all 1 i; j n. However, when propagated forward from time to t 1 é using æt 1 ;, the propagated vectors m i t 1 =æt 1 ;l i are in general no longer mutually orthogonal. Let us assume that hm i t 1 ;l j t 1 i6= 0 for all iéjand show that we arrive at a contradiction. Since vectors l j t 1 are linearly independent in fact orthogonal, we can express the vectors m i t 1 as linear combinations of l j t 1. For each i, let m i t 1 = nx j=1 c ij l j t 1 6
3 The above assumption then corresponds to c ij 6= 0 for all iéj. Eq. 6 implies that for each i, æt; t 1 m i t 1 = nx j=1 c ij æt; t 1 l j t 1 As t becomes large, the right hand side is dominated by the term corresponding to the vector l j ;j = j æ for which æt; t 1 l j æt 1 is the largest and c ij æ 6= 0. Since we have assumed the Lyapunov exponents which determine the rates of change in size of vectors when propagated by the linear æow to be arranged in descending order, the dominant term corresponds to when j æ is the smallest positive integer for which c ij æ 6= 0. Recall that the Lyapunov exponents i are constants along the reference trajectory ææ and independent of the starting point æ on the reference trajectory. Compute the Lyapunov exponents from the initial directions m i t 1 using 1 i = lim ln kæt; t 1 m i t 1 k t!1 t, t 1 7 In the light of the earlier argument, we can see that kæt; t 1 m i t 1 k!c ij ækæt; t 1 l j æk as t!1and that i = j æ. This contradicts the fact that Lyapunov exponents are constants independent of the starting point on the reference trajectory unless j æ = i. The condition j æ = i holds if and only if c ij = 0 for all i é j. Since this discussion is valid for any time t 1 é, we drop the subscript and conclude that hm i t;l j ti= for all times té. 0 if iéj; nonzero if i j A direct consequence of the above lemma is the following theorem. Theorem 1 For any 1 k n, the Lyapunov direction æelds corresponding to the k smallest Lyapunov exponents of the LTV system 2 deæne an invariant distribution ææt = spanfl n,k+1 t;::: ;l n tg. Proof: Consider the Lyapunov directions l n,k+1 ;::: ;l n at some initial time which correspond to the k smallest Lyapunov exponents. When propagated forward to time t, let m i t =æt; l i for all n, k + 1 i n. Lemma 1 implies that at time t, spanfm n,k+1 t;::: ;m n tg = spanfl n,k+1 t;::: ;l n tg. This is because, the directions m n,k+1 t;::: ;m n t cannot have components along any of the directions l 1 t;::: ;l n,k t. Consequently, ifv2ææ then æt; v 2 ææt for all té. Therefore æ is an invariant distribution. The above theorem and and lemma show that for regular LTV systems with distinct Lyapunov exponents the Lyapunov direction æelds deæne a collection of invariant distributions æ k æ = spanfl k ;::: ;l n g such that dimæ k = n, k +1 and æ 1 ::: æ n. When æxt 2 æ i then kæxtk Ke it where K is a positive constant. When these results are applied to an LTI system typically resulting from linearizing nonlinear dynamics about an equilibrium point æ 0, the distributions æ i become subspaces W i of T æ0 R n. So the Lyapunov directions reduce to Schur vectors that span the respective nested invariant subspaces W 1 ::: W n. This establishes that the Lyapunov direction æelds for an LTV system are time-varying analogs of Schur vectors for an LTI system. Since a special case of the above theorem where k =1 has many important consequences, we present it as a separate result. Theorem 2 The direction æeld l n æ corresponding to the smallest Lyapunov exponent n is invariant when propagated using the transition matrix, i.e., æt; l n = n t; l n t. Proof: For the LTV system 2 with Lyapunov exponents 1 é ::: é n, l n t is the direction æeld corresponding to the smallest Lyapunov exponent n. Applying Lemma 1 to the direction l n results in hæt; l n ;l j ti=0; 8j6=n 8 This implies that at time t, the direction l n t points in the same direction as m n t =æt; l n, i.e, æt; l n = n t; l n t 9 where n t; is the scaling factor that determines that rate of growth of the size of l n. Hence wehave shown that the Lyapunov direction æeld l n æ corresponding to the smallest Lyapunov exponent is invariant when propagated by the transition matrix. The following important corollaries and examples further justify the validity of the direction information provided by this Lyapunov direction æeld as an accurate generalization by showing that it is indeed consistent with the known results in the special cases of LTI systems and periodic LTV systems. Corollary 1 The Lyapunov direction corresponding to the smallest Lyapunov exponent, computed for an LTI
4 system æ _x = Aæx which is assumed to have real and distinct eigenvalues ordered from largest to smallest, points in the same direction as the eigenvector direction corresponding to the smallest eigenvalue. Proof: As stated earlier, it is already known that the Lyapunov exponents of an LTI system are the same as real parts of the eigenvalues of A. Theorem 2 states that the Lyapunov direction l n corresponding to the smallest Lyapunov exponent satisæes the condition of invariance when propagated using æt; = expa æ t,. Therefore æt; l n = n t, l n where n = exp n æ t,, i.e., l n is an eigenvector of æ with eigenvalue n which is the same as the eigenvector of A corresponding to the smallest eigenvalue. The following example shows that the Lyapunov directions computed for an LTI system reduce to the Schur vectors and in addition illustrates the above corollary. Example 1 Consider the LTI system æ _x = Aæx where,6 5 A = 10 4,5 The eigenvalues of A are 1 =,1 ; 2 =,10 and the corresponding eigenvectors are e 1 = ë1; 1ë and e 2 = ë5;,4ë respectively. The Schur vectors obtained by orthogonalizing the eigenvectors are s 2 = e 2 and s 1 = ë4; 5ë. The transition matrix of this LTI system is given by æt; = expa æ t, = 1=9ëæ ij ë where æ 11 = 4e,t, +5e,10t, æ 12 =,5e,10t, +5e,t, æ 21 =,4e,10t, +4e,t, æ 22 = 5e,t, +4e,10t, The eigenvalues of æ T t; æt; are q 1;2 t; =e,2t, w æ w 2, 6561e,18t, where w =41e,18t,, e,9t, +41. The Lyapunov exponents are given by 1 1;2 = lim t!1 2t, ln 1;2 11 As t!1, 1 t;! a exp,2t, and 2 t;! b exp,20t, where a and b are some constants whose true value is unimportant for our purpose. Therefore we ænd the Lyapunov exponents 1 =,1 and 2 =,10 which are the same as the eigenvalues 1 and 2 of A. The corresponding Lyapunov directions are given by the eigenvectors of æ T t; æt; as t!1. They are l 1;2 = lim t!1,9e,9t, +9æ p w 2,6561e,18t, 41e,18t,, e,9t,, 40 Therefore, l 2 = ë50;,40ë and l 1 = ë,32;,40ë which point in the same directions as the Schur vectors s 2 and s 1 respectively. Also note that l 2 points in the same direction as e 2. We shall now show through the following corollary that the result of Theorem 2 is also consistent with the results of Floquet theory when applied to periodic LTV systems. Corollary 2 The Lyapunov direction æeld corresponding to the smallest Lyapunov exponent, computed for a periodic LTV system æ _x = Atæx with At+T =At which is assumed to have real and distinct Floquet exponents ordered from largest to smallest points in the same direction as the Floquet direction æeld corresponding Proof: As stated earlier, it is already known that the Lyapunov exponents of a periodic LTV system are the same as real parts of the Floquet exponents. Theorem 2 states that the Lyapunov direction æeld l n æ corresponding to the smallest Lyapunov exponent satisæes the invariance condition. Therefore æt; l n = n t; l n t 12 where n = exp n æ t,. We know that the Lyapunov direction æelds are dependent only on the position in state space along the reference trajectory ææ. In this case, ææ is periodic, i.e., æ + T = æ. Consequently, l n + T =l n, i.e., l n æ is periodic. Together with the invariance condition, we get æ + T;l n = exp n T l n 13 This implies that l n is an eigenvector of æ + T; corresponding to Floquet exponent n = n. This shows that l n æ is a periodic direction æeld pointing in the same direction as the Floquet direction æeld corresponding We shall illustrate this corollary with the following example. Example 2 Consider the periodic LTV system æ _x =
5 Atæx where At =ëa ij të so that A 11 t =, 1 2 cos 2t sin 2t, 11 2 A 12 t = 1 2 sin 2t cos 2t A 21 t = 1 2 sin 2t cos 2t, 3 2 A 22 t =, 9 2 sin 2t cos 2t, 11 2 with period T =. This system was constructed from the LTI system in example 1represented here asæ_x= Bæx with a periodic transformation æx = Stæx where,6 5 cos t, sin t B = ; St = 14 4,5 sin t cos t so that At =S,1 tbst, S,1 t _St. The Floquet exponents turn out to be 1 =,1 and 2 =,10. The corresponding periodic Floquet direction æelds are given by e 1 t = cos t + sin t, sin t + cos t ; e 2 t = 5 cos t, 4 sin t,5 sin t, 4 cos t The transition matrix for this periodic LTV system is given by æt; =S T tæt; S 15 where æt; is the transition matrix in example 1. The Lyapunov exponents and associated directions are computed from the eigenvalues and eigenvectors of æ T t; æt; as t!1. Note that æ T t; æt; =S T æ T t; æt; S 16 because StS T t =I n. Consequently, the eigenvalues of æ T æ are the same as the eigenvalues of æ T æ and its eigenvectors are related to that of æ T æ by the transformation S T. Therefore, the Lyapunov exponents of the periodic LTV system turn out to be 1 =,1 and 2 =,10 which are the same as the Floquet exponents. The direction æeld corresponding to the smallest Lyapunov exponent 2 computed from the eigenvector of æ T t; æt; as t!1is given by 5 l 2 =S T =,4 5 cos, 4 sin,5 sin, 4 cos 17 which is the same as the Floquet direction æeld corresponding 3 Conclusions Lyapunov direction æelds span invariant distributions on which the average asymptotic evolution rates of solutions of the LTV dynamics are extremal. They reduce to Schur vectors when computed for an LTI system thereby establishing that they are appropriate timevarying analogs of Schur vectors. In particular, the Lyapunov direction æeld corresponding to the smallest Lyapunov exponent is invariant under the linear æow and this direction æeld is consistent with the corresponding eigenvector when computed for LTI dynamics and with the corresponding Floquet direction æeld when computed for periodic LTV dynamics. Illuminating examples were presented to illustrate these results. References ë1ë F. Colonius and W. Kliemann, ëthe Lyapunov Spectrum of Families of Time-Varying Matrices," Transactions of The American Mathematical Society, Vol. 348, No. 11, pp. 4389í4408, ë2ë L. Dieci, R. D. Russell and E. S. Van Vleck, ëon the Computation of Lyapunov Exponents for Continuous Dynamical Systems," SIAM Journal of Numerical Analysis, Vol. 34, No. 1, pp. 402í423, ë3ë M. Farkas, Periodic Motions, Springer-Verlag, New York, ë4ë J. M. Greene, and J.-S. Kim, ëthe Calculation of Lyapunov Spectra," Physica D, Vol. 24, pp. 213í225, ë5ë W. Hahn, Stability of Motion, Springer-Verlag, New York, ë6ë A. Isidori, ënonlinear Control Systems", Second edition, Springer, ë7ë A. M. Lyapunov, ëthe General Problem of Stability of Motion," International Journal of Control, Vol. 55, No. 3, pp , ë8ë R. J. Sacker and G. R. Sell, ëa Spectral Theory for Linear Diæerential Systems," Journal of Diæerential Equations, Vol. 27, pp , ë9ë G. Sansone and R. Conti, Nonlinear Diæerential Equations. Pergamon Press, New York, ë10ë W. E. Wiesel, ëoptimal Pole Placement in Time- Dependent Linear Systems," Journal of Guidance, Control and Dynamics, Vol. 18, No. 5, pp. 995í999, ë11ë M.-Y. Wu, ëa New Concept of Eigenvalues and Eigenvectors and its Applications," IEEE Transactions on Automatic Control, Vol. 25, No. 4, pp. 824í826, ë12ë J. J. Zhu, ëa Uniæed Spectral Theory for Linear Time-Varying Systems - Progress and Challenges," Proceedings of the 34th IEEE Conference on Decision and Control, New Orleans, IEEE, New York, Vol. 3, pp. 2540í2546, 1995.
Observations on the Stability Properties of Cooperative Systems
1 Observations on the Stability Properties of Cooperative Systems Oliver Mason and Mark Verwoerd Abstract We extend two fundamental properties of positive linear time-invariant (LTI) systems to homogeneous
More informationPassivity-based Stabilization of Non-Compact Sets
Passivity-based Stabilization of Non-Compact Sets Mohamed I. El-Hawwary and Manfredi Maggiore Abstract We investigate the stabilization of closed sets for passive nonlinear systems which are contained
More informationIn the previous chapters we have presented synthesis methods for optimal H 2 and
Chapter 8 Robust performance problems In the previous chapters we have presented synthesis methods for optimal H 2 and H1 control problems, and studied the robust stabilization problem with respect to
More informationFinding eigenvalues for matrices acting on subspaces
Finding eigenvalues for matrices acting on subspaces Jakeniah Christiansen Department of Mathematics and Statistics Calvin College Grand Rapids, MI 49546 Faculty advisor: Prof Todd Kapitula Department
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 informationj=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent.
Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. Let u = [u
More informationMath 5BI: Problem Set 6 Gradient dynamical systems
Math 5BI: Problem Set 6 Gradient dynamical systems April 25, 2007 Recall that if f(x) = f(x 1, x 2,..., x n ) is a smooth function of n variables, the gradient of f is the vector field f(x) = ( f)(x 1,
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 informationOutput Input Stability and Minimum-Phase Nonlinear Systems
422 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 47, NO. 3, MARCH 2002 Output Input Stability and Minimum-Phase Nonlinear Systems Daniel Liberzon, Member, IEEE, A. Stephen Morse, Fellow, IEEE, and Eduardo
More informationMa/CS 6b Class 23: Eigenvalues in Regular Graphs
Ma/CS 6b Class 3: Eigenvalues in Regular Graphs By Adam Sheffer Recall: The Spectrum of a Graph Consider a graph G = V, E and let A be the adjacency matrix of G. The eigenvalues of G are the eigenvalues
More informationJordan Normal Form. Chapter Minimal Polynomials
Chapter 8 Jordan Normal Form 81 Minimal Polynomials Recall p A (x) =det(xi A) is called the characteristic polynomial of the matrix A Theorem 811 Let A M n Then there exists a unique monic polynomial q
More informationLECTURE Review. In this lecture we shall study the errors and stability properties for numerical solutions of initial value.
LECTURE 24 Error Analysis for Multi-step Methods 1. Review In this lecture we shall study the errors and stability properties for numerical solutions of initial value problems of the form è24.1è dx = fèt;
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
More informationThis property turns out to be a general property of eigenvectors of a symmetric A that correspond to distinct eigenvalues as we shall see later.
34 To obtain an eigenvector x 2 0 2 for l 2 = 0, define: B 2 A - l 2 I 2 = È 1, 1, 1 Î 1-0 È 1, 0, 0 Î 1 = È 1, 1, 1 Î 1. To transform B 2 into an upper triangular matrix, subtract the first row of B 2
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationMinimum-Phase Property of Nonlinear Systems in Terms of a Dissipation Inequality
Minimum-Phase Property of Nonlinear Systems in Terms of a Dissipation Inequality Christian Ebenbauer Institute for Systems Theory in Engineering, University of Stuttgart, 70550 Stuttgart, Germany ce@ist.uni-stuttgart.de
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 informationJune 2011 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations
June 20 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations The topics covered in this exam can be found in An introduction to differential equations
More informationLOCALLY POSITIVE NONLINEAR SYSTEMS
Int. J. Appl. Math. Comput. Sci. 3 Vol. 3 No. 4 55 59 LOCALLY POSITIVE NONLINEAR SYSTEMS TADEUSZ KACZOREK Institute of Control Industrial Electronics Warsaw University of Technology ul. Koszykowa 75 66
More informationOn the simultaneous diagonal stability of a pair of positive linear systems
On the simultaneous diagonal stability of a pair of positive linear systems Oliver Mason Hamilton Institute NUI Maynooth Ireland Robert Shorten Hamilton Institute NUI Maynooth Ireland Abstract In this
More informationMATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors.
MATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors. Orthogonal sets Let V be a vector space with an inner product. Definition. Nonzero vectors v 1,v
More information4 Second-Order Systems
4 Second-Order Systems Second-order autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization
More informationFundamentals of Linear Algebra. Marcel B. Finan Arkansas Tech University c All Rights Reserved
Fundamentals of Linear Algebra Marcel B. Finan Arkansas Tech University c All Rights Reserved 2 PREFACE Linear algebra has evolved as a branch of mathematics with wide range of applications to the natural
More informationA note on linear differential equations with periodic coefficients.
A note on linear differential equations with periodic coefficients. Maite Grau (1) and Daniel Peralta-Salas (2) (1) Departament de Matemàtica. Universitat de Lleida. Avda. Jaume II, 69. 251 Lleida, Spain.
More informationFamily Feud Review. Linear Algebra. October 22, 2013
Review Linear Algebra October 22, 2013 Question 1 Let A and B be matrices. If AB is a 4 7 matrix, then determine the dimensions of A and B if A has 19 columns. Answer 1 Answer A is a 4 19 matrix, while
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationLINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM
LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM Unless otherwise stated, all vector spaces in this worksheet are finite dimensional and the scalar field F is R or C. Definition 1. A linear operator
More informationOn the convergence of normal forms for analytic control systems
Chapter 1 On the convergence of normal forms for analytic control systems Wei Kang Department of Mathematics, Naval Postgraduate School Monterey, CA 93943 wkang@nps.navy.mil Arthur J. Krener Department
More informationproblem of detection naturally arises in technical diagnostics, where one is interested in detecting cracks, corrosion, or any other defect in a sampl
In: Structural and Multidisciplinary Optimization, N. Olhoæ and G. I. N. Rozvany eds, Pergamon, 1995, 543í548. BOUNDS FOR DETECTABILITY OF MATERIAL'S DAMAGE BY NOISY ELECTRICAL MEASUREMENTS Elena CHERKAEVA
More informationThe University of Texas at Austin Department of Electrical and Computer Engineering. EE381V: Large Scale Learning Spring 2013.
The University of Texas at Austin Department of Electrical and Computer Engineering EE381V: Large Scale Learning Spring 2013 Assignment Two Caramanis/Sanghavi Due: Tuesday, Feb. 19, 2013. Computational
More informationChapter 2: Linear Independence and Bases
MATH20300: Linear Algebra 2 (2016 Chapter 2: Linear Independence and Bases 1 Linear Combinations and Spans Example 11 Consider the vector v (1, 1 R 2 What is the smallest subspace of (the real vector space
More informationLecture 1. Introduction. The importance, ubiquity, and complexity of embedded systems are growing
Lecture 1 Introduction Karl Henrik Johansson The importance, ubiquity, and complexity of embedded systems are growing tremendously thanks to the revolution in digital technology. This has created a need
More informationarxiv: v1 [math.cv] 14 Nov 2017
EXTENDING A FUNCTION JUST BY MULTIPLYING AND DIVIDING FUNCTION VALUES: SMOOTHNESS AND PRIME IDENTITIES arxiv:1711.07887v1 [math.cv] 14 Nov 2017 PATRICK ARTHUR MILLER ABSTRACT. We describe a purely-multiplicative
More information4 Vector Spaces. 4.1 Basic Definition and Examples. Lecture 10
Lecture 10 4 Vector Spaces 4.1 Basic Definition and Examples Throughout mathematics we come across many types objects which can be added and multiplied by scalars to arrive at similar types of objects.
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
More informationThere is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth)
82 Introduction Liapunov Functions Besides the Liapunov spectral theorem, there is another basic method of proving stability that is a generalization of the energy method we have seen in the introductory
More information0.1 Tangent Spaces and Lagrange Multipliers
01 TANGENT SPACES AND LAGRANGE MULTIPLIERS 1 01 Tangent Spaces and Lagrange Multipliers If a differentiable function G = (G 1,, G k ) : E n+k E k then the surface S defined by S = { x G( x) = v} is called
More informationMath 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.
Math 350 Fall 2011 Notes about inner product spaces In this notes we state and prove some important properties of inner product spaces. First, recall the dot product on R n : if x, y R n, say x = (x 1,...,
More informationW 1 æw 2 G + 0 e? u K y Figure 5.1: Control of uncertain system. For MIMO systems, the normbounded uncertainty description is generalized by assuming
Chapter 5 Robust stability and the H1 norm An important application of the H1 control problem arises when studying robustness against model uncertainties. It turns out that the condition that a control
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 information8 Periodic Linear Di erential Equations - Floquet Theory
8 Periodic Linear Di erential Equations - Floquet Theory The general theory of time varying linear di erential equations _x(t) = A(t)x(t) is still amazingly incomplete. Only for certain classes of functions
More informationPROBLEMS In each of Problems 1 through 12:
6.5 Impulse Functions 33 which is the formal solution of the given problem. It is also possible to write y in the form 0, t < 5, y = 5 e (t 5/ sin 5 (t 5, t 5. ( The graph of Eq. ( is shown in Figure 6.5.3.
More informationLinear Algebra, 4th day, Thursday 7/1/04 REU Info:
Linear Algebra, 4th day, Thursday 7/1/04 REU 004. Info http//people.cs.uchicago.edu/laci/reu04. Instructor Laszlo Babai Scribe Nick Gurski 1 Linear maps We shall study the notion of maps between vector
More informationChapter 6: Orthogonality
Chapter 6: Orthogonality (Last Updated: November 7, 7) These notes are derived primarily from Linear Algebra and its applications by David Lay (4ed). A few theorems have been moved around.. Inner products
More informationNonlinear Discrete-Time Observer Design with Linearizable Error Dynamics
622 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 48, NO. 4, APRIL 2003 Nonlinear Discrete-Time Observer Design with Linearizable Error Dynamics MingQing Xiao, Nikolaos Kazantzis, Costas Kravaris, Arthur
More informationElementary linear algebra
Chapter 1 Elementary linear algebra 1.1 Vector spaces Vector spaces owe their importance to the fact that so many models arising in the solutions of specific problems turn out to be vector spaces. The
More information322 HENDRA GUNAWAN AND MASHADI èivè kx; y + zk çkx; yk + kx; zk: The pair èx; kæ; ækè is then called a 2-normed space. A standard example of a 2-norme
SOOCHOW JOURNAL OF MATHEMATICS Volume 27, No. 3, pp. 321-329, July 2001 ON FINITE DIMENSIONAL 2-NORMED SPACES BY HENDRA GUNAWAN AND MASHADI Abstract. In this note, we shall study ænite dimensional 2-normed
More informationAsymptotic results for empirical measures of weighted sums of independent random variables
Asymptotic results for empirical measures of weighted sums of independent random variables B. Bercu and W. Bryc University Bordeaux 1, France Workshop on Limit Theorems, University Paris 1 Paris, January
More informationMemoryless output feedback nullification and canonical forms, for time varying systems
Memoryless output feedback nullification and canonical forms, for time varying systems Gera Weiss May 19, 2005 Abstract We study the possibility of nullifying time-varying systems with memoryless output
More information1. General Vector Spaces
1.1. Vector space axioms. 1. General Vector Spaces Definition 1.1. Let V be a nonempty set of objects on which the operations of addition and scalar multiplication are defined. By addition we mean a rule
More informationMATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.
MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation
More informationSIMPLE AND POSITIVE ROOTS
SIMPLE AND POSITIVE ROOTS JUHA VALKAMA MASSACHUSETTS INSTITUTE OF TECHNOLOGY Let V be a Euclidean space, i.e. a real finite dimensional linear space with a symmetric positive definite inner product,. We
More informationAlternative Characterization of Ergodicity for Doubly Stochastic Chains
Alternative Characterization of Ergodicity for Doubly Stochastic Chains Behrouz Touri and Angelia Nedić Abstract In this paper we discuss the ergodicity of stochastic and doubly stochastic chains. We define
More informationLimit Cycles in High-Resolution Quantized Feedback Systems
Limit Cycles in High-Resolution Quantized Feedback Systems Li Hong Idris Lim School of Engineering University of Glasgow Glasgow, United Kingdom LiHonIdris.Lim@glasgow.ac.uk Ai Poh Loh Department of Electrical
More informationContinuity of Bçezier patches. Jana Pçlnikovça, Jaroslav Plaçcek, Juraj ç Sofranko. Faculty of Mathematics and Physics. Comenius University
Continuity of Bezier patches Jana Plnikova, Jaroslav Placek, Juraj Sofranko Faculty of Mathematics and Physics Comenius University Bratislava Abstract The paper is concerned about the question of smooth
More informationBare-bones outline of eigenvalue theory and the Jordan canonical form
Bare-bones outline of eigenvalue theory and the Jordan canonical form April 3, 2007 N.B.: You should also consult the text/class notes for worked examples. Let F be a field, let V be a finite-dimensional
More informationNONCOMMUTATIVE POLYNOMIAL EQUATIONS. Edward S. Letzter. Introduction
NONCOMMUTATIVE POLYNOMIAL EQUATIONS Edward S Letzter Introduction My aim in these notes is twofold: First, to briefly review some linear algebra Second, to provide you with some new tools and techniques
More informationPractical Stabilization of Integrator Switched Systems
Practical Stabilization of Integrator Switched Systems Xuping Xu and Panos J. Antsaklis 2 Department of Electrical and Computer Engineering Penn State Erie, Erie, PA 6563, USA. E-mail: Xuping-Xu@psu.edu
More informationDetailed Proof of The PerronFrobenius Theorem
Detailed Proof of The PerronFrobenius Theorem Arseny M Shur Ural Federal University October 30, 2016 1 Introduction This famous theorem has numerous applications, but to apply it you should understand
More informationRecall the convention that, for us, all vectors are column vectors.
Some linear algebra Recall the convention that, for us, all vectors are column vectors. 1. Symmetric matrices Let A be a real matrix. Recall that a complex number λ is an eigenvalue of A if there exists
More informationTHEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)
4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M
More information18.06 Problem Set 8 - Solutions Due Wednesday, 14 November 2007 at 4 pm in
806 Problem Set 8 - Solutions Due Wednesday, 4 November 2007 at 4 pm in 2-06 08 03 Problem : 205+5+5+5 Consider the matrix A 02 07 a Check that A is a positive Markov matrix, and find its steady state
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More information4 Sums of Independent Random Variables
4 Sums of Independent Random Variables Standing Assumptions: Assume throughout this section that (,F,P) is a fixed probability space and that X 1, X 2, X 3,... are independent real-valued random variables
More informationTopic 2 Quiz 2. choice C implies B and B implies C. correct-choice C implies B, but B does not imply C
Topic 1 Quiz 1 text A reduced row-echelon form of a 3 by 4 matrix can have how many leading one s? choice must have 3 choice may have 1, 2, or 3 correct-choice may have 0, 1, 2, or 3 choice may have 0,
More informationSPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS
SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS G. RAMESH Contents Introduction 1 1. Bounded Operators 1 1.3. Examples 3 2. Compact Operators 5 2.1. Properties 6 3. The Spectral Theorem 9 3.3. Self-adjoint
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
More informationLyapunov-based control of quantum systems
Lyapunov-based control of quantum systems Symeon Grivopoulos Bassam Bamieh Department of Mechanical and Environmental Engineering University of California, Santa Barbara, CA 936-57 symeon,bamieh@engineering.ucsb.edu
More informationSolution via Laplace transform and matrix exponential
EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative
More informationDiagonalizing Matrices
Diagonalizing Matrices Massoud Malek A A Let A = A k be an n n non-singular matrix and let B = A = [B, B,, B k,, B n ] Then A n A B = A A 0 0 A k [B, B,, B k,, B n ] = 0 0 = I n 0 A n Notice that A i B
More informationConsensus of Information Under Dynamically Changing Interaction Topologies
Consensus of Information Under Dynamically Changing Interaction Topologies Wei Ren and Randal W. Beard Abstract This paper considers the problem of information consensus among multiple agents in the presence
More informationMATRICES ARE SIMILAR TO TRIANGULAR MATRICES
MATRICES ARE SIMILAR TO TRIANGULAR MATRICES 1 Complex matrices Recall that the complex numbers are given by a + ib where a and b are real and i is the imaginary unity, ie, i 2 = 1 In what we describe below,
More informationQUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY
QUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY D.D. Olesky 1 Department of Computer Science University of Victoria Victoria, B.C. V8W 3P6 Michael Tsatsomeros Department of Mathematics
More informationSTOCHASTIC PROCESSES Basic notions
J. Virtamo 38.3143 Queueing Theory / Stochastic processes 1 STOCHASTIC PROCESSES Basic notions Often the systems we consider evolve in time and we are interested in their dynamic behaviour, usually involving
More informationWhat is A + B? What is A B? What is AB? What is BA? What is A 2? and B = QUESTION 2. What is the reduced row echelon matrix of A =
STUDENT S COMPANIONS IN BASIC MATH: THE ELEVENTH Matrix Reloaded by Block Buster Presumably you know the first part of matrix story, including its basic operations (addition and multiplication) and row
More informationLecture Notes 1: Vector spaces
Optimization-based data analysis Fall 2017 Lecture Notes 1: Vector spaces In this chapter we review certain basic concepts of linear algebra, highlighting their application to signal processing. 1 Vector
More informationDefinition 1. A set V is a vector space over the scalar field F {R, C} iff. there are two operations defined on V, called vector addition
6 Vector Spaces with Inned Product Basis and Dimension Section Objective(s): Vector Spaces and Subspaces Linear (In)dependence Basis and Dimension Inner Product 6 Vector Spaces and Subspaces Definition
More informationSTABILITY OF 2D SWITCHED LINEAR SYSTEMS WITH CONIC SWITCHING. 1. Introduction
STABILITY OF D SWITCHED LINEAR SYSTEMS WITH CONIC SWITCHING H LENS, M E BROUCKE, AND A ARAPOSTATHIS 1 Introduction The obective of this paper is to explore necessary and sufficient conditions for stability
More informationarxiv: v1 [math.pr] 22 May 2008
THE LEAST SINGULAR VALUE OF A RANDOM SQUARE MATRIX IS O(n 1/2 ) arxiv:0805.3407v1 [math.pr] 22 May 2008 MARK RUDELSON AND ROMAN VERSHYNIN Abstract. Let A be a matrix whose entries are real i.i.d. centered
More information1. What is the determinant of the following matrix? a 1 a 2 4a 3 2a 2 b 1 b 2 4b 3 2b c 1. = 4, then det
What is the determinant of the following matrix? 3 4 3 4 3 4 4 3 A 0 B 8 C 55 D 0 E 60 If det a a a 3 b b b 3 c c c 3 = 4, then det a a 4a 3 a b b 4b 3 b c c c 3 c = A 8 B 6 C 4 D E 3 Let A be an n n matrix
More informationConsensus Seeking in Multi-agent Systems Under Dynamically Changing Interaction Topologies
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, SUBMITTED FOR PUBLICATION AS A TECHNICAL NOTE. 1 Consensus Seeking in Multi-agent Systems Under Dynamically Changing Interaction Topologies Wei Ren, Student Member,
More informationInner Product Spaces
Inner Product Spaces Introduction Recall in the lecture on vector spaces that geometric vectors (i.e. vectors in two and three-dimensional Cartesian space have the properties of addition, subtraction,
More informationAddendum: Lyapunov Exponent Calculation
Addendum: Lyapunov Exponent Calculation Experiment CP-Ly Equations of Motion The system phase u is written in vector notation as θ u = (1) φ The equation of motion can be expressed compactly using the
More informationZero controllability in discrete-time structured systems
1 Zero controllability in discrete-time structured systems Jacob van der Woude arxiv:173.8394v1 [math.oc] 24 Mar 217 Abstract In this paper we consider complex dynamical networks modeled by means of state
More informationMAT Linear Algebra Collection of sample exams
MAT 342 - Linear Algebra Collection of sample exams A-x. (0 pts Give the precise definition of the row echelon form. 2. ( 0 pts After performing row reductions on the augmented matrix for a certain system
More informationRemark By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.
Sec 6 Eigenvalues and Eigenvectors Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called an eigenvalue of A if there is a nontrivial
More informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More information(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);
STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend
More informationIr O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v )
Section 3.2 Theorem 3.6. Let A be an m n matrix of rank r. Then r m, r n, and, by means of a finite number of elementary row and column operations, A can be transformed into the matrix ( ) Ir O D = 1 O
More informationChapter Two Elements of Linear Algebra
Chapter Two Elements of Linear Algebra Previously, in chapter one, we have considered single first order differential equations involving a single unknown function. In the next chapter we will begin to
More information2. Laws of Exponents (1) b 0 1 (2) b x b y b x y (3) bx b y. b x y (4) b n (5) b r s b rs (6) n b b 1/n Example: Solve the equations (a) e 2x
7.1 Derivative of Exponential Function 1. Exponential Functions Let b 0 and b 1. f x b x is called an exponential function. Domain of f x :, Range of f x : 0, Example: lim x e x x2 2. Laws of Exponents
More information9 Facta Universitatis ser.: Elect. and Energ. vol. 11, No.3 è1998è this paper we have considered shaping gain for two interesting quantization procedu
FACTA UNIVERSITATIS èniçsè Series: Electronics and Energetics vol. 11, No.3 è1998è, 91-99 NONLINEAR TRANSFORMATION OF ONEíDIMENSIONAL CONSTELLATION POINTS IN ORDER TO ERROR PROBABILITY DECREASING Zoran
More informationHegselmann-Krause Dynamics: An Upper Bound on Termination Time
Hegselmann-Krause Dynamics: An Upper Bound on Termination Time B. Touri Coordinated Science Laboratory University of Illinois Urbana, IL 680 touri@illinois.edu A. Nedić Industrial and Enterprise Systems
More informationHamiltonian Dynamics In The Theory of Abstraction
Hamiltonian Dynamics In The Theory of Abstraction Subhajit Ganguly. email: gangulysubhajit@indiatimes.com Abstract: This paper deals with fluid flow dynamics which may be Hamiltonian in nature and yet
More informationModule 06 Stability of Dynamical Systems
Module 06 Stability of Dynamical Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha October 10, 2017 Ahmad F. Taha Module 06
More informationFundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad
Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian
More informationRow and Column Representatives in Qualitative Analysis of Arbitrary Switching Positive Systems
ROMANIAN JOURNAL OF INFORMATION SCIENCE AND TECHNOLOGY Volume 19, Numbers 1-2, 2016, 127 136 Row and Column Representatives in Qualitative Analysis of Arbitrary Switching Positive Systems Octavian PASTRAVANU,
More informationOn Linear Copositive Lyapunov Functions and the Stability of Switched Positive Linear Systems
1 On Linear Copositive Lyapunov Functions and the Stability of Switched Positive Linear Systems O. Mason and R. Shorten Abstract We consider the problem of common linear copositive function existence for
More informationConjugate Gradient (CG) Method
Conjugate Gradient (CG) Method by K. Ozawa 1 Introduction In the series of this lecture, I will introduce the conjugate gradient method, which solves efficiently large scale sparse linear simultaneous
More information