核心内容摘要
ACFun黄化龙年限定：当国潮遇上国漫，一场跨越次元的盛宴
行列式展开的定义拉普拉斯展开定理行列式的展开正式称为拉普拉斯展开Laplace expansion或余子式展开cofactor expansion是计算n阶行列式的一种递归方法。
它允许我们将高阶行列式转化为若干低阶行列式的加权和非常实用尤其在手算时。
核心定义对于一个n阶方阵其行列式 det(A) 可以按任意一行或任意一列展开。
按第i行展开其中 是余子式cofactor是子式minor即删除第i行和第j列后得到的(n-
阶子矩阵的行列式。
按第k列展开类似关键是符号项它决定正负号。
可以记忆为“棋盘格局”从左上角开始交替正负简单例子2×2行列式基础2阶行列式本身就是最简单的展开形式按第一行展开这正是定义本身。
详细例子3×3行列式取一个具体矩阵我们按第一行展开i1现在计算每个2阶子行列式第一项→ 贡献第二项→ 贡献第三项→ 贡献总和−312−90这个矩阵的行列式确实为0因为第三行 第一行×3 - 第二行×2行相关。
如果你按其他行或列展开会得到相同结果行列式性质保证。
图中展示了一个4×4 矩阵的行列式使用拉普拉斯展开Laplace expansion选择按第二行展开最终将 4 阶行列式转化为三个 3 阶行列式的加权和其中一项因元素为 0 而省略。
原矩阵为什么选择按第二行展开第二行元素
25, x₂, 0,
3只有两个非零元素第三项为 0可直接忽略展开后只有三项需要计算实际两项非零系数 一项 0。
这比按其他行/列展开更简洁例如第一行有三个非零元素第四行四个全非零。
展开后得到的 3×3 子矩阵结构较好便于进一步计算图中标注 “LU-decomposition applicable”可能适合 LU 分解来快速求行列式。
拉普拉斯展开公式按第 i2 行其中是删除第 2 行、第 j 列后的3×3 子式minor的行列式。
逐项计算符号和贡献j1元素符号贡献j2元素符号贡献j3元素贡献直接为 0图中省略j4元素符号贡献于是这正是图中所示的展开式。
每个子矩阵是如何得到的删除过程可视化第一项删除第2行、第1列保留列 2,3,4第二项删除第2行、第2列保留列 1,3,4 第一列全为 0除最后一行结构较“稀疏”。
第四项删除第2行、第4列保留列 1,2,3 第一列同样几乎为 0。
为什么标注 “LU-decomposition applicable”观察三个 3×3 子矩阵第
第三个子矩阵的第一列几乎全为 0只有最后一行非零这使得它们非常适合LU 分解无需或仅需少量行交换。
LU 分解后行列式 L 和 U 对角元素乘积计算极快尤其是手工计算时。
第一个子矩阵虽不那么稀疏但整体结构也利于消元。
通过对这三个 3×3 矩阵分别进行 LU 分解或进一步展开最终可得到关于 的线性表达式图中最下方所示的部分其他常数项图中标注 “a part of Expression”说明这只是完整表达式中的线性部分。
图中子矩阵的LU分解过程讲解在图中通过拉普拉斯展开得到的三个3×3子矩阵具有明显的稀疏结构许多零元素第一列往往只有底部非零这使得它们特别适合使用LU分解或带行交换的PLU分解来高效计算行列式。
LU分解的基本原理与计算行列式的优势LU分解将矩阵A分解为A LU或PA LUP为置换矩阵L单位下三角矩阵对角线为1下三角存放消元乘数U上三角矩阵计算行列式的关键无行交换时有s次行交换时过程本质是高斯消元从左到右、从上到下处理每一列k。
若当前主元pivot交换行使主元非零记录交换次数。
对下方的每一行i (k)计算乘数第i行第k行最终得到UL由乘数填充。
优势在图中子矩阵中零元素多交换后往往只需1~2步消元计算量极小非常适合手算。
示例1x₂对应的子矩阵最典型结构最稀疏矩阵B步骤1检查第一主元a₁₁ 0 → 需要行交换。
交换第1行和第3行一次交换符号因子-1交换后步骤2k1消元第一列下方两行第一列已为0无需任何操作l₂₁0, l₃₁0。
步骤3k2消元第二列主元a₂₂
3 ≠0l₃₂
5 /
3 5/3第3行 - (5/
× 第2行第三行新值第2列
5 - (5/
×
3
5 -
5 0第3列0 - (5/
×
5 -5/6得到U计算行列式对角积
15 ×
3 × (-5/
-
0375一次交换 → det(B) (-
× (-
0.
0375)
0375贡献项x₂ ×
0375图中
0375x₂正是由此而来示例
2
3对应的子矩阵结构类似含1-x₁矩阵C步骤1第一主元0→ 交换第1行和第3行一次交换交换后步骤2k1→ 下方已为0无消元。
步骤3k2主元
2 l₃₂ (1-x₁)/
2 5(1-x₁) 第3行 - 5(1-x₁) × 第2行第2列1-x₁ - 5(1-x₁)×
2 1-x₁ - (1-x₁) 0完美抵消第3列
5 - 5(1-x₁)×
3
5 -
5(1-x₁)
5x₁ - 1得到U对角
15,
2,
5x₁ - 1对角积
03(
5x₁ -
1)
045x₁ -
03det(C) (-
× (
045x₁ -
0.
-
045x₁
03贡献项
3 × det(C) -
0135x₁
009示例3-
25对应的子矩阵含1-x₁无需交换矩阵D步骤1第一主元1-x₁假设≠0→ 无需交换。
步骤2~3正常消元两列计算稍多但仍可手算。
消元后U的对角会包含(1-x₁)因子和常数项最终det(D)会产生另一部分x₁线性项。
贡献项-
25 × det(D)提供剩余的x₁系数和常数项。
为什么图中特别标注“LU-decomposition applicable”第
第三个子矩阵只需一次行交换 一到两次消元且常出现完美抵消如第二列直接变0。
第一子矩阵虽稍复杂但整体计算量远低于直接用定义或Sarrus法则展开3阶行列式。
合并三项后自然得到关于x₁、x₂的简单线性表达式图中最下方所示部分。
这种结构优化是典型的手算技巧稀疏零元素 → 少量运算即可得精确或高精度系数。
为什么有用递归性高阶行列式总能降到低阶。
实际计算如果某行/列有很多0选择含0最多的行/列展开能大幅简化很多项直接为0。
理论证明这是行列式公理化定义的等价形式之一。
总结图中巧妙利用第二行只有两个有效非零元素进行拉普拉斯展开将 4 阶问题降为三个 3 阶问题。
子矩阵结构稀疏便于后续 LU 分解或手工计算。
这种展开策略是典型的手算优化技巧总是优先选择含 0 最多或变量最少的行/列展开。