<<

Quadratic Forms

Math 422

Definition 1 A is a function f : Rn R of form → f (x)=xT Ax, where A is an n n symmetric . × 2 3 x Example 2 f (x, y)=2x2 +3xy 4y2 = xy 2 . − 3 4 y ∙ 2 − ¸ ∙ ¸ £ ¤ Note that the Euclidean inner product () of two (column) vectors a and b can be expressed in terms of matrix multiplication as a, b = bT a. h i Thus, a quadratic form can be expressed in terms of the Euclidean inner product as

xT Ax = Ax, x = x,Ax . h i h i n 1 n Let S − denote the (n 1)-dimensional sphere in R , i.e., relative to the Euclidean inner product − n 1 n S − = x R : x, x =1 . { ∈ h i } n 1 n n 1 Since S − is a closed and bounded subset of R , continuous functions on S − attain their maximum and minimum values.

n 1 T Question #1: For x S − , what are the maximum and minimum values of a quadratic form x Ax? ∈

Theorem 3 Let A be a symmetric n n matrix with eigenvalues λ1 λ2 λn. Then × ≥ ≥ ···≥ T n 1 1. λ1 x Ax λn for all x S − . ≥ ≥ ∈ n 1 T 2. If x1 S − is an eigenvalue associated with λ1, then λ1 = x Ax1. ∈ 1 n 1 T 3. If xn S − is an eigenvalue associated with λn, then λn = x Axn. ∈ n The maximum and minimum of a quadratic form xT Ax can be found by computing the largest and smallest eigenvalue of A. The maximum (respectively, minimum) will always be attained at diametrically opposite x points on the x ,wherex is any eigenvector associated with λ1 (respectively, λn). ± k k

01 x1 Example 4 Consider f (x1,x2)=2x1x2 = x1 x2 . Since the eigenvalues of A = 10 x2 ∙ ¸ ∙ ¸ 01 £ ¤ 1 are λ1 =1and λ2 = 1, the maximum and minimum values of f on the unit circle S are 10 − ∙ ¸ 1 and 1, respectively. Furthermore, the maximum value is attained at the eigenvectors on S1 associated − 1/√2 1 with λ1 =1, namely ; the minimum value is attained at the eigenvectors on S associated with ± 1/√2 ∙ ¸ 1/√2 λ2 = 1, namely − . − ± 1/√2 ∙ ¸

Question #2: Under what conditions is the quadratic form xT Ax >0 for all x = 0? 6

1 Definition 5 xT Ax is positive definite iff xT Ax >0 for all x = 0 and xT Ax =0iff x = 0. A symmetric 6 matrix is positive definite iff xT Ax is positive definite.

Example 6 The Euclidean inner product is a positive definite quadratic form since

x2 + + x2 = x, x = xT x = xT Ix. 1 ··· n h i Theorem 7 A A is positive definite iff all eigenvalues of A are positive.

422 Example 8 The matrix A = 242 is positive definite since the eigenvalues of A are λ1 =8,λ2 =2 ⎡ 224⎤ T 2 2 2 and λ3 =2. Note that if x = 0⎣, then x Ax⎦ =2x +4x +4x +4x1x2 +4x1x3 +4x2x3 > 0. 6 1 2 3 th Definition 9 For 1 k n, the k principal submatrix of an n n matrix A =[aij] is ≤ ≤ ×

a11 a1k ··· . . ⎡ . . ⎤ . ak1 akk ⎢ ··· ⎥ ⎣ ⎦ Theorem 10 A symmetric matrix A is positive definite iff every principal subdeterminant of A is positive.

2 1 3 − − 2 1 Example 11 The principal subdeterminants of the matrix A = 12 4 are det [2] = 2, det − = ⎡ − ⎤ 12 34 9 ∙ − ¸ − 3 and det A =1. Since all are positive, the quadratic form xT A⎣x is positive definite.⎦

Question #3: If xT Ax is a quadratic form with non-diagonal A, under what conditions does there exist an orthogonal change variables x =P y so that (P y)T A (P y)=yT P T AP y has no cross-terms? ¡ ¢ Definition 12 An n n matrix A is orthogonally diagonalizable iff there exists an P such × that P T AP is a .

Theorem 13 If A is an n n matrix, then the following are equivalent: × 1. A is orthogonally diagonalizable. 2. A has an orthonormal set of n eigenvectors. 3. A is symmetric.

Theorem 14 If A is symmetric, then

1. The eigenvalues of A are real . 2. Eigenvectors from different eigenspaces are orthogonal with respect to the Euclidean inner product.

Use the following procedure to orthogonally diagonalize A:

