LiuShuixiuのBlog
  • 学习笔记
  • My
  • 深度学习
LiuShuixiu
2025-03-15
目录

线性代数

# 线性代数

在介绍完如何存储和操作数据后，接下来将简要地回顾一下部分基本线性代数内容。 这些内容有助于读者了解和实现本书中介绍的大多数模型。 本节将介绍线性代数中的基本数学对象、算术和运算，并用数学符号和相应的代码实现来表示它们。

# 标量

如果你曾经在餐厅支付餐费，那么应该已经知道一些基本的线性代数，比如在数字间相加或相乘。 例如，北京的温度为

$52∘F$

（华氏度，除摄氏度外的另一种温度计量单位）。

严格来说，仅包含一个数值被称为_标量_（scalar）。 如果要将此华氏度值转换为更常用的摄氏度， 则可以计算表达式c=59(f−32)，并将f赋为52。 在此等式中，每一项（5、9和32）都是标量值。 符号c和f称为_变量_（variable），它们表示未知的标量值。

本书采用了数学表示法，其中标量变量由普通小写字母表示（例如，x、y和z）。 本书用R表示所有（连续）实数_标量的空间，之后将严格定义_空间（space）是什么， 但现在只要记住表达式x∈R是表示x是一个实值标量的正式形式。 符号∈称为“属于”，它表示“是集合中的成员”。 例如x,y∈{0,1}可以用来表明x和y是值只能为0或1的数字。

(标量由只有一个元素的张量表示)。 下面的代码将实例化两个标量，并执行一些熟悉的算术运算，即加法、乘法、除法和指数。

import torch

x = torch.tensor(3.0)
y = torch.tensor(2.0)

x + y, x * y, x / y, x**y
1
2
3
4
5
6

1
2
3
4
5
6

(tensor(5.), tensor(6.), tensor(1.5000), tensor(9.))
1
2

# 向量

[向量可以被视为标量值组成的列表]。 这些标量值被称为向量的_元素_（element）或_分量_（component）。 当向量表示数据集中的样本时，它们的值具有一定的现实意义。 例如，如果我们正在训练一个模型来预测贷款违约风险，可能会将每个申请人与一个向量相关联， 其分量与其收入、工作年限、过往违约次数和其他因素相对应。 如果我们正在研究医院患者可能面临的心脏病发作风险，可能会用一个向量来表示每个患者， 其分量为最近的生命体征、胆固醇水平、每天运动时间等。 在数学表示法中，向量通常记为粗体、小写的符号 （例如，x、y和z)）。

人们通过一维张量表示向量。一般来说，张量可以具有任意长度，取决于机器的内存限制。

x = torch.arange(4)
x
1
2

1
2

tensor([0, 1, 2, 3])
1
2

我们可以使用下标来引用向量的任一元素，例如可以通过xi来引用第i个元素。 注意，元素xi是一个标量，所以我们在引用它时不会加粗。 大量文献认为列向量是向量的默认方向，在本书中也是如此。 在数学中，向量x可以写为：

x=[x1x2⋮xn],

:eqlabel:eq_vec_def

其中x1,…,xn是向量的元素。在代码中，我们(通过张量的索引来访问任一元素)。

x[3]
1

1

tensor(3)
1
2

# 长度、维度和形状

向量只是一个数字数组，就像每个数组都有一个长度一样，每个向量也是如此。 在数学表示法中，如果我们想说一个向量x由n个实值标量组成， 可以将其表示为x∈Rn。 向量的长度通常称为向量的_维度_（dimension）。

与普通的Python数组一样，我们可以通过调用Python的内置len()函数来[访问张量的长度]。

len(x)
1

1

4
1
2

当用张量表示一个向量（只有一个轴）时，我们也可以通过.shape属性访问向量的长度。 形状（shape）是一个元素组，列出了张量沿每个轴的长度（维数）。 对于(只有一个轴的张量，形状只有一个元素。 )

x.shape
1

1

torch.Size([4])
1
2

请注意，维度（dimension）这个词在不同上下文时往往会有不同的含义，这经常会使人感到困惑。 为了清楚起见，我们在此明确一下： _向量_或_轴_的维度被用来表示_向量_或_轴_的长度，即向量或轴的元素数量。 然而，张量的维度用来表示张量具有的轴数。 在这个意义上，张量的某个轴的维数就是这个轴的长度。

