矩阵的代数运算是线性代数的核心计算基础。本节系统讨论矩阵的五大基本代数运算:加法与减法、数乘运算、矩阵乘法、转置运算(及对称矩阵与反对称矩阵)以及方阵的行列式,重点阐释矩阵乘法不满足交换律与消去律的本质特征及其在网络连通分析中的几何应用。
一、矩阵的加法
Section titled “一、矩阵的加法”设有两个同型的 m × n 矩阵 A = (aᵢj)ₘ × n 和 B = (bᵢj)ₘ × n,由 A 与 B 的对应元素相加所得到的一个 m × n 矩阵,称为 A 与 B 的和,记为 A + B,即
A + B = (aᵢj + bᵢj)ₘ × n.只有当两个矩阵 A 与 B 是同型矩阵时,它们才可以相加。不同型的矩阵之间没有加法运算。
矩阵加法的运算律
Section titled “矩阵加法的运算律”设 A, B, C 都是 m × n 矩阵,O = (0)ₘ × n 是同型零矩阵,则有:
- 交换律:A + B = B + A;
- 结合律:(A + B) + C = A + (B + C);
- 零元性质:A + O = O + A = A;
- 消去律:A + C = B + C iff A = B。
负矩阵与减法
Section titled “负矩阵与减法”设 A = (aᵢj)ₘ × n,称矩阵
-a₁1 & -a₁2 & … & -a₁n -a₂1 & -a₂2 & … & -a₂n vdots & vdots & ddots & vdots -aₘ1 & -aₘ2 & … & -aₘn为 A 的负矩阵,记为 -A。显然有
A + (-A) = (-A) + A = O.由此定义矩阵的减法为:
A - B = A + (-B) = (aᵢj - bᵢj)ₘ × n.二、数乘运算
Section titled “二、数乘运算”对于任意一个矩阵 A = (aᵢj)ₘ × n 和任意一个数 k,规定 k 与 A 的乘积为
kA = ka₁1 & ka₁2 & … & ka₁n ka₂1 & ka₂2 & … & ka₂n vdots & vdots & ddots & vdots kaₘ1 & kaₘ2 & … & kaₘn .数 k 与矩阵 A 的乘积是矩阵 A 中的所有元素都要乘 k。这与行列式提公因式仅针对单行(列)有着本质区别。
根据数乘矩阵运算的定义可知,数量矩阵 aEₙ 就是数 a 与单位矩阵 Eₙ 的乘积。而且对任意矩阵 A 有 -A = (-1)A = (-aᵢj)ₘ × n。
矩阵的加法与数乘运算合称为矩阵的线性运算,满足下面的运算律(设 A, B 为 m × n 矩阵,k, l 为常数):
- 结合律:(kl)A = k(lA) = l(kA);
- 分配律:k(A + B) = kA + kB,(k + l)A = kA + lA。
已知
A = -1 & 2 & 3 & 1 0 & 2 & -1 & 3 4 & 2 & 0 & 5 , B = 1 & 2 & -1 & 0 4 & -3 & 1 & 1 1 & 0 & 2 & 5 ,求 2A - 3B。
解
根据数乘与减法运算法则:
2A - 3B &= 2 -1 & 2 & 3 & 1 0 & 2 & -1 & 3 4 & 2 & 0 & 5 - 3 1 & 2 & -1 & 0 4 & -3 & 1 & 1 1 & 0 & 2 & 5 &= -2 - 3 & 4 - 6 & 6 + 3 & 2 - 0 0 - 12 & 4 + 9 & -2 - 3 & 6 - 3 8 - 3 & 4 - 0 & 0 - 6 & 10 - 15 = -5 & -2 & 9 & 2 -12 & 13 & -5 & 3 5 & 4 & -6 & -5 .已知
A = 3 & 0 & -1 & 2 2 & 8 & 3 & 1 , B = 5 & 6 & 3 & 2 2 & 4 & 7 & -1 ,且 A + 2X = B,求 X。
解
由方程可得 2X = B - A,因此
X = 1/2(B - A) = 1/2 2 & 6 & 4 & 0 0 & -4 & 4 & -2 = 1 & 3 & 2 & 0 0 & -2 & 2 & -1 .三、矩阵的乘法
Section titled “三、矩阵的乘法”矩阵的乘法是矩阵最具特色的运算。在给出它的定义之前,先看下面的变量代换模型。
设变量 y₁, y₂, y₃ 能用变量 x₁, x₂, x₃ 线性地表示为
y₁ = b₁1x₁ + b₁2x₂ + b₁3x₃, y₂ = b₂1x₁ + b₂2x₂ + b₂3x₃, y₃ = b₃1x₁ + b₃2x₂ + b₃3x₃, (2.4)其中 bᵢj 都是常数,式 (2.4) 是由变量 x₁, x₂, x₃ 到变量 y₁, y₂, y₃ 的线性变换。
如果另有变量 y₁, y₂, y₃ 到变量 z₁, z₂, z₃ 的线性变换
z₁ = a₁1y₁ + a₁2y₂ + a₁3y₃, z₂ = a₂1y₁ + a₂2y₂ + a₂3y₃, z₃ = a₃1y₁ + a₃2y₂ + a₃3y₃, (2.5)将 (2.4) 代入 (2.5),整理后可得由变量 x₁, x₂, x₃ 到变量 z₁, z₂, z₃ 的复合线性变换:
z₁ = (a₁1b₁1 + a₁2b₂1 + a₁3b₃1)x₁ + (a₁1b₁2 + a₁2b₂2 + a₁3b₃2)x₂ + (a₁1b₁3 + a₁2b₂3 + a₁3b₃3)x₃, z₂ = (a₂1b₁1 + a₂2b₂1 + a₂3b₃1)x₁ + (a₂1b₁2 + a₂2b₂2 + a₂3b₃2)x₂ + (a₂1b₁3 + a₂2b₂3 + a₂3b₃3)x₃, z₃ = (a₃1b₁1 + a₃2b₂1 + a₃3b₃1)x₁ + (a₃1b₁2 + a₃2b₂2 + a₃3b₃2)x₂ + (a₃1b₁3 + a₃2b₂3 + a₃3b₃3)x₃. (2.6)如果将线性变换 (2.4), (2.5), (2.6) 的系数矩阵分别记为 B = (bᵢj), A = (aᵢj) 和 C = (cᵢj),则矩阵元素之间的关系为:
cᵢj = aᵢ1b₁j + aᵢ2b₂j + aᵢ3b₃j = ∑ₖ=1³ aᵢkbₖj (i, j = 1, 2, 3).即矩阵 C 的第 i 行第 j 列位置的元素等于矩阵 A 的第 i 行元素与矩阵 B 的第 j 列对应元素乘积的和。
由此自然引出一般矩阵相乘的定义:
设矩阵 A = (aᵢj)ₘ × s, B = (bᵢj)ₛ × n。令 C = (cᵢj)ₘ × n 是由
cᵢj = aᵢ1b₁j + aᵢ2b₂j + … + aᵢsbₛj = ∑ₖ=1^s aᵢkbₖj (i = 1, 2, …, m; j = 1, 2, …, n)构成的 m 行 n 列矩阵,称矩阵 C 为矩阵 A 与矩阵 B 的乘积,记为 C = AB。
两个矩阵 A 和 B 可以相乘,当且仅当 A 的列数与 B 的行数相等。当 C = AB 时,C 的行数等于 A 的行数,C 的列数等于 B 的列数。
设矩阵
A = 1 & 0 & -1 2 & 1 & 0 3 & 2 & -1 , B = 1 & 0 3 & 1 0 & 2 ,求 AB。
解
根据乘法定义:
AB &= 1 & 0 & -1 2 & 1 & 0 3 & 2 & -1 1 & 0 3 & 1 0 & 2 &= 1 × 1 + 0 × 3 + (-1) × 0 & 1 × 0 + 0 × 1 + (-1) × 2 2 × 1 + 1 × 3 + 0 × 0 & 2 × 0 + 1 × 1 + 0 × 2 3 × 1 + 2 × 3 + (-1) × 0 & 3 × 0 + 2 × 1 + (-1) × 2 = 1 & -2 5 & 1 9 & 0 .这里矩阵 A 是 3 × 3 矩阵,而矩阵 B 是 3 × 2 矩阵,由于 B 的列数与 A 的行数不相等,所以 BA 没有意义。
在矩阵的乘法中,必须特别注意矩阵相乘的顺序!若 A 是 m × n 矩阵,B 是 n × m 矩阵,则 AB 与 BA 均有意义,但 AB 是 m 阶方阵,而 BA 是 n 阶方阵。当 m ≠ n 时,AB ≠ BA。即使 m = n,即 A, B 是同阶方阵,AB 与 BA 仍然通常不相等。
设矩阵
A = 1 & 0 1 & 0 , B = 0 & 0 1 & 1 ,求 AB 和 BA。
解
直接计算可得:
AB = 1 & 0 1 & 0 0 & 0 1 & 1 = 0 & 0 0 & 0 ,BA = 0 & 0 1 & 1 1 & 0 1 & 0 = 0 & 0 2 & 0 .此例明确揭示了矩阵乘法的两个根本事实:
- 矩阵乘法一般不满足交换律,即一般有 AB ≠ BA;
- 当 AB = O 时,一般不能推出 A = O 或 B = O。这说明矩阵乘法不满足消去律。若 ab = 0,在实数范围内必有 a = 0 或 b = 0;但在矩阵世界中,两个非零矩阵的乘积完全可以是零矩阵。
矩阵乘法的运算律
Section titled “矩阵乘法的运算律”当两个方阵满足 AB = BA 时,称 A 与 B 是可交换的(此时 A 与 B 必为同阶方阵)。
矩阵乘法满足下列运算律(假设运算都是有意义的):
- 结合律:(AB)C = A(BC);
- 分配律:(A + B)C = AC + BC,A(B + C) = AB + AC;
- 数乘结合律:k(AB) = (kA)B = A(kB)(k 为任意实数);
- 单位阵性质:Eₘ Aₘ × n = Aₘ × n Eₙ = Aₘ × n。
关于矩阵乘法应高度警惕:当 A = B 时,在可乘条件下必可推出 AC = BC 和 CA = CB;但反之由 AC = BC 且 C ≠ O,未必能推出 A = B。
利用矩阵乘法,式 (2.2) 中的线性变换可紧凑表示为:
y = Ax,其中 x = (x₁, x₂, …, xₙ)^(T), y = (y₁, y₂, …, yₙ)^(T), A = (aᵢj)ₙ × n,A 称为该线性变换的系数矩阵。
方阵的幂与多项式
Section titled “方阵的幂与多项式”设 A 为 n 阶方阵,由于矩阵乘法满足结合律,所以 underbraceAA…Aₖ 个 可以不加括号而有完全确定的意义。规定:
A⁰ = E, A¹ = A, A² = AA, …, A^k = underbraceAA…Aₖ 个.方阵的幂满足指数运算法则:
A^k A^l = A^(k+l), (A^k)^l = A^(kl) (k, l ∈ ℕ⁺).由于矩阵乘法一般不满足交换律,对于两个同阶方阵 A 与 B,一般有 (AB)^k ≠ A^k B^k。只有当 AB = BA 时该等式才成立。
设 x 的 m 次多项式 f(x) = a₀ x^m + a₁ x^(m-1) + … + aₘ-1x + aₘ。定义 A 的 m 次多项式为:
f(A) = a₀ A^m + a₁ A^(m-1) + … + aₘ-1A + aₘ E.设矩阵 A = 1 & 0 2 & 1,求所有与 A 可交换的矩阵。
解
因为与 A 可交换的矩阵必为二阶矩阵,设 X = x₁1 & x₁2 x₂1 & x₂2。由 AX = XA:
AX = 1 & 0 2 & 1 x₁1 & x₁2 x₂1 & x₂2 = x₁1 & x₁2 2x₁1 + x₂1 & 2x₁2 + x₂2 ,XA = x₁1 & x₁2 x₂1 & x₂2 1 & 0 2 & 1 = x₁1 + 2x₁2 & x₁2 x₂1 + 2x₂2 & x₂2 .对比对应元素:x₁1 = x₁1 + 2x₁2 implies x₁2 = 0;2x₁1 + x₂1 = x₂1 + 2x₂2 implies x₁1 = x₂2。
而 x₁1, x₂1 可取任意值。因此所有与 A 可交换的矩阵具有形式:
X = x₁1 & 0 x₂1 & x₁1 .用数学归纳法证明以下矩阵等式:
- 1 & 1 0 & 1 ⁿ = 1 & n 0 & 1;
- 1 & 1 1 & 1 ⁿ = 2ⁿ-1 1 & 1 1 & 1。
证明
(1) 当 n = 1 时显然成立。假设当 n = k 时等式成立,则当 n = k+1 时:
1 & 1 0 & 1 ^(k+1) = 1 & 1 0 & 1 ^k 1 & 1 0 & 1 = 1 & k 0 & 1 1 & 1 0 & 1 = 1 & k+1 0 & 1 .由数学归纳法,对任意正整数 n 等式成立。
(2) 当 n = 1 时等式两端均为 1 & 1 1 & 1,成立。假设当 n = k 时成立,则当 n = k+1 时:
1 & 1 1 & 1 ^(k+1) = 1 & 1 1 & 1 ^k 1 & 1 1 & 1 = 2^(k-1) 1 & 1 1 & 1 1 & 1 1 & 1 = 2^(k-1) 2 & 2 2 & 2 = 2^k 1 & 1 1 & 1 .由数学归纳法,结论对任意正整数 n 成立。
设 A = 1 & -2 -2 & 4 1 & -2 , B = 1 & 2 & 3 3 & 2 & 1。令 C = AB, D = BA。求 C, D, C², C³ 和 Cⁿ(n ≥ 2)。
解
先计算 D:
D = BA = 1 & 2 & 3 3 & 2 & 1 1 & -2 -2 & 4 1 & -2 = 0 & 0 0 & 0 = O.再计算 C = AB:
C = 1 & -2 -2 & 4 1 & -2 1 & 2 & 3 3 & 2 & 1 = -5 & -2 & 1 10 & 4 & -2 -5 & -2 & 1 .利用结合律与 BA = O:
C² = (AB)² = A(BA)B = AOB = O,C³ = (AB)³ = A(BA)(BA)B = O.一般地,对任意 n ≥ 2:
Cⁿ = (AB)ⁿ = Aunderbrace(BA)…(BA)ₙ-1 个B = O.若 f(x) = x² - 2x + 2,A = 2 & 0 & 1 -3 & 1 & 2 1 & -1 & 0,求 f(A)。
解
计算 A²:
A² = 5 & -1 & 2 -7 & -1 & -1 5 & -1 & -1 .代入求得:
f(A) = 5 & -1 & 2 -7 & -1 & -1 5 & -1 & -1 - 4 & 0 & 2 -6 & 2 & 4 2 & -2 & 0 + 2 & 0 & 0 0 & 2 & 0 0 & 0 & 2 = 3 & -1 & 0 -1 & -1 & -5 3 & 1 & 1 .矩阵乘法在网络连通分析中的应用
Section titled “矩阵乘法在网络连通分析中的应用”回顾例2 中的通信网络邻接矩阵:
A = 0 & 1 & 1 & 1 1 & 0 & 1 & 0 1 & 1 & 0 & 1 1 & 0 & 1 & 0 , A² = 3 & 1 & 2 & 1 1 & 2 & 1 & 2 2 & 1 & 3 & 1 1 & 2 & 1 & 2 .由矩阵乘法定义可知,A² 的第 i 行第 j 列元素为:
[A²]ᵢj = aᵢ1a₁j + aᵢ2a₂j + aᵢ3a₃j + aᵢ4a₄j.由于 aᵢtaₜj = 1 iff aᵢt = aₜj = 1(即点 i 与点 t、点 t 与点 j 均连通),因此 [A²]ᵢj 正好表示从通信点 i 经由一个中间中转点 t 传输信息至通信点 j 的通路数目!例如 [A²]₂4 = 2,表示第 2 点经由一个中间点到第 4 点有 2 条通道(即 2 → 1 → 4 和 2 → 3 → 4)。
进而矩阵
M₂ = A + A² = 3 & 2 & 3 & 2 2 & 2 & 2 & 2 3 & 2 & 3 & 2 2 & 2 & 2 & 2的元素 [M₂]ᵢj 表示从点 i 出发,直接或者经由一个中转点传达至点 j 的总通路数目。
四、转置矩阵、对称矩阵和反对称矩阵
Section titled “四、转置矩阵、对称矩阵和反对称矩阵”把矩阵 A 的所有行换成相应的列所得到的矩阵,称为矩阵 A 的转置矩阵,记为 A^(T) 或 A'。若
A = a₁1 & a₁2 & … & a₁n a₂1 & a₂2 & … & a₂n vdots & vdots & ddots & vdots aₘ1 & aₘ2 & … & aₘn ₘ × n, 则 A^(T) = a₁1 & a₂1 & … & aₘ1 a₁2 & a₂2 & … & aₘ2 vdots & vdots & ddots & vdots a₁n & a₂n & … & aₘn ₙ × m.显然 A 与 A^(T) 互为转置矩阵。特别地,n 维行向量的转置为 n 维列向量。
转置矩阵的运算律
Section titled “转置矩阵的运算律”设运算均有意义,k 为实数:
- (A^(T))^(T) = A;
- (A + B)^(T) = A^(T) + B^(T);
- (kA)^(T) = kA^(T);
- 反序律:(AB)^(T) = B^(T)A^(T)。推广到多个矩阵:(A₁ A₂ … Aₖ)^(T) = Aₖ^(T) … A₂^(T) A₁^(T)。
反序律证明
设 A = (aᵢj)ₘ × s, B = (bᵢj)ₛ × n。显然 (AB)^(T) 和 B^(T)A^(T) 都是 n × m 矩阵。
(AB)^(T) 的第 (i, j) 元素为:
[(AB)^(T)]ᵢj = [AB]ⱼi = ∑ₜ=1^s aⱼtbₜi = ∑ₜ=1^s bₜiaⱼt.而 B^(T) 的第 i 行为 (b₁i, b₂i, …, bₛi),A^(T) 的第 j 列为 (aⱼ1, aⱼ2, …, aⱼs)^(T)。故 B^(T)A^(T) 的第 (i, j) 元素也是 ∑ₜ=1^s bₜiaⱼt。因此 (AB)^(T) = B^(T)A^(T)。
对称矩阵与反对称矩阵
Section titled “对称矩阵与反对称矩阵”设 A = (aᵢj) 为 n 阶方阵。
- 若 A 满足 A^(T) = A,即 aᵢj = aⱼi(i, j = 1, 2, …, n),则称 A 为对称矩阵;
- 若 A 满足 A^(T) = -A,即 aᵢj = -aⱼi(i, j = 1, 2, …, n),则称 A 为反对称矩阵。
对于反对称矩阵,当 i = j 时有 aᵢi = -aᵢi implies aᵢi = 0。即反对称矩阵的主对角元必全为零!
证明任意一个方阵 A 都可以唯一地表示为一个对称矩阵与一个反对称矩阵之和。
证明
设 A = X + Y,其中 X^(T) = X,Y^(T) = -Y。两端取转置:
A^(T) = X^(T) + Y^(T) = X - Y.联立方程:
X + Y = A, X - Y = A^(T).解得唯一解:
X = 1/2(A + A^(T)), Y = 1/2(A - A^(T)).易验:X^(T) = 1/2(A^(T) + A) = X,Y^(T) = 1/2(A^(T) - A) = -Y。
五、方阵的行列式
Section titled “五、方阵的行列式”由 n 阶方阵 A 的元素按照原来的顺序构成的行列式称为方阵 A 的行列式,记作 |A| 或 det(A)。
当且仅当 A 为方阵时,才可以取行列式。对于不是方阵的长方矩阵,绝对不存在行列式。
方阵行列式的运算性质
Section titled “方阵行列式的运算性质”设 A, B 为 n 阶方阵,k 为数:
- |A^(T)| = |A|;
- |kA| = kⁿ |A|;
- 行列式乘法定理:|AB| = |A| · |B|。
对于同阶方阵 A, B,虽然一般有 AB ≠ BA,但由乘法定理可知它们的行列式恒相等:
|AB| = |A| · |B| = |B| · |A| = |BA|.设 A = 1 & 3 2 & -2 , B = 2 & 5 3 & 4。
计算验证 |AB| = |BA| = |A| · |B|。
解
计算乘积:
AB = 1 & 3 2 & -2 2 & 5 3 & 4 = 11 & 17 -2 & 2 , |AB| = 11 × 2 - 17 × (-2) = 56.BA = 2 & 5 3 & 4 1 & 3 2 & -2 = 12 & -4 11 & 1 , |BA| = 12 × 1 - (-4) × 11 = 56.而 |A| = -2 - 6 = -8,|B| = 8 - 15 = -7。因此 |A| · |B| = (-8) × (-7) = 56。等式完全成立。
设 A 是 n 阶方阵,满足 AA^(T) = E,且 |A| = -1。求 |A + E|。
解
由于 AA^(T) = E:
|A + E| &= |A + AA^(T)| = |A(E + A^(T))| &= |A| · |E + A^(T)| &= -|(E + A)^(T)| &= -|E + A| = -|A + E|.因此 2|A + E| = 0 implies |A + E| = 0。