《线性代数及其应用》

引言

线性代数不仅是大学数学的重要组成部分,更是理解现代人工智能系统的核心语言。在向量空间中建模,在矩阵变换中理解结构,在特征分解中揭示本质,线代的每一个概念都与机器学习、深度学习和强化学习息息相关。

本系列笔记旨在对线性代数的基础概念、常用定理、典型方法与几何直觉进行系统梳理与逐章讲解,并融合我在人工智能学习中的理解与应用联系,为之后深入强化学习等课程打下坚实的数学基础。

在这里我也会给出一些优秀的学习资料参考,希望对大家有所帮助:

第一章 线性代数中的线性方程组

1.1 线性方程组

线性方程组是形如下式的方程:

a1x1+a2x2+...+anxn=ba_1x_1+a_2x_2+...+a_nx_n=b

其中 bb 与系数 aia_i 为实数或复数,通常是已知的,未知数 xix_i 为变量。

方程组所有可能的解的集合称为线性方程组的解集

若两个线性方程组有相同的解集,则这两个线性方程组称为等价的。

线性方程组的解有如下三种情况:

  1. 无解
  2. 有唯一解
  3. 有无穷多个解

我们称一个线性方程组是相容的,如果它至少有一个解。否则,称其为不相容的

这里还涉及到概念系数矩阵A增广矩阵C

A=[a11a12...a1na21a22...a2n............am1am2...amn]A=\begin{bmatrix} a_{11} & a_{12} & ... & a_{1n} \\ a_{21} & a_{22} & ... & a_{2n} \\ ... & ... & ... & ... \\ a_{m1} & a_{m2} & ... & a_{mn} \end{bmatrix}

C=[AB]=[a11a12...a1nb1a21a22...a2nb2...............am1am2...amnbm]C=[A\quad B]=\begin{bmatrix} a_{11} & a_{12} & ... & a_{1n} & b_1 \\ a_{21} & a_{22} & ... & a_{2n} & b_2 \\ ... & ... & ... & ... & ... \\ a_{m1} & a_{m2} & ... & a_{mn} & b_m \end{bmatrix}

常见的三种初等行变换:

  1. (倍加变换)把某一行换成它本身与另一行的倍数的和;
  2. (对调变换)交换矩阵的两行;
  3. (倍乘变换)用一个非零常数乘某一行。

如果其中一个矩阵可以经过一系列初等行变换称为另外一个矩阵,我们称这两个矩阵是行等价的,而且行变换是可逆的

如果两个线性方程组的增广矩阵行等价的,那么这两个线性方程组有相同的解集


1.2 行化简与阶梯形矩阵

1.2.1 阶梯形与行简化阶梯形概念

我们先引入如下两个重要概念:

  1. 非零行或列:指矩阵中至少包含一个非零元素的行或列;
  2. 非零行的先导元素:该行中最左边的非零元素。

接下来我们引入阶梯形,我们称一个矩阵是阶梯形矩阵,满足以下三条性质:

  1. 每一个非零行都在零行之上;
  2. 非零行的先导元素所在列的下面元素全为零;
  3. 某一行行的先导元素所在的列位于前一行先导元素所在列的右侧。

形如下式:

[2321014800052]\left[ \begin{array}{rrrr} 2 & -3 & 2 & 1 \\ 0 & 1 & -4 & 8 \\ 0 & 0 & 0 & \dfrac{5}{2} \end{array} \right]

再引入行简化阶梯形,在满足阶梯形的基础上,满足以下两条性质:

  1. 每一非零行的先导元素为1;
  2. 每一先导元素1是该元素所在列的唯一非零元素。

形如下式:

[10029010160013]\left[ \begin{array}{ccc|c} 1 & 0 & 0 & 29 \\ 0 & 1 & 0 & 16 \\ 0 & 0 & 1 & 3 \end{array} \right]

这里我们引入一条重要定理:每个矩阵行等价于唯一的简化阶梯形矩阵。

任何非零矩阵都可以进行行简化变为阶梯形矩阵,但用不同的方法可以得到不同的阶梯形矩阵;然而,一个矩阵只能化作唯一的行简化阶梯形矩阵

1.2.2 主元位置

这里我们定义如下:

  1. 主元:非零行的先导元素1;
  2. 主元位置:非零行的先导元素1所在位置;
  3. 主元列AA 的含有主元位置的列。

1.2.3 行化简算法

算法前四步骤产生一个阶梯形矩阵,第五步产生简化阶梯形矩阵。这里我们引入一个实例说明:

[03664537858939129615]\left[ \begin{array}{rrrrrr} 0 & 3 & -6 & 6 & 4 & -5 \\ 3 & -7 & 8 & -5 & 8 & 9 \\ 3 & -9 & 12 & -9 & 6 & 15 \end{array} \right]

第一步:由最左的非零行开始,这是一个主元列,主元位置在该列顶端。

[03664537858939129615]\left[ \begin{array}{rrrrrr} 0 & 3 & -6 & 6 & 4 & -5 \\ 3 & -7 & 8 & -5 & 8 & 9 \\ 3 & -9 & 12 & -9 & 6 & 15 \end{array} \right]

其中第 1 列为当前主元列,因为在该列中第 2 行或第 3 行的元素为非零,适合进行行交换以选取主元。

第二步:在主元列中选取一个非零元素作为主元,若有必要的话,对换两行使这个元素移动到主元位置上。

[39129615378589036645]\left[ \begin{array}{rrrrrr} \mathbf{3} & -9 & 12 & -9 & 6 & 15 \\ 3 & -7 & 8 & -5 & 8 & 9 \\ 0 & 3 & -6 & 6 & 4 & -5 \end{array} \right]

其中第 1 行第 1 列的元素 3\mathbf{3} 为主元(pivot)。

第三步:用倍加行变换将主元下面的元素变为0。

[39129615024426036645]\left[ \begin{array}{rrrrrr} \mathbf{3} & -9 & 12 & -9 & 6 & 15 \\ 0 & 2 & -4 & 4 & 2 & -6 \\ 0 & 3 & -6 & 6 & 4 & -5 \end{array} \right]

第四步:暂时不管包含主元位置的行以及它上面的各行,对剩下的子矩阵使用上述的三个步骤直到没有非零行需要处理为止。

[39129615024426000014]\left[ \begin{array}{rrrrrr} \mathbf{3} & -9 & 12 & -9 & 6 & 15 \\ 0 & 2 & -4 & 4 & 2 & -6 \\ 0 & 0 & 0 & 0 & \mathbf{1} & 4 \end{array} \right]

第五步:由最右边的主元开始,把每个主元上方的各元素变为0,若某个主元不是1,用倍乘变换将其变为1。

[1023024012207000014]\left[ \begin{array}{rrrrrr} 1 & 0 & -2 & 3 & 0 & -24 \\ 0 & 1 & -2 & 2 & 0 & -7 \\ 0 & 0 & 0 & 0 & 1 & 4 \end{array} \right]

1.2.4 线性方程组的解

行简化算法应用于方程组的增广矩阵的时候,可以得到线性方程组解集的一种显示表示法

我们先引入基本变量自由变量的概念:

  1. 基本变量:对应于主元列的变量;
  2. 自由变量:对应于非主元列的变量。

这里我们如下示例:

[105101140000]\left[ \begin{array}{rrrr} 1 & 0 & -5 & 1 \\ 0 & 1 & 1 & 4 \\ 0 & 0 & 0 & 0 \end{array} \right]

这里我们转换为方程组来表示为:

x15x3=1x2+x3=4\begin{aligned} x_1-5x_3&=1 \\ x_2+x_3&=4 \end{aligned}

这里 x1x_1x2x_2基本变量x3x_3自由变量

只要一个线性方程组是相容的,其解集就可以显示表示,可以如下形式表示通解

