Hybrid Attitude Control and Estimation on SO(3)

Total Page:16

File Type:pdf, Size:1020Kb

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. A smoothing mechanism is proposed to filter out the discrete behaviour of the hybrid switching mechanism leading to control torques that are continuous. Next, the problem of attitude estimation from continuous body-frame vector mea- surements of known inertial directions is considered. Two hybrid attitude and gyro bias observers designed directly on SO(3) × R3 are proposed. The first observer uses a set of innovation terms and a switching mechanism that selects the appropriate innovation term. The second observer uses a fixed innovation term and allows the attitude state to be reset (experience discrete transition or jump) to an adequately chosen value on SO(3). Both hybrid observers guarantee global exponential stability of the zero estimation errors. Finally, in the case where the body-frame vector measurements are intermittent, an event-triggered attitude estimation scheme on SO(3) is proposed. The observer consists in integrating the continuous angular velocity during the interval of time where the vector measurements are not available, and updating the attitude state upon the arrival of the vector measurements. Both cases of synchronous and asynchronous vector measurements with possible irregular sampling periods are considered. Moreover, some modifications to the intermittent observer are developed to handle different practical issues such as discrete-time implementation, noise filtering and gyro bias compensation. ii To all those whom I love... iii Acknowledgement I would like to express my deepest gratitude to my advisor, Professor Abdelhamid Tayebi, for his support, motivation and guidance throughout the four years I have spent at Western University. He taught me what it meant to be a high quality researcher who carries fundamental values such as honesty, vision, passion and rigorousness. His timely feedback and suggestions helped me to smoothly conduct my research and write my thesis. I am grateful to my friend, colleague and collaborator, Dr. Abdelkader Abdessameud, with whom I have completed several parts of this research and shared countless discus- sions in the research field. I am also thankful to my alternate supervisor at Western, Dr. Ilia Polushin, who have helped me in several occasions. I would also like to sincerely thank Prof. Jin Jiang and Dr. Lyndon Brown for taking the time to serve as my PhD thesis examiners and for their constructive comments and feedback. I would especially like to thank Prof. Andrew R. Teel from the University of California Santa Barbara for agreeing to review and examine my thesis. I have had several fruitful discussions with him during our meetings at CDC and ACC conferences. Lastly, I would like to extend my sincere appreciation to my family. To my mother, Chahrazed, my sisters and my friends for their support and prayers in all my pursuits. And of course, to my wife, Intissar, for her love, support and understanding. iv List of Symbols • N and N>0 denote the natural and strictly positive natural numbers, respectively. • R, R≥0 and R>0 denote the real, nonnegative and positive real numbers, respec- tively. • Rn is the n-dimensional Euclidean space . • Rn×m is the set of real-valued n × m matrices. • Given A 2 Rn×m, A> denotes its transpose. • Given A 2 Rn×n, det(A) denotes its determinant. • Given A 2 Rn×n, tr(A) denotes the sum of its diagonal entries (trace). • Given A; B 2 Rn×m, their Euclidean inner product is defined as hhA; Bii = tr(A>B). n×m p • Given A 2 R , its Frobenius norm is kAkF = hhA; Aii. p • Given x 2 Rn, its Euclidean norm (2-norm) is given by kxk = x>x. • Sn = fx 2 Rn+1 : kxk = 1g is the unit n-sphere embedded in Rn+1. • B = fx 2 Rn : kxk ≤ 1g is the closed unit ball in Euclidean space. • SO(3) = fR 2 R3×3 : det(R) = 1; RR> = R>R = Ig is the Special Orthogonal group of order 3 where I denotes the three-dimensional identity matrix. • so(3) = fΩ 2 R3×3 :Ω> = −Ωg is the set of all skew-symmetric 3 × 3 matrices and defines also the Lie algebra of SO(3). • Q = f(η; ) 2 R × R3 : η2 + > = 1g is the set of unit quaternions. p • Given a rotation matrix R 2 SO(3), let us define jRjI = kI − Rk= 8. 3×3 3×3 1 > • The map E : R ! R is defined as E(M) = 2 (tr(M)I − M ). 2 • The map Ra : R × S ! SO(3) is defined as Ra(θ; u) = I + sin(θ)[u]× + (1 − 2 cos(θ))[u]×. 2 • The map Ru : Q ! SO(3) is defined as Ru(η; ) = I + 2[]× + 2η[]×. v 3 1 2 > • The map Rr : R ! SO(3) is defined as Rr(z) = 1+kzk2 (1 − kzk )I + 2zz + 2[z]× . 3×3 3×3 • The map PSO(3) : R ! SO(3) is defined as the projection of X 2 R on SO(3) (i.e. closest rotation matrix to X). 3×3 > • The map Pso(3) : R ! so(3) is defined as Pso(3)(M) = (M − M )=2. 3 3 • The map [·]× : R ! so(3) is defined as [x]×y = x × y for all x; y 2 R where × denotes the cross product on R3. 3 • The map vex : so(3) ! R is the inverse map of [·]×. 3×3 3 • The map : R ! R is defined as = vex ◦ Pso(3) where the symbol ◦ is used to denote function composition. • The set ΠSO(3) = fR 2 SO(3) : tr(R) 6= −1g is the set of all rotations with an angle ◦ 2 different than 180 . In other words, ΠSO(3) = SO(3) n Ra(π; S ). • The set ΠQ = f(η; ) 2 Q : η 6= 0g is the double cover of ΠSO(3) in Q. 3 3 3 • The map Pc : R × R ! R , for a given c 2 R≥0, defines the parameter projection function which is given in (2.3). • Given a set S, cl(S) denotes its closure (S together with all of its limit points). • A set-valued map F : M ⇒ N assigns to every x 2 M a set of values F(x) ⊆ N . • Let F : M ⇒ N be a set-valued function. The domain of F is defined as the set domF = fx 2 M : F(x) 6= ;g where ; denotes the empty set. • Cn(M; N ) is the set of all functions f : M!N such that the first n 2 N derivatives of each function f 2 Cn(M; N ) exist and are continuous. vi Contents Abstract ii Acknowlegements iv List of Symbols v List of Figures xi 1 Introduction 1 1.1 General Introduction . 1 1.2 Attitude Control . 2 1.3 Attitude Estimation . 5 1.4 Thesis Contributions . 7 1.5 List of Publications . 10 1.6 Thesis Outline . 12 2 Background and Preliminaries 14 2.1 General Notations . 14 2.2 Rigid Body Attitude . 16 2.2.1 Attitude Parametrizations . 18 2.2.1.1 Exponential Coordinates Representation . 18 2.2.1.2 Angle-Axis Representation . 19 2.2.1.3 Unit Quaternions Representation . 20 2.2.1.4 Rodrigues Vector Representation . 21 2.2.2 Metrics on SO(3) . 21 2.2.3 Attitude Visualization . 23 2.2.4 Useful Identities and Lemmas . 24 2.3 Hybrid Systems Framework . 28 2.3.1 Exponential Stability for Hybrid Systems . 30 2.4 Numerical Integration Tools . 32 vii 2.4.1 Numerical Integration on Rn ..................... 33 2.4.2 Numerical Integration on SO(3) . 34 2.4.3 Numerical Integration of Hybrid Systems . 35 3 Hybrid Attitude Control on SO(3) 37 3.1 Introduction . 37 3.2 Motivation Using Planar Rotations on S1 . 38 3.3 Attitude Stabilization on SO(3) . 47 3.3.1 Problem Formulation . 47 3.3.2 Smooth Attitude Stabilization on SO(3) . 48 3.3.3 Synergistic and Exp-Synergistic Potential Functions . 51 3.3.4 Hybrid Attitude Stabilization on SO(3) . 53 3.3.5 Construction of Exp-Synergistic Potential Functions . 55 3.3.6 Illustration of the Switching Mechanism . 64 3.4 Attitude Tracking for Rigid Body Systems .
Recommended publications
  • Rotational Motion (The Dynamics of a Rigid Body)
    University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Robert Katz Publications Research Papers in Physics and Astronomy 1-1958 Physics, Chapter 11: Rotational Motion (The Dynamics of a Rigid Body) Henry Semat City College of New York Robert Katz University of Nebraska-Lincoln, [email protected] Follow this and additional works at: https://digitalcommons.unl.edu/physicskatz Part of the Physics Commons Semat, Henry and Katz, Robert, "Physics, Chapter 11: Rotational Motion (The Dynamics of a Rigid Body)" (1958). Robert Katz Publications. 141. https://digitalcommons.unl.edu/physicskatz/141 This Article is brought to you for free and open access by the Research Papers in Physics and Astronomy at DigitalCommons@University of Nebraska - Lincoln. It has been accepted for inclusion in Robert Katz Publications by an authorized administrator of DigitalCommons@University of Nebraska - Lincoln. 11 Rotational Motion (The Dynamics of a Rigid Body) 11-1 Motion about a Fixed Axis The motion of the flywheel of an engine and of a pulley on its axle are examples of an important type of motion of a rigid body, that of the motion of rotation about a fixed axis. Consider the motion of a uniform disk rotat­ ing about a fixed axis passing through its center of gravity C perpendicular to the face of the disk, as shown in Figure 11-1. The motion of this disk may be de­ scribed in terms of the motions of each of its individual particles, but a better way to describe the motion is in terms of the angle through which the disk rotates.
    [Show full text]
  • The Rotation Group in Plate Tectonics and the Representation of Uncertainties of Plate Reconstructions
    Geophys. J. In?. (1990) 101, 649-661 The rotation group in plate tectonics and the representation of uncertainties of plate reconstructions T. Chang’ J. Stock2 and P. Molnar3 ‘Department of Mathematics, University of Virginia, Charlottesville, Virginia 22903-3199, USA ‘Department of Earth and Planetary Sciences, Harvard University, Cambridge, MA 02138, USA ’Department of Earth, Atmospheric, and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139, USA Accepted 1989 December 27. Received 1989 December 27; in original form 1989 August 11 Downloaded from SUMMARY The calculation of the uncertainty in an estimated rotation requires a parametriza- tion of the rotation group; that is, a unique mapping of the rotation group to a point in 3-D Euclidean space, R3. Numerous parametrizations of a rotation exist, http://gji.oxfordjournals.org/ including: (1) the latitude and longitude of the axis of rotation and the angle of rotation; (2) a representation as a Cartesian vector with length equal to the rotation angle and direction parallel to the rotation axis; (3) Euler angles; or (4) unit length quaternions (or, equivalently, Cayley-Klein parameters). The uncertainty in a rotation is determined by the effect of nearby rotations on the rotated data. The uncertainty in a rotation is small, if rotations close to the best fitting rotation degrade the fit of the data by a large amount, and it is large, if only rotations differing by a large amount cause such a degradation. Ideally, we would at California Institute of Technology on August 29, 2014 like to parametrize the rotations in such a way so that their representation as points in R3 would have the property that the distance between two points in R3 reflects the effects of the corresponding rotations on the fit of the data.
    [Show full text]
  • 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]
  • 4-2 Degrees and Radians
    4-2 Degrees and Radians Write each degree measure in radians as a multiple of π and each radian measure in degrees. 10. 30° SOLUTION: To convert a degree measure to radians, multiply by ANSWER: 14. SOLUTION: To convert a radian measure to degrees, multiply by ANSWER: 120° Identify all angles that are coterminal with the given angle. Then find and draw one positive and one negative angle coterminal with the given angle. 20. 225° SOLUTION: All angles measuring are coterminal with a angle. Sample answer: Let n = 1 and −1. eSolutions Manual - Powered by Cognero Page 1 ANSWER: 225° + 360n°; Sample answer: 585°, −135° 24. SOLUTION: All angles measuring are coterminal with a angle. Sample answer: Let n = 1 and −1. ANSWER: Find the length of the intercepted arc with the given central angle measure in a circle with the given radius. Round to the nearest tenth. 29. , r = 4 yd SOLUTION: ANSWER: 5.2 yd 30. 105°, r = 18.2 cm SOLUTION: Method 1 Convert 105° to radian measure, and then use s = rθ to find the arc length. Substitute r = 18.2 and θ = . Method 2 Use s = to find the arc length. ANSWER: 33.4 cm Find the rotation in revolutions per minute given the angular speed and the radius given the linear speed and the rate of rotation. 36. = 104π rad/min SOLUTION: The angular speed is 104π radians per minute. Each revolution measures 2π radians 104π ÷ 2π = 52 The angle of rotation is 52 revolutions per minute. ANSWER: 52 rev/min 37. v = 82.3 m/s, 131 rev/min SOLUTION: 131 × 2π = 262π The linear speed is 82.3 meters per second with an angle of rotation of 262π radians per minute.
    [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]
  • Rotational Motion CONTENTS CHAPTER-OPENING QUESTION—Guess Now! 8–1 Angular Quantities a Solid Ball and a Solid Cylinder Roll Down a Ramp
    You too can experience rapid rotation—if your stomach can take the high angular velocity and centripetal acceleration of some of the faster amusement park rides. If not, try the slower merry-go-round or Ferris wheel. Rotating carnival rides have rotational kinetic energy as well as angular momentum. Angular acceleration is produced by a net torque, and rotating objects have rotational kinetic energy. P T A E H R C 8 Rotational Motion CONTENTS CHAPTER-OPENING QUESTION—Guess now! 8–1 Angular Quantities A solid ball and a solid cylinder roll down a ramp. They both start from rest at the 8–2 Constant Angular Acceleration same time and place. Which gets to the bottom first? 8–3 Rolling Motion (a) They get there at the same time. (Without Slipping) (b) They get there at almost exactly the same time except for frictional differences. 8–4 Torque (c) The ball gets there first. 8–5 Rotational Dynamics; (d) The cylinder gets there first. Torque and Rotational Inertia (e) Can’t tell without knowing the mass and radius of each. 8–6 Solving Problems in Rotational Dynamics ntil now, we have been concerned mainly with translational motion. We 8–7 Rotational Kinetic Energy discussed the kinematics and dynamics of translational motion (the role 8–8 Angular Momentum and of force). We also discussed the energy and momentum for translational Its Conservation U *8–9 Vector Nature of motion. In this Chapter we will deal with rotational motion. We will discuss the Angular Quantities kinematics of rotational motion and then its dynamics (involving torque), as well as rotational kinetic energy and angular momentum (the rotational analog of linear momentum).
    [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]
  • ROTATION: a Review of Useful Theorems Involving Proper Orthogonal Matrices Referenced to Three- Dimensional Physical Space
    Unlimited Release Printed May 9, 2002 ROTATION: A review of useful theorems involving proper orthogonal matrices referenced to three- dimensional physical space. Rebecca M. Brannon† and coauthors to be determined †Computational Physics and Mechanics T Sandia National Laboratories Albuquerque, NM 87185-0820 Abstract Useful and/or little-known theorems involving33× proper orthogonal matrices are reviewed. Orthogonal matrices appear in the transformation of tensor compo- nents from one orthogonal basis to another. The distinction between an orthogonal direction cosine matrix and a rotation operation is discussed. Among the theorems and techniques presented are (1) various ways to characterize a rotation including proper orthogonal tensors, dyadics, Euler angles, axis/angle representations, series expansions, and quaternions; (2) the Euler-Rodrigues formula for converting axis and angle to a rotation tensor; (3) the distinction between rotations and reflections, along with implications for “handedness” of coordinate systems; (4) non-commu- tivity of sequential rotations, (5) eigenvalues and eigenvectors of a rotation; (6) the polar decomposition theorem for expressing a general deformation as a se- quence of shape and volume changes in combination with pure rotations; (7) mix- ing rotations in Eulerian hydrocodes or interpolating rotations in discrete field approximations; (8) Rates of rotation and the difference between spin and vortici- ty, (9) Random rotations for simulating crystal distributions; (10) The principle of material frame indifference (PMFI); and (11) a tensor-analysis presentation of classical rigid body mechanics, including direct notation expressions for momen- tum and energy and the extremely compact direct notation formulation of Euler’s equations (i.e., Newton’s law for rigid bodies).
    [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]