跳到主要内容

京都大学 情報学研究科 システム科学専攻 2019年8月実施 数学 I

Author

思齐塾, 祭音Myyura

Description

日本語版

問1

R3\mathbb{R}^3 から R3\mathbb{R}^3 への線形写像 ff の表現行列が AA であるとして, 以下の設問に答えよ。 ただし, aa は実数とする。

A=(12101010a)A = \begin{pmatrix} 1 & 2 & -1 \\ 0 & 1 & 0 \\ 1 & 0 & a \end{pmatrix}

(i) 行列 AA の階数(ランク)が最小になる aa の値を求めよ。また, このときの行列 AA の階数を求めよ。

(ii) 行列 AA の階数が最小になるとき, 線形写像 ff の核(カーネル)を求めよ。また, このとき, ff の像の正規直交基底を求めよ。

(iii) 行列 AA が対角化できなくなる aa の値を求めよ。

(iv) a=3a=3 のとき, 行列 AA は行列 BB と相似であることを示せ。

B=(100021002)B = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 2 & 1 \\ 0 & 0 & 2 \end{pmatrix}

問2

以下の設問に答えよ。ただし, det\det は行列式, TT は転置, 行列の右上の 1-1 は逆行列を意味するものとし, 問題文中の行列およびベクトルの成分, スカラーはすべて実数とする。また, n,m,ln, m, l は正の整数とする。さらに, m×mm \times m 行列 SS , m×lm \times l 行列 TT , l×ml \times m 行列 UU , l×ll \times l 正則行列 VV について, (m+l)×(m+l)(m+l) \times (m+l) 行列 [STUV]\begin{bmatrix} S & T \\ U & V \end{bmatrix}STV1US-TV^{-1}U が正則であれば,

[STUV]1=[(STV1U)1(STV1U)1TV1V1U(STV1U)1V1+V1U(STV1U)1TV1]\begin{bmatrix} S & T \\ U & V \end{bmatrix}^{-1} = \begin{bmatrix} (S-TV^{-1}U)^{-1} & -(S-TV^{-1}U)^{-1}TV^{-1} \\ -V^{-1}U(S-TV^{-1}U)^{-1} & V^{-1} + V^{-1}U(S-TV^{-1}U)^{-1}TV^{-1} \end{bmatrix}

となる。

(i) 正則な n×nn \times n 行列 AA , nn 次元列ベクトル b,cb, c , スカラー dd について, 以下が成り立つことを示せ。

det[AbcTd]=(detA)×(dcTA1b)\det \begin{bmatrix} A & b \\ c^T & d \end{bmatrix} = (\det A) \times (d - c^T A^{-1} b)

(ii) n2n \geq 2 とする。 AAn×nn \times n 正定値対称行列とする。このとき, A1A^{-1}n×nn \times n 正定値対称行列となり, スカラー α>0\alpha > 0 , (n1)(n-1) 次元列ベクトル β\beta , (n1)×(n1)(n-1) \times (n-1) 行列 Δ\Delta を用いて,

A1=[αβTβΔ]A^{-1} = \begin{bmatrix} \alpha & \beta^T \\ \beta & \Delta \end{bmatrix}

と表すことができる。 A~\tilde{A} を行列 AA から最初の行と列を除いた (n1)×(n1)(n-1) \times (n-1) 小行列とするとき,

A~1=ΔββTα\tilde{A}^{-1} = \Delta - \frac{\beta\beta^T}{\alpha}

となることを示せ。

(iii) 設問 (ii) の条件に加えて, x=[x1,x2,,xn]Tx = [x_1, x_2, \dots, x_n]^T , また, xx(n1)(n-1) 次元部分ベクトルを x~=[x2,x3,,xn]T\tilde{x} = [x_2, x_3, \dots, x_n]^T とする。このとき, 二次形式 xTA1xx^T A^{-1} xx1x_1 について二次式となるが, その二次式の x1x_1 に関する最小値は x~TA~1x~\tilde{x}^T \tilde{A}^{-1} \tilde{x} となることを示せ。

(iv) 正方行列 AnA_n を以下のように定義する。

A1=[a0],A2=[a0a1a1a0],,An+1=[a0a1a2ana1a0a1an1a2a1a0an2anan1an2a0]A_1 = [a_0], \quad A_2 = \begin{bmatrix} a_0 & a_1 \\ a_1 & a_0 \end{bmatrix}, \quad \dots, \quad A_{n+1} = \begin{bmatrix} a_0 & a_1 & a_2 & \cdots & a_n \\ a_1 & a_0 & a_1 & \cdots & a_{n-1} \\ a_2 & a_1 & a_0 & \cdots & a_{n-2} \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ a_n & a_{n-1} & a_{n-2} & \cdots & a_0 \end{bmatrix}