Example 15 1. Find a basis for each eigenspace of A. 2. Apply Gram-Schmidt and obtain an orthonormal basis for each eigenspace. 3. Form the matrix P whose columns are the basis vectors constructed in step 2.

2 422 Example 16 Consider the matrix A = 242 whose eigenvalues are λ1 =8,λ2 =2and λ3 =2. ⎡ 224⎤ ⎣ ⎦ 1 Canonical bases for the eigenspaces are x1 = 1 ⎧ ⎡ ⎤⎫ ⎨ 1 ⎬ 1 1 ⎣ ⎦ − − ⎩ ⎭ 8 and x2 = 0 , x3 = 1 2. Note that x1, x2 = x1, x3 =0. Applying Gram-Schmidt ↔ ⎧ ⎡ ⎤ ⎡ ⎤⎫ ↔ h i h i ⎨ 1 0 ⎬ ⎣ ⎦ ⎣ ⎦ 1/√3 gives ⎩ ⎭v = x1 = 1/√3 and 1 x1 ⎧ k k ⎡ ⎤⎫ ⎨ 1/√3 ⎬ 1/√2 ⎣1/√6 ⎦ 1/√3 1/√2 1/√6 x2 − x3 x⎩3,v2 v2 − ⎭ − − v2 = = 0 , v3 = −h i = 2/√6 . The matrix P = 1/√30 2/√6 ⎧ x2 ⎡ ⎤ x3 x3,v2 v2 ⎡ ⎤⎫ ⎡ ⎤ k k 1/√2 k −h i k 1/√6 1/√31/√2 1/√6 ⎨ − ⎬ − ⎣ 800⎦ ⎣ ⎦ ⎣ ⎦ and⎩ P T AP = 020 is a diagonal matrix. ⎭ ⎡ 002⎤ ⎣ ⎦ T Theorem 17 Let x Ax be a quadratic form in variables x1,...,xn. Let P be an orthogonal matrix that orthogonally diagonalizes A. If λ1,...,λn are the eigenvalues of A and y1,...,yn are new variables such that x = P y, then T T T 2 2 x Ax = y P AP y = λ y + + λny 1 1 ··· n and ¡ ¢ λ1 0 0 ··· 0 λ2 0 T P AP = ⎡ . . . ⎤ . . .. . ⎢ ⎥ ⎢ 00 λn ⎥ ⎢ ··· ⎥ ⎣ ⎦ 422 2 2 2 Example 18 Let A = 242 . By the calculations in Example 16, 2x1 +4x2 +4x3 +4x1x2 +4x1x3 + ⎡ 224⎤ T T T 2 2 2 4x2x3 = x Ax = y P⎣ AP y =8⎦y1 +2y2 +2y3. ¡ ¢ Question #4: If xT Ax is a quadratic form in two or three variables and c is a constant, what does the graph of the level set xT Ax =c look like?

Theorem 19 If xT Ax is a quadratic form in two variables and c is a constant, the level curve given by xT Ax =c is a conic. If xT Ax is a quadratic form in three variables and c is a constant, the level surface given by xT Ax =c is a .

1/√2 1/√2 Example 20 In Example 4, let P = − ; then P T AP = 1/√21/√2 ∙ ¸ 10 2 2 . The level curve given by 2x1x2 =1is the y y =1since 0 1 1 − 2 ∙ − ¸

T T 10 2 2 2x1x2 = x Ax = y y = y y . 0 1 1 − 2 ∙ − ¸ 2 2 2 From Example 18 we observe that the level surface 2x1 +4x2 +4x3 +4x1x2 +4x1x3 +4x2x3 =1is the 2 2 2 8y1 +2y2 +2y3 =1.

3 Exercise 21 Since the quadratic form in Example 11 is positive definite, the quadric given by xT Ax =1 is an ellipsoid. Eliminate the cross-terms by performing an orthogonal change of variables. Express this 2 2 2 y1 y2 y3 ellipsoid in the standard form a2 + b2 + c2 =1. Definition 22 A quadratic form xT Ax is non-degenerate if all eigenvalues of A are non-zero.

Definition 23 The signature of a non-degenerate quadratic form xT Ax, denoted by sig (A) , is the of negative eigenvalues of A.

Theorem 24 Let xT Ax be a non-degenerate quadratic form in two variables.

1. If sig (A)=0, then xT Ax =1 is an . 2. If sig (A)=1,thenxT Ax =1 is an hyperbola.

Theorem 25 Let xT Ax be a non-degenerate quadratic form in three variables.

1. If sig (A)=0, then xT Ax =1 is an ellipsoid. 2. If sig (A)=1,thenxT Ax =1 is an of one sheet. 3. If sig (A)=2,thenxT Ax =1 is an hyperboloid of two sheets.

11-17-08

4