跳到主要内容

東京大学 情報理工学系研究科 創造情報学専攻 2012年8月実施 筆記試験 第3問

Author

itsuitsuki

Description

In kinematic calculations and visual computations in robotics, vector operations are expressed with matrices sometimes. Answer the following questions on inner-product, outer-product, projection and rotation of three dimensional vectors. I\boldsymbol{I} is the 3×33\times3 identity matrix. The three-dimensional vectors x\boldsymbol{x}, a\boldsymbol{a}, b\boldsymbol{b} and n\boldsymbol{n} are 3×13\times1 column vectors: x=[xxxyxz],a=[axayaz],b=[bxbybz],n=[nxnynz].\boldsymbol{x}=\begin{bmatrix}x_x\\x_y\\x_z\end{bmatrix},\quad\boldsymbol{a}=\begin{bmatrix}a_x\\a_y\\a_z\end{bmatrix},\quad\boldsymbol{b}=\begin{bmatrix}b_x\\b_y\\b_z\end{bmatrix},\quad\boldsymbol{n}=\begin{bmatrix}n_x\\n_y\\n_z\end{bmatrix}. xT\boldsymbol{x}^T, which is a 1×31\times3 row vector, shows the transpose of x\boldsymbol{x}.

(1) On the inner product (a,b)(\boldsymbol{a},\boldsymbol{b}) between vectors a\boldsymbol{a} and b\boldsymbol{b}, describe i) the value of (a,b)(\boldsymbol{a},\boldsymbol{b}), ii) a 3×33\times3 matrix A\boldsymbol{A} which satisfies (a,b)a=Ab(\boldsymbol{a},\boldsymbol{b})\boldsymbol{a}=\boldsymbol{A}\boldsymbol{b}, and iii) A\boldsymbol{A} with vector a\boldsymbol{a} and its transpose aT\boldsymbol{a}^T.

(2) On the outer product a×b\boldsymbol{a}\times\boldsymbol{b} from a\boldsymbol{a} to b\boldsymbol{b}, describe i) 3×13\times1 expression of a×b\boldsymbol{a}\times\boldsymbol{b}, ii) 3×33\times3 expression of the matrix A\boldsymbol{A} which satisfies a×b=Ab\boldsymbol{a}\times\boldsymbol{b}=\boldsymbol{A}\boldsymbol{b}, iii) the 3×33\times3 matrix Q\boldsymbol{Q} which satisfies A=a×Q\boldsymbol{A}=\boldsymbol{a}\times\boldsymbol{Q} where x×Q\boldsymbol{x}\times\boldsymbol{Q} between a vector x\boldsymbol{x} and a 3×33\times3 matrix Q\boldsymbol{Q} means a 3×33\times3 matrix whose column vectors are three outer products from the vector x\boldsymbol{x} to each column vector in the matrix Q\boldsymbol{Q} respectively.

(3) As Figure 1 shows, a vector x\boldsymbol{x} is vertically projected to a vector y\boldsymbol{y} on a plane whose normal vector is a unit vector n\boldsymbol{n}. If the vector y\boldsymbol{y} is described as y=Px\boldsymbol{y}=\boldsymbol{P}\boldsymbol{x}, show that the matrix P\boldsymbol{P} becomes P=InnT\boldsymbol{P}=\boldsymbol{I}-\boldsymbol{n}\boldsymbol{n}^T.

(4) Three rotational matrices Rx(θx)\boldsymbol{R}_x(\theta_x), Ry(θy)\boldsymbol{R}_y(\theta_y) and Rz(θz)\boldsymbol{R}_z(\theta_z) are rotational matrices which rotate a vector x\boldsymbol{x} around the X-axis, Y-axis and Z-axis with θx\theta_x, θy\theta_y and θz\theta_z respectively, where the direction of the rotation for plus is clock-wise around the axis from the origin to infinity.

