xy0v0's Blog

Back

1 从排列展开式理解代数余子式

设 A=(aij)A=(a_{ij}) 为 nn 阶方阵,AijA_{ij} 表示元素 aija_{ij} 的代数余子式。 固定第二列到第 nn 列,把第一列换成变量 x1,…,xnx_1,\ldots,x_n,记 Δ(x1,…,xn)=∣x1a12⋯a1nx2a22⋯a2n⋮⋮⋮xnan2⋯ann∣.\Delta(x_1,\ldots,x_n)= \begin{vmatrix}x_1&a_{12}&\cdots&a_{1n}\\x_2&a_{22}&\cdots&a_{2n}\\ \vdots&\vdots&&\vdots\\x_n&a_{n2}&\cdots&a_{nn}\end{vmatrix}. 按列的次序写出排列定义,有 Δ(x1,…,xn)=∑(i1,…,in)(−1)τ(i1,…,in)xi1ai2,2⋯ain,n,\Delta(x_1,\ldots,x_n) =\sum_{(i_1,\ldots,i_n)}(-1)^{\tau(i_1,\ldots,i_n)} x_{i_1}a_{i_2,2}\cdots a_{i_n,n}, 其中求和遍历 1,…,n1,\ldots,n 的所有排列,τ\tau 表示逆序数。 每一项恰含一个 xkx_k。把含同一个 xkx_k 的项合并,就得到 Δ(x1,…,xn)=∑k=1nxkAk1.\Delta(x_1,\ldots,x_n)=\sum_{k=1}^n x_k A_{k1}. 因此,代数余子式 Ak1A_{k1} 就是把第一列看作变量时,xkx_k 的系数。 它不依赖第一列中的元素,只由其余列决定。换成任意一列或一行,道理相同。

不同行(列)的对应乘积之和