# 矩阵

正如向量将标量从零阶推广到一阶，矩阵将向量从一阶推广到二阶。 矩阵，我们通常用粗体、大写字母来表示 （例如，X、Y和Z）， 在代码中表示为具有两个轴的张量。

数学表示法使用A∈Rm×n 来表示矩阵A，其由m行和n列的实值标量组成。 我们可以将任意矩阵A∈Rm×n视为一个表格， 其中每个元素aij属于第i行第j列：

A=[a11a12⋯a1na21a22⋯a2n⋮⋮⋱⋮am1am2⋯amn].

:eqlabel:eq_matrix_def

对于任意A∈Rm×n， A的形状是（m,n）或m×n。 当矩阵具有相同数量的行和列时，其形状将变为正方形； 因此，它被称为_方阵_（square matrix）。

当调用函数来实例化张量时， 我们可以[通过指定两个分量m和n来创建一个形状为m×n的矩阵]。

A = torch.arange(20).reshape(5, 4)
A
1
2

1
2

tensor([[ 0,  1,  2,  3],
        [ 4,  5,  6,  7],
        [ 8,  9, 10, 11],
        [12, 13, 14, 15],
        [16, 17, 18, 19]])
1
2
3
4
5
6

我们可以通过行索引（i）和列索引（j）来访问矩阵中的标量元素aij， 例如[A]ij。 如果没有给出矩阵A的标量元素，如在 :eqref:eq_matrix_def那样， 我们可以简单地使用矩阵A的小写字母索引下标aij 来引用[A]ij。 为了表示起来简单，只有在必要时才会将逗号插入到单独的索引中， 例如a2,3j和[A]2i−1,3。

当我们交换矩阵的行和列时，结果称为矩阵的_转置_（transpose）。 通常用a⊤来表示矩阵的转置，如果B=A⊤， 则对于任意i和j，都有bij=aji。 因此，在 :eqref:eq_matrix_def中的转置是一个形状为n×m的矩阵：

A⊤=[a11a21…am1a12a22…am2⋮⋮⋱⋮a1na2n…amn].

现在在代码中访问(矩阵的转置)。

A.T
1

1

tensor([[ 0,  4,  8, 12, 16],
        [ 1,  5,  9, 13, 17],
        [ 2,  6, 10, 14, 18],
        [ 3,  7, 11, 15, 19]])
1
2
3
4
5

作为方阵的一种特殊类型，[对称矩阵（symmetric matrix）A等于其转置：A=A⊤]。 这里定义一个对称矩阵B：

B = torch.tensor([[1, 2, 3], [2, 0, 4], [3, 4, 5]])
B
1
2

1
2

tensor([[1, 2, 3],
        [2, 0, 4],
        [3, 4, 5]])
1
2
3
4

现在我们将B与它的转置进行比较。

B == B.T
1

1

tensor([[True, True, True],
        [True, True, True],
        [True, True, True]])
1
2
3
4

矩阵是有用的数据结构：它们允许我们组织具有不同模式的数据。 例如，我们矩阵中的行可能对应于不同的房屋（数据样本），而列可能对应于不同的属性。 曾经使用过电子表格软件或已阅读过 :numref:sec_pandas的人，应该对此很熟悉。 因此，尽管单个向量的默认方向是列向量，但在表示表格数据集的矩阵中， 将每个数据样本作为矩阵中的行向量更为常见。 后面的章节将讲到这点，这种约定将支持常见的深度学习实践。 例如，沿着张量的最外轴，我们可以访问或遍历小批量的数据样本。

# 张量

[就像向量是标量的推广，矩阵是向量的推广一样，我们可以构建具有更多轴的数据结构]。 张量（本小节中的“张量”指代数对象）是描述具有任意数量轴的n维数组的通用方法。 例如，向量是一阶张量，矩阵是二阶张量。 张量用特殊字体的大写字母表示（例如，X、Y和Z）， 它们的索引机制（例如xijk和[X]1,2i−1,3）与矩阵类似。

当我们开始处理图像时，张量将变得更加重要，图像以n维数组形式出现， 其中3个轴对应于高度、宽度，以及一个_通道_（channel）轴， 用于表示颜色通道（红色、绿色和蓝色）。 现在先将高阶张量暂放一边，而是专注学习其基础知识。