i) Describe 3×33\times3 expression of the matrix Rx(θx)\boldsymbol{R}_x(\theta_x), ii) As Figure 2 shows, the 3×33\times3 matrix Rn(θn)\boldsymbol{R}_n(\theta_n) is defined as the rotation matrix around a unit orientation vector n\boldsymbol{n} with θn\theta_n. Rn(θn)\boldsymbol{R}_n(\theta_n) is described as Rn(θn)=Rx(α)Ry(β)Rz(θn)Ry(β)Rx(α).\boldsymbol{R}_n(\theta_n)=\boldsymbol{R}_x(-\alpha)\boldsymbol{R}_y(\beta)\boldsymbol{R}_z(\theta_n)\boldsymbol{R}_y(-\beta)\boldsymbol{R}_x(\alpha). Explain what the variables α\alpha and β\beta become and explain why the expression is satisfied.

题目描述

机器人运动学和视觉计算常用矩阵表示向量运算。回答三维向量的内积、外积、投影与旋转问题。令 I\boldsymbol I3×33\times3 单位矩阵,并令

x=[xxxyxz],a=[axayaz],b=[bxbybz],n=[nxnynz]\boldsymbol{x}=\begin{bmatrix}x_x\\x_y\\x_z\end{bmatrix},\quad \boldsymbol{a}=\begin{bmatrix}a_x\\a_y\\a_z\end{bmatrix},\quad \boldsymbol{b}=\begin{bmatrix}b_x\\b_y\\b_z\end{bmatrix},\quad \boldsymbol{n}=\begin{bmatrix}n_x\\n_y\\n_z\end{bmatrix}

均为 3×13\times1 列向量;xT\boldsymbol x^T 是其 1×31\times3 转置行向量。

  1. a,b\boldsymbol a,\boldsymbol b 的内积 (a,b)(\boldsymbol a,\boldsymbol b),分别写出:a. 内积数值;b. 满足

    (a,b)a=Ab(\boldsymbol a,\boldsymbol b)\boldsymbol a=\boldsymbol A\boldsymbol b

    3×33\times3 矩阵 A\boldsymbol A;c. 用 a,aT\boldsymbol a,\boldsymbol a^T 表示该 A\boldsymbol A

  2. 对从 a\boldsymbol ab\boldsymbol b 的外积 a×b\boldsymbol a\times\boldsymbol b,分别写出:a. 其 3×13\times1 分量表达式;b. 满足 a×b=Ab\boldsymbol a\times\boldsymbol b=\boldsymbol A\boldsymbol b3×33\times3 矩阵 A\boldsymbol A;c. 满足 A=a×Q\boldsymbol A=\boldsymbol a\times\boldsymbol Q3×33\times3 矩阵 Q\boldsymbol Q。这里 x×Q\boldsymbol x\times\boldsymbol Q 定义为 3×33\times3 矩阵,其每一列分别是 x\boldsymbol xQ\boldsymbol Q 对应列向量的外积。

  3. 如图 1,把 x\boldsymbol x 正交投影到法向量为单位向量 n\boldsymbol n 的平面,得到 y\boldsymbol y。若 y=Px\boldsymbol y=\boldsymbol P\boldsymbol x,证明

    P=InnT.\boldsymbol P=\boldsymbol I-\boldsymbol n\boldsymbol n^T.
  4. Rx(θx),Ry(θy),Rz(θz)\boldsymbol R_x(\theta_x),\boldsymbol R_y(\theta_y),\boldsymbol R_z(\theta_z) 分别把向量绕 X,Y,ZX,Y,Z 轴旋转 θx,θy,θz\theta_x,\theta_y,\theta_z;从原点沿轴正方向看去,顺时针规定为正。

    1. 写出 Rx(θx)\boldsymbol R_x(\theta_x)3×33\times3 形式。

    2. 如图 2,Rn(θn)\boldsymbol R_n(\theta_n) 表示绕单位方向向量 n\boldsymbol n 旋转 θn\theta_n,且

      Rn(θn)=Rx(α)Ry(β)Rz(θn)Ry(β)Rx(α).\boldsymbol R_n(\theta_n)= \boldsymbol R_x(-\alpha)\boldsymbol R_y(\beta)\boldsymbol R_z(\theta_n) \boldsymbol R_y(-\beta)\boldsymbol R_x(\alpha).

      α,β\alpha,\beta,并解释该矩阵分解成立的原因。