京都大学 情報学研究科 システム科学専攻 2020年8月実施 線形代数
Author
思齐塾, 祭音Myyura
Description
正方行列 An を以下のように定義する。
A1=[a0],A2=[a0a1a1a0],…,An+1=a0a1a2⋮ana1a0a1⋮an−1a2a1a0⋮an−2⋯⋯⋯⋱⋯anan−1an−2⋮a0
いま, すべての n について An を正定値対称行列とするとき, n≥2 について,
detAn+1≤detAn−1(detAn)2
となることを示せ。
题目描述
定义实对称 Toeplitz 方阵序列
A1=[a0],A2=[a0a1a1a0],
以及一般的
An+1=a0a1a2⋮ana1a0a1⋮an−1a2a1a0⋮an−2⋯⋯⋯⋱⋯anan−1an−2⋮a0.
假设对每个 n,An 都是正定对称矩阵。证明对所有 n≥2,
detAn+1≤detAn−1(detAn)2.
Kai
解答
设 Dn=detAn 。根据题意, An 对所有 n 都是正定矩阵,因此其所有主子式都为正。特别地,行列式 Dn=detAn>0 对所有 n 成立。
需要证明的不等式为:
Dn+1≤Dn−1Dn2
因为 Dn−1>0 ,该不等式等价于:
Dn+1Dn−1≤Dn2
我们将使用 Desnanot-Jacobi 恒等式(也称为 Sylvester 行列式恒等式或 Dodgson 凝聚法)。对于任意一个 m×m 矩阵 M ,该恒等式为:
det(M)det(M1,m1,m)=det(M11)det(Mmm)−det(M1m)det(Mm1)
其中, Mij 表示从 M 中移除第 i 行和第 j 列得到的子矩阵, Mi,jk,l 表示移除第 i,j 行和第 k,l 列得到的子矩阵。
我们将此恒等式应用于 (n+1)×(n+1) 矩阵 M=An+1 。此时 m=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)
接下来,我们根据 An+1 的结构来确定上式中各项的行列式:
-
An+1 本身: det(An+1)=Dn+1 。
-
(An+1)11 是移除第一行和第一列得到的子矩阵。根据 An+1 的定义,这是一个 n×n 的主子矩阵,其结构与 An 完全相同。
(An+1)11=An⟹det((An+1)11)=Dn
- (An+1)n+1n+1 是移除最后一行和最后一列得到的子矩阵。同样,这也是 An 。
(An+1)n+1n+1=An⟹det((An+1)n+1n+1)=Dn
- (An+1)1,n+11,n+1 是移除第一行、最后一行、第一列和最后一列得到的子矩阵。这是一个 (n−1)×(n−1) 的中心子矩阵,其结构与 An−1 相同(此步骤要求 n≥2 )。
(An+1)1,n+11,n+1=An−1⟹det((An+1)1,n+11,n+1)=Dn−1
- (An+1)n+11 和 (An+1)1n+1 是非主子矩阵。由于 An+1 是一个对称矩阵,其元素满足 (An+1)ij=(An+1)ji 。我们来比较这两个子矩阵。
令 B=(An+1)1n+1 和 C=(An+1)n+11 。
C 的 (i,j) 元素为 (An+1)i,j+1 。
B 的 (j,i) 元素为 (An+1)j+1,i 。
因为 An+1 是对称的, (An+1)i,j+1=(An+1)j+1,i 。所以 Cij=Bji 。
这意味着 C=BT 。因此,它们的行列式相等:
det((An+1)n+11)=det(((An+1)1n+1)T)=det((An+1)1n+1)
将以上结果代入 Desnanot-Jacobi 恒等式:
Dn+1⋅Dn−1=Dn⋅Dn−det((An+1)1n+1)⋅det((An+1)n+11)
Dn+1Dn−1=Dn2−(det((An+1)1n+1))2
令 Cn=det((An+1)1n+1) 。由于 An+1 的元素是实数, Cn 是一个实数,所以 Cn2≥0 。
因此,我们得到:
Dn+1Dn−1=Dn2−Cn2≤Dn2
我们已经证明了 Dn+1Dn−1≤Dn2 。因为 An−1 是正定矩阵,所以 Dn−1=detAn−1>0 。两边同时除以 Dn−1 ,不等号方向不变:
Dn+1≤Dn−1Dn2
证明完毕。