跳转到内容

3.3 克拉默法则、体积和线性变换

[美] David C. Lay、Steven R. Lay 等著,刘深泉等译机械工业出版社《线性代数及其应用(原书第5版)》大学本科 · 教材双语助读 (Alt+Y)

本节应用前面几节的理论得到一些重要且有理论意义的公式以及行列式的几何解释.

克拉默法则在各种理论计算中是必需的.例如,它被用来研究 Ax = b 的解受 ±bb 中元素的变化而受到什么样的影响.然而,这个公式对手工运算是没有多大效果的,除非是 2 × 2 或 3 × 3 矩阵: 对任意 n × n 矩阵 A 和任意的 ℝⁿ 中向量 ±bb,令 Aᵢ(±bb) 表示 A 中第 i 列由向量 ±bb 替换得到的矩阵:

A _ i (±b b) = [ ±b a _ 1 … ±b b … ±b a _ n ]
定理 7 克拉默法则

设 A 是一个可逆的 n × n 矩阵,对 Rⁿ 中任意向量 b,方程 Ax = b 的唯一解可由下式给出:

x _ i = frac det A _ i (boldsymbol b)det A, i = 1, 2, … , n (1)
证明

用 a₁, …, aₙ 表示 A 的列,用 e₁, …, eₙ 表示 n × n 单位阵 I 的列。若 Ax = b ,则由矩阵乘法的定义有

r l A · I _ i (boldsymbol x) & = A [ boldsymbol e _ 1 … boldsymbol x … boldsymbol e _ n ] = [ A boldsymbol e _ 1 … A boldsymbol x … A boldsymbol e _ n ] & = [ boldsymbol a _ 1 … boldsymbol b … boldsymbol a _ n ] = A _ i (boldsymbol b)

由行列式的乘法性质,

(det A) (det I _ i (boldsymbol x)) = det A _ i (boldsymbol b)

左边第二个行列式为 xᵢ (沿第 i 行作余因子展开),从而 (det A)xᵢ = det Aᵢ(±bb) 由于 A 可逆,从而 operatornamedet A ≠ 0 ,于是证明了(1)式.

例 1

利用克拉默法则解方程组

r l & 3 x _ 1 - 2 x _ 2 = 6 & - 5 x _ 1 + 4 x _ 2 = 8
解

视此方程组为 Ax=b 型,利用上面引入的记号,

A = [ c c 3 & - 2 - 5 & 4 ], A _ 1 (boldsymbol b) = [ c c 6 & - 2 8 & 4 ], A _ 2 (boldsymbol b) = [ c c 3 & 6 - 5 & 8 ]

由于 operatornamedet A = 2 ,故此方程组有唯一解。由克拉默法则,有

x _ 1 = frac det A _ 1 (±b b)det A = frac 2 4 + 1 62 = 2 0x _ 2 = frac det A _ 2 (±b b)det A = frac 2 4 + 3 02 = 2 7

许多工程上的问题(特别是在电子工程和控制论中)能用拉普拉斯变换进行分析。这种技巧将一个适当的线性微分方程组转变为一个线性代数方程组,它的系数含有一个参数。下一个例子说明可能出现的代数方程组类型。

例 2

考虑下列方程组,其中 s 是一个未定的参数。确定 s 的值,使得这个方程组有唯一解,并利用克拉默法则写出这个解。

l 3 s x _ 1 - 2 x _ 2 = 4 - 6 x _ 1 + s x _ 2 = 1
解

视此方程组为 Ax = b 型,则

A = [ c c 3 s & - 2 - 6 & s ], A _ 1 (boldsymbol b) = [ c c 4 & - 2 1 & s ], A _ 2 (boldsymbol b) = [ c c 3 s & 4 - 6 & 1 ]

由于 operatornamedet A = 3s² - 12 = 3(s + 2)(s - 2) ,当 s ≠ ± 2 时,这个方程组有唯一解。对这样的 s ,方程组的解为 (x₁, x₂) ,其中

