跳到主要内容

明治大学 先端数理科学研究科 現象数理学専攻 2021年8月実施 线性代数

Author​

思齐塾, 祭音Myyura

Description​

R3\mathbb{R}^3 の座標を (x,y,z)(x, y, z) として、 x=rcos⁡θ,y=rsin⁡θx = r \cos \theta, y = r \sin \theta とおく。 R3\mathbb{R}^3 の点 (rcos⁡θ,rsin⁡θ,z)(r \cos \theta, r \sin \theta, z) を (rcos⁡(θ+α),rsin⁡(θ+α),z)(r \cos (\theta + \alpha), r \sin (\theta + \alpha), z) にうつす変換を zz 軸のまわりの α\alpha 回転とよび fαf_\alpha で表す。次の問に答えよ。

(1) fαf_\alpha は R3\mathbb{R}^3 の線形変換であることを示し、これを3次の正方行列で表せ。

(2) xx 軸のまわりの π/2\pi / 2 回転と zz 軸のまわりの π/2\pi / 2 回転を、それぞれ行列 A1A_1 と A2A_2 で表す。行列 A1A_1 および A2A_2 を求めよ。

(3) 行列の積 A1A2A_1 A_2 を MM とおく。 MM の固有値をすべて求めよ。また、行列 MM は R3\mathbb{R}^3 のあるベクトルのまわりの回転を表すことを示し、回転軸となるベクトルと回転角の大きさを求めよ。

题目描述​

在 R3\mathbb{R}^3 中使用坐标 (x,y,z)(x,y,z),并令

x=rcos⁡θ,y=rsin⁡θ.x=r\cos\theta,\qquad y=r\sin\theta.

把 R3\mathbb{R}^3 中的点

(rcos⁡θ,rsin⁡θ,z)(r\cos\theta,r\sin\theta,z)

映射为

(rcos⁡(θ+α),rsin⁡(θ+α),z)(r\cos(\theta+\alpha),r\sin(\theta+\alpha),z)

的变换称为绕 zz 轴旋转 α\alpha,记作 fαf_\alpha。回答下列问题。

(1) 证明 fαf_\alpha 是 R3\mathbb{R}^3 上的线性变换,并用三阶方阵表示该变换。

(2) 分别用矩阵 A1,A2A_1,A_2 表示绕 xx 轴旋转 π/2\pi/2 与绕 zz 轴旋转 π/2\pi/2。求 A1A_1 和 A2A_2。

(3) 令

M=A1A2.M=A_1A_2.

求 MM 的全部特征值。再证明矩阵 MM 表示 R3\mathbb{R}^3 中绕某个向量的旋转,并求作为旋转轴的向量以及旋转角的大小。

Kai​

(1) fα(rcos⁡θ,rsin⁡θ,z)=(rcos⁡(θ+α),rsin⁡(θ+α),z)f_\alpha(r\cos\theta, r\sin\theta, z) = (r\cos(\theta+\alpha), r\sin(\theta+\alpha), z) . Let v1=(r1cos⁡θ1,r1sin⁡θ1,z1)v_1 = (r_1\cos\theta_1, r_1\sin\theta_1, z_1) and v2=(r2cos⁡θ2,r2sin⁡θ2,z2)v_2 = (r_2\cos\theta_2, r_2\sin\theta_2, z_2) . Then fα(v1+v2)=fα(r1cos⁡θ1+r2cos⁡θ2,r1sin⁡θ1+r2sin⁡θ2,z1+z2)f_\alpha(v_1 + v_2) = f_\alpha(r_1\cos\theta_1 + r_2\cos\theta_2, r_1\sin\theta_1 + r_2\sin\theta_2, z_1+z_2) is difficult to analyze. Instead, consider fα(x,y,z)=(xcos⁡α−ysin⁡α,xsin⁡α+ycos⁡α,z)f_\alpha(x,y,z) = (x\cos\alpha - y\sin\alpha, x\sin\alpha + y\cos\alpha, z) . fα(x1+x2,y1+y2,z1+z2)=((x1+x2)cos⁡α−(y1+y2)sin⁡α,(x1+x2)sin⁡α+(y1+y2)cos⁡α,z1+z2)=(x1cos⁡α−y1sin⁡α,x1sin⁡α+y1cos⁡α,z1)+(x2cos⁡α−y2sin⁡α,x2sin⁡α+y2cos⁡α,z2)=fα(x1,y1,z1)+fα(x2,y2,z2)f_\alpha(x_1+x_2, y_1+y_2, z_1+z_2) = ((x_1+x_2)\cos\alpha - (y_1+y_2)\sin\alpha, (x_1+x_2)\sin\alpha + (y_1+y_2)\cos\alpha, z_1+z_2) = (x_1\cos\alpha - y_1\sin\alpha, x_1\sin\alpha + y_1\cos\alpha, z_1) + (x_2\cos\alpha - y_2\sin\alpha, x_2\sin\alpha + y_2\cos\alpha, z_2) = f_\alpha(x_1,y_1,z_1) + f_\alpha(x_2, y_2, z_2) fα(cx,cy,cz)=(cxcos⁡α−cysin⁡α,cxsin⁡α+cycos⁡α,cz)=c(xcos⁡α−ysin⁡α,xsin⁡α+ycos⁡α,z)=cfα(x,y,z)f_\alpha(cx, cy, cz) = (cx\cos\alpha - cy\sin\alpha, cx\sin\alpha + cy\cos\alpha, cz) = c(x\cos\alpha - y\sin\alpha, x\sin\alpha + y\cos\alpha, z) = c f_\alpha(x,y,z) So it is a linear transformation. The matrix representation is

