分块矩阵法是高等代数中处理矩阵相关问题的最重要的方法之一, 其重要性可以说怎么强调都不过分. 分块矩阵法的核心思想是根据具体问题构造适当的分块矩阵, 然后运用广义初等变换, 将某些子块消为零块, 得到特殊的分块矩阵从而解决问题. 该方法几乎贯穿了线性代数的始终, 在矩阵求逆、矩阵秩不等式、行列式、线性方程组、线性变换、二次型等方面有着广泛的应用.

例如, 证明行列式的乘法公式|AB|=|A||B|时, 我们就构造了分块矩阵

(AOEB).

作广义初等变换

(AOEB)c1B+c2(AABEO)r1r2(EOAAB).

取行列式, 得

|A||B|=|AOEB|=(1)n|EOAAB|=(1)n|En||AB|,

|AB|=|A||B|.

再例如, 若子矩阵块A可逆,利用初等变换将分块矩阵(ABCD)
的子块C消成O, 即

(ABCD)r1CA1+r2(ABODCA1B).

得到行列式

|ABCD|=|ABODCA1B|=|A||DCA1B|.

该式有时也称为行列式的降阶公式.

§1 由于广义初等变换不改变矩阵的秩,所以 关于矩阵秩的等式或不等式的证明, 分块矩阵法可以大显神通. 使用该方法的主要思路是根据等式或不等式的一边,利用下面例2中的等式或不等式,构造合适的分块矩阵,对其实施分块矩阵的初等变换,得到等式或不等式的另一边的秩关系.

例1Am×n,Bn×l, 则r(AB)min{r(A),r(B)}.

证明 由于

(A,AB)(EBOE)=(A,O),(EOAE)(BAB)=(BO),

所以

r(AB)r(A,AB)=r(A,O)=r(A),r(AB)r(BAB)=r(BO)=r(B),

r(AB)min{r(A),r(B)}.

例2 r(A)+r(B)=r(AOOB)r(AOCB).

证明r(A)=r1, r(B)=r2, 则存在可逆矩阵P1,Q1P2,Q2 使得

P1AQ1=(Er1O),P2BQ2=(Er2O).

从而

(P1P2)(AB)(Q1Q2)=(Er1OEr2O).

所以由(1)得

r(AB)=r1+r2=r(A)+r(B).

另一方面, 由于通过初等变换,

(AOCB)(Er1OD11D12Er2D21D22O)(Er1OOOEr2ODO)(Er1Er2DO),

所以

r(AOCB)=r1+r2+r(D)r1+r2=r(A)+r(B).

例3A,B是同阶矩阵,则
|r(A)r(B)|r(A±B)r(A)+r(B).

证明 由于

(AB)(AOAB)(AAAA+B),

所以

r(A+B)r(AAAA+B)=r(AB)=r(A)+r(B).

另一方面, 不妨设r(A)r(B), 则由

(A+BB)(A+BBOB)(ABBB),

可得

r(A)r(ABBB)=r(A+BB)=r(A+B)+r(B),

r(A)r(B)r(A+B).

而且,

r(A)=r(AB+B)r(AB)+r(B),

所以 r(A)r(B)r(AB),.
类似地,

r(B)r(A)r(BA)=r(AB).

所以
|r(A)r(B)|r(AB),,

|r(A)r(B)|r(A±B).

例4 [第一降阶定理]
A是可逆矩阵, (ABCD)(r,mr)×(r,mr)分块矩阵. 证明r(ABCD)=r(A)+r(DCA1B).

证明 因为

(ErOCA1Emr)(ABCD)(ErA1BOEmr)=(AOODCA1B),

(ErOCA1Emr)(ErA1BOEmr)都是非异矩阵, 所以

r(ABCD)=r(A)+r(DCA1B).

例5 [第二降阶定理]
A,D分别是r,s阶非奇异矩阵, BFr×s, CFs×r, 则

r(DCA1B)=r(D)r(A)+r(ABD1C).

证明 由第一讲解定理得

r(ABCD)=r(A)+r(DCA1B).

另一方面, 由D非奇异知

(ErBD1OEmr)(ABCD)(ErOD1CEmr)=(ABD1COOD),

所以

r(ABCD)=r(D)+r(ABD1C).

