Quadratic forms. Defn

Similar documents
Non-square diagonal matrices

Rotations & reflections

Rotations & reflections

Old & new co-ordinates

Jordan blocks. Defn. Let λ F, n Z +. The size n Jordan block with e-value λ is the n n upper triangular matrix. J n (λ) =

Data fitting. Ideal equation for linear relation

Rank & nullity. Defn. Let T : V W be linear. We define the rank of T to be rank T = dim im T & the nullity of T to be nullt = dim ker T.

Matrix-valued functions

Matrix-valued functions

Matrix of linear maps. Matrix-vector product

Definition of adjoint

A A x i x j i j (i, j) (j, i) Let. Compute the value of for and

Chapter 7: Symmetric Matrices and Quadratic Forms

Kernel. Prop-Defn. Let T : V W be a linear map.

i.e. the i-th column of A is the value of T at the i-th standard basis vector e i R n.

Orthogonal complement

Economics 204 Summer/Fall 2010 Lecture 10 Friday August 6, 2010

Here are some additional properties of the determinant function.

Week Quadratic forms. Principal axes theorem. Text reference: this material corresponds to parts of sections 5.5, 8.2,

Eigenvectors. Prop-Defn

Jordan chain. Defn. E.g. T = J 3 (λ) Daniel Chan (UNSW) Lecture 29: Jordan chains & tableaux Semester / 10

Chapter 12 Review Vector. MATH 126 (Section 9.5) Vector and Scalar The University of Kansas 1 / 30

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

Linear Algebra 20 Symmetric Matrices and Quadratic Forms

7. Symmetric Matrices and Quadratic Forms

Lecture 7: Positive Semidefinite Matrices

Chapter 5: Matrices. Daniel Chan. Semester UNSW. Daniel Chan (UNSW) Chapter 5: Matrices Semester / 33

Chapter 2: Vector Geometry

1 Last time: least-squares problems

Data for a vector space

Bilinear and quadratic forms

Final Exam, Linear Algebra, Fall, 2003, W. Stephen Wilson

Algebra Workshops 10 and 11

APPLICATIONS The eigenvalues are λ = 5, 5. An orthonormal basis of eigenvectors consists of

Linear Algebra II. Ulrike Tillmann. January 4, 2018

Lecture 1: Systems of linear equations and their solutions

8. Diagonalization.

MATH 20F: LINEAR ALGEBRA LECTURE B00 (T. KEMP)

Econ 204 Supplement to Section 3.6 Diagonalization and Quadratic Forms. 1 Diagonalization and Change of Basis

Math 3191 Applied Linear Algebra

Lecture 2: Linear Algebra Review

Math Matrix Algebra

MATH 1020 WORKSHEET 12.1 & 12.2 Vectors in the Plane

Review problems for MA 54, Fall 2004.

Lecture 15, 16: Diagonalization

complex dot product x, y =

Linear Algebra: Matrix Eigenvalue Problems

Transpose & Dot Product

ICS 6N Computational Linear Algebra Symmetric Matrices and Orthogonal Diagonalization

ANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3

Transpose & Dot Product

Exercise Sheet 1.

Chapter 1: Introduction to Vectors (based on Ian Doust s notes)

Section 7.3: SYMMETRIC MATRICES AND ORTHOGONAL DIAGONALIZATION

Symmetric matrices and dot products

Maths for Signals and Systems Linear Algebra in Engineering

Linear Algebra. Session 12

Chapter 6: Orthogonality

Section 6.2, 6.3 Orthogonal Sets, Orthogonal Projections

Fall River Joint Unified School District

