重要矩阵、满秩条件与矩阵幂

汇总实矩阵、零矩阵、正交矩阵、实对称矩阵、秩一矩阵、行列满秩矩阵及矩阵幂的核心结论与使用方法。

来源:邂逅遗憾 26 考研数学思维课(数一)整理。

本文将重要矩阵、满秩条件与矩阵幂合并整理,保留原稿中的核心结论、证明过程和帮助理解结论的小例子,不展开整理完整例题。

一、实矩阵、零矩阵与可逆矩阵

1. 实矩阵的定义与特征值

所有元素均为实数的矩阵称为实矩阵。考研数学讨论的矩阵通常都是实矩阵,不考虑矩阵元素本身为复数的情况。

实矩阵的特征值不一定都是实数,但实矩阵的非实特征值一定成共轭对出现。若

λ=a+bi\lambda=a+bi

为特征值,则

λ=abi\overline\lambda=a-bi

也必为特征值。

证明: 实矩阵 AA 的特征多项式 λEA|\lambda E-A| 的系数均为实数。实系数多项式的非实根必成共轭对出现,因此 a+bia+biabia-bi 同为 AA 的特征值。

2. 零矩阵的秩判定

所有元素均为零的矩阵称为零矩阵,记作 OO。对于任意矩阵 AA

A=O    r(A)=0.A=O\iff r(A)=0.

证明: A=OA=O 时不存在非零行或非零列,所以 r(A)=0r(A)=0;反之,r(A)=0r(A)=0 说明 AA 中没有任何非零元素,故 A=OA=O

3. 零矩阵的迹判定

AA 为任意实矩阵,则有两个重要的充要条件:

A=O    tr(AAT)=0,A=O\iff \operatorname{tr}(AA^T)=0,

以及

A=O    tr(ATA)=0.A=O\iff \operatorname{tr}(A^TA)=0.

证明:

tr(AAT)=tr(ATA)=i,jaij2.\operatorname{tr}(AA^T) =\operatorname{tr}(A^TA) =\sum_{i,j}a_{ij}^2.

由于 AA 是实矩阵,每个 aij20a_{ij}^2\ge0。若平方和为零,则每一项都必须为零,所以 A=OA=O。这也可统一写成

A=O    AF=0.A=O\iff \lVert A\rVert_F=0.

必须注意:tr(A)=0\operatorname{tr}(A)=0 不能推出 A=OA=OAB=OAB=O 也不能直接推出 A=OA=OB=OB=O

4. 可逆矩阵的等价条件

nn 阶方阵 AA,下列条件等价:

A 可逆    A0    r(A)=n    Ax=0 只有零解.A\text{ 可逆} \iff |A|\ne0 \iff r(A)=n \iff Ax=0\text{ 只有零解}.

AA 可逆,则

A1=1AA,A^{-1}=\frac1{|A|}A^*,

并且

AA1=A1A=E,AA=AA=AE.AA^{-1}=A^{-1}A=E, \qquad AA^*=A^*A=|A|E.

证明: A0|A|\ne0 时,伴随矩阵公式给出 A1=A1AA^{-1}=|A|^{-1}A^*,所以 AA 可逆。反过来,若 AA 可逆,则 AA1=1|A||A^{-1}|=1,故 A0|A|\ne0。而 A0    r(A)=n    Ax=0|A|\ne0\iff r(A)=n\iff Ax=0 只有零解,分别来自行列式、秩与齐次方程组的基本结论。

二、可交换矩阵

1. 可交换矩阵的定义

一般而言,矩阵乘法不满足交换律,通常有

ABBA.AB\ne BA.

例如

(A+B)(AB)=A2AB+BAB2,(A+B)(A-B)=A^2-AB+BA-B^2,

只有在 AB=BAAB=BA 时,才能化为

(A+B)(AB)=A2B2.(A+B)(A-B)=A^2-B^2.

AB=BAAB=BA,则称矩阵 A,BA,B 可交换。考研真题不仅可能直接给出 AB=BAAB=BA,也可能给出其他条件,要求反推出 AB=BAAB=BA

2. 可交换矩阵的条件

2.1 零矩阵、数量矩阵与对角矩阵

A,BA,B 均为 nn 阶方阵,下列任一条件都能推出 AB=BAAB=BA

  1. 至少一个矩阵为零矩阵。
  2. 至少一个矩阵为数量矩阵。例如 B=kEB=kE 时,AB=BA=kAAB=BA=kA
  3. 两个矩阵均为对角矩阵。

