Products of Positive Definite Symplectic Matrices
Total Page:16
File Type:pdf, Size:1020Kb
Load more
Recommended publications
-
Modeling, Inference and Clustering for Equivalence Classes of 3-D Orientations Chuanlong Du Iowa State University
Iowa State University Capstones, Theses and Graduate Theses and Dissertations Dissertations 2014 Modeling, inference and clustering for equivalence classes of 3-D orientations Chuanlong Du Iowa State University Follow this and additional works at: https://lib.dr.iastate.edu/etd Part of the Statistics and Probability Commons Recommended Citation Du, Chuanlong, "Modeling, inference and clustering for equivalence classes of 3-D orientations" (2014). Graduate Theses and Dissertations. 13738. https://lib.dr.iastate.edu/etd/13738 This Dissertation is brought to you for free and open access by the Iowa State University Capstones, Theses and Dissertations at Iowa State University Digital Repository. It has been accepted for inclusion in Graduate Theses and Dissertations by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. Modeling, inference and clustering for equivalence classes of 3-D orientations by Chuanlong Du A dissertation submitted to the graduate faculty in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Major: Statistics Program of Study Committee: Stephen Vardeman, Co-major Professor Daniel Nordman, Co-major Professor Dan Nettleton Huaiqing Wu Guang Song Iowa State University Ames, Iowa 2014 Copyright c Chuanlong Du, 2014. All rights reserved. ii DEDICATION I would like to dedicate this dissertation to my mother Jinyan Xu. She always encourages me to follow my heart and to seek my dreams. Her words have always inspired me and encouraged me to face difficulties and challenges I came across in my life. I would also like to dedicate this dissertation to my wife Lisha Li without whose support and help I would not have been able to complete this work. -
Chapter 0 Preliminaries
Chapter 0 Preliminaries Objective of the course Introduce basic matrix results and techniques that are useful in theory and applications. It is important to know many results, but there is no way I can tell you all of them. I would like to introduce techniques so that you can obtain the results you want. This is the standard teaching philosophy of \It is better to teach students how to fish instead of just giving them fishes.” Assumption You are familiar with • Linear equations, solution sets, elementary row operations. • Matrices, column space, row space, null space, ranks, • Determinant, eigenvalues, eigenvectors, diagonal form. • Vector spaces, basis, change of bases. • Linear transformations, range space, kernel. • Inner product, Gram Schmidt process for real matrices. • Basic operations of complex numbers. • Block matrices. Notation • Mn(F), Mm;n(F) are the set of n × n and m × n matrices over F = R or C (or a more general field). n • F is the set of column vectors of length n with entries in F. • Sometimes we write Mn;Mm;n if F = C. A1 0 • If A = 2 Mm+n(F) with A1 2 Mm(F) and A2 2 Mn(F), we write A = A1 ⊕ A2. 0 A2 • xt;At denote the transpose of a vector x and a matrix A. • For a complex matrix A, A denotes the matrix obtained from A by replacing each entry by its complex conjugate. Furthermore, A∗ = (A)t. More results and examples on block matrix multiplications. 1 1 Similarity, Jordan form, and applications 1.1 Similarity Definition 1.1.1 Two matrices A; B 2 Mn are similar if there is an invertible S 2 Mn such that S−1AS = B, equivalently, A = T −1BT with T = S−1. -
Sometimes Only Square Matrices Can Be Diagonalized 19
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 52, October 1975 SOMETIMESONLY SQUARE MATRICESCAN BE DIAGONALIZED LAWRENCE S. LEVY ABSTRACT. It is proved that every SQUARE matrix over a serial ring is equivalent to some diagonal matrix, even though there are rectangular matri- ces over these rings which cannot be diagonalized. Two matrices A and B over a ring R are equivalent if B - PAQ for in- vertible matrices P and Q over R- Following R. B. Warfield, we will call a module serial if its submodules are totally ordered by inclusion; and we will call a ring R serial it R is a direct sum of serial left modules, and also a direct sum of serial right mod- ules. When R is artinian, these become the generalized uniserial rings of Nakayama LE-GJ- Note that, in serial rings, finitely generated ideals do not have to be principal; for example, let R be any full ring of lower triangular matrices over a field. Moreover, it is well known (and easily proved) that if aR + bR is a nonprincipal right ideal in any ring, the 1x2 matrix [a b] is not equivalent to a diagonal matrix. Thus serial rings do not have the property that every rectangular matrix can be diagonalized. 1. Background. Serial rings are clearly semiperfect (P/radR is semisim- ple artinian and idempotents can be lifted modulo radP). (1.1) Every indecomposable projective module over a semiperfect ring R is S eR for some primitive idempotent e of R. (See [M].) A double-headed arrow -» will denote an "onto" homomorphism. -
Workshop on Applied Matrix Positivity International Centre for Mathematical Sciences 19Th – 23Rd July 2021
Workshop on Applied Matrix Positivity International Centre for Mathematical Sciences 19th { 23rd July 2021 Timetable All times are in BST (GMT +1). 0900 1000 1100 1600 1700 1800 Monday Berg Bhatia Khare Charina Schilling Skalski Tuesday K¨ostler Franz Fagnola Apanasovich Emery Cuevas Wednesday Jain Choudhury Vishwakarma Pascoe Pereira - Thursday Sharma Vashisht Mishra - St¨ockler Knese Titles and abstracts 1. New classes of multivariate covariance functions Tatiyana Apanasovich (George Washington University, USA) The class which is refereed to as the Cauchy family allows for the simultaneous modeling of the long memory dependence and correlation at short and intermediate lags. We introduce a valid parametric family of cross-covariance functions for multivariate spatial random fields where each component has a covariance function from a Cauchy family. We present the conditions on the parameter space that result in valid models with varying degrees of complexity. Practical implementations, including reparameterizations to reflect the conditions on the parameter space will be discussed. We show results of various Monte Carlo simulation experiments to explore the performances of our approach in terms of estimation and cokriging. The application of the proposed multivariate Cauchy model is illustrated on a dataset from the field of Satellite Oceanography. 1 2. A unified view of covariance functions through Gelfand pairs Christian Berg (University of Copenhagen) In Geostatistics one examines measurements depending on the location on the earth and on time. This leads to Random Fields of stochastic variables Z(ξ; u) indexed by (ξ; u) 2 2 belonging to S × R, where S { the 2-dimensional unit sphere { is a model for the earth, and R is a model for time. -
Doubling Algorithms with Permuted Lagrangian Graph Bases
DOUBLING ALGORITHMS WITH PERMUTED LAGRANGIAN GRAPH BASES VOLKER MEHRMANN∗ AND FEDERICO POLONI† Abstract. We derive a new representation of Lagrangian subspaces in the form h I i Im ΠT , X where Π is a symplectic matrix which is the product of a permutation matrix and a real orthogonal diagonal matrix, and X satisfies 1 if i = j, |Xij | ≤ √ 2 if i 6= j. This representation allows to limit element growth in the context of doubling algorithms for the computation of Lagrangian subspaces and the solution of Riccati equations. It is shown that a simple doubling algorithm using this representation can reach full machine accuracy on a wide range of problems, obtaining invariant subspaces of the same quality as those computed by the state-of-the-art algorithms based on orthogonal transformations. The same idea carries over to representations of arbitrary subspaces and can be used for other types of structured pencils. Key words. Lagrangian subspace, optimal control, structure-preserving doubling algorithm, symplectic matrix, Hamiltonian matrix, matrix pencil, graph subspace AMS subject classifications. 65F30, 49N10 1. Introduction. A Lagrangian subspace U is an n-dimensional subspace of C2n such that u∗Jv = 0 for each u, v ∈ U. Here u∗ denotes the conjugate transpose of u, and the transpose of u in the real case, and we set 0 I J = . −I 0 The computation of Lagrangian invariant subspaces of Hamiltonian matrices of the form FG H = H −F ∗ with H = H∗,G = G∗, satisfying (HJ)∗ = HJ, (as well as symplectic matrices S, satisfying S∗JS = J), is an important task in many optimal control problems [16, 25, 31, 36]. -
The Intersection of Some Classical Equivalence Classes of Matrices Mark Alan Mills Iowa State University
Iowa State University Capstones, Theses and Retrospective Theses and Dissertations Dissertations 1999 The intersection of some classical equivalence classes of matrices Mark Alan Mills Iowa State University Follow this and additional works at: https://lib.dr.iastate.edu/rtd Part of the Mathematics Commons Recommended Citation Mills, Mark Alan, "The intersection of some classical equivalence classes of matrices " (1999). Retrospective Theses and Dissertations. 12153. https://lib.dr.iastate.edu/rtd/12153 This Dissertation is brought to you for free and open access by the Iowa State University Capstones, Theses and Dissertations at Iowa State University Digital Repository. It has been accepted for inclusion in Retrospective Theses and Dissertations by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. INFORMATION TO USERS This manuscript has been reproduced from the microfilm master. UMI films the text directly fi'om the original or copy submitted. Thus, some thesis and dissertation copies are in typewriter face, while others may be from any type of computer printer. The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleedthrough, substandard margins, and improper alignment can adversely affect reproduction. In the unlikely event that the author did not send UMI a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion. Oversize materials (e.g., maps, drawings, charts) are reproduced by sectioning the original, beginning at the upper left-hand comer and continuing firom left to right in equal sections with small overlaps. -
A Review of Matrix Scaling and Sinkhorn's Normal Form for Matrices
A review of matrix scaling and Sinkhorn’s normal form for matrices and positive maps Martin Idel∗ Zentrum Mathematik, M5, Technische Universität München, 85748 Garching Abstract Given a nonnegative matrix A, can you find diagonal matrices D1, D2 such that D1 AD2 is doubly stochastic? The answer to this question is known as Sinkhorn’s theorem. It has been proved with a wide variety of methods, each presenting a variety of possible generalisations. Recently, generalisations such as to positive maps between matrix algebras have become more and more interesting for applications. This text gives a review of over 70 years of matrix scaling. The focus lies on the mathematical landscape surrounding the problem and its solution as well as the generalisation to positive maps and contains hardly any nontrivial unpublished results. Contents 1. Introduction3 2. Notation and Preliminaries4 3. Different approaches to equivalence scaling5 3.1. Historical remarks . .6 3.2. The logarithmic barrier function . .8 3.3. Nonlinear Perron-Frobenius theory . 12 3.4. Entropy optimisation . 14 3.5. Convex programming and dual problems . 17 3.6. Topological (non-constructive) approaches . 20 arXiv:1609.06349v1 [math.RA] 20 Sep 2016 3.7. Other ideas . 21 3.7.1. Geometric proofs . 22 3.7.2. Other direct convergence proofs . 23 ∗[email protected] 1 4. Equivalence scaling 24 5. Other scalings 27 5.1. Matrix balancing . 27 5.2. DAD scaling . 29 5.3. Matrix Apportionment . 33 5.4. More general matrix scalings . 33 6. Generalised approaches 35 6.1. Direct multidimensional scaling . 35 6.2. Log-linear models and matrices as vectors . -
Structure-Preserving Flows of Symplectic Matrix Pairs
STRUCTURE-PRESERVING FLOWS OF SYMPLECTIC MATRIX PAIRS YUEH-CHENG KUO∗, WEN-WEI LINy , AND SHIH-FENG SHIEHz Abstract. We construct a nonlinear differential equation of matrix pairs (M(t); L(t)) that are invariant (Structure-Preserving Property) in the class of symplectic matrix pairs X12 0 IX11 (M; L) = S2; S1 X = [Xij ]1≤i;j≤2 is Hermitian ; X22 I 0 X21 where S1 and S2 are two fixed symplectic matrices. Furthermore, its solution also preserves de- flating subspaces on the whole orbit (Eigenvector-Preserving Property). Such a flow is called a structure-preserving flow and is governed by a Riccati differential equation (RDE) of the form > > W_ (t) = [−W (t);I]H [I;W (t) ] , W (0) = W0, for some suitable Hamiltonian matrix H . We then utilize the Grassmann manifolds to extend the domain of the structure-preserving flow to the whole R except some isolated points. On the other hand, the Structure-Preserving Doubling Algorithm (SDA) is an efficient numerical method for solving algebraic Riccati equations and nonlinear matrix equations. In conjunction with the structure-preserving flow, we consider two special classes of sym- plectic pairs: S1 = S2 = I2n and S1 = J , S2 = −I2n as well as the associated algorithms SDA-1 k−1 and SDA-2. It is shown that at t = 2 ; k 2 Z this flow passes through the iterates generated by SDA-1 and SDA-2, respectively. Therefore, the SDA and its corresponding structure-preserving flow have identical asymptotic behaviors. Taking advantage of the special structure and properties of the Hamiltonian matrix, we apply a symplectically similarity transformation to reduce H to a Hamiltonian Jordan canonical form J. -
Math 623: Matrix Analysis Final Exam Preparation
Mike O'Sullivan Department of Mathematics San Diego State University Spring 2013 Math 623: Matrix Analysis Final Exam Preparation The final exam has two parts, which I intend to each take one hour. Part one is on the material covered since the last exam: Determinants; normal, Hermitian and positive definite matrices; positive matrices and Perron's theorem. The problems will all be very similar to the ones in my notes, or in the accompanying list. For part two, you will write an essay on what I see as the fundamental theme of the course, the four equivalence relations on matrices: matrix equivalence, similarity, unitary equivalence and unitary similarity (See p. 41 in Horn 2nd Ed. He calls matrix equivalence simply \equivalence."). Imagine you are writing for a fellow master's student. Your goal is to explain to them the key ideas. Matrix equivalence: (1) Define it. (2) Explain the relationship with abstract linear transformations and change of basis. (3) We found a simple set of representatives for the equivalence classes: identify them and sketch the proof. Similarity: (1) Define it. (2) Explain the relationship with abstract linear transformations and change of basis. (3) We found a set of representatives for the equivalence classes using Jordan matrices. State the Jordan theorem. (4) The proof of the Jordan theorem can be broken into two parts: (1) writing the ambient space as a direct sum of generalized eigenspaces, (2) classifying nilpotent matrices. Sketch the proof of each. Explain the role of invariant spaces. Unitary equivalence (1) Define it. (2) Explain the relationship with abstract inner product spaces and change of basis. -
Linear Transformations on Matrices
TRANSACTIONS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 198, 1974 LINEAR TRANSFORMATIONSON MATRICES BY D. I. DJOKOVICO) ABSTRACT. The real orthogonal group 0(n), the unitary group U(n) and the symplec- tic group Sp (n) are embedded in a standard way in the real vector space of n x n real, complex and quaternionic matrices, respectively. Let F be a nonsingular real linear transformation of the ambient space of matrices such that F(G) C G where G is one of the groups mentioned above. Then we show that either F(x) ■ ao(x)b or F(x) = ao(x')b where a, b e G are fixed, x* is the transpose conjugate of the matrix x and o is an automorphism of reals, complexes and quaternions, respectively. 1. Introduction. In his survey paper [4], M. Marcus has stated seven conjec- tures. In this paper we shall study the last two of these conjectures. Let D be one of the following real division algebras R (the reals), C (the complex numbers) or H (the quaternions). As usual we assume that R C C C H. We denote by M„(D) the real algebra of all n X n matrices over D. If £ G D we denote by Ï its conjugate in D. Let A be the automorphism group of the real algebra D. If D = R then A is trivial. If D = C then A is cyclic of order 2 with conjugation as the only nontrivial automorphism. If D = H then the conjugation is not an isomorphism but an anti-isomorphism of H. -
The Sturm-Liouville Problem and the Polar Representation Theorem
THE STURM-LIOUVILLE PROBLEM AND THE POLAR REPRESENTATION THEOREM JORGE REZENDE Grupo de Física-Matemática da Universidade de Lisboa Av. Prof. Gama Pinto 2, 1649-003 Lisboa, PORTUGAL and Departamento de Matemática, Faculdade de Ciências da Universidade de Lisboa e-mail: [email protected] Dedicated to the memory of Professor Ruy Luís Gomes Abstract: The polar representation theorem for the n-dimensional time-dependent linear Hamiltonian system Q˙ = BQ + CP , P˙ = −AQ − B∗P , with continuous coefficients, states that, given two isotropic solu- tions (Q1,P1) and (Q2,P2), with the identity matrix as Wronskian, the formula Q2 = r cos ϕ, Q1 = r sin ϕ, holds, where r and ϕ are continuous matrices, det r =6 0 and ϕ is symmetric. In this article we use the monotonicity properties of the matrix ϕ arXiv:1006.5718v1 [math.CA] 29 Jun 2010 eigenvalues in order to obtain results on the Sturm-Liouville prob- lem. AMS Subj. Classification: 34B24, 34C10, 34A30 Key words: Sturm-Liouville theory, Hamiltonian systems, polar representation. 1. Introduction Let n = 1, 2,.... In this article, (., .) denotes the natural inner 1 product in Rn. For x ∈ Rn one writes x2 =(x, x), |x| =(x, x) 2 . If ∗ M is a real matrix, we shall denote M its transpose. Mjk denotes the matrix entry located in row j and column k. In is the identity n × n matrix. Mjk can be a matrix. For example, M can have the four blocks M11, M12, M21, M22. In a case like this one, if M12 = M21 =0, we write M = diag (M11, M22). -
DECOMPOSITION of SYMPLECTIC MATRICES INTO PRODUCTS of SYMPLECTIC UNIPOTENT MATRICES of INDEX 2∗ 1. Introduction. Decomposition
Electronic Journal of Linear Algebra, ISSN 1081-3810 A publication of the International Linear Algebra Society Volume 35, pp. 497-502, November 2019. DECOMPOSITION OF SYMPLECTIC MATRICES INTO PRODUCTS OF SYMPLECTIC UNIPOTENT MATRICES OF INDEX 2∗ XIN HOUy , ZHENGYI XIAOz , YAJING HAOz , AND QI YUANz Abstract. In this article, it is proved that every symplectic matrix can be decomposed into a product of three symplectic 0 I unipotent matrices of index 2, i.e., every complex matrix A satisfying AT JA = J with J = n is a product of three −In 0 T 2 matrices Bi satisfying Bi JBi = J and (Bi − I) = 0 (i = 1; 2; 3). Key words. Symplectic matrices, Product of unipotent matrices, Symplectic Jordan Canonical Form. AMS subject classifications. 15A23, 20H20. 1. Introduction. Decomposition of elements in a matrix group into products of matrices with a special nature (such as unipotent matrices, involutions and so on) is a popular topic studied by many scholars. In the n × n matrix k ring Mn(F ) over a field F , a unipotent matrix of index k refers to a matrix A satisfying (A − In) = 0. Fong and Sourour in [4] proved that every matrix in the group SLn(C) (the special linear group over complex field C) is a product of three unipotent matrices (without limitation on the index). Wang and Wu in [6] gave a further result that every matrix in the group SLn(C) is a product of four unipotent matrices of index 2. In particular, decompositions of symplectic matrices have drawn considerable attention from numerous 0 I scholars (see [1, 3]).