跳到主要内容

明治大学 先端数理科学研究科 現象数理学専攻 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_\alphaR3\mathbb{R}^3 の線形変換であることを示し、これを3次の正方行列で表せ。

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

(3) 行列の積 A1A2A_1 A_2MM とおく。 MM の固有値をすべて求めよ。また、行列 MMR3\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_\alphaR3\mathbb{R}^3 上的线性变换,并用三阶方阵表示该变换。

(2) 分别用矩阵 A1,A2A_1,A_2 表示绕 xx 轴旋转 π/2\pi/2 与绕 zz 轴旋转 π/2\pi/2。求 A1A_1A2A_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))=(100001010)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)=(010100001)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=(100001010)(010100001)=(010001100)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 . (MI)v=0(M - I)v = 0 , where v=(x,y,z)v = (x,y,z) .

(110011101)(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}

xy=0,yz=0,xz=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) .