克拉默法则、体积和线性变换
行列式可用于克拉默法则求解线性方程组,但仅适用于非零且规模较小的情况。行列式的几何意义体现在:二阶行列式可表示平行四边形的面积,三阶行列式可表示平行六面体的体积。
1. 克拉默法则
行列式不仅可以用于判断矩阵的可逆性,还可用于直接求解线性方程组的解。下面要介绍的克拉默法则 (Cramer’s Rule) 就是一种基于行列式求解线性方程组的方法,它可通过计算不同行列式的比值来确定方程组中的每个未知数。具体法则如下:
定理 7
克拉默法则
设
A 是一个
n×n 的可逆矩阵。对于任意向量
b∈Rn ,线性方程组
Ax=b 的唯一解的第
i 个分量由下式给出:
xi=detAdetAi(b),i=1,2,…,n(1) 其中,符号 Ai(b) 表示一个通过替换矩阵 A 的某一列而得到的新矩阵。若 n 阶矩阵 A=[a1,…,ai,…,an] ,那么 Ai(b)=[a1,…,b,…,an] 。
符号定义矩阵
A 的列记为
(a1,a2,…,an)单位矩阵
I 的列记为
(e1,e2,…,en)向量
x 是矩阵方程方程
Ax=b 的解
矩阵乘法的展开
将
A 乘以一个替换列后的矩阵
Ii(x) ,这个矩阵是单位矩阵
I 的第
i 列替换为向量
x :
A(Ii(x))=A[e1⋯x⋯en]=[Ae1…Ax…Aen]=[a1…b…an]=Ai(b)
行列式的乘法性质根据行列式的乘法性质,有:
det(A)⋅det(Ii(x))=det(Ai(b)) 由于
Ii(x) 是一个替换列后的单位矩阵,其行列式值为向量
xi ,上式可以进一步化简:
det(A)xi=det(Ai(b))⇒xi=det(A)det(Ai(b)) 定理证毕!
使用克拉默法则来解线性方程组,请看下面示例:
2. 克拉默法则求矩阵的逆
克拉默法则可以用于直接求解逆矩阵的
(i,j) 元素,其推导过程如下: 给定一个
n×n 的方阵
A ,其逆矩阵
A−1 的第
j 列可以看做是一个向量
x ,满足:
Ax=ej 其中
ej 是单位矩阵的第
j 列,
x 的第
i 个元素是
A−1 中的
(i,j) 元素。
从求解线性方程组的角度出发,根据克拉默法则可得:
(A−1)ij=xi=detAdetAi(ej)(2) Ai(ej) 是用单位矩阵
In 的第
j 列替换矩阵
A 的第
i 行得到的,它的行列式
detAi(ej) 是矩阵
A 按第
j 行第
i 列的余子展开式,有如下表达式:
detAi(ej)=(−1)i+jdetAji=Cji(3) 这个推导过程可以结合下面的动画演示来辅助理解:
由
(2),(3) 可得
A−1 :
A−1=detA1C11C12⋮C1nC21C22⋮C2n⋯⋯⋱⋯Cn1Cn2⋮Cnn(4) 公式
(4) 右边的代数余子式的矩阵又被称为
A 的
伴随矩阵 (adjugate) 或
经典伴随矩阵 (classical adjoint) ,记作
adjA 。对于
(4) 式,有如下定理:
对
(5) 式两边同时乘以矩阵
A ,再做简单推导,可得如下结论:
adj(A)⋅A=det(A)⋅In(6) 公式
(6) 可以在已知矩阵
A 和伴随矩阵
adjA 的情况计算行列式
detA 。 不过,如果使用定理
8 来求解矩阵的逆,计算复杂度较高,它只适合处理小规模矩阵。例如,对于
3×3 的矩阵,使用该方法所需的计算量远高于行化简算法。请观察下面的计算过程:
克拉姆法则的局限性
克拉默法则虽然直观,但它的计算复杂度高,而且只能处理非奇异矩阵,且不适用于大规模矩阵,因此在实际应用中较少用于高维或复杂的线性方程组求解或求逆矩阵。在工程应用中,常使用
LU 分解、高斯消元法或最小二乘法等更高效的方法。
3. 行列式的几何意义
将行列式与几何意义联系起来,能为我们提供更直观的理解:它不仅是一个代数工具,还能用于描述由矩阵列向量围成的区域面积或体积。有如下定理:
定理 9
行列式与面积或体积的关系
(1) 如果
A 是一个
2×2 矩阵,那么由矩阵
A 的列向量所确定的平行四边形的面积为
∣detA∣ 。
(2) 如果
A 是一个
3×3 矩阵,那么由矩阵
A 的列向量所确定的平行六面体的面积为
∣detA∣ 。
3.1 二维空间中的行列式与面积
二阶矩阵如果是三角矩阵,我们可以很容易看出由其列向量围成的平行四边形的面积就是矩阵对应的行列式,如下所示:
要证明定理 9(1) ,我们可以先把任意矩阵通过行变换(或列变换)化简为三角矩阵,只要能证明原矩阵的列向量和三角矩阵的列向量围成的平行四边形面积相同即可。对任意可逆的二阶矩阵 A=[a1a2] 经过列倍加变换后,总能得到三角矩阵 A′=[a1a2+ca1] 。下面只需要证明 A 和 A′ 的列向量围成的平行四边形面积是相同的即可,下面的动画过程展示了我们总能找到一个合适的标量 c ,可以使得 a2+ca1 垂直于 a1 ,从而很容易观察到两个平行四边形的面积相等:
证明命题给定两个非零向量
a1,a2 证明对于任意标量
c ,由
a1 和
a2+ca1 围成的平行四边形的面积与由
a1 和
a2 围成的平行四边形面积相等。
a1,a2 共线的情况 a1,a2 共线,平行四边形退化为一条直线,面积为
0 。
几何解释定义平行线
L: 假设
L 是过原点
0 和向量
a1 的直线。
向量
a2+ca1 是在
a2 的基础上沿
a1 方向平移。即,
a2+ca1 在平行于
L 的直线上。
与 L 的垂直距离
a2 和
a2+ca1 相对直线
L 的垂直距离相等。观察动画不难发现从
a2+ca1 到直线
L 的垂直距离不会因为沿
a1 方向的滑动而改变。
结论:两个平行四边形面积相等两个平行四边形的底边相同(即
a1 ),并且它们垂直高度相同。因此这两个平行四边形的面积相等。
3.2 三维空间中的行列式与体积
三阶三角矩阵的列向量围成的平行六面体的体积等于对角线元素的乘积(取绝对值)。只要三角矩阵的对角线元素保持不变,该平行六面体的体积就不会发生变化。请观察下面的动画过程:
v1v2v3 2.000.000.000.003.000.000.000.004.00在上面的示例中,向量的数值发生变化时,由向量 x,y 确定的平行四边形的面积(底)以及向量 z 到 xy 平面的垂直距离(高)始终不变,所以动画过程中平行六面体的体积并没有发生变化。证明定理 9(2) 和证明定理 9(1) 的思路是一样的,下面我们用动画来演示证明的过程:
v1v2v3 3.001.002.005.000.001.004.001.001.00任何 3×3 矩阵都可以通过不改变 ∣detA∣ 的列倍加变换转化为对角矩阵。在变换过程中由矩阵列向量构成的六面体的体积始终保持不变,上面的动画展示了一次列倍加变换的过程:[a1a2a3]→[a1+c1a2+c2a3a2a3]。当把矩阵化简为对角矩阵时,列向量构成的会是一个轴对齐的长方体。
4.1 行列式作为缩放因子
行列式还可以用来揭示了线性变换如何影响空间的面积、体积以及方向。我们在 1.8 线性变换的矩阵一节介绍过一些常见的变换矩阵( R2 空间)。我们可以观察到旋转矩阵、剪切矩阵的行列式都为 1 ,变换后图形的面积不会发生变化;缩放矩阵的行列式为 k (k∈R2) ,变换后图形的面积会缩放 ∣k∣ 倍(如果是投影矩阵, k=0 )。行列式的绝对值描述了线性变换对区域面积 (R2) 或体积 (R3) 的缩放比例,而符号则反映了变换的方向(如:是否翻转)。这种几何视角帮助我们从形状和大小的变化来理解矩阵变换的本质。注意观察下面变换矩阵的行列式和变换图形的面积变化:
给定一个线性变化 T ,定义在集合 S 上, T(S) 表示集合 S 中所有点在变换 T 下的像。我们关注的是变换后区域的面积或体积与原始区域的变化。有如下定理:
定义原点处的平行四边形考虑
2×2 的情况,设
A=[a1a2]。在
R2 中,由向量
b1 和
b2 确定的、位于原点的平行四边形
S 可以表示为:
S={s1b1+s2b2:0≤s1≤1,0≤s2≤1} 计算 S 在 T 下的像集合
S 在线性变换
T 下的像由以下形式的点组成:
T(s1b1+s2b2)=s1T(b1)+s2T(b2)=s1Ab1+s2Ab2 其中
0≤s1≤1 ,
0≤s2≤1 。
将 T(S) 表示为矩阵乘积由上一步可知,
T(S) 是由矩阵
[Ab1Ab2]的列向量确定的平行四边形。这个矩阵可以写成
AB 的形式,其中
B=[b1b2]。
应用行列式乘法定理根据定理
9 (行列式的乘法性质),有:
{T(S) 的面积}=∣detAB∣=∣detA∣∣detB∣=∣detA∣⋅{S 的面积} 推广到任意平行四边形任意平行四边形可以表示为
p+S 的形式,其中
p 是一个向量,
S 是位于原点的平行四边形。 由于平移不改变集合的面积,有:
{T(p+S) 的面积}={T(p)+T(S) 的面积}={T(S) 的面积}=∣detA∣⋅{S 的面积}=∣detA∣⋅{(p+S) 的面积}(平移不改变面积)(由第4步)(平移不改变面积) 结论这证明了定理对
R2 中所有平行四边形都成立。证毕!
定义原点处的平行六面体考虑
3×3 的情况,设
A=[a1a2a3]。在
R3 中,由向量
b1 、
b2 和
b3 确定的、位于原点的平行六面体
S 可以表示为:
S={s1b1+s2b2+s3b3:0≤s1,s2,s3≤1} 计算 S 在 T 下的像集合
S 在线性变换
T 下的像由以下形式的点组成:
T(s1b1+s2b2+s3b3)=s1T(b1)+s2T(b2)+s3T(b3)=s1Ab1+s2Ab2+s3Ab3 其中
0≤s1,s2,s3≤1 。
将 T(S) 表示为矩阵乘积由上一步可知,
T(S) 是由矩阵
[Ab1Ab2Ab3]的列向量确定的平行六面体。这个矩阵可以写成
AB 的形式,其中
B=[b1b2b3]。
应用行列式乘法定理根据定理
9 (行列式的乘法性质),有:
{T(S) 的体积}=∣detAB∣=∣detA∣∣detB∣=∣detA∣⋅{S 的体积} 推广到任意平行六面体任意平行六面体可以表示为
p+S 的形式,其中
p 是一个向量,
S 是位于原点的平行六面体。 由于平移不改变集合的体积,有:
{T(p+S) 的体积}={T(p)+T(S) 的体积}={T(S) 的体积}=∣detA∣⋅{S 的体积}=∣detA∣⋅{(p+S) 的体积}(平移不改变体积)(由第4步)(平移不改变体积) 结论这证明了定理对
R3 中所有平行六面体都成立。证毕!
4.2 定理 10 的推广
定理 10 告诉我们:线性变换对平行四边形或平行六面体的面积/体积的影响,就是简单地乘以变换矩阵行列式的绝对值。但如果我们想把这个结论推广到任意形状的区域(比如圆形、椭圆或其他不规则图形),该怎么办呢?
核心思想是用小方块去逼近任意区域。对于一个圆形区域 R ,我们可以用一堆小正方形去填充它。当这些小正方形足够小时,它们的面积之和就能无限接近圆的真实面积——这正是微积分中计算面积的基本思想。
现在,当我们对这个区域施加一个线性变换
T 时,每一个小正方形都会被变换成一个小平行四边形。根据定理
10 ,每个小平行四边形的面积都是原来小正方形面积的
∣detA∣ 倍。由于
每一个小块的面积都被缩放了相同的倍数 ∣detA∣ ,那么所有小块加起来的总面积自然也被缩放了
∣detA∣ 倍。当我们让小方块无限细分时,就得到了定理
10 的推广形式:
{T(R) 的面积}=∣detA∣⋅{R 的面积} 这个结论对于
R2 中任意有限面积的区域都成立,同样的思路也可以推广到
R3 中的体积计算。这就是为什么行列式在多元微积分的换元积分(雅可比行列式)中扮演着如此重要的角色。
定理
10 的结论对于
R2 中任意有限面积的区域,或
R3 中任意有限体积的区域都成立。
下面通过一个经典的例子来展示定理 10 推广形式的应用——计算椭圆的面积。
例5:设
a 和
b 是正数,求由椭圆方程
a2x12+b2x22=1 所围成的区域
E 的面积。
解:我们可以把椭圆
E 看作是单位圆盘
D 在某个线性变换
T 下的像。设变换矩阵为
A=[a00b] ,若
u=[u1u2] ,
x=[x1x2] , 且
x=Au ,那么:
u1=ax1,u2=bx2 当
u 在单位圆盘
D 内时(即
u12+u22≤1 ),对应的
x 满足
(x1/a)2+(x2/b)2≤1 ,正好落在椭圆
E 内。因此椭圆
E 就是单位圆盘
D 在变换
T 下的像,即
E=T(D) 。根据定理
10 的推广:
{椭圆 E 的面积}={T(D) 的面积}=∣detA∣⋅{D 的面积}=ab⋅π(1)2=πab