X = torch.arange(24).reshape(2, 3, 4)
X
1
2

1
2

tensor([[[ 0,  1,  2,  3],
         [ 4,  5,  6,  7],
         [ 8,  9, 10, 11]],

        [[12, 13, 14, 15],
         [16, 17, 18, 19],
         [20, 21, 22, 23]]])
1
2
3
4
5
6
7
8

# 张量算法的基本性质

标量、向量、矩阵和任意数量轴的张量（本小节中的“张量”指代数对象）有一些实用的属性。 例如，从按元素操作的定义中可以注意到，任何按元素的一元运算都不会改变其操作数的形状。 同样，[给定具有相同形状的任意两个张量，任何按元素二元运算的结果都将是相同形状的张量]。 例如，将两个相同形状的矩阵相加，会在这两个矩阵上执行元素加法。

A = torch.arange(20, dtype=torch.float32).reshape(5, 4)
B = A.clone()  # 通过分配新内存，将A的一个副本分配给B
A, A + B
1
2
3

1
2
3

(tensor([[ 0.,  1.,  2.,  3.],
         [ 4.,  5.,  6.,  7.],
         [ 8.,  9., 10., 11.],
         [12., 13., 14., 15.],
         [16., 17., 18., 19.]]),
 tensor([[ 0.,  2.,  4.,  6.],
         [ 8., 10., 12., 14.],
         [16., 18., 20., 22.],
         [24., 26., 28., 30.],
         [32., 34., 36., 38.]]))
1
2
3
4
5
6
7
8
9
10
11

具体而言，[两个矩阵的按元素乘法称为_Hadamard积_（Hadamard product）（数学符号⊙）]。 对于矩阵B∈Rm×n， 其中第i行和第j列的元素是bij。 矩阵A（在 :eqref:eq_matrix_def中定义）和B的Hadamard积为：

A⊙B=[a11b11a12b12…a1nb1na21b21a22b22…a2nb2n⋮⋮⋱⋮am1bm1am2bm2…amnbmn].

A * B
1

1

tensor([[  0.,   1.,   4.,   9.],
        [ 16.,  25.,  36.,  49.],
        [ 64.,  81., 100., 121.],
        [144., 169., 196., 225.],
        [256., 289., 324., 361.]])
1
2
3
4
5
6

将张量乘以或加上一个标量不会改变张量的形状，其中张量的每个元素都将与标量相加或相乘。

a = 2
X = torch.arange(24).reshape(2, 3, 4)
a + X, (a * X).shape
1
2
3

1
2
3

(tensor([[[ 2,  3,  4,  5],
          [ 6,  7,  8,  9],
          [10, 11, 12, 13]],
 
         [[14, 15, 16, 17],
          [18, 19, 20, 21],
          [22, 23, 24, 25]]]),
 torch.Size([2, 3, 4]))
1
2
3
4
5
6
7
8
9

# 降维

🏷subseq_lin-alg-reduction

我们可以对任意张量进行的一个有用的操作是[计算其元素的和]。 数学表示法使用∑符号表示求和。 为了表示长度为d的向量中元素的总和，可以记为∑i=1dxi。 在代码中可以调用计算求和的函数：

x = torch.arange(4, dtype=torch.float32)
x, x.sum()
1
2

1
2

(tensor([0., 1., 2., 3.]), tensor(6.))
1
2

我们可以(表示任意形状张量的元素和)。 例如，矩阵A中元素的和可以记为∑i=1m∑j=1naij。

A.shape, A.sum()
1

1

(torch.Size([5, 4]), tensor(190.))
1
2

默认情况下，调用求和函数会沿所有的轴降低张量的维度，使它变为一个标量。 我们还可以[指定张量沿哪一个轴来通过求和降低维度]。 以矩阵为例，为了通过求和所有行的元素来降维（轴0），可以在调用函数时指定axis=0。 由于输入矩阵沿0轴降维以生成输出向量，因此输入轴0的维数在输出形状中消失。

A_sum_axis0 = A.sum(axis=0)
A_sum_axis0, A_sum_axis0.shape
1
2

1
2

(tensor([40., 45., 50., 55.]), torch.Size([4]))
1
2

