TBP MATH33A Review Sheet. November 24, 2018

Transcription

1 TBP MATH33A Review Sheet November 24, 2018

2 General Transformation Matrices: Function Scaling by k Orthogonal projection onto line L Implementation If we want to scale I 2 by k, we use the following: [ k 0 k I 2 = 0 k Given a unit vector u that is parallel to L such that [ u u1 = u 2 The matrix [ 2 u1 u 1 u 2 u 1 u 2 2 u 2 Reflection about a line Rotation through angle θ generates an orthogonal projection onto line Lwhere u u 2 2 = 1. In order to reflect about a line, we utilize a matrix of the form [ a b b -a where a 2 + b 2 = 1. Rotation through an angle θ involves a matrix of the form [ cosθ -sinθ sinθ cosθ or, more generally, [ a -b b a Rotation through angle θ combined with scaling by r where a 2 + b 2 = 1. This matrix takes the form [ a -b b a which is equivalent to [ cosθ -sinθ r sinθ cosθ Horizontal shear Matrices of the form [ 1 k 0 1 where k is an arbitrary constant. 1

3 Vertical shear Matrices of the form [ 1 0 k 1 where k is an arbitrary constant. Chapter 1 Rank shows dimension of im(a) where A is an arbitrary matrix If rank(a ) = rank(a), where Aĩs the augmented matrix of A, then the system is solvable and the dimension of the solution space is equal to the number of free variables. Furthermore, if, for the mxn matrix, m < n and the system is solvable, there are infinitely many solutions A transformation T is linear iff Chapter 2 T( v + w) = T( v ) + T( w) for all vectors v and w in IR m T(k v ) = kt( v ) for all vectors v in IR m and all scalars k Earlier, we defined an orthogonal projection onto line L where the unit vector u is on the line. We now define the more general case with a vector u that is on the line but not a unit vector. In this case [ w w1 = w 2 and the matrix is 1 w w 2 2 [ 2 w1 w 1 w 2 w 1 w 2 2 w 2 Furthermore, a reflection with respect to line L with an arbirtrary vector w on L as defined before can be found using 1 [ 2 2 w1 - w 2 2w 1 w 2 w w 2 2w 1 w 2 w w 1 which more generally takes the form [ a b b a where, as before, a 2 + b 2 = 1. Generally, for subspace L and V s.t. V =L, Ref V ( x) = -Ref L ( x) 2

4 The projection matrix onto a space where u is a unit vector in the space of interest ()and defined in the same way as before) is found by computing the matrix [ u1 u u2 u for the 2D case and [ u1 u u2 u u3 u Finally, we define that the matrix for reflection about a subspace V is equal to the identity matrix of that dimension minus the projection matrix. Chapter 3 A R V = I d A P V We begin by defining a subspace as something that is closed under addition and allows for scalar multiplication. Furthermore, if T ( x) = A x is a linear transformation from IR m to IR n, then ker(t ) = ker(a) is a subspace of IR m im(t ) = im(a) is a subspace of IR n We define a basis to be a set of vectors in substance V of IR n s.t. B = ( v 1, v 2, v 3,... v m ) Any vector x in V can be written as x = c 1 v 1 + c 1 v 2 + c 1 v v m where c 1, c 2, c 3,...c m are the β coordinates of x and the vector c 1 c 2 c 3... c m is defined as the β coordinate vector of x = [ x β We see x = B [ x β To get a better idea of the relationship between transformations and bases, we refer to the following diagram from page 152 of the Bretscher textbook. 3

5 To explain this diagram further, we start at the top left corner and consider an arbitrary vector x with a basis as described above that is on an arbitrary line L. In order to describe the vector on the line, we need to use two basis vectors and two constants. Looking at the bottom left of the diagram, we see that the coordinate vector of x is defined with two different coefficients. However, let us consider now x the basis so that one of the basis vectors v 1 is on L. This reduces the coordinate vector so that only the coefficient c 1 is relevant. The equivalent transformation for this change in basis is shown so that T transforms [ c1 c 2 into [ c1 0 We observe a parallelism in the transformation matrix and the coordinate vector. Chapter 5 We define a set of orthonormal vectors to be a set of all unit vectors that are all perpendicular to each other and are linearly independent. We also define a transform of A with notation A T such that, if [ A = then A T = Consequently, v T = v w. The orthogonal projection of v T onto V is denoted as v 4

6 if V is a subspace of IR n with orthonormal basis u 1, u 2, u 3,... u m proj V x = v = ( u 1 x) u ( u m x) u m Cauchy-Schwartz inequality: x y x y Gram-Schmidt process: convert an old base to a new base First, we consider a basis with vectors v 1, v 2,... v m and a new basis with vectors found by where we define u 1 = 1 v 1 u 1, u 2,... u m v 1, u 2 = 1 v 2 v 2,... u m = 1 v m v j = v j - ( u 1 x j ) u ( u j - 1 x j ) u j - 1 Q-R factorization is a method we use to convert a matrix R from the old base to a new base. We use the equation M = QR where M is a matrix with the basis vectors v 1, v 2,... v m as its columns and Q is a matrix with the unit basis vectors u 1, u 2,... u m as its columns. This equation can also be rewritten as R = Q T A Note that, in using this, Q must have orthonormal columns and R is an upper triangular matrix with positive diagonals. A transformation is classified as an orthogonal transformation if it preserves length If T ( x) = A x is an orthogonal transformation, A is an orthogonal matrix. In general, an orthogonal matrix is invertable, its determinant is not equal to zero, the matrix multiplied by its transpose is the identity matrix, and the inverse of the matrix is equal to its transpose. The matrix of an orthogonal projection is P = QQ T where General properties of transpose: (A + B) T = A T + B T (ka) T = ka T Q = [ u1 u2... u m v m 5

7 (AB) T = B T A T rank(a T ) = rank(a) (A T ) -1 = (A -1 ) T (Im(A)) = ker(a T ) Least square solution (A T A) x = A T b If A T A is invertible/linearly ind, then x = (A T A) -1 A T b and is unique For A x = b, A x is the orthogonal projection of b onto the image of A. The matrix of an orthogonal projection is Chapter 6 A(A T A) -1 A T b The determinant of an (upper or lower) triangular matrix is the product of the diagonal entries of the matrix. det(a T ) = det(a) Determinant operations: Chapter 7 Swapping rows results in x ( 1) where x is the determinant Dividing a row by scalar k results in x ( 1 k ) det(a) A square matrix is invertivle if and only if det(a) = 0 To make it easier to solve for the determinant of a given matrix A, we turn A into some upper triangular matrix B and find that det(a) = ( 1) s k 1 k 2... k r det(b) where s is the number of row swaps, k is a scalar that we divide a row in A by, and r is the number of scalars we the rows of A by. det(ab) = (det(a)) (det(b)) det(a m ) = (det(a)) m det(a -1 ) = 1 det(a) = (det(a))-1 An eigenvector is a vector v such that where λ is an eigenvalue An eigenbasis is a set of eigenvectors A v = λ v 6

8 [ #» v1 #» v2 #» v3... v #» n B = S -1 AS or A = SBS -1 where The trace of a matrix is defined as S = [ #» v1 #» v2 #» v3... v #» n λ B = 0 λ λ n tr(a) = a 11 +a a nn where a ii is the value of matrix at row i and column i. The roots of f A (λ) = det(a λi). An eigenspace is defined as E λ = ker(a λi n ) or { v in IR n : A v = λ v } Consequently, dim(e λ ) = geometric multiplicity. Complex Numbers z = a + bi = r(cosθ + sintθi) e iθ = cosθ + sinθi Thus, z = re iθ z= a - bi z * z = a 2 + b 2 = z 2 = r 2 If A has eigenvalues a ± ib where b 0 and v + i w is an eigenvector of A with eigenvalue a ± ib, then S -1 AS = [ a b b a where S = [ w v 7

9 Alternatively, [ a b A = S b a where S = [ v1 v2 and v 1 + v 2 i is an eigenvector for a ± ib. Chapter 8 A symmetric matrix is defined as A = A T. If A nn is symmetric matrix, There exists an orthonormal eigenbases (and thus the matrix is orthogonally diagonalizable). All eigenvalues are real ( λ = λ), there are n real eigenvalues, and an algebraic multiplicity of n. x T (T y) = y T (T x) For any orthonormal basis β, [T β is symmetric. If v 1 and v 2 are eigenvectors, of A with distinct eigenvalues λ 1 and λ 2, then v 1 v 2 = 0 and v 1 is orthogonal to v 2. To perform an orthogonal diagonalization of a symmetric matrix A, Find the eigenvalues of A, and find a basis of each eigenspace. Use the Gram-Schmidt process, find an orthonormal basis of each eigenspace. Forma n orthonormal eigenbases v 1, v 2, v 3,..., v #» n for A by concatenating the orthonormal basis you found in the second part such that S = [ v1 v2 v3... v #» n S -1 From this, S is orthogonal and S -1 AS will be diagonal. Quadratic forms: q( x) = x A x = x T A x = λ 1 c λ 2 c λ n c 2 n where we define A to be a symmetric nxn matrix, β to be an orthonormal eigenbasis for A with associated λ 1 c 2 1, λ 2 c 2 2,..., λ n c 2 n, and c i to be the coordinates of x with respect to β. Orthonormal coordinate transformation such that x = Sy to y = S -1 x = S T x Q( #» y 1, #» y 2,... #» y n ) = λ 1 y λ 2 y λ n y 2 n 8

10 A is positive definite if q( x) is positive for all nonzero x in IR n and iff det(a) (m) 0 for all m = 1, 2,..., n. This implies that all eigenvalues are positive. positive semidefinite if q( x) 0 for all x in IR n. This implies that all eigenvalues are positive or zero. indefinite if q( x) takes positive as well as negative values. The principal axes for q( x) = x A x with A being a symmetric nxn matrix with n distinct eigenvalues are the eigenspaces of A. Conside the curve defined by q( x 1, x 2 ) = a x b x #» #» 1 x 2 + c x #» 2 = 1 with eigenvalues λ 1, λ 2 of matrix [ a b/2 b/2 c of q. If both λ 1, λ 2 are positive, the curve is an ellipse, if one is positive and one is negative, the curve is a hyperbola. A linear transformation is described as L( x) = A x and is an invertivle linear transformation from IR 2 to IR 2 where A T A is symmetric, the image of the unit circle under L is an ellipse E, and the lengths of semimajor and semiminor axes of E are σ 1 and σ 2 of A. We define the following equation: where A = UΣV T V is an orthonormal eigenbasis for A T A with respect to eigenvalues defined such that λ 1 λ 2... λ r 0 = λ r = λ m with corresponding singular values (defined such athat singular values are the square root of λ) σ 1 σ 2... σ r 0... We associate the columns of V to be the orthonormal vectors v 1, v 2, v 3,... v #» n in IR n such that L( v 1 ), L( v 2 ), L( v 3 ),...L( v #» n ) are orthogonal. Σ = σ σ σ r U = [ u #» 1... u #» 2 ur + 1 = where u i are ybut vectors defined by u #» 1 = 1 σ 1 A v 1... u #» r = 1 σ r A v r such that A v i = σ 1 ui for i = 1, 2, 3...r and A v i = 0 for i = r + 1,...m. Furthermore, if A is an n x m matrix of rank r, then singular values σ 1, σ 2,... σ r are nonzero while σ r + 1,... σ m are zero. 9