いま, すべての nn について AnA_n を正定値対称行列とするとき, n2n \geq 2 について,

detAn+1(detAn)2detAn1\det A_{n+1} \leq \frac{(\det A_n)^2}{\det A_{n-1}}

となることを示せ。

题目描述

問1

设从 R3\mathbb{R}^3R3\mathbb{R}^3 的线性映射 ff 在标准基下的表示矩阵为

A=(12101010a),A= \begin{pmatrix} 1&2&-1\\ 0&1&0\\ 1&0&a \end{pmatrix},

其中 aa 为实数。回答:

  1. 求使 rankA\operatorname{rank}A 最小的 aa,并求此时的最小秩。
  2. AA 的秩最小时,求线性映射 ff 的核,并求 imf\operatorname{im}f 的一组标准正交基。
  3. 求所有使 AA 在实数域上不可对角化的 aa
  4. a=3a=3 时,证明 AA 与矩阵
B=(100021002)B= \begin{pmatrix} 1&0&0\\ 0&2&1\\ 0&0&2 \end{pmatrix}

相似。

問2

以下 det\det 表示行列式,记号 ATA^{\mathrm T} 表示矩阵 AA 的转置, A1A^{-1} 表示其逆矩阵;题中所有矩阵、向量和标量的分量均为实数, n,m,ln,m,l 均为正整数。

可以使用如下分块逆矩阵公式:设 SSm×mm\times m 矩阵,TTm×lm\times l 矩阵, UUl×ml\times m 矩阵,VV 为可逆的 l×ll\times l 矩阵。若

[STUV]STV1U\begin{bmatrix} S&T\\ U&V \end{bmatrix} \quad\text{和}\quad S-TV^{-1}U

均可逆,则

[STUV]1=[(STV1U)1(STV1U)1TV1V1U(STV1U)1V1+V1U(STV1U)1TV1].\begin{bmatrix} S&T\\ U&V \end{bmatrix}^{-1} = \begin{bmatrix} (S-TV^{-1}U)^{-1} & -(S-TV^{-1}U)^{-1}TV^{-1} \\[1mm] -V^{-1}U(S-TV^{-1}U)^{-1} & V^{-1} +V^{-1}U(S-TV^{-1}U)^{-1}TV^{-1} \end{bmatrix}.

回答下列问题。

  1. AA 是可逆的 n×nn\times n 矩阵, b,c\boldsymbol{b},\boldsymbol{c}nn 维列向量,dd 是标量。证明
det[AbcTd]=(detA)(dcTA1b).\det \begin{bmatrix} A&\boldsymbol{b}\\ \boldsymbol{c}^{\mathrm T}&d \end{bmatrix} = (\det A) \left( d-\boldsymbol{c}^{\mathrm T}A^{-1}\boldsymbol{b} \right).
  1. n2n\geq2AAn×nn\times n 正定实对称矩阵。于是 A1A^{-1} 也是正定实对称矩阵,并可写为
A1=[αβTβΔ],A^{-1} = \begin{bmatrix} \alpha&\boldsymbol{\beta}^{\mathrm T}\\ \boldsymbol{\beta}&\Delta \end{bmatrix},

其中 α>0\alpha>0 是标量, β\boldsymbol{\beta}(n1)(n-1) 维列向量, Δ\Delta(n1)×(n1)(n-1)\times(n-1) 矩阵。令 A~\widetilde{A} 为从 AA 删除第一行和第一列所得的 (n1)×(n1)(n-1)\times(n-1) 主子矩阵。证明

A~1=ΔββTα.\widetilde{A}^{-1} = \Delta -\frac{\boldsymbol{\beta}\boldsymbol{\beta}^{\mathrm T}}{\alpha}.
  1. 在上一小问的条件下,令
x=[x1,x2,,xn]T,x~=[x2,x3,,xn]T.\boldsymbol{x} = [x_1,x_2,\ldots,x_n]^{\mathrm T}, \qquad \widetilde{\boldsymbol{x}} = [x_2,x_3,\ldots,x_n]^{\mathrm T}.

把二次型 xTA1x\boldsymbol{x}^{\mathrm T}A^{-1}\boldsymbol{x} 视为关于 x1x_1 的二次函数,并保持 x~\widetilde{\boldsymbol{x}} 固定。证明它关于 x1x_1 的最小值为

x~TA~1x~.\widetilde{\boldsymbol{x}}^{\mathrm T} \widetilde{A}^{-1} \widetilde{\boldsymbol{x}}.
  1. 定义实对称 Toeplitz 方阵序列
A1=[a0],A2=[a0a1a1a0],A_1=[a_0], \qquad A_2= \begin{bmatrix} a_0&a_1\\ a_1&a_0 \end{bmatrix},