{x1=1+5x3x2=4x3x3 是自由变量\left\{ \begin{aligned} x_1 &= 1 + 5x_3 \\ x_2 &= 4 - x_3 \\ x_3 &\text{ 是自由变量} \end{aligned} \right.

1.2.5 解集的参数表示

如式子 (1.14) 形式就是解集的参数表示,其中自由变量作为参数。解方程组就是要求出解集的这种参数表示或确定它无解。

当方程组不相容的时候,无论是否有自由变量,解集都是空集。

1.2.6 存在性与唯一性问题

线性方程组相容的充分必要条件是增广矩阵的最右列不是主元列,也就是说,增广矩阵的阶梯形没有形如:

[00b]b0[0 \quad \cdots \quad 0 \quad b],b\neq 0

的行。若线性方程组是相容的,则它的解集有如下两种情况:

  1. 当没有自由变量的时候,有唯一解;
  2. 当有自由变量的时候,有无穷多解。

以下是求解线性方程组的步骤:

  1. 写出方程组的增广矩阵;
  2. 用行化简法将增广矩阵化为阶梯形,确定方程组是否相容,如果没有解就停止;否则进行下一步;
  3. 继续行化简算法得到它的简化阶梯形;
  4. 写出由第 3 步所得矩阵对应的方程组;
  5. 把第 4 步所得的每个非零方程改写为用任意自由变量表示其基本变量的形式。

1.3 向量方程

1.3.1 向量的基本概念

在线性代数中,向量是一个有序数组,可以表示为:

v=[v1v2vn]\mathbf{v} = \begin{bmatrix} v_1 \\ v_2 \\ \vdots \\ v_n \end{bmatrix}

其中 v1,v2,,vnv_1, v_2, \ldots, v_n 是向量 v\mathbf{v} 的分量。在 Rn\mathbb{R}^n 中,向量有 nn 个分量。

这里解释一下 Rn\mathbb{R}^n 的含义:表示 nn 维实数空间,其中的向量有 nn 个分量,每个分量都是实数

所有元素都是零的向量称为零向量,记为 0\mathbf{0}

1.3.2 向量的运算

向量加法:两个向量相加,对应分量相加。

u+v=[u1u2un]+[v1v2vn]=[u1+v1u2+v2un+vn]\mathbf{u} + \mathbf{v} = \begin{bmatrix} u_1 \\ u_2 \\ \vdots \\ u_n \end{bmatrix} + \begin{bmatrix} v_1 \\ v_2 \\ \vdots \\ v_n \end{bmatrix} = \begin{bmatrix} u_1 + v_1 \\ u_2 + v_2 \\ \vdots \\ u_n + v_n \end{bmatrix}

标量乘法:向量乘以标量,每个分量都乘以该标量。

cv=c[v1v2vn]=[cv1cv2cvn]c\mathbf{v} = c\begin{bmatrix} v_1 \\ v_2 \\ \vdots \\ v_n \end{bmatrix} = \begin{bmatrix} cv_1 \\ cv_2 \\ \vdots \\ cv_n \end{bmatrix}

向量加法满足以下性质:

  1. 交换律:u+v=v+u\mathbf{u} + \mathbf{v} = \mathbf{v} + \mathbf{u}
  2. 结合律:(u+v)+w=u+(v+w)(\mathbf{u} + \mathbf{v}) + \mathbf{w} = \mathbf{u} + (\mathbf{v} + \mathbf{w})

标量乘法满足以下性质:

  1. c(dv)=(cd)vc(d\mathbf{v}) = (cd)\mathbf{v}
  2. (c+d)v=cv+dv(c+d)\mathbf{v} = c\mathbf{v} + d\mathbf{v}
  3. c(u+v)=cu+cvc(\mathbf{u} + \mathbf{v}) = c\mathbf{u} + c\mathbf{v}

1.3.3 线性组合

给定向量 v1,v2,,vp\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_p 和标量 c1,c2,,cpc_1, c_2, \ldots, c_p,表达式

c1v1+c2v2++cpvpc_1\mathbf{v}_1 + c_2\mathbf{v}_2 + \cdots + c_p\mathbf{v}_p

称为向量 v1,v2,,vp\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_p线性组合,其中 c1,c2,,cpc_1, c_2, \ldots, c_p权重系数

这里我们引入一道例题来分析存在性问题:

例题

a1=[125],a2=[256],b=[743]a_1 = \begin{bmatrix} 1 \\ -2 \\ -5 \end{bmatrix}, \quad a_2 = \begin{bmatrix} 2 \\ 5 \\ 6 \end{bmatrix}, \quad b = \begin{bmatrix} 7 \\ 4 \\ -3 \end{bmatrix}

确定 bb 是否能写成 a1\mathbf{a_1}a2\mathbf{a_2} 的线性组合,也就是说,确定是否存在权 x1x_1x2x_2,使得

x1a1+x2a2=bx_1 \mathbf{a_1} + x_2 \mathbf{a_2} = b

若向量方程有解,求它的解。

解:

首先我们将向量方程写成线性方程组的形式:

x1+2x2=72x1+5x2=45x1+6x2=3\begin{aligned} x_1 + 2x_2 &= 7 \\ -2x_1 + 5x_2 &= 4 \\ -5x_1 + 6x_2 &= -3 \end{aligned}

然后写出对应的增广矩阵:

[127254563]\begin{bmatrix} 1 & 2 & 7 \\ -2 & 5 & 4 \\ -5 & 6 & -3 \end{bmatrix}

我们使用行简化算法将其化为简化阶梯形:

[103012000]\begin{bmatrix} 1 & 0 & 3 \\ 0 & 1 & 2 \\ 0 & 0 & 0 \end{bmatrix}

这里说明有唯一解为:

x1=3,x2=2x_1 = 3, \quad x_2 = 2

因此我们可以得出以下结论:

向量方程

x1a1+x2a2++xnan=bx_1\mathbf{a}_1 + x_2\mathbf{a}_2 + \cdots + x_n\mathbf{a}_n = \mathbf{b}

和增广矩阵为

[a1a2anb]\begin{bmatrix} \mathbf{a}_1 & \mathbf{a}_2 & \cdots & \mathbf{a}_n & \mathbf{b} \end{bmatrix}

的线性方程组有相同的解集。特别地,b\mathbf{b} 可以表示为 a1,a2,,an\mathbf{a}_1, \mathbf{a}_2, \ldots, \mathbf{a}_n 的线性组合且当且仅当对应于上式的线性方程组有解。

1.3.4 张成空间

给定向量 v1,v2,,vp\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_p,所有这些向量的线性组合的集合称为由 v1,v2,,vp\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_p 张成张成空间,记作 Span{v1,v2,,vp}\text{Span}\{\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_p\}。也就是说,Span{v1,v2,,vp}\text{Span}\{\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_p\}是所有形如

c1v1+c2v2++cpvpc_1\mathbf{v}_1 + c_2\mathbf{v}_2 + \cdots + c_p\mathbf{v}_p

的向量的集合,其中 c1,c2,,cpc_1, c_2, \ldots, c_p 为标量。

要判断向量 b\mathbf{b} 是否在张成空间 Span{v1,v2,,vp}\text{Span}\{\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_p\} 中,只需要判断向量方程

x1v1+x2v2++xpvp=bx_1\mathbf{v}_1 + x_2\mathbf{v}_2 + \cdots + x_p\mathbf{v}_p = \mathbf{b}

是否有解即可,或者等价地,我们判断增广矩阵

[v1v2vpb]\begin{bmatrix} \mathbf{v}_1 & \mathbf{v}_2 & \cdots & \mathbf{v}_p & \mathbf{b} \end{bmatrix}

线性方程组是否有解即可。

注意:Span{v1,v2,,vp}\text{Span}\{\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_p\} 包含了 vi\mathbf{v_i} 的所有倍数,因为

civi=civi+0v1++0vi1+0vi+1++0vpc_i\mathbf{v_i}=c_i\mathbf{v_i} + 0\mathbf{v_1} + \cdots + 0\mathbf{v_{i-1}} + 0\mathbf{v_{i+1}} + \cdots + 0\mathbf{v_p}

特别地,一定包含零向量

这里我们还给出Span{v}\text{Span}\{\mathbf{v}\}Span{u,v}\text{Span}\{\mathbf{u,v}\}的几何解释:

v\mathbf{v}R3\mathbb{R}^3中的一个非零向量,则Span{v}\text{Span}\{\mathbf{v}\}是通过原点和v\mathbf{v}的直线上所有点的集合,如下图所示:

alt text

u\mathbf{u}v\mathbf{v}R3\mathbb{R}^3中的两个不共线的非零向量,则Span{u,v}\text{Span}\{\mathbf{u,v}\}是通过原点、u\mathbf{u}v\mathbf{v}的平面上所有点的集合,如下图所示:

alt text

例题

a1=[123]\mathbf{a}_1=\begin{bmatrix}1\\-2\\3\end{bmatrix}a2=[5133]\mathbf{a}_2=\begin{bmatrix}5\\-13\\-3\end{bmatrix}b=[381]\mathbf{b}=\begin{bmatrix}-3\\8\\1\end{bmatrix},则 Span{a1,a2}\operatorname{Span}\left\{\mathbf{a}_1,\mathbf{a}_2\right\}R3\mathbb{R}^{3} 中通过原点的一个平面,问 b\mathbf{b} 是否在该平面上?

解:

我们判断向量方程

x1a1+x2a2=bx_1\mathbf{a}_1 + x_2\mathbf{a}_2 = \mathbf{b}

是否有解即可,或者等价地,我们判断增广矩阵为

[a1a2b]\begin{bmatrix} \mathbf{a}_1 & \mathbf{a}_2 & \mathbf{b} \end{bmatrix}

的线性方程组是否有解即可。

这里我们代入增广矩阵并进行化简:

[1532138331]\begin{bmatrix} 1 & 5 & -3 \\ -2 & -13 & 8 \\ 3 & -3 & 1 \end{bmatrix}

化简为阶梯形如下:

[153032002]\begin{bmatrix} 1 & 5 & -3 \\ 0 & -3 & 2 \\ 0 & 0 & -2 \end{bmatrix}

明显看出向量方程无解,因此 b\mathbf{b} 不在 Span{a1,a2}\operatorname{Span}\left\{\mathbf{a}_1,\mathbf{a}_2\right\} 中。

1.3.5 向量方程与线性方程组的等价性

考虑向量方程:

x1a1+x2a2++xnan=bx_1\mathbf{a}_1 + x_2\mathbf{a}_2 + \cdots + x_n\mathbf{a}_n = \mathbf{b}

其中 a1,a2,,an\mathbf{a}_1, \mathbf{a}_2, \ldots, \mathbf{a}_nb\mathbf{b}Rm\mathbb{R}^m 中的向量,x1,x2,,xnx_1, x_2, \ldots, x_n 是未知数。

这个向量方程等价于一个 m×nm \times n 的线性方程组。如果我们将向量 aj\mathbf{a}_j 表示为:

aj=[a1ja2jamj]\mathbf{a}_j = \begin{bmatrix} a_{1j} \\ a_{2j} \\ \vdots \\ a_{mj} \end{bmatrix}

那么向量方程可以写成:

x1[a11a21am1]+x2[a12a22am2]++xn[a1na2namn]=[b1b2bm]x_1\begin{bmatrix} a_{11} \\ a_{21} \\ \vdots \\ a_{m1} \end{bmatrix} + x_2\begin{bmatrix} a_{12} \\ a_{22} \\ \vdots \\ a_{m2} \end{bmatrix} + \cdots + x_n\begin{bmatrix} a_{1n} \\ a_{2n} \\ \vdots \\ a_{mn} \end{bmatrix} = \begin{bmatrix} b_1 \\ b_2 \\ \vdots \\ b_m \end{bmatrix}

这等价于线性方程组:

a11x1+a12x2++a1nxn=b1a21x1+a22x2++a2nxn=b2am1x1+am2x2++amnxn=bm\begin{aligned} a_{11}x_1 + a_{12}x_2 + \cdots + a_{1n}x_n &= b_1 \\ a_{21}x_1 + a_{22}x_2 + \cdots + a_{2n}x_n &= b_2 \\ \vdots \\ a_{m1}x_1 + a_{m2}x_2 + \cdots + a_{mn}x_n &= b_m \end{aligned}

1.3.6 向量方程的几何解释

向量方程 x1a1+x2a2++xnan=bx_1\mathbf{a}_1 + x_2\mathbf{a}_2 + \cdots + x_n\mathbf{a}_n = \mathbf{b} 有解的充分必要条件是 b\mathbf{b} 可以表示为 a1,a2,,an\mathbf{a}_1, \mathbf{a}_2, \ldots, \mathbf{a}_n 的线性组合,即 b\mathbf{b} 在由 a1,a2,,an\mathbf{a}_1, \mathbf{a}_2, \ldots, \mathbf{a}_n 张成的向量空间中。

从几何角度看:

  • R2\mathbb{R}^2 中,两个不共线的向量可以张成整个平面
  • R3\mathbb{R}^3 中,三个不共面的向量可以张成整个空间

1.3.7 齐次向量方程

b=0\mathbf{b} = \mathbf{0} 时,向量方程 x1a1+x2a2++xnan=0x_1\mathbf{a}_1 + x_2\mathbf{a}_2 + \cdots + x_n\mathbf{a}_n = \mathbf{0} 称为齐次向量方程

齐次向量方程总是有平凡解 x1=x2==xn=0x_1 = x_2 = \cdots = x_n = 0。如果存在非零解,则称向量 a1,a2,,an\mathbf{a}_1, \mathbf{a}_2, \ldots, \mathbf{a}_n线性相关的;否则,称它们是线性无关的。


1.4 矩阵方程 Ax=b

1.4.1 向量方程和矩阵方程

定义:

AA 是一个 m×nm \times n 矩阵,它的各列为 a1,a2,,an\mathbf{a}_1, \mathbf{a}_2, \ldots, \mathbf{a}_n,则 AAx\mathbf{x} 的积(记为 AxA\mathbf{x})就是 AA 的各列以 x\mathbf{x} 中对应元素为权的线性组合,即

Ax=[a1a2an][x1x2xn]=x1a1+x2a2++xnanA\mathbf{x} = \begin{bmatrix} \mathbf{a}_1 & \mathbf{a}_2 & \cdots & \mathbf{a}_n \end{bmatrix} \begin{bmatrix} x_1 \\ x_2 \\ \vdots \\ x_n \end{bmatrix} = x_1\mathbf{a}_1 + x_2\mathbf{a}_2 + \cdots + x_n\mathbf{a}_n

注意这里 AxA\mathbf{x} 当且仅当 AA 的列数等于 x\mathbf{x} 的行数时才有定义。

形式如 x1a1+x2a2++xnan=bx_1\mathbf{a}_1 + x_2\mathbf{a}_2 + \cdots + x_n\mathbf{a}_n = \mathbf{b} 的方程我们称为向量方程

形式如 Ax=bA\mathbf{x}=\mathbf{b} 的方程我们称为矩阵方程。’

定理1:
AAm×nm \times n 矩阵,它的各列为 a1,,an\boldsymbol{a}_{1}, \cdots, \boldsymbol{a}_{n},而 b\boldsymbol{b} 属于 Rm\mathbb{R}^{m},则矩阵方程

Ax=bA\boldsymbol{x} = \boldsymbol{b}

向量方程

x1a1+x2a2++xnan=bx_{1}\boldsymbol{a}_{1} + x_{2}\boldsymbol{a}_{2} + \cdots + x_{n}\boldsymbol{a}_{n} = \boldsymbol{b}

有相同的解集。它又与增广矩阵

[a1a2anb] \left[\begin{array}{llll} \boldsymbol{a}_{1} & \boldsymbol{a}_{2} & \cdots & \boldsymbol{a}_{n} & \boldsymbol{b} \end{array}\right]

线性方程组有相同的解集。

1.4.2 解的存在性

方程 Ax=bA\mathbf{x}=\mathbf{b} 有解当且仅当 b\mathbf{b} 可以表示为 AA 的各列的线性组合。

定理2:

AAm×nm \times n 矩阵,则下列命题是逻辑上等价的(即对某个 AA,它们都成立或都不成立):

  1. Rm\mathbb{R}^{m} 中每个 b\boldsymbol{b},方程 Ax=bA\boldsymbol{x} = \boldsymbol{b} 有解。
  2. Rm\mathbb{R}^{m} 中的每个 b\boldsymbol{b} 都是 AA 的列的一个线性组合。
  3. AA 的各列生成 Rm\mathbb{R}^{m}
  4. AA 在每一行都有一个主元位置。

注意:这里讨论的是系数矩阵而非增广矩阵。如果是对应增广矩阵,如果每一行都有一个主元位置,并不一定有解。

1.4.3 Ax 的计算

若乘积 AxA\mathbf{x} 有定义,则 AxA\mathbf{x} 的第 ii 个分量是矩阵 AA 的第 ii 行与列向量 x\mathbf{x} 的点积。

这里需要注意,如果 A=IA=I,即为单位阵,这时候 Ax=xA\mathbf{x} = \mathbf{x}

1.4.4 矩阵-向量积 Ax 的性质

AAm×nm \times n 矩阵,u\mathbf{u}v\mathbf{v}nn 维向量,cc 是标量,则

A(u+v)=Au+AvA(cu)=cAuA(\mathbf{u} + \mathbf{v}) = A\mathbf{u} + A\mathbf{v} \\ A(c\mathbf{u}) = cA\mathbf{u}


1.5 线性方程组的解集

1.5.1 齐次线性方程组

齐次线性方程组是形如 Ax=0A\mathbf{x} = \mathbf{0} 的方程组,其中 AA 是一个 m×nm \times n 矩阵,x\mathbf{x}nn 维向量,0\mathbf{0}Rm\mathbb{R}^m 中的零向量。

齐次线性方程组的一个重要性质是它总是有解,至少有零解(或称为平凡解x=0\mathbf{x} = \mathbf{0}

定理3: 齐次线性方程组 Ax=0A\mathbf{x} = \mathbf{0} 有非平凡解的充分必要条件是方程组中的自由变量至少有一个,或等价地,AA列数大于行秩(即 n>rn > r,其中 rrAA 的秩)。

注意:非平凡解向量 x\mathbf{x} 可能会有一些零元素,但只要不是全为零,就称为非平凡解。

强调:齐次方程 Ax=0A\mathbf{x} = \mathbf{0} 总可以表示为张成空间的形式,即

Nul(A)=span{v1,v2,,vp}\text{Nul}(A) = \text{span}\{\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_p\}

这里 v1,v2,,vp\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_pAx=0A\mathbf{x} = \mathbf{0} 的一组基,ppAA 的零空间的维度,也称为 AA零度(nullity),满足 p=nrp = n - r,其中 nnAA 的列数,rrAA 的行秩。

齐次方程通解的参数向量形式为:

x=c1v1+c2v2++cpvp\mathbf{x} = c_1\mathbf{v}_1 + c_2\mathbf{v}_2 + \cdots + c_p\mathbf{v}_p

其中 c1,c2,,cpc_1, c_2, \ldots, c_p 是任意常数。

1.5.2 非齐次线性方程组

非齐次线性方程组是形如 Ax=bA\mathbf{x} = \mathbf{b} 的方程组,其中 b0\mathbf{b} \neq \mathbf{0}

非齐次线性方程组的解集可以表示为:

{p+v:vNul(A)}\{\mathbf{p} + \mathbf{v} : \mathbf{v} \in \text{Nul}(A)\}

其中 p\mathbf{p} 是非齐次方程组的一个特解,Nul(A)\text{Nul}(A) 是对应齐次方程组 Ax=0A\mathbf{x} = \mathbf{0} 的解集,也称为 AA零空间

定理4: 如果 p\mathbf{p} 是非齐次线性方程组 Ax=bA\mathbf{x} = \mathbf{b} 的一个特解,而 v1,v2,,vp\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_p 构成对应齐次方程组 Ax=0A\mathbf{x} = \mathbf{0} 解空间的一组基,则 Ax=bA\mathbf{x} = \mathbf{b} 的通解可以表示为:

x=p+c1v1+c2v2++cpvp\mathbf{x} = \mathbf{p} + c_1\mathbf{v}_1 + c_2\mathbf{v}_2 + \cdots + c_p\mathbf{v}_p

其中 c1,c2,,cpc_1, c_2, \ldots, c_p 是任意常数。

1.5.3 解集的参数表示

线性方程组的解集可以用参数表示。对于齐次线性方程组,解集可以表示为:

x=c1v1+c2v2++cpvp\mathbf{x} = c_1\mathbf{v}_1 + c_2\mathbf{v}_2 + \cdots + c_p\mathbf{v}_p

其中 v1,v2,,vp\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_p 是零空间的一组基,c1,c2,,cpc_1, c_2, \ldots, c_p 是任意常数。

对于非齐次线性方程组,解集可以表示为:

x=p+c1v1+c2v2++cpvp\mathbf{x} = \mathbf{p} + c_1\mathbf{v}_1 + c_2\mathbf{v}_2 + \cdots + c_p\mathbf{v}_p

其中 p\mathbf{p} 是一个特解。

1.5.4 解集的几何解释

线性方程组的解集在几何上有直观的解释:

  • 对于齐次线性方程组 Ax=0A\mathbf{x} = \mathbf{0},解集是一个过原点的子空间。

    • 如果只有零解,则解集仅为原点。
    • 如果有一个自由变量,解集是一条过原点的直线。
    • 如果有两个自由变量,解集是一个过原点的平面。
    • 依此类推。
  • 对于非齐次线性方程组 Ax=bA\mathbf{x} = \mathbf{b},解集是一个不过原点的仿射子空间(即平行于某个子空间的平面)。

    • 如果对应的齐次方程组只有零解,则解集是一个点。
    • 如果对应的齐次方程组有一个自由变量,解集是一条不过原点的直线。
    • 如果对应的齐次方程组有两个自由变量,解集是一个不过原点的平面。
    • 依此类推。

1.5.5 例题

例题1: 求解线性方程组

x12x2+x3=42x15x2+2x3=7x14x2+3x3=10\begin{aligned} x_1 - 2x_2 + x_3 &= 4 \\ 2x_1 - 5x_2 + 2x_3 &= 7 \\ x_1 - 4x_2 + 3x_3 &= 10 \end{aligned}

解: 首先写出增广矩阵并进行行简化:

[1214252714310]\begin{bmatrix} 1 & -2 & 1 & 4 \\ 2 & -5 & 2 & 7 \\ 1 & -4 & 3 & 10 \end{bmatrix}

通过行变换将其化为简化阶梯形:

[121401010013]\begin{bmatrix} 1 & -2 & 1 & 4 \\ 0 & -1 & 0 & -1 \\ 0 & 0 & 1 & 3 \end{bmatrix}

从简化阶梯形可以得到:

x12x2+x3=4x2=1x3=3\begin{aligned} x_1 - 2x_2 + x_3 &= 4 \\ -x_2 &= -1 \\ x_3 &= 3 \end{aligned}

解得:x2=1x_2 = 1x3=3x_3 = 3,代入第一个方程得 x121+3=4x_1 - 2 \cdot 1 + 3 = 4,即 x1=3x_1 = 3

因此,方程组的唯一解为 x1=3x_1 = 3x2=1x_2 = 1x3=3x_3 = 3,或写成向量形式 x=[313]\mathbf{x} = \begin{bmatrix} 3 \\ 1 \\ 3 \end{bmatrix}

例题2: 求解齐次线性方程组

x13x2+x3=02x16x2+2x3=0x1+3x2x3=0\begin{aligned} x_1 - 3x_2 + x_3 &= 0 \\ 2x_1 - 6x_2 + 2x_3 &= 0 \\ -x_1 + 3x_2 - x_3 &= 0 \end{aligned}

解: 首先写出增广矩阵并进行行简化:

[131026201310]\begin{bmatrix} 1 & -3 & 1 & 0 \\ 2 & -6 & 2 & 0 \\ -1 & 3 & -1 & 0 \end{bmatrix}

通过行变换将其化为简化阶梯形:

[131000000000]\begin{bmatrix} 1 & -3 & 1 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{bmatrix}

从简化阶梯形可以得到:

x13x2+x3=0x_1 - 3x_2 + x_3 = 0

这里 x2x_2x3x_3 是自由变量,可以任意取值。令 x2=sx_2 = sx3=tx_3 = t,则 x1=3stx_1 = 3s - t

因此,方程组的通解对应的参数向量形式为:

x=[3stst]=s[310]+t[101]\mathbf{x} = \begin{bmatrix} 3s - t \\ s \\ t \end{bmatrix} = s\begin{bmatrix} 3 \\ 1 \\ 0 \end{bmatrix} + t\begin{bmatrix} -1 \\ 0 \\ 1 \end{bmatrix}

其中 sstt 是任意实数。

这里我们令 u=[310]\mathbf{u}=\begin{bmatrix} 3 \\ 1 \\ 0 \end{bmatrix}v=[101]\mathbf{v}=\begin{bmatrix} -1 \\ 0 \\ 1 \end{bmatrix},平凡解只需要 s=0s=0t=0t=0 即可。

几何意义下,解集是R3\mathbb{R}^3中过原点的平面,即为Span{u,v}\text{Span}\{\mathbf{u},\mathbf{v}\}u\mathbf{u}v\mathbf{v} 是该平面的一组基向量。


1.6 线性无关

1.6.1 线性无关的定义

定义:

v1,v2,,vp\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_pRn\mathbb{R}^n 中的一组向量,如果方程 c1v1+c2v2++cpvp=0c_1\mathbf{v}_1 + c_2\mathbf{v}_2 + \cdots + c_p\mathbf{v}_p = \mathbf{0} 的唯一解是 c1=c2==cp=0c_1 = c_2 = \cdots = c_p = 0(即仅有平凡解),则称向量 v1,v2,,vp\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_p线性无关的;否则,称它们是线性相关的

这里也可以理解为,如果矩阵方程 Ax=0A\mathbf{x}=\mathbf{0} 只有平凡解,那么矩阵 AA 的各列是线性无关的。

1.6.2 矩阵各列的线性无关

定理5:

AAm×nm \times n 矩阵,它的各列为 a1,a2,,an\mathbf{a}_1, \mathbf{a}_2, \ldots, \mathbf{a}_n,则下列命题是逻辑上等价的(即对某个 AA,它们都成立或都不成立):

  1. 矩阵方程 Ax=0A\mathbf{x} = \mathbf{0} 仅有平凡解。
  2. 向量方程 x1a1+x2a2++xnan=0x_1\mathbf{a}_1 + x_2\mathbf{a}_2 + \cdots + x_n\mathbf{a}_n = \mathbf{0} 仅有平凡解。
  3. 向量 a1,a2,,an\mathbf{a}_1, \mathbf{a}_2, \ldots, \mathbf{a}_n 是线性无关的。

1.6.3 一个或两个向量的集合

定理6:

仅含一个向量的集合线性无关当且仅当该向量非零。因为当 v0\mathbf{v} \neq \mathbf{0} 时,向量方程 x1v=0x_1\mathbf{v} = \mathbf{0} 的唯一解是 x1=0x_1 = 0。而零向量是线性相关的 ,因为向量方程 00=00\mathbf{0} = \mathbf{0} 有非平凡解 x1=1x_1 = 1

定理7:

两个向量 v1,v2\mathbf{v}_1, \mathbf{v}_2 是线性无关的当且仅当 v1\mathbf{v}_1 不是 v2\mathbf{v}_2 的倍数。

我们从几何意义上看,两个向量线性相关,当且仅当它们落在通过原点的同一直线上,如下图所示:

alt text

1.6.4 两个或更多向量的集合

定理8:

v1,v2,,vp\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_pRn\mathbb{R}^n 中的一组向量,集合 S={v1,v2,,vp}S=\{\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_p\} 线性相关,当且仅当 SS至少有一个向量是其他向量的线性组合。事实上,若 SS 线性相关,且 v10\mathbf{v_1} \neq \mathbf{0},则某个向量 vj(j>1)\mathbf{v_j}(j>1) 是它前面向量 v1,v2,,vj1\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_{j-1} 的线性组合。

证明:

S={v1,v2,,vp}S = \{\boldsymbol{v}_1, \boldsymbol{v}_2, \ldots, \boldsymbol{v}_p\} 为向量集合,则:

  1. 充分性证明
    若存在 vj\boldsymbol{v}_j 可表示为其他向量的线性组合,即

    vj=ijcivi(ciR)\boldsymbol{v}_j = \sum_{i \neq j} c_i \boldsymbol{v}_i \quad (c_i \in \mathbb{R})

    通过移项可得非平凡线性关系:

    (1)vj+ijcivi=0(-1)\boldsymbol{v}_j + \sum_{i \neq j} c_i \boldsymbol{v}_i = \boldsymbol{0}

    SS 线性相关。

  2. 必要性证明
    SS 线性相关,则存在不全为零的标量 c1,,cpc_1, \ldots, c_p 使得:

    i=1pcivi=0\sum_{i=1}^p c_i \boldsymbol{v}_i = \boldsymbol{0}

    • 情形1​:若 v1=0\boldsymbol{v}_1 = \boldsymbol{0},则 v1\boldsymbol{v}_1 可平凡表示为其他向量的线性组合(系数全零)且此时 c10c_1 \neq 0
    • 情形2​:若 v10\boldsymbol{v}_1 \neq \boldsymbol{0},设 jj 是最大下标使得 cj0c_j \neq 0
      • j=1j=1 时,c1v1=0c_1 \boldsymbol{v}_1 = \boldsymbol{0}v10\boldsymbol{v}_1 \neq \boldsymbol{0} 矛盾;
      • j>1j>1 时,可解出 vj\boldsymbol{v}_j 为前 j1j-1 个向量的线性组合:

        vj=i=1j1cicjvi\boldsymbol{v}_j = -\sum_{i=1}^{j-1} \frac{c_i}{c_j} \boldsymbol{v}_i

定理9:

若一个向量组的向量个数超过每个向量的元素个数,那么这个向量组线性相关,就是说,Rn\mathbb{R}^n 中任意向量组 {v1,,vn+1}\{ \mathbf{v_1, \ldots, v_{n+1}} \} 当满足 p>np>n 时线性相关。

这里理解为:如果满足定理要求,则线性方程组中未知数个数大于方程个数,因此有无数个非平凡解。

定理10:

Rn\mathbb{R}^n 中向量组 {v1,,vp}\{ \mathbf{v_1, \ldots, v_p} \} 包含零向量,则它线性相关。


1.7 线性变换介绍

1.7.1 线性变换的定义

线性变换是一种特殊的函数,它将一个向量空间的元素映射到另一个向量空间,并保持向量加法和标量乘法的性质。

T:RnRmT: \mathbb{R}^n \rightarrow \mathbb{R}^m 是一个函数,如果对于 Rn\mathbb{R}^n 中的任意向量 u\mathbf{u}v\mathbf{v} 以及任意标量 cc,满足以下两个条件:

  1. T(u+v)=T(u)+T(v)T(\mathbf{u} + \mathbf{v}) = T(\mathbf{u}) + T(\mathbf{v})(加法保持性)
  2. T(cu)=cT(u)T(c\mathbf{u}) = cT(\mathbf{u})(标量乘法保持性)

则称 TT 为从 Rn\mathbb{R}^nRm\mathbb{R}^m线性变换

这里我们称 Rn\mathbb{R}^n定义域,而 Rm\mathbb{R}^m余定义域(也称为取值空间)。

这里也很容易推出以下两条重要性质

  1. TT 是线性变换,则 T(0)=0T(\mathbf{0}) = \mathbf{0}
  2. 且对 TT 定义域中一切向量 u\mathbf{u}v\mathbf{v} 以及任意标量 ccdd,有 T(cu+dv)=cT(u)+dT(v)T(c\mathbf{u} + d\mathbf{v}) = cT(\mathbf{u}) + dT(\mathbf{v})

进一步推广我们可以得到叠加原理

T:RnRmT: \mathbb{R}^n \rightarrow \mathbb{R}^m 是线性变换,则对 Rn\mathbb{R}^n 中的任意向量 u1,u2,,up\mathbf{u}_1, \mathbf{u}_2, \ldots, \mathbf{u}_p,以及标量 c1,c2,,cpc_1, c_2, \ldots, c_p,有

T(c1u1+c2u2++cpup)=c1T(u1)+c2T(u2)++cpT(up)T(c_1\mathbf{u}_1 + c_2\mathbf{u}_2 + \cdots + c_p\mathbf{u}_p) = c_1T(\mathbf{u}_1) + c_2T(\mathbf{u}_2) + \cdots + c_pT(\mathbf{u}_p)

1.7.2 线性变换的矩阵表示

定理11:

每个从 Rn\mathbb{R}^nRm\mathbb{R}^m 的线性变换 TT 都可以用一个 m×nm \times n 的矩阵 AA 来表示,使得对于 Rn\mathbb{R}^n 中的任意向量 x\mathbf{x},都有:

T(x)=AxT(\mathbf{x}) = A\mathbf{x}

矩阵 AA 的列向量可以通过计算标准基向量的像来确定:

A=[T(e1)T(e2)T(en)]A = \begin{bmatrix} T(\mathbf{e}_1) & T(\mathbf{e}_2) & \cdots & T(\mathbf{e}_n) \end{bmatrix}

其中 ej\mathbf{e}_jRn\mathbb{R}^n 中的第 jj标准基向量(第 jj 个分量为 1,其余分量为 0),即对应于单位矩阵 InI_n 的第 jj 列。

证明
x=Inx=[e1en]x=x1e1++xnen\boldsymbol{x} = I_n \boldsymbol{x} = [\boldsymbol{e}_1 \cdots \boldsymbol{e}_n] \boldsymbol{x} = x_1 \boldsymbol{e}_1 + \cdots + x_n \boldsymbol{e}_n,由于 TT 是线性变换,可得:

T(x)=T(x1e1++xnen)=x1T(e1)++xnT(en)=[T(e1)T(en)][x1xn]=Ax\begin{aligned} T(\boldsymbol{x}) &= T(x_1 \boldsymbol{e}_1 + \cdots + x_n \boldsymbol{e}_n) \\ &= x_1 T(\boldsymbol{e}_1) + \cdots + x_n T(\boldsymbol{e}_n) \\ &= \left[ T(\boldsymbol{e}_1) \cdots T(\boldsymbol{e}_n) \right] \begin{bmatrix} x_1 \\ \vdots \\ x_n \end{bmatrix} \\ &= A \boldsymbol{x} \end{aligned}

例题:

例3
T:R2R2T:\mathbb{R}^2 \rightarrow \mathbb{R}^2 为将 R2\mathbb{R}^2 中每个点绕原点逆时针旋转角度 φ\varphi 的变换。通过几何分析可证明这是线性变换(参见图1-39)。求该变换的标准矩阵。

  1. 基向量旋转结果:
    • [10]\begin{bmatrix}1\\0\end{bmatrix} 旋转后为 [cosφsinφ]\begin{bmatrix}\cos\varphi\\\sin\varphi\end{bmatrix}
    • [01]\begin{bmatrix}0\\1\end{bmatrix} 旋转后为 [sinφcosφ]\begin{bmatrix}-\sin\varphi\\\cos\varphi\end{bmatrix}

alt text

  1. 定理11​(线性变换的标准矩阵构造法),直接得到:

    A=[cosφsinφsinφcosφ]A = \begin{bmatrix} \cos\varphi & -\sin\varphi \\ \sin\varphi & \cos\varphi \end{bmatrix}

1.7.3 线性变换的几何解释

线性变换可以从几何角度理解为对向量空间的一种变形,它保持了向量间的线性关系。以下是一些常见线性变换的几何解释:

  1. 伸缩变换:将向量沿着某个方向拉伸或压缩。
    例如,矩阵 [2003]\begin{bmatrix} 2 & 0 \\ 0 & 3 \end{bmatrix} 表示在 xx 方向上拉伸为原来的 2 倍,在 yy 方向上拉伸为原来的 3 倍。

  2. 旋转变换:将向量绕原点旋转一定角度。
    例如,矩阵 [cosθsinθsinθcosθ]\begin{bmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{bmatrix} 表示将向量逆时针旋转 θ\theta 角度。

  3. 投影变换:将向量投影到某个子空间上。
    例如,矩阵 [1000]\begin{bmatrix} 1 & 0 \\ 0 & 0 \end{bmatrix} 表示将向量投影到 xx 轴上。

  4. 反射变换:将向量关于某个子空间反射。
    例如,矩阵 [1001]\begin{bmatrix} 1 & 0 \\ 0 & -1 \end{bmatrix} 表示将向量关于 xx 轴反射。

  5. 剪切变换:将向量沿着某个方向剪切。
    例如,矩阵 [1k01]\begin{bmatrix} 1 & k \\ 0 & 1 \end{bmatrix} 表示将向量沿 xx 方向剪切,剪切因子为 kk

更详细的描述见下表展示:

alt text

alt text

alt text

alt text

1.7.4 线性变换的存在性与唯一性问题

定义:

映射 TT:若 Rm\mathbb{R}^m 中每个 b\mathbf{b}Rn\mathbb{R}^n至少一个的像,则RnRm\mathbb{R}^n \rightarrow \mathbb{R}^m 称为到 Rm\mathbb{R}^m 上的映射(也称为满射)。

我们也根据判断 TT 的值域是否是整个余定义域 R\mathbb{R} 来判断是否是满射,如下图所示更清晰一点:

alt text

定义:

映射 TT:若 Rm\mathbb{R}^m 中每个 b\mathbf{b}Rn\mathbb{R}^n至多一个的像,则RnRm\mathbb{R}^n \rightarrow \mathbb{R}^m 称为到 Rm\mathbb{R}^m 的单射(也称为一对一映射)。

alt text

定理12:

T:RnRmT: \mathbb{R}^n \rightarrow \mathbb{R}^m 是线性变换,则 TT单射当且仅当方程 Ax=0A\mathbf{x}=\mathbf{0} 仅有平凡解

​**证明:**​

  1. 充分性:因 TT 是线性的,故 T(0)=T(0x)=0T(x)=0T(\boldsymbol{0}) = T(0\mathbf{x}) = 0T(\mathbf{x}) = \boldsymbol{0}

  2. 必要性:

    • TT 是一对一的
      • 则方程 T(x)=0T(\boldsymbol{x}) = \boldsymbol{0} 至多有一个解,因此仅有平凡解 x=0\boldsymbol{x} = \boldsymbol{0}
    • TT 不是一对一的
      • 存在 Rn\mathbb{R}^n 中两个相异向量 u\boldsymbol{u}v\boldsymbol{v},使得 T(u)=bT(\boldsymbol{u}) = \boldsymbol{b}T(v)=bT(\boldsymbol{v}) = \boldsymbol{b}
      • 由线性性可得:

      T(uv)=T(u)T(v)=bb=0T(\boldsymbol{u} - \boldsymbol{v}) = T(\boldsymbol{u}) - T(\boldsymbol{v}) = \boldsymbol{b} - \boldsymbol{b} = \boldsymbol{0}

      这里又因为 TT 是线性的且方程 Ax=0A\mathbf{x}=\mathbf{0} 仅有平凡解,所以向量 uv=0\mathbf{u-v} = \mathbf{0},即 u=v\mathbf{u} = \mathbf{v},这与假设矛盾。

例题
TT 是线性变换,它的标准矩阵为

A=[148102130005]A = \begin{bmatrix} 1 & -4 & 8 & 1 \\ 0 & 2 & -1 & 3 \\ 0 & 0 & 0 & 5 \end{bmatrix}

  1. ​​TT 是否把 R4\mathbb{R}^{4} 映上到 R3\mathbb{R}^{3}?​
  2. ​​TT 是否是一对一映射?


AA 已是阶梯形,可直接观察到:

  • AA 在每一行都有主元位置(由1.4.2节的定理2
  • 对任意 bR3\boldsymbol{b} \in \mathbb{R}^{3},方程 Ax=bA\boldsymbol{x} = \boldsymbol{b} 相容
    TTR4\mathbb{R}^{4} 映射到 R3\mathbb{R}^{3}

但方程 Ax=bA\boldsymbol{x} = \boldsymbol{b} 存在自由变量(4变量中仅有3个基本变量),故:

  • 每个 b\boldsymbol{b} 对应多个 x\boldsymbol{x}
    TT 不是一对一映射

定理13:

T:RnRmT:\mathbb{R}^{n}\rightarrow\mathbb{R}^{m} 是线性变换,AATT 的标准矩阵,则:

  1. TTRn\mathbb{R}^{n} 映上到 Rm\mathbb{R}^{m},当且仅当 AA 的列生成 Rm\mathbb{R}^{m}
  2. TT 是一对一的,当且仅当 AA 的列线性无关。

例题:

例5
设线性变换 T(x1,x2)=(3x1+x2,5x1+7x2,x1+3x2)T(x_1,x_2) = (3x_1 + x_2, 5x_1 + 7x_2, x_1 + 3x_2),证明 TT 是一对一线性变换。TT 是否将 R2\mathbb{R}^2 映射到 R3\mathbb{R}^3


T(x)T(\boldsymbol{x}) 表示为矩阵形式:

T(x)=[3x1+x25x1+7x2x1+3x2]=[315713][x1x2]=AxT(\boldsymbol{x}) = \begin{bmatrix} 3x_1 + x_2 \\ 5x_1 + 7x_2 \\ x_1 + 3x_2 \end{bmatrix} = \begin{bmatrix} 3 & 1 \\ 5 & 7 \\ 1 & 3 \end{bmatrix} \begin{bmatrix} x_1 \\ x_2 \end{bmatrix} = A\boldsymbol{x}

  1. 一对一证明

    • 矩阵 AA 的列向量 [351]\begin{bmatrix}3\\5\\1\end{bmatrix}[173]\begin{bmatrix}1\\7\\3\end{bmatrix} 线性无关(非倍数关系);
    • 定理13TT 是一对一的。
  2. 映射分析

    • AA3×23 \times 2 矩阵,最多有2个主元位置;
    • 定理2AA 的列无法生成 R3\mathbb{R}^3(需要3个主元) ;
    • TT 不是将 R2\mathbb{R}^2 映射到 R3\mathbb{R}^3 上的。

下图展示更为清晰一些:

alt text

这里很容易知到 T(e1)=[351]T(\mathbf{e_1}) = \begin{bmatrix}3\\5\\1\end{bmatrix}T(e2)=[173]T(\mathbf{e_2}) = \begin{bmatrix}1\\7\\3\end{bmatrix},即二者对应于线性变换后的两个基向量,Span{T(e1),T(e2)}\text{Span}\{T(\mathbf{e_1}), T(\mathbf{e_2})\}只能张成平面,而非三维空间。



第二章 矩阵代数

2.1 矩阵运算

矩阵是线性代数中最基本的数学对象之一,本节将介绍矩阵的基本运算及其性质。

2.1.1 矩阵加法与标量乘法

矩阵加法:若 AABB 是同型矩阵(即具有相同的行数和列数),则它们的和 A+BA+B 是一个同型矩阵,其中每个元素是 AABB 对应元素的和。

(A+B)ij=Aij+Bij(A+B)_{ij} = A_{ij} + B_{ij}

标量乘法:若 cc 是标量,AA 是矩阵,则 cAcA 是将 AA 的每个元素都乘以 cc 得到的矩阵。

(cA)ij=cAij(cA)_{ij} = c \cdot A_{ij}

性质

  1. 交换律:A+B=B+AA + B = B + A
  2. 结合律:(A+B)+C=A+(B+C)(A + B) + C = A + (B + C)
  3. 分配律:c(A+B)=cA+cBc(A + B) = cA + cB(c+d)A=cA+dA(c + d)A = cA + dA

2.1.2 矩阵乘法

矩阵乘法:若 AAm×nm \times n 矩阵,BBn×pn \times p 矩阵,则它们的乘积 ABAB 是一个 m×pm \times p 矩阵,其中:

(AB)ij=k=1nAikBkj(AB)_{ij} = \sum_{k=1}^{n} A_{ik} \cdot B_{kj}

也就是说,ABAB 的第 ii 行第 jj 列元素是 AA 的第 ii 行与 BB 的第 jj 列的点积。

例题:计算下列矩阵乘积:

A=[123456],B=[789101112]A = \begin{bmatrix} 1 & 2 & 3 \\ 4 & 5 & 6 \end{bmatrix}, \quad B = \begin{bmatrix} 7 & 8 \\ 9 & 10 \\ 11 & 12 \end{bmatrix}

AA2×32 \times 3 矩阵,BB3×23 \times 2 矩阵,所以 ABAB2×22 \times 2 矩阵。

AB=[123456][789101112]=[17+29+31118+210+31247+59+61148+510+612]=[5864139154]\begin{aligned} AB &= \begin{bmatrix} 1 & 2 & 3 \\ 4 & 5 & 6 \end{bmatrix} \begin{bmatrix} 7 & 8 \\ 9 & 10 \\ 11 & 12 \end{bmatrix} \\ &= \begin{bmatrix} 1 \cdot 7 + 2 \cdot 9 + 3 \cdot 11 & 1 \cdot 8 + 2 \cdot 10 + 3 \cdot 12 \\ 4 \cdot 7 + 5 \cdot 9 + 6 \cdot 11 & 4 \cdot 8 + 5 \cdot 10 + 6 \cdot 12 \end{bmatrix} \\ &= \begin{bmatrix} 58 & 64 \\ 139 & 154 \end{bmatrix} \end{aligned}

矩阵乘法的性质

  1. 一般情况下,矩阵乘法不满足交换律,即 ABBAAB \neq BA
  2. 结合律:(AB)C=A(BC)(AB)C = A(BC)
  3. 分配律:A(B+C)=AB+ACA(B+C) = AB + AC(A+B)C=AC+BC(A+B)C = AC + BC

关于结合律的解释:

这里我们可以理解为复合映射,具体见下图所示:

这里我们取矩阵 CC 为列向量 x\mathbf{x},且 xRp\mathbf{x} \in \mathbb{R}^p

alt text

关于结合律的简单证明:

AAm×nm \times n 矩阵,BBn×pn \times p 矩阵,xRp\mathbf{x} \in \mathbb{R}^p,用 b1,,bp\mathbf{b_1}, \cdots, \mathbf{b_p} 表示 BB 的各列,而 xx 的元素为 x1,,xpx_1, \cdots, x_p,则:

Bx=x1b1++xpbpB\mathbf{x} = x_1 \mathbf{b_1} + \cdots + x_p \mathbf{b_p}

由乘以 AA 的线性性质,

A(Bx)=A(x1b1++xpbp)=x1Ab1++xpAbpA(B\mathbf{x}) = A(x_1 \mathbf{b_1} + \cdots + x_p \mathbf{b_p}) = x_1 A\mathbf{b_1} + \cdots + x_p A\mathbf{b_p}

向量 A(Bx)A(B\mathbf{x}) 是向量 Ab1,,AbpA\mathbf{b_1}, \cdots, A\mathbf{b_p} 的线性组合,以 xx 的元素为权。若我们把这些向量表示成一个矩阵的各列,就有:

A(Bx)=[Ab1Ab2Abp]xA(B\mathbf{x}) = [A\mathbf{b_1} \quad A\mathbf{b_2} \quad \cdots \quad A\mathbf{b_p}] \mathbf{x}

于是乘以矩阵 [Ab1Ab2Abp][A\mathbf{b_1} \quad A\mathbf{b_2} \quad \cdots \quad A\mathbf{b_p}]x\mathbf{x} 变换为 A(Bx)A(B\mathbf{x}),这便找到了所需要的矩阵。

定义:

AAm×nm \times n 的矩阵,BBn×pn \times p 的矩阵,BB 的各列为 b1,,bp\mathbf{b_1}, \cdots, \mathbf{b_p},则 ABAB 是一个 m×pm \times p 的矩阵,它的各列为 Ab1,,AbpA\mathbf{b_1}, \cdots, A\mathbf{b_p},即

AB=[Ab1Ab2Abp]AB = [A\mathbf{b_1} \quad A\mathbf{b_2} \quad \cdots \quad A\mathbf{b_p}]

计算AB的行列法则:

若乘积 ABAB 有定义,则 ABAB 的第ii行第jj列的元素是AA的第ii行与BB的第jj列的点积,即

(AB)ij=k=1nAikBkj(AB)_{ij} = \sum_{k=1}^{n} A_{ik} \cdot B_{kj}

2.1.3 矩阵乘法的性质

定理14:

AAm×nm \times n 矩阵,BBCC 的维数使下列各式的乘积有定义:

  1. A(BC)=(AB)C(乘法结合律)A(BC) = (AB)C \quad \text{(乘法结合律)}

  2. A(B+C)=AB+AC(乘法左分配律)A(B + C) = AB + AC \quad \text{(乘法左分配律)}

  3. (B+C)A=BA+CA(乘法右分配律)(B + C)A = BA + CA \quad \text{(乘法右分配律)}

  4. r(AB)=(rA)B=A(rB)r 为任意数r(AB) = (rA)B = A(rB),\quad r \text{ 为任意数}

  5. ImA=A=AIn(矩阵乘法的恒等式)I_m A = A = A I_n \quad \text{(矩阵乘法的恒等式)}

2.1.4 矩阵的乘幂

定义:

AA 是一个 n×nn \times n方阵(只有方阵才定义乘幂),则称:

  • A0=InA^0 = I_n,单位矩阵;
  • A1=AA^1 = A
  • A2=AAA^2 = AA,即矩阵与自身相乘;
  • Ak=AAAkA^k = \underbrace{A A \cdots A}_{k \text{个}}

即:

Ak={In,k=0AAk1,k>0A^k = \begin{cases} I_n, & k = 0 \\ A \cdot A^{k-1}, & k > 0 \end{cases}

运算规则:

AAn×nn \times n 矩阵,r,sNr, s \in \mathbb{N},则有:

  • 幂的加法法则

    ArAs=Ar+sA^r A^s = A^{r+s}

  • 幂的乘法法则

    (Ar)s=Ars\left( A^r \right)^s = A^{rs}

  • 与标量的结合律

    (rA)k=rkAk(当 rR)(rA)^k = r^k A^k \quad\text{(当 \( r \in \mathbb{R} \))}

  • 单位矩阵的乘幂恒等式:

    Ink=InI_n^k = I_n

注意事项

  • 非方阵没有幂的定义,即只有方阵才能进行乘幂运算;
  • 一般来说,矩阵幂并不具可交换性,即 AkBk(AB)kA^k B^k \neq (AB)^k
  • 幂运算常用于研究差分方程、幂级数展开、Markov链等问题。

2.1.5 矩阵转置

矩阵 AA转置记为 ATA^T,是将 AA 的行与列互换得到的矩阵。如果 AAm×nm \times n 矩阵,则 ATA^Tn×mn \times m 矩阵,且 (AT)ij=Aji(A^T)_{ij} = A_{ji}

转置的性质

  1. (AT)T=A(A^T)^T = A
  2. (A+B)T=AT+BT(A + B)^T = A^T + B^T
  3. (cA)T=cAT(cA)^T = cA^T,其中 cc 是标量
  4. (AB)T=BTAT(AB)^T = B^T A^T

例题:验证 (AB)T=BTAT(AB)^T = B^T A^T

:设 AAm×nm \times n 矩阵,BBn×pn \times p 矩阵。

((AB)T)ij=(AB)ji=k=1nAjkBki=k=1nBkiAjk=k=1nBikTAkjT=(BTAT)ij\begin{aligned} ((AB)^T)_{ij} &= (AB)_{ji} \\ &= \sum_{k=1}^{n} A_{jk} \cdot B_{ki} \\ &= \sum_{k=1}^{n} B_{ki} \cdot A_{jk} \\ &= \sum_{k=1}^{n} B^T_{ik} \cdot A^T_{kj} \\ &= (B^T A^T)_{ij} \end{aligned}

因此,(AB)T=BTAT(AB)^T = B^T A^T

2.1.6 特殊矩阵

零矩阵:所有元素都为零的矩阵,记为 OO

单位矩阵:主对角线上的元素都为1,其余元素都为0的方阵,记为 InI_nnn 表示阶数)。

In=[100010001]I_n = \begin{bmatrix} 1 & 0 & \cdots & 0 \\ 0 & 1 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & 1 \end{bmatrix}

单位矩阵的重要性质:对任意 n×nn \times n 矩阵 AA,有 AIn=InA=AAI_n = I_n A = A

对角矩阵:非主对角线元素都为0的方阵。

D=[d1000d2000dn]D = \begin{bmatrix} d_1 & 0 & \cdots & 0 \\ 0 & d_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & d_n \end{bmatrix}

对称矩阵:满足 A=ATA = A^T 的方阵,即 Aij=AjiA_{ij} = A_{ji}

反对称矩阵:满足 A=ATA = -A^T 的方阵,即 Aij=AjiA_{ij} = -A_{ji}


2.2 矩阵的逆

2.2.1 可逆矩阵的定义

定义:

AA 是一个 n×nn \times n 的方阵,如果存在另一个 n×nn \times n 的方阵 BB,使得:

AB=BA=InAB = BA = I_n

则称 AA可逆矩阵(或非奇异矩阵),BBAA逆矩阵,记为 A1A^{-1}

一般来说,我们称可逆矩阵为非奇异矩阵不可逆矩阵为奇异矩阵

定理15:

如果方阵 AA 是可逆的,则它的逆矩阵是唯一的。

证明:

假设 BBCC 都是 AA 的逆矩阵,则:

AB=BA=InAC=CA=InAB = BA = I_n \quad \text{且} \quad AC = CA = I_n

于是有:

B=BIn=B(AC)=(BA)C=InC=CB = BI_n = B(AC) = (BA)C = I_nC = C

因此,AA 的逆矩阵是唯一的。

2.2.2 逆矩阵的性质

AABB 都是可逆的 n×nn \times n 矩阵,则:

  1. (A1)1=A(A^{-1})^{-1} = A
  2. (AB)1=B1A1(AB)^{-1} = B^{-1}A^{-1}
  3. (AT)1=(A1)T(A^T)^{-1} = (A^{-1})^T
  4. 对任意非零标量 cc(cA)1=1cA1(cA)^{-1} = \frac{1}{c}A^{-1}
  5. 对任意正整数 kk(Ak)1=(A1)k(A^k)^{-1} = (A^{-1})^k

证明:

对于性质2,我们需要验证 (AB)(B1A1)=(B1A1)(AB)=In(AB)(B^{-1}A^{-1}) = (B^{-1}A^{-1})(AB) = I_n

(AB)(B1A1)=A(BB1)A1=AInA1=AA1=In\begin{aligned} (AB)(B^{-1}A^{-1}) &= A(BB^{-1})A^{-1} \\ &= AI_nA^{-1} \\ &= AA^{-1} \\ &= I_n \end{aligned}

同理可证 (B1A1)(AB)=In(B^{-1}A^{-1})(AB) = I_n

2.2.3 矩阵可逆的充要条件

定理16:

n×nn \times n 矩阵 AA 可逆的充要条件是 AA 行等价于 InI_n

进一步说明为:把 AA 经过一系列初等行变换化简为 InI_n 的同时,可以将 InI_n 化简为 A1A^{-1}

定理17:

n×nn \times n 矩阵 AA 可逆的充要条件是 AA 的简化行阶梯形为 InI_n

定理18(可逆矩阵定理):

AAn×nn \times n 矩阵,则以下命题等价:

  1. AA 是可逆矩阵。
  2. AA 行等价于 InI_n
  3. AAnn 个主元位置。
  4. 齐次方程 Ax=0A\mathbf{x} = \mathbf{0} 仅有平凡解。
  5. 方程 Ax=bA\mathbf{x} = \mathbf{b} 对任意 bRn\mathbf{b} \in \mathbb{R}^n 有唯一解。
  6. AA 的列向量线性无关。
  7. AA 的列向量构成 Rn\mathbb{R}^n 的一组基。
  8. AA 的行向量线性无关。
  9. AA 的行向量构成 Rn\mathbb{R}^n 的一组基。
  10. det(A)0\text{det}(A) \neq 0(行列式不为零)。
  11. 0 不是 AA 的特征值。
  12. rank(A)=n\text{rank}(A) = n(满秩)。

2.2.4 求逆矩阵的方法

方法一:初等行变换法

将矩阵 AA 与单位矩阵 InI_n 并排写成增广矩阵 [AIn][A|I_n],然后通过初等行变换将左侧的 AA 化为 InI_n,此时右侧得到的矩阵就是 A1A^{-1}

这里我们做出如下解释:

原始矩阵方程为 Ax=bA\mathbf{x} = \mathbf{b},我们同时在等号两侧左乘 A1A^{-1} 后,我们可以得到 x=A1b\mathbf{x} = A^{-1}\mathbf{b},这两个矩阵方程是等价的,因此对应的增广矩阵也是等价的

例题:

求矩阵 A=[1237]A = \begin{bmatrix} 1 & 2 \\ 3 & 7 \end{bmatrix} 的逆矩阵。

解:

构造增广矩阵 [AI2][A|I_2]

[12103701]\begin{bmatrix} 1 & 2 & 1 & 0 \\ 3 & 7 & 0 & 1 \end{bmatrix}

通过初等行变换将左侧变为单位矩阵:

  1. R23R1R2R_2 - 3R_1 \rightarrow R_2

[12100131]\begin{bmatrix} 1 & 2 & 1 & 0 \\ 0 & 1 & -3 & 1 \end{bmatrix}

  1. R12R2R1R_1 - 2R_2 \rightarrow R_1

[10720131]\begin{bmatrix} 1 & 0 & 7 & -2 \\ 0 & 1 & -3 & 1 \end{bmatrix}

因此,A1=[7231]A^{-1} = \begin{bmatrix} 7 & -2 \\ -3 & 1 \end{bmatrix}

方法二:伴随矩阵法

对于 n×nn \times n 矩阵 AA,其逆矩阵可以表示为:

A1=1det(A)adj(A)A^{-1} = \frac{1}{\text{det}(A)} \text{adj}(A)

其中,det(A)\text{det}(A)AA 的行列式,adj(A)\text{adj}(A)AA 的伴随矩阵。

伴随矩阵我们在后续章节中再介绍。

2.2.5 逆矩阵的应用

  1. 求解线性方程组:对于方程组 Ax=bA\mathbf{x} = \mathbf{b},若 AA 可逆,则解为 x=A1b\mathbf{x} = A^{-1}\mathbf{b}

  2. 线性变换的逆变换:若线性变换 TT 由矩阵 AA 表示,且 AA 可逆,则 TT 的逆变换由 A1A^{-1} 表示。

  3. 矩阵方程求解:对于矩阵方程 AX=BAX = BXA=BXA = B,若 AA 可逆,则解分别为 X=A1BX = A^{-1}BX=BA1X = BA^{-1}


2.3 可逆矩阵的特征

2.3.1 可逆矩阵的特征

AAn×nn \times n 矩阵,则下列 ​12 个命题等价​(对任意给定的 AA,它们同时为真或同时为假):

标记 命题描述
a. AA 是可逆矩阵
b. AA 行等价于 n×nn \times n 单位矩阵
c. AAnn 个主元位置
d. 方程 Ax=0A\mathbf{x} = \mathbf{0} 仅有平凡解
e. AA 的各列线性无关
f. 线性变换 xAx\mathbf{x} \mapsto A\mathbf{x} 是一对一的
g. Rn\mathbb{R}^n 中任意 b\mathbf{b},方程 Ax=bA\mathbf{x} = \mathbf{b} 至少有一个解
h. AA 的各列生成 Rn\mathbb{R}^n
i. 线性变换 xAx\mathbf{x} \mapsto A\mathbf{x}Rn\mathbb{R}^n 映上到 Rn\mathbb{R}^n
j. 存在 n×nn \times n 矩阵 CC 使得 CA=ICA = I
k. 存在 n×nn \times n 矩阵 DD 使得 AD=IAD = I
1. ATA^{T} 是可逆矩阵

核心等价关系图解

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
graph LR
a[a.可逆矩阵] --> b[b.行等价单位矩阵]
b --> c[c.有n主元位置]
c --> d[d.仅有零解]
d --> e[e.列线性无关]
e --> f[f.一对一变换]
f --> g[g.解存在性]
g --> h[h.列生成Rⁿ]
h --> i[i.映上变换]
i --> j[j.左逆存在]
j --> k[k.右逆存在]
k --> 1[1.转置可逆]
1 --> a
classDef logic stroke:#FF6B6B,stroke-width:2px;
class a,b,c,d,e,f,g,h,i,j,k,1 logic;

2.3.2 可逆线性变换

线性变换 T:RnRnT: \mathbb{R}^n \to \mathbb{R}^n 称为可逆的,若存在函数 S:RnRnS: \mathbb{R}^n \to \mathbb{R}^n 使得:

  1. 对所有 xRnx \in \mathbb{R}^n,有

    S(T(x))=x(1)S(T(x)) = x \tag{1}

  2. 对所有 xRnx \in \mathbb{R}^n,有

    T(S(x))=x(2)T(S(x)) = x \tag{2}

下列定理说明:若这样的 SS 存在,则它是唯一的且必是线性变换。我们称 SSTT 的逆,记作 T1T^{-1}

定理19

T:RnRnT: \mathbb{R}^n \to \mathbb{R}^n 为线性变换,AATT 的标准矩阵。则:

TT 可逆 当且仅当 AA 是可逆矩阵。
此时由 S(x)=A1xS(x) = A^{-1}x 定义的线性变换 SS,是满足 (1) 和 (2) 的唯一函数

证明

(⇒)设 TT 是可逆的:

由条件 (2) 可知,TT 是从 Rn\mathbb{R}^n 映到 Rn\mathbb{R}^n 的映射。

因为 bRn\mathbf{b} \in \mathbb{R}^n,设 x=S(b)\mathbf{x} = S(\mathbf{b}),则:

T(x)=T(S(b))=bT(\mathbf{x}) = T(S(\mathbf{b})) = \mathbf{b}

所以每个 b\mathbf{b} 属于 TT 的值域。

因此可知矩阵 AA 是可逆的。

(⇐)反之,若 AA 是可逆的:

令:

S(x)=A1xS(x) = A^{-1}x

显然 SS 是线性变换,且满足条件 (1) 和 (2)。例如:

S(T(x))=S(Ax)=A1(Ax)=xS(T(x)) = S(Ax) = A^{-1}(Ax) = x

于是 TT 是可逆的。


2.4 分块矩阵

2.4.1 分块矩阵的概念

分块矩阵(或称块矩阵分部矩阵)是将一个大矩阵按照行和列划分成若干个小矩阵(块)的表示方法。这种表示方法在处理大型矩阵时非常有用,可以简化计算并揭示矩阵的结构特性。

例如,矩阵 AA 可以按如下方式划分为四个块:

A=[A11A12A21A22]A = \begin{bmatrix} A_{11} & A_{12} \\ A_{21} & A_{22} \end{bmatrix}

其中 A11,A12,A21,A22A_{11}, A_{12}, A_{21}, A_{22} 是适当维度的子矩阵。

2.4.2 分块矩阵的运算

1. 分块矩阵的加法

如果矩阵 AABB 具有相同的分块方式,则它们的和是对应块的和:

A+B=[A11A12A21A22]+[B11B12B21B22]=[A11+B11A12+B12A21+B21A22+B22]A + B = \begin{bmatrix} A_{11} & A_{12} \\ A_{21} & A_{22} \end{bmatrix} + \begin{bmatrix} B_{11} & B_{12} \\ B_{21} & B_{22} \end{bmatrix} = \begin{bmatrix} A_{11}+B_{11} & A_{12}+B_{12} \\ A_{21}+B_{21} & A_{22}+B_{22} \end{bmatrix}

2. 分块矩阵的乘法

AABB 是两个分块矩阵:

A=[A11A12A21A22],B=[B11B12B21B22]A = \begin{bmatrix} A_{11} & A_{12} \\ A_{21} & A_{22} \end{bmatrix}, \quad B = \begin{bmatrix} B_{11} & B_{12} \\ B_{21} & B_{22} \end{bmatrix}

如果各个块的维度使得乘法有定义,则它们的乘积为:

AB=[A11B11+A12B21A11B12+A12B22A21B11+A22B21A21B12+A22B22]AB = \begin{bmatrix} A_{11}B_{11} + A_{12}B_{21} & A_{11}B_{12} + A_{12}B_{22} \\ A_{21}B_{11} + A_{22}B_{21} & A_{21}B_{12} + A_{22}B_{22} \end{bmatrix}

这种乘法规则类似于普通矩阵乘法,但操作对象是矩阵块而非单个元素。

注意:为使分块矩阵乘法有定义,必须保证相邻块的维度匹配,即 AA列分块方式必须与 BB行分块方式一致。

定理20

这里我们令符号 colk(A)col_k(A) 表示 AA 的第 kk 列块,rowk(B)row_k(B) 表示 BB 的第 kk 行块。

AAm×nm \times n 矩阵,BBn×pn \times p 矩阵,则:

AB=[col1(A)col2(A)coln(A)][row1(B)row2(B)rown(B)](1)AB = \begin{bmatrix} \text{col}_1(A) & \text{col}_2(A) & \cdots & \text{col}_n(A) \end{bmatrix} \begin{bmatrix} \text{row}_1(B) \\ \text{row}_2(B) \\ \vdots \\ \text{row}_n(B) \end{bmatrix} \tag{1}

展开为:

AB=col1(A)row1(B)++coln(A)rown(B)AB = \text{col}_1(A)\text{row}_1(B) + \cdots + \text{col}_n(A)\text{row}_n(B)

证明

对每个行指标 ii 和列指标 jj,乘积 colk(A)rowk(B)\text{col}_k(A)\text{row}_k(B)(i,j)(i, j) 元素是:

  • colk(A)\text{col}_k(A) 中的元素 aika_{ik}
  • rowk(B)\text{row}_k(B) 中的元素 bkjb_{kj} 的乘积

因此,在等式 (1) 的和中,(i,j)(i,j) 元素为:

ai1b1j+ai2b2j++ainbnj{(k=1)(k=2)(k=n)a_{i1}b_{1j} + a_{i2}b_{2j} + \cdots + a_{in}b_{nj} \quad \begin{cases} (k=1) \\ (k=2) \\ \vdots \\ (k=n) \end{cases}

根据行列法则,该和正好是矩阵 ABAB 的第 (i,j)(i,j) 元素。

3. 分块矩阵的转置

分块矩阵的转置是将每个块转置,并且交换块的位置:

[A11A12A21A22]T=[A11TA21TA12TA22T]\begin{bmatrix} A_{11} & A_{12} \\ A_{21} & A_{22} \end{bmatrix}^T = \begin{bmatrix} A_{11}^T & A_{21}^T \\ A_{12}^T & A_{22}^T \end{bmatrix}

2.4.3 特殊的分块矩阵

1. 分块对角矩阵

当非对角块都为零矩阵时,称为分块对角矩阵:

[A11OOA22]\begin{bmatrix} A_{11} & O \\ O & A_{22} \end{bmatrix}

分块对角矩阵的行列式等于对角块行列式的乘积:

det[A11OOA22]=det(A11)det(A22)\det\begin{bmatrix} A_{11} & O \\ O & A_{22} \end{bmatrix} = \det(A_{11}) \cdot \det(A_{22})

2. 分块三角矩阵

分块上三角矩阵:

[A11A12OA22]\begin{bmatrix} A_{11} & A_{12} \\ O & A_{22} \end{bmatrix}

分块下三角矩阵:

[A11OA21A22]\begin{bmatrix} A_{11} & O \\ A_{21} & A_{22} \end{bmatrix}

2.4.4 分块矩阵的逆

对于 2×22 \times 2 分块矩阵,如果满足特定条件,其逆矩阵可以通过以下方式计算:

1. 当 A22A_{22} 可逆时

[A11A12A21A22]1=[(A11A12A221A21)1(A11A12A221A21)1A12A221A221A21(A11A12A221A21)1A221+A221A21(A11A12A221A21)1A12A221]\begin{bmatrix} A_{11} & A_{12} \\ A_{21} & A_{22} \end{bmatrix}^{-1} = \begin{bmatrix} (A_{11} - A_{12}A_{22}^{-1}A_{21})^{-1} & -(A_{11} - A_{12}A_{22}^{-1}A_{21})^{-1}A_{12}A_{22}^{-1} \\ -A_{22}^{-1}A_{21}(A_{11} - A_{12}A_{22}^{-1}A_{21})^{-1} & A_{22}^{-1} + A_{22}^{-1}A_{21}(A_{11} - A_{12}A_{22}^{-1}A_{21})^{-1}A_{12}A_{22}^{-1} \end{bmatrix}

其中 S=A11A12A221A21S = A_{11} - A_{12}A_{22}^{-1}A_{21} 称为 A22A_{22}舒尔补

2. 当 A11A_{11} 可逆时

[A11A12A21A22]1=[A111+A111A12(A22A21A111A12)1A21A111A111A12(A22A21A111A12)1(A22A21A111A12)1A21A111(A22A21A111A12)1]\begin{bmatrix} A_{11} & A_{12} \\ A_{21} & A_{22} \end{bmatrix}^{-1} = \begin{bmatrix} A_{11}^{-1} + A_{11}^{-1}A_{12}(A_{22} - A_{21}A_{11}^{-1}A_{12})^{-1}A_{21}A_{11}^{-1} & -A_{11}^{-1}A_{12}(A_{22} - A_{21}A_{11}^{-1}A_{12})^{-1} \\ -(A_{22} - A_{21}A_{11}^{-1}A_{12})^{-1}A_{21}A_{11}^{-1} & (A_{22} - A_{21}A_{11}^{-1}A_{12})^{-1} \end{bmatrix}

其中 T=A22A21A111A12T = A_{22} - A_{21}A_{11}^{-1}A_{12} 称为 A11A_{11}舒尔补

2.4.5 分块矩阵的应用

  1. 大型线性系统求解:将大型线性方程组分解为较小的子系统求解。

  2. 矩阵分解:许多矩阵分解方法(如LU分解、QR分解)可以用分块矩阵形式表示。

  3. 控制理论:在状态空间表示中,系统矩阵常常以分块形式出现。

  4. 图像处理:大型图像可以分块处理,提高计算效率。

例题

计算下列分块矩阵的乘积:

[I2OAI3][I2BOI3]\begin{bmatrix} I_2 & O \\ A & I_3 \end{bmatrix} \begin{bmatrix} I_2 & B \\ O & I_3 \end{bmatrix}

其中 AA3×23 \times 2 矩阵,BB2×32 \times 3 矩阵,I2I_2I3I_3 分别是 2×22 \times 23×33 \times 3 单位矩阵,OO 表示适当维度的零矩阵。

根据分块矩阵乘法规则:

[I2OAI3][I2BOI3]=[I2I2+OOI2B+OI3AI2+I3OAB+I3I3]=[I2BAAB+I3]\begin{bmatrix} I_2 & O \\ A & I_3 \end{bmatrix} \begin{bmatrix} I_2 & B \\ O & I_3 \end{bmatrix} = \begin{bmatrix} I_2 \cdot I_2 + O \cdot O & I_2 \cdot B + O \cdot I_3 \\ A \cdot I_2 + I_3 \cdot O & A \cdot B + I_3 \cdot I_3 \end{bmatrix} = \begin{bmatrix} I_2 & B \\ A & AB + I_3 \end{bmatrix}