Arxiv:1612.07142V1 [Math.AT] 21 Dec 2016 and Agn Spaces

Total Page:16

File Type:pdf, Size:1020Kb

Arxiv:1612.07142V1 [Math.AT] 21 Dec 2016 and Agn Spaces CAYLEY TRANSFORM ON STIEFEL MANIFOLDS ENRIQUE MAC´IAS-VIRGOS,´ MAR´IA JOSE´ PEREIRA-SAEZ,´ AND DANIEL TANRE´ Abstract. We define a Cayley transform on Stiefel manifolds. Applications to the Lusternik-Schnirelmann category and optimisation problems are presented. Introduction Denote by K the algebra of either the real numbers R, the complex numbers C or the quaternions H. Let G(n) = O(n, K) be the Lie group of matrices A ∈ Kn×n such ∗ ∗ ¯t that AA = In, where A = A is the conjugate transpose. Depending on K this group corresponds to the orthogonal group O(n), the unitary group U(n) or the symplectic group Sp(n). n×n Let I ∈ K be the identity matrix. The classical Cayley transform, cI : Ω(I) → Ω(I) −1 n×n −1 is defined by cI (X) = (I − X)(I + X) with Ω(I) = X ∈ K | (I + X) exists . 2 This map satisfies the equality cI = id. Moreover cI induces a diffeomorphism between n×n ∗ the tangent space TI G(n) = {X ∈ K | X + X = 0} and Ω(I) ∩ G(n). This construction was generalized by A. G´omez-Tato and the first two authors [4] to any ∗ A ∈ G(n) as a map cA : Ω(A) → Ω(A ) defined by ∗ −1 ∗ ∗ cA(X) = (I − A X)(A + X) = cI (A X)A , (1) with Ω(A)= X ∈ Kn×n | (A + X)−1 exists . −1 ∗ In this case we have cA = cA and there is a diffeomorphism between the tangent space n×n ∗ ∗ ∗ TAG(n)= {X ∈ K | A X + X A = 0} and Ω(A ) ∩ G(n). In this work, we construct Cayley transforms on Stiefel manifolds. We first specify some conventions and notations in use in this paper and state our main results. n Let 0 ≤ k ≤ n. The compact Stiefel manifold On,k of orthonormal k-frames in K n×k ∗ is the set of matrices x ∈ K such that x x = Ik. This manifold appears also as the basis of the principal fibration ι ρ G(n − k) −→ G(n) −→ On,k, arXiv:1612.07142v1 [math.AT] 21 Dec 2016 B 0 where ι(B) = and ρ is the projection onto the last k columns. If A ∈ G(n) 0 Ik and x = ρ(A), we denote by ρ∗A : TAG(n) → TxOn,k the map induced between the tangent spaces. Date: March 25, 2018. 2010 Mathematics Subject Classification. Primary 43A85; Secondary 58C99; 55M30. Key words and phrases. Quaternionic Stiefel manifold; Cayley transform; Lusternik-Schnirelmann category; Optimisation. 1 2 E.MAC´IAS-VIRGOS,´ M.-J. PEREIRA-SAEZ,´ AND D. TANRE´ The next statement contains the existence and the main properties of a Cayley trans- form in Stiefel manifolds. T Theorem A. Let 0 ≤ k ≤ n and x = ∈ O with P ∈ Kk×k. We choose P n,k α T A = ∈ G(n). Then there exists a map β P A γ : TxOn,k → On,k A such that γ ◦ ρ∗A = ρ ◦ cA. Moreover we have the following properties. 1) The map γA is injective on the open subset x X −1 Γ = v = A ∈ TxO | (βX + P ) exists . Y n,k This subset Γx does not depend on the choice of A such that ρ(A)= x. Furthermore, A x if γ is injective on an open subset U ⊂ TxOn,k then we have U ⊂ Γ . A x 2) The map γ induces a diffeomorphism between Γ ⊂ TxOn,k and the open subset τ Ωx = ∈ O | (π + P ∗)−1 exists . π n,k An explicit formula for γA is given in Definition 1.3. Also, the expression of the inverse A −1 map (γ |Γx ) appears in Equations (12) and (13). ∗ As we said before for the group G(n), the Cayley transform cA : TAG(n) → Ω(A ) ∩ G(n) is a diffeomorphism. Therefore the Cayley open subset Ω(A∗)∩G(n) is contractible. This property cannot be extended as it stands in the case of a Stiefel manifold. However the image of the injectivity domain of a Cayley transform in On,k is contractible in On,k. x Theorem B. For every x ∈ On,k the open subset Ω is contractible in On,k. This property is a consequence of the existence of a local section (see Proposition 3.3) A x s : Ω → G(n) of the projection ρ: G(n) → On,k and the contractibility of the Cayley open subsets Ω(A∗) ∩ G(n). The contents of the paper are as follows. Section 1 contains the construction of the A Cayley transform γ : TxOn,k → On,k. The study of the injectivity of its derivative is done in Section 2 and the proofs of Theorems A and B occupy Section 3. Finally, Section 4 is devoted to applications of this construction to Lusternik-Schnirelmann cat- egory of the quaternionic Stiefel manifolds and to optimisation problems on real Stiefel manifolds. 1. Construction Let Kn be either the real vector space Rn, the complex vector space Cn or the quater- nionic vector space Hn (with the structure of a right H-vector space) endowed with the ∗ inner product hu, vi = u v. Let 0 ≤ k ≤ n. The compact Stiefel manifold On,k of n n×k ∗ orthonormal k-frames in K is the set of matrices x ∈ K such that x x = Ik. It is standard to denote On,k by Vn,k in the real case, Wn,k in the complex case and Xn,k in the quaternionic case. CAYLEYTRANSFORMONSTIEFELMANIFOLDS 3 T Usually we write x = ∈ O , with T ∈ K(n−k)×k and P ∈ Kk×k. The linear left P n,k 0 action of G(n) on On,k is transitive and the isotropy group of x0 = is isomorphic Ik to G(n−k). Therefore On,k is diffeomorphic to G(n)/G(n−k) and we have the principal ι ρ fibration G(n − k) −→ G(n) −→ On,k. Let x ∈ On,k. We complete x to a matrix A ∈ G(n) such that ρ(A) = x. The tangent space TI G(n) of the group G(n) at the identity is the set of skew-Hermitian (skew-symmetric in the real case) matrices and the tangent space TAG(n) at A equals A · TI G(n). On the other side, recall that Kn×k ∗ ∗ TxOn,k = {v ∈ | v x + x v = 0} = A · Tx0 On,k, 0 where x0 = = ρ(In). So each tangent vector v ∈ TxOn,k can be written as Ik X v = A , with X ∈ K(n−k)×k, Y ∈ Kk×k and Y + Y ∗ = 0. Grants to the principal Y fibration defining the Stiefel manifold, the tangent space TxOn,k can be identified to the ⊥ orthogonal (TAG(n − k)) of the image of the inclusion of G(n − k) in G(n), ⊥ 0 X ∗ (TAG(n − k)) = A | Y + Y = 0 =∼ TxO . −X∗ Y n,k X ∗ With this identification, the tangent space TxO = A | Y + Y = 0 is consid- n,k Y ered as a subspace of TAG(n) and we may apply the Cayley map cA of G(n) on it. From Equation (1) we have cA = RA∗ ◦ cI ◦ LA∗ , (2) where LA∗ and RA∗ denote as usual the left and right multiplications in a Lie group. ⊥ Thus, we have first to determine cI on the elements of (TI G(n − k)) . Lemma 1.1. Let X ∈ K(n−k)×k and Y ∈ Kk×k such that Y + Y ∗ = 0. Then the matrix ∗ Ik + X X + Y is invertible. 0 X Proof. The skew-symmetric matrix M = cannot have real eigenvalues, −X∗ Y In−k X then In + M = ∗ is invertible. In the following product, −X Ik + Y In−k 0 In−k X ∗ · (In + M)= ∗ , (3) X Ik 0 X X + Ik + Y the two factors on the left-hand side admit an inverse. So, the matrix on the right-hand side admits an inverse, and the results follows. We denote ∗ −1 b = (Ik + X X + Y ) . (4) 4 E.MAC´IAS-VIRGOS,´ M.-J. PEREIRA-SAEZ,´ AND D. TANRE´ 0 X Proposition 1.2. Let M = ∈ TI G(n), then we have −X∗ Y ∗ In−k − 2XbX −2Xb cI (M)= ∗ . (5) 2bX −Ik + 2b Proof. From Equation (3) we have −1 −1 In−k X In−k 0 In−k −Xb In−k 0 (In + M) = −1 ∗ = ∗ . 0 b X Ik 0 b X Ik By applying the definition of cI : TI G(n) → G(n), we get: −1 cI (M) = (In − M) (In + M) ∗ In−k −X −1 In−k − 2XbX −2Xb = ∗ (In + M) = ∗ . X Ik − Y 2bX −Ik + 2b A computation from (2) and (5) gives directly: ∗ ∗ ∗ ∗ ∗ ∗ (In−k − 2XbX )α − 2XbT (In−k − 2XbX )β − 2XbP cA (AM)= ∗ ∗ ∗ ∗ ∗ ∗ . 2bX α + (−Ik + 2b)T 2bX β + (−Ik + 2b)P The Cayley transform is now obtained by projecting this expression on On,k. T α T Definition 1.3. Let x = ∈ O . We choose A = ∈ G(n), and consider P n,k β P X A v = A ∈ TxO . The Cayley transform on the Stiefel manifold, γ : TxO → Y n,k n,k On,k, is defined by ∗ ∗ ∗ ∗ A (In−k − 2XbX )β − 2XbP −Xb ∗ β γ (v)= ∗ ∗ ∗ = 2 (βX + P ) + ∗ , (6) 2bX β + (−Ik + 2b)P b −P where b is given in Equation (4). Remark 1.4. The map γA depends on the choice of A such that ρ(A) = x. With the previous notation, the elements of G(n) that are sent on x are the matrices A•E := E 0 αE T X E∗X A = with E ∈ G(n−k).
Recommended publications
  • Continuous Cayley Transform
    Continuous Cayley Transform A. Jadczyk In complex analysis, the Cayley transform is a mapping w of the complex plane to itself, given by z − i w(z) = : z + i This is a linear fractional transformation of the type az + b z 7! cz + d given by the matrix a b 1 −i C = = : c d 1 i The determinant of this matrix is 2i; and, noticing that (1 + i)2 = 2i; we can 1 as well replace C by 1+i C in order to to have C of determinant one - thus in the group SL(2; C): The transformation itself does not change after such a replacement. The Cayley transform is singular at the point z = −i: It maps the upper complex half{plane fz : Im(z) > 0g onto the interior of the unit disk D = fw : jwj < 1g; and the lower half{plane fz : Im(z) < 0g; except of the singular point, onto the exterior of the closure of D: The real line fz : Im(z) = 0g is mapped onto the unit circle fw : jwj = 1g - the boundary of D: The singular point z = −i can be thought of as being mapped to the point at infinity of the Riemann sphere - the one{point compactification of the complex plane. By solving the equation w(z) = z we easily find that there are two fixed points (z1; z2) of the Cayley transform: p p ! 1 3 1 3 z = + + i − − ; 1 2 2 2 2 p p ! 1 3 1 3 z = − + i − + : 2 2 2 2 2 In order to interpolate between the identity mapping and the Cayley trans- form we will choose a one{parameter group of transformations in SL(2; C) shar- ing these two fixed points.
    [Show full text]
  • Characterization and Computation of Matrices of Maximal Trace Over Rotations
    Characterization and Computation of Matrices of Maximal Trace over Rotations Javier Bernal1, Jim Lawrence1,2 1National Institute of Standards and Technology, Gaithersburg, MD 20899, USA 2George Mason University, 4400 University Dr, Fairfax, VA 22030, USA javier.bernal,james.lawrence @nist.gov [email protected] { } Abstract The constrained orthogonal Procrustes problem is the least-squares problem that calls for a rotation matrix that optimally aligns two cor- responding sets of points in d dimensional Euclidean space. This − problem generalizes to the so-called Wahba’s problem which is the same problem with nonnegative weights. Given a d d matrix M, × solutions to these problems are intimately related to the problem of finding a d d rotation matrix U that maximizes the trace of UM, × i.e., that makes UM a matrix of maximal trace over rotation matrices, and it is well known this can be achieved with a method based on the computation of the singular value decomposition (SVD) of M. As the main goal of this paper, we characterize d d matrices of maximal trace × over rotation matrices in terms of their eigenvalues, and for d = 2, 3, we show how this characterization can be used to determine whether a matrix is of maximal trace over rotation matrices. Finally, although depending only slightly on the characterization, as a secondary goal of the paper, for d = 2, 3, we identify alternative ways, other than the SVD, of obtaining solutions to the aforementioned problems. MSC: 15A18, 15A42, 65H17, 65K99, 93B60 Keywords: Eigenvalues, orthogonal, Procrustes, rotation, singular value decomposition, trace, Wahba 1 Contents 1 Introduction 2 2 Reformulation of Problems as Maximizations of the Trace of Matrices Over Rotations 4 3 Characterization of Matrices of Maximal Trace Over Rota- tions 6 4 The Two-Dimensional Case: Computation without SVD 17 5 The Three-Dimensional Case: Computation without SVD 22 6 The Three-Dimensional Case Revisited 28 Summary 33 References 34 1 Introduction Suppose P = p1,...,pn and Q = q1,...,qn are each sets of n points in Rd.
    [Show full text]
  • Rotation Matrix - Wikipedia, the Free Encyclopedia Page 1 of 22
    Rotation matrix - Wikipedia, the free encyclopedia Page 1 of 22 Rotation matrix From Wikipedia, the free encyclopedia In linear algebra, a rotation matrix is a matrix that is used to perform a rotation in Euclidean space. For example the matrix rotates points in the xy -Cartesian plane counterclockwise through an angle θ about the origin of the Cartesian coordinate system. To perform the rotation, the position of each point must be represented by a column vector v, containing the coordinates of the point. A rotated vector is obtained by using the matrix multiplication Rv (see below for details). In two and three dimensions, rotation matrices are among the simplest algebraic descriptions of rotations, and are used extensively for computations in geometry, physics, and computer graphics. Though most applications involve rotations in two or three dimensions, rotation matrices can be defined for n-dimensional space. Rotation matrices are always square, with real entries. Algebraically, a rotation matrix in n-dimensions is a n × n special orthogonal matrix, i.e. an orthogonal matrix whose determinant is 1: . The set of all rotation matrices forms a group, known as the rotation group or the special orthogonal group. It is a subset of the orthogonal group, which includes reflections and consists of all orthogonal matrices with determinant 1 or -1, and of the special linear group, which includes all volume-preserving transformations and consists of matrices with determinant 1. Contents 1 Rotations in two dimensions 1.1 Non-standard orientation
    [Show full text]
  • Generalization of Rotational Mechanics and Application to Aerospace Systems
    CORE Metadata, citation and similar papers at core.ac.uk Provided by Texas A&amp;M Repository GENERALIZATION OF ROTATIONAL MECHANICS AND APPLICATION TO AEROSPACE SYSTEMS A Dissertation by ANDREW JAMES SINCLAIR Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY May 2005 Major Subject: Aerospace Engineering GENERALIZATION OF ROTATIONAL MECHANICS AND APPLICATION TO AEROSPACE SYSTEMS A Dissertation by ANDREW JAMES SINCLAIR Submitted to Texas A&M University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Approved as to style and content by: John L. Junkins John E. Hurtado (Co-Chair of Committee) (Co-Chair of Committee) Srinivas R. Vadali John H. Painter (Member) (Member) Helen L. Reed (Head of Department) May 2005 Major Subject: Aerospace Engineering iii ABSTRACT Generalization of Rotational Mechanics and Application to Aerospace Systems. (May 2005) Andrew James Sinclair, B.S., University of Florida; M.S., University of Florida Co–Chairs of Advisory Committee: Dr. John L. Junkins Dr. John E. Hurtado This dissertation addresses the generalization of rigid-body attitude kinematics, dynamics, and control to higher dimensions. A new result is developed that demon- strates the kinematic relationship between the angular velocity in N-dimensions and the derivative of the principal-rotation parameters. A new minimum-parameter de- scription of N-dimensional orientation is directly related to the principal-rotation parameters. The mapping of arbitrary dynamical systems into N-dimensional rotations and the merits of new quasi velocities associated with the rotational motion are studied.
    [Show full text]
  • Medical Issue2012 4
    Journal of Science and Technology Vol. 13 June /2012 ISSN -1605 – 427X Natural and Medical Sciences (NMS No. 1) www.sustech.edu On the Development and Recent Work of Cayley Transforms [I] Adam Abdalla Abbaker Hassan Nile Valley University, Faculty of Education- Department of Mathematics ABSTRACT: This paper deals with the recent work and developments of Cayely transforms which had been rapidly developed. ا : : ﺘﺘﻌﻠﻕ ﻫﺫﻩ ﺍﻝﻭﺭﻗﺔ ﺍﻝﻌﻠﻤﻴﺔ ﺒﺘﻘﺩﻴﻡ ﺍﻝﻌﻤل ﺍﻝﺤﺩﻴﺙ ﻭ ﺍﻝﺘﻁﻭﺭﺍﺕ ﺍﻝﺴﺭﻴﻌﺔ ﻝﺘﺤﻭﻴﻼﺕ ﻜﻤ ﺎﻴﻠﻲ . INTRODUCTION In this study we shall give the basic is slightly simplified variant of that of important transforms and historical Langer (1, 2) . background of the subject, then introduce -1 the accretive operator of Cayley transform Cayley Transforms V = ( A – iI)(A+iI) and a maximal accretive and dissipative First we introduce the transformation B = operator, and relate this with Cayley -1 ( I + T*T) and Transform which yield the homography in function calculus. After this we will study C = T ( I + T*T) -1 in order to study the the fractional power and maximal Cayely transforms. accretive operators; indeed all this had been developed by Sz.Nagy and If the linear transform T is bounded , it is C. Foias (1) . The method of extending and clear that the transformation B appearing accretive (or dissipative) operator to above is also bounded , symmetric , and a maximal one via Cayley transforms such that 0 ≤ B ≤ I ; (3) . modeled on Von Neumann’s theory on symmetric operators, is due to Philips (1) C = TB is then bounded too. If T is a fractional powers of operators A in Hilbert linear transformation with dense domain, space , or even in Banachh space, such as we know that T*, and consequently also T*T, exist, but we know nothing of their – A is infinitesimal generator of a conts (4) one parameter sem–group of contractions, domains of definition .
    [Show full text]
  • Complex Variables Visualized
    Diplomarbeit Revised Version - June 17, 2014 Computer Algebra and Analysis: Complex Variables Visualized Ausgef¨uhrtam Research Institute for Symbolic Computation (RISC) der Johannes Kepler Universit¨atLinz unter Anleitung von Univ.-Prof. Dr. Peter Paule durch Thomas Ponweiser Stockfeld 15 4283 Bad Zell Datum Unterschrift ii Kurzfassung Ziel dieser Diplomarbeit ist die Visualisierung einiger grundlegender Ergeb- nisse aus dem Umfeld der Theorie der modularen Gruppe sowie der modularen Funktionen unter Zuhilfenahme der Computer Algebra Software Mathematica. Die Arbeit gliedert sich in drei Teile. Im ersten Kapitel werden f¨urdiese Ar- beit relevante Begriffe aus der Gruppentheorie zusammengefasst. Weiters wer- den M¨obiusTransformationen eingef¨uhrtund deren grundlegende geometrische Eigenschaften untersucht. Das zweite Kapitel ist der Untersuchung der modularen Gruppe aus al- gebraischer und geometrischer Sicht gewidmet. Der kanonische Fundamen- talbereich der modularen Gruppe sowie die daraus abgeleitete Kachelung der oberen Halbebene wird eingef¨uhrt.Weiters wird eine generelle Methode zum Auffinden von Fundamentalbereichen f¨urUntergruppen der modularen Gruppe vorgestellt, welche sich zentral auf die Konzepte der 2-dimensionalen hyper- bolischen Geometrie st¨utzt. Im dritten Kapitel geben wir einige konkrete Beispiele, wie die aufge- baute Theorie f¨urdie Visualisierung bestimmter mathematischer Sachverhalte angewendet werden kann. Neben der Visualisierung von Graphen modularer Funktionen stellt sich auch der Zusammenhang zwischen modularen Transfor- mationen und Kettenbr¨uchen als besonders sch¨onesErgebnis dar. iii iv Abstract The aim of this diploma thesis is the visualization of some fundamental results in the context of the theory of the modular group and modular functions. For this purpose the computer algebra software Mathematica is utilized. The thesis is structured in three parts.
    [Show full text]
  • Metaplectic Group, Symplectic Cayley Transform, and Fractional Fourier Transforms
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector J. Math. Anal. Appl. 416 (2014) 947–968 Contents lists available at ScienceDirect Journal of Mathematical Analysis and Applications www.elsevier.com/locate/jmaa Metaplectic group, symplectic Cayley transform, and fractional Fourier transforms MauriceA.deGossona,∗,FranzLuefb a University of Vienna, Faculty of Mathematics (NuHAG), 1090 Vienna, Austria b Department of Mathematics, Norwegian University of Science and Technology, 7491 Trondheim, Norway article info abstract Article history: We begin with a survey of the standard theory of the metaplectic group with some Received 15 February 2013 emphasis on the associated notion of Maslov index. We thereafter introduce the Available online 14 March 2014 Cayley transform for symplectic matrices, which allows us to study in detail the Submitted by P.G. Lemarie-Rieusset spreading functions of metaplectic operators, and to prove that they are basically Keywords: quadratic chirps. As a non-trivial application we give new formulae for the fractional Metaplectic and symplectic group Fourier transform in arbitrary dimension. We also study the regularity of the Fractional Fourier transform solutions to the Schrödinger equation in the Feichtinger algebra. Weyl operator © 2014 The Authors. Published by Elsevier Inc. This is an open access article Feichtinger algebra under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/). Contents 1. Introduction.....................................................................
    [Show full text]
  • Superconnection Character Forms and the Cayley Transform
    0040 9283 RX so3 00 + 00 Pergamon Press plc SUPERCONNECTION CHARACTER FORMS AND THE CAYLEY TRANSFORM DANIEL QUILLEN (Received June 1987) $0. INTRODUCTION LET E = E+ @E- be a super vector bundle equipped with inner product, and let D be a connection on E compatible with the grading and inner product. Let x= where T:E’-+E- is a vector bundle map. In [l l] we introduced the superconnection character form (1) O(E, D,X)=tr,exp u(X’+[D,X]~+D*) which generalizes the well-known definition of the Chern character in de Rham co- homology. This construction was motivated by the problem of proving a local version of the Atiyah-Singer index theorem for families of Dirac operators. The formula (1) is used,in an infinite dimensional context in Bismut’s solution of this problem [3]. From the cohomological viewpoint the differential form (1) is not very interesting. Its cohomology class remains unchanged as X is contracted to zero by tX, so the class doesn’t depend on X. In the first part of this paper we remedy this defect by showing that (1) has a meaning when X is allowed to become infinite in a certain sense. We show that this formula can be written in terms of the (ungraded) subbundle of E given by the graph of T. Moreover the resulting expression makes sense for an arbitrary subbundle W. Thus superconnection character forms can be defined for any correspondence between E+ and E-. We also show that the corresponding cohomology class is the Chern character of the difference W-E-.
    [Show full text]
  • Cayley Transformation and Numerical Stability of Calibration Equation
    Int J Comput Vis (2009) 82: 156–184 DOI 10.1007/s11263-008-0193-x Cayley Transformation and Numerical Stability of Calibration Equation F.C. Wu · Z.H. Wang · Z.Y. Hu Received: 25 May 2008 / Accepted: 22 October 2008 / Published online: 5 December 2008 © US Government 2008 Abstract The application of Cayley transformation to en- it is proved that the standard calibration equation, the Cay- hance the numerical stability of camera calibration is inves- ley calibration family and the S-Cayley calibration family tigated. First, a new calibration equation, called the standard are all some special cases of this generic calibration family. calibration equation, is introduced using the Cayley transfor- mation and its analytical solution is obtained. The standard Keywords Camera calibration · The absolute conic · calibration equation is equivalent to the classical calibration Calibration equation · Cayley transformation · Numerical equation, but it exhibits remarkable better numerical stabil- stability ity. Second, a one-parameter calibration family, called the Cayley calibration family which is equivalent to the stan- dard calibration equation, is obtained using also the Cayley 1 Introduction transformation and it is found that this family is composed Camera calibration is a necessary step to recover 3D met- of those infinite homographies whose rotation has the same ric information from 2D images. Many methods and tech- axis with the rotation between the two given views. The con- niques for camera calibration have been proposed in the last dition number of equations in the Cayley calibration family decades. The existing methods can be roughly classified into varies with the parameter value, and an algorithm to deter- two categories: methods based on a reference calibration mine the best parameter value is provided.
    [Show full text]
  • Spectral Theory of Unbounded Self-Adjoint Operators in Hilbert Spaces
    University of Helsinki Department of Mathematics and Statistics Pro Gradu Thesis Spectral Theory of Unbounded Self-adjoint Operators in Hilbert spaces Timo Weckman 013626675 Instructor: Jari Taskinen April 30, 2020 HELSINGIN YLIOPISTO HELSINGFORS UNIVERSITET UNIVERSITY OF HELSINKI Tiedekunta/Osasto Fakultet/Sektion Faculty Laitos Institution Department Faculty of Science Department of Mathematics and Statistics Tekijä Författare Author Timo Weckman Työn nimi Arbetets titel Title Spectral Theory of Unbounded Self-adjoint Operators in Hilbert spaces Oppiaine Läroämne Subject Mathematics Työn laji Arbetets art Level Aika Datum Month and year Sivumäärä Sidoantal Number of pages Pro gradu -thesis April 2020 44 p. Tiivistelmä Referat Abstract In this work we present a derivation of the spectral theorem of unbounded spectral operators in a Hilbert space. The spectral theorem has several applications, most notably in the theory of quantum mechanics. The theorem allows a self-adjoint linear operator on a Hilbert space to be represented in terms of simpler operators, projections. The focus of this work are the self-adjoint operators in Hilbert spaces. The work begins with the int- roduction of vector and inner product spaces and the denition of the complete inner product space, the Hilbert space. Three classes of bounded linear operators relevant for this work are introduced: self-adjoint, unitary and projection operators. The unbounded self-adjoint operator and its proper- ties are also discussed. For the derivation of the spectral theorem, the basic spectral properties of operators in Hilbert space are presented. The spectral theorem is rst derived for bounded operators. With the denition of basic spectral properties and the introduction of the spectral family, the spectral theorem for bounded self-adjoint operators is presented with a proof.
    [Show full text]
  • Hybrid Attitude Control and Estimation on SO(3)
    Western University Scholarship@Western Electronic Thesis and Dissertation Repository 11-29-2017 4:00 PM Hybrid Attitude Control and Estimation On SO(3) Soulaimane Berkane The University of Western Ontario Supervisor Tayebi, Abdelhamid The University of Western Ontario Graduate Program in Electrical and Computer Engineering A thesis submitted in partial fulfillment of the equirr ements for the degree in Doctor of Philosophy © Soulaimane Berkane 2017 Follow this and additional works at: https://ir.lib.uwo.ca/etd Part of the Acoustics, Dynamics, and Controls Commons, Controls and Control Theory Commons, and the Navigation, Guidance, Control and Dynamics Commons Recommended Citation Berkane, Soulaimane, "Hybrid Attitude Control and Estimation On SO(3)" (2017). Electronic Thesis and Dissertation Repository. 5083. https://ir.lib.uwo.ca/etd/5083 This Dissertation/Thesis is brought to you for free and open access by Scholarship@Western. It has been accepted for inclusion in Electronic Thesis and Dissertation Repository by an authorized administrator of Scholarship@Western. For more information, please contact [email protected]. Abstract This thesis presents a general framework for hybrid attitude control and estimation de- sign on the Special Orthogonal group SO(3). First, the attitude stabilization problem on SO(3) is considered. It is shown that, using a min-switch hybrid control strategy de- signed from a family of potential functions on SO(3), global exponential stabilization on SO(3) can be achieved when this family of potential functions satisfies certain properties. Then, a systematic methodology to construct these potential functions is developed. The proposed hybrid control technique is applied to the attitude tracking problem for rigid body systems.
    [Show full text]
  • Math 512 / Problem Set 11 (Three Pages) • Study/Read: Inner Product
    Due: Wednesday, December 3, at 4 PM Math 512 / Problem Set 11 (three pages) • Study/read: Inner product spaces & Operators on inner product spaces - Ch. 5 & 6 from LADW by Treil. - Ch. 8 & 9 from Linear Algebra by Hoffman & Kunze. (!) Make sure that you understand perfectly the definitions, examples, theorems, etc. R 1 Legendre Polynomials Recall that (f, g) = −1 f(x)g(x)dx is an inner product on C(I, R), where I = [−1, 1]. The system of polynomials (L0(t), L1(t), ... , Ln−1(t)) resulting from the n−1 Gram-Schmidt orthonormalization of the standard basis E := (1, t, ... , t ) of Pn(R) are called the Legendre polynomials. Google the term “Legendre polynomials” and learn about their applications. 1) Answer / do the following: a) Compute L0, L1, L2, L3. 3 b) Write t − t + 1 as a linear combination of L0, L1, L2, L3. c) Show that in general, one has deg(Lk) = k. k d) What is the coefficient of t in Lk(t)? 2) Which of the following pairs of matrices are unitarily equivalent: 0 1 0 2 0 0 0 1 0 1 0 0 1 0! 0 1! a) & b) −1 0 0 & 0 −1 0 c) −1 0 0 & 0 −i 0 0 1 1 0 0 0 1 0 0 0 0 0 1 0 0 i [Hint: Unitarily equivalent matrices have the same characteristic polynomial (WHY), thus the same eigenvalues including multiplicities, the same trace, determinant, etc.] 3) For each of the following matrices over F = R decide whether they are unitarily diago- nalizable, and if so, find the unitary matrix P over R such that P AP −1 is diagonal.
    [Show full text]