以及一般的

An+1=[a0a1a2ana1a0a1an1a2a1a0an2anan1an2a0].A_{n+1} = \begin{bmatrix} a_0&a_1&a_2&\cdots&a_n\\ a_1&a_0&a_1&\cdots&a_{n-1}\\ a_2&a_1&a_0&\cdots&a_{n-2}\\ \vdots&\vdots&\vdots&\ddots&\vdots\\ a_n&a_{n-1}&a_{n-2}&\cdots&a_0 \end{bmatrix}.

假设对每个 nnAnA_n 都是正定对称矩阵。证明对所有 n2n\geq2

detAn+1(detAn)2detAn1.\det A_{n+1} \leq \frac{(\det A_n)^2}{\det A_{n-1}}.

Kai

問1

(i) 解答 行列 AA の階数(ランク)は, AA の行列式 det(A)\det(A) が 0 でないときは 3, 0 のときは 3 未満となります。 まず, AA の行列式を計算します。

det(A)=12101010a=1(1a00)2(0a01)+(1)(0011)=a+1\det(A) = \begin{vmatrix} 1 & 2 & -1 \\ 0 & 1 & 0 \\ 1 & 0 & a \end{vmatrix} = 1(1 \cdot a - 0 \cdot 0) - 2(0 \cdot a - 0 \cdot 1) + (-1)(0 \cdot 0 - 1 \cdot 1) = a + 1

det(A)=0\det(A) = 0 となるのは a+1=0a+1=0 , すなわち a=1a=-1 のときです。このとき, 階数が最小になる可能性があります。

  • a1a \neq -1 のとき, det(A)0\det(A) \neq 0 なので, rank(A)=3\mathrm{rank}(A) = 3 です。
  • a=1a = -1 のとき, det(A)=0\det(A) = 0 なので, rank(A)<3\mathrm{rank}(A) < 3 です。 このときの行列 AA は,
A=(121010101)A = \begin{pmatrix} 1 & 2 & -1 \\ 0 & 1 & 0 \\ 1 & 0 & -1 \end{pmatrix}

左上の 2×22 \times 2 小行列式は 1201=10\begin{vmatrix} 1 & 2 \\ 0 & 1 \end{vmatrix} = 1 \neq 0 であるため, rank(A)2\mathrm{rank}(A) \geq 2 です。 したがって, a=1a=-1 のとき, rank(A)=2\mathrm{rank}(A)=2 となります。 これが最小の階数です。

答え: 階数が最小になる aa の値は a=1a=-1 で, そのときの階数は 2 です。

(ii) 解答 a=1a=-1 のとき, 行列 AA の階数は最小になります。

  1. 線形写像 ff の核(カーネル) Ker(f)\mathrm{Ker}(f)

Ker(f)\mathrm{Ker}(f) は, 方程式 Ax=0A\mathbf{x} = \mathbf{0} の解空間です。

(121010101)(xyz)=(000)\begin{pmatrix} 1 & 2 & -1 \\ 0 & 1 & 0 \\ 1 & 0 & -1 \end{pmatrix} \begin{pmatrix} x \\ y \\ z \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}

これを連立方程式として解くと,