l x _ 1 = frac det A _ 1 (±b b)det A = frac 4 s + 23 (s + 2) (s - 2) x _ 2 = frac det A _ 2 (±b b)det A = frac 3 s + 2 43 (s + 2) (s - 2) = frac s + 8(s + 2) (s - 2)

克拉默法则可以容易地导出一个求 n × n 矩阵 A 的逆的一般公式.A⁻1 的第 j 列是一个向量 x,满足

A boldsymbol x = boldsymbol e _ j

其中 eⱼ 是单位矩阵的第 j 列,x 的第 i 个元素是 A⁻1 中的 (i, j) 元素。由克拉默法则,

{A ^ - 1 中 (i, j) 元素 } = x _ i = frac det A _ i (±b e _ j)det A (2)

回顾 Aⱼi 表示 A 的子矩阵,它由 A 去掉第 j 行和第 i 列得到。 Aᵢ(±beⱼ) 按第 i 列的余因子展开式为

det A _ i (boldsymbol e _ j) = (- 1) ^ i + j det A _ j i = C _ j i (3)

其中 Cⱼi 是 A 的一个余因子.由(2), A⁻1 的 (i, j) 元素等于余因子 Cⱼi 除以 det A 。(注意: Cⱼi 的下标是 (i, j) 的颠倒。)于是

A ^ - 1 = frac 1det A [ c c c c C _ 1 1 & C _ 2 1 & … & C _ n 1 C _ 1 2 & C _ 2 2 & … & C _ n 2 vdots & vdots & & vdots C _ 1 n & C _ 2 n & … & C _ n n ] (4)

公式(4)右边的余因子的矩阵称为 A 的伴随矩阵,记为 adj A。(伴随这个术语在后面的线性变换课程中还有另一层意思。)下一个定理简单复述(4)式.

定理 8 逆矩阵公式

设 A 是一个可逆的 n × n 矩阵,则 A⁻1 = (1)/(det A) operatornameadj A .

例 3

求矩阵 A = 2 & 1 & 3 1 & -1 & 1 1 & 4 & -2 的逆.

解

9 个余因子为

C _ 1 1 = + | c c - 1 & 1 4 & - 2 | = - 2, C _ 1 2 = - | c c 1 & 1 1 & - 2 | = 3, C _ 1 3 = + | c c 1 & - 1 1 & 4 | = 5C _ 2 1 = - | c c 1 & 3 4 & - 2 | = 1 4, C _ 2 2 = + | c c 2 & 3 1 & - 2 | = - 7, C _ 2 3 = - | c c 2 & 1 1 & 4 | = - 7C _ 3 1 = + | c c 1 & 3 - 1 & 1 | = 4,C _ 3 2 = - | l l 2 & 3 1 & 1 | = 1,C _ 3 3 = + | c c 2 & 1 1 & - 1 | = - 3

伴随矩阵是余因子的矩阵的转置(例如, C₁2 去到(2,1)位置),从而

operatornameadj A = [ c c c - 2 & 1 4 & 4 3 & - 7 & 1 5 & - 7 & - 3 ]

我们可以直接计算 det A ,但下列计算对上面的结果提供了一个检验并计算出 det A 的方法:

(adj A) · A = [ r r r - 2 & 1 4 & 4 3 & - 7 & 1 5 & - 7 & - 3 ] [ r r r 2 & 1 & 3 1 & - 1 & 1 1 & 4 & - 2 ] = [ r r r 1 4 & 0 & 0 0 & 1 4 & 0 0 & 0 & 1 4 ] = 1 4 I

由于 (adjA)· A = 14I ,由定理8得 operatorname detA = 14 ,且

A ^ - 1 = frac 11 4 [ c c c - 2 & 1 4 & 4 3 & - 7 & 1 5 & - 7 & - 3 ] = [ c c c - 1 / 7 & 1 & 2 / 7 3 / 1 4 & - 1 / 2 & 1 / 1 4 5 / 1 4 & - 1 / 2 & - 3 / 1 4 ]

数值计算的注解 定理8主要用于理论上的计算.A⁻1 的公式使我们不用实际计算出 A⁻1 就可以推导出 A⁻1 的性质.如果确实需要求 A⁻1,除了特殊情形外, 2.2节中的算法给出了一个计算 A⁻1 的很好的方法.

