Orthogonal Matrices

Orthogonal Matrices

Part 2. 1 Part 2. 2 Orthogonal Matrices • If u and v are nonzero vectors then u · v = kuk kvk cos(θ) ◦ £¥¤ ¢ is 0 if and only if cos(θ) = 0, i.e., θ = 90 . ¡ Hence, we say that two vectors u and v are perpendicular or orthogonal (in symbols u ⊥ v) if u · v = 0. • A vector u is a unit vector if kuk = 1. Now rotate the vectors, i.e., apply T . • Let T : R2 → R2 be rotation around the £ ¥ origin by θ degrees. This is a linear transformation. To see this, consider ¡£¢¥¤ the addition of u and v as in the §£¨ © following diagram. ¦ This diagram shows that T (u) + T (v) = T (u + v). It’s even easier to Part 2. 3 Part 2. 4 check that T (sv) = sT (v), so T is linear. • Recall the addition laws for sin and cos. We need to find T (e1), T (e2). If we cos(A + B) = cos(A) cos(B) − sin(A) sin(B) rotate e by θ, we get (cos(θ), sin(θ)) by 1 cos(A − B) = cos(A) cos(B) + sin(A) sin(B) the unit circle definition of sin and cos. sin(A + B) = sin(A) cos(B) + cos(A) sin(B) sin(A − B) = sin(A) cos(B) − cos(A) sin(B) £¥¤§¦¢¨ © Recall also that ¢¡ cos(−θ) = cos(θ) sin(−θ) = − sin(θ) The vector w = (− sin(θ), cos(θ)) is a unit • Exercise: Rotating by θ and then by ϕ vector and w · T (e1) = 0, so T (e2) must should be the same as rotating θ + ϕ. be either w or −w. Considering the sign Thus, we have of the first component shows that R(θ)R(ϕ) = R(θ + ϕ) = R(ϕ)R(θ). T (e2) = w. Use this and matrix multiplication to Thus, the matrix of T is derive the addition laws for sin and cos. −1 cos θ − sin(θ) What is R(θ) ?. R(θ) = . sin(θ) cos(θ) • Exercise: What is the matrix of reflection through the x-axis? Part 2. 5 Part 2. 6 • A matrix A is orthogonal if kAvk = kvk Hence, if A is orthogonal, we have for all vectors v. If A and B are 2Au · Av = kAuk2 + kAvk2 − kAu − Avk2 orthogonal, so is AB. If A is orthogonal, = kAuk2 + kAvk2 − kA(u − v)k2 so is A−1. Clearly I is orthogonal.. Rotation matrices are orthogonal. = kuk2 + kvk2 − ku − vk2 The set of orthogonal 2 × 2 matrices is = 2u · v denoted by O(2). so Au · Av = u · v. Classification of Orthogonal • If A is orthogonal, we must have Matrices 2 kAe1k = Ae1 · Ae1 = e1 · e1 = 1 2 • If A is orthogonal, then Au · Av = u · v (∗) kAe2k = Ae1 · Ae2 = e2 · e2 = 1 for all vectors u and v. (It follows that Ae1 · Ae2 = e1 · e2 = 0. A preserves the angles between vectors). To see this, recall that • Exercise: If A satisfies (∗), A is orthogonal. 2u · v = kuk2 + kv2k − ku − vk2. • Equations (∗) say that the columns of A are orthogonal unit vectors. Thus, if the first column of A is a 2 2 , a + b = 1 b Part 2. 7 Part 2. 8 the possibilities for the second column Call this matrix S(θ). What does this are represent geometrically? −b b , • We claim that there is a line L through a −a the origin so that A = S(θ) leaves every So A must look like one of the matrices point on L fixed. To see this, it will suffice to find a unit vector a −b a b , . cos(ϕ) b a b −a u = sin(ϕ) • Since a2 + b2 = 1, we can find θ so that so that Au = u. If we do this, any other a = cos(θ) and b = sin(θ). Thus, in one vector along L will have the form tu for case, we get some scalar t, and we will have A(tu) = tAu = tu. Since we only need to a −b cos θ − sin(θ) consider each line through the origin A = = = R(θ), b a sin(θ) cos(θ) once, we can suppose 0 ≤ ϕ < 180◦. We may as well suppose that 0 ≤ θ < 360◦. so A is rotation by θ degrees. In the Thus, we are trying to find ϕ so that other case we have Au = u, i.e. a b cos θ sin(θ) cos θ sin(θ) cos(ϕ) cos(ϕ) A = = = b −a sin(θ) − cos(θ) sin(θ) − cos(θ) sin(ϕ) sin(ϕ) Part 2. 9 Part 2. 10 If we carry out the multiplication, this which is a unit vector orthogonal to u. becomes Calculate Av as follows cos(θ) cos(ϕ) + sin(θ) sin(ϕ) cos(ϕ) cos θ sin(θ) − sin(θ/2) = Av = sin(θ) cos(ϕ) − cos(θ) sin(ϕ) sin(ϕ) sin(θ) − cos(θ) cos(θ/2) Using the addition laws, this is the same − cos(θ) sin(θ/2) + sin(θ) cos(θ/2) = as − sin(θ) sin(θ/2) − cos(θ) cos(θ/2) cos(θ − ϕ) cos(ϕ) = sin(θ − θ/2) sin(θ − ϕ) sin(ϕ) = − cos(θ − θ/2) Because of our angle restrictions, we must have θ − ϕ = ϕ, so ϕ = θ/2. Thus, sin(θ/2) = Au = u for − cos(θ/2) cos(θ/2) = −v, u = sin(θ/2). So Av = −v. A little thought shows that A is reflection through the line L that Consider the vector passes through the origin and is parallel to u. Any vector w can be written as − sin(θ/2) v = , w = su + tv for scalars s and t and cos(θ/2) A(su + tv) = sAu + tAv = su − tv. Part 2. 11 Part 2. 12 Isometries of the Plane £¥¤ ¦¨§¥© • If a point has coordinates (x, y), the ¢¡ vector with components (x, y) is the vector from the origin to the point. This is called the position vector of the point. Generally, we use the same notation for the point and its position vector. • Exercise: Verify the following • The distance between two points p and equations. q in R2 is given by d(p, q) = kp − qk. R(θ)S(ϕ) = S(θ + ϕ) • A transformation T : R2 → R2 is called S(ϕ)R(θ) = S(ϕ − θ) an isometry if it is 1-1 and onto and d(T (p),T (q)) = d(p, q) for all points p and S(θ)S(ϕ) = R(θ − ϕ). q. • One kind of isometry we have discovered is T (p) = Ap, where A is an orthogonal matrix. • Another kind of isometry is translation. Let v be a fixed vector and define Tv by Part 2. 13 Part 2. 14 Tv(p) = p + v. Translations are not linear. • If x, y and z are three noncollinear points and T is an isometry such that • The identity mapping id: 2 → 2 R R T (x) = x, T (y) = y and T (z) = z, then defined by id(p) = p is obviously an T = id. isometry. The composition of two To see this, suppose that is any point. isometries is an isometry. The inverse p Then is determined by the numbers of an isometry is an isometry. p d(x, p), d(y, p) and d(z, p). Since T is an • Let x, y and z be noncollinear points. isometry Then a point p is uniquely determined d(x, p) = d(T (x),T (p)) = d(x, T (p)). by the three numbers d(x, p), d(y, p) and Similarly, d(y, T (p)) = d(y, p) and d(z, p). d(z, T (p)) = d(z, p). Thus, we must have T (p) = p. Problem: If points , and form a ¢ • x y z right triangle with the right angle at x, £ then T (x), T (y) and T (z) form a right ¡ triangle with the right angle at T (x). Part 2. 15 Part 2. 16 • Theorem If T is an isometry, there is a let M be the matrix of reflection vector v and an orthogonal matrix A so through the x axis. In either case M is 00 that T (p) = Ap + v for all p, i.e., T is the an orthogonal matrix and MT (e2) = e2. 000 000 composition of a rotation or reflection Let T = MRTvT . Then T (0) = 0, 000 000 000 at the origin and a translation. T (e1) = e1 and T (e2) = e2, so T must be the identity. If we set , is Proof: The points 0, e1 and e2 are A = MR A 000 noncollinear and form a right triangle, orthogonal and we have T = ATvT = id. so the points T (0), T (e1) and T (e2) form Thus, for any point point p, −1 a right triangle. ATvT (p) = p. Multiplying by A we get T T (p) = A−1p. Apply T to both sides Let v = −T (0). Let Tv be translation by v −v 0 to get T (p) = T A−1p = A−1p − v. Since v. Then T = TvT is an isometry and −v −1 T 0(0) = 0. We must have A is orthogonal, the proof is 0 complete. d(0,T (e1)) = d(0, e1) = 1, so we can 0 rotate T (e1) around the origin so it • Every isometry can be written in the coincides with e1. Let R be the rotation form T (x) = Ax + u. To abbreviate this, matrix for this rotation. Then we’ll simply denote this isometry by 00 00 T = RTvT is an isometry with T (0) = 0 (A | u), so the rule is 00 00 and T (e1) = e1.

View Full Text

Details

  • File Type
    pdf
  • Upload Time
    -
  • Content Languages
    English
  • Upload User
    Anonymous/Not logged-in
  • File Pages
    12 Page
  • File Size
    -

Download

Channel Download Status
Express Download Enable

Copyright

We respect the copyrights and intellectual property rights of all users. All uploaded documents are either original works of the uploader or authorized works of the rightful owners.

  • Not to be reproduced or distributed without explicit permission.
  • Not used for commercial purposes outside of approved use cases.
  • Not used to infringe on the rights of the original creators.
  • If you believe any content infringes your copyright, please contact us immediately.

Support

For help with questions, suggestions, or problems, please contact us