前两章把矩阵从「数表」讲成了「线性映射的坐标版」:取定基,映射与矩阵一一对应(2.3),而且同一个映射换一组基就换一个矩阵(B=P−1AQ,定理 2.3.5)。于是「矩阵有多少种样子」这个问题变成了「同一个映射能变成多少种样子」——本章换一个方向问:能不能把矩阵拆开,拆成几个结构简单、一看就懂的矩阵的乘积或和?这就是矩阵分解。
分解要从最简单的一种开始。任给一个矩阵,无论它多大多怪,总有一件事一定能做到:用可逆矩阵从两边夹它,把它压成一个「前面一块单位矩阵、其余全零」的标准形——相抵标准形 Γ。这条结论写下来只有一行(引理 3.1.1),却把「秩」这个数彻底说清了:任何矩阵都相抵于且只相抵于由秩决定的标准形(本编补注)。本节随后把这条分解用到两个地方,都是考试与后续章节反复出现的老朋友:解线性方程组(例 3.1.1:相抵标准形把「有解没解、有几个自由度」变成看几个分量是否为零)与长方阵的单侧逆(定理 3.1.1:行满秩就有右逆、列满秩就有左逆)。最后给出本节最重要的一条秩等式 rank(A)=rank(AHA)(命题 3.1.1),它给出的右逆公式 AH(AAH)−1 是「最小二乘」等一大类应用的起点。
素材与说明(本节无录播)
本页依据课程第三章课件 3.1 相抵分解(15 页)整理。本节无对应录播,故不引用课堂口述与时间戳:课件给出的引理、定理、命题、注记与例题照录并保留原编号(引理 3.1.1、例 3.1.1、定理 3.1.1、注 1、命题 3.1.1);课件省略的证明、只写一半的推导、未给出答案的练习,由本编补出并逐处标注「本编补」。编号与表述问题见勘误与补注。
图解读图:任何矩阵都能被两个可逆矩阵「夹」成一个只由秩决定的标准形 Γ。左乘 P⁻¹ 换到达端的基,右乘 Q 换出发端的基——矩阵在换基,映射本身没变。这也解释了为什么相抵标准形唯一、分解式不唯一。
先交代「分解」这件事本身。课件开篇只用了三句话和一个类比的引子:把矩阵分解为较为简单或具有某种特性的若干矩阵的乘积或和,就像 36=4×9=22×32——同一个数换几种乘法写法,每种写法都能显出它的部分性质;分解的意义在于清晰地反映出原矩阵的某些特征,提供有效的理论分析或数值计算依据。
那么任给一个矩阵,能进行何种分解?最朴素、对任何矩阵都可行的答案是:从左右两边乘可逆矩阵,把它化成只有一个参数(秩)的标准形。这就是矩阵相抵。
引理 3.1.1 · 相抵分解式与相抵标准形
设 A∈Fm×n 的秩为 r,则存在两个可逆矩阵 P∈Fm×m 和 Q∈Fn×n 使得
A=P−1ΓQ,Γ=[Ir000]m×n.
上式称为 A 的相抵分解式,式中 Γ 称为 A 的相抵标准形。
课件把这条结论列为引理、只述不证——它正是 2.3 留下的一句承诺(「相抵的标准形只由秩决定,将在讨论矩阵的秩时给出」,见 2.3.4 的注 1)。本编用「选基 + 换基公式」把证明补出;你会看到它把前三章各一件工具串在了一起。
证明(课件未给,本编补:选一组好基,再用换基公式)
把 A 读成映射。取 Fn 与 Fm 的标准基,设线性映射 T:Fn→Fm 在这两组标准基下的矩阵为 A。(A 本来就是某个映射的矩阵,2.3 的 定义 2.3.1 读法。)
选出发端的基。由亏加秩定理(定理 2.2.5),dimN(T)=n−r,取 N(T) 的一组基 ξ1,⋯,ξn−r,并用基扩充定理(推论 1.3.1)扩充为 Fn 的一组基
ξ1,⋯,ξn−r, ξn−r+1,⋯,ξn.
像的后 r 个基向量独立。记 fi=T(ξn−r+i)(i=1,⋯,r)。若 ∑i=1rkifi=0,即 T(∑ikiξn−r+i)=0,则 ∑ikiξn−r+i∈N(T)=span{ξ1,⋯,ξn−r};而这两个向量组同属一组基,由表示唯一性(定理 1.3.2)得全部系数为零。故 f1,⋯,fr 线性无关;又任一 T(x) 只是 ξn−r+1,⋯,ξn 的像的组合(ξ1,⋯,ξn−r 被映成零),故 f1,⋯,fr 张成 R(T)——它们是 R(T) 的一组基,个数 r=dimR(T)。
选到达端的基。把 f1,⋯,fr 扩充为 Fm 的一组基
η1=f1, ⋯, ηr=fr, ηr+1,⋯,ηm.
读出矩阵。在基 ξ1,⋯,ξn 与 η1,⋯,ηm 下,映射 T 的矩阵恰为 Γ:
T(ξ1)=⋯=T(ξn−r)=0,T(ξn−r+i)=ηi(i=1,⋯,r),
前 n−r 列全零、后 r 列依次是 e1,⋯,er,按 定义 2.3.1 逐列对照正是 Γ。
用换基公式收尾。设由标准基到 ξ 的过渡矩阵为 S、由标准基到 η 的过渡矩阵为 R。换基公式(定理 2.3.5)说的是:旧基下矩阵 A 与新基下矩阵 Γ 满足 Γ=R−1AS,即
A=RΓS−1.
取 P=R−1∈Fm×m、Q=S−1∈Fn×n(过渡矩阵可逆,命题 1.3.1),两者都可逆,得 A=P−1ΓQ。∎
证明的要点:Γ 之所以可能出现,是因为「秩」把坐标的 n 个自由度分成了两半——n−r 个被映射压死(进零空间)、r 个原样传下去(映成值域的基)。把这两半分别选成出发端与到达端的基,映射就只剩下「前 n−r 列归零、后 r 列照抄」这一副骨架。
∎
本编补注 · 「相抵」这个名字与它刚兑现的承诺
相抵(等价) 的定义在 2.3.4 的注 1:A 经有限次初等变换变成 B,即存在可逆 P,Q 使 B=P−1AQ,记 A≅B(矩阵相抵,与空间的同构 V≅W 同名不同义)。相抵是等价关系,而 引理 3.1.1 把每一类都选出了一个代表元 Γ,所以:
两个同型矩阵相抵 ⟺ 它们的秩相等。(本编补:⟹ 用初等变换不改变秩;⟸ 设秩同为 r,由 引理 3.1.1 两者都相抵于同一个 Γ,再用传递性。)
这正是课件那句「相抵分解适于任何矩阵,常解决矩阵秩相关的命题」的准确含义:关于秩的命题都可以先化到 Γ 上验证,再搬回来。后续几节会反复使用这一套路(例如 3.2 把满秩分解的第一种构造直接搭在相抵分解上)。
两个容易混的同名概念再区分一次:相抵允许左右用各自独立的两组基(P,Q 无关),相似(B=P−1AP,推论 2.3.4)只允许同一空间换一组基、两侧用同一个 P——相似必相抵,相抵未必相似(只有同阶方阵才谈相似)。
分解唯一性问题(课件原文)
课件紧接着问:分解式唯一吗?答案是不唯一。设 Γ=[Ir000],对任意适当阶数的矩阵 ∗,有
A=P−1ΓQ=P−1[Ir0∗I]Γ[Ir∗0I]Q.
为什么(本编补上这一步的验算):夹在中间的乘法之所以能把 Γ 换出来,是因为
[Ir0∗I]Γ=Γ=Γ[Ir∗0I]
(直接按分块乘法核对:Γ 的后 n−r 列全零、后 m−r 行全零,∗ 乘上去也归零)。而 [Ir0∗I] 与 [Ir∗0I] 都是可逆矩阵(对角块为单位阵的三角块,行列式为 1),乘到 P−1、Q 上仍是可逆矩阵。因此矩阵 A 的相抵分解是不唯一的——Γ 两边的可逆矩阵可以有无穷多种取法。
不唯一的是「分解式」,唯一的是「标准形」:无论怎么写分解式,中间那个 Γ 只由 r 决定(因为秩在相抵下不变,rankΓ=rankA=r)。这一对「式不唯一、标准形唯一」后面还会以同样的面貌出现(3.2 的满秩分解不唯一,但相差一个可逆矩阵 D)。
课件给相抵分解配的第一个应用,是它最经典的用武之地——解方程组 Ax=b。道理是:相抵变换是「换基」,换基不改变方程的解集(它只是用可逆矩阵 P、Q 把等号两边同时翻译一遍)。翻译之后,方程组的系数矩阵变成 Γ,于是「有解没解」「有几个自由参数」都变成数分量的问题。
例 3.1.1(课件原文)· 考查线性方程组 Ax=b 的求解问题
设 A∈Rm×n、b∈Rm 已知,x∈Rn 是未知向量。将矩阵 A 作相抵分解,讨论方程组的解。
解答(课件给出前两步,本编补完解的结构与可解性判据)
第一步(换元):把 A=P−1ΓQ 代入,得
P−1ΓQx=b.
两边左乘 P(可逆):
Γ(Qx)=Pb.
令 y=Qx(这是一次可逆的未知量替换:Q 可逆,y↔x 是一一对应),上式改写为
Γy=Pb,Γ=[Ir000]m×n.
第二步(读出全部信息):把 Pb 按前 r 个、后 m−r 个分量分块,写成 Pb=[cd],c∈Fr、d∈Fm−r;把 y 按前 r 个、后 n−r 个分量分块为 y=[y1y2]。方程 Γy=Pb 逐分量读出来是
y1=c(前 r 个方程),0=d(后 m−r 个方程).
于是结论有三条(本编补全):
- 可解判据:方程组有解 ⟺d=0,即 Pb 的后 m−r 个分量全为零。这是「rank(A∣b)=rankA」的坐标版:那个条件说的正是 b 落进 A 的列空间,而这里的 d=0 就是它的显式算法形式。
- 有解时的通解:y1 被唯一确定(=c),y2∈Fn−r 任意;再乘回 x=Q−1y 得
x=Q−1[cy2],y2 任取,
即解集是 Fn 中一个 n−r 维的仿射子空间——自由参数的个数恰是 n−r,与亏加秩定理(齐次解空间维数 =n−r)完全一致。
3. 无解时的信息:d=0 时,越界的分量就是「方程之间的矛盾」,m−r 个多余方程里只要有矛盾,整个方程组就无解。
两处记号提醒:① 解出 y 之后必须乘回 Q−1 才是原未知量 x(与 2.4 求特征向量时「坐标要乘回基」是同一个陷阱);② P 与 Q 的角色不对称:P 管等号左边的行(它把 A 的行压到标准位置),Q 管未知量的换元。∎
剩下一个绕不开的实际问题:P、Q 怎么求?课件给出的答案是一张分块矩阵——把两个可逆矩阵「挂」在 A 的两侧,与 A 同时接受变换,A 化到 Γ 时它们也自动积好了。本编把它的记号与 引理 3.1.1 统一(课件原页的记号见紧随的勘误框)。
相抵变换求 P 与 Q−1(本编整理,记号与引理 3.1.1 统一)
构造分块矩阵(∗ 是任意合适阶数的块,只作占位)
M=[AInIm∗].
操作规则:只对前 m 行做行变换、只对前 n 列做列变换,直到左上块变成 Γ。此时 M 自动变成
M ⇝ [ΓQ−1P∗].
为什么(一行验算):把所做的行变换总积记为 E、列变换总积记为 F,则整套操作就是左乘 [E00I]、右乘 [F00I]:
[E00I][AInIm∗][F00I]=[EAFFE∗].
左上块 EAF=Γ 说明 A=E−1ΓF−1;与 A=P−1ΓQ 对照,右上块 E 就是 P、左下块 F 就是 Q−1(等价地:P=E、Q=F−1,而 P−1=E−1)。
记忆:A 左边的行变换记录在右上,右边的列变换记录在左下——「变换乘在哪侧,就积在哪个对角块的对面」。算完务必回代验算 EAF=Γ(比验 P−1ΓQ 更省事,也不容易弄错栏目)。
课件第 7 页的记号与引理 3.1.1 不同调(本编勘误 · E3.1)
课件第 7 页的分解式写作 A=P−1ΓQ(引理 3.1.1),但紧接着的分块矩阵一页把读出的两块标成 P−1 与 Q−1:按该页自己的式子 [P−100I][AII∗][Q−100I]=Ir0Q−100P−1∗ 展开,得到的是 P−1AQ−1=Γ,即 A=PΓQ——与引理的 A=P−1ΓQ 差一个位置。两处只要取同一组字母就对不上。
本编处理:正文的 相抵变换求 P 与 Q−1 已改成自洽写法(读出的两块是真 P=E 与 Q−1=F,且 EAF=Γ 可直接验算);引用本方法时请以正文为准。这类「不同页对同一字母取不同角色」的记号循环在 3.1 内还有一处(右逆充分性证明法 1 的 A=P[Im,0]Q 与引理的 A=P−1ΓQ,本编正文已按引理记号统一),在 3.2 也出现过(同样按引理记号调顺),不影响结论,只影响抄写。
下面把整套流程在一个非平凡例子上走完。例题的每一步都做过回代验算,你可以照着复现。
例 3.1.A(新增)· 相抵分解全流程:一个秩 2 的三阶方阵
求 A=1232573710 的相抵分解式。
解答 · 例 3.1.A
先定秩:第三行 = 第一行 + 第二行,且前两行不成比例,故 r=2。
挂上分块矩阵(右下块只占位,记作 0):
M=1231002570103710001100010001.
行变换(只动前 3 行):r3←r3−r1−r2(第三行清零),r2←r2−2r1:
1001002100103100011−2−101−1001.
列变换(只动前 3 列):c2←c2−2c1,c3←c3−3c1,再 c3←c3−c2:
M ⇝ 100100010−210000−1−111−2−101−1001.
读出结果(按 相抵变换求 P 与 Q−1:左上块 =Γ、右上块 =E=P、左下块 =F=Q−1):
Γ=100010000,E=P=1−2−101−1001,F=Q−1=100−210−1−11.
回代验算(用 EAF=Γ 最省事):
EAF=1−2−101−10011232573710100−210−1−11=100210310100−210−1−11=100010000=Γ ✓
把分解式写成 A=P−1ΓQ:P 已在手上(右上块),只需把左下块求逆得 Q=F−1,
P=1−2−101−1001,P−1=123011001,Q=F−1=100210311,
于是 A=P−1ΓQ 即所求的相抵分解式。顺手取出的两块又给出了 A 的「骨架」:P−1 的前 r 列 123011 是 R(A) 的一组基,Q 的前 r 行 [102131] 是各列在这组基下的坐标——两者相乘恰好还原 A。这不是巧合,3.2 的满秩分解(方法 1)正建立在它上面。
例 3.1.B(新增 · 承接例 3.1.A)· 用相抵标准形解方程组
设 A 同 例 3.1.A。(a) 解方程组 Ax=b1,b1=(1,2,3)T;(b) 判断 Ax=b2(b2=(1,2,4)T)是否有解。
解答 · 例 3.1.B(新增
用 例 3.1.A 的 P、Γ、Q:方程 Γy=Pb,x=Q−1y。
(a) Pb1=1−2−101−1001123=100。后 m−r=1 个分量为 0 ✓ 有解。由 Γy=Pb1 读得 y1=1、y2=0、y3 任意,记 y3=t:
x=Q−1y=100−210−1−1110t=1−t−tt=100+t−1−11.
验算:A(1,0,0)T=(1,2,3)T=b1 ✓(原方程组的一个特解);A(−1,−1,1)T=0 ✓(齐次解)。通解 = 特解 + 齐次通解,自由参数个数 1=n−r ✓。
(b) Pb2=101,后 1 个分量 =1=0,故无解。对照经典判据:rank(A∣b2)=3>2=rankA ✓(第三个方程 = 前两个之和,而右端 4=1+2=3,矛盾)。
要点:相抵标准形把「解方程组」化成了「看 Pb 的后 m−r 个分量是否为零」——可解性是一个分量计数问题,自由度数是一个维数计数问题。∎
相抵分解与高斯消元是什么关系?(本编)
高斯消元(行变换)把 A 化成行阶梯形,用于判断解与求解;相抵分解则允许同时使用行变换与列变换(等价于左右各乘可逆矩阵),于是能把 A 压得更彻底——一路压到只由秩决定的标准形 Γ。两者的关系可以这样说:
- 行变换(左乘 P)在不改变未知量的前提下化简方程,是「求解」的工具;
- 列变换(右乘 Q)等价于对未知量换元(x=Q−1y),是「换坐标」的工具;
- 行阶梯形只把 A 化到「阶梯」为止(保留了解的信息),相抵标准形则把列方向的冗余也清掉(牺牲了直接可解性,换来了唯一的代表元)。
用哪个取决于目的:要解方程,行变换就够;要研究「秩」本身、或要构造后续的分解,标准形更方便。
第二个应用从一个看似天真的问题开始:什么样的矩阵乘积能为单位阵? 换成等式就是问 AB=I 能对哪些 A 成立——若是方阵,AB=I 已经足够推出 B=A−1、两者都可逆;但若 A 是长方阵,结论完全不同。
问题 1 与问题 2(课件原文)
问题 1:何种矩阵乘积为单位阵?初看几乎要脱口而出「乘积为单位阵的两矩阵必可逆」——错。反例:
AAH=[1 0][10]=1.
这里 A=[1 0] 不是可逆矩阵(连方阵都不是),却与 AH 乘出了单位阵(1×1 的单位阵)。
小结(修正后的正确说法):若两方阵乘积为单位阵,则这两矩阵必可逆。前提「方阵」一个字都不能少。
问题 2:何种长方阵乘积为单位阵?——这就是下文定理 3.1.1 回答的问题。
定理 3.1.1 · 右逆与左逆
设矩阵 A∈Crm×n(r>0),则
- A 有右逆(即存在矩阵 B 使得 AB=I)的充分必要条件是 A 为行满秩矩阵(r=m);
- A 有左逆(即存在矩阵 B 使得 BA=I)的充分必要条件是 A 为列满秩矩阵(r=n)。
记号说明(本编):Crm×n 表示全体「m×n、秩恰为 r」的复矩阵;AB=I 中的 I 按上下文分别是 Im(右逆)、In(左逆)。下面先证必要性(课件给出右逆情形),再给出两条独立的充分性路线。
必要性(课件证右逆,左逆对称,本编补)
右逆:设 AB=Im。由 定理 2.3.6(矩阵的秩 = 映射的秩)与秩的基本不等式 rank(AB)≤rankA(乘积的秩不超过任一因子的秩,本编注:R(AB)⊆R(A)),得
rank(AB)=m ⟹ m≤rankA≤min{m,n},
故 rankA=m,即 A 行满秩。
左逆(本编补:对 AH 用右逆结论):设 BA=In,两边取共轭转置得 AHBH=In,即 AH 有右逆。刚刚证过的右逆情形给出 rank(AH)=n,而共轭转置不改变秩(命题 3.1.1 的链条),故 rankA=n,即 A 列满秩。∎
一条推论(本编):AB=I 的等式还说明 A 不能比 B 更「宽」:A 行满秩要求 m≤n,A 列满秩要求 n≤m——单侧逆只在长方阵的「长边」方向上存在。
∎
充分性 · 法 1(相抵分解;课件原文,记号与引理 3.1.1 调顺)
设 A 行满秩,即 r=m。此时相抵标准形是横长条 Γ=[Im, 0](m×n),引理 3.1.1 给出
A=P−1[Im, 0]Q.
Γ 的右侧还有一根竖条可以「配平」——注意 [Im, 0][Im0]=Im。于是定义一个候选的右逆:
B=Q−1[Im0]P.
代入验算:
AB=P−1[Im, 0]Q⋅Q−1[Im0]P=P−1[Im, 0][Im0]P=P−1ImP=Im.
故 B 是 A 的右逆,充分性得证。∎
方法的要点:相抵标准形把「求右逆」变成了「求 Γ 的右逆」——而 Γ=[Im,0] 的右逆就摆在眼前:[Im0]。两端的可逆矩阵只需原样搬运过去。法 1 只用相抵分解(加、乘、求逆),所以对任意数域 F 都成立,不依赖内积。
∎
注 1 · 矩阵 A 右(左)逆存在性小结(课件原文)
- 当 A 为可逆矩阵时,其右(左)逆唯一存在,即为 A 的逆矩阵;
- 当 A 为行(列)满秩矩阵(非可逆阵)时,其右(左)逆存在且不唯一;
- 当 A 为其他情况(既不方阵满秩、也不相应满秩)时,其右(左)逆不存在。
本编补注 · 左逆 / 右逆在映射语言里是什么?
把 A 读成映射 T:Fn→Fm,x↦Ax(定理 2.3.4 的坐标公式),则:
- BA=In 说明 B 能把 T 的输出还原:若 T(x1)=T(x2),两边作用 B 得 x1=x2——这正是单射,而 定理 2.2.2 说单射 ⟺ 列满秩。左逆 = 单射的「反向还原器」;
- AB=Im 说明 B 能为每个输出定一个原像:任给 y,取 x=By 便有 T(x)=y——这正是满射。右逆 = 满射的「选原像器」(每个 y 选一个原像,选法即 B,所以不唯一)。
于是 定理 3.1.1 的三条与 2.2 的判据一一对应,只差一层翻译。这也解释了「不唯一」:单射的原像本来只有一个,所以右逆未必唯一(给值域外的每个点随便选原像);而列满秩时 A 作为映射单射,其左逆把像上的每个点送回唯一原像——不同之处在于左逆在像之外可以任意取值,所以同样不唯一(见 例 3.1.C)。
对照记忆:方阵情形两者合一(可逆),这是注 1 第一句;长方阵只能要一半,这是注 1 第二句。
例 3.1.C(新增)· 求全部右逆、全部左逆
(a) A=[1 0]∈C1×2,求它的所有右逆;(b) A=[10]∈C2×1,求它的所有左逆;(c) A=100010,写出一个左逆,并说明左逆有无穷多个。
解答 · 例 3.1.C
(a) A 行满秩(1×2,秩 1=m),右逆存在。设 B=[b1b2],则
AB=[1 0][b1b2]=b1 ⟹ AB=1⟺b1=1.
故全部右逆为 B=[1t](t∈C 任意)——无穷多个,确实不唯一(也印证注 1)。
(b) A 列满秩(2×1,秩 1=n),左逆存在。设 B=[b1 b2],则
BA=[b1 b2][10]=b1 ⟹ BA=1⟺b1=1.
全部左逆为 B=[1 t](t 任意)。
(c) 与 (b) 同理:设 B=[b11b21b12b22b13b23],则 BA=[b11b21b12b22](第 3 列乘零消失)。故 BA=I2⟺ 前两列为 I2 的两列,第三列 (b13,b23)T 任意。例如
B0=[100100],以及 [10017−5], ⋯
——左逆把「多余的一列」留给自由参数。几何上:A 把 C2 嵌入 C3(列向量张成 xy 平面),左逆只需把平面上的点送回去,平面外的分量怎么定都不影响 BA=I。
(a)(b)(c) 的共同点:三个例子都验证了注 1 的第二句——满秩非方阵的单侧逆有无穷多个,多出来的自由度来自「方阵之外的维数」:A 是 m×n 满秩时,差 ∣n−m∣ 个方向本就属于选法的自由。∎
法 1 用相抵分解造右逆,好处是不挑数域,代价是要先算相抵分解。还有一条更短的路线:直接写出一个公式当作右逆——比如 AH(AAH)−1——只需验证 AAH 可逆。而「可逆」是一个秩问题,于是问题归结为下面这条秩等式。课件把它列为命题,并用方程组给出证明,它是本节最重要的结论。
命题 3.1.1 · 秩的等式
rank(A)=rank(AAH),rank(AH)=rank(AHA),
从而(A 与 AH 的秩相等,故四个量相等)
rank(A)=rank(AHA)=rank(AAH)=rank(AH).
这是一个与矩阵秩相关的重要结论。
证明(课件两半,本编整理为「同解 ⟹ 同秩」三步)
第一步(证明 N(A)=N(AHA)):考查两个齐次方程组
Ax=0,AHAx=0.
- N(A)⊆N(AHA):若 z 满足 Az=0,两边左乘 AH 得 AHAz=0 ✓(这一步只用乘法,不需要任何内积)。
- N(AHA)⊆N(A):若 z 满足 AHAz=0,两边左乘 zH 得 zHAHAz=0,即
(Az,Az)=0.
由内积的正定性(定义 1.4.1(4):(u,u)=0⟺u=0)得 Az=0 ✓。
故两个方程组同解:N(A)=N(AHA)。
第二步(同解 ⟹ 同秩):两个矩阵同型(都是 n 列),用亏加秩定理(定理 2.2.5)
rank(A)=n−dimN(A)=n−dimN(AHA)=rank(AHA).
第三步(得到四个量相等):对 AH 重做第一步与第二步(AH 是 n×m 矩阵,其零空间在 Cm 中):
rank(AH)=rank((AH)HAH)=rank(AAH).
又转置(共轭转置)不改变秩:rank(A)=rank(AH)。串起来即得命题的四项等式链。∎
证明的要点:关键一步是正定性——(Az,Az)=0 只有正定内积才允许推出 Az=0(若换成一般二次型,「自己要内积为零」不等于「自己是零」)。这正是课件把话说死的那句「这是一个与矩阵秩相关的重要结论」的分量所在:它把「秩」这一个数,与内积全部挂钩。
∎
本编补注 · 为什么是 AHA 而不是 A2?
AHA 的两个好性质:① 一定是 Hermite 半正定矩阵(xHAHAx=(Ax,Ax)≥0,定义 1.4.3);② 它与 A 同秩(命题 3.1.1)。于是
AHA 可逆 ⟺ rank(AHA)=n ⟺ rankA=n ⟺ A 列满秩;
AAH 可逆 ⟺ rank(AAH)=m ⟺ A 行满秩。
这两条是本课程后面反复引用的「可逆性转换器」:长方阵的满秩问题,转成方阵 AHA(或 AAH)的可逆问题——后者可以用行列式、可以用正定性、可以直接求逆。(换成 A2 就全塌了:A 非方阵时 A2 根本没法定义;就算方阵也不例外,取 A=[0010] 便有 A2=O,秩从 1 掉到 0。)
充分性 · 法 2(用命题 3.1.1;课件原文)
设 A∈Cm×n 为行满秩矩阵(rankA=m)。由 命题 3.1.1,
rank(AAH)=rank(A)=m,
即 m 阶方阵 AAH 满秩,可逆。于是
AAH(AAH)−1=Im,
即
B=AH(AAH)−1
是矩阵 A 的一个右逆(A⋅AH(AAH)−1=Im ✓)。∎
与法 1 的分工(本编):法 1(相抵)给出「存在、且可机械地算出」,用不着内积,对数域 F 通用;法 2(本节)给出闭式公式,只需求逆一次,但依赖 AH(用于实/复内积空间)。考试时能写出来就是好办法:法 2 的公式短、验证快,是首选;法 1 是「万一忘了公式」时的保底路线。对称地,列满秩时 B=(AHA)−1AH 是 A 的一个左逆。
∎
课程预告 · 这两个公式后面会一直出现
B=AH(AAH)−1(右逆)与 B=(AHA)−1AH(左逆)不是两个临时的技巧——它们是广义逆(伪逆)的雏形,也是「最小二乘」的核心:
- 对列满秩的超定方程组 Ax=b(方程个数多于未知数,通常无解),最小的「残差平方和」解是
x=(AHA)−1AHb,
其中 (AHA)−1AH 正是 A 的一个左逆(对称地,行满秩时 AH(AAH)−1 是右逆);等式 (AHA)x=AHb 称为正规方程;
- 对任何矩阵(包括不满秩、不对称的情形),上述两条路线可以被统一成同一个对象 A+(Moore–Penrose 伪逆),它满足 AA+A=A——这正是 3.2 的例 3.2.2 将要构造的 X。
本编此处只用到 定理 3.1.1 与 命题 3.1.1 的结论;广义逆的完整理论与最小二乘的推导留待后续课程内容展开,这里先埋一个锚点。
命题 3.1.1 两页的表述次序不一致(本编勘误 · E3.2)
课件第 12 页把两条等式写成 rank(A)=rank(AAH) 与 rank(AH)=rank(AHA),第 13 页又把它们写成 rank(A)=rank(AHA) 与 rank(AH)=rank(AAH)——两页的配对次序相反。这不是矛盾:四个量本来就相等(见上方证明第三步),两种写法都成立;但抄写时容易抄错配对,建议直接背四项等式链
rank(A)=rank(AHA)=rank(AAH)=rank(AH),
它不依赖先写哪一对,也把两页的写法一并涵盖。课件所证的主线是 N(A)=N(AHA) 与由此得到的 rank(A)=rank(AHA)。
例 3.1.D(新增)· 用公式造右逆并验证
设 A=[102101]。(a) 判定 A 是否有右逆;(b) 若有,用 定理 3.1.1 法 2 的公式造一个并验算。
解答 · 例 3.1.D
(a) A 是 2×3 矩阵,秩为 2=m(两行不成比例),行满秩,故右逆存在(且由注 1 不唯一)。
(b) 先算
AAH=[102101]120011=[5222],det(AAH)=10−4=6=0,
可逆(也验证了 命题 3.1.1:rank(AAH)=rank(A)=2 ✓)。求逆:
(AAH)−1=61[2−2−25].
于是
B=AH(AAH)−1=61120011[2−2−25]=6122−2−215.
验算:
AB=61[102101]22−2−215=61[6006]=I2 ✓
顺带核对不唯一性:由注 1,右逆不唯一——本例中 A 的右逆构成一个 2 参数的仿射族(每个输出的原像可加 N(A) 中任意向量的坐标,n−r=1,两列各自贡献一个参数),完整写法见习题 C4。∎
- 一条主线:任何矩阵都能相抵到标准形(引理 3.1.1:A=P−1ΓQ,证明 = 选基 + 换基公式)⇒ 分解式不唯一、标准形唯一(本编补注:相抵 ⟺ 同秩)⇒ 用它解方程组(例 3.1.1:可解判据 = 看 Pb 的后 m−r 个分量,自由度数 =n−r)⇒ 用它造单侧逆、并给出满秩判别(定理 3.1.1:行满秩 ⟺ 右逆,列满秩 ⟺ 左逆)⇒ 最后通过秩等式 rank(A)=rank(AHA)(命题 3.1.1)拿到闭式的单侧逆公式 AH(AAH)−1、(AHA)−1AH。
- 一条对应:分解 ↔ 秩。相抵标准形由秩唯一决定;满秩分解(3.2)的中间因子个数是秩;单侧逆的存在性由「是否达到 m 或 n」决定。本章所有结论都是「用结构把秩显出来」。
常见错误清单 · 自查五条
- 把相抵分解当成唯一的:不唯一的是分解式(P、Q 有无穷多种),唯一的是标准形 Γ;别把两句话写反;
- 分块矩阵记错栏目:右上块积的是行变换(P−1=E),左下块积的是列变换(Q−1=F);算完一定用 EAF=Γ 回代;
- 忘记「方阵」前提:AB=I 只在方阵情形推出 BA=I;长方阵只能得单侧逆,且右逆要求行满秩、左逆要求列满秩,方向别写反;
- x 与 y 混用:用相抵分解解方程组时,Γy=Pb 解出的是 y,必须乘回 Q−1 才是 x;可解判据看的是 Pb(不是 b)的后 m−r 个分量;
- 把 rank(A)=rank(AHA) 当默认事实:它的证明用了内积正定性,是「实 / 复内积空间」的结论;写成 ATA 或一般域上的 A∗A 时要先确认前提(本课程一律在实 / 复内积空间中使用)。
即时自测(不看上文,5 分钟)
- 写出相抵分解式与相抵标准形;「分解不唯一」与「标准形唯一」各指什么?
- A=[1224]:写出 A 的相抵标准形与一个分解式(P−1、Q 都给出)。
- A∈C3×2 且 rankA=2:A 有右逆还是左逆?AAH 与 AHA 哪个可逆?另一个的秩是多少?
- 默写 命题 3.1.1 的证明主线,并指出哪一步用了内积的哪条公理。
解答 · 即时自测
- A=P−1ΓQ,Γ=[Ir000];不唯一指 P,Q 的取法无穷(对任意 ∗,P−1[Ir0∗I] 与 [Ir∗0I]Q 仍给出分解);唯一指 r 决定 Γ,无论哪种分解式中间都是同一个 Γ。
- 秩 1,Γ=[1000]。一个分解:第一行减第二行的一半方向——取 P=[1−201](r2←r2−2r1),F=Q−1=[10−21](c2←c2−2c1),则 PAQ−1=[1000]=Γ,即 A=P−1ΓQ。验算:P−1ΓQ=[1201][1000][1021]=[1224] ✓。
- 列满秩(rank=2=n):有左逆,没有右逆(m=3>n=2 不满足行满秩);AHA(2×2)可逆;AAH(3×3)不可逆,由 命题 3.1.1 其秩为 2。
- 主线:N(A)⊆N(AHA)(两边左乘 AH)与 N(AHA)⊆N(A)(左乘 zH 得 (Az,Az)=0,由正定性推出 Az=0)⇒ 同解 ⇒ 亏加秩给出 rank(A)=rank(AHA);对 AH 重复即得四项等式链。用到的公理是 定义 1.4.1(4)正定性。