克拉默法则也是一个理论工具,它可用来研究当 b 或 A 中某元素改变时,Ax = b 的解如何变化。(可能由于取得 b 或 A 中的元素时存在误差。)当 A 是一个具有复数元素的 3 × 3 矩阵时,因为涉及复数的四则计算, [A b] 的行化简可能很麻烦,而行列式的计算相对容易,故克拉默法则有时在手算时被选用。对一个很大的 n × n (实数或复数)矩阵,克拉默法则是无效的。仅计算一个行列式就大约与用行化简解 Ax = b 有相同的工作量。

在下一个应用中,我们证明在本章介绍中阐述的行列式的几何解释。尽管 ℝⁿ 中长度和距离的一般讨论直到第6章才给出,但我们在这里假设 ℝ² 和 ℝ³ 中通常的欧几里得长度、面积、体积已经是清楚的。

定理 9

若 A 是一个 2 × 2 矩阵,则由 A 的列确定的平行四边形的面积为 |det A| 。若 A 是一个 3 × 3 矩阵,则由 A 的列确定的平行六面体的体积为 |det A| 。

证明

若 A 为 2 阶对角矩阵,定理显然成立.

| det [ c c a & 0 0 & d ] | = | a d | = {text 矩阵的面积 }

见图 3-2. 若 A 不为对角情形,只需证 A = [a₁, a₂] 能变换成一个对角矩阵,变换时既不改变相应的平行四边形面积又不改变 |det A|.由3.2节我们知

图 3-2 面积=|ad|
道,当行列式的两列交换或一列的倍数加到另一列上时,行列式的绝对值不改变。同时容易看到,这样的运算足以将 A 变换成对角矩阵。由于列交换不会改变对应的平行四边形,所以只需证明下列应用于 ℝ² 和 ℝ³ 中的向量的简单几何现象就足够了。

设 a₁ 和 a₂ 为非零向量,则对任意数c,由 a₁ 和 a₂ 确定的平行四边形的面积等于由 a₁ 和 a₂+ca₁ 确定的平行四边形的面积.

为了证明这个结论,我们可以假设 a₂ 不是 a₁ 的倍数,否则这两个平行四边形将退化成面积为0.若 L 是通过0和 a₁ 的直线,则 a₂ + L 是通过 a₂ 且平行于 L 的直线, a₂ + ca₁ 在此直线上,见图 3-3.点 a₂ 和 a₂ + ca₁ 具有相同的垂直距离,因此图 3-3中的两个平行四边形具有相同的底边,即由0到 L 具有相同的垂直距离,因此图3-3中的两个平行四边形具有相同的底边,即由0到 a₁ 的线段,所以这两个平行四边形具有相同的面积,这就完成了 ℝ² 的情形的证明.

图 3-3 两个等面积的平行四边形

类似可证明 ℝ³ 的情形.定理对 3 × 3 对角阵显然成立,见图 3-4. 任意一个 3 × 3 矩阵 A 均可用不改变 |det A| 的列变换变换成对角矩阵 (考虑在 A^(T) 上做行变换), 所以只需证明这些变换不影响由 A 的列确定的平行六面体的体积就行了.

图 3-5中用阴影给出一个平行六面体,它具有两个倾斜的侧面,其体积等于在平面 operatornameSpan{±ba₁, ±ba₃} 上的底面的面积乘以 ±ba₂ 到 operatornameSpan{±ba₁, ±ba₃} 的垂直距离。任意向量 ±ba₂ + c±ba₁ 到 operatornameSpan{±ba₁, ±ba₃} 的距离与 ±ba₂ 到 operatornameSpan{±ba₁, ±ba₃} 的距离相同,这是因为 ±ba₂ + c±ba₁ 位于平行于 operatornameSpan{±ba₁, ±ba₃} 的平面 ±ba₂ + operatornameSpan{±ba₁, ±ba₃} 上,所以当 [±ba₁, ±ba₂, ±ba₃] 变为 [±ba₁, ±ba₂ + c±ba₁, ±ba₃] 时,平行六面体的体积不变。于是列的倍加变换不影响平行六面体的体积。由于列的交换也不影响体积,所以定理证

