# Closed-Form Solution Of Absolute Orientation Using Unit Quaternions

Size: px
Start display at page:

Download "Closed-Form Solution Of Absolute Orientation Using Unit Quaternions"

Transcription

1 Closed-Form Solution Of Absolute Orientation Using Unit Berthold K. P. Horn Department of Computer and Information Sciences November 11, 2004

2 Outline 1 Introduction 2 3

3 The Problem Given: two sets of corresponding points in different coordinate systems Task: find the absolute orientation of the two systems (scale, rotation, translation) Previous approaches: Only iterative solutions Only close to least-squares solution Only handle three points Selectively neglect constraints This paper: Closed-form, exact least-squares solution No constraints neglected

4 Quaternion Notation q = q 0 + iq x + jq y + kq z q = (q 0, q x, q y, q z ) T : Basic Properties i 2 = 1 j 2 = 1 k 2 = 1 ij = k jk = i ki = j ji = k kj = i ik = j

5 Products of Quaternion Product as a Matrix Vector Product r 0 r x r y r z r q = r x r 0 r z r y r y r z r 0 r x q = R q q r = r z r y r x r 0 r 0 r x r y r z r x r 0 r z r y r y r z r 0 r x q = R q r z r y r x r 0 Watch Out! The quaternion product is not commutative!

6 Dot Product Dot Product of Two p q = p 0 q 0 + p x q x + p y q y + p z q z q 2 = q q Conjugate and Inverse p = q 0 iq x jq y kq z q q = q0 2 + qx 2 + qy 2 + qz 2 = q q A non-zero quaternion has an inverse q 1 = q q q

7 Useful Properties of Dot Products Representing Vectors Let r = (x, y, z) T R 3, then ( q p) ( q r) = ( q q)( p r) ( p q) r = p ( r q ) r = 0 + ix + iy + iz. Matrices associated with such purely imaginary quaternions are skew symmetric: (go back to matrices) R T = R R T = R

