:label:sec_linear-algebra
在介绍完如何存储和操作数据后,接下来将简要地回顾一下部分基本线性代数内容。
这些内容有助于读者了解和实现本书中介绍的大多数模型。
本节将介绍线性代数中的基本数学对象、算术和运算,并用数学符号和相应的代码实现来表示它们。
标量
如果你曾经在餐厅支付餐费,那么应该已经知道一些基本的线性代数,比如在数字间相加或相乘。 例如,北京的温度为52∘𝐹(华氏度,除摄氏度外的另一种温度计量单位)。
严格来说,仅包含一个数值被称为标量(scalar)。
如果要将此华氏度值转换为更常用的摄氏度,则可以计算表达式𝑐=59(𝑓−32),并将𝑓赋为52。在此等式中,每一项(5、9和32)都是标量值。 符号𝑐和𝑓称为变量(variable),它们表示未知的标量值。
本书采用了数学表示法,其中标量变量由普通小写字母表示(例如,𝑥、𝑦和𝑧)。
本书用𝑅表示所有(连续)实数标量的空间,之后将严格定义空间(space)是什么, 但现在只要记住表达式𝑥∈𝑅是表示𝑥是一个实值标量的正式形式。 符号∈称为“属于”,它表示“是集合中的成员”。 例如𝑥,𝑦∈{0,1}可以用来表明𝑥和𝑦是值只能为0或1的数字。
标量由只有一个元素的张量表示。
下面的代码将实例化两个标量,并执行一些熟悉的算术运算,即加法、乘法、除法和指数。
向量
[向量可以被视为标量值组成的列表]。
这些标量值被称为向量的元素(element)或分量(component)。
当向量表示数据集中的样本时,它们的值具有一定的现实意义。 例如,如果我们正在训练一个模型来预测贷款违约风险,可能会将每个申请人与一个向量相关联, 其分量与其收入、工作年限、过往违约次数和其他因素相对应。 如果我们正在研究医院患者可能面临的心脏病发作风险,可能会用一个向量来表示每个患者, 其分量为最近的生命体征、胆固醇水平、每天运动时间等。
在数学表示法中,向量通常记为粗体、小写的符号 (例如,𝑥、𝑦和𝑧))。
人们通过一维张量表示向量。一般来说,张量可以具有任意长度,取决于机器的内存限制。
我们可以使用下标来引用向量的任一元素,例如可以通过𝑥𝑖来引用第𝑖个元素。
注意,元素𝑥𝑖是一个标量,所以我们在引用它时不会加粗。 大量文献认为列向量是向量的默认方向,在本书中也是如此。 在数学中,向量𝑥可以写为:
:eqlabel:eq_vec_def
其中𝑥1,…,𝑥𝑛是向量的元素。
在代码中,我们通过张量的索引来访问任一元素。
长度、维度和形状
向量只是一个数字数组,就像每个数组都有一个长度一样,每个向量也是如此。
在数学表示法中,如果我们想说一个向量𝑥由𝑛个实值标量组成, 可以将其表示为𝑥∈𝑅𝑛。
向量的长度通常称为向量的维度(dimension)。
与普通的Python数组一样,我们可以通过调用Python的内置len()
函数来[访问张量的长度]。
当用张量表示一个向量(只有一个轴)时,我们也可以通过.shape
属性访问向量的长度。
形状(shape)是一个元素组,列出了张量沿每个轴的长度(维数)。
对于只有一个轴的张量,形状只有一个元素。
请注意,维度(dimension)这个词在不同上下文时往往会有不同的含义,这经常会使人感到困惑。
为了清楚起见,我们在此明确一下:
向量或轴的维度 = 向量或轴的长度 = 向量或轴的元素数量
然而,张量的维度用来表示张量具有的轴数。
在这个意义上,张量的某个轴的维数就是这个轴的长度。
矩阵
正如向量将标量从零阶推广到一阶,矩阵将向量从一阶推广到二阶。
矩阵,我们通常用粗体、大写字母来表示,在代码中表示为具有两个轴的张量。
数学表示法使用𝐴∈𝑅𝑚×𝑛 来表示矩阵𝐴,其由𝑚行和𝑛列的实值标量组成。 我们可以将任意矩阵𝐴∈𝑅𝑚×𝑛视为一个表格, 其中每个元素𝑎𝑖𝑗属于第𝑖行第𝑗列:
:eqlabel:eq_matrix_def
对于任意𝐴∈𝑅𝑚×𝑛, 𝐴的形状是(𝑚,𝑛)或𝑚×𝑛。 当矩阵具有相同数量的行和列时,其形状将变为正方形; 因此,它被称为方阵(square matrix)。
当调用函数来实例化张量时, 我们可以[通过指定两个分量𝑚和𝑛来创建一个形状为𝑚×𝑛的矩阵]。
我们可以通过行索引(𝑖)和列索引(𝑗)来访问矩阵中的标量元素𝑎𝑖𝑗, 例如[𝐴]𝑖𝑗。
如果没有给出矩阵𝐴的标量元素,如在 :eqref:eq_matrix_def
那样, 我们可以简单地使用矩阵𝐴的小写字母索引下标𝑎𝑖𝑗 来引用[𝐴]𝑖𝑗。 为了表示起来简单,只有在必要时才会将逗号插入到单独的索引中, 例如𝑎2,3𝑗和[𝐴]2𝑖−1,3。
当我们交换矩阵的行和列时,结果称为矩阵的转置(transpose)。
通常用𝑎⊤来表示矩阵的转置,如果𝐵=𝐴⊤, 则对于任意𝑖和𝑗,都有𝑏𝑖𝑗=𝑎𝑗𝑖。
因此,在 :eqref:eq_matrix_def
中的转置是一个形状为𝑛×𝑚的矩阵:
现在在代码中访问(矩阵的转置)。
作为方阵的一种特殊类型,[对称矩阵(symmetric matrix)𝐴等于其转置:𝐴=𝐴⊤]。
这里定义一个对称矩阵𝐵:
现在我们将B
与它的转置进行比较。
矩阵是有用的数据结构:它们允许我们组织具有不同模式的数据。
例如,我们矩阵中的行可能对应于不同的房屋(数据样本),而列可能对应于不同的属性。 曾经使用过电子表格软件或已阅读过 :numref:sec_pandas
的人,应该对此很熟悉。 因此,尽管单个向量的默认方向是列向量,但在表示表格数据集的矩阵中, 将每个数据样本作为矩阵中的行向量更为常见。 后面的章节将讲到这点,这种约定将支持常见的深度学习实践。 例如,沿着张量的最外轴,我们可以访问或遍历小批量的数据样本。
张量
[就像向量是标量的推广,矩阵是向量的推广一样,我们可以构建具有更多轴的数据结构]。
张量(本小节中的“张量”指代数对象)是描述具有任意数量轴的𝑛维数组的通用方法。
向量是一阶张量,矩阵是二阶张量。
张量用特殊字体的大写字母表示, 它们的索引机制(例如𝑥𝑖𝑗𝑘和[𝑋]1,2𝑖−1,3)与矩阵类似。
当我们开始处理图像时,张量将变得更加重要,图像以𝑛维数组形式出现, 其中3个轴对应于高度、宽度,以及一个通道(channel)轴, 用于表示颜色通道(红色、绿色和蓝色)。现在先将高阶张量暂放一边,而是专注学习其基础知识。
张量算法的基本性质
标量、向量、矩阵和任意数量轴的张量(本小节中的“张量”指代数对象)有一些实用的属性。 例如,从按元素操作的定义中可以注意到,任何按元素的一元运算都不会改变其操作数的形状。
[给定具有相同形状的任意两个张量,任何按元素二元运算的结果都将是相同形状的张量]。
例如,将两个相同形状的矩阵相加,会在这两个矩阵上执行元素加法。
[两个矩阵的按元素乘法称为Hadamard积(Hadamard product)(数学符号⊙)]。
对于矩阵𝐵∈𝑅𝑚×𝑛, 其中第𝑖行和第𝑗列的元素是𝑏𝑖𝑗。 矩阵𝐴(在 :eqref:eq_matrix_def
中定义)和𝐵的Hadamard积为:
将张量乘以或加上一个标量不会改变张量的形状,其中张量的每个元素都将与标量相加或相乘。
降维
:label:subseq_lin-alg-reduction
我们可以对任意张量进行的一个有用的操作是[计算其元素的和]。
数学表示法使用∑符号表示求和。
为了表示长度为𝑑的向量中元素的总和,可以记为
在代码中可以调用计算求和的函数:
我们可以(表示任意形状张量的元素和)。 例如,矩阵𝐴中元素的和可以记为
默认情况下,调用求和函数会沿所有的轴降低张量的维度,使它变为一个标量。
我们还可以[指定张量沿哪一个轴来通过求和降低维度]。
以矩阵为例,为了通过求和所有行的元素来降维(轴0),可以在调用函数时指定axis=0
。 由于输入矩阵沿0轴降维以生成输出向量,因此输入轴0的维数在输出形状中消失。
矩阵A:
指定axis=1
将通过汇总所有列的元素降维(轴1)。因此,输入轴1的维数在输出形状中消失。
沿着行和列对矩阵求和,等价于对矩阵的所有元素进行求和。
[一个与求和相关的量是平均值(mean或average)]。
我们通过将总和除以元素总数来计算平均值。
在代码中,我们可以调用函数来计算任意形状张量的平均值。
同样,计算平均值的函数也可以沿指定轴降低张量的维度。
非降维求和
:label:subseq_lin-alg-non-reduction
但是,有时在调用函数来[计算总和或均值时保持轴数不变]会很有用。
例如,由于sum_A
在对每行进行求和后仍保持两个轴,我们可以通过广播将A
除以sum_A
。
如果我们想沿[某个轴计算A
元素的累积总和], 比如axis=0
(按行计算),可以调用cumsum
函数。 此函数不会沿任何轴降低输入张量的维度。
点积(Dot Product)
我们已经学习了按元素操作、求和及平均值。 另一个最基本的操作之一是点积。
给定两个向量𝑥,𝑦∈𝑅𝑑, 它们的点积(dot product)𝑥⊤𝑦 (或〈𝑥,𝑦〉) 是相同位置的按元素乘积的和:
[点积是相同位置的按元素乘积的和]
注意,(我们可以通过执行按元素乘法,然后进行求和来表示两个向量的点积):
点积在很多场合都很有用。 例如,给定一组由向量𝑥∈𝑅𝑑表示的值, 和一组由𝑤∈𝑅𝑑表示的权重。 𝑥中的值根据权重𝑤的加权和, 可以表示为点积𝑥⊤𝑤。 当权重为非负数且和为1(即(∑𝑖=1𝑑𝑤𝑖=1))时, 点积表示加权平均(weighted average)。 将两个向量规范化得到单位长度后,点积表示它们夹角的余弦。 本节后面的内容将正式介绍长度(length)的概念。
矩阵-向量积
现在我们知道如何计算点积,可以开始理解矩阵-向量积(matrix-vector product)。
回顾分别在 :eqref:eq_matrix_def
和 :eqref:eq_vec_def
中定义的矩阵𝐴∈𝑅𝑚×𝑛和向量𝑥∈𝑅𝑛。
让我们将矩阵𝐴用它的行向量表示:
其中每个𝑎𝑖⊤∈𝑅𝑛都是行向量,表示矩阵的第𝑖行。
[矩阵向量积𝐴𝑥是一个长度为𝑚的列向量, 其第𝑖个元素是点积𝑎𝑖⊤𝑥]:
我们可以把一个矩阵𝐴∈𝑅𝑚×𝑛乘法看作一个从𝑅𝑛到𝑅𝑚向量的转换。
这些转换是非常有用的,例如可以用方阵的乘法来表示旋转。 后续章节将讲到,我们也可以使用矩阵-向量积来描述在给定前一层的值时, 求解神经网络每一层所需的复杂计算。
在代码中使用张量表示矩阵-向量积,我们使用mv
函数。
当我们为矩阵A
和向量x
调用torch.mv(A, x)
时,会执行矩阵-向量积。
注意,A
的列维数(沿轴1的长度)必须与x
的维数(其长度)相同。
矩阵-矩阵乘法
在掌握点积和矩阵-向量积的知识后, 那么矩阵-矩阵乘法(matrix-matrix multiplication)应该很简单。
假设有两个矩阵𝐴∈𝑅𝑛×𝑘和𝐵∈𝑅𝑘×𝑚:
用行向量𝑎𝑖⊤∈𝑅𝑘表示矩阵𝐴的第𝑖行,并让列向量𝑏𝑗∈𝑅𝑘作为矩阵𝐵的第𝑗列。
要生成矩阵积𝐶=𝐴𝐵,最简单的方法是考虑𝐴的行向量和𝐵的列向量:
当我们简单地将每个元素𝑐𝑖𝑗计算为点积𝑎𝑖⊤𝑏𝑗:
[我们可以将矩阵-矩阵乘法AB看作简单地执行𝑚次矩阵-向量积,并将结果拼接在一起,形成一个𝑛×𝑚矩阵]。
在下面的代码中,我们在A
和B
上执行矩阵乘法。 这里的A
是一个5行4列的矩阵,B
是一个4行3列的矩阵。 两者相乘后,我们得到了一个5行3列的矩阵。
矩阵-矩阵乘法可以简单地称为矩阵乘法,不应与"Hadamard积"混淆
范数
:label:subsec_lin-algebra-norms
线性代数中最有用的一些运算符是范数(norm)。
非正式地说,向量的范数是表示一个向量有多大。 这里考虑的大小(size)概念不涉及维度,而是分量的大小。
在线性代数中,向量范数是将向量映射到标量的函数𝑓。 给定任意向量𝑥,向量范数要满足一些属性。
第一个性质:
如果我们按常数因子𝛼缩放向量的所有元素, 其范数也会按相同常数因子的绝对值缩放
𝑓(𝛼𝑥)=|𝛼|𝑓(𝑥).
第二个性质:
熟悉的三角不等式 𝑓(𝑥+𝑦)≤𝑓(𝑥)+𝑓(𝑦).
第三个性质:
范数必须是非负的 𝑓(𝑥)≥0.
这是有道理的。因为在大多数情况下,任何东西的最小的大小是0。
最后一个性质要求范数最小为0,当且仅当向量全由0组成。
∀𝑖,[𝑥]𝑖=0⇔𝑓(𝑥)=0.
范数听起来很像距离的度量。 欧几里得距离和毕达哥拉斯定理中的非负性概念和三角不等式可能会给出一些启发。
欧几里得距离是一个𝐿2范数:
假设𝑛维向量𝑥中的元素是𝑥1,…,𝑥𝑛,其[**𝐿2范数是向量元素平方和的平方根**]
其中,在𝐿2范数中常常省略下标2,也就是说∥𝑥∥等同于∥𝑥∥2。 在代码中,我们可以按如下方式计算向量的𝐿2范数。
深度学习中更经常地使用𝐿2范数的平方,也会经常遇到[𝐿1范数,它表示为向量元素的绝对值之和:]
与𝐿2范数相比,𝐿1范数受异常值的影响较小。
为了计算𝐿1范数,我们将绝对值函数和按元素求和组合起来。
𝐿2范数和𝐿1范数都是更一般的𝐿𝑝范数的特例:
类似于向量的𝐿2范数,[矩阵]𝑋∈𝑅𝑚×𝑛的Frobenius范数(Frobenius norm)是矩阵元素平方和的平方根:
Frobenius范数满足向量范数的所有性质,它就像是矩阵形向量的𝐿2范数。
调用以下函数将计算矩阵的Frobenius范数。
范数和目标
:label:subsec_norms_and_objectives
在深度学习中,我们经常试图解决优化问题: 最大化分配给观测数据的概率; 最小化预测和真实观测之间的距离。 用向量表示物品(如单词、产品或新闻文章),以便最小化相似项目之间的距离,最大化不同项目之间的距离。 目标,或许是深度学习算法最重要的组成部分(除了数据),通常被表达为范数。
小结
- 标量、向量、矩阵和张量是线性代数中的基本数学对象。
- 向量泛化自标量,矩阵泛化自向量。
- 标量、向量、矩阵和张量分别具有零、一、二和任意数量的轴。
- 一个张量可以通过
sum
和mean
沿指定的轴降低维度。 - 两个矩阵的按元素乘法被称为他们的Hadamard积。它与矩阵乘法不同。
- 在深度学习中,我们经常使用范数,如𝐿1范数、𝐿2范数和Frobenius范数。
- 我们可以对标量、向量、矩阵和张量执行各种操作。
练习
- 证明一个矩阵𝐴的转置的转置是𝐴,即(𝐴⊤)⊤=𝐴。
- 给出两个矩阵𝐴和𝐵,证明“它们转置的和”等于“它们和的转置”,即𝐴⊤+𝐵⊤=(𝐴+𝐵)⊤。
- 给定任意方阵𝐴,𝐴+𝐴⊤总是对称的吗?为什么?
- 本节中定义了形状(2,3,4)的张量
X
。len(X)
的输出结果是什么? - 对于任意形状的张量
X
,len(X)
是否总是对应于X
特定轴的长度?这个轴是什么?轴0
- 运行
A/A.sum(axis=1)
,看看会发生什么。请分析一下原因? - 考虑一个具有形状(2,3,4)的张量,在轴0、1、2上的求和输出是什么形状?
- 为
linalg.norm
函数提供3个或更多轴的张量,并观察其输出。对于任意形状的张量这个函数计算得到什么?