图3-4 体积 = |a b c|

图3-5 两个等体积的平行六面体

例 4

计算由点 (-2,-2) , (0,3) , (4,-1) 和 (6,4) 确定的平行四边形的面积.见图 3-6a.

图3-6 平移一个平行四边形不改变其面积

解

先将此平行四边形平移到使原点作为其某一顶点的情形.例如,将每个顶点坐标减去顶点 (-2, -2).这样,新的平行四边形面积与原平行四边形面积相同,其顶点为

(0, 0), (2, 5), (6, 1) text 和 (8, 6)

见图 3-6b. 此平行四边形由 A=2 & 6 5 & 1 的列所确定.由于 |det A|=-28 ,所以所求的平行四边形的面积为 28.

行列式可用于描述平面和 ℝ³ 中线性变换的一个重要的几何性质.若 T 是一个线性变换,S 是 T 的定义域内的一个集合,用 T(S) 表示 S 中点的像集.我们对 T(S) 的面积(体积)与原来的集合 S 的面积(体积)相对比有何变化这件事感兴趣.为了方便,当 S 是一个边界为平行四边形的区域时,我们就用 S 表示一个平行四边形.

定理 10

设 T: ℝ² → ℝ² 是由一个 2 × 2 矩阵 A 确定的线性变换,若 S 是 ℝ² 中一个平行四边形,则

{T (S) text 的面积 } = | det A | · {S text 的面积 } (5)

若 T 是一个由 3 × 3 矩阵 A 确定的线性变换,而 S 是 R³ 中的一个平行六面体,则

{T (S) text 的体积 } = | det A | · {S text 的体积 } (6)
证明

考虑 2 × 2 的情形,A = [a₁ a₂].一个顶点在 ℝ² 中原点的平行四边形由向量 ±bb₁ 和 ±bb₂ 确定,具有以下形式:

S = {s _ 1 ±b b _ 1 + s _ 2 ±b b _ 2: 0 ≤slant s _ 1 ≤slant 1, 0 ≤slant s _ 2 ≤slant 1 }

S 在 T 下的像由以下形式的点组成:

T (s _ 1 boldsymbol b _ 1 + s _ 2 boldsymbol b _ 2) = s _ 1 T (boldsymbol b _ 1) + s _ 2 T (boldsymbol b _ 2) = s _ 1 A boldsymbol b _ 1 + s _ 2 A boldsymbol b _ 2

其中 0 ≤slant s₁ ≤slant 1, 0 ≤slant s₂ ≤slant 1 ,从而 T(S) 是一个由矩阵 [Ab₁, Ab₂] 的列确定的平行四边形。这个矩阵可写成 AB ,其中 B = [b₁, b₂] 。由定理9和行列式的乘积定理,

r l {T (S) text 的面积 } & = | det A B | = | det A | · | det B | & = | det A | · {S text 的面积 } (7)

任意一个平行四边形都具有形式 p + S ,其中 ±bp 是一个向量, S 是一个和上面一样有一个顶点在原点的平行四边形。容易看出 T 将 ±bp + S 变换为 T(±bp) + T(S) (见习题26)。由于平移不改变一个集合的面积,所以

r l {T (±b p + S) text 的面积 } & = {T (±b p) + T (S) text 的面积 } & = {T (S) text 的面积 } text 平移 & = | det A | · {S text 的面积 } text 由(7) & = | det A | · {±b p + S text 的面积 } text 平移

这表明(5)式对 R² 中任意的平行四边形均成立。(6)式在 3 × 3 情形的证明是类似的.

当我们尝试把定理10推广到 ℝ² 和 ℝ³ 中不是由直线或平面所围的区域时,必须面对如何定义和计算这个区域的面积或体积的问题。这是一个在微积分中研究的问题,我们仅仅对 ℝ² 情形给出基本思想。若 R 是一个有有限面积的平面区域,则 R 的面积可由这些小正方形的面积之和充分地逼近。见图 3-7。R 的面积可由这些小正方形的面积之和充分地逼近。见图3-7。