上述观点也解释了上一讲的零和公式。若把第 jj 列换成第 ll 列(l≠jl\ne j), 新行列式有两列相同,故为零;第 jj 列的代数余子式却不变。因此 ∑i=1nailAij=0(l≠j).\sum_{i=1}^n a_{il}A_{ij}=0\quad(l\ne j). 同理,按行有 ∑j=1nakjAij=0(k≠i).\sum_{j=1}^n a_{kj}A_{ij}=0\quad(k\ne i). 与按行(列)展开公式合在一起,即为 ∑r=1nahrAir={det⁡A,h=i,0,h≠i,∑r=1narkArj={det⁡A,k=j,0,k≠j.\sum_{r=1}^n a_{hr}A_{ir} =\begin{cases}\det A,&h=i,\\0,&h\ne i,\end{cases} \qquad \sum_{r=1}^n a_{rk}A_{rj} =\begin{cases}\det A,&k=j,\\0,&k\ne j.\end{cases} 同一行(列)对应相乘求和得到原行列式,换成另一行(列)则得到零。

2 行列式的降阶与递推计算

例题 2.1. 计算 D=∣53−120172520−23100−4−14002350∣.D=\begin{vmatrix} 5&3&-1&2&0\\1&7&2&5&2\\0&-2&3&1&0\\0&-4&-1&4&0\\0&2&3&5&0 \end{vmatrix}.

解. 先按第五列展开,再按所得四阶行列式的第一列展开,得 D=2(−1)2+5∣53−120−2310−4−140235∣=−10∣−231−4−14235∣.D=2(-1)^{2+5} \begin{vmatrix}5&3&-1&2\\0&-2&3&1\\0&-4&-1&4\\0&2&3&5\end{vmatrix} =-10\begin{vmatrix}-2&3&1\\-4&-1&4\\2&3&5\end{vmatrix}. 对右侧三阶行列式作 r3←r3+r1r_3\leftarrow r_3+r_1,再从第三行提出 66,得 D=−60∣−231−4−14011∣=−60[(−2)(−1−4)−3(−4)+(−4)]=−1080.D=-60\begin{vmatrix}-2&3&1\\-4&-1&4\\0&1&1\end{vmatrix} =-60\bigl[(-2)(-1-4)-3(-4)+(-4)\bigr] =\boxed{-1080}. 每次展开都选择非零元素最少的一行或一列;行倍加用于进一步制造零元素。 ◻

例题 2.2. 求主对角线为 a+ba+b、上副对角线为 aa、下副对角线为 bb 的三对角行列式 Dn=∣a+ba0⋯0ba+ba⋱⋮0ba+b⋱0⋮⋱⋱⋱a0⋯0ba+b∣(n≥1).D_n=\begin{vmatrix} a+b&a&0&\cdots&0\\b&a+b&a&\ddots&\vdots\\ 0&b&a+b&\ddots&0\\\vdots&\ddots&\ddots&\ddots&a\\0&\cdots&0&b&a+b \end{vmatrix}\quad(n\ge1).

解. 约定 D0=1D_0=1,则 D1=a+bD_1=a+b。按第一列展开,第二项的余子式再按第一行展开,得 Dn=(a+b)Dn−1−abDn−2(n≥2).D_n=(a+b)D_{n-1}-abD_{n-2}\quad(n\ge2). 令 En=Dn−bDn−1E_n=D_n-bD_{n-1},则 En=aEn−1E_n=aE_{n-1},且 E1=aE_1=a,故 En=anE_n=a^n。 于是 Dn=an+bDn−1D_n=a^n+bD_{n-1},逐次代入得 Dn=∑k=0nan−kbk={an+1−bn+1a−b,a≠b,(n+1)an,a=b.\boxed{D_n=\sum_{k=0}^n a^{n-k}b^k =\begin{cases}\dfrac{a^{n+1}-b^{n+1}}{a-b},&a\ne b,\\[6pt] (n+1)a^n,&a=b.\end{cases}} 有限和的形式不需要除以 a−ba-b,因而也适用于 a=ba=b 或参数为零的情形。 ◻

3 多行(列)的拉普拉斯展开

子式、余子式与展开公式

设 AA 为 nn 阶方阵,1≤r<n1\le r<n。选取递增的行指标集 I={i1<⋯<ir}I=\{i_1<\cdots<i_r\} 和列指标集 J={j1<⋯<jr}J=\{j_1<\cdots<j_r\}。 记 A[I,J]A[I,J] 为这些行、列交叉位置构成的子矩阵,Ic,JcI^c,J^c 为相应的补集。 所有子矩阵中的行、列均保持原来的顺序。

定义 3.1. rr 阶行列式 det⁡A[I,J]\det A[I,J] 称为一个 rr 阶子式。 删去所选行、列后所得的行列式,称为它的余子式,即 det⁡A[Ic,Jc]\det A[I^c,J^c]。 其代数余子式为 (−1)∑I+∑Jdet⁡A[Ic,Jc],∑I=i1+⋯+ir,∑J=j1+⋯+jr.(-1)^{\sum I+\sum J}\det A[I^c,J^c], \qquad \sum I=i_1+\cdots+i_r,\quad\sum J=j_1+\cdots+j_r.

定理 3.2 (拉普拉斯展开). 固定 rr 列 JJ,对所有 rr 行的选择 II 求和,有 det⁡A=∑|I|=r(−1)∑I+∑Jdet⁡A[I,J]det⁡A[Ic,Jc].\boxed{\det A=\sum_{|I|=r}(-1)^{\sum I+\sum J} \det A[I,J]\det A[I^c,J^c].} 固定 rr 行 II 时,则对所有 rr 列的选择 JJ 求和: det⁡A=∑|J|=r(−1)∑I+∑Jdet⁡A[I,J]det⁡A[Ic,Jc].\det A=\sum_{|J|=r}(-1)^{\sum I+\sum J} \det A[I,J]\det A[I^c,J^c].

每次固定一组行或列,只遍历另一组,共有 (nr)\binom nr 项。 当 r=1r=1 时,这正是上一讲的按一行或一列展开公式。

按前 rr 列展开时,符号从哪里来

固定 J={1,…,r}J=\{1,\ldots,r\}。排列展开式的每一项在前 rr 列取出的元素, 必定来自 rr 个不同的行。按这组行 I={k1<⋯<kr}I=\{k_1<\cdots<k_r\} 把所有项分组。

先把所选行排在前面、其余行排在后面,且两组内部均递增。 第 ksk_s 行前面有 ks−1k_s-1 个较小的行指标,其中 s−1s-1 个已在所选组内, 故它与另一组产生 ks−sk_s-s 个逆序。组间的总逆序数为 ∑s=1r(ks−s)=∑s=1rks−r(r+1)2.\sum_{s=1}^r(k_s-s)=\sum_{s=1}^r k_s-\frac{r(r+1)}2. 组内的排列分别构成两个行列式,因此这一组的和为 (−1)∑I+r(r+1)/2det⁡A[I,{1,…,r}]det⁡A[Ic,{r+1,…,n}].(-1)^{\sum I+r(r+1)/2} \det A[I,\{1,\ldots,r\}]\det A[I^c,\{r+1,\ldots,n\}]. 这里把减号改成加号不改变奇偶性。对所有 II 求和即得前 rr 列的公式。 一般列集 JJ 先移到前面,另产生 ∑J−r(r+1)/2\sum J-r(r+1)/2 次对换, 故总符号为 (−1)∑I+∑J(-1)^{\sum I+\sum J}。按行展开可由转置得到。

4 拉普拉斯展开的应用

分块三角行列式

命题 4.1. 设 AA、DD 分别是 pp 阶、qq 阶方阵,BB 为 p×qp\times q 矩阵,则 det⁡(AB0D)=det⁡Adet⁡D.\det\begin{pmatrix}A&B\\0&D\end{pmatrix}=\det A\det D. 同理,det⁡(A0CD)=det⁡Adet⁡D\det\begin{pmatrix}A&0\\C&D\end{pmatrix}=\det A\det D。

证明. 按前 pp 列展开。这些列在后 qq 行全为零,因此只有选取前 pp 行时子式才可能非零。 该项的符号为 (−1)2(1+⋯+p)=1(-1)^{2(1+\cdots+p)}=1,子式为 det⁡A\det A,余子式为 det⁡D\det D。 下分块三角的情形由转置得到。 ◻

这里要求对角块是方阵,但不要求它们可逆。一般的分块矩阵不能直接套用下面的二阶形式: det⁡(ABCD)=det⁡Adet⁡D−det⁡Bdet⁡C.\det\begin{pmatrix}A&B\\C&D\end{pmatrix} =\det A\det D-\det B\det C. 例如四个块均为 2×22\times2 矩阵,取 A=D=0A=D=0、B=C=I2B=C=I_2。 整体矩阵对应两次对换,行列式为 11,而右端为 −1-1。

一个四阶行列式

例题 4.2. 从矩阵 (a1b1c1a2b2c2)\begin{pmatrix}a_1&b_1&c_1\\a_2&b_2&c_2\end{pmatrix} 中分别删去第 1,2,31,2,3 列, 所得二阶行列式依次记为 A=∣b1c1b2c2∣,B=∣a1c1a2c2∣,C=∣a1b1a2b2∣.A=\begin{vmatrix}b_1&c_1\\b_2&c_2\end{vmatrix},\quad B=\begin{vmatrix}a_1&c_1\\a_2&c_2\end{vmatrix},\quad C=\begin{vmatrix}a_1&b_1\\a_2&b_2\end{vmatrix}. 证明 H=∣a1b1c10a2b2c200a1b1c10a2b2c2∣=AC−B2.H=\begin{vmatrix}a_1&b_1&c_1&0\\a_2&b_2&c_2&0\\ 0&a_1&b_1&c_1\\0&a_2&b_2&c_2\end{vmatrix}=AC-B^2.

证明. 按前两行展开。含第四列的二阶子式均为零;选取第 2,32,3 列时,余子式有一列为零。 故六种列选择中只有 {1,2}\{1,2\}、{1,3}\{1,3\} 两项可能非零,分别为 (−1)1+2+1+2CA=CA,(−1)1+2+1+3BB=−B2.(-1)^{1+2+1+2}CA=CA,\qquad (-1)^{1+2+1+3}BB=-B^2. 相加即得 H=AC−B2H=AC-B^2。此处 A,B,CA,B,C 均是二阶行列式的值,即三个数。 ◻

5 代数余子式与三维叉积

设 𝜶=(a11,a21,a31)𝖳\boldsymbol\alpha=(a_{11},a_{21},a_{31})^{\mathsf T}、 𝜷=(a12,a22,a32)𝖳\boldsymbol\beta=(a_{12},a_{22},a_{32})^{\mathsf T} 为 ℝ3\mathbb R^3 中的向量。 把它们作为前两列,第三列记为 𝒙=(x,y,z)𝖳\boldsymbol x=(x,y,z)^{\mathsf T}。按第三列展开: det⁡(𝜶,𝜷,𝒙)=∣a11a12xa21a22ya31a32z∣=x∣a21a22a31a32∣−y∣a11a12a31a32∣+z∣a11a12a21a22∣.\begin{align*} \det(\boldsymbol\alpha,\boldsymbol\beta,\boldsymbol x) &=\begin{vmatrix}a_{11}&a_{12}&x\\a_{21}&a_{22}&y\\a_{31}&a_{32}&z\end{vmatrix}\\ &=x\begin{vmatrix}a_{21}&a_{22}\\a_{31}&a_{32}\end{vmatrix} -y\begin{vmatrix}a_{11}&a_{12}\\a_{31}&a_{32}\end{vmatrix} +z\begin{vmatrix}a_{11}&a_{12}\\a_{21}&a_{22}\end{vmatrix}. \end{align*} 因此,若令 𝜸=(a21a32−a31a22a31a12−a11a32a11a22−a21a12),\boldsymbol\gamma= \begin{pmatrix} a_{21}a_{32}-a_{31}a_{22}\\ a_{31}a_{12}-a_{11}a_{32}\\ a_{11}a_{22}-a_{21}a_{12} \end{pmatrix}, 就有 det⁡(𝜶,𝜷,𝒙)=𝜸⋅𝒙.\det(\boldsymbol\alpha,\boldsymbol\beta,\boldsymbol x) =\boldsymbol\gamma\cdot\boldsymbol x. 取 𝒙=𝜶\boldsymbol x=\boldsymbol\alpha 或 𝜷\boldsymbol\beta,左侧都有两列相同,因而 𝜸⋅𝜶=0,𝜸⋅𝜷=0.\boldsymbol\gamma\cdot\boldsymbol\alpha=0, \qquad \boldsymbol\gamma\cdot\boldsymbol\beta=0. 这解释了为什么由这三个代数余子式组成的向量,同时与原来的两个向量正交。

定义 5.1. 上述向量 𝜸\boldsymbol\gamma 称为 𝜶\boldsymbol\alpha 与 𝜷\boldsymbol\beta 的叉积,记为 𝜶×𝜷\boldsymbol\alpha\times\boldsymbol\beta。

例如 (1,0,0)𝖳×(0,1,0)𝖳=(0,0,1)𝖳.(1,0,0)^{\mathsf T}\times(0,1,0)^{\mathsf T}=(0,0,1)^{\mathsf T}. 当 𝜶,𝜷\boldsymbol\alpha,\boldsymbol\beta 线性无关时,叉积非零,给出它们张成平面的一个法向量。 当两向量线性相关时,这些二阶子式均为零,叉积为零向量,不能用来指定法向方向。 交换两向量,相当于交换行列式的前两列,故 𝜷×𝜶=−𝜶×𝜷.\boldsymbol\beta\times\boldsymbol\alpha=-\boldsymbol\alpha\times\boldsymbol\beta. 所以叉积的三个坐标并非另行记忆的一组公式,它们就是第三列的代数余子式,符号依次为 +,−,++,-,+。

阅读 / 下载 PDF

高等代数(一)1.5:行列式计算与拉普拉斯展开
https://www.xy0v0.top/archives/higher-algebra-i-01-05-determinant-calculation-and-laplace-expansion
Author xy0v0
Published at 九月 30, 2026