x i e i ) + Q(x n e n ) + ( i<n c ij x i x j

Mathematical foundations - linear algebra

Math 520 Exam 2 Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

Announce Statistical Motivation Properties Spectral Theorem Other ODE Theory Spectral Embedding. Eigenproblems I

Linear Algebra, part 2 Eigenvalues, eigenvectors and least squares solutions

Throughout these notes we assume V, W are finite dimensional inner product spaces over C.

Chapter 13: Vectors and the Geometry of Space

Chapter 13: Vectors and the Geometry of Space

Chapter 6. Eigenvalues. Josef Leydold Mathematical Methods WS 2018/19 6 Eigenvalues 1 / 45

Linear Algebra Review. Vectors

Spectral Theorem for Self-adjoint Linear Operators

Quantum Computing Lecture 2. Review of Linear Algebra

We wish the reader success in future encounters with the concepts of linear algebra.

Algebra II. Paulius Drungilas and Jonas Jankauskas

The Singular Value Decomposition and Least Squares Problems

The value of a problem is not so much coming up with the answer as in the ideas and attempted ideas it forces on the would be solver I.N.

Math 290-2: Linear Algebra & Multivariable Calculus Northwestern University, Lecture Notes

MATH 16C: MULTIVARIATE CALCULUS

Here each term has degree 2 (the sum of exponents is 2 for all summands). A quadratic form of three variables looks as

Mathematical Methods wk 2: Linear Operators

OHSx XM521 Multivariable Differential Calculus: Homework Solutions 13.1

Part IA. Vectors and Matrices. Year

SYLLABUS. 1 Linear maps and matrices

Problem # Max points possible Actual score Total 120

Mathematical foundations - linear algebra

Computational math: Assignment 1

Ir O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v )

22m:033 Notes: 7.1 Diagonalization of Symmetric Matrices

Lecture 7. Econ August 18

Linear Algebra: Characteristic Value Problem

8 Eigenvectors and the Anisotropic Multivariate Gaussian Distribution

1 Planar rotations. Math Abstract Linear Algebra Fall 2011, section E1 Orthogonal matrices and rotations

EECS 275 Matrix Computation

Linear Least Square Problems Dr.-Ing. Sudchai Boonto

j=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent.

Algebra II. A2.1.1 Recognize and graph various types of functions, including polynomial, rational, and algebraic functions.

CS168: The Modern Algorithmic Toolbox Lecture #8: How PCA Works

Geometry and Motion, MA 134 Week 1

Announce Statistical Motivation ODE Theory Spectral Embedding Properties Spectral Theorem Other. Eigenproblems I

Algebra 2A Unit 1 Week 1 Day Activity Unit 1 Week 2 Day Activity Unit 1 Week 3 Day Activity Unit 2 Week 1 Day Activity

Transcription:

Quadratic forms Aim lecture: We use the spectral thm for self-adjoint operators to study solns to some multi-variable quadratic eqns. In this lecture, we work over the real field F = R. We use the notn x = x 1,..., x n ), y = y 1,..., y n ) R n Defn A real) quadratic form is a function Q : R n R of the form Qx) = i j β ijx i x j for some β ij R. A quadric in R n is a non-empty set of the form V Q) = {x R n Qx) = 1} for some given quadratic form Q i.e. V Q) is the set of solns to Qx) = 1. Rem If Q is a quadratic form & d R 0 such that the solns V to Qx) = d is non-empty, then V is the quadric V 1 d Q). E.g. For a ij R, consider the function Q x 1 ) x 2 ) = a11 a12 x1 x 2 ) a 21 a 22 ) x 1 ) x 2 = Daniel Chan UNSW) Lecture 40: Quadrics Semester 2 2013 1 / 10

Quadratic forms via symmetric matrices Quadratic forms can be studied via matrices because of the following Prop-Defn 1 Let A M nn R). Then Q A : R n R : x x T Ax is a quadratic form. 2 Let Q : R n R be a quadratic form. Then there is a unique symmetric matrix A M nn R) such that Q = Q A. In fact if A = a ij ) ij then Q A x) = i a ii x 2 i + a ij + a ji )x i x j. i<j Proof. Just generalise the calculation in the previous example. E.g. Daniel Chan UNSW) Lecture 40: Quadrics Semester 2 2013 2 / 10

Non-degenerate forms Recall that the e-values of a real symmetric matrix are all real even when considered as a complex matrix). Defn Let Q : R n R be a quadratic form & A M nn R) be the symmetric matrix with Q = Q A. Let τ +, τ be the number of positive resp negative e-values of A, counted with multiplicity alg = geom). We say Q or V Q)) is non-degenerate if A is invertible. Otherwise we say it is degenerate. The rank of Q is rankq = rank A = τ + + τ. The signature of Q is τ + τ. E.g. A degenerate quadric in R 2 looks like Daniel Chan UNSW) Lecture 40: Quadrics Semester 2 2013 3 / 10