证明: 零矩阵与任意同阶矩阵相乘仍为零矩阵;数量矩阵 kEkE 与任意同阶矩阵相乘均为 kAkA;两个对角矩阵相乘时,只需逐个相乘主对角元素,交换次序不改变结果。

2.2 一个矩阵是另一个矩阵的函数

B=f(A),B=f(A),

其中 f(A)f(A) 表示 AkA^kA1A^{-1}AA^* 等矩阵的线性组合,则 AABB 可交换。

证明: AAAkA^k 显然可交换;当 AA 可逆时,AA1=A1A=EAA^{-1}=A^{-1}A=E;对于伴随矩阵,无论 AA 是否可逆,都有 AA=AA=AEAA^*=A^*A=|A|E。所以 AkA^kA1A^{-1}AA^* 都与 AA 可交换。与 AA 可交换的矩阵作线性组合后仍与 AA 可交换,因此

Af(A)=f(A)A.Af(A)=f(A)A.

2.3 乘积为非零数量矩阵

AB=kE,k0,AB=kE,\qquad k\ne0,

AB=BAAB=BA

证明:AB=kEAB=kE 可知 A,BA,B 均可逆,并且 B=kA1B=kA^{-1}。因此

BA=kA1A=kE=AB.BA=kA^{-1}A=kE=AB.

2.4 乘积是两个矩阵的线性组合

AB=aA+bB,ab0,AB=aA+bB,\qquad ab\ne0,

AB=BAAB=BA

证明: 将已知条件移项并配成

(AbE)(BaE)=abE,(A-bE)(B-aE)=abE,

乘积为非零数量矩阵,所以两个因子可交换,即

(AbE)(BaE)=(BaE)(AbE).(A-bE)(B-aE)=(B-aE)(A-bE).

展开并约去相同项,得到 AB=BAAB=BA

2.5 二次式与混合乘积的关系

A2+aAB=bE,a0,b0,A^2+aAB=bE, \qquad a\ne0,\quad b\ne0,

AB=BAAB=BA

证明: 原式可写成

A(A+aB)=bE.A(A+aB)=bE.

乘积为非零数量矩阵,所以 AAA+aBA+aB 可交换:

A(A+aB)=(A+aB)A.A(A+aB)=(A+aB)A.

展开并利用 a0a\ne0,得到 AB=BAAB=BA

三、正交矩阵

1. 正交矩阵的定义与等价条件

AAnn 阶实矩阵。若

ATA=E,A^TA=E,

则称 AA 为正交矩阵。nn 阶矩阵 AA 是正交矩阵的充要条件是

AAT=E,A1=AT,ATA=E.AA^T=E, \qquad A^{-1}=A^T, \qquad A^TA=E.

证明:ATA=EA^TA=E 可知 ATA^TAA 的逆矩阵,所以 A1=ATA^{-1}=A^T,进而 AAT=EAA^T=E。反过来,任取其中任一等式,也能得到其余等式。ATA=EA^TA=E 表明 AA 的列向量两两正交且模长为 11AAT=EAA^T=E 表明 AA 的行向量也具有同样性质。

2. 正交矩阵的封闭性

AA 为正交矩阵,则 ATA^TA1A^{-1}AA^*A-A 都是正交矩阵;若 A,BA,B 均为同阶正交矩阵,则 ABAB 为正交矩阵,而 A+BA+B 不一定为正交矩阵。

证明: 因为 A1=ATA^{-1}=A^T,所以 ATA^TA1A^{-1} 实为同一矩阵,并且仍满足正交矩阵的定义。又由

A=AA1,A=±1,A^*=|A|A^{-1},\qquad |A|=\pm1,

可知 A=±A1A^*=\pm A^{-1},故 AA^* 也正交。对于 A-A,有 (A)T(A)=ATA=E(-A)^T(-A)=A^TA=E。对于乘积,

(AB)(AB)T=ABBTAT=AAT=E.(AB)(AB)^T=ABB^TA^T=AA^T=E.

(A+B)T(A+B)(A+B)^T(A+B) 中还含有交叉项,一般不等于 EE

3. 正交矩阵的行列式

ATA=EA^TA=E 取行列式,得到

A2=1,|A|^2=1,

因此

A=1 或 1.|A|=1\text{ 或 }-1.

证明:

ATA=ATA=A2=E=1,|A^TA|=|A^T||A|=|A|^2=|E|=1,

所以正交矩阵的行列式只能为 111-1

