Skip to content

算法Algorithm

行化简

施工中 · 内容逐步修订

Row reduction

用初等行变换把矩阵化为阶梯形以求解线性方程组和判定秩。

形式陈述 ​

域 F 上矩阵的三种初等行变换是交换两行、把一行乘以非零标量、把一行的标量倍加到另一行。每步等价于左乘一个可逆初等矩阵。把同一行变换同时作用于系数矩阵与右端项,也就是作用于增广矩阵 [A∣b],便得到与 Ax=b 解集相同的方程组。

反复变换可得到行阶梯形:零行位于底部,每个非零行的首个非零元称为主元,下一行的主元位于上一行主元的右侧。继续使主元为 1,并将每个主元列的其他位置化为 0,就得到简化行阶梯形(RREF)。每个矩阵的 RREF 唯一,主元数等于矩阵的秩;对相容方程组,系数矩阵的非主元列对应自由变量。

两个同型矩阵行等价,当且仅当它们有相同 RREF,也等价于行空间相同。

一个有限的 Gauss–Jordan 过程 ​

输入是可以执行精确域运算的 m×n 矩阵 B,输出是它的 RREF。维护下一主元行 r=1,从左到右扫描列 c=1,…,n:

  1. 在第 r 行及以下寻找 Bic≠0。若没有,就跳过本列;若有,可固定选择行号最小的一项
  2. 交换第 i,r 行,再把第 r 行除以它在第 c 列的非零主元
  3. 对每个 j≠r,从第 j 行减去 Bjc 倍的主元行,使本列除主元外全为零
  4. 令 r←r+1;所有行都已成为主元行时结束,否则继续扫描下一列

每个主元行在此前主元列上已经为零,所以后续操作不会破坏已完成的主元列。列数和行数有限,过程必终止;零矩阵直接跳过全部列。输入增广矩阵时,要保留系数列与右端列的分界:右端列上的孤立主元表示矛盾行,不是一个未知量的主元。

直觉

行变换是在不改变方程约束集合的前提下重写方程:把冗余和依赖逐步显露,最终让主变量与自由变量一目了然。

行化简的三类变换与主元阶梯

例如“第二个方程减去第一个方程的两倍”可以用相反的加法还原。每一步都可逆,原方程组与化简后的方程组便表达同一组约束。出现零行意味着某个方程可由其他方程组合得到;出现主元则留下一个独立约束。

例子与边界

从主元读出解 ​

对增广矩阵

(12|324|6)⟶(12|300|0),

第二行减去第一行的两倍后成为零行,系数矩阵只有一个主元,秩为 1。第二个变量自由:令 y=t,第一行给出 x=3−2t,全部解为 (3−2t,t)。

若把原来第二行的右端项由 6 改为 7,同一步消元就得到 (0,0∣1),即 0=1,所以无解。一般地,增广矩阵中的行 [0 ⋯ 0∣c] 在 c≠0 时给出这种矛盾。

行变换如何改变矩阵 ​

行变换保持行空间,也保持列之间的线性依赖关系,因此保持秩。列向量本身则会随左乘的初等矩阵移动。例如把 (1010) 的第二行减去第一行,列空间就从 span{(1,1)} 变为 span{(1,0)},维数仍为 1。

对方阵,交换两行使行列式变号,一行乘以 c≠0 使行列式乘以 c,一行加上另一行的倍数则保持行列式。这样可以在消元过程中跟踪行列式。相似变换 P−1AP 则表示同一算子的换基,保持特征值;单独做行变换没有这项性质,例如在实数域上,把单位矩阵的一行乘以 2 就改变了特征值。

列变换对应未知变量的坐标替换。例如交换系数矩阵的两列,就要同时交换这两个变量的位置。用列变换求解时,把这些替换记录下来,才能将最后的坐标还原为原未知量。

精确消元与浮点消元 ​

精确代数中,每个非零主元都可以用来做除法。浮点消元还要考虑尺度:除以很小的主元可能产生很大的消元倍数,放大舍入误差。选主元与增长因子用于分析并控制这种误差传播。

推论与应用

行化简统一处理线性方程组的可解性、参数解与秩。把自由变量依次取为标准坐标,可以从齐次方程的参数解中读出核的一组基;把可逆矩阵的增广块 [A∣I] 化成 [I∣B],则得到逆矩阵 B=A−1。

完整 RREF 适合展示规范形和参数解。大型浮点方阵求解通常使用带主元的 Gaussian 消元或 LU 因子,再以回代求解;同一系数矩阵对应多个右端项时,保存的因子还可以重复使用。

在Horton圈基算法中,每条候选圈是一行二元边向量,消元只需异或主元行。剩余为零表示该圈已在旧张成内;非零剩余给新主元。用于消元的剩余向量可能是多个圈的并,输出时仍须保留被接受的原简单圈及其边ID。

参考资料
  • Gilbert Strang, Introduction to Linear Algebra, 6th ed., Wellesley-Cambridge Press, 2023,Chs. 2–3, elimination, echelon forms, and complete solutions of Ax=b。
  • Michael Artin, Algebra, 2nd ed., Pearson, 2011,Ch. 1, row operations and reduced echelon matrices。
关系图谱30 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系