r(DCA1B)=r(D)r(A)+r(ABD1C).

例6Am×n矩阵, 证明

r(EmAAT)r(EnATA)=mn.

*证明???

r(EmAAT)=r(EmAEnAT)=r(Em)r(En)+r(EnATEmA)=mn+r(EnATA).

r(EmAAT)r(EnATA)=mn.

例7 [Sylvester不等式]
A,Bn阶方阵,则r(A)+r(B)r(AB)+n.

** 证明** 对分块矩阵作广义初等变换, (AO\EB)(OAB EB)(OAB EO)(EO\OAB),.
由(3)得

r(A)+r(B)r(AOEB)=r(EOOAB)=r(AB)+n,

r(A)+r(B)nr(AB).

例[8 复旦大学1999,武汉大学2001,厦门大学2001,北京理工大学2003,西安电子科技大学2004,浙江大学2007,北京大学2009,中南大学2010,兰州大学2020]
A,B,C分别是m×n,n×s,s×t矩阵. 证明Frobenius不等式

r(ABC)+r(B)r(AB)+r(BC).

证明 (方法一) 利用分块初等变换,考虑

(ABCOOB)(ABCABOB)(OABBCB)(ABOBBC),

于是

r(ABC)+r(B)=r(ABCOOB)=r(ABOBBC)r(ABOOBC)=r(AB)+r(BC).

(方法二) 利用满秩分解与Sylvester不等式. 设r(B)=r, 则存在行、列满秩矩阵H,L使得B=LH. 则

r(ABC)=r(ALHC)r(AL)+r(HC)rr(AB)+r(BC)r(B).

下面的例子中, 由于A+B未必可逆,所以下面的一步使用的不是(广义)初等变换

(A+BE)(AAAA+B)=左乘(A+B)((A+B)A(A+B)AAA+B)

由于矩阵乘积的秩不大于每一个因子的秩,所以

r(AAAA+B)((A+B)A(A+B)AAA+B)

例9 [武汉大学2008]
A,B是数域\mF上的n阶方阵,AB=BA. 证明:
r(A+B)r(A)+r(B)r(AB).

证明 (方法一) 利用线性方程组与维数公式,参见??????.

(方法二)
因为

(AOOB)(AOAB)(AAAA+B)左乘(A+B)((A+B)A(A+B)AAA+B)(O(A+B)ABA+B)(ABOBA+B)

由于A+B未必可逆,所以左乘A+B不是初等变换,从而

r(A)+r(B)=r(AOOB)=r(AAAA+B)r((A+B)A(A+B)AAA+B)=r(ABOBA+B)(ABOOA+B)=r(AB)+r(A+B).

例10A,Bm×n,n×m矩阵,证明

r(AABA)=r(A)+r(EnBA)n.

证明 因为

(AOOEnBA)(AOBAEnBA)(AABAEn)(AABAOBAEn)(AABAOOEn)

所以

r(A)+r(EnBA)=r(AOOEnBA)=r(AABAOOEn)=r(AABA)+n.

§2 由于矩阵的标准形(等价标准形、Jordan 标准形、Frobenius标准形)都是分块矩阵,因此运用标准形解决矩阵问题常常结合分块矩阵的方法.

例11A,B分别是复数域上的m,n阶矩阵,则A,B没有公共特征值当且仅当矩阵方程AX=XB只有零解.

** 证明**
必要性:(反证法)假设X0AX=XB的非零解,设存在可逆矩阵P,Q使得

PX0Q=(ErOOO),r1.

AX0=X0B得到

(PAP1)(PX0Q)=(PX0Q)(Q1BQ).

PAP1=(A1A2A3A4),Q1BQ=(B1B2B3B4),

(A1A2A3A4)(ErOOO)=(ErOOO)(B1B2B3B4),

(A1OA3O)=(B1B2OO).

从而可得A1=B1,A3=O,B2=O. 从而A1=B1的特征值是A,B的公共特征值,矛盾!所以AX=XB只有零解.

充分性:(反证法)假设λ0A,B的公共特征值,因为fB(λ)=fBT(λ), 所以λ0也是A,BT的公共特征值, 故存在可逆矩阵P,Q使得

P1AP=(λ0A20A4),Q1BTQ=(λ0B20B4),