4. 正交矩阵的特征值

正交矩阵的特征值模长为 11。其实际特征值只能为 111-1;非实特征值 a+bia+bi 满足

a2+b2=1,a^2+b^2=1,

并且一定成共轭对出现。若正交矩阵的阶数为奇数且行列式小于零,则 1-1 必为其特征值。

证明:Ax=λxAx=\lambda x,其中 x0x\ne0。正交变换保持向量长度,因此

x=Ax=λx=λx,\lVert x\rVert=\lVert Ax\rVert =\lVert\lambda x\rVert =|\lambda|\lVert x\rVert,

从而 λ=1|\lambda|=1。实特征值只能为 111-1;非实特征值满足 a2+b2=1a^2+b^2=1,且因 AA 为实矩阵而成共轭对出现。共轭对的乘积为 11,所以奇数阶正交矩阵若行列式为负,实特征值中必有 1-1

四、实对称矩阵

1. 实对称矩阵的定义

AT=A,A^T=A,

aij=ajia_{ij}=a_{ji},则称 AA 为实对称矩阵。

2. 特征值为实数

实对称矩阵的特征值均为实数,并且可以取实特征向量。

证明:Ax=λxAx=\lambda x,其中 xx 可以是复向量。因为 AA 为实对称矩阵,

xAx=λxxx^*Ax=\lambda x^*x

是实数,而 xx>0x^*x>0,所以 λ\lambda 必为实数。由于 AλEA-\lambda E 为实矩阵,其齐次方程组可以取非零实解,因此可以取实特征向量。

3. 不同特征值的特征向量正交

Aα=λ1α,Aβ=λ2β,λ1λ2.A\alpha=\lambda_1\alpha, \qquad A\beta=\lambda_2\beta, \qquad \lambda_1\ne\lambda_2.

αTβ=0\alpha^T\beta=0

证明: 一方面,

αTAβ=λ2αTβ;\alpha^TA\beta=\lambda_2\alpha^T\beta;

另一方面,由 AT=AA^T=A

αTAβ=(Aα)Tβ=λ1αTβ.\alpha^TA\beta=(A\alpha)^T\beta =\lambda_1\alpha^T\beta.

(λ1λ2)αTβ=0(\lambda_1-\lambda_2)\alpha^T\beta=0。因为 λ1λ2\lambda_1\ne\lambda_2,所以 αTβ=0\alpha^T\beta=0

4. 正交相似对角化

实对称矩阵可用正交矩阵相似对角化:存在正交矩阵 QQ,使

QTAQ=Λ.Q^TAQ=\Lambda.

这里 QQ 是由 AA 的单位正交化特征向量构成的,Λ\Lambda 是由 AA 的特征值构成的对角矩阵,且 QQ 的每一列与 Λ\Lambda 的每一个主对角元素对应。

证明: 实对称矩阵的 Rayleigh 商在单位球面上能取到最大值,取到最大值的向量就是一个实特征向量。这个特征向量的正交补空间在 AA 的作用下保持不变;把 AA 限制在该正交补空间上,仍得到实对称线性变换。对阶数作归纳,可得到 nn 个两两正交的实特征向量。再把每个特征子空间内的基单位正交化,并将所得单位特征向量依次作为列向量组成 QQ,则 QTQ=EQ^TQ=E,并且 AQ=QΛAQ=Q\Lambda。两端左乘 QTQ^T,得到 QTAQ=ΛQ^TAQ=\Lambda

5. 用正交特征向量判定实对称矩阵

nn 阶实矩阵 AA 是实对称矩阵的充要条件是:AAnn 个两两正交的实特征向量。

证明: 必要性由正交相似对角化结论得到。充分性方面,将这 nn 个特征向量单位化并组成正交矩阵 QQ,则

A=QΛQT.A=Q\Lambda Q^T.

由于 ΛT=Λ\Lambda^T=\Lambda,所以

AT=(QΛQT)T=QΛQT=A,A^T=(Q\Lambda Q^T)^T=Q\Lambda Q^T=A,

AA 为实对称矩阵。

五、秩一矩阵

AAnn 阶方阵,若 r(A)=1r(A)=1,则有以下结论。

1. 秩一矩阵的分解形式

r(A)=1    A=αβT,α0, β0.r(A)=1 \iff A=\alpha\beta^T, \qquad \alpha\ne0,\ \beta\ne0.

即:秩一矩阵一定可以写成非零列向量乘以非零行向量的形式,非零列向量乘以非零行向量得到的结果一定是一个秩一矩阵。

