0.7 分块矩阵
类似于集合的划分,一个矩阵的分块是把该矩阵完全地分成一些互不相交的子矩阵。使得原矩阵的每一个元落到日只落到一个分块子矩阵中。为识别一些有用的结构,矩阵分块往往是一个方便的方法。
0.7.1子矩阵设 Λ∈Mm,n(F) ,对于指标集 α⊆{1,…,m} 和 β⊆{1,…,n} ,把 A 中位于标号 α 的诸行和标号为 β 的诸列的(子)矩阵记作 A(α,β) 例如
147258369({1,3},{1,2,3})=[172839]. 如果 m=n 且 β=α ,了矩阵 A(α,α) 称为 Λ 的主子矩阵,简记为 A(α) 。说明一个子矩阵或一个主子矩阵是经划去那些行或列得来的,常常要比说它们包括那些行或列要方便些。这可以通
16
过补充指标集来完成。例如, A(α′,β′) 是划去标号为 α 的诸行与标号为 β 的诸列后的结果。
矩阵 A 的子方阵的行列式称为 A 的一个子式。如果子矩阵是一个主子矩阵,那么其子式称为主子式。那些出现在 Laplace 展开式 (0.3.1) 中的带正负号的子式 [(−1)i−1detAij] 称为 A 的代数余子式。约定,空主子式是 1,即 detA(ϕ)=1 。
0.7.2 分块与乘法 如果 α1,⋯,αs 组成 {1,⋯,m} 的一个划分,而 β1,⋯,βs 组成 {1,⋯,n} 的一个划分,那么诸矩阵 A(αi,βj) , 1⩽i⩽t , 1⩽j⩽s ,构成 A∈Mm,n(F) 的一个分块。如果 A∈Mm,n(F) 和 B∈Mn,p(F) 已经分块,且使 {1,⋯,n} 的两个划分重合,那么就说这两个矩阵分块是共形的。这时,
17
[AB](αi,γj)=k∑sA(αi,βk)B(βk,γj), 其中 A(αt,βk) 和 B(βk,γj) 是 A 和 B 的共形分块.等式左边是,乘积 AB 的一个子矩阵(按通常方式计算),而右边的每一个被加项是一个标准的矩阵乘积.因此,共形的分块矩阵的乘法与通常的矩阵乘法相仿.当各被加项有相同的分块时,分块矩阵的加法也是有意义的.
0.7.3 分块矩阵的逆 求出非奇异分块矩阵 A 的逆的相应子块,即用相应的分块形式给出分块矩阵的逆,有时很有用。可以采用各种不同的,但彼此等价的方式来求分块矩阵的逆——假定 A∈Mn(F) 的某些子矩阵及 A−1 也是非奇异的。为简单起见,设 A 是如下的分块
A=[A11A21A12A22], 其中, An∈Mn1(F) , i=1 ,2日 n1+n2=n ,关于 A−1 的相应分块形式有一个有用的表示式
[[A11−A12A22−1A21]−1[A21A11−1A12−A22]−1A21A11−1A11−1A12[A21A11−1A12−A22]−1[A22−A21A11−1A12]−1], 其中,假定所有有关的逆存在。或者,用一般的指标集记号,可以记
A−1(α)=[A(α)−A(α,α′)A(α′)1A(α′,α)]−1, 以及
A′(α,α′)=A(α)−1A(α,α′)[A(α′,α)A(α)−1A(α,α′)−A(α′)]−1, 仍假定有关的逆存在。还可以写出其余的表示式。注意, A−1(α) 是 A−1 的子矩阵,而 A(α)−1 是 A 的一个子矩阵的逆,并且这两个矩阵一般不相同。
0.7.4 小秩修正矩阵的逆 如果已知某个矩阵的逆,了解该矩阵再加上一个“小”秩矩阵时,其逆如何变化,这同样是一个有意义的问题。有这样的简便公式,只要修正矩阵的形式足够简单,就可以使新逆的计算比从头开始计算要简便。设非奇异矩阵 A∈Mn(F) 的逆 A−1 已知,考虑
B=A+XRY, 其中, X 是 n×r 矩阵,而 R 是 r×r 非奇异矩阵.如果 B 是非奇异的,那么
B1=A−1⋅A1X(R−1+YA1X)−1YA−1. 如果 r 比 n 小得多,那么求 R 和 R−1+YA−1X 的逆可能要比求 B 的逆更为容易,并且,如果
A 是容易求逆的,且有使乘法简化的形式,那么,采用这个公式可能要胜过直接求 B 的逆。例如,如果修正矩阵有秩 1, X 是 n×1 的, Y 是 1×n 的,且 R=[1] ,上述公式就变成
B−1=A−1−1+YA−1X1A−1XYA−1, (注意,此时 XY=B−A ).特别地,如果
B=I+xy1 对于 x,y∈Fn,l∈Mn(F) 成立,那么,只要 yix=−1 ,就有
B′=I−1+ylx1xyl. 0.8 行列式(续)
关于行列式的一些补充材料和恒等式对于某些论题的阐述很有用。其中大部分内容在基础线性代数中不易找到。
0.8.1 复合矩阵 设矩阵 Λ∈Mm,n(F) ,其某个阶数的所有子式组成的阵列称为 Λ 的复合矩阵。特别地,它的 α,β 元为 detA(α,β) 的 (km)×(kn) 矩阵叫做 A 的第 k 次复合矩阵,记作 Ck(A) 。这里, α⊆{1,⋯,m} 和 β⊆{1,⋯,n} 都是基数为 k⩽min{m,n} 的指标集,按通常的辞典式次序排列,即 {1,2,4} 排在 {1,2,5} 之前, {1,2,5} 排在 {1,3,4} 之前,等等。例如,如果
A=147258369. 则
C2(A)=det[1425]det[1728]det[4758]det[1436]det[1739]det[1769]det[2536]det[2839]det[−586−9]=−3−6−3−612−6−3−6−3, 如果 A∈Mk,n(F) ,且 B∈Mk,n(F) ,那么
C,(AB)=C,(A)C,(B),r⋮1minmk,n. 另外,还有
Ci(tA)=t′Ci(A),i∈F. 若 I∈Mn , 则 Ck(I)=I∈M(1) .
若 A∈Mn 是非奇异的,则 Ck(A)−1=Ck(A1)
若 A∈Mm,n(F) ,则 Ck(AT)=Ck(A)T .
及
若 A∈Mm,n(C) ,则 Ck(A∗)=Ck(A)∗ 。
0.8.2 经典伴随和逆 如果 Λ∈Mn(F) ,由诸代数余子式
bij=(−1)r+jdetA({j}′,{i}′) 组成的转置矩阵 B=[bij]∈Mn(F) 称为 Λ 的(经典)伴随,常记作 adjA 。有时用转置伴随来代替伴随,以免与Hermite伴随 A∗ 相混淆。注意,
adjΛ=ECn−1(Λ)′E0, 其中
E=1−101−10⋱±1. 用Laplace展开式计算行列式的公式说明.
(adjΛ)A=Λ(adjΛ)=(detA)I. 因而,如果 A 是非奇异的( detA=0 ),那么
A−1=detA1adjA. 一般说来,用伴随来数值计算矩阵的逆不是可取的方法,但是伴随对于给出逆的解析表达式是有用的.
⊖ 原书给出的这个关系式是错误的(但书中尚未用到这个关系式),adj A 与 (A) 应是如下关系:
adjA=J0Cn−1(A)tJ1=J1Cn−1(A)tJ(n为 偶 数) adjA=J0Cn−1(A)TJ1=J1Cn−1(A)TJ1(n为 奇 数). 其中 n 阶方阵 J0 和 J1 分别为:
J0…10⋱−110−1,J1=11⋱−11−1. 0.8.3 Cramer法则 当 A∈Mn(F) 非奇异时,Cramer法则是求线性方程组 Ax=b 唯一解的一种方法。它和逆的伴随表示有相同的计算手续,一般,只有在需要解析地给出解向量的个别解析分量时,这种方法才有用。如果 xi 是解向量 x∈Fn 的第 i 个分量,那么,Cramer法则可述为公式
xi=detAdet(A⋅b). 记号 A←b 表示 Mn 中第 i 列是 b ,其余各列与 A 的对应列相同的矩阵。Cramer 法则可直接由行列式的乘法性质推出。把方程组 Ai=b 改写成
A(I←r)=A←b, 然后两边取行列式(利用乘法性质)可得
(detA)det(I←x)=det(A←b). 而 det(I−x)=x ,因而公式得证。
0.8.4 逆的子式 推广非奇异矩阵的逆的伴随公式,有如下重要公式:
detA1(α′,β′)−(−1)(∑i=1n∑j=1n)detAdetA(β,α), 它把 Λ1 的诸子式与 A∈Mλ(F) 的诸子式联系起来。对于主子阵,这个公式有简单的形式
detA′(α′)=detAdetA(α). 0.8.5 Schur补和行列式公式 对于给定的矩阵 A∈Mn(F) ,设 α⊆{1,…,n} 是使 A(α) 非奇异的指标集.记 A(α) 的逆为 A(α)−1 .用 α 和 α′ 对 A 作 2×2 分块,据此,关于det A 的重要公式是
detA−detA(α)det[A(α′)−A(α′,α)A(α)′A(α,α′)]. 注意,这个公式推广了(0.3.1)中的关于 2×2 矩阵行列式的普通公式。称特殊矩阵
A(α′)−A(α′,α)A(α)−1A(α,α′) 为 A 的Schur补,将
[Λ11Λ21Λ12Λ22]右 乘 以[I0−A11−1A12I], 然后把 A:1 和 A(α) 等同起来,就可验证 detA 和 Schur 补公式成立。注意,Schur 补已在 A−1 的分块中出现过[见(0.7.3)]。
0.8.6 Sylvester 恒等式 设 α⊆{1,…,n} 是固定的指标集,设 B=[bij]∈Mn−k(F) 由
bij=detA(α⋃{i},α⋃{j}) 所定义,其中 k 是 α 的基数, i,j∈{1,…,n} 是不包含在 α 中的指标, A∈Mn(F) 。另一个有用的行列式恒等式是
detB⋅[detA(u)]−nkdetA 0.8.7 Cauchy-Binet 公式 这个有用的公式是可以想起来的,这是因为它看上去与矩阵的乘法公式相类似。它等价于复合矩阵的乘法性质(见 0.8.1),所以这不是偶然的巧合。设 A∈Mm,k(F) , B∈Mk,p(F) 和 C=AB 。再设 1⩽r⩽min{m,k,n} , α⊆{1,⋯,m} 和 β⊆{1,⋯,n} 都是基
数为 r 的指标集。关于 C 的 α,β 子式的表示式是
detC(α,β)=γ∑detA(α,γ)detB(γ,β), 其中和式取遍基数为 r 的所有指标集 γ⊆{1,…,k} .
0.8.8子式间的关系已知 A∈Mm,n(F) ,给定基数为 k 的固定指标集 α⊆{1,…,m} ,当 ω⊆{1,…,n} 取遍基数为 k 的有序指标集时,诸子式
detA(α,ω) 不是代数无关的,因为在各子矩阵中诸子式多于诸子矩阵中的各不相同的元,在这些子式当中,二次关系是知道的。设 i1,i2,…,ik∈{1,…,n} 是 k 个互异的指标,不一定取自然顺序,又设
A(α;i1,…,ik) 表示这样一个矩阵,它的各行用 α 标号,而它的第 j 列是 A(α,{1,…,n}) 的 ij 列。这与前述记号的差别是,列可以不按自然顺序,如在 A({1,3};4,2) 中,它的第 i 列有 A 中的1,4元和3,4元。于是有关系式
detA(α;i1,…,ik)detA(α;j1,…,jk)=∑i=1kdetA(α;i1,…,ir−1,jt,is+1,…,ir)detA(α;j2,…,jt−1,i,jt+1,…,jk). 它对于每个 s=1,…,k 和互异指标的所有序列
i1,…,ik∈{1,…,n}和j1,…,jk∈{1,…,n} 成立.