指定axis=1将通过汇总所有列的元素降维（轴1）。因此，输入轴1的维数在输出形状中消失。

A_sum_axis1 = A.sum(axis=1)
A_sum_axis1, A_sum_axis1.shape
1
2

1
2

(tensor([ 6., 22., 38., 54., 70.]), torch.Size([5]))
1
2

沿着行和列对矩阵求和，等价于对矩阵的所有元素进行求和。

A.sum(axis=[0, 1])  # 结果和A.sum()相同
1

1

tensor(190.)
1
2

[一个与求和相关的量是_平均值_（mean或average）]。 我们通过将总和除以元素总数来计算平均值。 在代码中，我们可以调用函数来计算任意形状张量的平均值。

A.mean(), A.sum() / A.numel()
1

1

(tensor(9.5000), tensor(9.5000))
1
2

同样，计算平均值的函数也可以沿指定轴降低张量的维度。

A.mean(axis=0), A.sum(axis=0) / A.shape[0]
1

1

(tensor([ 8.,  9., 10., 11.]), tensor([ 8.,  9., 10., 11.]))
1
2

# 非降维求和

🏷subseq_lin-alg-non-reduction

但是，有时在调用函数来[计算总和或均值时保持轴数不变]会很有用。

sum_A = A.sum(axis=1, keepdims=True)
sum_A
1
2

1
2

tensor([[ 6.],
        [22.],
        [38.],
        [54.],
        [70.]])
1
2
3
4
5
6

例如，由于sum_A在对每行进行求和后仍保持两个轴，我们可以(通过广播将A除以sum_A)。

A / sum_A
1

1

tensor([[0.0000, 0.1667, 0.3333, 0.5000],
        [0.1818, 0.2273, 0.2727, 0.3182],
        [0.2105, 0.2368, 0.2632, 0.2895],
        [0.2222, 0.2407, 0.2593, 0.2778],
        [0.2286, 0.2429, 0.2571, 0.2714]])
1
2
3
4
5
6

如果我们想沿[某个轴计算A元素的累积总和]， 比如axis=0（按行计算），可以调用cumsum函数。 此函数不会沿任何轴降低输入张量的维度。

A.cumsum(axis=0)
1

1

tensor([[ 0.,  1.,  2.,  3.],
        [ 4.,  6.,  8., 10.],
        [12., 15., 18., 21.],
        [24., 28., 32., 36.],
        [40., 45., 50., 55.]])
1
2
3
4
5
6

# 点积（Dot Product）

我们已经学习了按元素操作、求和及平均值。 另一个最基本的操作之一是点积。 给定两个向量 x,y∈Rd， 它们的_点积_（dot product）x⊤y （或⟨x,y⟩） 是相同位置的按元素乘积的和：x⊤y=∑i=1dxiyi。

[~点积是相同位置的按元素乘积的和~]

y = torch.ones(4, dtype = torch.float32)
x, y, torch.dot(x, y)
1
2

1
2

(tensor([0., 1., 2., 3.]), tensor([1., 1., 1., 1.]), tensor(6.))
1
2

注意，(我们可以通过执行按元素乘法，然后进行求和来表示两个向量的点积)：

torch.sum(x * y)
1

1

tensor(6.)
1
2

点积在很多场合都很有用。 例如，给定一组由向量 x∈Rd 表示的值， 和一组由w∈Rd表示的权重。 x中的值根据权重w的加权和， 可以表示为点积x⊤w。 当权重为非负数且和为1（即(∑i=1dwi=1)）时， 点积表示_加权平均_（weighted average）。 将两个向量规范化得到单位长度后，点积表示它们夹角的余弦。 本节后面的内容将正式介绍_长度_（length）的概念。

# 矩阵-向量积

现在我们知道如何计算点积，可以开始理解_矩阵-向量积_（matrix-vector product）。 回顾分别在 :eqref:eq_matrix_def和 :eqref:eq_vec_def中定义的矩阵A∈Rm×n和向量x∈Rn。 让我们将矩阵A用它的行向量表示：

A=[a1⊤a2⊤⋮am⊤],

其中每个ai⊤∈Rn都是行向量，表示矩阵的第i行。 [矩阵向量积Ax是一个长度为m的列向量， 其第i个元素是点积ai⊤x]：