证明:r(A)=1r(A)=1,任取 AA 的一个非零列向量 α\alpha,其余每一列都是 α\alpha 的倍数。把这些倍数依次排成行向量 βT\beta^T,便有 A=αβTA=\alpha\beta^T。反过来,αβT\alpha\beta^T 的每一列都是 α\alpha 的倍数,且 α,β\alpha,\beta 均非零,所以其秩恰为 11

2. 秩一矩阵的迹

秩一矩阵的迹等于行向量乘以列向量:

tr(A)=αTβ=βTα.\operatorname{tr}(A)=\alpha^T\beta=\beta^T\alpha.

证明:α=(a1,,an)T\alpha=(a_1,\ldots,a_n)^Tβ=(b1,,bn)T\beta=(b_1,\ldots,b_n)^T,则 A=αβTA=\alpha\beta^T 的主对角元素依次为 a1b1,,anbna_1b_1,\ldots,a_nb_n,所以

tr(A)=i=1naibi=αTβ=βTα.\operatorname{tr}(A)=\sum_{i=1}^na_ib_i =\alpha^T\beta=\beta^T\alpha.

3. 秩一矩阵的高次幂

l=tr(A),l=\operatorname{tr}(A),

A2=αβTαβT=(βTα)A=lA,A^2=\alpha\beta^T\alpha\beta^T =(\beta^T\alpha)A =lA,

因而

An=ln1A.A^n=l^{n-1}A.

证明:A=αβTA=\alpha\beta^TβTα=l\beta^T\alpha=l,直接得到 A2=lAA^2=lA。若 Ak=lk1AA^k=l^{k-1}A,则

Ak+1=AkA=lk1A2=lkA.A^{k+1}=A^kA=l^{k-1}A^2=l^kA.

由数学归纳法可得 An=ln1AA^n=l^{n-1}A

4. 秩一矩阵的特征值

AA 的特征值为

λ1=λ2==λn1=0,λn=tr(A).\lambda_1=\lambda_2=\cdots=\lambda_{n-1}=0, \qquad \lambda_n=\operatorname{tr}(A).

证明: 因为 r(A)=1r(A)=1,所以 Ax=0Ax=0 的基础解系含有 n1n-1 个线性无关向量,故 00 至少是 n1n-1 重特征值。又因为全部特征值之和等于矩阵的迹,剩余一个特征值只能是 tr(A)\operatorname{tr}(A)

5. 秩一矩阵相似对角化的条件

AA 能相似对角化的充要条件是

A 能相似对角化    tr(A)0.A\text{ 能相似对角化} \iff \operatorname{tr}(A)\ne0.

tr(A)0\operatorname{tr}(A)\ne0,则

Adiag(0,,0,tr(A)).A\sim \operatorname{diag}(0,\ldots,0,\operatorname{tr}(A)).

tr(A)=0\operatorname{tr}(A)=0,则 A2=OA^2=O,且 AA 不能相似对角化。

证明:tr(A)0\operatorname{tr}(A)\ne0,则 00 对应的特征子空间维数为 nr(A)=n1n-r(A)=n-1,而 tr(A)\operatorname{tr}(A) 还对应一个与它们线性无关的特征向量,所以 AAnn 个线性无关的特征向量,可以相似对角化。若 tr(A)=0\operatorname{tr}(A)=0,则所有特征值均为 00。假如 AA 可以相似对角化,就只能相似于零矩阵,从而 A=OA=O,这与 r(A)=1r(A)=1 矛盾。

6. 秩一矩阵的特征向量

tr(A)=0\operatorname{tr}(A)=0AA 的特征值全为 00,则只需解

Ax=0,Ax=0,

即可得到特征值 00 对应的特征向量。

tr(A)0\operatorname{tr}(A)\ne0,则

λ1=λ2==λn1=0,λn=tr(A)=αTβ=βTα.\lambda_1=\lambda_2=\cdots=\lambda_{n-1}=0, \qquad \lambda_n=\operatorname{tr}(A)=\alpha^T\beta=\beta^T\alpha.

特征值 00 对应的特征向量仍只需解 Ax=0Ax=0

证明: 特征值 00 对应的特征向量按定义就是齐次方程组 Ax=0Ax=0 的非零解。又因为

Aα=αβTα=(βTα)α=tr(A)α,A\alpha=\alpha\beta^T\alpha =(\beta^T\alpha)\alpha =\operatorname{tr}(A)\alpha,

