跳到主要内容

京都大学 情報学研究科 数理工学専攻 2012年8月実施 基礎数学 II

Author​

思齐塾, 祭音Myyura

Description​

C n×m\mathbb{C}^{\,n \times m} で、 n×mn \times m 複素行列の全体を表す。 行列 A∈C n×mA \in \mathbb{C}^{\,n \times m} のエルミート共役を A∗A^* で表す。 H∈C n×nH \in \mathbb{C}^{\,n \times n} を正定値エルミート行列とする。 以下の問いに答えよ。 (i) HH は逆行列をもつが, H−1H^{-1} もまた正定値エルミート行列であることを示せ.

(ii) nn の分割 n=r+s,r>0,s>0n = r + s, r > 0, s > 0 , に応じて, HH を次のようにブロック行列に分割する.

H=(ABB∗D)(∗)H = \begin{pmatrix} A & B \\ B^* & D \end{pmatrix} \qquad (*)

ただし, A∈Cr×r,B∈Cr×s,B∗∈Cs×r,D∈Cs×sA \in C^{r \times r}, B \in C^{r \times s}, B^* \in C^{s \times r}, D \in C^{s \times s} . このとき, AA も DD も正定値エルミート行列であることを示せ.

(iii) 適当な下三角行列を左から, 上三角行列を右からそれぞれかけることにより, (*) の行列は, AA を変化させることなくブロック対角行列 (左下と右上のブロックが零行列の形の行列) に変形できることを示せ.

(iv) 次に, H∈Cn×nH \in C^{n \times n} の (i,j)(i, j) 成分を hijh_{ij} で表すとき,

det⁡H≤h11h22⋯hnn\det H \leq h_{11} h_{22} \cdots h_{nn}

の成り立つことを nn に関する数学的帰納法で証明せよ.

(v) A∈Cn×nA \in C^{n \times n} を A=(a1,…,an)A = (\mathbf{a}_1, \dots, \mathbf{a}_n) の形に表す. ただし, ak∈Cn,k=1,…,n\mathbf{a}_k \in C^n, k = 1, \dots, n . このとき

∣det⁡(A)∣2≤∣∣a1∣∣2⋯∣∣an∣∣2|\det(A)|^2 \leq ||\mathbf{a}_1||^2 \cdots ||\mathbf{a}_n||^2

の成り立つことを示せ. ただし, ∣∣ak∣∣||\mathbf{a}_k|| は ak\mathbf{a}_k の標準的ノルムを表す.

题目描述​

以 Cn×m\mathbb C^{n\times m} 表示全体 n×mn\times m 复矩阵,以 A∗A^* 表示矩阵 AA 的 Hermite 共轭。设 H∈Cn×nH\in\mathbb C^{n\times n} 是正定 Hermite 矩阵。完成以下各问:

  1. 已知 HH 可逆,证明 H−1H^{-1} 也是正定 Hermite 矩阵。

  2. 对分拆 n=r+sn=r+s(r>0, s>0r>0,\ s>0),将 HH 写成

    H=(ABB∗D),(*)H= \begin{pmatrix} A&B\\ B^*&D \end{pmatrix},\tag{*}

    其中

    A∈Cr×r,B∈Cr×s,B∗∈Cs×r,D∈Cs×s.A\in\mathbb C^{r\times r},\quad B\in\mathbb C^{r\times s},\quad B^*\in\mathbb C^{s\times r},\quad D\in\mathbb C^{s\times s}.

    证明 AA 与 DD 都是正定 Hermite 矩阵。

  3. 证明可以在 (∗)(*) 左侧乘某个适当的下三角分块矩阵、右侧乘某个适当的上三角分块矩阵,在保持左上角分块 AA 不变的同时,把 HH 化为上下两个非对角分块均为零的分块对角矩阵。

  4. 记 HH 的 (i,j)(i,j) 元为 hijh_{ij}。用关于 nn 的数学归纳法证明

    det⁡H≤h11h22⋯hnn.\det H\leq h_{11}h_{22}\cdots h_{nn}.
  5. 将任意 A∈Cn×nA\in\mathbb C^{n\times n} 按列写成

    A=(a1,…,an),ak∈Cn.A=(\boldsymbol a_1,\ldots,\boldsymbol a_n), \qquad \boldsymbol a_k\in\mathbb C^n.

    证明

    ∣det⁡A∣2≤∥a1∥2⋯∥an∥2,|\det A|^2 \leq\|\boldsymbol a_1\|^2\cdots\|\boldsymbol a_n\|^2,

    其中 ∥ak∥\|\boldsymbol a_k\| 为标准范数。

