跳到主要内容

京都大学 情報学研究科 システム科学専攻 2020年8月実施 線形代数

Author

思齐塾, 祭音Myyura

Description

正方行列 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}}

となることを示せ。

题目描述

定义实对称 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

解答

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}}

证明完毕。