所以 α\alpha 是特征值

λn=tr(A)=αTβ=βTα\lambda_n=\operatorname{tr}(A)=\alpha^T\beta=\beta^T\alpha

对应的特征向量。

7. 数量矩阵与秩一矩阵之和

若矩阵

B=A+kE,B=A+kE,

则可利用秩一矩阵的结论,快速计算矩阵 BB 的特征值;BB 的特征值等于 AA 的特征值全部加上 kk

证明:Ax=λxAx=\lambda x,则

Bx=(A+kE)x=(λ+k)x,Bx=(A+kE)x=(\lambda+k)x,

所以 BB 的特征值是在 AA 的每个特征值上加 kk

8. 两个秩一矩阵相加

见到两个秩一矩阵相加,条件反射想到

r(A+B)r(A)+r(B).r(A+B)\le r(A)+r(B).

证明: A+BA+B 的列空间包含于 AABB 的列空间之和,因此其维数不超过两个列空间维数之和。特别地,当 r(A)=r(B)=1r(A)=r(B)=1 时,r(A+B)2r(A+B)\le2

六、二次多项式矩阵

AAnn 阶方阵,且

(AaE)(AbE)=O,ab.(A-aE)(A-bE)=O, \qquad a\ne b.

这是一个重要的秩推导过程,可依次得到以下结论。

1. 特征值只能为 aabb

证明:

(AaE)(AbE)=O(A-aE)(A-bE)=O

可知特征值满足

(λa)(λb)=0,(\lambda-a)(\lambda-b)=0,

所以

λ=aλ=b.\lambda=a\quad\text{或}\quad\lambda=b.

即矩阵 AA 的特征值只能为 aabb

2. 两个矩阵的秩之和为 nn

证明:

(AaE)(AbE)=O(A-aE)(A-bE)=O

可得

r(AaE)+r(AbE)n.r(A-aE)+r(A-bE)\le n.

另一方面,

r(AaE)+r(AbE)=r(aEA)+r(AbE)r(aEA+AbE)=r((ab)E)=n.\begin{aligned} r(A-aE)+r(A-bE) &=r(aE-A)+r(A-bE)\\ &\ge r(aE-A+A-bE)\\ &=r((a-b)E)=n. \end{aligned}

因此

r(AaE)+r(AbE)=n.r(A-aE)+r(A-bE)=n.

3. 特征值 aa 对应的特征向量

证明: 仍由

(AaE)(AbE)=O(A-aE)(A-bE)=O

可知,(AbE)(A-bE) 的每一个列向量都是齐次方程组

(AaE)x=0(A-aE)x=0

的解。因此,(AbE)(A-bE) 的非零列向量都是特征值 λ=a\lambda=a 对应的特征向量,从而 λ=a\lambda=a 对应的线性无关特征向量的个数至少为

r(AbE).r(A-bE).

4. 特征值 bb 对应的特征向量

证明: 展开两个因子的乘积,有

(AaE)(AbE)=A2(a+b)A+abE=(AbE)(AaE)=O.\begin{aligned} (A-aE)(A-bE) &=A^2-(a+b)A+abE\\ &=(A-bE)(A-aE)=O. \end{aligned}

所以,(AaE)(A-aE) 的每一个列向量都是齐次方程组

(AbE)x=0(A-bE)x=0

的解。因此,(AaE)(A-aE) 的非零列向量都是特征值 λ=b\lambda=b 对应的特征向量,从而 λ=b\lambda=b 对应的线性无关特征向量的个数至少为

r(AaE).r(A-aE).

5. 矩阵 AA 可以相似对角化

证明: 综合前面的结论,矩阵 AA 的线性无关特征向量至少有

r(AbE)+r(AaE)=nr(A-bE)+r(A-aE)=n

个。又因为 nn 阶矩阵最多有 nn 个线性无关的特征向量,所以 AA 恰有 nn 个线性无关的特征向量,因而 AA 可以相似对角化。

6. 常见二次关系

幂等矩阵和满足二次关系的矩阵都可以按照这一思路处理。例如:

A2=AA(AE)=O(A0E)(AE)=O,A^2=A \Rightarrow A(A-E)=O \Rightarrow (A-0E)(A-E)=O,

所以 AA 可以相似对角化;

A2=2E(A+2E)(A2E)=O,A^2=2E \Rightarrow (A+\sqrt2E)(A-\sqrt2E)=O,

所以 AA 可以相似对角化;