从而

A=P(λ0A20A4)P1,B=(QT)1(λ00B2TB4T)QT.

X0=P(100O)QTO, 则AX)=X0B, 矛盾!所以A,B没有公共特征值.

该例的几何解法参见???. 该例的结论可用于解决如下问题:

1、A,B没有公共特征值 (fA(λ),fB(λ))=1 fB(A)可逆 矩阵方程AX=XB只有零解 线性变换σ:Cm×nCm×n,XAXXB可逆.

2、(中国科学技术大学,2013)设A,Bn阶复方阵,Cn×n上的线性变换\A:XAXXB. 证明:\A可逆的充分必要条件是AB无公共特征值.

3、(武汉大学,2013)设λ1,λ2,,λnn阶实方阵A的全部特征值,但λi(i=1,2,,n)不是A的特征值. 定义Rn×n的线性变换

σ(X)=ATX+XA,XRn×n.

证明:(1) σ是可逆线性变换;

(2) 对任意实对称矩阵C, 存在唯一的实对称矩阵B使得ATB+BA=C.

证明
(1) 只需证明σ是单射. 设σ(X)=ATX+XA=O, 则XA=ATX, 于是

fA(AT)X=XfA(A)=O.

由于λi(i=1,2,,n)不是A的特征值, 所以fA(AT)可逆,从而X=O, σ是单射,从而是双射,故可逆.

(2) 由(1)知存在唯一的矩阵B使得σ(B)=C. 由于

σ(BT)=ATBT+BTA=(ATB+BA)T=CT=C,

所以BT=B, 即B是对称矩阵.

4、设A,B分别是m,n阶方阵,Cm×n阶矩阵. 证明:AX+XB=C有唯一解的充要条件是对A,B的任一特征值λi,μj, 都有λi+μj0.

(提示:只需证AX+XB=O只有零解当且仅当λi+μj0, 这等价于AXXB=O只有零解当且仅当A,B没有公共特征值.)

矩阵方程AXXB=C称为Sylvester方程,在控制理论中有着重要应用.

例12
A,B分别是复数域上的m,n阶矩阵,且A,B没有公共特征值,对任意的CCm×n, 若矩阵方程AXXB=C有解,则有唯一解.

证明 (反证法) 假设矩阵方程AXXB=C有两个不同的解X1,X2, 则X1X2矩阵方程AXXB=O的非零解,与??? 矛盾!所以 矩阵方程AXXB=C有唯一解.

[Roth第二定理]
矩阵方程AXXB=C有解当且仅当(ACOB)(AOOB)相似,见???.

**例13 ** 设A,B分别是复数域上的m,n阶矩阵,若A,B没有公共特征值,则D=(ACOB)可对角化当且仅当A,B可对角化.

** 证明** 因为A,B没有公共特征值,所以矩阵方程AXXB=C有唯一解,设为X0. 由于

(EX0OE)(ACOB)(EX0OE)1=(AAX0+X0B+COB)=(AOOB),

所以D可对角化当且仅当A,B可对角化.

例14A,B分别是复数域C上的n阶与m阶矩阵, 0<rmin{n,m}. 证明:如果AXXB=O有秩为r 的矩阵解, 则AB至少有r个公共特征值(重根按重数计算).

证明CAXXB=O的矩阵解, 且r(C)=r. 则存在可逆矩阵PCn×n, QCm×m, 使得

C=P(ErOOO)Q.

由于AC=CB, 所以

AP(ErOOO)Q=P(ErOOO)QB,

从而

P1AP(ErOOO)=(ErOOO)QBQ1.

P1AP,QBQ1写成分块矩阵的形式

P1AP=(A11A12A21A22),QBQ1=(B11B12B21B22)

(A11OA21O)=(B11B12OO).

由此得A11=B11,A21=O,B12=O.
所以

P1AP=(A11A12OA22),QBQ1=(A11OB12B22).

所以A,B的特征多项式分别为

fA(λ)=|λErA11||λEA22|,fB(λ)=|λErA11||λEB22|.

它们有公因式|λErA11|, 因而AB至少有r 个公共特征值(重根按重数计算).

例15 [Roth第一定理]
A,B,C分别为m×n, n×mn×n 矩阵, 则\r(AOOB)=r(AOCB)的充要条件为存在矩阵X,Y使得
XABY=C.

