跳到主要内容

東京大学 新領域創成科学研究科 メディカル情報生命専攻 2016年8月実施 問題8

Author

zephyr

Description

Answer the following questions about linear algebra.

(1) Denote by o\mathbf{o} the zero vector. Let a\mathbf{a} denote a two-dimensional vector that is not o\mathbf{o}. Ta(x)T_{\mathbf{a}}(\mathbf{x}) is the orthogonal projection of a point x\mathbf{x} on a\mathbf{a}. Prove the following propositions.

  • (1.1) Ta(Ta(x))=Ta(x)T_{\mathbf{a}}(T_{\mathbf{a}}(\mathbf{x})) = T_{\mathbf{a}}(\mathbf{x}) for any two-dimensional point x\mathbf{x}.

  • (1.2) Tb(Ta(x))=oT_{\mathbf{b}}(T_{\mathbf{a}}(\mathbf{x})) = \mathbf{o} for any non-zero two-dimensional vector b\mathbf{b} orthogonal to a\mathbf{a}.

(2) Assume that a real symmetric matrix P\mathbf{P} satisfies P2=P\mathbf{P}^2 = \mathbf{P}. Prove that the eigenvalues of P\mathbf{P} are either 0 or 1.

(3) Denote by a1,a2\mathbf{a_1}, \mathbf{a_2} the column vectors corresponding to the bases of a two-dimensional subspace of the three dimensional space. Describe the projection matrix to the subspace using A=[a1,a2]\mathbf{A} = [\mathbf{a_1}, \mathbf{a_2}].


回答以下有关线性代数的问题。

(1) 用 o\mathbf{o} 表示零向量。设 a\mathbf{a} 表示一个二维向量,它不是 o\mathbf{o}Ta(x)T_{\mathbf{a}}(\mathbf{x}) 是点 x\mathbf{x}a\mathbf{a} 上的正交投影。证明以下命题。

  • (1.1) 对于任意二维点 x\mathbf{x}Ta(Ta(x))=Ta(x)T_{\mathbf{a}}(T_{\mathbf{a}}(\mathbf{x})) = T_{\mathbf{a}}(\mathbf{x})

  • (1.2) 对于任意非零二维向量 b\mathbf{b},它与 a\mathbf{a} 正交,Tb(Ta(x))=oT_{\mathbf{b}}(T_{\mathbf{a}}(\mathbf{x})) = \mathbf{o}

(2) 假设一个实对称矩阵 P\mathbf{P} 满足 P2=P\mathbf{P}^2 = \mathbf{P}。证明 P\mathbf{P} 的特征值要么是 0,要么是 1。

(3) 用 a1,a2\mathbf{a_1}, \mathbf{a_2} 表示对应于三维空间的二维子空间基的列向量。用 A=[a1,a2]\mathbf{A} = [\mathbf{a_1}, \mathbf{a_2}] 描述该子空间的投影矩阵。

Kai

(1)

(1.1)

The orthogonal projection of x\mathbf{x} on a\mathbf{a} is given by:

Ta(x)=axaaaT_{\mathbf{a}}(\mathbf{x}) = \frac{\mathbf{a} \cdot \mathbf{x}}{\mathbf{a} \cdot \mathbf{a}} \mathbf{a}

To prove the proposition, we apply TaT_{\mathbf{a}} again on Ta(x)T_{\mathbf{a}}(\mathbf{x}):

Ta(Ta(x))=Ta(axaaa)T_{\mathbf{a}}(T_{\mathbf{a}}(\mathbf{x})) = T_{\mathbf{a}}\left( \frac{\mathbf{a} \cdot \mathbf{x}}{\mathbf{a} \cdot \mathbf{a}} \mathbf{a} \right)

Using the definition of orthogonal projection:

Ta(axaaa)=a(axaaa)aaaT_{\mathbf{a}}\left( \frac{\mathbf{a} \cdot \mathbf{x}}{\mathbf{a} \cdot \mathbf{a}} \mathbf{a} \right) = \frac{\mathbf{a} \cdot \left( \frac{\mathbf{a} \cdot \mathbf{x}}{\mathbf{a} \cdot \mathbf{a}} \mathbf{a} \right)}{\mathbf{a} \cdot \mathbf{a}} \mathbf{a}

Simplifying the dot products:

=(ax)(aa)(aa)aaa=axaaa= \frac{\frac{(\mathbf{a} \cdot \mathbf{x})(\mathbf{a} \cdot \mathbf{a})}{(\mathbf{a} \cdot \mathbf{a})}}{\mathbf{a} \cdot \mathbf{a}} \mathbf{a} = \frac{\mathbf{a} \cdot \mathbf{x}}{\mathbf{a} \cdot \mathbf{a}} \mathbf{a}

Thus:

Ta(Ta(x))=Ta(x)T_{\mathbf{a}}(T_{\mathbf{a}}(\mathbf{x})) = T_{\mathbf{a}}(\mathbf{x})

(1.2)

Given ba=0\mathbf{b} \cdot \mathbf{a} = 0, we need to show:

Tb(Ta(x))=Tb(axaaa)T_{\mathbf{b}}(T_{\mathbf{a}}(\mathbf{x})) = T_{\mathbf{b}}\left( \frac{\mathbf{a} \cdot \mathbf{x}}{\mathbf{a} \cdot \mathbf{a}} \mathbf{a} \right)

Using the definition of orthogonal projection:

Tb(axaaa)=b(axaaa)bbbT_{\mathbf{b}}\left( \frac{\mathbf{a} \cdot \mathbf{x}}{\mathbf{a} \cdot \mathbf{a}} \mathbf{a} \right) = \frac{\mathbf{b} \cdot \left( \frac{\mathbf{a} \cdot \mathbf{x}}{\mathbf{a} \cdot \mathbf{a}} \mathbf{a} \right)}{\mathbf{b} \cdot \mathbf{b}} \mathbf{b}