A2+2A=3E(A+3E)(AE)=O,A^2+2A=3E \Rightarrow (A+3E)(A-E)=O,

所以 AA 也可以相似对角化。

证明方法: 只要把矩阵满足的二次关系分解为

(AaE)(AbE)=O,ab,(A-aE)(A-bE)=O,\qquad a\ne b,

就可以直接套用前面的证明。两个一次因子对应两个不同的特征值,因此矩阵可以相似对角化。

七、列满秩矩阵

设矩阵 Am×nA_{m\times n} 列满秩,即

r(A)=n,r(A)=n,

此时隐含条件为 mnm\ge n。列满秩矩阵有以下重要结论。

r(A)=n    Ar(A)=n\iff A 可经有限次初等行变换化为 (EnO)\begin{pmatrix}E_n\\O\end{pmatrix}

r(A)=nr(A)=n 当且仅当 AA 可经过有限次初等行变换化为

(EnO),\begin{pmatrix}E_n\\O\end{pmatrix},

即存在可逆矩阵 PP,使得

PA=(EnO).PA=\begin{pmatrix}E_n\\O\end{pmatrix}.

证明: 初等行变换不改变矩阵的列秩。r(A)=nr(A)=n 表明 AA 的每一列都有主元,所以其行最简形必为 (EnO)\begin{pmatrix}E_n\\O\end{pmatrix}。有限次初等行变换等价于左乘可逆矩阵 PP,因此得到上述矩阵等式;反向推导同理。

记住这一结论的最好方式是举例并手推。例如

A=(142536),r(A)=2,A=\begin{pmatrix}1&4\\2&5\\3&6\end{pmatrix}, \qquad r(A)=2,

AA 作初等行变换,可以将其化为

(100100).\begin{pmatrix}1&0\\0&1\\0&0\end{pmatrix}.

由此还可得到一个推论:若矩阵 AABB 的列数相同,且均列满秩,则矩阵 AA 和矩阵 BB 的行向量组等价。因为二者都可以通过初等行变换化为

(EnO),\begin{pmatrix}E_n\\O\end{pmatrix},

以它作为桥梁即可。

r(A)=n    Ax=0r(A)=n\iff Ax=0 只有零解

r(A)=nr(A)=n 当且仅当齐次方程组

Ax=0Ax=0

只有零解。其逆否命题是:若 x0x\ne0,则 Ax0Ax\ne0。因此还可以写成

r(A)=n    对于任意 α0, 均有 Aα0.r(A)=n \iff \text{对于任意 }\alpha\ne0,\text{ 均有 }A\alpha\ne0.

证明: r(A)=nr(A)=n 表明齐次方程组的未知数个数等于秩,自由未知数个数为零,因此只有零解。反过来,若 Ax=0Ax=0 只有零解,则 nr(A)=0n-r(A)=0,所以 r(A)=nr(A)=n

r(A)=n    r(A)=n\iffAB=OAB=O,则 B=OB=O

r(A)=nr(A)=n 当且仅当:对任意适配的矩阵 BB,若 AB=OAB=O,则 B=OB=O

证明: 把矩阵 BB 按列分块,每一列都满足 Ax=0Ax=0。因为该齐次方程组只有零解,所以 BB 的每一列都是零向量,即 B=OB=O

r(A)=n    r(A)=n\iffAB=ACAB=AC,则 B=CB=C

r(A)=nr(A)=n 当且仅当:对任意适配的矩阵 B,CB,C,若 AB=ACAB=AC,则 B=CB=C

证明: 这是列满秩矩阵的“左消去律”,因为

A(BC)=OBC=O.A(B-C)=O\Rightarrow B-C=O.

r(A)=n    ATAr(A)=n\iff A^TA 为正定矩阵

r(A)=nr(A)=n 当且仅当 ATAA^TA 为正定矩阵,即

r(A)=n    ATA 为正定矩阵.r(A)=n\iff A^TA\text{ 为正定矩阵}.

先说明 ATAA^TA 是对称矩阵:

(ATA)T=ATA.(A^TA)^T=A^TA.

证明: 必要性:若 r(A)=nr(A)=n,则对任意 x0x\ne0,都有 Ax0Ax\ne0,从而

xTATAx=(Ax)T(Ax)>0.x^TA^TAx=(Ax)^T(Ax)>0.

由正定矩阵的定义可知,ATAA^TA 正定。

充分性:若 ATAA^TA 正定,则对任意 x0x\ne0,都有

xTATAx>0,x^TA^TAx>0,