8 Representing s With Representation Using a Unit Quaternion q = cos φ 2 + sin φ 2 (iω x + jω y + kω z ) }{{} unit vector ω The imaginary part gives the direction of the axis of rotation. The angle is encoded into the real part and the magnitude of the imaginary part. Rotating a Vector Note that a vector is represented using an imaginary quaternion! r = q r q r = p r p = p( q r q ) p = p q r( p q)

9 Given: n points, each measured in a left and right coordinate system {r l,i } and {r r,i } Try to find a transformation of the form r r = sr(r l ) + r 0 from the left to the right coordinate system. There will be errors e i = r r,i sr(r l,i ) r 0 Minimize the sum of squares of errors e i 2

10 Introduction Working Relative to the Centroids r l = 1 r l,i n r l,i = r l,i r l r r = 1 n r r,i r r,i = r r,i r r Rewrite the error term: e i = r r,i sr(r l,i ) r 0 where r 0 = r 0 r r + sr( r l ) The sum of squares of errors becomes e i 2 = r r,i sr(r l,i ) r 0 2

11 Introduction Decompose the sum: = r r,i sr(r l,i ) r 0 2 ( ) r r,i sr(r l,i ) 2 2r 0 r r,i sr(r l,i ) +n r 0 2 }{{} independent of translation } {{ } =0 The total error is minimized with r 0 = 0, i.e. r 0 = r r sr( r l ). Remaining error term: e i = r r,i sr(r l,i). Now minimize r r,i sr(r l,i ) 2.

12 Finding the Decompose the sum: = r r,i sr(r l,i ) 2 r r,i 2 2s r r,i R(r l,i ) + s2 r l,i 2 = S r 2sD + s 2 S l Complete the square: Best scale is s = D/S l. (s S l D/ S l ) 2 + (S }{{} r S l D 2 )/S }{{} l =0 independent of scale

13 Finding the s = D/S l = r r,i R(r l,i ) r l,i 2 This scale factor is not symmetric! When going from the right to the left system, we get s 1 s Use a symmetrical expression for the error term instead: e i = 1 s r r,i sr(r l,i ) = e i s

14 Finding the Total error is then 1 r s r,i 2 2 r r,i R(r l,i )+s r l,i 2 = 1 s S r 2D ss l Complete the square (slightly different from Horn p. 632) ( s S l S r / s) 2 + 2( S }{{} l S r D) }{{} =0 independent of scale Best scale is now independent of rotation: s = S r /S l = n / r r,i 2 r l,i 2

15 Introduction Remaining error term to be minimized: (S r S l D 2 )/S l or 2( S l S r D) (asymmetric) (symmetric) I.e., maximize D = r r,i R(r l,i ) Use the quaternion representation: find the unit quaternion q that maximizes D = ( q r l,i q ) r r,i

16 Introduction Using ( p q) r = p ( r q ) from earlier, rewrite the term as ( q r l,i q ) r r,i = ( q r l,i ) ( r r,i q) = = ( R l,i q) (R r,i q) q T R T l,i R r,i q ) ( = q T R T l,i R r,i q = q T ( N i ) q = q T N q

17 Introduction Utilize the 3 3 matrix M := S xx S xy S xz r l,i r T r,i := S yx S yy S yz S zx S zy S zz for a convenient representation of ( Sxx + S yy + S zz S yz S zy S zx S xz S xy S yx ) N = S yz S zy S zx S xz S xx S yy S zz S xy + S yx S xy + S yx S xx + S yy S zz S zx + S xz S yz + S zy S xy S yx S zx + S xz S yz + S zy S xx S yy + S zz Note that N is symmetric, and trace(n) = 0. This contains all information to solve the least-squares problem for rotation!

18 Introduction The unit quaternion that maximizes q T N q is the eigenvector to the most positive eigenvalue of N To find this eigenvalue solve the quartic obtained from Use e.g. Ferrari s method. det(n λi) = 0 For the eigenvalue λ m, the eigenvector e m is found by solving (N λ m I) e m = 0 A lot easier nowadays using SVD...

19 The Algorithm 1 Find the centroids r l, r r, substract from measurements 2 For each pair ( r l = (x l, y l, z l ), r r = (x r, y r, z r) ) compute all possible products x l x r, x l y r,..., z l z r and add up to obtain S xx, S xy,..., S zz 3 Compute elements of N 4 Calculate the coefficients of the quartic and solve quartic 5 Pick the most positive root and obtain corresponding eigenvector. The quaternion representing the rotation is a unit vector in the same direction 6 Compute the scale 7 Compute the translation

20 Special Cases and Extensions If points are coplanar (e.g. only three points given), the calculation simplifies greatly Can also use weighted errors: Minimize weighted sum of errors: w i e i 2 Calculate weighted centroids: r l = n / w i r l,i w i etc. Change scale factor calculation: S r = n w i r r,i 2 etc. Change components of matrix for rotation recovery: S xx = n w ix l,i x r,i

21 MatLab Implementation

22 I presented a closed-form solution to the absolute orientation problem Given a mechanism for SVD or eigenvalue computation, the solution is straightforward Non-Gaussian noise, statistical outliers are not handled well Numerical stability? Project Combine with registration? Use for non-rigid motion? E.g., by doing absolute orientation on local point sets?

23 Appendix References References Berthold K. P. Horn. Closed-form solution of absolute orientation using unit quaternions. Journal of the Optical Society of America A, 4(4): , April Ross J. Micheals, Terrance E. Boult. On the Robustness of Absolute Orientation F. Sansò. An exact solution of the roto-translation problem. Photogrammetria 29: , 1973.

### A DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any

Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#\$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y \$ Z % Y Y x x } / % «] «] # z» & Y X»

More information

### Two Posts to Fill On School Board

Y Y 9 86 4 4 qz 86 x : ( ) z 7 854 Y x 4 z z x x 4 87 88 Y 5 x q x 8 Y 8 x x : 6 ; : 5 x ; 4 ( z ; ( ) ) x ; z 94 ; x 3 3 3 5 94 ; ; ; ; 3 x : 5 89 q ; ; x ; x ; ; x : ; ; ; ; ; ; 87 47% : () : / : 83

More information

### Course 2BA1: Hilary Term 2007 Section 8: Quaternions and Rotations

Course BA1: Hilary Term 007 Section 8: Quaternions and Rotations David R. Wilkins Copyright c David R. Wilkins 005 Contents 8 Quaternions and Rotations 1 8.1 Quaternions............................ 1 8.

More information

### A. H. Hall, 33, 35 &37, Lendoi

7 X x > - z Z - ----»»x - % x x» [> Q - ) < % - - 7»- -Q 9 Q # 5 - z -> Q x > z»- ~» - x " < z Q q»» > X»? Q ~ - - % % < - < - - 7 - x -X - -- 6 97 9

More information

### LOWELL WEEKLY JOURNAL

Y -» \$ 5 Y 7 Y Y -Y- Q x Q» 75»»/ q } # ]»\ - - \$ { Q» / X x»»- 3 q \$ 9 ) Y q - 5 5 3 3 3 7 Q q - - Q _»»/Q Y - 9 - - - )- [ X 7» -» - )»? / /? Q Y»» # X Q» - -?» Q ) Q \ Q - - - 3? 7» -? #»»» 7 - / Q

More information

### Rigid body simulation. Once we consider an object with spatial extent, particle system simulation is no longer sufficient

