DOCSLIB.ORG

MATH 532: Linear Algebra Chapter 5: Norms, Inner Products and Orthogonality

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
MATH 532: Linear Algebra Chapter 5: Norms, Inner Products and Orthogonality MATH 532: Linear Algebra Chapter 5: Norms, Inner Products and Orthogonality Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Spring 2015 [email protected] MATH 532 1 Outline 1 Vector Norms 2 Matrix Norms 3 Inner Product Spaces 4 Orthogonal Vectors 5 Gram–Schmidt Orthogonalization & QR Factorization 6 Unitary and Orthogonal Matrices 7 Orthogonal Reduction 8 Complementary Subspaces 9 Orthogonal Decomposition 10 Singular Value Decomposition 11 Orthogonal Projections [email protected] MATH 532 2 Vector Norms Outline 1 Vector Norms 2 Matrix Norms 3 Inner Product Spaces 4 Orthogonal Vectors 5 Gram–Schmidt Orthogonalization & QR Factorization 6 Unitary and Orthogonal Matrices 7 Orthogonal Reduction 8 Complementary Subspaces 9 Orthogonal Decomposition 10 Singular Value Decomposition 11 Orthogonal Projections [email protected] MATH 532 3 Vector Norms Vector Norms Definition Let x; y 2 Rn (Cn). Then n T X x y = xi yi 2 R i=1 n ∗ X x y = x¯i yi 2 C i=1 is called the standard inner product for Rn (Cn). [email protected] MATH 532 4 Vector Norms Definition Let V be a vector space. A function k · k : V! R≥0 is called a norm provided for any x; y 2 V and α 2 R 1 kxk ≥ 0 and kxk = 0 if and only if x = 0, 2 kαxk = jαj kxk, 3 kx + yk ≤ kxk + kyk. Remark The inequality in (3) is known as the triangle inequality. [email protected] MATH 532 5 We will show that the standard inner product induces the Euclidean norm (cf. length of a vector). Remark Inner products let us define angles via xT y cos θ = : kxkkyk In particular, x; y are orthogonal if and only if xT y = 0. Vector Norms Remark Any inner product h·; ·i induces a norm via (more later) p kxk = hx; xi: [email protected] MATH 532 6 Remark Inner products let us define angles via xT y cos θ = : kxkkyk In particular, x; y are orthogonal if and only if xT y = 0. Vector Norms Remark Any inner product h·; ·i induces a norm via (more later) p kxk = hx; xi: We will show that the standard inner product induces the Euclidean norm (cf. length of a vector). [email protected] MATH 532 6 In particular, x; y are orthogonal if and only if xT y = 0. Vector Norms Remark Any inner product h·; ·i induces a norm via (more later) p kxk = hx; xi: We will show that the standard inner product induces the Euclidean norm (cf. length of a vector). Remark Inner products let us define angles via xT y cos θ = : kxkkyk [email protected] MATH 532 6 Vector Norms Remark Any inner product h·; ·i induces a norm via (more later) p kxk = hx; xi: We will show that the standard inner product induces the Euclidean norm (cf. length of a vector). Remark Inner products let us define angles via xT y cos θ = : kxkkyk In particular, x; y are orthogonal if and only if xT y = 0. [email protected] MATH 532 6 We show that k · k2 is a norm. We do this for the real case, but the complex case goes analogously. 1 Clearly, kxk2 ≥ 0. Also, 2 kxk2 = 0 () kxk2 = 0 n X 2 () xi = 0 () xi = 0; i = 1;:::; n; i=1 () x = 0: Vector Norms Example Let x 2 Rn and consider the Euclidean norm p T kxk2 = x x n !1=2 X 2 = xi : i=1 [email protected] MATH 532 7 1 Clearly, kxk2 ≥ 0. Also, 2 kxk2 = 0 () kxk2 = 0 n X 2 () xi = 0 () xi = 0; i = 1;:::; n; i=1 () x = 0: Vector Norms Example Let x 2 Rn and consider the Euclidean norm p T kxk2 = x x n !1=2 X 2 = xi : i=1 We show that k · k2 is a norm. We do this for the real case, but the complex case goes analogously. [email protected] MATH 532 7 Vector Norms Example Let x 2 Rn and consider the Euclidean norm p T kxk2 = x x n !1=2 X 2 = xi : i=1 We show that k · k2 is a norm. We do this for the real case, but the complex case goes analogously. 1 Clearly, kxk2 ≥ 0. Also, 2 kxk2 = 0 () kxk2 = 0 n X 2 () xi = 0 () xi = 0; i = 1;:::; n; i=1 () x = 0: [email protected] MATH 532 7 n !1=2 X 2 = jαj xi = jαj kxk2: i=1 3 To establish (3) we need Lemma Let x; y 2 Rn. Then T jx yj ≤ kxk2kyk2: (Cauchy–Schwarz–Bunyakovsky) Moreover, equality holds if and only if y = αx with xT y α = : 2 kxk2 Vector Norms Example (cont.) 2 We have n !1=2 X 2 kαxk2 = (αxi ) i=1 [email protected] MATH 532 8 3 To establish (3) we need Lemma Let x; y 2 Rn. Then T jx yj ≤ kxk2kyk2: (Cauchy–Schwarz–Bunyakovsky) Moreover, equality holds if and only if y = αx with xT y α = : 2 kxk2 Vector Norms Example (cont.) 2 We have n !1=2 n !1=2 X 2 X 2 kαxk2 = (αxi ) = jαj xi = jαj kxk2: i=1 i=1 [email protected] MATH 532 8 Vector Norms Example (cont.) 2 We have n !1=2 n !1=2 X 2 X 2 kαxk2 = (αxi ) = jαj xi = jαj kxk2: i=1 i=1 3 To establish (3) we need Lemma Let x; y 2 Rn. Then T jx yj ≤ kxk2kyk2: (Cauchy–Schwarz–Bunyakovsky) Moreover, equality holds if and only if y = αx with xT y α = : 2 kxk2 [email protected] MATH 532 8 aT b cos θ = : kakkbk Using trigonometry as in the figure, the projection of a onto b is then b aT b b aT b kak cos θ = kak = b: kbk kakkbk kbk kbk2 Now, we let y = a and x = b, so that the projection of y onto x is given by xT y αx; where α = : kxk2 Vector Norms Motivation for Proof of Cauchy–Schwarz–Bunyakovsky As already alluded to above, the angle θ between two vectors a and b is related to the inner product by [email protected] MATH 532 9 Using trigonometry as in the figure, the projection of a onto b is then b aT b b aT b kak cos θ = kak = b: kbk kakkbk kbk kbk2 Now, we let y = a and x = b, so that the projection of y onto x is given by xT y αx; where α = : kxk2 Vector Norms Motivation for Proof of Cauchy–Schwarz–Bunyakovsky As already alluded to above, the angle θ between two vectors a and b is related to the inner product by aT b cos θ = : kakkbk [email protected] MATH 532 9 the projection of a onto b is then b aT b b aT b kak cos θ = kak = b: kbk kakkbk kbk kbk2 Now, we let y = a and x = b, so that the projection of y onto x is given by xT y αx; where α = : kxk2 Vector Norms Motivation for Proof of Cauchy–Schwarz–Bunyakovsky As already alluded to above, the angle θ between two vectors a and b is related to the inner product by aT b cos θ = : kakkbk Using trigonometry as in the figure, [email protected] MATH 532 9 aT b b aT b = kak = b: kakkbk kbk kbk2 Now, we let y = a and x = b, so that the projection of y onto x is given by xT y αx; where α = : kxk2 Vector Norms Motivation for Proof of Cauchy–Schwarz–Bunyakovsky As already alluded to above, the angle θ between two vectors a and b is related to the inner product by aT b cos θ = : kakkbk Using trigonometry as in the figure, the projection of a onto b is then b kak cos θ kbk [email protected] MATH 532 9 aT b = b: kbk2 Now, we let y = a and x = b, so that the projection of y onto x is given by xT y αx; where α = : kxk2 Vector Norms Motivation for Proof of Cauchy–Schwarz–Bunyakovsky As already alluded to above, the angle θ between two vectors a and b is related to the inner product by aT b cos θ = : kakkbk Using trigonometry as in the figure, the projection of a onto b is then b aT b b kak cos θ = kak kbk kakkbk kbk [email protected] MATH 532 9 Now, we let y = a and x = b, so that the projection of y onto x is given by xT y αx; where α = : kxk2 Vector Norms Motivation for Proof of Cauchy–Schwarz–Bunyakovsky As already alluded to above, the angle θ between two vectors a and b is related to the inner product by aT b cos θ = : kakkbk Using trigonometry as in the figure, the projection of a onto b is then b aT b b aT b kak cos θ = kak = b: kbk kakkbk kbk kbk2 [email protected] MATH 532 9 Vector Norms Motivation for Proof of Cauchy–Schwarz–Bunyakovsky As already alluded to above, the angle θ between two vectors a and b is related to the inner product by aT b cos θ = : kakkbk Using trigonometry as in the figure, the projection of a onto b is then b aT b b aT b kak cos θ = kak = b: kbk kakkbk kbk kbk2 Now, we let y = a and x = b, so that the projection of y onto x is given by xT y αx; where α = : kxk2 [email protected] MATH 532 9 = y T y − 2αxT y + α2xT x 2 xT y xT y = y T y − 2 xT y + xT x kxk2 kxk4 |{z} 2 =kxk2 2 xT y = k k2 − : y 2 2 kxk2 This implies 2 T 2 2 x y ≤ kxk2kyk2; and the Cauchy–Schwarz–Bunyakovsky inequality follows by taking square roots.
Load more

