把行变换记录成一个矩阵
初等行变换最初是矩阵上的操作步骤:交换两行、把一行乘以非零标量、把一行的倍数加到另一行。在增广矩阵上,只有把操作施加到受影响的整行——包括增广线右边的常数项——才是在等价地改写方程。此时操作可以逆转,因而保留原方程组的解集。
但同一件事还有一个更有结构的读法:每一个初等行变换,都等同于左乘某个特别的方阵。这个观点很重要,因为它把一串初等行变换转成普通的矩阵等式。之后讨论行化简、可逆性、秩、行列式和基底时,这正是程序性操作和代数结构之间的桥。
这个结构观点也回答一个实际问题:若同一个行变换要同时作用于许多列,怎样只记录一次规则,而不逐项重做?单位矩阵恰好提供这种编码。它的每一行会记录新行由哪些旧行组成;随后,矩阵乘法便把这套规则同时施加到相容矩阵的每一列。
从单位矩阵建立初等矩阵
假设一个初等行变换 ρ 作用于有 p 行的矩阵。先从单位矩阵 Ip
开始,对 Ip 做同一个初等行变换,所得矩阵记作 Eρ。所有标量都属于矩阵所采用的底域;特别地,每个非零倍数 c 在该域中都有倒数
1/c。
与 p 行矩阵上的初等行变换 ρ 对应的初等行变换矩阵,就是把
ρ 作用在 Ip 上所得的 p×p 矩阵。
这种矩阵也简称初等矩阵;较长的名称可以明确指出它与指定行变换的关系。
换句话说,
Eρ=ρ(Ip).
这个定义有用,是因为以下定理。
设 A 是任何有 p 行的矩阵,ρ 是一个初等行变换,并以
ρ(A) 表示把它作用于 A 后得到的矩阵。则
ρ(A)=EρA.所以,对 A 做这个初等行变换,等同于用同一个操作在 Ip 上得到的初等行变换矩阵左乘 A。
这里必须是左乘。初等行变换改变的是行,会把新行写成旧行的线性组合;左乘正是把
A 的新行写成旧行的线性组合。右乘则会混合列,不是同一件事。
维数也说明了同一点。若 A 是 p×n 矩阵,则 Eρ 是
p×p 矩阵,所以 EρA 有定义,而且仍是 p×n。这里完全不要求 A 是方阵。相反,对这 n 列进行列变换的矩阵大小为 n×n,并从右边相乘。
为什么左乘会执行指定行变换
把任意 p×n 矩阵 A 的各行记作
a1,a2,…,ap,并以 ekT 表示 Ip 的第 k 行。基本等式
ekTA=ak表示左乘矩阵的一行会选取 A 的相应行。现在逐类检查。
-
对 Ri↔Rj,Eρ 的第 i 行是
ejT,第 j 行是 eiT,其余第 k 行仍是
ekT。因此 EρA 恰好交换 A 的第 i,j 行。
-
对 Ri←cRi(c=0),Eρ 的第 i 行是
ceiT,其余行不变。因此乘积的第 i 行为
(ceiT)A=cai,
而每个 k=i 的第 k 行仍为 ak。
-
对 Rj←Rj+cRi(i=j),Eρ 的第 j 行是
ejT+ceiT,其余行不变,所以新的第 j 行为
(ejT+ceiT)A=aj+cai,
正是指定的行加法。
三类情形穷尽了初等行变换,因此对每个相容矩阵 A 都有
EρA=ρ(A)。这个证明没有声称任意改动行的方式都对应一个初等矩阵;
ρ 必须是刚才检查的三类操作之一。
作用于增广矩阵
假设 Ax=b 由增广矩阵 [A∣b] 表示。初等行变换必须作用于整个增广矩阵,不能只改其中一边。相应的乘法等式是
Eρ[A∣b]=[EρA∣Eρb].
因此,系数和右端常数会接受完全相同的行组合。向量 x 满足 Ax=b,当且仅当它满足 EρAx=Eρb:正向可把原等式左乘 Eρ,反向则施加逆向初等行变换。这才是解集保持不变的严格理由。
若只改 A 而固定 b,或只改 b,一般会得到另一个方程组。增广线只是视觉上的分隔,不是行变换可以中途停止的边界。
分别记录替换、缩放与交换
三种初等行变换,对应三类初等行变换矩阵。
对有三行的矩阵,考虑
ρ:R2←R2+3R1.把它作用在 I3 上:
I3=100010001⟼Eρ=130010001.因此,对任何有三行的矩阵 A,
EρA就是把 A 的第 2 行改成“第 2 行加上 3 倍第 1 行”后所得的矩阵。
对
ρ:R3←−2R3,对应的初等行变换矩阵是
Eρ=10001000−2.行倍乘要求标量非零,这一点在矩阵上也看得出来:如果倍数是 0,所得矩阵会出现零行,不能反向恢复。
对
ρ:R1↔R3,我们得到
Eρ=001010100.左乘这个矩阵,就会交换任何相容矩阵的第 1 行和第 3 行。
逐步理解初等矩阵的左乘
对单位矩阵施加行变换,便记录了该变换的系数;再用所得矩阵左乘,就把同一个行变换施加到每一列。
把行变换读成左乘看清初等行变换如何变成左乘矩阵:先把操作作用在单位矩阵上,再用所得矩阵改变任何相容矩阵的行。
由 I 开始
对 ρ: R₂ ← R₂ + 3R₁,先把操作作用在 I₃ 上。
建立 Eρ
只有 I₃ 的第 2 行改变,得到第 2 行为 [3,1,0] 的 Eρ。
放在左边相乘
同一个矩阵会对 A 执行该操作:ρ(A)=Eρ A。
读新行
新的第 2 行等于旧第 2 行加三份旧第 1 行。
保留乘积顺序
在 A₃=E₂E₁A₁ 中,E₁ 最靠近 A₁,因为 ρ₁ 最先作用。
反向得到逆矩阵
R₂ ← R₂ + 3R₁ 的反向操作是 R₂ ← R₂ − 3R₁,所以相应矩阵相乘得到 I。
行变换矩阵是把行变换先作用在单位矩阵上得到的。之后用它左乘任何相容矩阵,就会执行同一个行变换;反向行变换则给出逆矩阵。
比较行变换与矩阵乘积
用下面的步骤工具重温一串行变换。每一步都要同时读成:在左侧乘上对应的初等矩阵。
互动 · 逐步演示第 1 步,共 5 步
跟着走完一条行化简路径
每一步做一组初等行变换。留意主元上下的元素怎样逐个变成 0。
先在第 1 列选主元。
1正在处理的主元0本步改动竖线右边:常数
先由增广矩阵开始。第一个主元的工作,是帮我们把它下面的元素清掉。
还未到行阶梯形
第 2 行第 1 列仍是 1,没有在第 1 行主元的右方。
要留意第 1 列的第一行已经有方便的主元 1,所以暂时不用换行。
一串初等行变换就是一个矩阵乘积
真正的用途不只是表示一个初等步骤,而是把整串初等行变换写成一个矩阵乘积。
假设
A1ρ1A2ρ2A3ρ3⋯ρkAk+1.则
Ak+1=EρkEρk−1⋯Eρ2Eρ1A1.这里每个 ρj 都必须是作用于同样行数矩阵的初等行变换。这个乘积表示整串操作,但乘积本身未必仍是单一个初等行变换矩阵。
这个公式中的顺序很重要。第一个行变换最靠近 A1,因为它最先作用:
A2=Eρ1A1,A3=Eρ2A2=Eρ2Eρ1A1.
令
A=101020110.依次做
ρ1: R2←R2+R1,ρ2: R1←R1+2R2.对应的行变换矩阵是
Eρ1=110010001,Eρ2=100210001.两个操作完成后,结果就是
Eρ2Eρ1A=310210001101020110=311420520.若不想直接相乘 Eρ2Eρ1,也可以把同一串行变换按同一顺序作用在
I3 上,所得矩阵就是合并后的左乘矩阵。最后的数值结果也确认:第二步使用的是已经更新后的第 2 行。
读懂较长的行变换乘积
把一长串初等行变换写成矩阵乘积时,需要先分清三种对象:
- 正在被变换的矩阵 A1,A2,…,Ak+1;
- 每一步对应的初等行变换矩阵 H1,H2,…,Hk;
- 合并后的单一左乘矩阵 J=Hk⋯H2H1。
以下是一个四行矩阵的典型例子。假设
ρ1ρ2ρ3ρ4ρ5ρ6:R2←−21R2,:R1↔R2,:R3←R3−2R1,:R4←R4−2R3,:R1←R1−R3,:R1←R1−R2.
如果
A1ρ1A2ρ2A3ρ3A4ρ4A5ρ5A6ρ6A7,
那么
A7=H6H5H4H3H2H1A1.
相应的行变换矩阵是
H1H3H5=10000−210000100001,=10−20010000100001,=10000100−10100001,H2H4H6=0100100000100001,=10000100001−20001,=1000−110000100001.
它们的乘积是
J=H6H5H4H3H2H1=−1100−2301−2−101−20001.
有效率的做法,是把六个行变换依次作用在 I4 上,而不是手动展开六个因子。矩阵 J 记录了整串步骤对行的总作用:
A7=JA1.
这个等式也可用来检查乘积顺序。如果把第一步放在最左边,所得乘积会描述另一串不同的操作。
快速检查
(在上面六步行变换中,若 K 是由 A7 反向回到 A1 的行变换矩阵乘积,J 和 K 应满足什么等式?)
解答 · 答案
矩阵 J 和 K 编码的是每一个四行矩阵上的变换,不只是题中展示的
A1。对任意四行矩阵 X,先做正向步骤,再做反向步骤,会得到
K(JX)=X.取 X=I4,便有 KJ=I4。反过来,先做反向步骤,再做正向步骤,也会恢复任意四行矩阵 Y,所以 J(KY)=Y;取 Y=I4,便有 JK=I4。因此
KJ=I4JK=I4.换句话说,K=J−1。
反向行变换与逆矩阵
每个初等行变换都有反向操作:
- Rj←Rj+cRi 的反向操作是
Rj←Rj−cRi;
- Ri←cRi(其中 c=0)的反向操作是
Ri←(1/c)Ri;
- 交换两行的反向操作仍然是同一次交换。
这给出一个精确的矩阵表述。
每个初等行变换矩阵都是可逆的。它的逆矩阵,就是反向初等行变换所对应的初等行变换矩阵。
设 Eρ 表示初等行变换 ρ,而 Eρ−1 表示它的反向初等操作。对每个有 p 行的矩阵 X,先做 ρ 再撤销,都会恢复 X。由左乘定理,
Eρ−1EρX=X.这个等式对每个这样的 X 都成立;特别取 X=Ip,得到
Eρ−1Eρ=Ip。先做反向操作再做 ρ,同理得到
EρEρ−1=Ip。因此
Eρ−1=Eρ−1.具体来说,交换两行的逆仍是同一次交换;乘以 c=0 的逆是乘以 1/c;把第 i 行的 c 倍加到第 j 行,则由加上 −c 倍第 i 行来撤销。
例如,若
E=130010001
表示 R2←R2+3R1,则
E−1=1−30010001
表示 R2←R2−3R1。
这就是初等行变换可逆的代数原因,也解释了为什么行化简和可逆矩阵有密切关系。
设 E1,…,Ek 按时间顺序编码同一组 p 行矩阵上的操作,合并左乘矩阵就是 J=Ek⋯E1。每个 Ej−1 表示对应初等行变换的反向操作。要恢复输入,首先必须撤销 Ek,因为它是最后一次对行所作的改变;然后依次撤销 Ek−1,…,E1。因此,反向链的合并矩阵为
K=E1−1E2−1⋯Ek−1.这里必须区分两种次序:反向操作的时间顺序从 Ek−1 开始;写成乘积时,最先作用的因子却放在最靠近待恢复矩阵的位置。以三步为例,由结合律可得
KJ=E1−1E2−1(E3−1E3)E2E1=E1−1(E2−1E2)E1=Ip.另一个乘积 JK 则从相邻的 E1E1−1 开始消去,最后同样得到 Ip。任意有限链都可这样逐对消去。整个论证没有交换因子,也不要求被变换的矩阵本身是方阵。
把这个读法用于前面的两步例子:合并矩阵是 Eρ2Eρ1,所以逆矩阵是 Eρ1−1Eρ2−1。从最终矩阵开始,先从第一行减去当前第二行的两倍,第一行便恢复为 [1 0 1];再从第二行减去刚恢复的第一行,第二行便恢复为 [0 2 1]。撤销操作是在逐个恢复中间状态,乘积次序也就随之确定,无须另背一个口诀。
为什么这个观点之后会有用?
如果 B 与 A 行等价,则存在有限多个初等行变换把 A 变成 B。因此存在某个初等行变换矩阵的乘积 E,使得
B=EA.
由于每个初等行变换矩阵都可逆,乘积 E 也可逆。因此行等价既可以用初等操作步骤来描述,也可以用左边乘上一个可逆矩阵的等式来描述。反向叙述则要小心:任意可逆左乘矩阵未必本身只代表一步初等操作。
这个观点会在多个地方使用:
- 一个方阵若行等价于 In,它就是若干初等行变换矩阵的乘积;
- 初等行变换保留齐次方程组的解信息,因为它等同于左乘可逆矩阵;
- 行列式关于初等行变换的规则,可以通过初等矩阵来表述;
- 秩和基底的论证可以使用行化简,同时清楚知道原矩阵的列本身通常已被改变。
常见错误
初等行变换由左乘表示。右乘会混合列,而不是混合行。
如果 ρ1 先于 ρ2 作用,合并矩阵是
Eρ2Eρ1,不是 Eρ1Eρ2。
行倍乘要求标量非零。把一行乘以 0 不能反向恢复,也不会产生可逆的初等行变换矩阵。
从某个特定矩阵 A 满足 KJA=A,一般不能推出 KJ=I,因为 A 的各列未必能检测所有可能输入。应证明复合行变换会恢复每个相容矩阵,或直接把它们作用于 Ip。所以上面对 KJ=JK=Ip 的论证先使用任意 p 行矩阵,再代入单位矩阵。
按作用次序理解乘积
要编码有 p 行矩阵上的一个初等行变换,只需把该操作施加到 Ip。所得初等矩阵 Eρ 会以左乘方式作用于每个相容矩阵;对交换、非零倍乘和行加法三类情形逐行检查,就得到 EρA=ρ(A)。若对象是增广矩阵,操作必须横跨整行,包括右端常数,方程组才会保持完全相同的解集。
依次施加 ρ1,…,ρk 时,乘积按作用顺序写成
Eρk⋯Eρ1。每个因子都可逆,因为反向初等操作给出它的双侧逆矩阵。因此,有限的行化简既可以看作逐步操作,也可以看作一次可逆的左乘;各个因子则继续记录每一步的准确内容。
快速检查
快速检查
(对有三行的矩阵,E=105010001 表示哪个行变换?)
解答 · 答案
只有第 3 行被改变,而且它变成“第 3 行加上 5 倍第 1 行”。因此
R3←R3+5R1.
快速检查
(R2←R2−4R1 的反向行变换是什么?)
解答 · 答案
反向操作是
R2←R2+4R1.
练习
快速检查
(写出有三行矩阵上 R1↔R2 的行变换矩阵。)
解答 · 引导解答
交换 I3 的前两行,得到
010100001.这就是左乘时交换第 1 行和第 2 行的矩阵。
快速检查
(假设 B 由 A 经以下两步得到:先做 R2←R2+R1,再做 R3←2R3。把 B 写成涉及 A 的矩阵乘积。)
解答 · 引导解答
令 E1 是 R2←R2+R1 的行变换矩阵,令 E2 是
R3←2R3 的行变换矩阵。因为 E1 先作用,而 E2 后作用,
B=E2E1A.
快速检查
(令 β1,β2=0。假设五行矩阵 A 经以下行变换链得到 B:先做 R3←R3+α1R1,再做 R2←β1R2,再做 R1↔R4,再做 R3←R3+α2R2,最后做 R1←β2R1。写出单一矩阵 G,使得 B=GA。)
把同一串操作作用在 I5 上,并记住后面的操作使用的是当时的行,不是原来的行。
解答 · 引导解答
把这些操作作用在 I5 上,得到
G=00α1100β1α2β10000100β2000000001.重点是 α2β1 这一项。操作
R3←R3+α2R2 发生在第 2 行已经被缩放成
β1R2 之后,所以来自原来第 2 行的贡献是
α2β1R2。
快速检查
(沿用上一题的矩阵 G,假设 D=GC。写出矩阵 H,使得 C=HD。)
解答 · 引导解答
反向链是
GR1←β21R1R3←R3−α2R2R1↔R4R2←β11R2R3←R3−α1R1I5.因此
H=000β2100β11−α2000010010−α10000001.这就是 G−1,所以 D=GC 会推出 C=HD。
相关笔记
这一页建立在
2.2 增广矩阵与行变换
和
3.2 矩阵乘法、单位矩阵与线性方程组
之上。接着可以读
3.6 分块矩阵,完成第三章。它也准备了
5.1 可逆矩阵
和
7.2 行变换、乘积与可逆性
中使用的行化简代数观点。
补充练习与完整解答
题 1
对(a)–(d),以下按次序的行变换把A送到B。写出从B返回A的操作链,并给出满足B=HA的方阵H。大小分别为(a)、(b):3×4;(c):3×6;(d):4×5。
(a)
1. R2←R2+(−1)R12. R3←R3+(−2)R13. R3←R3+(−2)R24. R3←(−1)R35. R1←R1+(−2)R26. R2←R2+(−1)R3
(b)
1. R1↔R22. R1←(−1)R13. R3←R3+(−2)R14. R3←R3+(−3)R25. R1←R1+(−2)R2
(c)
1. R1↔R22. R3←R3+(2)R13. R3←R3+(−3)R24. R1←R1+(−2)R25. R1←R1+(2)R36. R2←R2+(−2)R3
(d)
1. R2←R2+(−1)R12. R3←R3+(−3)R13. R4←R4+(−1)R14. R2←(−1)R25. R3←R3+(4)R26. R4←R4+(3)R27. R3↔R48. R4←R4+(4)R39. R1←R1+(−2)R210. R1←R1+(−1)R3
完整解答
以下在实数上讨论。行替换的源行与目标行不同,缩放系数非零;操作作用于整行,包括所有增广块。有限链中第一步先作用,所以合并左乘矩阵中的初等矩阵按相反次序排列。
把正向操作逐一作用于Ir。以下Hj记录每个中间乘数;对步数归纳,当前矩阵就是HjA,最后一个乘数即所求H。倒序撤销各步:加行倍数改为减同一倍数,缩放改为乘倒数,交换仍作同一交换。每个逆操作对任意相容矩阵都撤销原操作,因此所写逆链确实返回A。
(a)
H0=100010001
H1=1−10010001
H2=1−1−2010001
H3=1−1001−2001
H4=1−1001200−1
H5=3−10−21200−1
H6=3−10−2−1201−1
逆链按以下次序:
1. R2←R2+(1)R32. R1←R1+(2)R23. R3←(−1)R34. R3←R3+(2)R25. R3←R3+(2)R16. R2←R2+(1)R1
(b)
H0=100010001
H1=010100001
H2=010−100001
H3=010−102001
H4=01−3−102001
H5=−21−3−102001
逆链按以下次序:
1. R1←R1+(2)R22. R3←R3+(3)R23. R3←R3+(2)R14. R1←(−1)R15. R1↔R2
(c)
H0=100010001
H1=010100001
H2=010102001
H3=01−3102001
H4=−21−3102001
H5=−81−3502201
H6=−87−35−422−21
逆链按以下次序:
1. R2←R2+(2)R32. R1←R1+(−2)R33. R1←R1+(2)R24. R3←R3+(3)R25. R3←R3+(−2)R16. R1↔R2
(d)
H0=1000010000100001
H1=1−100010000100001
H2=1−1−30010000100001
H3=1−1−3−1010000100001
H4=11−3−10−10000100001
H5=111−10−1−4000100001
H6=11120−1−4−300100001
H7=11210−1−3−400010010
H8=11290−1−3−1600010014
H9=−11292−1−3−1600010014
H10=−31295−1−3−160001−1014
逆链按以下次序:
1. R1←R1+(1)R32. R1←R1+(2)R23. R4←R4+(−4)R34. R3↔R45. R4←R4+(−3)R26. R3←R3+(−4)R27. R2←(−1)R28. R4←R4+(1)R19. R3←R3+(3)R110. R2←R2+(1)R1
题 3
设A,A′为m×p矩阵,B为p×q矩阵,且A′=ρ(A)。(a) 原题要求直接由定义证明AB=ρ(A′B),其中(i)Rk←Rk+αRi、i=k;(ii)Rk←βRk、β=0;(iii)Ri↔Rk。(b) 用归纳法证明:若A,A′行等价,则AB,A′B在同一有限操作序列下行等价。
完整解答
以下在实数上讨论。行替换的源行与目标行不同,缩放系数非零;操作作用于整行,包括所有增广块。有限链中第一步先作用,所以合并左乘矩阵中的初等矩阵按相反次序排列。
(a):(i)、(ii)原印方向不成立
对(i),取A=(1,0)T、B=[1],并作R2←R2+R1。则A′B=(1,1)T,再作同一操作得到(1,2)T=AB。对(ii),取A=B=[1]、β=2,则A′B=[2],而ρ(A′B)=[4]=[1]。所以原印的普遍命题不能被证明。正确方向是A′B=ρ(AB),等价地AB=ρ−1(A′B)。交换操作满足ρ−1=ρ,所以(iii)原印方向成立。
直接逐项证明正向恒等式:在行替换情形,只有输出第k行改变:
(A′B)kj=r=1∑p(akr+αair)brj=(AB)kj+α(AB)ij.
其余输出行逐项相同。缩放时(A′B)kj=β(AB)kj,其余行不变。交换A的第i,k行,会恰好交换AB的第i,k行。由此直接证明三种正向传递恒等式;再作逆操作,便得到反向恒等式。
(b):对行变换序列长度归纳
长度为零时,两矩阵相同,乘积也相同。假设结论对长度n成立。对长度n+1的序列,令An为前n步后的矩阵。归纳假设表明,同样的前n步把AB送到AnB;再由(a)已逐项证明的正向恒等式,最后一步把AnB送到A′B。所以整个n+1步序列把AB送到A′B。归纳得任意有限长度的结论。所有乘积大小均为m×q,因此始终是同尺寸矩阵上的合法行变换。
题 4
3×5矩阵依次满足A1R3+2R1A2R2↔R3A33R1A4。(a)–(c) 给出满足Aj+1=JjAj的J1,J2,J3;(d) 给出A4=JA1中的J;(e) 给出A1=J′A4中的J′。
完整解答
以下在实数上讨论。行替换的源行与目标行不同,缩放系数非零;操作作用于整行,包括所有增广块。有限链中第一步先作用,所以合并左乘矩阵中的初等矩阵按相反次序排列。
每个Jj由对应的一步操作作用于I3得到。连续代入可得A4=J3J2J1A1。具体为:
J1=102010001
J2=100001010
J3=300010001
J=J3J2J1=320001010
J′=J−1=310−32001010
逆过程先把第1行乘1/3,再交换第2、3行,最后从第3行减去第1行的两倍。把这些操作作用于所写J得到I3;因此J′J=JJ′=I3,对任意A1都得到要求的反向等式。
题 5
下面六个4×6矩阵是连续阶段A,A(1),…,A(5)=A′。(a) 指出五步行变换ρ1,…,ρ5;(b) 给出满足A′=HA的明确4×4矩阵H。
A0stage=210−1131−2−101348403−4−308−4−4−1
A1stage=120−1311−20−1138440−43−30−48−4−1
A2stage=100−13−51−20−1138−1240−411−30−416−4−1
A3stage=10003−5110−1138−1248−411−3−4−416−4−5
A4stage=100031−5101−1384−128−4−311−4−4−416−5
A5stage=1000311101513845128−4−3−511−4−4−4−516−5
完整解答
以下在实数上讨论。行替换的源行与目标行不同,缩放系数非零;操作作用于整行,包括所有增广块。有限链中第一步先作用,所以合并左乘矩阵中的初等矩阵按相反次序排列。
(a) 比较整行,包括最后两列,依次得到:
1. R1↔R22. R2←R2+(−2)R13. R4←R4+(1)R14. R2↔R35. R3←(−51)R3
例如,交换后第2行减第1行的两倍,得到(0,−5,−1,−12,11,16);第4行加第1行,得到(0,1,3,8,−4,−5)。第四步交换第2、3行;第五步把当时的第3行乘−1/5,最后两项也同时缩放。(b) 把同样操作作用于I4,逐步得到:
H0=1000010000100001
H1=0100100000100001
H2=01001−20000100001
H3=01001−20100100001
H4=001010−2101000001
H5=00−5101052101000001
取H=H5。每个记录行恰为相应原行的线性组合,所以逐项证明了A′=HA。
题 6
考虑下面一串5×7矩阵,并将中间阶段记为U0,U1,U2,U3。箭头为Aρ1U0→⋯→U1ρ2U2→⋯→U3ρ3B1,其后是所示五步B1→⋯→B6。(a) 求ρ1,ρ2,ρ3;(b) 填写B3(I)(II)(III)B1的三处空白;(c) 求满足B5=KB4的K;(d) 是否存在从I5到某个满足B6=HB1的H的行变换链?若有,须写出全部操作和每个中间矩阵,包括H;若无,说明不存在。
A=311−14212−13113−12−20−1−2−2132−94−2−1−22−3−311−2−2
U0=131−14122−13113−120−2−1−2−2312−94−1−2−22−31−31−2−2
U1=100041−11031−22020−2−1−2−23−8−1−64−11−11−31−60−1−2
U2=100001−110−11−220−20−2−1−2−23−8−1−6−8−11−1111−60−1−6
U3=10000110001200002−3−2038−9−60−1−101016−6−10
B1=100001100012000021−20383−60−1−1010162−10
B2=1000011000120000210038300−1−101016230
B3=1000001000−12000−22100−583000−1010−56230
B4=1000001000−1200002100183000−1010−16230
B5=1000001000−1200000100123000−1010−12230
B6=1000001000−12000001001230000010−15230
从B1开始的五个箭头依次为:
1. R4←R4+(2)R32. R1←R1+(−1)R23. R1←R1+(2)R34. R2←R2+(−2)R35. R2←R2+(1)R4
完整解答
以下在实数上讨论。行替换的源行与目标行不同,缩放系数非零;操作作用于整行,包括所有增广块。有限链中第一步先作用,所以合并左乘矩阵中的初等矩阵按相反次序排列。
(a) 比较三个标记箭头两侧的整行:
ρ1:R1↔R2,ρ2:R5←R5−4R1,ρ3:R3←−31R3.
对ρ2,从(4,3,2,−2,4,−3,−2)减去(1,1,1,0,3,−1,1)的四倍,得到(0,−1,−2,−2,−8,1,−6)。对ρ3,缩放(0,0,0,−3,−9,0,−6)后得到(0,0,0,1,3,0,2)。(b) 把B1→B2→B3的两步倒序撤销:
(I):R1←R1+R2,(II)=B2,(III):R4←R4−2R3.
(c) 只有第2行改为第2行减第3行的两倍,故:
K=10000010000−21000001000001
(d) 存在。对I5作以下同一组五步操作:
1. R4←R4+(2)R32. R1←R1+(−1)R23. R1←R1+(2)R34. R2←R2+(−2)R35. R2←R2+(1)R4
H0=1000001000001000001000001
H1=1000001000001200001000001
H2=10000−11000001200001000001
H3=10000−11000201200001000001
H4=10000−110002−21200001000001
H5=10000−11000201200101000001
取H=H5。上面已经写出I5=H0→H1→⋯→H5的每个矩阵。每一步都以同一方式改变原输入行的组合,因此把操作链作用于B1时,第j步有HjB1=Bj+1;特别地HB1=B6。这对矩形B1同样成立,完全不需要把B1求逆。
题 7
3×4矩阵A1,…,A5按下面四步操作相连。(a) 求满足Aj+1=JjAj的J1,…,J4;(b) 若增广矩阵[A1∣B1],…,[A4∣B4],[A5∣I3]也按完全相同的四步相连,各Bj均为3×3,求全部四个Bj;(c) B1A1是否与A5行等价?若是,写出操作链的每一步;若否,解释原因。
1. R2←R2+(−1)R12. R3←R3+(−2)R13. R3←R3+(−1)R24. R1←R1+(−1)R2
完整解答
以下在实数上讨论。行替换的源行与目标行不同,缩放系数非零;操作作用于整行,包括所有增广块。有限链中第一步先作用,所以合并左乘矩阵中的初等矩阵按相反次序排列。
(a) 把各步指定操作作用于I3,得到:
J1=1−10010001
J2=10−2010001
J3=10001−1001
J4=100−110001
(b) 每个行变换也作用于右块,因此J1B1=B2、J2B2=B3、J3B3=B4、J4B4=I3。从最后一个等式开始,依次施行逆操作:
B4=J4−1,B3=J3−1J4−1,B2=J2−1J3−1J4−1,B1=J1−1J2−1J3−1J4−1.
B1=112123001
B2=102113001
B3=100111001
B4=100110001
把这些明确矩阵代入上述四个右块等式,也可逐一核对;并没有对矩形Aj作非法消去。(c) 是。令P=J4J3J2J1。增广操作链同时给出A5=PA1及PB1=I3。所以对B1A1施行一次四步操作,得到P(B1A1)=A1;再作一次四步,得到PA1=A5。因此完整操作链依次为:
1. R2←R2+(−1)R12. R3←R3+(−2)R13. R3←R3+(−1)R24. R1←R1+(−1)R25. R2←R2+(−1)R16. R3←R3+(−2)R17. R3←R3+(−1)R28. R1←R1+(−1)R2
题 8
设α1,α2,α3,β1,β2,β3为任意实数。四行矩阵A,B经下面六步操作相连。给出满足B=GA及A=HB的明确4×4矩阵G,H。
1. R2←R2+(α1)R12. R3←R3+(α2)R13. R4←R4+(α3)R14. R1←R1+(β1)R25. R1←R1+(β2)R36. R1←R1+(β3)R4
完整解答
以下在实数上讨论。行替换的源行与目标行不同,缩放系数非零;操作作用于整行,包括所有增广块。有限链中第一步先作用,所以合并左乘矩阵中的初等矩阵按相反次序排列。
把原行记为a1,a2,a3,a4。前三步后各行为a1、a2+α1a1、a3+α2a1、a4+α3a1。最后三步使用这些已经更新的行,因此最终首行为(1+α1β1+α2β2+α3β3)a1+β1a2+β2a3+β3a4,其余三行保持上述形式。这四条原行组合给出G;把六步倒序撤销,则得到H:
G=α1β1+α2β2+α3β3+1α1α2α3β1100β2010β3001
H=1−α1−α2−α3−β1α1β1+1α2β1α3β1−β2α1β2α2β2+1α3β2−β3α1β3α2β3α3β3+1
具体地,依次从第1行减去β3R4、β2R3、β1R2,恢复a1;再从第4、3、2行分别减去α3R1、α2R1、α1R1,恢复全部原行。若最终各行为bj,这就是说a1=b1−β1b2−β2b3−β3b4,且ai+1=bi+1−αia1;展开便恰好得到所写H。所以HG=GH=I4。这些操作全是行替换,没有任何参数必须非零,即使G左上角等于零也不影响结论。
题 9
设实数α1,α2,β1,β2满足β1β2=0。四行矩阵A,B经下面五步操作相连。(a) 给出一个逐项明确的4×4矩阵G,使B=GA,并写出从I4到G的操作链;不必写中间矩阵。(b) 写出从G返回I4的操作链。
1. R3←R3+(α1)R22. R1←(β1)R13. R1↔R44. R1←R1+(α2)R35. R2←(β2)R2
完整解答
以下在实数上讨论。行替换的源行与目标行不同,缩放系数非零;操作作用于整行,包括所有增广块。有限链中第一步先作用,所以合并左乘矩阵中的初等矩阵按相反次序排列。
(a) 设原行为a1,a2,a3,a4。第一步使第3行成为a3+α1a2;缩放第1行再交换第1、4行,使第1行为a4、第4行为β1a1。把当前第3行的α2倍加到第1行,得到a4+α2a3+α1α2a2。最后第2行变为β2a2。因此:
G=000β1α1α2β2α10α20101000
从I4出发的所求操作链,就是按原次序施行题目五步;它们作用于标准行后,逐项得到所写G。(b) 倒序撤销:
1. R2←(β21)R22. R1←R1+(−α2)R33. R1↔R44. R1←(β11)R15. R3←R3+(−α1)R2
前两步逆操作后,第2行为e2、第1行为e4。交换后第4行恢复为e4,第1行成为β1e1;乘倒数恢复e1,最后一步把第3行从e3+α1e2恢复为e3。终点因此为I4。β1,β2非零的条件,恰好用于这两次倒数缩放。
题 10
设m×n矩阵A,A′行等价。证明对每个n×p矩阵B,AB与A′B行等价。
完整解答
以下在实数上讨论。行替换的源行与目标行不同,缩放系数非零;操作作用于整行,包括所有增广块。有限链中第一步先作用,所以合并左乘矩阵中的初等矩阵按相反次序排列。
证明中把内维数和输出列数分别记为p和q;这只是重新命名题目的维数符号。
(a):单步正向传递
直接逐项证明正向恒等式:在行替换情形,只有输出第k行改变:
(A′B)kj=r=1∑p(akr+αair)brj=(AB)kj+α(AB)ij.
其余输出行逐项相同。缩放时(A′B)kj=β(AB)kj,其余行不变。交换A的第i,k行,会恰好交换AB的第i,k行。由此直接证明三种正向传递恒等式;再作逆操作,便得到反向恒等式。
(b):对行变换序列长度归纳
长度为零时,两矩阵相同,乘积也相同。假设结论对长度n成立。对长度n+1的序列,令An为前n步后的矩阵。归纳假设表明,同样的前n步把AB送到AnB;再由(a)已逐项证明的正向恒等式,最后一步把AnB送到A′B。所以整个n+1步序列把AB送到A′B。归纳得任意有限长度的结论。所有乘积大小均为m×q,因此始终是同尺寸矩阵上的合法行变换。