← Back to Math 010

7.1 Orthogonal Matrices: Class Notes

Definition

A square matrix A is said to be orthogonal if its transpose is the same as its inverse, that is, if A−1=AT or, equivalently, if AAT=ATA=I. A matrix transformation TA:Rn→Rn is said to be an orthogonal transformation or an orthogonal operator if A is an orthogonal matrix.

Consider the matrix in Exercise 6.

A=[1/32/32/32/3−2/31/3−2/3−1/32/3]

Multiply A by its transpose. Each row of A becomes the matching column of AT, shown in the same color:

AAT=[1/32/32/32/3−2/31/3−2/3−1/32/3][1/32/3−2/32/3−2/3−1/32/31/32/3]=[100010001]

A row of A times the matching column of AT (the same row again) gives each 1 on the diagonal, and every other product is 0. So AAT=I, and A is orthogonal.

Theorem 7.1.1

The following are equivalent for an n×n matrix A.

a) A is orthogonal.

b) The row vectors of A form an orthonormal set in Rn with the Euclidean inner product.

c) The column vectors of A form an orthonormal set in Rn with the Euclidean inner product.

Tip: The rows (or columns) must be orthonormal, not just orthogonal.

Theorem 7.1.2

a) The transpose of an orthogonal matrix is orthogonal.

b) The inverse of an orthogonal matrix is orthogonal.

c) A product of orthogonal matrices is orthogonal.

d) If A is orthogonal, then det(A)=1 or det(A)=−1.

Proof of (d)

A orthogonal ⇒AAT=I, so

det(AAT)=det(I)⇒det(A)det(AT)=1

Since det(A)=det(AT),

[det(A)]2=1⇒det(A)=±1

#8

Let TA:R3→R3 be multiplication by the orthogonal matrix in Exercise 6. Find TA(𝐱) for the vector 𝐱=(0,1,4), and confirm ‖TA(𝐱)‖=‖𝐱‖ relative to the Euclidean inner product on R3.

A=[1/32/32/32/3−2/31/3−2/3−1/32/3]

Solution

‖𝐱‖=02+12+42=17

TA(𝐱)=[1/32/32/32/3−2/31/3−2/3−1/32/3][014]=[10/32/37/3]

‖TA(𝐱)‖=13100+4+49=13153=139·17=17

So ‖TA(𝐱)‖ matches ‖𝐱‖=17. This is always the case if A is orthogonal.

Answer: TA(𝐱)=[10/32/37/3],‖TA(𝐱)‖=‖𝐱‖=17

Theorem 7.1.3

If A is an n×n matrix, then the following are equivalent.

a) A is orthogonal.

b) ‖A𝐱‖=‖𝐱‖ for all 𝐱 in Rn [length is preserved].

c) A𝐱·A𝐲=𝐱·𝐲 for all 𝐱 and 𝐲 in Rn [the dot product is preserved].

Proof of (a) ⇒ (b)

A orthogonal ⇒AAT=ATA=I, so

‖A𝐱‖=(A𝐱·A𝐱)1/2=(𝐱TATA⏟I𝐱)1/2=(𝐱·𝐱)1/2=‖𝐱‖

Example 4: Rotations in the xy-plane

Sometimes a problem is more easily approached in R2 by considering a set of rotated axes, x′y′. Then a fixed point (x,y) has coordinates (x′,y′) relative to the rotated axes.

We start with 𝐮1, 𝐮2 (in particular, 𝐞1, 𝐞2), unit vectors along the x- and y-axes, and consider 𝐮1′, 𝐮2′ along the x′- and y′-axes. We'll treat the rotation as a change from an old basis B={𝐮1,𝐮2} to B′={𝐮1′,𝐮2′}.

Unit vectors along the x, y axes and along the x prime, y prime axes rotated by theta

Unit circle with the x′- and y′-axes drawn horizontal and vertical. The x- and y-axes sit at angle θ clockwise from them, since the x′y′-axes are the xy-axes rotated counterclockwise through θ. Unit vectors 𝐮1 and 𝐮2 point along the x- and y-axes, and 𝐮1′ and 𝐮2′ point along the x′- and y′-axes.

From 4.7,

[x′y′]=PB→B′[xy],wherePB→B′=[[𝐮1]B′[𝐮2]B′]

Relative to x′ and y′,

𝐮1=(cos⁡θ,−sin⁡θ),𝐮2=(sin⁡θ,cos⁡θ)

so