证明 (EOXE)(AOCB)(EOYE)=(AOOB)即得.

P1AQ1=(ErOOO)=ΔH1,P2BQ2=(EsOOO)=ΔH2.
(1)(P1O\OP2)(AO\OB)(Q1O\OQ2)=(Er O Es O),
(2)(P1O\OP2)(AO CB)(Q1O\OQ2)=(P1AQ1O P2CQ1P2BQ2) =(ErO OO D1D2EsO D3D4OO)=(H1O DH2).
从(2.5)中消去D1, 即

(ErOED1OEsOOOOE)(ErOOOD1D2EsOD3D4OO)=(ErOOOOD2EsOD3D4OO).

此过程可简记为

(D1OOO)(ErOOO)+(D1D2D3D4)=(OD2D3D4).

继续消去D3,D2的过程可简记为

(OOD3O)(ErOOO)+(OD2D3D4)=(OD2OD4).

(EsOOO)(OD2OO)+(OD2OD4)=(OOOD4).

由题设, (2.4)与(2.5)式右端的秩相等, 所以D4=O. 从而该消去过程相当于存在矩阵U=(OOD3O)+(D1OOO),V=(OD2OO)使得

UH1+H2V+D=O.

U(P1AQ1)+(P2BQ2)V=P2CQ1,

或等价地,

(P21UP1)AB(Q2VQ11)=C.

X=P21UP1, Y=Q2VQ11, 即证.

例16A,B分别为m×nn×m矩阵, 则r(A)+r(B)=r(AB)+n 的充要条件为存在矩阵X,Y使得
XABY=En.

证明 由Sylvester不等式的证明过程知只需证明

r(AOOB)=r(AOEB)X,Y,s.t.XABY=En.

这由上题即得.

§3 通过构造分块矩阵,再利用广义初等变换将其化为分块上(下)三角矩阵是解决行列式等式的常用方法.

例17A,B,C,Dn阶方阵, |A|0. 证明

(1) 若AC=CA,\ 则
|ABCD|=|ADCB|;

(2) 若 AB=BA,\ 则
|ABCD|=|DACB|.

证明 (1),若A可逆, 则
\begin{eqnarray}(EO\-A1CE)(AB\CD)&=&\begin{pmatrix}A&B\
-A{-1}CA+C&-ACB+D\end{pmatrix}$$12pt]
&\begin{tabular}{c}
\scriptsize{AC=CA}\
\hline\hline
\end{tabular}&(AB\OA1CB+D),,\end{eqnarray
}
|AB\CD|=|A||A1CB+D|=|ADCB|,.

(2) 类似可证.

例18 [中国科学院大学2010]
A,B分别是n×mm×n 矩阵. 证明

|EnAB|=|EmBAEn|=|EmBA|.

证明 一方面,

(EOAE)(EmBAEn)=(EmBOEnAB)

|EmB AEn|=|Em||EnAB|=|EnAB|,.
另一方面,

(EBOE)(EmBAEn)=(EmBAOOEn)

|EmB AEn|=|EmBA||En|=|EmBA|,.
所以

|EnAB|=|EmBAEn|=|EmBA|.

[上海交通大学2005]
A,Bn阶方阵,满足
r(EAB)+r(E+BA)=n. 求证:r(A)=n.

(提示: 先证r(E+BA)=r(E+AB), 从而r(EAB)+r(E+AB)=n, 故(AB)2=E, A可逆从而r(A)=n, 参见???)
\end{remark}

下面的例子称为特征多项式的降阶公式.

例19 [上海大学2012]
A,B分别是n×mm×n 矩阵, λ0. 证明

|λEnAB|=λnm|λEmBA|.

*证明 一方面,

(EnOλ1AE)(λEmBAEn)=(λEmBOEnλ1AB)

故$$|λEmB AEn|
=|\lambda E_m||E_n-\lambda{-1}AB|=\lambda |\lambda E_n-AB|,.$$
另一方面,

(EBOE)(λEmBAEn)=(λEmBAOOEn)

|λEmB AEn|=|λEmBA||En|=|λEmBA|,.
所以

|λEnAB|=λnm|λEmBA|.