Since ba=0\mathbf{b} \cdot \mathbf{a} = 0:

=(ax)(ba)(aa)bbb=(ax)0(aa)(bb)b=0= \frac{\frac{(\mathbf{a} \cdot \mathbf{x})(\mathbf{b} \cdot \mathbf{a})}{(\mathbf{a} \cdot \mathbf{a})}}{\mathbf{b} \cdot \mathbf{b}} \mathbf{b} = \frac{(\mathbf{a} \cdot \mathbf{x}) \cdot 0}{(\mathbf{a} \cdot \mathbf{a}) \cdot (\mathbf{b} \cdot \mathbf{b})} \mathbf{b} = 0

Thus:

Tb(Ta(x))=oT_{\mathbf{b}}(T_{\mathbf{a}}(\mathbf{x})) = \mathbf{o}

(2)

Assume that a real symmetric matrix P\mathbf{P} satisfies P2=P\mathbf{P}^2 = \mathbf{P}. Prove that the eigenvalues of P\mathbf{P} are either 0 or 1.

Let P\mathbf{P} be a real symmetric matrix. Therefore, it is diagonalizable. Let v\mathbf{v} be an eigenvector of P\mathbf{P} with eigenvalue λ\lambda:

Pv=λv\mathbf{P}\mathbf{v} = \lambda \mathbf{v}

Applying P\mathbf{P} again:

P2v=P(Pv)=P(λv)=λPv=λ(λv)=λ2v\mathbf{P}^2 \mathbf{v} = \mathbf{P} (\mathbf{P} \mathbf{v}) = \mathbf{P} (\lambda \mathbf{v}) = \lambda \mathbf{P} \mathbf{v} = \lambda (\lambda \mathbf{v}) = \lambda^2 \mathbf{v}

Since P2=P\mathbf{P}^2 = \mathbf{P}, we have:

P2v=Pv=λv\mathbf{P}^2 \mathbf{v} = \mathbf{P} \mathbf{v} = \lambda \mathbf{v}

Thus:

λ2v=λv\lambda^2 \mathbf{v} = \lambda \mathbf{v}

Since v\mathbf{v} is a non-zero vector, we can conclude:

λ2=λ\lambda^2 = \lambda

Thus, the eigenvalues λ\lambda must satisfy:

λ(λ1)=0\lambda (\lambda - 1) = 0

Therefore:

λ=0orλ=1\lambda = 0 \quad \text{or} \quad \lambda = 1

(3)

Given the matrix A\mathbf{A} formed by two column vectors a1\mathbf{a_1} and a2\mathbf{a_2}, which represent the basis of a two-dimensional subspace in three-dimensional space, we want to find the projection matrix P\mathbf{P} that projects any vector in R3\mathbb{R}^3 onto this subspace.

Matrix A\mathbf{A} is:

A=[a1,a2]\mathbf{A} = [\mathbf{a_1}, \mathbf{a_2}]

where A\mathbf{A} is a 3×23 \times 2 matrix.

Derivation of the Projection Matrix

1. Projection of a Vector

The projection of a vector x\mathbf{x} onto the subspace spanned by the columns of A\mathbf{A} can be expressed as a linear combination of the columns of A\mathbf{A}:

Px=c1a1+c2a2\mathbf{P}\mathbf{x} = c_1 \mathbf{a_1} + c_2 \mathbf{a_2}

In matrix form, we write:

Px=Ac\mathbf{P}\mathbf{x} = \mathbf{A}\mathbf{c}

where c\mathbf{c} is a column vector of coefficients:

c=[c1c2]\mathbf{c} = \begin{bmatrix} c_1 \\ c_2 \end{bmatrix}
2. Finding the Coefficients

To determine the coefficients c\mathbf{c}, we use the property that the projection minimizes the distance to the subspace. This can be formulated as:

AT(xAc)=0\mathbf{A}^T (\mathbf{x} - \mathbf{A}\mathbf{c}) = \mathbf{0}

This equation implies:

ATx=ATAc\mathbf{A}^T \mathbf{x} = \mathbf{A}^T \mathbf{A} \mathbf{c}

Assuming ATA\mathbf{A}^T \mathbf{A} is invertible, we solve for c\mathbf{c}:

c=(ATA)1ATx\mathbf{c} = (\mathbf{A}^T \mathbf{A})^{-1} \mathbf{A}^T \mathbf{x}
3. Constructing the Projection Matrix

Substituting c\mathbf{c} back into the projection formula, we have:

Px=A(ATA)1ATx\mathbf{P}\mathbf{x} = \mathbf{A} (\mathbf{A}^T \mathbf{A})^{-1} \mathbf{A}^T \mathbf{x}

Since this holds for any vector x\mathbf{x}, the projection matrix P\mathbf{P} can be identified as:

P=A(ATA)1AT\mathbf{P} = \mathbf{A} (\mathbf{A}^T \mathbf{A})^{-1} \mathbf{A}^T

Knowledge

对称矩阵 特征值和特征向量 投影矩阵

重点词汇

  • Orthogonal projection 正交投影
  • Symmetric matrix 对称矩阵
  • Eigenvalue 特征值
  • Column vector 列向量
  • Subspace 子空间
  • Projection matrix 投影矩阵

参考资料

  1. Gilbert Strang, "Linear Algebra and Its Applications," Chap. 3, 5.
  2. David C. Lay, "Linear Algebra and Its Applications," Chap. 6.