Ax=[a1⊤a2⊤⋮am⊤]x=[a1⊤xa2⊤x⋮am⊤x].

我们可以把一个矩阵A∈Rm×n乘法看作一个从Rn到Rm向量的转换。 这些转换是非常有用的，例如可以用方阵的乘法来表示旋转。 后续章节将讲到，我们也可以使用矩阵-向量积来描述在给定前一层的值时， 求解神经网络每一层所需的复杂计算。

在代码中使用张量表示矩阵-向量积，我们使用mv函数。 当我们为矩阵A和向量x调用torch.mv(A, x)时，会执行矩阵-向量积。 注意，A的列维数（沿轴1的长度）必须与x的维数（其长度）相同。

A.shape, x.shape, torch.mv(A, x)
1

1

(torch.Size([5, 4]), torch.Size([4]), tensor([ 14.,  38.,  62.,  86., 110.]))
1
2

# 矩阵-矩阵乘法

在掌握点积和矩阵-向量积的知识后， 那么矩阵-矩阵乘法（matrix-matrix multiplication）应该很简单。

假设有两个矩阵A∈Rn×k和B∈Rk×m：

A=[a11a12⋯a1ka21a22⋯a2k⋮⋮⋱⋮an1an2⋯ank],B=[b11b12⋯b1mb21b22⋯b2m⋮⋮⋱⋮bk1bk2⋯bkm].

用行向量ai⊤∈Rk表示矩阵A的第i行，并让列向量bj∈Rk作为矩阵B的第j列。要生成矩阵积C=AB，最简单的方法是考虑A的行向量和B的列向量:

A=[a1⊤a2⊤⋮an⊤],B=[b1b2⋯bm].

当我们简单地将每个元素cij计算为点积ai⊤bj:

C=AB=[a1⊤a2⊤⋮an⊤][b1b2⋯bm]=[a1⊤b1a1⊤b2⋯a1⊤bma2⊤b1a2⊤b2⋯a2⊤bm⋮⋮⋱⋮an⊤b1an⊤b2⋯an⊤bm].

[我们可以将矩阵-矩阵乘法AB看作简单地执行m次矩阵-向量积，并将结果拼接在一起，形成一个n×m矩阵]。 在下面的代码中，我们在A和B上执行矩阵乘法。 这里的A是一个5行4列的矩阵，B是一个4行3列的矩阵。 两者相乘后，我们得到了一个5行3列的矩阵。

B = torch.ones(4, 3)
torch.mm(A, B)
1
2

1
2

tensor([[ 6.,  6.,  6.],
        [22., 22., 22.],
        [38., 38., 38.],
        [54., 54., 54.],
        [70., 70., 70.]])
1
2
3
4
5
6

矩阵-矩阵乘法可以简单地称为矩阵乘法，不应与"Hadamard积"混淆。

# 范数

🏷subsec_lin-algebra-norms

线性代数中最有用的一些运算符是_范数_（norm）。 非正式地说，向量的_范数_是表示一个向量有多大。 这里考虑的_大小_（size）概念不涉及维度，而是分量的大小。

在线性代数中，向量范数是将向量映射到标量的函数f。 给定任意向量x，向量范数要满足一些属性。 第一个性质是：如果我们按常数因子α缩放向量的所有元素， 其范数也会按相同常数因子的_绝对值_缩放：

f(αx)=|α|f(x).

第二个性质是熟悉的三角不等式:

f(x+y)≤f(x)+f(y).

第三个性质简单地说范数必须是非负的:

f(x)≥0.

这是有道理的。因为在大多数情况下，任何东西的最小的_大小_是0。 最后一个性质要求范数最小为0，当且仅当向量全由0组成。

∀i,[x]i=0⇔f(x)=0.

范数听起来很像距离的度量。 欧几里得距离和毕达哥拉斯定理中的非负性概念和三角不等式可能会给出一些启发。 事实上，欧几里得距离是一个L2范数： 假设n维向量x中的元素是x1,…,xn，其[L2_范数_是向量元素平方和的平方根： ]

($$|\mathbf{x}|2 = \sqrt{\sum{i=1}^n x_i^2},$$)

其中，在L2范数中常常省略下标2，也就是说∥x∥等同于∥x∥2。 在代码中，我们可以按如下方式计算向量的L2范数。

u = torch.tensor([3.0, -4.0])
torch.norm(u)
1
2