图 3-7 由正方形近似一个平面区域.近似的程度随格子变细小而改进

若 T 是一个由 2 × 2 矩阵 A 确定的线性变换,则平面区域 R 在 T 下的像由 R 中小正方形的像来近似.定理10的证明说明每一个这样的像是一个平行四边形,它的面积是 |det A| 乘以这个正方形的面积.若 R' 是 R 中这些正方形的并集,则 T(R') 的面积是 |det A| 的面积,见图 3-8. R' 的面积.当这个过程无限进行下去时,就可以得到下面推广的定理10的证明.T(R') 的面积也逼近 T(R) 的面积.当这个过程无限进行下去时,就可以得到下面推广的定理10的证明.

图3-8 T(R) 由平行四边形的并集近似

定理10 的结论对 R² 中任意具有有限面积的区域或 R³ 中具有有限体积的区域均成立.

例 5

若 a, b 是正数,求由方程为 (x₁²)/(a²) + (x₂²)/(b²) = 1 的椭圆为边界的区域 E 的面积.

解

我们断言 E 是单位圆盘 D 在线性变换 T 下的像.这里 T 由矩阵 A = a & 0 0 & b 确定,这是因为若 ±bu = u₁ u₂ , ±bx = x₁ x₂ ,且 ±bx = A±bu ,则

u _ 1 = frac x _ 1a, u _ 2 = frac x _ 2b

从而得 u 在此单位圆盘内,即满足 u₁² + u₂² ≤slant 1 ,当且仅当 x 在 E 内,即满足 (x₁/a)² + (x₂/b)² ≤slant 1 。由定理10 的推广,

r l & {text 椭圆的面积 } = {T (D) text 的面积 } & = | det A | · {D text 的面积 } & = a · b · π · (1) ^ 2 = π a b
练习题

设 S 为由向量 ±bb₁ = 1 3 和 ±bb₂ = 5 1 确定的平行四边形,令 A = 1 & -0.1 0 & 2.计算 S 在映射 ±bx → A±bx 下的像的面积.

习题3.3

利用克拉默法则求习题1~6中方程组的解.

  1. 5x₁ + 7x₂ = 3 2. 4x₁ + x₂ = 6 2x₁ + 4x₂ = 1 3x₁ + 2x₂ = 7 3. 3x₁ - 2x₂ = 3 4. -5x₁ + 2x₂ = 9 -4x₁ + 6x₂ = -5 3x₁ - x₂ = -4 5. x₁ + x₂ = 3 6. x₁ + 3x₂ + x₃ = 4 -3x₁ + 2x₃ = 0 -x₁ + 2x₃ = 2 x₂ - 2x₃ = 2 3x₁ + x₂ = 2

在习题 7~10 中,确定参数 s 的值,使方程组有唯一解,并求出该解. 7. 6sx₁ + 4x₂ = 5 8. 3sx₁ + 5x₂ = 3 9x₁ + 2sx₂ = -2 12x₁ + 5sx₂ = 2 9. sx₁ + 2sx₂ = -1 10. sx₁ - 2x₂ = 1 3x₁ + 6sx₂ = 4 4sx₁ + 4sx₂ = 2

在习题 11~16 中,求给出矩阵的伴随矩阵,然后利用定理8 给出矩阵的逆矩阵. 11. 0&-2&-15&0&0-1&1&1 12. 1&1&3-2&2&10&1&1 13. 3&5&41&0&12&1&1 14. 1&-1&20&2&12&0&4 15. 5&0&0-1&1&0-2&3&-1 16. 1&2&40&-3&10&0&-2

  1. 证明:若 A 是 2 × 2 矩阵,则对 A⁻1 ,定理8与