Rigid body dynamics Rigid body simulation Once we consider an object with spatial extent, particle system simulation is no longer sufficient Rigid body simulation Unconstrained system no contact Constrained

More information

### Lesson Rigid Body Dynamics

Lesson 8 Rigid Body Dynamics Lesson 8 Outline Problem definition and motivations Dynamics of rigid bodies The equation of unconstrained motion (ODE) User and time control Demos / tools / libs Rigid Body

More information

### ' Liberty and Umou Ono and Inseparablo "

3 5? #< q 8 2 / / ) 9 ) 2 ) > < _ / ] > ) 2 ) ) 5 > x > [ < > < ) > _ ] ]? <

More information

### Why Transforms? Want to animate objects and camera Translations Rotations Shears And more.. Want to be able to use projection transforms

Why Transforms? Want to animate objects and camera Translations Rotations Shears And more.. Want to be able to use projection transforms ITCS 3050:Game Engine Programming 1 Geometric Transformations Implementing

More information

### LOWELL WEEKLY JOURNAL

Y G y G Y 87 y Y 8 Y - \$ X ; ; y y q 8 y \$8 \$ \$ \$ G 8 q < 8 6 4 y 8 7 4 8 8 < < y 6 \$ q - - y G y G - Y y y 8 y y y Y Y 7-7- G - y y y ) y - y y y y - - y - y 87 7-7- G G < G y G y y 6 X y G y y y 87 G

More information

### Lecture 8 Analyzing the diffusion weighted signal. Room CSB 272 this week! Please install AFNI

Lecture 8 Analyzing the diffusion weighted signal Room CSB 272 this week! Please install AFNI http://afni.nimh.nih.gov/afni/ Next lecture, DTI For this lecture, think in terms of a single voxel We re still

More information

### Assignment 11 (C + C ) = (C + C ) = (C + C) i(c C ) ] = i(c C) (AB) = (AB) = B A = BA 0 = [A, B] = [A, B] = (AB BA) = (AB) AB

Arfken 3.4.6 Matrix C is not Hermition. But which is Hermitian. Likewise, Assignment 11 (C + C ) = (C + C ) = (C + C) [ i(c C ) ] = i(c C ) = i(c C) = i ( C C ) Arfken 3.4.9 The matrices A and B are both

More information

### Group, Rings, and Fields Rahul Pandharipande. I. Sets Let S be a set. The Cartesian product S S is the set of ordered pairs of elements of S,

Group, Rings, and Fields Rahul Pandharipande I. Sets Let S be a set. The Cartesian product S S is the set of ordered pairs of elements of S, A binary operation φ is a function, S S = {(x, y) x, y S}. φ

More information

### Homework 1/Solutions. Graded Exercises

MTH 310-3 Abstract Algebra I and Number Theory S18 Homework 1/Solutions Graded Exercises Exercise 1. Below are parts of the addition table and parts of the multiplication table of a ring. Complete both

More information

### OWELL WEEKLY JOURNAL

Y \»< - } Y Y Y & #»»» q ] q»»»>) & - - - } ) x ( - { Y» & ( x - (» & )< - Y X - & Q Q» 3 - x Q Y 6 \Y > Y Y X 3 3-9 33 x - - / - -»- --

More information

### County Council Named for Kent

