Solutions Q31 - 34#

# import all python add-ons etc that will be needed later on
%matplotlib inline
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.patches as patches
from sympy import *
init_printing()                      # allows printing of SymPy results in typeset maths format
plt.rcParams.update({'font.size': 14})  # set font size for plots

Q31 answer#

The coordinates of vector V1 are x1=rcos(α) and y1=rsin(α) and the second vector is rotated by θ from the first, so that V2 has x-coordinate

x2=rcos(αθ)=rcos(α)cos(θ)+rsin(α)sin(θ)

and by substitution x2=x1cos(θ)+y1sin(θ). The y-coordinate is

y2=rsin(αθ)=rcos(α)sin(θ)+rsin(α)cos(θ)

and y2=x1sin(θ)+y1cos(θ). Combining the two formulae produces the rotation matrix, equation 11.

[x2y2]=[cos(θ)sin(θ)sin(θ)cos(θ)][x1y1]

Q32 answer#

The rotation matrix is RRθ=[cos(θ)sin(θ)sin(θ)cos(θ)] and its square is

RRθ=[cos(θ)sin(θ)sin(θ)cos(θ)][cos(θ)sin(θ)sin(θ)cos(θ)]=[cos2(θ)sin2(θ)2sin(θ)cos(θ)2sin(θ)cos(θ)cos2(θ)sin2(θ)]

Substituting the double angle relationships into the last matrix shows that R(2θ)=R(θ)2.

Q33 answer#

There is no set answer to this question.

Q34 answer#

As rotation only occurs about the z-axis and by 90o,the product of the three rotation matrices becomes

[cos(θ)sin(θ)0sin(θ)cos(θ)0001]=[010100001]

An inversion is the matrix [100010001] and the rotation-inversion operation is

[100010001][010100001]=[010100001]

(b) Performing the matrix multiplication the other way round produces the same result, therefore Rot×Inv=Inv×Rot and the matrices commute.