九州大学 システム情報科学府 共通 2014年8月実施 线性代数
Author
思齐塾 , 祭音Myyura
Description
実ベクトル空間 R N \mathbb{R}^N R N にベクトル x 1 , x 2 , . . . , x k x_1, x_2, ..., x_k x 1 , x 2 , ... , x k が与えられているとする。 A = ∑ i = 1 k x i x i T ∈ R N × N A = \sum_{i=1}^k x_i x_i^T \in \mathbb{R}^{N \times N} A = ∑ i = 1 k x i x i T ∈ R N × N とおく。ただし、 x i T x_i^T x i T は x i x_i x i の転置を表す。また、実ベクトル u ∈ R N u \in \mathbb{R}^N u ∈ R N に対し、
σ u = ∑ i = 1 k ( u T x i ) 2 \sigma_u = \sqrt{\sum_{i=1}^k (u^T x_i)^2} σ u = i = 1 ∑ k ( u T x i ) 2
と定義する。このとき、以下の各問に答えよ。
(1) 任意の実ベクトル u ∈ R N u \in \mathbb{R}^N u ∈ R N に対し、 σ u 2 = u T A u \sigma_u^2 = u^T A u σ u 2 = u T A u が成り立つことを示せ。
(2) A A A の固有値がすべて非負であることを示せ。
(3) ある実行列 B B B が存在して、 A = B T B A = B^T B A = B T B が成り立つことを示せ。
ヒント: A A A は実対称行列なので対角化可能。
(4) 任意の u , v ∈ R N u, v \in \mathbb{R}^N u , v ∈ R N に対し、 σ u + σ v ≥ σ u + v \sigma_u + \sigma_v \geq \sigma_{u+v} σ u + σ v ≥ σ u + v が成り立つことを示せ。
ヒント:問(3) の行列 B B B による変換 z ↦ z ~ = B z z \mapsto \tilde{z} = Bz z ↦ z ~ = B z ( z ∈ R N z \in \mathbb{R}^N z ∈ R N ) を用いよ。
题目描述
在实向量空间 R N \mathbb R^N R N 中给定向量 x 1 , x 2 , … , x k x_1,x_2,\ldots,x_k x 1 , x 2 , … , x k ,定义
A = ∑ i = 1 k x i x i T ∈ R N × N , A=\sum_{i=1}^k x_ix_i^T\in\mathbb R^{N\times N}, A = i = 1 ∑ k x i x i T ∈ R N × N ,
其中 x i T x_i^T x i T 表示 x i x_i x i 的转置。对任意实向量 u ∈ R N u\in\mathbb R^N u ∈ R N ,再定义
σ u = ∑ i = 1 k ( u T x i ) 2 . \sigma_u=\sqrt{\sum_{i=1}^k(u^Tx_i)^2}. σ u = i = 1 ∑ k ( u T x i ) 2 .
回答下列问题:
证明对任意 u ∈ R N u\in\mathbb R^N u ∈ R N ,
σ u 2 = u T A u . \sigma_u^2=u^TAu. σ u 2 = u T A u .
证明 A A A 的全部特征值均非负。
证明存在实矩阵 B B B 使
可使用提示:A A A 是可对角化的实对称矩阵。
证明对任意 u , v ∈ R N u,v\in\mathbb R^N u , v ∈ R N ,
σ u + σ v ≥ σ u + v . \sigma_u+\sigma_v\ge\sigma_{u+v}. σ u + σ v ≥ σ u + v .
可使用提示:利用第 3 问的矩阵 B B B 所定义的变换 z ↦ z ~ = B z z\mapsto\widetilde z=Bz z ↦ z = B z 。
Kai
(1) σ u 2 = ∑ i = 1 k ( u T x i ) 2 = ∑ i = 1 k ( u T x i ) ( x i T u ) = ∑ i = 1 k u T x i x i T u = u T ( ∑ i = 1 k x i x i T ) u = u T A u \sigma_u^2 = \sum_{i=1}^k (u^T x_i)^2 = \sum_{i=1}^k (u^T x_i)(x_i^T u) = \sum_{i=1}^k u^T x_i x_i^T u = u^T (\sum_{i=1}^k x_i x_i^T) u = u^T A u σ u 2 = ∑ i = 1 k ( u T x i ) 2 = ∑ i = 1 k ( u T x i ) ( x i T u ) = ∑ i = 1 k u T x i x i T u = u T ( ∑ i = 1 k x i x i T ) u = u T A u
(2) A A A は実対称行列なので、直交行列 P P P により A = P D P − 1 = P D P T A = PDP^{-1} = PDP^T A = P D P − 1 = P D P T と対角化可能。ここで D D D は対角行列であり、対角成分は A A A の固有値である。任意のベクトル x x x に対して、 x T A x = x T ( ∑ i = 1 k x i x i T ) x = ∑ i = 1 k x T x i x i T x = ∑ i = 1 k ( x T x i ) 2 ≥ 0 x^T A x = x^T (\sum_{i=1}^k x_i x_i^T) x = \sum_{i=1}^k x^T x_i x_i^T x = \sum_{i=1}^k (x^T x_i)^2 \geq 0 x T A x = x T ( ∑ i = 1 k x i x i T ) x = ∑ i = 1 k x T x i x i T x = ∑ i = 1 k ( x T x i ) 2 ≥ 0 。したがって、 A A A は半正定値行列であり、固有値はすべて非負である。
(3) (2)より、 A A A は半正定値対称行列であり、固有値はすべて非負である。 A = P D P T A = PDP^T A = P D P T と対角化する。 D D D の対角成分は A A A の固有値 λ i \lambda_i λ i であり、 λ i ≥ 0 \lambda_i \geq 0 λ i ≥ 0 。そこで、 D 1 / 2 D^{1/2} D 1/2 を対角成分が λ i \sqrt{\lambda_i} λ i である対角行列とする。すると、 D = D 1 / 2 D 1 / 2 D = D^{1/2} D^{1/2} D = D 1/2 D 1/2 。 A = P D P T = P D 1 / 2 D 1 / 2 P T = ( D 1 / 2 P T ) T ( D 1 / 2 P T ) A = PDP^T = P D^{1/2} D^{1/2} P^T = (D^{1/2} P^T)^T (D^{1/2} P^T) A = P D P T = P D 1/2 D 1/2 P T = ( D 1/2 P T ) T ( D 1/2 P T ) 。 B = D 1 / 2 P T B = D^{1/2} P^T B = D 1/2 P T とおくと、 A = B T B A = B^T B A = B T B 。
(4) σ u = ∑ i = 1 k ( u T x i ) 2 = u T A u = u T B T B u = ( B u ) T ( B u ) = ∣ ∣ B u ∣ ∣ \sigma_u = \sqrt{\sum_{i=1}^k (u^T x_i)^2} = \sqrt{u^T A u} = \sqrt{u^T B^T B u} = \sqrt{(Bu)^T (Bu)} = ||Bu|| σ u = ∑ i = 1 k ( u T x i ) 2 = u T A u = u T B T B u = ( B u ) T ( B u ) = ∣∣ B u ∣∣ ここで、 ∣ ∣ ⋅ ∣ ∣ || \cdot || ∣∣ ⋅ ∣∣ はユークリッドノルムである。すると、 σ u + σ v = ∣ ∣ B u ∣ ∣ + ∣ ∣ B v ∣ ∣ ≥ ∣ ∣ B u + B v ∣ ∣ = ∣ ∣ B ( u + v ) ∣ ∣ = σ u + v \sigma_u + \sigma_v = ||Bu|| + ||Bv|| \geq ||Bu + Bv|| = ||B(u+v)|| = \sigma_{u+v} σ u + σ v = ∣∣ B u ∣∣ + ∣∣ B v ∣∣ ≥ ∣∣ B u + B v ∣∣ = ∣∣ B ( u + v ) ∣∣ = σ u + v 。これは三角不等式による。