Kai​

(i) H−1H^{-1} の正定値性​

HH は正定値なので ker⁡H={0}\ker H=\{0\} であり、正則である。また

(H−1)∗=(H∗)−1=H−1(H^{-1})^*=(H^*)^{-1}=H^{-1}

だから H−1H^{-1} は Hermite 行列である。任意の y≠0y\neq0 に対し x=H−1yx=H^{-1}y とおけば x≠0x\neq0 であり、

y∗H−1y=(Hx)∗H−1(Hx)=x∗Hx>0.y^*H^{-1}y=(Hx)^*H^{-1}(Hx)=x^*Hx>0.

よって H−1H^{-1} も正定値 Hermite 行列である。

(ii) 主対角ブロック​

u∈Cr∖{0}u\in\mathbb C^r\setminus\{0\} に対して (u0)≠0\binom u0\neq0 だから、

(u∗0)H(u0)=u∗Au>0.\begin{pmatrix}u^*&0\end{pmatrix} H \begin{pmatrix}u\\0\end{pmatrix} =u^*Au>0.

同様に、 v∈Cs∖{0}v\in\mathbb C^s\setminus\{0\} に対して

(0v∗)H(0v)=v∗Dv>0.\begin{pmatrix}0&v^*\end{pmatrix} H \begin{pmatrix}0\\v\end{pmatrix} =v^*Dv>0.

H=H∗H=H^* から A=A∗A=A^* 、 D=D∗D=D^* でもあるので、 A,DA,D はともに正定値 Hermite 行列である。

(iii) Schur 補行列​

(ii) より AA は正則である。次の下三角行列と上三角行列を用いる。

L=(Ir0−B∗A−1Is),U=(Ir−A−1B0Is).L= \begin{pmatrix} I_r&0\\ -B^*A^{-1}&I_s \end{pmatrix}, \qquad U= \begin{pmatrix} I_r&-A^{-1}B\\ 0&I_s \end{pmatrix}.

直接乗法すると、

LH=(AB0D−B∗A−1B)LH= \begin{pmatrix} A&B\\ 0&D-B^*A^{-1}B \end{pmatrix}

および

LHU=(A00D−B∗A−1B)LHU= \begin{pmatrix} A&0\\ 0&D-B^*A^{-1}B \end{pmatrix}

を得る。なお U=L∗U=L^* なので、これは合同変換でもある。したがって Schur 補行列

D−B∗A−1BD-B^*A^{-1}B

も正定値である。

(iv) Hadamard の不等式​

n=1n=1 の場合は等号である。 n−1n-1 次まで成立すると仮定し、 r=1r=1 、 s=n−1s=n-1 として

H=(h11BB∗D)H= \begin{pmatrix} h_{11}&B\\ B^*&D \end{pmatrix}

と分割する。(iii) から

det⁡H=h11det⁡(D−B∗h11−1B).\det H =h_{11}\det\left(D-B^*h_{11}^{-1}B\right).

S=B∗h11−1BS=B^*h_{11}^{-1}B とおくと SS は半正定値であり、 D−SD-S は正定値である。 D−1/2SD−1/2D^{-1/2}SD^{-1/2} の固有値を μi\mu_i とすれば 0≤μi<10\leq\mu_i<1 だから、

det⁡(D−S)=det⁡D∏i=1n−1(1−μi)≤det⁡D.\det(D-S) =\det D\prod_{i=1}^{n-1}(1-\mu_i) \leq\det D.

帰納法の仮定を DD に適用して、

det⁡H≤h11det⁡D≤h11h22⋯hnn\det H \leq h_{11}\det D \leq h_{11}h_{22}\cdots h_{nn}

を得る。

(v) 一般の行列に対する不等式​

AA が正則でなければ左辺は 00 なので明らかである。 AA が正則なら H=A∗AH=A^*A は正定値 Hermite 行列であり、

hii=ai∗ai=∥ai∥2h_{ii}=\boldsymbol a_i^*\boldsymbol a_i=\|\boldsymbol a_i\|^2

かつ

det⁡H=det⁡(A∗)det⁡A=∣det⁡A∣2.\det H=\det(A^*)\det A=|\det A|^2.

(iv) を適用すれば、

∣det⁡A∣2≤∥a1∥2⋯∥an∥2|\det A|^2 \leq\|\boldsymbol a_1\|^2\cdots\|\boldsymbol a_n\|^2

となる。