{x+2yz=0y=0xz=0\begin{cases} x + 2y - z = 0 \\ y = 0 \\ x - z = 0 \end{cases}

第2式より y=0y=0 。これを第1式に代入すると xz=0x-z=0 となり, 第3式と同じです。 x=zx=z となります。 z=tz=t ( tt は任意の実数)とおくと, x=t,y=0,z=tx=t, y=0, z=t となります。 よって, 解ベクトルは x=t(101)\mathbf{x} = t \begin{pmatrix} 1 \\ 0 \\ 1 \end{pmatrix} と表せます。 したがって, Ker(f)\mathrm{Ker}(f) は以下のように表せます。

Ker(f)={c(101)cR}\mathrm{Ker}(f) = \left\{ c \begin{pmatrix} 1 \\ 0 \\ 1 \end{pmatrix} \mid c \in \mathbb{R} \right\}

2. 線形写像 ff の像 Im(f)\mathrm{Im}(f) の正規直交基底

Im(f)\mathrm{Im}(f)AA の列ベクトルで張られる空間(列空間)です。 rank(A)=2\mathrm{rank}(A)=2 なので, 基底は2つの線形独立な列ベクトルで構成されます。 AA の第1列と第2列は線形独立なので, 基底として {u1=(101),u2=(210)}\left\{ \mathbf{u}_1 = \begin{pmatrix} 1 \\ 0 \\ 1 \end{pmatrix}, \mathbf{u}_2 = \begin{pmatrix} 2 \\ 1 \\ 0 \end{pmatrix} \right\} を取ることができます。

この基底にグラム・シュミットの直交化法を適用して正規直交基底 {w1,w2}\left\{ \mathbf{w}_1, \mathbf{w}_2 \right\} を求めます。

  • w1\mathbf{w}_1 の計算:
w1=u1u1=112+02+12(101)=12(101)\mathbf{w}_1 = \frac{\mathbf{u}_1}{\|\mathbf{u}_1\|} = \frac{1}{\sqrt{1^2+0^2+1^2}} \begin{pmatrix} 1 \\ 0 \\ 1 \end{pmatrix} = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ 0 \\ 1 \end{pmatrix}
  • w2\mathbf{w}_2 の計算: まず, u2\mathbf{u}_2 から w1\mathbf{w}_1 方向の成分を引きます。
u2=u2(u2w1)w1=(210)((210)12(101))12(101)\mathbf{u}_2' = \mathbf{u}_2 - (\mathbf{u}_2 \cdot \mathbf{w}_1)\mathbf{w}_1 = \begin{pmatrix} 2 \\ 1 \\ 0 \end{pmatrix} - \left( \begin{pmatrix} 2 \\ 1 \\ 0 \end{pmatrix} \cdot \frac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ 0 \\ 1 \end{pmatrix} \right) \frac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ 0 \\ 1 \end{pmatrix}
u2=(210)2212(101)=(210)(101)=(111)\mathbf{u}_2' = \begin{pmatrix} 2 \\ 1 \\ 0 \end{pmatrix} - \frac{2}{\sqrt{2}} \cdot \frac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ 0 \\ 1 \end{pmatrix} = \begin{pmatrix} 2 \\ 1 \\ 0 \end{pmatrix} - \begin{pmatrix} 1 \\ 0 \\ 1 \end{pmatrix} = \begin{pmatrix} 1 \\ 1 \\ -1 \end{pmatrix}

次に, u2\mathbf{u}_2' を正規化します。

w2=u2u2=112+12+(1)2(111)=13(111)\mathbf{w}_2 = \frac{\mathbf{u}_2'}{\|\mathbf{u}_2'\|} = \frac{1}{\sqrt{1^2+1^2+(-1)^2}} \begin{pmatrix} 1 \\ 1 \\ -1 \end{pmatrix} = \frac{1}{\sqrt{3}} \begin{pmatrix} 1 \\ 1 \\ -1 \end{pmatrix}

答え: ff の像の正規直交基底は {12(101),13(111)}\left\{ \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ 0 \\ 1 \end{pmatrix}, \frac{1}{\sqrt{3}} \begin{pmatrix} 1 \\ 1 \\ -1 \end{pmatrix} \right\} です。

(iii) 解答 ここでは f:R3R3f:\mathbb R^3\to\mathbb R^3 の対角化なので,実数体上で考える。 特性多項式は

det(AλI)=(1λ)[λ2(a+1)λ+a+1]\det(A - \lambda I) = (1-\lambda) [\lambda^2 - (a+1)\lambda + a+1]

である。二次因子に λ=1\lambda=1 を代入した値は 11 なので, 11 が二次因子の根と重なることはない。二次因子の判別式は

D=(a+1)24(a+1)=(a+1)(a3)D = (a+1)^2 - 4(a+1) = (a+1)(a-3)

である。

  • a<1a<-1 または a>3a>3 なら D>0D>0 であり,相異なる3個の実固有値を もつので実対角化できる。
  • 1<a<3-1<a<3 なら D<0D<0 であり,非実共役な固有値をもつので 実数体上では対角化できない。
  • a=1a=-1 では固有値 00 の代数的多重度は2であるが,(i) より dimkerA=1\dim\ker A=1 なので対角化できない。
  • a=3a=3 では固有値 22 の代数的多重度は2である。一方,
rank(A2I)=2,dimker(A2I)=1,\operatorname{rank}(A-2I)=2,\qquad \dim\ker(A-2I)=1,

なので対角化できない。

したがって,実数体上で対角化できないための必要十分条件は

1a3.\boxed{-1\leq a\leq3}.

なお,複素数体上の対角化を意味する場合には, 1<a<3-1<a<3 の3固有値は 相異なるため対角化でき,対角化できない値は端点 a=1,3a=-1,3 だけである。

(iv) 解答 a=3a=3 のとき, 行列 AA が行列 BB と相似であることを示すには, AA のジョルダン標準形が BB と一致することを示せばよいです。

A=(121010103),B=(100021002)A = \begin{pmatrix} 1 & 2 & -1 \\ 0 & 1 & 0 \\ 1 & 0 & 3 \end{pmatrix}, \quad B = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 2 & 1 \\ 0 & 0 & 2 \end{pmatrix}

(iii)の計算より, a=3a=3 のときの AA の固有値とその多重度は以下の通りです。

  • 固有値 λ1=1\lambda_1 = 1 : 代数的多重度 1, 幾何学的多重度 1
  • 固有値 λ2=2\lambda_2 = 2 : 代数的多重度 2, 幾何学的多重度 1

ジョルダン標準形の構造は, これらの多重度によって決まります。

  • λ1=1\lambda_1 = 1 について: 幾何学的多重度が1なので, 1×11 \times 1 のジョルダンブロックが1つです: (1)\begin{pmatrix} 1 \end{pmatrix}
  • λ2=2\lambda_2 = 2 について: 幾何学的多重度が1なので, ジョルダンブロックは1つです。代数的多重度が2なので, ブロックのサイズは 2×22 \times 2 になります: (2102)\begin{pmatrix} 2 & 1 \\ 0 & 2 \end{pmatrix}

したがって, AA のジョルダン標準形 JAJ_A はこれらのブロックを対角に並べたものになります。

JA=(100021002)J_A = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 2 & 1 \\ 0 & 0 & 2 \end{pmatrix}

これは行列 BB と一致します。行列はそのジョルダン標準形と相似であるため, AAJAJ_A と相似であり, JA=BJ_A = B なので, AABB と相似です。 (証明終)

問2

解答

(i) の証明

与えられた行列を M=[AbcTd]M = \begin{bmatrix} A & b \\ c^T & d \end{bmatrix} とする。 AA は正則(可逆)であるため、 A1A^{-1} が存在する。 行列 MM を次のようにブロック行列の積で分解できる。

[AbcTd]=[A0cT1][IA1b0dcTA1b]\begin{bmatrix} A & b \\ c^T & d \end{bmatrix} = \begin{bmatrix} A & 0 \\ c^T & 1 \end{bmatrix} \begin{bmatrix} I & A^{-1}b \\ 0 & d - c^T A^{-1} b \end{bmatrix}

この分解が正しいことを確認する:

[A0cT1][IA1b0dcTA1b]=[AI+00A(A1b)+0(dcTA1b)cTI+10cT(A1b)+1(dcTA1b)]=[AbcTd]\begin{bmatrix} A & 0 \\ c^T & 1 \end{bmatrix} \begin{bmatrix} I & A^{-1}b \\ 0 & d - c^T A^{-1} b \end{bmatrix} = \begin{bmatrix} A \cdot I + 0 \cdot 0 & A(A^{-1}b) + 0 \cdot (d - c^T A^{-1} b) \\ c^T \cdot I + 1 \cdot 0 & c^T(A^{-1}b) + 1 \cdot (d - c^T A^{-1} b) \end{bmatrix} = \begin{bmatrix} A & b \\ c^T & d \end{bmatrix}

行列の積の行列式は、各行列の行列式の積に等しいので、

det(M)=det([A0cT1])×det([IA1b0dcTA1b])\det(M) = \det \begin{pmatrix} \begin{bmatrix} A & 0 \\ c^T & 1 \end{bmatrix} \end{pmatrix} \times \det \begin{pmatrix} \begin{bmatrix} I & A^{-1}b \\ 0 & d - c^T A^{-1} b \end{bmatrix} \end{pmatrix}

ブロック三角行列の行列式は、対角ブロックの行列式の積となる。

det[A0cT1]=(detA)×(det1)=detA\det \begin{bmatrix} A & 0 \\ c^T & 1 \end{bmatrix} = (\det A) \times (\det 1) = \det A
det[IA1b0dcTA1b]=(detI)×det(dcTA1b)=1×(dcTA1b)\det \begin{bmatrix} I & A^{-1}b \\ 0 & d - c^T A^{-1} b \end{bmatrix} = (\det I) \times \det(d - c^T A^{-1} b) = 1 \times (d - c^T A^{-1} b)

ここで、 dcTA1bd - c^T A^{-1} b はスカラー( 1×11 \times 1 行列)なので、その行列式は値自身である。 したがって、

det[AbcTd]=(detA)×(dcTA1b)\det \begin{bmatrix} A & b \\ c^T & d \end{bmatrix} = (\det A) \times (d - c^T A^{-1} b)

が示された。

(ii) の証明

AAn×nn \times n 正定値対称行列とする。 AA を次のように分割する。

A=[abTbA~]A = \begin{bmatrix} a & b^T \\ b & \tilde{A} \end{bmatrix}

ここで、 aa はスカラー、 bb(n1)(n-1) 次元の列ベクトル、 A~\tilde{A}(n1)×(n1)(n-1) \times (n-1) の小行列である。 AA が正定値であるため、その主小行列 A~\tilde{A} も正定値であり、したがって正則である。

問題の冒頭で与えられたブロック行列の逆行列の公式を用いる。ここで S=aS=a , T=bTT=b^T , U=bU=b , V=A~V=\tilde{A} と対応させる。

A1=[abTbA~]1=[(abTA~1b)1(abTA~1b)1bTA~1A~1b(abTA~1b)1A~1+A~1b(abTA~1b)1bTA~1]A^{-1} = \begin{bmatrix} a & b^T \\ b & \tilde{A} \end{bmatrix}^{-1} = \begin{bmatrix} (a-b^T\tilde{A}^{-1}b)^{-1} & -(a-b^T\tilde{A}^{-1}b)^{-1}b^T\tilde{A}^{-1} \\ -\tilde{A}^{-1}b(a-b^T\tilde{A}^{-1}b)^{-1} & \tilde{A}^{-1} + \tilde{A}^{-1}b(a-b^T\tilde{A}^{-1}b)^{-1}b^T\tilde{A}^{-1} \end{bmatrix}

与えられた A1A^{-1} の分割形式と比較する。

A1=[αβTβΔ]A^{-1} = \begin{bmatrix} \alpha & \beta^T \\ \beta & \Delta \end{bmatrix}

各ブロックを比較すると、

  1. α=(abTA~1b)1\alpha = (a-b^T\tilde{A}^{-1}b)^{-1}
  2. β=A~1b(abTA~1b)1=αA~1b\beta = -\tilde{A}^{-1}b(a-b^T\tilde{A}^{-1}b)^{-1} = -\alpha \tilde{A}^{-1}b
  3. Δ=A~1+A~1b(abTA~1b)1bTA~1=A~1+α(A~1b)(bTA~1)\Delta = \tilde{A}^{-1} + \tilde{A}^{-1}b(a-b^T\tilde{A}^{-1}b)^{-1}b^T\tilde{A}^{-1} = \tilde{A}^{-1} + \alpha (\tilde{A}^{-1}b)(b^T\tilde{A}^{-1})

式(2)より、 A~1b=1αβ\tilde{A}^{-1}b = -\frac{1}{\alpha}\beta となる。これを式(3)に代入する。

Δ=A~1+α(1αβ)(1αβ)T=A~1+α(1αβ)(1αβT)\Delta = \tilde{A}^{-1} + \alpha \left(-\frac{1}{\alpha}\beta\right) \left(-\frac{1}{\alpha}\beta\right)^T = \tilde{A}^{-1} + \alpha \left(-\frac{1}{\alpha}\beta\right) \left(-\frac{1}{\alpha}\beta^T\right)
Δ=A~1+α1α2ββT=A~1+ββTα\Delta = \tilde{A}^{-1} + \alpha \frac{1}{\alpha^2} \beta\beta^T = \tilde{A}^{-1} + \frac{\beta\beta^T}{\alpha}

この式を A~1\tilde{A}^{-1} について解くと、

A~1=ΔββTα\tilde{A}^{-1} = \Delta - \frac{\beta\beta^T}{\alpha}

が示された。

(iii) の証明

二次形式 Q(x)=xTA1xQ(x) = x^T A^{-1} xx1x_1x~\tilde{x} を用いて展開する。

x=[x1x~],A1=[αβTβΔ]x = \begin{bmatrix} x_1 \\ \tilde{x} \end{bmatrix}, \quad A^{-1} = \begin{bmatrix} \alpha & \beta^T \\ \beta & \Delta \end{bmatrix}
Q(x)=[x1x~T][αβTβΔ][x1x~]Q(x) = \begin{bmatrix} x_1 & \tilde{x}^T \end{bmatrix} \begin{bmatrix} \alpha & \beta^T \\ \beta & \Delta \end{bmatrix} \begin{bmatrix} x_1 \\ \tilde{x} \end{bmatrix}
=[x1x~T][αx1+βTx~βx1+Δx~]= \begin{bmatrix} x_1 & \tilde{x}^T \end{bmatrix} \begin{bmatrix} \alpha x_1 + \beta^T \tilde{x} \\ \beta x_1 + \Delta \tilde{x} \end{bmatrix}
=x1(αx1+βTx~)+x~T(βx1+Δx~)= x_1(\alpha x_1 + \beta^T \tilde{x}) + \tilde{x}^T(\beta x_1 + \Delta \tilde{x})
=αx12+x1βTx~+x~Tβx1+x~TΔx~= \alpha x_1^2 + x_1\beta^T \tilde{x} + \tilde{x}^T\beta x_1 + \tilde{x}^T\Delta\tilde{x}

x1βTx~x_1\beta^T \tilde{x} はスカラーなので、その転置 x~Tβx1\tilde{x}^T\beta x_1 と等しい。したがって、

Q(x)=αx12+2(βTx~)x1+x~TΔx~Q(x) = \alpha x_1^2 + 2(\beta^T \tilde{x})x_1 + \tilde{x}^T\Delta\tilde{x}

この式は、 x1x_1 に関する二次関数である。 AA が正定値なので、 A1A^{-1} も正定値であり、その主小行列である α\alphaα>0\alpha > 0 である。したがって、この二次関数は下に凸の放物線であり、最小値を持つ。

最小値は、この二次関数を x1x_1 で微分して 0 とおくことで見つけられる。

Q(x)x1=2αx1+2(βTx~)=0\frac{\partial Q(x)}{\partial x_1} = 2\alpha x_1 + 2(\beta^T \tilde{x}) = 0

これを解くと、最小値を与える x1x_1 の値 x1x_1^* は、

x1=βTx~αx_1^* = -\frac{\beta^T \tilde{x}}{\alpha}

となる。 この x1x_1^*Q(x)Q(x) の式に代入して最小値を計算する。

minx1Q(x)=α(βTx~α)2+2(βTx~)(βTx~α)+x~TΔx~\min_{x_1} Q(x) = \alpha \left(-\frac{\beta^T \tilde{x}}{\alpha}\right)^2 + 2(\beta^T \tilde{x})\left(-\frac{\beta^T \tilde{x}}{\alpha}\right) + \tilde{x}^T\Delta\tilde{x}
=α(βTx~)2α22(βTx~)2α+x~TΔx~= \alpha \frac{(\beta^T \tilde{x})^2}{\alpha^2} - 2\frac{(\beta^T \tilde{x})^2}{\alpha} + \tilde{x}^T\Delta\tilde{x}
=(βTx~)2α2(βTx~)2α+x~TΔx~= \frac{(\beta^T \tilde{x})^2}{\alpha} - 2\frac{(\beta^T \tilde{x})^2}{\alpha} + \tilde{x}^T\Delta\tilde{x}
=(βTx~)2α+x~TΔx~= -\frac{(\beta^T \tilde{x})^2}{\alpha} + \tilde{x}^T\Delta\tilde{x}

ここで、 (βTx~)2=(x~Tβ)(βTx~)=x~T(ββT)x~(\beta^T \tilde{x})^2 = (\tilde{x}^T \beta)(\beta^T \tilde{x}) = \tilde{x}^T (\beta\beta^T) \tilde{x} と書けるので、

minx1Q(x)=x~TΔx~x~TββTx~α=x~T(ΔββTα)x~\min_{x_1} Q(x) = \tilde{x}^T\Delta\tilde{x} - \frac{\tilde{x}^T\beta\beta^T\tilde{x}}{\alpha} = \tilde{x}^T \left( \Delta - \frac{\beta\beta^T}{\alpha} \right) \tilde{x}

設問 (ii) の結果から、 A~1=ΔββTα\tilde{A}^{-1} = \Delta - \frac{\beta\beta^T}{\alpha} である。 したがって、二次形式の最小値は、

minx1xTA1x=x~TA~1x~\min_{x_1} x^T A^{-1} x = \tilde{x}^T \tilde{A}^{-1} \tilde{x}

となり、題意は示された。

(iv)

Dn=detAnD_n = \det A_n 。根据题意, AnA_n 对所有 nn 都是正定矩阵,因此其所有主子式都为正。特别地,行列式 Dn=detAn>0D_n = \det A_n > 0 对所有 nn 成立。

需要证明的不等式为:

Dn+1Dn2Dn1D_{n+1} \leq \frac{D_n^2}{D_{n-1}}

因为 Dn1>0D_{n-1} > 0 ,该不等式等价于:

Dn+1Dn1Dn2D_{n+1} D_{n-1} \leq D_n^2

我们将使用 Desnanot-Jacobi 恒等式(也称为 Sylvester 行列式恒等式或 Dodgson 凝聚法)。对于任意一个 m×mm \times m 矩阵 MM ,该恒等式为:

det(M)det(M1,m1,m)=det(M11)det(Mmm)det(M1m)det(Mm1)\det(M) \det(M_{1,m}^{1,m}) = \det(M_1^1) \det(M_m^m) - \det(M_1^m) \det(M_m^1)

其中, MijM_i^j 表示从 MM 中移除第 ii 行和第 jj 列得到的子矩阵, Mi,jk,lM_{i,j}^{k,l} 表示移除第 i,ji,j 行和第 k,lk,l 列得到的子矩阵。

我们将此恒等式应用于 (n+1)×(n+1)(n+1) \times (n+1) 矩阵 M=An+1M = A_{n+1} 。此时 m=n+1m=n+1

det(An+1)det((An+1)1,n+11,n+1)=det((An+1)11)det((An+1)n+1n+1)det((An+1)1n+1)det((An+1)n+11)\det(A_{n+1}) \det((A_{n+1})_{1,n+1}^{1,n+1}) = \det((A_{n+1})_1^1) \det((A_{n+1})_{n+1}^{n+1}) - \det((A_{n+1})_1^{n+1}) \det((A_{n+1})_{n+1}^1)

接下来,我们根据 An+1A_{n+1} 的结构来确定上式中各项的行列式:

  1. An+1A_{n+1} 本身: det(An+1)=Dn+1\det(A_{n+1}) = D_{n+1}

  2. (An+1)11(A_{n+1})_1^1 是移除第一行和第一列得到的子矩阵。根据 An+1A_{n+1} 的定义,这是一个 n×nn \times n 的主子矩阵,其结构与 AnA_n 完全相同。

(An+1)11=An    det((An+1)11)=Dn(A_{n+1})_1^1 = A_n \implies \det((A_{n+1})_1^1) = D_n
  1. (An+1)n+1n+1(A_{n+1})_{n+1}^{n+1} 是移除最后一行和最后一列得到的子矩阵。同样,这也是 AnA_n
(An+1)n+1n+1=An    det((An+1)n+1n+1)=Dn(A_{n+1})_{n+1}^{n+1} = A_n \implies \det((A_{n+1})_{n+1}^{n+1}) = D_n
  1. (An+1)1,n+11,n+1(A_{n+1})_{1,n+1}^{1,n+1} 是移除第一行、最后一行、第一列和最后一列得到的子矩阵。这是一个 (n1)×(n1)(n-1) \times (n-1) 的中心子矩阵,其结构与 An1A_{n-1} 相同(此步骤要求 n2n \geq 2 )。
(An+1)1,n+11,n+1=An1    det((An+1)1,n+11,n+1)=Dn1(A_{n+1})_{1,n+1}^{1,n+1} = A_{n-1} \implies \det((A_{n+1})_{1,n+1}^{1,n+1}) = D_{n-1}
  1. (An+1)n+11(A_{n+1})_{n+1}^1(An+1)1n+1(A_{n+1})_1^{n+1} 是非主子矩阵。由于 An+1A_{n+1} 是一个对称矩阵,其元素满足 (An+1)ij=(An+1)ji(A_{n+1})_{ij} = (A_{n+1})_{ji} 。我们来比较这两个子矩阵。 令 B=(An+1)1n+1B = (A_{n+1})_1^{n+1}C=(An+1)n+11C = (A_{n+1})_{n+1}^1CC(i,j)(i,j) 元素为 (An+1)i,j+1(A_{n+1})_{i, j+1}BB(j,i)(j,i) 元素为 (An+1)j+1,i(A_{n+1})_{j+1, i} 。 因为 An+1A_{n+1} 是对称的, (An+1)i,j+1=(An+1)j+1,i(A_{n+1})_{i, j+1} = (A_{n+1})_{j+1, i} 。所以 Cij=BjiC_{ij} = B_{ji} 。 这意味着 C=BTC = B^T 。因此,它们的行列式相等:
det((An+1)n+11)=det(((An+1)1n+1)T)=det((An+1)1n+1)\det((A_{n+1})_{n+1}^1) = \det(((A_{n+1})_1^{n+1})^T) = \det((A_{n+1})_1^{n+1})

将以上结果代入 Desnanot-Jacobi 恒等式:

Dn+1Dn1=DnDndet((An+1)1n+1)det((An+1)n+11)D_{n+1} \cdot D_{n-1} = D_n \cdot D_n - \det((A_{n+1})_1^{n+1}) \cdot \det((A_{n+1})_{n+1}^1)
Dn+1Dn1=Dn2(det((An+1)1n+1))2D_{n+1} D_{n-1} = D_n^2 - (\det((A_{n+1})_1^{n+1}))^2

Cn=det((An+1)1n+1)C_n = \det((A_{n+1})_1^{n+1}) 。由于 An+1A_{n+1} 的元素是实数, CnC_n 是一个实数,所以 Cn20C_n^2 \geq 0 。 因此,我们得到:

Dn+1Dn1=Dn2Cn2Dn2D_{n+1} D_{n-1} = D_n^2 - C_n^2 \leq D_n^2

我们已经证明了 Dn+1Dn1Dn2D_{n+1} D_{n-1} \leq D_n^2 。因为 An1A_{n-1} 是正定矩阵,所以 Dn1=detAn1>0D_{n-1} = \det A_{n-1} > 0 。两边同时除以 Dn1D_{n-1} ,不等号方向不变:

Dn+1Dn2Dn1D_{n+1} \leq \frac{D_n^2}{D_{n-1}}

证明完毕。