東京大学 新領域創成科学研究科 メディカル情報生命専攻 2016年8月実施 問題8
Author
zephyr
Description
Answer the following questions about linear algebra.
(1) Denote by o \mathbf{o} o the zero vector. Let a \mathbf{a} a denote a two-dimensional vector that is not o \mathbf{o} o . T a ( x ) T_{\mathbf{a}}(\mathbf{x}) T a ( x ) is the orthogonal projection of a point x \mathbf{x} x on a \mathbf{a} a . Prove the following propositions.
(1.1) T a ( T a ( x ) ) = T a ( x ) T_{\mathbf{a}}(T_{\mathbf{a}}(\mathbf{x})) = T_{\mathbf{a}}(\mathbf{x}) T a ( T a ( x )) = T a ( x ) for any two-dimensional point x \mathbf{x} x .
(1.2) T b ( T a ( x ) ) = o T_{\mathbf{b}}(T_{\mathbf{a}}(\mathbf{x})) = \mathbf{o} T b ( T a ( x )) = o for any non-zero two-dimensional vector b \mathbf{b} b orthogonal to a \mathbf{a} a .
(2) Assume that a real symmetric matrix P \mathbf{P} P satisfies P 2 = P \mathbf{P}^2 = \mathbf{P} P 2 = P . Prove that the eigenvalues of P \mathbf{P} P are either 0 or 1.
(3) Denote by a 1 , a 2 \mathbf{a_1}, \mathbf{a_2} a 1 , 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 = [ a 1 , a 2 ] \mathbf{A} = [\mathbf{a_1}, \mathbf{a_2}] A = [ a 1 , a 2 ] .
回答以下有关线性代数的问题。
(1) 用 o \mathbf{o} o 表示零向量。设 a \mathbf{a} a 表示一个二维向量,它不是 o \mathbf{o} o 。T a ( x ) T_{\mathbf{a}}(\mathbf{x}) T a ( x ) 是点 x \mathbf{x} x 在 a \mathbf{a} a 上的正交投影。证明以下命题。
(1.1) 对于任意二维点 x \mathbf{x} x ,T a ( T a ( x ) ) = T a ( x ) T_{\mathbf{a}}(T_{\mathbf{a}}(\mathbf{x})) = T_{\mathbf{a}}(\mathbf{x}) T a ( T a ( x )) = T a ( x ) 。
(1.2) 对于任意非零二维向量 b \mathbf{b} b ,它与 a \mathbf{a} a 正交,T b ( T a ( x ) ) = o T_{\mathbf{b}}(T_{\mathbf{a}}(\mathbf{x})) = \mathbf{o} T b ( T a ( x )) = o 。
(2) 假设一个实对称矩阵 P \mathbf{P} P 满足 P 2 = P \mathbf{P}^2 = \mathbf{P} P 2 = P 。证明 P \mathbf{P} P 的特征值要么是 0,要么是 1。
(3) 用 a 1 , a 2 \mathbf{a_1}, \mathbf{a_2} a 1 , a 2 表示对应于三维空间的二维子空间基的列向量。用 A = [ a 1 , a 2 ] \mathbf{A} = [\mathbf{a_1}, \mathbf{a_2}] A = [ a 1 , a 2 ] 描述该子空间的投影矩阵。
题目描述
回答下列线性代数问题:
记零向量为 o \mathbf o o ,取非零二维向量 a \mathbf a a ,并令 T a ( x ) T_{\mathbf a}(\mathbf x) T a ( x ) 表示点 x \mathbf x x 在 a \mathbf a a 所张成直线上的正交投影。结合上图证明:
对任意二维点 x \mathbf x x ,
T a ( T a ( x ) ) = T a ( x ) ; T_{\mathbf a}(T_{\mathbf a}(\mathbf x))=T_{\mathbf a}(\mathbf x); T a ( T a ( x )) = T a ( x ) ;
对任意与 a \mathbf a a 正交的非零二维向量 b \mathbf b b ,
T b ( T a ( x ) ) = o . T_{\mathbf b}(T_{\mathbf a}(\mathbf x))=\mathbf o. T b ( T a ( x )) = o .
设实对称矩阵 P \mathbf P P 满足 P 2 = P \mathbf P^2=\mathbf P P 2 = P ,证明其特征值只能是 0 0 0 或 1 1 1 。
在三维空间中,列向量 a 1 , a 2 \mathbf a_1,\mathbf a_2 a 1 , a 2 构成某二维子空间的一组基。令
A = [ a 1 , a 2 ] , \mathbf A=[\mathbf a_1,\mathbf a_2], A = [ a 1 , a 2 ] ,
用 A \mathbf A A 写出到该子空间的正交投影矩阵。
正交投影算子 :用投影的几何意义或内积公式证明幂等性,以及正交方向上的再次投影为零。
投影矩阵的谱性质 :把 P 2 = P \mathbf P^2=\mathbf P P 2 = P 作用于特征向量,推出特征值满足 λ 2 = λ \lambda^2=\lambda λ 2 = λ 。
非正交基下的投影矩阵 :由基矩阵 A \mathbf A A 与 Gram 矩阵构造到其列空间的正交投影。
Kai
(1)
(1.1)
The orthogonal projection of x \mathbf{x} x on a \mathbf{a} a is given by:
T a ( x ) = a ⋅ x a ⋅ a a T_{\mathbf{a}}(\mathbf{x}) = \frac{\mathbf{a} \cdot \mathbf{x}}{\mathbf{a} \cdot \mathbf{a}} \mathbf{a} T a ( x ) = a ⋅ a a ⋅ x a
To prove the proposition, we apply T a T_{\mathbf{a}} T a again on T a ( x ) T_{\mathbf{a}}(\mathbf{x}) T a ( x ) :
T a ( T a ( x ) ) = T a ( a ⋅ x a ⋅ a a ) 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) T a ( T a ( x )) = T a ( a ⋅ a a ⋅ x a )
Using the definition of orthogonal projection:
T a ( a ⋅ x a ⋅ a a ) = a ⋅ ( a ⋅ x a ⋅ a a ) a ⋅ a a T_{\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} T a ( a ⋅ a a ⋅ x a ) = a ⋅ a a ⋅ ( a ⋅ a a ⋅ x a ) a
Simplifying the dot products:
= ( a ⋅ x ) ( a ⋅ a ) ( a ⋅ a ) a ⋅ a a = a ⋅ x a ⋅ a a = \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} = a ⋅ a ( a ⋅ a ) ( a ⋅ x ) ( a ⋅ a ) a = a ⋅ a a ⋅ x a
Thus:
T a ( T a ( x ) ) = T a ( x ) T_{\mathbf{a}}(T_{\mathbf{a}}(\mathbf{x})) = T_{\mathbf{a}}(\mathbf{x}) T a ( T a ( x )) = T a ( x )
(1.2)
Given b ⋅ a = 0 \mathbf{b} \cdot \mathbf{a} = 0 b ⋅ a = 0 , we need to show:
T b ( T a ( x ) ) = T b ( a ⋅ x a ⋅ a a ) 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) T b ( T a ( x )) = T b ( a ⋅ a a ⋅ x a )
Using the definition of orthogonal projection:
T b ( a ⋅ x a ⋅ a a ) = b ⋅ ( a ⋅ x a ⋅ a a ) b ⋅ b b T_{\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} T b ( a ⋅ a a ⋅ x a ) = b ⋅ b b ⋅ ( a ⋅ a a ⋅ x a ) b
Since b ⋅ a = 0 \mathbf{b} \cdot \mathbf{a} = 0 b ⋅ a = 0 :
= ( a ⋅ x ) ( b ⋅ a ) ( a ⋅ a ) b ⋅ b b = ( a ⋅ x ) ⋅ 0 ( a ⋅ a ) ⋅ ( b ⋅ b ) 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 = b ⋅ b ( a ⋅ a ) ( a ⋅ x ) ( b ⋅ a ) b = ( a ⋅ a ) ⋅ ( b ⋅ b ) ( a ⋅ x ) ⋅ 0 b = 0
Thus:
T b ( T a ( x ) ) = o T_{\mathbf{b}}(T_{\mathbf{a}}(\mathbf{x})) = \mathbf{o} T b ( T a ( x )) = o
(2)
Assume that a real symmetric matrix P \mathbf{P} P satisfies P 2 = P \mathbf{P}^2 = \mathbf{P} P 2 = P . Prove that the eigenvalues of P \mathbf{P} P are either 0 or 1.
Let P \mathbf{P} P be a real symmetric matrix. Therefore, it is diagonalizable. Let v \mathbf{v} v be an eigenvector of P \mathbf{P} P with eigenvalue λ \lambda λ :
P v = λ v \mathbf{P}\mathbf{v} = \lambda \mathbf{v} Pv = λ v
Applying P \mathbf{P} P again:
P 2 v = P ( P v ) = P ( λ v ) = λ P v = λ ( λ v ) = λ 2 v \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} P 2 v = P ( Pv ) = P ( λ v ) = λ Pv = λ ( λ v ) = λ 2 v
Since P 2 = P \mathbf{P}^2 = \mathbf{P} P 2 = P , we have:
P 2 v = P v = λ v \mathbf{P}^2 \mathbf{v} = \mathbf{P} \mathbf{v} = \lambda \mathbf{v} P 2 v = Pv = λ v
Thus:
λ 2 v = λ v \lambda^2 \mathbf{v} = \lambda \mathbf{v} λ 2 v = λ v
Since v \mathbf{v} v is a non-zero vector, we can conclude:
λ 2 = λ \lambda^2 = \lambda λ 2 = λ
Thus, the eigenvalues λ \lambda λ must satisfy:
λ ( λ − 1 ) = 0 \lambda (\lambda - 1) = 0 λ ( λ − 1 ) = 0
Therefore:
λ = 0 or λ = 1 \lambda = 0 \quad \text{or} \quad \lambda = 1 λ = 0 or λ = 1
(3)
Given the matrix A \mathbf{A} A formed by two column vectors a 1 \mathbf{a_1} a 1 and a 2 \mathbf{a_2} 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} P that projects any vector in R 3 \mathbb{R}^3 R 3 onto this subspace.
Matrix A \mathbf{A} A is:
A = [ a 1 , a 2 ] \mathbf{A} = [\mathbf{a_1}, \mathbf{a_2}] A = [ a 1 , a 2 ]
where A \mathbf{A} A is a 3 × 2 3 \times 2 3 × 2 matrix.
Derivation of the Projection Matrix
1. Projection of a Vector
The projection of a vector x \mathbf{x} x onto the subspace spanned by the columns of A \mathbf{A} A can be expressed as a linear combination of the columns of A \mathbf{A} A :
P x = c 1 a 1 + c 2 a 2 \mathbf{P}\mathbf{x} = c_1 \mathbf{a_1} + c_2 \mathbf{a_2} Px = c 1 a 1 + c 2 a 2
In matrix form, we write:
P x = A c \mathbf{P}\mathbf{x} = \mathbf{A}\mathbf{c} Px = Ac
where c \mathbf{c} c is a column vector of coefficients:
c = [ c 1 c 2 ] \mathbf{c} = \begin{bmatrix}
c_1 \\
c_2
\end{bmatrix} c = [ c 1 c 2 ]
2. Finding the Coefficients
To determine the coefficients c \mathbf{c} c , we use the property that the projection minimizes the distance to the subspace. This can be formulated as:
A T ( x − A c ) = 0 \mathbf{A}^T (\mathbf{x} - \mathbf{A}\mathbf{c}) = \mathbf{0} A T ( x − Ac ) = 0
This equation implies:
A T x = A T A c \mathbf{A}^T \mathbf{x} = \mathbf{A}^T \mathbf{A} \mathbf{c} A T x = A T Ac
Assuming A T A \mathbf{A}^T \mathbf{A} A T A is invertible, we solve for c \mathbf{c} c :
c = ( A T A ) − 1 A T x \mathbf{c} = (\mathbf{A}^T \mathbf{A})^{-1} \mathbf{A}^T \mathbf{x} c = ( A T A ) − 1 A T x
3. Constructing the Projection Matrix
Substituting c \mathbf{c} c back into the projection formula, we have:
P x = A ( A T A ) − 1 A T x \mathbf{P}\mathbf{x} = \mathbf{A} (\mathbf{A}^T \mathbf{A})^{-1} \mathbf{A}^T \mathbf{x} Px = A ( A T A ) − 1 A T x
Since this holds for any vector x \mathbf{x} x , the projection matrix P \mathbf{P} P can be identified as:
P = A ( A T A ) − 1 A T \mathbf{P} = \mathbf{A} (\mathbf{A}^T \mathbf{A})^{-1} \mathbf{A}^T P = A ( A T A ) − 1 A T
Knowledge
对称矩阵 特征值和特征向量 投影矩阵
重点词汇
Orthogonal projection 正交投影
Symmetric matrix 对称矩阵
Eigenvalue 特征值
Column vector 列向量
Subspace 子空间
Projection matrix 投影矩阵
参考资料
Gilbert Strang, "Linear Algebra and Its Applications," Chap. 3, 5.
David C. Lay, "Linear Algebra and Its Applications," Chap. 6.