1
2

tensor(5.)
1
2

深度学习中更经常地使用L2范数的平方，也会经常遇到[L1范数，它表示为向量元素的绝对值之和： ]

($$|\mathbf{x}|1 = \sum{i=1}^n \left|x_i \right|.$$)

与L2范数相比，L1范数受异常值的影响较小。 为了计算L1范数，我们将绝对值函数和按元素求和组合起来。

torch.abs(u).sum()
1

1

tensor(7.)
1
2

L2范数和L1范数都是更一般的Lp范数的特例：

∥x∥p=(∑i=1n|xi|p)1/p.

类似于向量的L2范数，[矩阵]X∈Rm×n(的_Frobenius范数_（Frobenius norm）是矩阵元素平方和的平方根： )

($$|\mathbf{X}|F = \sqrt{\sum{i=1}^m \sum_{j=1}^n x_{ij}^2}.$$)

Frobenius范数满足向量范数的所有性质，它就像是矩阵形向量的L2范数。 调用以下函数将计算矩阵的Frobenius范数。

torch.norm(torch.ones((4, 9)))
1

1

tensor(6.)
1
2

# 范数和目标

🏷subsec_norms_and_objectives

在深度学习中，我们经常试图解决优化问题： _最大化_分配给观测数据的概率; _最小化_预测和真实观测之间的距离。 用向量表示物品（如单词、产品或新闻文章），以便最小化相似项目之间的距离，最大化不同项目之间的距离。 目标，或许是深度学习算法最重要的组成部分（除了数据），通常被表达为范数。

# 关于线性代数的更多信息

仅用一节，我们就教会了阅读本书所需的、用以理解现代深度学习的线性代数。 线性代数还有很多，其中很多数学对于机器学习非常有用。 例如，矩阵可以分解为因子，这些分解可以显示真实世界数据集中的低维结构。 机器学习的整个子领域都侧重于使用矩阵分解及其向高阶张量的泛化，来发现数据集中的结构并解决预测问题。 当开始动手尝试并在真实数据集上应用了有效的机器学习模型，你会更倾向于学习更多数学。 因此，这一节到此结束，本书将在后面介绍更多数学知识。

如果渴望了解有关线性代数的更多信息，可以参考线性代数运算的在线附录 (opens new window) (opens new window)或其他优秀资源 :cite:Strang.1993,Kolter.2008,Petersen.Pedersen.ea.2008。

# 小结

  • 标量、向量、矩阵和张量是线性代数中的基本数学对象。
  • 向量泛化自标量，矩阵泛化自向量。
  • 标量、向量、矩阵和张量分别具有零、一、二和任意数量的轴。
  • 一个张量可以通过sum和mean沿指定的轴降低维度。
  • 两个矩阵的按元素乘法被称为他们的Hadamard积。它与矩阵乘法不同。
  • 在深度学习中，我们经常使用范数，如L1范数、L2范数和Frobenius范数。
  • 我们可以对标量、向量、矩阵和张量执行各种操作。

# 练习

  1. 证明一个矩阵A的转置的转置是A，即(A⊤)⊤=A。
  2. 给出两个矩阵A和B，证明“它们转置的和”等于“它们和的转置”，即A⊤+B⊤=(A+B)⊤。
  3. 给定任意方阵A，A+A⊤总是对称的吗?为什么?
  4. 本节中定义了形状(2,3,4)的张量X。len(X)的输出结果是什么？
  5. 对于任意形状的张量X,len(X)是否总是对应于X特定轴的长度?这个轴是什么?
  6. 运行A/A.sum(axis=1)，看看会发生什么。请分析一下原因？
  7. 考虑一个具有形状(2,3,4)的张量，在轴0、1、2上的求和输出是什么形状?
  8. 为linalg.norm函数提供3个或更多轴的张量，并观察其输出。对于任意形状的张量这个函数计算得到什么?
上次更新: 2025/03/15, 18:48:43
小程序笔记
数值稳定性和模型初始化

← 数值稳定性和模型初始化→

最近更新
01
Suwayomi漫画（电脑版）
03-23
02
常用软件的安装
03-23
03
Mihon漫画
03-16
更多文章>
Theme by Vdoing | Copyright © 2023-2025 LiuShuixiu | MIT License