2.2节中定理4给出相同的公式.

  1. 假设 A 中所有元素均为整数,且 operatornamedet A = 1 。解释为什么 A⁻1 中所有元素也是整数。在习题19~22中,分别求出给定顶点的平行四边形的面积。

  2. (0,0),(5,2),(6,4),(11,6)

  3. (0,0),(-2,4),(4,-5),(2,-1)

  4. (-2,0),(0,3),(1,3),(-1,0)

  5. (0, -2) , (5, -2) , (-3, 1) , (2, 1)

  6. 求一个顶点在原点且相邻顶点在 (1,0,-3) , (1,2,4) , (5,1,0) 的平行六面体的体积.

  7. 求一个顶点在原点且相邻顶点在 (1,3,0) , (-2,0,2) , (-1,3,-1) 的平行六面体的体积.

  8. 利用体积的概念解释为什么 3 × 3 矩阵 A 的行列式为零,当且仅当 A 不可逆。要求不利用 3.2 节中的定理4。(提示:考虑 A 的列。)

  9. 令 T: ℝ^m → ℝⁿ 为线性变换, ±bp 为 ℝ^m 中一个向量, S 为 ℝ^m 中集合,证明: ±bp + S 在 T 下的像等于 ℝⁿ 中平移的集合 T(±bp) + T(S) .

  10. 令 S 是由向量 ±bb₁ = -2 3 和 ±bb₂ = -2 5 确定的平行四边形, A = 6 & -3 -3 & 2 ,求 S 在映射 ±bx → Ax

下像的面积.

  1. 对 b₁=4 -7, b₂=0 1, A=5&2 1&1 ,重复 27 题.

  2. 求 R² 中顶点为 0, v₁ , v₂ 的三角形的面积公式.

  3. 令 R 是顶点为 (x₁, y₁), (x₂, y₂), (x₃, y₃) 的三角形.证明: {三角形的面积} = 1/2 | det x₁ & y₁ & 1 x₂ & y₂ & 1 x₃ & y₃ & 1 | . (提示:通过减去一个顶点向量,将 R 平移到原点,再利用习题 29.)

  4. 令 T: ℝ³ → ℝ³ 是由矩阵 A = a & 0 & 0 0 & b & 0 0 & 0 & c 确定的线性变换,其中 a, b, c 为正数。令 S 为单位球体,其边界曲面方程为 x₁² + x₂² + x₃² = 1 。 a. 证明: T(S) 的边界为椭球面,方程为 (x₁²)/(a²) + (x₂²)/(b²) + (x₃²)/(c²) = 1 。 b. 利用单位球体体积为 4/3π 的事实确定边界为(a)中椭球面的区域的体积。

  5. 令 S 是 R³ 中具有顶点 0, e₁, e₂, e₃ 的四面体, S' 为具有顶点 0, v₁, v₂, v₃ 的四面体,见下面的图形.

a.描述一个将 S 映射到 S' 上的线性变换.

b.利用 {S' 的体积 } = 1/3{底面积} · {高} 这个事实给出一个四面体 S' 的体积公式.

  1. [M]任取一个 4 × 4 矩阵,检验定理8中的逆矩阵公式.利用你的矩阵程序计算所有的 3 × 3 子矩阵的余因子,构造伴随矩阵,令 B = (operatornameadj A) / (det A).计算 B - operatornameinv(A),其中 operatornameinv(A) 是由矩阵程序计算出的 A 的逆.使用浮点运算保留最大可能的小数位.报告你的结论.

  2. [M]任取一个 4 × 4 矩阵 A 和一个 4 × 1 向量 ±bb,检验克拉默法则.求 Ax = b 的解并与 A⁻1b 中元素相比较.利用克拉默法则,对你的计算机程序,写出输出 ±bx 中第二个元素的指令.

  3. [M]在 MATLAB 中,对一个任意取的 30 × 30 的矩阵 A ,使用 flops 命令统计出计算 A⁻1 所需要的浮点运算的次数。将这个数与计算 (adj A)/(det A) 所需的浮点运算次数相比较。

练习题答案

S 的面积为 |det 1 & 5 3 & 1 | = 14, det A = 2 。由定理10, S 在映射 x → Ax 下的像的面积为 |det A| · {S 的面积 } = 2 · 14 = 28 .

共 0 题