拉普拉斯展开法以余子式与代数余子式为基本构件,将行列式的计算化为低阶行列式的运算。
设 A 是一个 n 阶方阵, 记其行列式为 D=detA:
A=a11a21⋮an1a12a22⋮an2⋯⋯⋱⋯a1na2n⋮ann
按第 k 行展开时,其公式为:
D=j=1∑n(−1)k+jakjMkj
按第 l 列展开时,其公式为:
D=i=1∑n(−1)i+lailMil
即某一行(列)中所有元素与其对应的代数余子式 (−1)i+jMij 的乘积之和。
Mij 代表将 A 的第 i 行,第 j 列元素 aij 所在的行和列删去后,留下的 (n−1) 阶余子式。
(−1)i+j 代表行列式的位置系数。
嵌套展开
拉普拉斯展开可以递归进行:第一次展开得到 n 个 (n−1) 阶余子式,每个余子式本身又是一个行列式,可再按某一行或一列展开,阶数逐次降低,直到化为 2×2 行列式 ad−bc 直接求值。
以三阶行列式为例,先按第一行展开,再把每个二阶余子式写成 ad−bc:
a11a21a31a12a22a32a13a23a33=a11a22a32a23a33−a12a21a31a23a33+a13a21a31a22a32=a11(a22a33−a23a32)−a12(a21a33−a23a31)+a13(a21a32−a22a31).
按列展开同理,先沿第 l 列展开得到 n 个 (n−1) 阶余子式,再继续嵌套降阶。
上三角行列式
拉普拉斯展开可以直接推出行列式性质中上(下)三角行列式的结论:其值等于主对角线上各元素的乘积。
对 n 阶上三角行列式按第一列展开,由于第一列只有 a11 非零,展开只剩一项:
a110⋮0a12a22⋮0⋯⋯⋱⋯a1na2n⋮ann=a11a220⋮0a23a33⋮0⋯⋯⋱⋯a2na3n⋮ann
余子式 M11 仍是 (n−1) 阶上三角行列式,反复按第一列展开即得
a11a22⋯ann=i=1∏naii.
下三角行列式同理,按第一行展开得到同一结果。
步骤
- 选择某一行/列,通常选0最多的那一行/列
- 将该行/列中的每个元素乘以它所对应的余子式,并附以恰当的符号即:(−1)行号+列号
- 将所有同行/列元素的乘积之和作为所求行列式的值
Example
对于一个4阶行列式:
15613−28240579136
按第2行展开,则结果为:
15613−28240579136=−5382457936+(−2)161457936−0161382936+1161382457
拉普拉斯展开法的关键在于将高阶行列式化为低阶行列式的运算,可以很大程度减少计算量。
这种方法在手算时很有用,但在数值计算时更多采用高斯消元法。不过拉普拉斯展开还是揭示了行列式与子式的递推关系,体现了行列式理论的内在结构。
相关内容