跳到主要内容

九州大学 システム情報科学府 共通 2014年8月実施 线性代数

Author

思齐塾, 祭音Myyura

Description

実ベクトル空間 RN\mathbb{R}^N にベクトル x1,x2,...,xkx_1, x_2, ..., x_k が与えられているとする。 A=i=1kxixiTRN×NA = \sum_{i=1}^k x_i x_i^T \in \mathbb{R}^{N \times N} とおく。ただし、 xiTx_i^Txix_i の転置を表す。また、実ベクトル uRNu \in \mathbb{R}^N に対し、

σu=i=1k(uTxi)2\sigma_u = \sqrt{\sum_{i=1}^k (u^T x_i)^2}

と定義する。このとき、以下の各問に答えよ。

(1) 任意の実ベクトル uRNu \in \mathbb{R}^N に対し、 σu2=uTAu\sigma_u^2 = u^T A u が成り立つことを示せ。

(2) AA の固有値がすべて非負であることを示せ。

(3) ある実行列 BB が存在して、 A=BTBA = B^T B が成り立つことを示せ。 ヒント: AA は実対称行列なので対角化可能。

(4) 任意の u,vRNu, v \in \mathbb{R}^N に対し、 σu+σvσu+v\sigma_u + \sigma_v \geq \sigma_{u+v} が成り立つことを示せ。 ヒント:問(3) の行列 BB による変換 zz~=Bzz \mapsto \tilde{z} = Bz ( zRNz \in \mathbb{R}^N ) を用いよ。

题目描述

在实向量空间 RN\mathbb R^N 中给定向量 x1,x2,,xkx_1,x_2,\ldots,x_k,定义

A=i=1kxixiTRN×N,A=\sum_{i=1}^k x_ix_i^T\in\mathbb R^{N\times N},

其中 xiTx_i^T 表示 xix_i 的转置。对任意实向量 uRNu\in\mathbb R^N,再定义

σu=i=1k(uTxi)2.\sigma_u=\sqrt{\sum_{i=1}^k(u^Tx_i)^2}.

回答下列问题:

  1. 证明对任意 uRNu\in\mathbb R^N

    σu2=uTAu.\sigma_u^2=u^TAu.
  2. 证明 AA 的全部特征值均非负。

  3. 证明存在实矩阵 BB 使

    A=BTB.A=B^TB.

    可使用提示:AA 是可对角化的实对称矩阵。

  4. 证明对任意 u,vRNu,v\in\mathbb R^N

    σu+σvσu+v.\sigma_u+\sigma_v\ge\sigma_{u+v}.

    可使用提示:利用第 3 问的矩阵 BB 所定义的变换 zz~=Bzz\mapsto\widetilde z=Bz

Kai

(1) σu2=i=1k(uTxi)2=i=1k(uTxi)(xiTu)=i=1kuTxixiTu=uT(i=1kxixiT)u=uTAu\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

(2) AA は実対称行列なので、直交行列 PP により A=PDP1=PDPTA = PDP^{-1} = PDP^T と対角化可能。ここで DD は対角行列であり、対角成分は AA の固有値である。任意のベクトル xx に対して、 xTAx=xT(i=1kxixiT)x=i=1kxTxixiTx=i=1k(xTxi)20x^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 。したがって、 AA は半正定値行列であり、固有値はすべて非負である。

(3) (2)より、 AA は半正定値対称行列であり、固有値はすべて非負である。 A=PDPTA = PDP^T と対角化する。 DD の対角成分は AA の固有値 λi\lambda_i であり、 λi0\lambda_i \geq 0 。そこで、 D1/2D^{1/2} を対角成分が λi\sqrt{\lambda_i} である対角行列とする。すると、 D=D1/2D1/2D = D^{1/2} D^{1/2}A=PDPT=PD1/2D1/2PT=(D1/2PT)T(D1/2PT)A = PDP^T = P D^{1/2} D^{1/2} P^T = (D^{1/2} P^T)^T (D^{1/2} P^T)B=D1/2PTB = D^{1/2} P^T とおくと、 A=BTBA = B^T B

(4) σu=i=1k(uTxi)2=uTAu=uTBTBu=(Bu)T(Bu)=Bu\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|| ここで、 || \cdot || はユークリッドノルムである。すると、 σu+σv=Bu+BvBu+Bv=B(u+v)=σu+v\sigma_u + \sigma_v = ||Bu|| + ||Bv|| \geq ||Bu + Bv|| = ||B(u+v)|| = \sigma_{u+v} 。これは三角不等式による。