(Ax)T(Ax)>0,(Ax)^T(Ax)>0,

所以 Ax0Ax\ne0。由前面的逆否命题可知 r(A)=nr(A)=n,充分性得证。

r(A)=n    r(A)=n\iff 存在矩阵 BB,使得 BA=EBA=E

r(A)=nr(A)=n 当且仅当存在矩阵 BB,使得

BA=E.BA=E.

等价地,存在矩阵 BB,使得

ATBT=E.A^TB^T=E.

证明: 必要性:若 AA 列满秩,则 ATA^T 行满秩。根据行满秩矩阵的基本结论,方程 ATx=bA^Tx=b 一定有解,因此矩阵方程

ATX=EA^TX=E

也一定有解。将 XX 按列分块即可知,存在矩阵 BTB^T,使得

ATBT=E,A^TB^T=E,

即存在矩阵 BB,使得 BA=EBA=E

充分性:若 BA=EBA=E,则由矩阵秩的强弱关系可知

r(A)r(BA)=r(E)=n.r(A)\ge r(BA)=r(E)=n.

AA 只有 nn 列,所以 r(A)=nr(A)=n,即 AA 列满秩。充分性得证。

r(A)=nBx=0r(A)=n\Rightarrow Bx=0ABx=0ABx=0 同解

r(A)=nr(A)=n,则方程组

Bx=0Bx=0

ABx=0ABx=0

同解。

证明:Bx=0Bx=0,则等式两端左乘 AA,可得 ABx=0ABx=0,所以 Bx=0Bx=0 的解都是 ABx=0ABx=0 的解。反过来,若 ABx=0ABx=0,因为 r(A)=nr(A)=n,齐次方程组 Ay=0Ay=0 只有零解,令 y=Bxy=Bx,便有

A(Bx)=0Bx=0.A(Bx)=0\Rightarrow Bx=0.

所以 ABx=0ABx=0 的解也是 Bx=0Bx=0 的解,二者同解。

方程组同解的充要条件是行向量组等价,因此 BB 的行向量组与 ABAB 的行向量组等价,从而

r(B)=r(AB)=r(BAB).r(B)=r(AB)=r\begin{pmatrix}B\\AB\end{pmatrix}.

常用的“左乘列满秩,秩不变”,就是上述结论的一部分:

r(A)=nr(B)=r(AB).r(A)=n\Rightarrow r(B)=r(AB).

八、行满秩矩阵

设矩阵 Am×nA_{m\times n} 行满秩,即

r(A)=m,r(A)=m,

隐含条件为 mnm\le n。行满秩矩阵独有的性质并不多,最常用的是:若 AA 行满秩,则 Ax=bAx=b 一定有解。

① 若 AA 行满秩,则 Ax=bAx=b 一定有解,AX=BAX=B 也一定有解

Ax=bAx=b

一定有解;矩阵方程

AX=BAX=B

也一定有解。

证明:r(A)=mr(A)=m,则

m=r(A)r(A,b)m,m=r(A)\le r(A,b)\le m,

所以 r(A)=r(A,b)=mr(A)=r(A,b)=m,即 Ax=bAx=b 一定有解。将矩阵方程 AX=BAX=B 中的 XXBB 按列分块,XX 的每一列都对应一个有解的非齐次方程组,因此 AX=BAX=B 一定有解。

r(A)=m    Ar(A)=m\iff A 可经有限次初等列变换化为 (Em,O)(E_m,O)

即存在可逆矩阵 PP,使得

AP=(Em,O).AP=(E_m,O).

证明: r(A)=mr(A)=m 表明 AA 的每一行都有主元。初等列变换不改变行秩,因此其列最简形为 (Em,O)(E_m,O);有限次初等列变换等价于右乘可逆矩阵 PP。反过来,若 AP=(Em,O)AP=(E_m,O)PP 可逆,则初等列变换不改变秩,所以 r(A)=mr(A)=m

行满秩矩阵其余的更多性质,基本上都是因为“AA 行满秩等价于 ATA^T 列满秩”,然后利用上一节列满秩矩阵的相关结论进行推导。

r(A)=mr(BA)=r(B)r(A)=m\Rightarrow r(BA)=r(B),即右乘行满秩,秩不变

r(A)=mr(BA)=r(B),r(A)=m\Rightarrow r(BA)=r(B),

证明: ATA^T 列满秩,由“左乘列满秩,秩不变”可知

r(BA)=r(ATBT)=r(BT)=r(B).r(BA)=r(A^TB^T)=r(B^T)=r(B).