Recommended publications
  • Introduction to Linear Bialgebra
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by University of New Mexico University of New Mexico UNM Digital Repository Mathematics and Statistics Faculty and Staff Publications Academic Department Resources 2005 INTRODUCTION TO LINEAR BIALGEBRA Florentin Smarandache University of New Mexico, [email protected] W.B. Vasantha Kandasamy K. Ilanthenral Follow this and additional works at: https://digitalrepository.unm.edu/math_fsp Part of the Algebra Commons, Analysis Commons, Discrete Mathematics and Combinatorics Commons, and the Other Mathematics Commons Recommended Citation Smarandache, Florentin; W.B. Vasantha Kandasamy; and K. Ilanthenral. "INTRODUCTION TO LINEAR BIALGEBRA." (2005). https://digitalrepository.unm.edu/math_fsp/232 This Book is brought to you for free and open access by the Academic Department Resources at UNM Digital Repository. It has been accepted for inclusion in Mathematics and Statistics Faculty and Staff Publications by an authorized administrator of UNM Digital Repository. For more information, please contact [email protected], [email protected], [email protected]. INTRODUCTION TO LINEAR BIALGEBRA W. B. Vasantha Kandasamy Department of Mathematics Indian Institute of Technology, Madras Chennai – 600036, India e-mail: [email protected] web: http://mat.iitm.ac.in/~wbv Florentin Smarandache Department of Mathematics University of New Mexico Gallup, NM 87301, USA e-mail: [email protected] K. Ilanthenral Editor, Maths Tiger, Quarterly Journal Flat No.11, Mayura Park, 16, Kazhikundram Main Road, Tharamani, Chennai – 600 113, India e-mail: [email protected] HEXIS Phoenix, Arizona 2005 1 This book can be ordered in a paper bound reprint from: Books on Demand ProQuest Information & Learning (University of Microfilm International) 300 N.
    [Show full text]
  • A New Description of Space and Time Using Clifford Multivectors
    A new description of space and time using Clifford multivectors James M. Chappell† , Nicolangelo Iannella† , Azhar Iqbal† , Mark Chappell‡ , Derek Abbott† †School of Electrical and Electronic Engineering, University of Adelaide, South Australia 5005, Australia ‡Griffith Institute, Griffith University, Queensland 4122, Australia Abstract Following the development of the special theory of relativity in 1905, Minkowski pro- posed a unified space and time structure consisting of three space dimensions and one time dimension, with relativistic effects then being natural consequences of this space- time geometry. In this paper, we illustrate how Clifford’s geometric algebra that utilizes multivectors to represent spacetime, provides an elegant mathematical framework for the study of relativistic phenomena. We show, with several examples, how the application of geometric algebra leads to the correct relativistic description of the physical phenomena being considered. This approach not only provides a compact mathematical representa- tion to tackle such phenomena, but also suggests some novel insights into the nature of time. Keywords: Geometric algebra, Clifford space, Spacetime, Multivectors, Algebraic framework 1. Introduction The physical world, based on early investigations, was deemed to possess three inde- pendent freedoms of translation, referred to as the three dimensions of space. This naive conclusion is also supported by more sophisticated analysis such as the existence of only five regular polyhedra and the inverse square force laws. If we lived in a world with four spatial dimensions, for example, we would be able to construct six regular solids, and in arXiv:1205.5195v2 [math-ph] 11 Oct 2012 five dimensions and above we would find only three [1].
    [Show full text]
  • Orthogonal Projection, Low Rank Approximation, and Orthogonal Bases
    Week 11 Orthogonal Projection, Low Rank Approximation, and Orthogonal Bases 11.1 Opening Remarks 11.1.1 Low Rank Approximation * View at edX 463 Week 11. Orthogonal Projection, Low Rank Approximation, and Orthogonal Bases 464 11.1.2 Outline 11.1. Opening Remarks..................................... 463 11.1.1. Low Rank Approximation............................. 463 11.1.2. Outline....................................... 464 11.1.3. What You Will Learn................................ 465 11.2. Projecting a Vector onto a Subspace........................... 466 11.2.1. Component in the Direction of ............................ 466 11.2.2. An Application: Rank-1 Approximation...................... 470 11.2.3. Projection onto a Subspace............................. 474 11.2.4. An Application: Rank-2 Approximation...................... 476 11.2.5. An Application: Rank-k Approximation...................... 478 11.3. Orthonormal Bases.................................... 481 11.3.1. The Unit Basis Vectors, Again........................... 481 11.3.2. Orthonormal Vectors................................ 482 11.3.3. Orthogonal Bases.................................. 485 11.3.4. Orthogonal Bases (Alternative Explanation).................... 488 11.3.5. The QR Factorization................................ 492 11.3.6. Solving the Linear Least-Squares Problem via QR Factorization......... 493 11.3.7. The QR Factorization (Again)........................... 494 11.4. Change of Basis...................................... 498 11.4.1. The Unit Basis Vectors,
    [Show full text]
  • Lecture 4: April 8, 2021 1 Orthogonality and Orthonormality
    Mathematical Toolkit Spring 2021 Lecture 4: April 8, 2021 Lecturer: Avrim Blum (notes based on notes from Madhur Tulsiani) 1 Orthogonality and orthonormality Definition 1.1 Two vectors u, v in an inner product space are said to be orthogonal if hu, vi = 0. A set of vectors S ⊆ V is said to consist of mutually orthogonal vectors if hu, vi = 0 for all u 6= v, u, v 2 S. A set of S ⊆ V is said to be orthonormal if hu, vi = 0 for all u 6= v, u, v 2 S and kuk = 1 for all u 2 S. Proposition 1.2 A set S ⊆ V n f0V g consisting of mutually orthogonal vectors is linearly inde- pendent. Proposition 1.3 (Gram-Schmidt orthogonalization) Given a finite set fv1,..., vng of linearly independent vectors, there exists a set of orthonormal vectors fw1,..., wng such that Span (fw1,..., wng) = Span (fv1,..., vng) . Proof: By induction. The case with one vector is trivial. Given the statement for k vectors and orthonormal fw1,..., wkg such that Span (fw1,..., wkg) = Span (fv1,..., vkg) , define k u + u = v − hw , v i · w and w = k 1 . k+1 k+1 ∑ i k+1 i k+1 k k i=1 uk+1 We can now check that the set fw1,..., wk+1g satisfies the required conditions. Unit length is clear, so let’s check orthogonality: k uk+1, wj = vk+1, wj − ∑ hwi, vk+1i · wi, wj = vk+1, wj − wj, vk+1 = 0. i=1 Corollary 1.4 Every finite dimensional inner product space has an orthonormal basis.
    [Show full text]
  • Orthogonality Handout
    3.8 (SUPPLEMENT) | ORTHOGONALITY OF EIGENFUNCTIONS We now develop some properties of eigenfunctions, to be used in Chapter 9 for Fourier Series and Partial Differential Equations. 1. Definition of Orthogonality R b We say functions f(x) and g(x) are orthogonal on a < x < b if a f(x)g(x) dx = 0 . [Motivation: Let's approximate the integral with a Riemann sum, as follows. Take a large integer N, put h = (b − a)=N and partition the interval a < x < b by defining x1 = a + h; x2 = a + 2h; : : : ; xN = a + Nh = b. Then Z b f(x)g(x) dx ≈ f(x1)g(x1)h + ··· + f(xN )g(xN )h a = (uN · vN )h where uN = (f(x1); : : : ; f(xN )) and vN = (g(x1); : : : ; g(xN )) are vectors containing the values of f and g. The vectors uN and vN are said to be orthogonal (or perpendicular) if their dot product equals zero (uN ·vN = 0), and so when we let N ! 1 in the above formula it makes R b sense to say the functions f and g are orthogonal when the integral a f(x)g(x) dx equals zero.] R π 1 2 π Example. sin x and cos x are orthogonal on −π < x < π, since −π sin x cos x dx = 2 sin x −π = 0. 2. Integration Lemma Suppose functions Xn(x) and Xm(x) satisfy the differential equations 00 Xn + λnXn = 0; a < x < b; 00 Xm + λmXm = 0; a < x < b; for some numbers λn; λm. Then Z b 0 0 b (λn − λm) Xn(x)Xm(x) dx = [Xn(x)Xm(x) − Xn(x)Xm(x)]a: a Proof.
    [Show full text]
  • Inner Products and Orthogonality
    Advanced Linear Algebra – Week 5 Inner Products and Orthogonality This week we will learn about: • Inner products (and the dot product again), • The norm induced by the inner product, • The Cauchy–Schwarz and triangle inequalities, and • Orthogonality. Extra reading and watching: • Sections 1.3.4 and 1.4.1 in the textbook • Lecture videos 17, 18, 19, 20, 21, and 22 on YouTube • Inner product space at Wikipedia • Cauchy–Schwarz inequality at Wikipedia • Gram–Schmidt process at Wikipedia Extra textbook problems: ? 1.3.3, 1.3.4, 1.4.1 ?? 1.3.9, 1.3.10, 1.3.12, 1.3.13, 1.4.2, 1.4.5(a,d) ??? 1.3.11, 1.3.14, 1.3.15, 1.3.25, 1.4.16 A 1.3.18 1 Advanced Linear Algebra – Week 5 2 There are many times when we would like to be able to talk about the angle between vectors in a vector space V, and in particular orthogonality of vectors, just like we did in Rn in the previous course. This requires us to have a generalization of the dot product to arbitrary vector spaces. Definition 5.1 — Inner Product Suppose that F = R or F = C, and V is a vector space over F. Then an inner product on V is a function h·, ·i : V × V → F such that the following three properties hold for all c ∈ F and all v, w, x ∈ V: a) hv, wi = hw, vi (conjugate symmetry) b) hv, w + cxi = hv, wi + chv, xi (linearity in 2nd entry) c) hv, vi ≥ 0, with equality if and only if v = 0.
    [Show full text]
  • Introduction to Clifford's Geometric Algebra
    解説：特集 コンピューテーショナル・インテリジェンスの新展開 —クリフォード代数表現など高次元表現を中心として— Introduction to Clifford’s Geometric Algebra Eckhard HITZER＊ ＊Department of Applied Physics, University of Fukui, Fukui, Japan ＊E-mail: [email protected] Key Words：hypercomplex algebra, hypercomplex analysis, geometry, science, engineering. JL 0004/12/5104–0338 C 2012 SICE erty 15), 21) that any isometry4 from the vector space V Abstract into an inner-product algebra5 A over the field6 K can Geometric algebra was initiated by W.K. Clifford over be uniquely extended to an isometry7 from the Clifford 130 years ago. It unifies all branches of physics, and algebra Cl(V )intoA. The Clifford algebra Cl(V )is has found rich applications in robotics, signal process- the unique associative and multilinear algebra with this ing, ray tracing, virtual reality, computer vision, vec- property. Thus if we wish to generalize methods from tor field processing, tracking, geographic information algebra, analysis, calculus, differential geometry (etc.) systems and neural computing. This tutorial explains of real numbers, complex numbers, and quaternion al- the basics of geometric algebra, with concrete exam- gebra to vector spaces and multivector spaces (which ples of the plane, of 3D space, of spacetime, and the include additional elements representing 2D up to nD popular conformal model. Geometric algebras are ideal subspaces, i.e. plane elements up to hypervolume ele- to represent geometric transformations in the general ments), the study of Clifford algebras becomes unavoid- framework of Clifford groups (also called versor or Lip- able. Indeed repeatedly and independently a long list of schitz groups). Geometric (algebra based) calculus al- Clifford algebras, their subalgebras and in Clifford al- lows e.g., to optimize learning algorithms of Clifford gebras embedded algebras (like octonions 17))ofmany neurons, etc.
    [Show full text]
  • A Bit About Hilbert Spaces
    A Bit About Hilbert Spaces David Rosenberg New York University October 29, 2016 David Rosenberg (New York University ) DS-GA 1003 October 29, 2016 1 / 9 Inner Product Space (or “Pre-Hilbert” Spaces) An inner product space (over reals) is a vector space V and an inner product, which is a mapping h·,·i : V × V ! R that has the following properties 8x,y,z 2 V and a,b 2 R: Symmetry: hx,yi = hy,xi Linearity: hax + by,zi = ahx,zi + b hy,zi Postive-definiteness: hx,xi > 0 and hx,xi = 0 () x = 0. David Rosenberg (New York University ) DS-GA 1003 October 29, 2016 2 / 9 Norm from Inner Product For an inner product space, we define a norm as p kxk = hx,xi. Example Rd with standard Euclidean inner product is an inner product space: hx,yi := xT y 8x,y 2 Rd . Norm is p kxk = xT y. David Rosenberg (New York University ) DS-GA 1003 October 29, 2016 3 / 9 What norms can we get from an inner product? Theorem (Parallelogram Law) A norm kvk can be generated by an inner product on V iff 8x,y 2 V 2kxk2 + 2kyk2 = kx + yk2 + kx - yk2, and if it can, the inner product is given by the polarization identity jjxjj2 + jjyjj2 - jjx - yjj2 hx,yi = . 2 Example d `1 norm on R is NOT generated by an inner product. [Exercise] d Is `2 norm on R generated by an inner product? David Rosenberg (New York University ) DS-GA 1003 October 29, 2016 4 / 9 Pythagorean Theroem Definition Two vectors are orthogonal if hx,yi = 0.
    [Show full text]
  • True False Questions from 4,5,6
    (c)2015 UM Math Dept licensed under a Creative Commons By-NC-SA 4.0 International License. Math 217: True False Practice Professor Karen Smith 1. For every 2 × 2 matrix A, there exists a 2 × 2 matrix B such that det(A + B) 6= det A + det B. Solution note: False! Take A to be the zero matrix for a counterexample. 2. If the change of basis matrix SA!B = ~e4 ~e3 ~e2 ~e1 , then the elements of A are the same as the element of B, but in a different order. Solution note: True. The matrix tells us that the first element of A is the fourth element of B, the second element of basis A is the third element of B, the third element of basis A is the second element of B, and the the fourth element of basis A is the first element of B. 6 6 3. Every linear transformation R ! R is an isomorphism. Solution note: False. The zero map is a a counter example. 4. If matrix B is obtained from A by swapping two rows and det A < det B, then A is invertible. Solution note: True: we only need to know det A 6= 0. But if det A = 0, then det B = − det A = 0, so det A < det B does not hold. 5. If two n × n matrices have the same determinant, they are similar. Solution note: False: The only matrix similar to the zero matrix is itself. But many 1 0 matrices besides the zero matrix have determinant zero, such as : 0 0 6.
    [Show full text]
  • Eigenvalues, Eigenvectors and the Similarity Transformation
    Eigenvalues, Eigenvectors and the Similarity Transformation Eigenvalues and the associated eigenvectors are ‘special’ properties of square matrices. While the eigenvalues parameterize the dynamical properties of the system (timescales, resonance properties, amplification factors, etc) the eigenvectors define the vector coordinates of the normal modes of the system. Each eigenvector is associated with a particular eigenvalue. The general state of the system can be expressed as a linear combination of eigenvectors. The beauty of eigenvectors is that (for square symmetric matrices) they can be made orthogonal (decoupled from one another). The normal modes can be handled independently and an orthogonal expansion of the system is possible. The decoupling is also apparent in the ability of the eigenvectors to diagonalize the original matrix, A, with the eigenvalues lying on the diagonal of the new matrix, . In analogy to the inertia tensor in mechanics, the eigenvectors form the principle axes of the solid object and a similarity transformation rotates the coordinate system into alignment with the principle axes. Motion along the principle axes is decoupled. The matrix mechanics is closely related to the more general singular value decomposition. We will use the basis sets of orthogonal eigenvectors generated by SVD for orbit control problems. Here we develop eigenvector theory since it is more familiar to most readers. Square matrices have an eigenvalue/eigenvector equation with solutions that are the eigenvectors xand the associated eigenvalues : Ax = x The special property of an eigenvector is that it transforms into a scaled version of itself under the operation of A. Note that the eigenvector equation is non-linear in both the eigenvalue () and the eigenvector (x).
    [Show full text]
  • A Guided Tour to the Plane-Based Geometric Algebra PGA
    A Guided Tour to the Plane-Based Geometric Algebra PGA Leo Dorst University of Amsterdam Version 1.15{ July 6, 2020 Planes are the primitive elements for the constructions of objects and oper- ators in Euclidean geometry. Triangulated meshes are built from them, and reflections in multiple planes are a mathematically pure way to construct Euclidean motions. A geometric algebra based on planes is therefore a natural choice to unify objects and operators for Euclidean geometry. The usual claims of `com- pleteness' of the GA approach leads us to hope that it might contain, in a single framework, all representations ever designed for Euclidean geometry - including normal vectors, directions as points at infinity, Pl¨ucker coordinates for lines, quaternions as 3D rotations around the origin, and dual quaternions for rigid body motions; and even spinors. This text provides a guided tour to this algebra of planes PGA. It indeed shows how all such computationally efficient methods are incorporated and related. We will see how the PGA elements naturally group into blocks of four coordinates in an implementation, and how this more complete under- standing of the embedding suggests some handy choices to avoid extraneous computations. In the unified PGA framework, one never switches between efficient representations for subtasks, and this obviously saves any time spent on data conversions. Relative to other treatments of PGA, this text is rather light on the mathematics. Where you see careful derivations, they involve the aspects of orientation and magnitude. These features have been neglected by authors focussing on the mathematical beauty of the projective nature of the algebra.
    [Show full text]
  • The Determinant Inner Product and the Heisenberg Product of Sym(2)
    INTERNATIONAL ELECTRONIC JOURNAL OF GEOMETRY VOLUME 14 NO. 1 PAGE 1–12 (2021) DOI: HTTPS://DOI.ORG/10.32323/IEJGEO.602178 The determinant inner product and the Heisenberg product of Sym(2) Mircea Crasmareanu (Dedicated to the memory of Prof. Dr. Aurel BEJANCU (1946 - 2020)Cihan Ozgur) ABSTRACT The aim of this work is to introduce and study the nondegenerate inner product < ·; · >det induced by the determinant map on the space Sym(2) of symmetric 2 × 2 real matrices. This symmetric bilinear form of index 2 defines a rational symmetric function on the pairs of rays in the plane and an associated function on the 2-torus can be expressed with the usual Hopf bundle projection 3 2 1 S ! S ( 2 ). Also, the product < ·; · >det is treated with complex numbers by using the Hopf invariant map of Sym(2) and this complex approach yields a Heisenberg product on Sym(2). Moreover, the quadratic equation of critical points for a rational Morse function of height type generates a cosymplectic structure on Sym(2) with the unitary matrix as associated Reeb vector and with the Reeb 1-form being half of the trace map. Keywords: symmetric matrix; determinant; Hopf bundle; Hopf invariant AMS Subject Classification (2020): Primary: 15A15 ; Secondary: 15A24; 30C10; 22E47. 1. Introduction This short paper, whose nature is mainly expository, is dedicated to the memory of Professor Dr. Aurel Bejancu, who was equally a gifted mathematician and a dedicated professor. As we hereby present an in memoriam paper, we feel it would be better if we choose an informal presentation rather than a classical structured discourse, where theorems and propositions are followed by their proofs.
    [Show full text]