线性代数在各大理工科,乃至经济金融领域的使用之广泛,毋庸置疑。 一直以来,我虽也知道线性代数的重要,但从内心上其实一直是犯怵的(尤其是学习论文、算法中,基本只要看到对方把算法向量化之后就蒙圈了),当年在学校学习过程中很多也是靠着死记硬背过来的,对它的直观意义一直都没能有很好的理解。
最近,这么一本书进入了我的视线:《线性代数及其应用》,听书名感觉平平,但只翻了几页就感觉十分过瘾,仿佛打通了任督二脉。以往很多死记硬背的知识点在这本书的解释下,变成了可以直观推导出来的结果。这本书不仅对线性代数的基本概念阐述地很直观形象,而且还有许多现实生活中的应用,特别是经济、物理、计算机领域,真正让人领略到线性代数作为现代数学的魅力。
我特将自己的读书总结和体会记录于此,也是希望借此加深自己的理解。
注意,这个系列假设你已经有了线性代数基础,像是行变换、将矩阵转换为行阶梯形式这种基本技巧已经掌握。本文不再赘述具体操作步骤,主要关注于概念的直观理解。
一、线性方程组
线性代数,最基本的问题,就是解线性方程组了。线性方程组就是一组形如 a1x1+a2x2+⋯+anxn=b 的方程。一个线性方程组中的变量是相同的,如果第一个方程是关于 x1⋯xn 的,那么其他的也都应该如此(这些变量不一定都出现,因为系数可以有 0)。
1.1 线性方程组的矩阵形式
方程组
⎩⎨⎧x1−4x1−+2x22x25x2+−+x38x39x3=0=8=−9
可以通过増广矩阵形式描述:
10−4−2251−8908−9
增广矩阵去掉最后一列,就是该方程组的系数矩阵。
矩阵形式只是线性方程组的一种表示形式。今后的很多关于线性方程组的计算,都将在矩阵形式上进行操作,然而你也需要知道,在这些操作进行的同时,线性方程组也在进行类似的变换。比如,将增广矩阵的第一、二行对换,那么同时,它所代表的线性方程组中,第一、二个方程也进行了对调。
1.2 线性方程组的解
解一个线性方程组,就是通过对其矩阵形式行变换(三种方式:交换方程的先后顺序,一个方程左右同乘以某数,和两个方程相加) 转换为行阶梯形式。比如
10−4−2251−8908−9转化为行阶梯形式:100−2101−41043和最简形式:10001000129163
上面最简形式的矩阵对应的线性方程组是
⎩⎨⎧x1x2x3=29=16=3
这个线性方程组和一开始的方程组是等价的,只是处于不同的状态,它们的解也是相同的,而显然行最简形式的方程组最容易解,所以我们一般都将线性方程组的増广矩阵转化为行最简形式继而求解。
1.3 解的存在性和唯一性
还记得线性代数时经常讨论的「无解「「唯一解」「无穷多解」吧?
首先来看刚才的方程组,经过行变换后,方程组的解已经很显然了:x1=29,x2=16,x3=3。这个方程组的解就只有一个,是唯一解。
1.3.1 无解
我们再来看一个方程组:
⎩⎨⎧2x15x1−−x23x28x2−++4x32x37x3=8=1=1
它的增广矩阵
0251−3−8−427811 行变换得到:200−3102−40185/2
变换后的矩阵所对应的方程组为
⎩⎨⎧2x1−3x2x2+−2x34x30=1=8=5/2
显然,第三个方程 0=5/2 是无解的。对比这个方程组和它对应的增广矩阵,我们可以发现,当增广矩阵的行阶梯形式存在 [0 ⋯ 0 b] 形式的行时,方程组无解。
1.3.2 有解
当增广矩阵变换为行阶梯形式后,不存在 [0 ⋯ 0 b] 形式的行,则说明方程有解。我们接下来讨论下它的解具体会是怎么样的。
假设现在有这样一个已经化为行最简形式的增广矩阵:
100010−510140
这个矩阵有 4 列,故而有 3 个变量。相对应的方程组为:
⎩⎨⎧x1x2−+5x3x30=1=4=0
观察这个方程组,x1 和 x2 只存在于一个方程中(对应行最简形式中的主元位置),x3 存在于两个方程中。那么我们可以通过 x3 来表示 x1 和 x2:
⎩⎨⎧x1x2x3=5x3+1=x3+4是自由变量
上面列出来的实际上就是这个方程组的解集了。x1 和 x2 被称为「基本变量」;x3被称为「自由变量」,因为它在解集里不受任何约束,而基本变量需要自由变量来表示; 也就是说,自由变量确定了一个值,基本变量也就随之确定了一个值。上面这个解集形式也被称为方程组的「通解」,因为它给出了方程组所有解的显式表示。
需要注意的是,我们需要先将增广矩阵变换为行最简形式,才能知道谁是自由变量,谁是基本变量。
因为自由变量能取任意值,所以,存在自由变量的线性方程组有无穷多解,而没有自由变量的线性方程组则只有一个唯一解(就像本文第一个方程组那样)。
总结一下:
- 当增广矩阵的行阶梯形式(当然行最简形式也可以)存在[0 ⋯ 0 b] 形式时,方程组无解;
- 当增广矩阵的行最简形式不存在自由变量时,方程组有唯一解;
- 当增广矩阵的行最简形式存在自由变量时,方程组有无穷多解;
二、向量方程
n 维空间中的点可用 n 维向量表示。
向量之间可以线性组合:
y=c1v1+⋯+cpvp
那么,假如有三个向量:a1=[1,−2,−5]T,a2=[2,5,6]T,b=[7,4,−3]T,想要知道 b 是否能通过 a1 和 a2 线性表示,实际上就是求线性方程 x1a1+x2a2=b 是否有解的问题。
把这个方程展开来看,就是:
x11−2−5+x2256=74−3
等同于
x1−2x1−5x1+2x25x26x2=74−3
和
x1+2x2−2x1+5x2−5x1+6x2=74−3
所以这个问题其实和一个线性方程组是等价的,这个线性方程组对应的増广矩阵就是([a1,a2,b]):
1−2−525−674−3
化简为行最简形式就是:
100010320
可以看出,这个线性方程组的解为 x1=3 和 x2=2。继而我们就知道了 b 和 a1, a2 的关系:
31−2−5+2256=74−3
我们反过来回顾这一过程,可以发现,之前我们线性方程组的的增广矩阵表示形式,其实也可以看做是列向量组成的形式,在这个例子中,增广矩阵可以表示为
[a1,a2,b]
。把增广矩阵按列拆开看,我们就可以得到线性方程组的向量方程表示形式。
向量方程是线性方程组另一种重要的表现形式,它能帮助我们将矩阵、线性方程组的抽象概念同几何的直观联系起来。
在几何中,n 个向量 v1,v2,⋯,vp 的所有线性组合 c1v1+c2v2+⋯+cpvp 成为一个空间,称作由 v1,v2,⋯,vp 张成的 Rn 的子空间,记作 Span{v1,⋯,vp}。
一个向量张成的空间是一根直线,两个向量张成的空间是一个平面。
三、矩阵方程
向量的线性组合可以看作向量与矩阵的乘积,比如一个 m×n 的矩阵 A,各列为 a1,⋯,an,而 x 为 n 维向量,则有:
Ax=[a1 a2 ⋯an]x1⋮xn=x1a1+x2a2+⋯+xnan
这种形如 Ax=b 的形式,就称为矩阵方程。
由矩阵方程的定义,我们可以得出:**方程Ax=b有解当且仅当b为A中列的线性组合。**又因为我们之前提到,这些列向量的所有线性组合构成了Span{a1,⋯,an},向量 b 是否存在于这个空间,就等价于 Ax=b 有解。
下面我们来讨论下任意 b∈Rm 的情况。
设
A=1−4−332−24−6−7,b=b1b2b3
求方程 Ax=b 是否对 b1,b2,b3 的所有取值都有解?
我们首先对增广矩阵化简:
1−4−332−24−6−7b1b2b3∼10031404100b1b2+4b1b1−21b2+b3
可以看出,当b 取某些值时,b1−21b2+b3 不等于0,于是就会有无解的情况。只有当
b1−21b2+b3=0
时方程才有解。注意,这个式子在几何中表示三维中的一个平面, 结合Ax=b,这个平面就是A 中列向量线性组合构成的集合。
本来 b 是三维的向量,如果没有限制的话它可以表示整个三维空间,然而,在这个空间中,一大部分都不满足使 Ax=b 有解。这仅剩的一个平面就是 A 的列向量所能张成的全部空间。 这些三维列向量最终张成了一个二维平面。
观察行最简形式矩阵,可以知道,之所以 b 的一些取值造成矩阵方程无解,是因为系数矩阵 A 中最后一行没有主元,在行最简形式中变成了形如 [0 0 0 b] 的行。如果系数矩阵 A 中每一行都有主元的话,那么就不会出现无解的情况。
反过来看,当 n 个 m 维列向量能张成 Rm 时,就说明对任意 b∈Rm,方程 Ax=b 都有解,也就是说,Rm 空间中的任意向量,都可以由 A 的列线性表示。
总结一下,就是以下四点相互等价。
- 对任意 b∈Rm,方程 Ax=b 都有解。
- 任意 b∈Rm 都是A 中列的一个线性组合。
- A 的列张成 Rm。
- A 中每一行都有主元位置。
四、三种等价形式
矩阵方程
Ax=b
和向量方程
x1a1+x2a2+⋯+xnan=b
以及下列增广矩阵对应的线性方程组具有相同的解集
[a1 a2 ⋯ an b]
矩阵方程、向量方程和线性方程组是三种不同但却相互等价的形式。在现实生活中构造一个数学模型时,我们可以在任何情况下自由选择其中任何一种最自然、最便利的陈述形式。
以上三种形式就是我们在解线性方程组时的三个工具,结合具体问题,我们可以通过不同角度观察问题,进而求解。 另外,这三种形式的求解,都是对增广矩阵进行行化简,因此,増广矩阵的行变换是一切的基础。
本文中所有文字、图片版权均属本人所有,如需转载请注明来源。