(cos⁡α−sin⁡α0sin⁡αcos⁡α0001)\begin{pmatrix} \cos\alpha & -\sin\alpha & 0 \\ \sin\alpha & \cos\alpha & 0 \\ 0 & 0 & 1 \end{pmatrix}

(2) Rotation about x-axis by π/2\pi/2 is

A1=(1000cos⁡(π/2)−sin⁡(π/2)0sin⁡(π/2)cos⁡(π/2))=(10000−1010)A_1 = \begin{pmatrix} 1 & 0 & 0 \\ 0 & \cos(\pi/2) & -\sin(\pi/2) \\ 0 & \sin(\pi/2) & \cos(\pi/2) \end{pmatrix} = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & -1 \\ 0 & 1 & 0 \end{pmatrix}

Rotation about z-axis by π/2\pi/2 is

A2=(cos⁡(π/2)−sin⁡(π/2)0sin⁡(π/2)cos⁡(π/2)0001)=(0−10100001)A_2 = \begin{pmatrix} \cos(\pi/2) & -\sin(\pi/2) & 0 \\ \sin(\pi/2) & \cos(\pi/2) & 0 \\ 0 & 0 & 1 \end{pmatrix} = \begin{pmatrix} 0 & -1 & 0 \\ 1 & 0 & 0 \\ 0 & 0 & 1 \end{pmatrix}

(3) M=A1A2=(10000−1010)(0−10100001)=(0−1000−1100)M = A_1A_2 = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & -1 \\ 0 & 1 & 0 \end{pmatrix} \begin{pmatrix} 0 & -1 & 0 \\ 1 & 0 & 0 \\ 0 & 0 & 1 \end{pmatrix} = \begin{pmatrix} 0 & -1 & 0 \\ 0 & 0 & -1 \\ 1 & 0 & 0 \end{pmatrix} Eigenvalues: det⁡(M−λI)=0\det(M - \lambda I) = 0

∣−λ−100−λ−110−λ∣=−λ(λ2)+(−1)(−1)=−λ3+1=0\begin{vmatrix} -\lambda & -1 & 0 \\ 0 & -\lambda & -1 \\ 1 & 0 & -\lambda \end{vmatrix} = -\lambda(\lambda^2) + (-1)(-1) = -\lambda^3 + 1 = 0
λ3=1\lambda^3 = 1
λ=1,ω,ω2\lambda = 1, \omega, \omega^2

where ω=e2πi/3=−12+i32\omega = e^{2\pi i / 3} = -\frac{1}{2} + i \frac{\sqrt{3}}{2} M is a rotation matrix. Tr(M) = 0. The rotation angle θ\theta satisfies 1+2cos⁡θ=01 + 2\cos\theta = 0 , so cos⁡θ=−12\cos\theta = -\frac{1}{2} , and θ=2π3\theta = \frac{2\pi}{3} . To find the rotation axis, we need to find the eigenvector corresponding to λ=1\lambda = 1 . (M−I)v=0(M - I)v = 0 , where v=(x,y,z)v = (x,y,z) .

(−1−100−1−110−1)(xyz)=(000)\begin{pmatrix} -1 & -1 & 0 \\ 0 & -1 & -1 \\ 1 & 0 & -1 \end{pmatrix} \begin{pmatrix} x \\ y \\ z \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}

−x−y=0,−y−z=0,x−z=0-x - y = 0, -y - z = 0, x - z = 0 x=z,y=−xx = z, y = -x So the eigenvector is (x,−x,x)=x(1,−1,1)(x, -x, x) = x(1, -1, 1) . So the rotation axis is (1,−1,1)(1, -1, 1) .