\ Y Y 8 9 69 6» > 69 ««] 6 : 8 «V z 9 8 x 9 8 8 8?? 9 V q» :: q;; 8 x () «; 8 x ( z x 9 7 ; x >«\ 8 8 ; 7 z x [ q z «z : > ; ; ; ( 76 x ; x z «7 8 z ; 89 9 z > q _ x 9 : ; 6? ; ( 9 [ ) 89 _ ;»» «; x V

More information

### ON A NEW SPECIES OF IMAGINARY QUANTITIES CONNECTED WITH A THEORY OF QUATERNIONS. William Rowan Hamilton

ON A NEW SPECIES OF IMAGINARY QUANTITIES CONNECTED WITH A THEORY OF QUATERNIONS By William Rowan Hamilton (Proceedings of the Royal Irish Academy, 2 (1844), 424 434.) Edited by David R. Wilkins 1999 On

More information

### MECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso

MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Preliminary Math Concept of Stress Stress Components Equilibrium

More information

### CS 246 Review of Linear Algebra 01/17/19

1 Linear algebra In this section we will discuss vectors and matrices. We denote the (i, j)th entry of a matrix A as A ij, and the ith entry of a vector as v i. 1.1 Vectors and vector operations A vector

More information

### Useful Formulae ( )

Appendix A Useful Formulae (985-989-993-) 34 Jeremić et al. A.. CHAPTER SUMMARY AND HIGHLIGHTS page: 35 of 536 A. Chapter Summary and Highlights A. Stress and Strain This section reviews small deformation

More information

### GG303 Lecture 6 8/27/09 1 SCALARS, VECTORS, AND TENSORS

GG303 Lecture 6 8/27/09 1 SCALARS, VECTORS, AND TENSORS I Main Topics A Why deal with tensors? B Order of scalars, vectors, and tensors C Linear transformation of scalars and vectors (and tensors) II Why

More information

### oenofc : COXT&IBCTOEU. AU skaacst sftwer thsa4 aafcekr will be ehat«s«ai Bi. C. W. JUBSSOS. PERFECT THBOUGH SDFFEBISG. our

x V - --- < x x 35 V? 3?/ -V 3 - ) - - [ Z8 - & Z - - - - - x 0-35 - 3 75 3 33 09 33 5 \ - - 300 0 ( -? 9 { - - - -- - < - V 3 < < - - Z 7 - z 3 - [ } & _ 3 < 3 ( 5 7< ( % --- /? - / 4-4 - & - % 4 V 2

More information

### Neatest and Promptest Manner. E d i t u r ami rul)lihher. FOIt THE CIIILDIIES'. Trifles.

» ~ \$ ) 7 x X ) / ( 8 2 X 39 ««x» ««! «! / x? \» «({? «» q «(? (?? x! «? 8? ( z x x q? ) «q q q ) x z x 69 7( X X ( 3»«! ( ~«x ««x ) (» «8 4 X «4 «4 «8 X «x «(» X) ()»» «X «97 X X X 4 ( 86) x) ( ) z z

More information

### Lecture 7. Quaternions

Matthew T. Mason Mechanics of Manipulation Spring 2012 Today s outline Motivation Motivation have nice geometrical interpretation. have advantages in representing rotation. are cool. Even if you don t

More information

### LOWELL WEEKLY JOURNAL

G \$ G 2 G ««2 ««q ) q «\ { q «««/ 6 «««««q «] «q 6 ««Z q «««Q \ Q «q «X ««G X G ««? G Q / Q Q X ««/«X X «««Q X\ «q «X \ / X G XX «««X «x «X «x X G X 29 2 ««Q G G «) 22 G XXX GG G G G G G X «x G Q «) «G

More information

### FMIA. Fluid Mechanics and Its Applications 113 Series Editor: A. Thess. Moukalled Mangani Darwish. F. Moukalled L. Mangani M.

FMIA F. Moukalled L. Mangani M. Darwish An Advanced Introduction with OpenFOAM and Matlab This textbook explores both the theoretical foundation of the Finite Volume Method (FVM) and its applications in

More information

### Introduction to quaternions

. Introduction Introduction to uaternions Invented and developed by William Hamilton in 843, uaternions are essentially a generalization of complex numbers to four dimensions (one real dimension, three

More information

### L bor y nnd Union One nnd Inseparable. LOW I'LL, MICHIGAN. WLDNHSDA Y. JULY ), I8T. liuwkll NATIdiNAI, liank

G k y \$5 y / >/ k «««# ) /% < # «/» Y»««««?# «< >«>» y k»» «k F 5 8 Y Y F G k F >«y y

More information

### Vision 3D articielle Session 2: Essential and fundamental matrices, their computation, RANSAC algorithm

Vision 3D articielle Session 2: Essential and fundamental matrices, their computation, RANSAC algorithm Pascal Monasse monasse@imagine.enpc.fr IMAGINE, École des Ponts ParisTech Contents Some useful rules

More information

### P A L A C E P IE R, S T. L E O N A R D S. R a n n o w, q u a r r y. W WALTER CR O TC H, Esq., Local Chairman. E. CO O PER EVANS, Esq.,.

? ( # [ ( 8? [ > 3 Q [ ««> » 9 Q { «33 Q> 8 \ \ 3 3 3> Q»«9 Q ««« 3 8 3 8 X \ [ 3 ( ( Z ( Z 3( 9 9 > < < > >? 8 98 ««3 ( 98 < # # Q 3 98? 98 > > 3 8 9 9 ««««> 3 «>

More information

### A Memorial. Death Crash Branch Out. Symbol The. at Crossing Flaming Poppy. in Belding

- G Y Y 8 9 XXX G - Y - Q 5 8 G Y G Y - - * Y G G G G 9 - G - - : - G - - ) G G- Y G G q G G : Q G Y G 5) Y : z 6 86 ) ; - ) z; G ) 875 ; ) ; G -- ) ; Y; ) G 8 879 99 G 9 65 q 99 7 G : - G G Y ; - G 8

More information

382 4 7 q X

More information

### Problem Set 2 Due Tuesday, September 27, ; p : 0. (b) Construct a representation using five d orbitals that sit on the origin as a basis: 1

Problem Set 2 Due Tuesday, September 27, 211 Problems from Carter: Chapter 2: 2a-d,g,h,j 2.6, 2.9; Chapter 3: 1a-d,f,g 3.3, 3.6, 3.7 Additional problems: (1) Consider the D 4 point group and use a coordinate

More information

### THE I Establiifrad June, 1893

89 : 8 Y Y 2 96 6 - - : - 2 q - 26 6 - - q 2 2 2 4 6 4«4 ' V () 8 () 6 64-4 '2" () 6 ( ) * 'V ( 4 ) 94-4 q ( / ) K ( x- 6% j 9*V 2'%" 222 27 q - - K 79-29 - K x 2 2 j - -% K 4% 2% 6% ' K - 2 47 x - - j

More information

### Covariance to PCA. CS 510 Lecture #14 February 23, 2018

Covariance to PCA CS 510 Lecture 14 February 23, 2018 Overview: Goal Assume you have a gallery (database) of images, and a probe (test) image. The goal is to find the database image that is most similar

More information

### Group Theory: Math30038, Sheet 6

Group Theory: Math30038, Sheet 6 Solutions GCS 1. Consider the group D ofrigidsymmetriesofaregularn-gon (which may be turned over). Prove that this group has order 2n, is non-abelian, can be generated

More information

### Quantum Information & Quantum Computing

Math 478, Phys 478, CS4803, February 9, 006 1 Georgia Tech Math, Physics & Computing Math 478, Phys 478, CS4803 Quantum Information & Quantum Computing Problems Set 1 Due February 9, 006 Part I : 1. Read

More information

### GG612 Lecture 3. Outline

GG61 Lecture 3 Strain and Stress Should complete infinitesimal strain by adding rota>on. Outline Matrix Opera+ons Strain 1 General concepts Homogeneous strain 3 Matrix representa>ons 4 Squares of line

More information

### LOWELL JOURNAL. MUST APOLOGIZE. such communication with the shore as Is m i Boimhle, noewwary and proper for the comfort

- 7 7 Z 8 q ) V x - X > q - < Y Y X V - z - - - - V - V - q \ - q q < -- V - - - x - - V q > x - x q - x q - x - - - 7 -» - - - - 6 q x - > - - x - - - x- - - q q - V - x - - ( Y q Y7 - >»> - x Y - ] [

More information

### Stress, Strain, Mohr s Circle

Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected

More information

### Applications of Eigenvalues & Eigenvectors

Applications of Eigenvalues & Eigenvectors Louie L. Yaw Walla Walla University Engineering Department For Linear Algebra Class November 17, 214 Outline 1 The eigenvalue/eigenvector problem 2 Principal

More information

### CS286.2 Lecture 8: A variant of QPCP for multiplayer entangled games

CS286.2 Lecture 8: A variant of QPCP for multiplayer entangled games Scribe: Zeyu Guo In the first lecture, we saw three equivalent variants of the classical PCP theorems in terms of CSP, proof checking,

More information

### 6. 3D Kinematics DE2-EA 2.1: M4DE. Dr Connor Myant

DE2-EA 2.1: M4DE Dr Connor Myant 6. 3D Kinematics Comments and corrections to connor.myant@imperial.ac.uk Lecture resources may be found on Blackboard and at http://connormyant.com Contents Three-Dimensional

More information

### forms Christopher Engström November 14, 2014 MAA704: Matrix factorization and canonical forms Matrix properties Matrix factorization Canonical forms

Christopher Engström November 14, 2014 Hermitian LU QR echelon Contents of todays lecture Some interesting / useful / important of matrices Hermitian LU QR echelon Rewriting a as a product of several matrices.

More information

### M E 320 Professor John M. Cimbala Lecture 10

M E 320 Professor John M. Cimbala Lecture 10 Today, we will: Finish our example problem rates of motion and deformation of fluid particles Discuss the Reynolds Transport Theorem (RTT) Show how the RTT

More information

### Demonstration of the Coupled Evolution Rules 163 APPENDIX F: DEMONSTRATION OF THE COUPLED EVOLUTION RULES

Demonstration of the Coupled Evolution Rules 163 APPENDIX F: DEMONSTRATION OF THE COUPLED EVOLUTION RULES Before going into the demonstration we need to point out two limitations: a. It assumes I=1/2 for

More information

### A b r i l l i a n t young chemist, T h u r e Wagelius of N e w Y o r k, ac. himself with eth

6 6 0 x J 8 0 J 0 z (0 8 z x x J x 6 000 X j x "" "" " " x " " " x " " " J " " " " " " " " x : 0 z j ; J K 0 J K q 8 K K J x 0 j " " > J x J j z ; j J q J 0 0 8 K J 60 : K 6 x 8 K J :? 0 J J K 0 6% 8 0

More information

### Complex Numbers and Quaternions for Calc III

Complex Numbers and Quaternions for Calc III Taylor Dupuy September, 009 Contents 1 Introduction 1 Two Ways of Looking at Complex Numbers 1 3 Geometry of Complex Numbers 4 Quaternions 5 4.1 Connection

More information

### Knowledge Discovery and Data Mining 1 (VO) ( )

Knowledge Discovery and Data Mining 1 (VO) (707.003) Review of Linear Algebra Denis Helic KTI, TU Graz Oct 9, 2014 Denis Helic (KTI, TU Graz) KDDM1 Oct 9, 2014 1 / 74 Big picture: KDDM Probability Theory

More information

### Tensor Visualization. CSC 7443: Scientific Information Visualization

Tensor Visualization Tensor data A tensor is a multivariate quantity Scalar is a tensor of rank zero s = s(x,y,z) Vector is a tensor of rank one v = (v x,v y,v z ) For a symmetric tensor of rank 2, its

More information

### Linear Algebra. Chapter 8: Eigenvalues: Further Applications and Computations Section 8.2. Applications to Geometry Proofs of Theorems.

Linear Algebra Chapter 8: Eigenvalues: Further Applications and Computations Section 8.2. Applications to Geometry Proofs of Theorems May 1, 2018 () Linear Algebra May 1, 2018 1 / 8 Table of contents 1

More information

«4 [< «

More information

### CS 143 Linear Algebra Review

CS 143 Linear Algebra Review Stefan Roth September 29, 2003 Introductory Remarks This review does not aim at mathematical rigor very much, but instead at ease of understanding and conciseness. Please see

More information

### ECS 178 Course Notes QUATERNIONS

ECS 178 Course Notes QUATERNIONS Kenneth I. Joy Institute for Data Analysis and Visualization Department of Computer Science University of California, Davis Overview The quaternion number system was discovered

More information

### Problem Set 2 Due Thursday, October 1, & & & & # % (b) Construct a representation using five d orbitals that sit on the origin as a basis:

Problem Set 2 Due Thursday, October 1, 29 Problems from Cotton: Chapter 4: 4.6, 4.7; Chapter 6: 6.2, 6.4, 6.5 Additional problems: (1) Consider the D 3h point group and use a coordinate system wherein

More information

### LOWELL WEEKLY JOURNAL.

Y 5 ; ) : Y 3 7 22 2 F \$ 7 2 F Q 3 q q 6 2 3 6 2 5 25 2 2 3 \$2 25: 75 5 \$6 Y q 7 Y Y # \ x Y : { Y Y Y : ( \ _ Y ( ( Y F [ F F ; x Y : ( : G ( ; ( ~ x F G Y ; \ Q ) ( F \ Q / F F \ Y () ( \ G Y ( ) \F

More information

### V o l u m e 5, N u m b e r 5 2, 1 6 P a g e s. Gold B e U ClUt Stamps Double Stamp D a y E v e r y Wednesday

1 6 5 J 9 6 " " z k ; k x k k k z z k j " " ( k " " k 8 1959 " " x k j 5 25 ; ; k k qz ; x 13 x k * k ( ) k k : qz 13 k k k j ; q k x ; x 615 26 ( : k z 113 99751 z k k q ; 15 k k k j q " " k j x x ( *»

More information

### The Product Operator Formalism

2 The Product Operator Formalism 1. INTRODUCTION In this section we will see that the density matrix at equilibrium can be expressed in terms of the spin angular momentum component I z of each nucleus.

More information

### MANY BILLS OF CONCERN TO PUBLIC

- 6 8 9-6 8 9 6 9 XXX 4 > -? - 8 9 x 4 z ) - -! x - x - - X - - - - - x 00 - - - - - x z - - - x x - x - - - - - ) x - - - - - - 0 > - 000-90 - - 4 0 x 00 - -? z 8 & x - - 8? > 9 - - - - 64 49 9 x - -

More information

### Lagrange Multipliers

Optimization with Constraints As long as algebra and geometry have been separated, their progress have been slow and their uses limited; but when these two sciences have been united, they have lent each

More information

### Rotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep.

Rotational Motion Chapter 4 P. J. Grandinetti Chem. 4300 Sep. 1, 2017 P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. 1, 2017 1 / 76 Angular Momentum The angular momentum of a particle with respect

More information

### Assignment 10. Arfken Show that Stirling s formula is an asymptotic expansion. The remainder term is. B 2n 2n(2n 1) x1 2n.

Assignment Arfken 5.. Show that Stirling s formula is an asymptotic expansion. The remainder term is R N (x nn+ for some N. The condition for an asymptotic series, lim x xn R N lim x nn+ B n n(n x n B

More information

### Modeling photon generation

Modeling photon generation C. J. McKinstrie Bell Laboratories, Alcatel Lucent, Holmdel, New Jersey 07733 mckinstrie@alcatel-lucent.com / 732-888-7275 Abstract The photons required for a variety of quantum

More information

### ..«W- tn^zmxmmrrx/- NEW STORE. Popular Goods at Popular D. E. SPRING, Mas just opened a large fdo.k of DRY GOODS & GROCERIES,

B y «X }() z zxx/ X y y y y )3 y «y

More information

### Classical Mechanics. Luis Anchordoqui

1 Rigid Body Motion Inertia Tensor Rotational Kinetic Energy Principal Axes of Rotation Steiner s Theorem Euler s Equations for a Rigid Body Eulerian Angles Review of Fundamental Equations 2 Rigid body

More information

### CSE 554 Lecture 7: Alignment

CSE 554 Lecture 7: Alignment Fall 2012 CSE554 Alignment Slide 1 Review Fairing (smoothing) Relocating vertices to achieve a smoother appearance Method: centroid averaging Simplification Reducing vertex

More information

### Lecture 4: Least Squares (LS) Estimation

ME 233, UC Berkeley, Spring 2014 Xu Chen Lecture 4: Least Squares (LS) Estimation Background and general solution Solution in the Gaussian case Properties Example Big picture general least squares estimation:

More information

### Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik. Robot Dynamics. Dr.-Ing. John Nassour J.

Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik Robot Dynamics Dr.-Ing. John Nassour 25.1.218 J.Nassour 1 Introduction Dynamics concerns the motion of bodies Includes Kinematics

More information

### Lecture 5: Random Walks and Markov Chain

Spectral Graph Theory and Applications WS 20/202 Lecture 5: Random Walks and Markov Chain Lecturer: Thomas Sauerwald & He Sun Introduction to Markov Chains Definition 5.. A sequence of random variables

More information

### 19. Principal Stresses

19. Principal Stresses I Main Topics A Cauchy s formula B Principal stresses (eigenvectors and eigenvalues) C Example 10/24/18 GG303 1 19. Principal Stresses hkp://hvo.wr.usgs.gov/kilauea/update/images.html

More information

### Dynamics. 1 Copyright c 2015 Roderic Grupen

Dynamics The branch of physics that treats the action of force on bodies in motion or at rest; kinetics, kinematics, and statics, collectively. Websters dictionary Outline Conservation of Momentum Inertia

More information

### On complexified quantum mechanics and space-time

On complexified quantum mechanics and space-time Dorje C. Brody Mathematical Sciences Brunel University, Uxbridge UB8 3PH dorje.brody@brunel.ac.uk Quantum Physics with Non-Hermitian Operators Dresden:

More information

### APPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.

APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product

More information

### Elementary Linear Algebra

Matrices J MUSCAT Elementary Linear Algebra Matrices Definition Dr J Muscat 2002 A matrix is a rectangular array of numbers, arranged in rows and columns a a 2 a 3 a n a 2 a 22 a 23 a 2n A = a m a mn We

More information

### ELA

Volume 5 pp 97-33 November 6 http://mathtechnionacil/iic/ela ZEROS OF UNILATERAL QUATERNIONIC POLYNOMIALS STEFANO DE LEO GISELE DUCATI AND VINICIUS LEONARDI Abstract The purpose of this paper is to show

More information

### Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition

Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2016 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing

More information

### Bidiagonal pairs, Tridiagonal pairs, Lie algebras, and Quantum Groups

Bidiagonal pairs, Tridiagonal pairs, Lie algebras, and Quantum Groups Darren Funk-Neubauer Department of Mathematics and Physics Colorado State University - Pueblo Pueblo, Colorado, USA darren.funkneubauer@colostate-pueblo.edu

More information

### Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Rigid body physics Particle system Most simple instance of a physics system Each object (body) is a particle Each particle

More information

### Quaternions. R. J. Renka 11/09/2015. Department of Computer Science & Engineering University of North Texas. R. J.

Quaternions R. J. Renka Department of Computer Science & Engineering University of North Texas 11/09/2015 Definition A quaternion is an element of R 4 with three operations: addition, scalar multiplication,

More information

### Next is material on matrix rank. Please see the handout

B90.330 / C.005 NOTES for Wednesday 0.APR.7 Suppose that the model is β + ε, but ε does not have the desired variance matrix. Say that ε is normal, but Var(ε) σ W. The form of W is W w 0 0 0 0 0 0 w 0

More information

### General Physics I. Lecture 10: Rolling Motion and Angular Momentum.

General Physics I Lecture 10: Rolling Motion and Angular Momentum Prof. WAN, Xin (万歆) 万歆 ) xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ Outline Rolling motion of a rigid object: center-of-mass motion

More information

### Properties of Matrices and Operations on Matrices

Properties of Matrices and Operations on Matrices A common data structure for statistical analysis is a rectangular array or matris. Rows represent individual observational units, or just observations,

More information

### Rigid body dynamics. Basilio Bona. DAUIN - Politecnico di Torino. October 2013

Rigid body dynamics Basilio Bona DAUIN - Politecnico di Torino October 2013 Basilio Bona (DAUIN - Politecnico di Torino) Rigid body dynamics October 2013 1 / 16 Multiple point-mass bodies Each mass is

More information

### COMP 175 COMPUTER GRAPHICS. Lecture 04: Transform 1. COMP 175: Computer Graphics February 9, Erik Anderson 04 Transform 1

Lecture 04: Transform COMP 75: Computer Graphics February 9, 206 /59 Admin Sign up via email/piazza for your in-person grading Anderson@cs.tufts.edu 2/59 Geometric Transform Apply transforms to a hierarchy

More information

### Basic Equations of Elasticity

A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ

More information

### A Probability Review

A Probability Review Outline: A probability review Shorthand notation: RV stands for random variable EE 527, Detection and Estimation Theory, # 0b 1 A Probability Review Reading: Go over handouts 2 5 in

More information

### Determination of Locally Varying Directions through Mass Moment of Inertia Tensor

Determination of Locally Varying Directions through Mass Moment of Inertia Tensor R. M. Hassanpour and C.V. Deutsch Centre for Computational Geostatistics Department of Civil and Environmental Engineering

More information

### Macroscopic theory Rock as 'elastic continuum'

Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave

More information

### Math 5630: Iterative Methods for Systems of Equations Hung Phan, UMass Lowell March 22, 2018

1 Linear Systems Math 5630: Iterative Methods for Systems of Equations Hung Phan, UMass Lowell March, 018 Consider the system 4x y + z = 7 4x 8y + z = 1 x + y + 5z = 15. We then obtain x = 1 4 (7 + y z)

More information

### (d) Since we can think of isometries of a regular 2n-gon as invertible linear operators on R 2, we get a 2-dimensional representation of G for

Solutions to Homework #7 0. Prove that [S n, S n ] = A n for every n 2 (where A n is the alternating group). Solution: Since [f, g] = f 1 g 1 fg is an even permutation for all f, g S n and since A n is

More information

### MATH 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 information

### PEAT SEISMOLOGY Lecture 2: Continuum mechanics

PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a

More information

### An eightfold path to E 8

An eightfold path to E 8 Robert A. Wilson First draft 17th November 2008; this version 29th April 2012 Introduction Finite-dimensional real reflection groups were classified by Coxeter [2]. In two dimensions,

More information

### Solving Linear Systems of Equations

November 6, 2013 Introduction The type of problems that we have to solve are: Solve the system: A x = B, where a 11 a 1N a 12 a 2N A =.. a 1N a NN x = x 1 x 2. x N B = b 1 b 2. b N To find A 1 (inverse

More information

### Course MA2C02, Hilary Term 2010 Section 4: Vectors and Quaternions

Course MA2C02, Hilary Term 2010 Section 4: Vectors and Quaternions David R. Wilkins Copyright c David R. Wilkins 2000 2010 Contents 4 Vectors and Quaternions 47 4.1 Vectors...............................

More information

### Lowell Dam Gone Out. Streets Turned I n t o Rivers. No Cause For Alarm Now However As This Happened 35 Years A&o

V ()\\ ))? K K Y 6 96 Y - Y Y V 5 Z ( z x z \ - \ - - z - q q x x - x 5 9 Q \ V - - Y x 59 7 x x - Y - x - - x z - z x - ( 7 x V 9 z q &? - 9 - V ( x - - - V- [ Z x z - -x > -) - - > X Z z ( V V V

More information

### Linear Algebra Practice Final

. Let (a) First, Linear Algebra Practice Final Summer 3 3 A = 5 3 3 rref([a ) = 5 so if we let x 5 = t, then x 4 = t, x 3 =, x = t, and x = t, so that t t x = t = t t whence ker A = span(,,,, ) and a basis

More information

### LINEAR 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 information

### COMS 4721: Machine Learning for Data Science Lecture 19, 4/6/2017

COMS 4721: Machine Learning for Data Science Lecture 19, 4/6/2017 Prof. John Paisley Department of Electrical Engineering & Data Science Institute Columbia University PRINCIPAL COMPONENT ANALYSIS DIMENSIONALITY

More information

### 1 First Theme: Sums of Squares

I will try to organize the work of this semester around several classical questions. The first is, When is a prime p the sum of two squares? The question was raised by Fermat who gave the correct answer

More information