Conics Defn A non-degenerate quadric in R 2 is called a conic. Hopefully, you already know the following examples. Below we let a 1, a 2 > 0. E.g. 1 If D = 1) 1) = I 2 then V Q D ) is the circle x1 2 + x 2 2 = 1 & Q D has signature 2. E.g. 2 If D = 1 ) 1 ) then V Q a1 2 a2 2 D ) is the ellipse x1 + x2 = 1 & Q D has signature 2. a 1 a 2 N.B. For a 1 a 2, the x 1 and x 2 -intercepts give the points closest and furthest from 0, 0). Daniel Chan UNSW) Lecture 40: Quadrics Semester 2 2013 4 / 10

Conics cont d E.g. 3 If D = 1 ) 1 ) then V Q a1 2 a2 2 D ) is the hyperbola x1 x2 = 1 a 1 a 2 & Q D has signature 0. Note the following geometric features: Closest points to 0,0): Asymptotes: Daniel Chan UNSW) Lecture 40: Quadrics Semester 2 2013 5 / 10

Quadric surfaces A quadric surface is a non-degenerate quadric in R 3. Let a 1, a 2, a 3 > 0. Prop If D = 1 ) 1 ) 1 ) then V Q a1 2 a2 2 a3 2 D ) is the ellipsoid x1 + x2 + x3 = 1. a 1 a 2 a 3 The points closest resp furthest) from 0, 0) found by setting x i = 0 for all non-minimal resp non-maximal) a i. The signature is 3. Prop If D = 1 ) 1 ) 1 ) then V Q a1 2 a2 2 a3 2 D ) is the hyperboloid of 1-sheet x1 + x2 x3 = 1. a 1 a 2 a 3 The pts closest to 0, 0) must lie on the ellipse x 3 = 0. The signature = 1. Daniel Chan UNSW) Lecture 40: Quadrics Semester 2 2013 6 / 10

Quadric surfaces cont d A quadric surface is a non-degenerate quadric in R 3. Let a 1, a 2, a 3 > 0. Prop If D = 1 ) 1 ) 1 ) then V Q a1 2 a2 2 a3 2 D ) is the hyperboloid of 2-sheets x1 x2 x3 = 1. a 1 a 2 a 3 The points closest to 0, 0) are the x 1 -intercepts. The signature is -1. Sketch Let s re-scale axes by changing variables to y i = x i a i. Then in R 3 y -space the surface is just the surface of revolution of Daniel Chan UNSW) Lecture 40: Quadrics Semester 2 2013 7 / 10

Principal axis theorem Let A M nn R) be a symm matrix. We apply the spectral thm for self-adjoint operators to orthog diagonalise D = U T AU where U O n & D = i λ i ). Let v i be the i-th column of U, a unit norm e-vector of A with e-value λ i. Consider now the orthonormal change of co-ord system U : R n y R n x & introduce new co-ords y = U 1 x = U T x. In our new co-ords, our quadratic form Q A becomes Q A x) = x T Ax = Uy) T UDU T )Uy) = y T U T UDU T Uy = y T Dy = Q D y). We obtain Theorem Principal axis) There is an orthonormal change of co-ords such that the quadric V Q A ) has the form V Q D ) for D the diagonal matrix of e-values. The new co-ord axes called principal axes go through the corresponding e-vectors v i. In particular, any conic is either an ellipse or an hyperbola depending on whether the signature of A is 2 or 0. Any quadric surface is either an ellipsoid, an hyperboloid of 1-sheet, or an hyperboloid of 2-sheets, depending on whether the signature is 3,1 or -1. Daniel Chan UNSW) Lecture 40: Quadrics Semester 2 2013 8 / 10

Example E.g. Describe the quadric 2x 2 1 + x 2 2 + x 2 3 + 4x 1 x 2 + 4x 1 x 3 + 8x 2 x 3 = 3 Find the closest point on the quadric to the 0, 0, 0) T. A We may re-write the quadric as x T Ax = 3 where 2 2 2 A = 2 1 4 2 4 1 We orthogonally diagonalised this last lecture to see A = UDU T where D = 6) 3) 3) & U O 3. Also the e-spaces were E 6 = Span1, 2, 2) T, E 3 = E6. Hence there s an orthonormal change of co-ords U : R 3 y R 3 x so that Daniel Chan UNSW) Lecture 40: Quadrics Semester 2 2013 9 / 10

Example cont d Daniel Chan UNSW) Lecture 40: Quadrics Semester 2 2013 10 / 10