[𝐮1]B′=[cos⁡θ−sin⁡θ],[𝐮2]B′=[sin⁡θcos⁡θ]

So

[x′y′]=[cos⁡θsin⁡θ−sin⁡θcos⁡θ][xy]

⇒x′=xcos⁡θ+ysin⁡θ,y′=−xsin⁡θ+ycos⁡θ

These are the rotation equations for R2.

For R3 (rotating about the z-axis):

P=[cos⁡θsin⁡θ0−sin⁡θcos⁡θ0001]

This differs from the Rθ operator in 1.8 because Rθ rotates a vector to a new location. Here, the point is fixed; the axes are rotated.

Note that the rotation matrices are orthogonal matrices:

AT=A−1 ⇒ 𝐜i·𝐜j={1if i=j0if i≠j

AAT=I ⇒ 𝐫i·𝐫j={1if i=j0if i≠j

#14

Consider the rectangular x′y′-coordinate system obtained by rotating a rectangular xy-coordinate system counterclockwise through the angle θ=3π/4.

a. Find the x′y′-coordinates of the point whose xy-coordinates are (−2,6).

b. Find the xy-coordinates of the point whose x′y′-coordinates are (5,2).

Solution

The x y axes and the x prime y prime axes rotated by 3 pi over 4, with the point (minus 2, 6)

The xy-axes with the x′-axis (dashed) at angle θ=3π/4 counterclockwise from the positive x-axis, pointing up and to the left, and the y′-axis (dashed) pointing down and to the left. The point (−2,6) is in quadrant II of the xy-system.

With θ=3π/4, cos⁡θ=−1/2 and sin⁡θ=1/2.

a.

[x′y′]=[cos⁡θsin⁡θ−sin⁡θcos⁡θ][xy]

[x′y′]=[−1/21/2−1/2−1/2][−26]=[42−22]

Answer: (x′,y′)=(42,−22)

b. Solve for x and y with the inverse. The rotation matrix is orthogonal, so its inverse is its transpose: its first column becomes the first row and its second column becomes the second row.

[xy]=[cos⁡θsin⁡θ−sin⁡θcos⁡θ]−1[x′y′]=[cos⁡θ−sin⁡θsin⁡θcos⁡θ][x′y′]

[xy]=[−1/2−1/21/2−1/2][52]=[−7/23/2]

Answer: (x,y)=(−72,32)

Theorem 7.1.5

Let V be a finite-dimensional inner product space. If P is the transition matrix from one orthonormal basis for V to another orthonormal basis for V, then P is an orthogonal matrix.

Theorem 7.1.4

If S is an orthonormal basis for an n-dimensional inner product space V, and if (𝐮)S=(u1,u2,…,un) and (𝐯)S=(v1,v2,…,vn) then:

a) ‖𝐮‖=u12+u22+⋯+un2

b) d(𝐮,𝐯)=(u1−v1)2+(u2−v2)2+⋯+(un−vn)2

c) ⟨𝐮,𝐯⟩=u1v1+u2v2+⋯+unvn

#23

The set S={13, 12x, 32x2−23} is an orthonormal basis for P2 with respect to the evaluation inner product at the points x0=−1, x1=0, x2=1. Let 𝐩=p(x)=1+x+x2 and 𝐪=q(x)=2x−x2.

a. Find (𝐩)S and (𝐪)S.

b. Use Theorem 7.1.4 to compute ‖𝐩‖, d(𝐩,𝐪), and ⟨𝐩,𝐪⟩.

Solution

Call the basis vectors, in order, 𝐩1, 𝐩2, 𝐩3.

Part of (a) Since S is orthonormal,

(𝐩)S=(⟨𝐩,𝐩1⟩, ⟨𝐩,𝐩2⟩, ⟨𝐩,𝐩3⟩)

⟨𝐩,𝐩1⟩=p(−1)p1(−1)+p(0)p1(0)+p(1)p1(1)=(1)(13)+(1)(13)+(3)(13)=53

⟨𝐩,𝐩2⟩=p(−1)p2(−1)+p(0)p2(0)+p(1)p2(1)=(1)(−12)+(1)(0)+(3)(12)=2

⟨𝐩,𝐩3⟩=(1)(32−23)+(1)(−23)+(3)(32−23)=432−523=12−106=26=23

Answer: (𝐩)S=(53, 2, 23)

Part of (b)

‖𝐩‖=‖(𝐩)S‖=253+2+23=11

Answer: ‖𝐩‖=11