r(A)=m    r(A)=m\iffBA=CABA=CA,则 B=CB=C

BA=CAB=C.BA=CA\Rightarrow B=C.

证明:BA=CABA=CA 两端同时取转置,得到

ATBT=ATCT.A^TB^T=A^TC^T.

因为 ATA^T 列满秩,由列满秩矩阵的左消去律可得 BT=CTB^T=C^T,即 B=CB=C

r(A)=m    AATr(A)=m\iff AA^T 为正定矩阵

r(A)=m    AAT 为正定矩阵.r(A)=m\iff AA^T\text{ 为正定矩阵}.

证明: 因为 r(A)=m    r(AT)=mr(A)=m\iff r(A^T)=m,可把列满秩矩阵的正定判定应用于 ATA^T

r(AT)=m    (AT)TAT=AAT 为正定矩阵.r(A^T)=m \iff (A^T)^TA^T=AA^T\text{ 为正定矩阵}.

其余更多性质也可以由 ATA^T 列满秩推出。

九、矩阵的高次幂

计算矩阵的 nn 次幂,应优先识别结构。常用方法如下。

1. 秩一矩阵的二级结论

优先看是否为秩一矩阵。若 A=αβTA=\alpha\beta^T,令 l=tr(A)l=\operatorname{tr}(A),则

An=ln1A.A^n=l^{n-1}A.

证明: 因为 A2=(βTα)A=tr(A)A=lAA^2=(\beta^T\alpha)A=\operatorname{tr}(A)A=lA,连续左乘 AA 即得 An=ln1AA^n=l^{n-1}A

2. 二项式展开

必须满足可交换。一般写成

A=B+kE,A=B+kE,

其中 BB 在进行几次乘法运算后为零矩阵,例如严格上三角矩阵或严格下三角矩阵。若 Bs=OB^s=O,则

(kE+B)n=j=0s1(nj)knjBj.(kE+B)^n =\sum_{j=0}^{s-1}\binom{n}{j}k^{n-j}B^j.

证明: 数量矩阵 kEkEBB 可交换,所以可以使用二项式定理。原展开式本应从 j=0j=0 写到 nn;由于 Bs=OB^s=O,所有 jsj\ge s 的项均为零,只需保留到 s1s-1 项。

严格三角矩阵的非零元素都在对角线下方或上方,非零元素构成的三角形边长为 kk 时,常有

BkO,Bk+1=O.B^k\ne O, \qquad B^{k+1}=O.

这是因为严格三角矩阵每乘一次,可能出现非零元素的位置就向远离主对角线的方向移动一层;移动超过矩阵的有效三角区域后,所有元素都变为零。

3. 分块矩阵

分块矩阵可结合前面几种方法分别处理各个对角块;块对角结构满足

(BOOC)n=(BnOOCn).\begin{pmatrix}B&O\\O&C\end{pmatrix}^{n} =\begin{pmatrix}B^n&O\\O&C^n\end{pmatrix}.

证明: 两个块对角矩阵相乘时,非对角块始终为零,对角块分别相乘。对幂次作数学归纳即可得到上述公式。

4. 相似对角化

这是通法,也是最重要的方法。若

P1AP=Λ,P^{-1}AP=\Lambda,

An=PΛnP1.A^n=P\Lambda^nP^{-1}.

高次幂问题通常将矩阵对角化,然后利用对角矩阵的幂容易计算这一点处理。

证明:P1AP=ΛP^{-1}AP=\LambdaA=PΛP1A=P\Lambda P^{-1},所以

An=(PΛP1)n=PΛnP1,A^n =(P\Lambda P^{-1})^n =P\Lambda^nP^{-1},

中间相邻的 P1PP^{-1}P 全部约为单位矩阵。

5. 初等变换与找规律

若矩阵由简单的初等变换得到,或无法相似对角化,可从初等变换的角度理解它的作用;也可以先计算 A2,A3A^2,A^3 寻找递推规律。找规律可以使用,但仍建议优先以结构和相似对角化为主。

依据: 若计算出较低次幂后出现 Am=c0E+c1A++cm1Am1A^m=c_0E+c_1A+\cdots+c_{m-1}A^{m-1},后续高次幂都可以用这个关系递推降次。其本质仍是寻找矩阵满足的多项式关系;若能进一步分解或对角化,应优先改用前面的结构化方法。

Back to archive

Discussion

Comments

Post

Share questions, corrections, or extra notes about this post.