torch.nn 到底是什么?

译者:lhc741

作者:Jeremy Howard,fast.ai。感谢Rachel Thomas和Francisco Ingham的帮助和支持。

我们推荐使用notebook来运行这个教程,而不是脚本,点击这里下载notebook(.ipynb)文件。

Pytorch提供了torch.nntorch.optimDatasetDataLoader这些设计优雅的模块和类以帮助使用者创建和训练神经网络。 为了最大化利用这些模块和类的功能,并使用它们做出适用于你所研究问题的模型,你需要真正理解他们是如何工作的。 为了做到这一点,我们首先基于MNIST数据集训练一个没有任何特征的简单神经网络。 最开始我们只会用到PyTorch中最基本的tensor功能,然后我们将会逐渐的从torch.nntorch.optimDatasetDataLoader中选择一个特征加入到模型中,来展示新加入的特征会对模型产生什么样的效果,以及它是如何使模型变得更简洁或更灵活。

在这个教程中,我们假设你已经安装好了PyTorch,并且已经熟悉了基本的tensor运算。(如果你熟悉Numpy的数组运算,你将会发现这里用到的PyTorch tensor运算和numpy几乎是一样的)

MNIST数据安装

我们将要使用经典的MNIST数据集,这个数据集由手写数字(0到9)的黑白图片组成。

我们将使用pathlib来处理文件路径的相关操作(python3中的一个标准库),使用request来下载数据集。 我们只会在用到相关库的时候进行引用,这样你就可以明确在每个操作中用到了哪些库。

  1. from pathlib import Path
  2. import requests
  3. DATA_PATH = Path("data")
  4. PATH = DATA_PATH / "mnist"
  5. PATH.mkdir(parents=True, exist_ok=True)
  6. URL = "http://deeplearning.net/data/mnist/"
  7. FILENAME = "mnist.pkl.gz"
  8. if not (PATH / FILENAME).exists():
  9. content = requests.get(URL + FILENAME).content
  10. (PATH / FILENAME).open("wb").write(content)`

该数据集采用numpy数组格式,并已使用pickle存储,pickle是一个用来把数据序列化为python特定格式的库。

  1. import pickle
  2. import gzip
  3. with gzip.open((PATH / FILENAME).as_posix(), "rb") as f:
  4. ((x_train, y_train), (x_valid, y_valid), _) = pickle.load(f, encoding="latin-1")

每一幅图像都是28 x 28的,并被拉平成长度为784(=28x28)的一行。 我们以其中一个为例展示一下,首先需要将这个一行的数据重新变换为一个2d的数据。

  1. from matplotlib import pyplot
  2. import numpy as np
  3. pyplot.imshow(x_train[0].reshape((28, 28)), cmap="gray")
  4. print(x_train.shape)

image

输出:

  1. (50000, 784)

PyTorch使用torch.tensor,而不是numpy数组,所以我们需要将数据转换。

  1. import torch
  2. x_train, y_train, x_valid, y_valid = map(
  3. torch.tensor, (x_train, y_train, x_valid, y_valid)
  4. )
  5. n, c = x_train.shape
  6. x_train, x_train.shape, y_train.min(), y_train.max()
  7. print(x_train, y_train)
  8. print(x_train.shape)
  9. print(y_train.min(), y_train.max())

输出:

  1. tensor([[0., 0., 0., ..., 0., 0., 0.],
  2. [0., 0., 0., ..., 0., 0., 0.],
  3. [0., 0., 0., ..., 0., 0., 0.],
  4. ...,
  5. [0., 0., 0., ..., 0., 0., 0.],
  6. [0., 0., 0., ..., 0., 0., 0.],
  7. [0., 0., 0., ..., 0., 0., 0.]]) tensor([5, 0, 4, ..., 8, 4, 8])
  8. torch.Size([50000, 784])
  9. tensor(0) tensor(9)

神经网络从零开始(不使用torch.nn)

我们先来建立一个只使用PyTorch张量运算的模型。 我们假设你已经熟悉神经网络的基础。(如果你还不熟悉,可以访问course.fast.ai进行学习)。

PyTorch提供创建随机数填充或全零填充张量的方法,我们使用该方法初始化一个简单线性模型的权重和偏置。 这两个都是普通的张量,但它们有一个特殊的附加条件:设置需要计算梯度的参数为True。这样PyTorch就会记录所有与这个张量相关的运算,使其能在反向传播阶段自动计算梯度。

对于weights而言,由于我们希望初始化张量过程中存在梯度,所以我们在初始化之后设置requires_grad。(注意:尾缀为_的方法在PyTorch中表示这个操作会被立即被执行。)

  • 注意:

我们以Xavier初始化方法(每个元素都除以1/sqrt(n))为例来对权重进行初始化。

  1. import math
  2. weights = torch.randn(784, 10) / math.sqrt(784)
  3. weights.requires_grad_()
  4. bias = torch.zeros(10, requires_grad=True)

多亏了PyTorch具有自动梯度计算功能,我们可以使用Python中任何标准函数(或者可调用对象)来创建模型! 因此,让我们编写一个普通的矩阵乘法和广播加法建立一个简单的线性模型。 我们还需要一个激活函数,所以我们编写并使用一个log_softmax函数。 请记住:尽管Pytorch提供了许多预先编写好的损失函数、激活函数等等,你仍然可以使用纯python轻松实现你自己的函数。 Pytorch甚至可以自动地为你的函数创建快速的GPU代码或向量化的CPU代码。

  1. def log_softmax(x):
  2. return x - x.exp().sum(-1).log().unsqueeze(-1)
  3. def model(xb):
  4. return log_softmax(xb @ weights + bias)

在上面的一段代码中,@表示点积运算符。我们将调用我们的函数计算一个批次的数据(本例中为64幅图像)。 这是一次模型前向传递的过程。 请注意,因为我们使用了随机数来初始化权重,所以在这个阶段我们的预测值并不会比随机的更好。

  1. bs = 64 # 一批数据个数
  2. xb = x_train[0:bs] # 从x获取一小批数据
  3. preds = model(xb) # 预测值
  4. preds[0], preds.shape
  5. print(preds[0], preds.shape)

输出:

  1. tensor([-2.4513, -2.5024, -2.0599, -3.1052, -3.2918, -2.2665, -1.9007, -2.2588,
  2. -2.0149, -2.0287], grad_fn=<SelectBackward>) torch.Size([64, 10])

可以从上面的结果不难看出,张量preds不仅包括了张量值,还包括了梯度函数。这个梯度函数我们可以在后面的反向传播阶段用到。

下面我们来实现一个负的对数似然函数(Negative log-likehood)作为损失函数(同样也使用纯python实现):

  1. def nll(input, target):
  2. return -input[range(target.shape[0]), target].mean()
  3. loss_func = nll

让我们查看下随机模型的损失值,这样我们就可以确认在执行反向传播的步骤后,模型的预测效果是否有了改进。

  1. yb = y_train[0:bs]
  2. print(loss_func(preds, yb))

输出:

  1. tensor(2.3620, grad_fn=<NegBackward>)

我们再来实现一个用来计算模型准确率的函数。对于每次预测,我们规定如果预测结果中概率最大的数字和图片实际对应的数字是相同的,那么这次预测就是正确的。

  1. def accuracy(out, yb):
  2. preds = torch.argmax(out, dim=1)
  3. return (preds == yb).float().mean()

我们先来看一下被随机初始化的模型的准确率,这样我们就可以看到损失值降低的时候准确率是否提高了。

  1. print(accuracy(preds, yb))

输出:

  1. tensor(0.0938)

现在我们可以运行一个完整的训练步骤了,每次迭代,我们会进行以下几个操作:

  • 从全部数据中选择一小批数据(大小为bs
  • 使用模型进行预测
  • 计算当前预测的损失值
  • 使用loss.backward()更新模型中的梯度,在这个例子中,更新的是weightsbias

现在,我们来利用计算出的梯度对权值和偏置项进行更新,因为我们不希望这一步的操作被用于下一次迭代的梯度计算,所以我们在torch.no_grad()这个上下文管理器中完成。想要了解更多PyTorch Autograd记录操作现,可以点击这里

接下来,我们将梯度设置为0,来为下一次循环做准备。否则我们的梯度将会记录所有已经执行过的运算(如,loss.backward()会将梯度变化值直接与变量已有值进行累加,而不是替换变量原有的值)。

小贴士

您可以使用标准python调试器对PyTorch代码进行单步调试,从而在每一步检查不同的变量值。取消下面的set_trace()来尝试该功能。

  1. from IPython.core.debugger import set_trace
  2. lr = 0.5 # 学习率
  3. epochs = 2 # 训练的轮数
  4. for epoch in range(epochs):
  5. for i in range((n - 1) // bs + 1):
  6. # set_trace()
  7. start_i = i * bs
  8. end_i = start_i + bs
  9. xb = x_train[start_i:end_i]
  10. yb = y_train[start_i:end_i]
  11. pred = model(xb)
  12. loss = loss_func(pred, yb)
  13. loss.backward()
  14. with torch.no_grad():
  15. weights -= weights.grad * lr
  16. bias -= bias.grad * lr
  17. weights.grad.zero_()
  18. bias.grad.zero_()

目前为止,我们已经从零开始完成了建立和训练一个最小的神经网络(因为我们建立的logistic回归模型不包含隐层)!

现在,我们来看一下模型的损失值和准确率,并于我们之前输出的值进行比较。结果正如我们预期的,损失值下降,准确率提高。

  1. print(loss_func(model(xb), yb), accuracy(model(xb), yb))

输出:

  1. tensor(0.0822, grad_fn=<NegBackward>) tensor(1.)

torch.nn.functional的使用

现在,我们要对前面的代码进行重构,使代码在完成相同功能的同时,用PyTorch的nn来使代码变得更加简洁和灵活。 从现在开始,接下来的每一步我们都会使代码变得更短,更好理解或更灵活。

要进行的第一步也是最简单的一步,是使用torch.nn.functional(通过会在引用时用F表示)中的函数替换我们自己的激活函数和损失函数使代码变得更短。 这个模块包含了torch.nn库中的所有函数(这个库的其它部分是各种类),所以在这个模块中还会找到其它便于建立神经网络的函数,比如池化函数。(模块中还包含卷积函数,线性函数等等,不过在后面的内容中我们会看到,这些操作使用库中的其它部分会更好。)

如果你使用负对数似然损失和对数柔性最大值(softmax)激活函数,PyTorch有一个结合了这两个函数的简单函数F.cross_entropy供你使用,这样我们就可以删掉模型中的激活函数。

  1. import torch.nn.functional as F
  2. loss_func = F.cross_entropy
  3. def model(xb):
  4. return xb @ weights + bias

注意在model函数中我们不再调用log_softmax。现在我们来确认一下损失值和准确率与之前相同。

  1. print(loss_func(model(xb), yb), accuracy(model(xb), yb))

输出:

  1. tensor(0.0822, grad_fn=<NllLossBackward>) tensor(1.)

使用nn.Module进行重构

接下来,我们将会用到nn.Modelnn.Parameter来完成一个更加清晰简洁的训练循环。我们继承nn.Module(它是一个能够跟踪状态的类)。在这个例子中,我们想要新建一个类,实现存储权重,偏置和前向传播步骤中所有用到方法。nn.Module包含了许多属性和方法(比如.parameters().zero_grad()),我们会在后面用到。

注意

nn.Module(M大写)是一个PyTorch中特有的概念,它是一个会经常用到的类。不要和Python中module(m小写)混淆,module是一个可以被引入的Python代码文件。

  1. from torch import nn
  2. class Mnist_Logistic(nn.Module):
  3. def __init__(self):
  4. super().__init__()
  5. self.weights = nn.Parameter(torch.randn(784, 10) / math.sqrt(784))
  6. self.bias = nn.Parameter(torch.zeros(10))
  7. def forward(self, xb):
  8. return xb @ self.weights + self.bias

既然现在我们要使用一个对象而不是函数,我们要先对模型进行实例化。

  1. model = Mnist_Logistic()

现在我们可以像之前那样计算损失值了。注意nn.Module对象的使用方式很像函数(例如它们是可调用的),但是PyTorch将会自动调用我们的forward函数

  1. print(loss_func(model(xb), yb))

输出:

  1. tensor(2.2082, grad_fn=<NllLossBackward>)

之前在每个训练循环中,我们通过变量名对每个变量的值进行更新,并手动的将每个变量的梯度置为0,像这样:

  1. with torch.no_grad():
  2. weights -= weights.grad * lr
  3. bias -= bias.grad * lr
  4. weights.grad.zero_()
  5. bias.grad.zero_()

现在我们可以利用model.parameters()model.zero_grad()(这两个都是PyTorch定义在nn.Module中的)使这些步骤变得更加简洁并且更不容易忘记更新部分参数,尤其是模型很复杂的情况:

  1. with torch.no_grad():
  2. for p in model.parameters(): p -= p.grad * lr
  3. model.zero_grad()

下面我们把训练循环封装进fit函数中,这样就能在后面再次运行。

  1. def fit():
  2. for epoch in range(epochs):
  3. for i in range((n - 1) // bs + 1):
  4. start_i = i * bs
  5. end_i = start_i + bs
  6. xb = x_train[start_i:end_i]
  7. yb = y_train[start_i:end_i]
  8. pred = model(xb)
  9. loss = loss_func(pred, yb)
  10. loss.backward()
  11. with torch.no_grad():
  12. for p in model.parameters():
  13. p -= p.grad * lr
  14. model.zero_grad()
  15. fit()

我们来再次检查一下我们的损失值是否下降。

  1. print(loss_func(model(xb), yb))

输出:

  1. tensor(0.0812, grad_fn=<NllLossBackward>)

使用nn.Linear进行重构

我们继续对代码进行重构。我们将用PyTorch中的nn.Linear代替手动定义和初始化self.weightsself.bias以及计算xb @ self.weights + self.bias, 因为nn.Linear可以完成这些操作。 PyTorch中预设了很多类型的神经网络层,使用它们可以极大的简化我们的代码,通常还会带来速度上的提升。

  1. class Mnist_Logistic(nn.Module):
  2. def __init__(self):
  3. super().__init__()
  4. self.lin = nn.Linear(784, 10)
  5. def forward(self, xb):
  6. return self.lin(xb)

我们初始化模型并像之前那样计算损失值:

  1. model = Mnist_Logistic()
  2. print(loss_func(model(xb), yb))

输出:

  1. tensor(2.2731, grad_fn=<NllLossBackward>)

我们仍然可以像之前那样使用fit函数:

  1. fit()
  2. print(loss_func(model(xb), yb))

输出:

  1. tensor(0.0820, grad_fn=<NllLossBackward>)

使用optim进行重构

PyTorch还有一个包含很多优化算法的包————torch.optim。我们可以使用优化器中的step方法执行前向传播过程中的步骤来替换手动更新每个参数。

这个方法将允许我们替换之前手动编写的优化步骤:

  1. with torch.no_grad():
  2. for p in model.parameters(): p -= p.grad * lr
  3. model.zero_grad()

替换后如下:

  1. opt.step()
  2. opt.zero_grad()

(optim.zero_grad()将梯度重置为0,我们需要在计算下一次梯度之前调用它)

  1. from torch import optim

我们将建立模型和优化器的步骤定义为一个小函数方便将来复用。

  1. def get_model():
  2. model = Mnist_Logistic()
  3. return model, optim.SGD(model.parameters(), lr=lr)
  4. model, opt = get_model()
  5. print(loss_func(model(xb), yb))
  6. for epoch in range(epochs):
  7. for i in range((n - 1) // bs + 1):
  8. start_i = i * bs
  9. end_i = start_i + bs
  10. xb = x_train[start_i:end_i]
  11. yb = y_train[start_i:end_i]
  12. pred = model(xb)
  13. loss = loss_func(pred, yb)
  14. loss.backward()
  15. opt.step()
  16. opt.zero_grad()
  17. print(loss_func(model(xb), yb))

输出:

  1. tensor(2.3785, grad_fn=<NllLossBackward>)
  2. tensor(0.0802, grad_fn=<NllLossBackward>)

使用Dataset进行重构

Pytorch包含一个Dataset抽象类。Dataset可以是任何东西,但它始终包含一个__len__函数(通过Python中的标准函数len调用)和一个用来索引到内容中的__getitem__函数。 这篇教程以创建Dataset的自定义子类FacialLandmarkDataset为例进行介绍。

PyTorch中的TensorDataset是一个封装了张量的Dataset。通过定义长度和索引的方式,是我们可以对张量的第一维进行迭代,索引和切片。这将使我们在训练中,获取同一行中的自变量和因变量更加容易。

  1. from torch.utils.data import TensorDataset

可以把x_trainy_train中的数据合并成一个简单的TensorDataset,这样就可以方便的进行迭代和切片操作。

  1. train_ds = TensorDataset(x_train, y_train)

之前,我们不得不分别对x和y的值进行迭代循环。

  1. xb = x_train[start_i:end_i]
  2. yb = y_train[start_i:end_i]

现在我们可以将这两步合二为一了。

  1. xb,yb = train_ds[i*bs : i*bs+bs]
  1. model, opt = get_model()
  2. for epoch in range(epochs):
  3. for i in range((n - 1) // bs + 1):
  4. xb, yb = train_ds[i * bs: i * bs + bs]
  5. pred = model(xb)
  6. loss = loss_func(pred, yb)
  7. loss.backward()
  8. opt.step()
  9. opt.zero_grad()
  10. print(loss_func(model(xb), yb))

输出:

  1. tensor(0.0817, grad_fn=<NllLossBackward>)

使用DataLoader进行重构

PyTorch的DataLoader负责批量数据管理,你可以使用任意的Dataset创建一个DataLoaderDataLoader使得对批量数据的迭代更容易。DataLoader自动的为我们提供每一小批量的数据来代替切片的方式train_ds[i*bs : i*bs+bs]

  1. from torch.utils.data import DataLoader
  2. train_ds = TensorDataset(x_train, y_train)
  3. train_dl = DataLoader(train_ds, batch_size=bs)

之前我们像下面这样按批(xb,yb)对数据进行迭代:

  1. for i in range((n-1)//bs + 1):
  2. xb,yb = train_ds[i*bs : i*bs+bs]
  3. pred = model(xb)

现在我们的循环变得更加简洁,因为使用了data loader来自动获取数据。

  1. for xb,yb in train_dl:
  2. pred = model(xb)
  1. model, opt = get_model()
  2. for epoch in range(epochs):
  3. for xb, yb in train_dl:
  4. pred = model(xb)
  5. loss = loss_func(pred, yb)
  6. loss.backward()
  7. opt.step()
  8. opt.zero_grad()
  9. print(loss_func(model(xb), yb))

输出:

  1. tensor(0.0817, grad_fn=<NllLossBackward>)

多亏PyTorch中的nn.Modulenn.ParameterDatasetDataLoader,我们的训练代码变得非常简洁易懂。下面我们来试着增加一些用于提高模型效率所必需的的基本特征。

增加验证集

在第一部分,我们仅仅是试着为我们的训练集构建一个合理的训练步骤,但实际上,我们始终应该有一个验证集来确认模型是否过拟合。

打乱训练数据的顺序通常是避免不同批数据中存在相关性和过拟合的重要步骤。另一方面,无论是否打乱顺序计算出的验证集损失值都是一样的。鉴于打乱顺序还会消耗额外的时间,所以打乱验证集数据是没有任何意义的。

我们在验证集上用到的每批数据的数量是训练集的两倍,这是因为在验证集上不需要进行反向传播,这样就会占用较小的内存(因为它并不需要储存梯度)。我们利用了这一点,使用了更大的batchsize,更快的计算出了损失值。

  1. train_ds = TensorDataset(x_train, y_train)
  2. train_dl = DataLoader(train_ds, batch_size=bs, shuffle=True)
  3. valid_ds = TensorDataset(x_valid, y_valid)
  4. valid_dl = DataLoader(valid_ds, batch_size=bs * 2)

我们将会在每轮(epoch)结束后计算并输出验证集上的损失值。

(注意:在训练前我们总是会调用model.train()函数,在推断之前调用model.eval()函数,因为这些会被nn.BatchNorm2dnn.Dropout等层使用,确保在不同阶段的准确性。)

  1. model, opt = get_model()
  2. for epoch in range(epochs):
  3. model.train()
  4. for xb, yb in train_dl:
  5. pred = model(xb)
  6. loss = loss_func(pred, yb)
  7. loss.backward()
  8. opt.step()
  9. opt.zero_grad()
  10. model.eval()
  11. with torch.no_grad():
  12. valid_loss = sum(loss_func(model(xb), yb) for xb, yb in valid_dl)
  13. print(epoch, valid_loss / len(valid_dl))

输出:

  1. 0 tensor(0.2999)
  2. 1 tensor(0.2742)

编写fit()和get_data()函数

现在我们来重构一下我们自己的函数。 我们在计算训练集和验证集上的损失值时执行了差不多的过程两次,因此我们将这部分代码提炼成一个函数loss_batch,用来计算每个批的损失值。

我们为训练集传递一个优化器参数来执行反向传播。对于验证集我们不传优化器参数,这样就不会执行反向传播。

  1. def loss_batch(model, loss_func, xb, yb, opt=None):
  2. loss = loss_func(model(xb), yb)
  3. if opt is not None:
  4. loss.backward()
  5. opt.step()
  6. opt.zero_grad()
  7. return loss.item(), len(xb)

fit执行了训练模型的必要操作,并在每轮(epoch)结束后计算模型在训练集和测试集上的损失。

  1. import numpy as np
  2. def fit(epochs, model, loss_func, opt, train_dl, valid_dl):
  3. for epoch in range(epochs):
  4. model.train()
  5. for xb, yb in train_dl:
  6. loss_batch(model, loss_func, xb, yb, opt)
  7. model.eval()
  8. with torch.no_grad():
  9. losses, nums = zip(
  10. *[loss_batch(model, loss_func, xb, yb) for xb, yb in valid_dl]
  11. )
  12. val_loss = np.sum(np.multiply(losses, nums)) / np.sum(nums)
  13. print(epoch, val_loss)

get_data返回训练集和验证集需要使用到的dataloaders。

  1. def get_data(train_ds, valid_ds, bs):
  2. return (
  3. DataLoader(train_ds, batch_size=bs, shuffle=True),
  4. DataLoader(valid_ds, batch_size=bs * 2),
  5. )

现在,我们只需要三行代码就可以获取数据、拟合模型了。

  1. train_dl, valid_dl = get_data(train_ds, valid_ds, bs)
  2. model, opt = get_model()
  3. fit(epochs, model, loss_func, opt, train_dl, valid_dl)

输出:

  1. 0 0.2961075816631317
  2. 1 0.28558296990394594

现在你能用这三行代码训练各种各样的模型。我们来看一下能否用它们来训练一个卷积神经网络(CNN)吧!

应用到卷积神经网络

我们现在将要创建一个包含三个卷积层的神经网络。因为前面章节中没有一个函数涉及到模型的具体形式,所以我们不需要对它们进行任何修改就可以训练一个卷积神经网络。

我们将会使用PyTorch中预先定义好的Conv2d类作为我们的卷积层。我们定义一个有三个卷积层的卷积神经网络。每个卷积层之后会执行ReLu。在最后,我们会执行一个平均池化操作。 (注意:view是PyTorch版的numpy reshape

  1. class Mnist_CNN(nn.Module):
  2. def __init__(self):
  3. super().__init__()
  4. self.conv1 = nn.Conv2d(1, 16, kernel_size=3, stride=2, padding=1)
  5. self.conv2 = nn.Conv2d(16, 16, kernel_size=3, stride=2, padding=1)
  6. self.conv3 = nn.Conv2d(16, 10, kernel_size=3, stride=2, padding=1)
  7. def forward(self, xb):
  8. xb = xb.view(-1, 1, 28, 28)
  9. xb = F.relu(self.conv1(xb))
  10. xb = F.relu(self.conv2(xb))
  11. xb = F.relu(self.conv3(xb))
  12. xb = F.avg_pool2d(xb, 4)
  13. return xb.view(-1, xb.size(1))
  14. lr = 0.1

Momentum是随机梯度下降的一个变型,它将前面步骤的更新也考虑在内,通常能够加快训练速度。

  1. model = Mnist_CNN()
  2. opt = optim.SGD(model.parameters(), lr=lr, momentum=0.9)
  3. fit(epochs, model, loss_func, opt, train_dl, valid_dl)

输出:

  1. 0 0.3829730714321136
  2. 1 0.2258522843360901

nn.Sequential

torch.nn中还有另一个类可以方便的用来简化我们的代码:Sequential。一个Sequential对象可以序列化运行它包含的模块。这是一个更简单的搭建神经网络的方式。

想要充分利用这一优势,我们要能够使用给定的函数轻松的定义一个自定义层。比如说,PyTorch中没有view层,我们需要为我们的网络定义一个。 Lambda函数将会创建一个层,并在后面使用Sequential定义神经网络的时候用到。

  1. class Lambda(nn.Module):
  2. def __init__(self, func):
  3. super().__init__()
  4. self.func = func
  5. def forward(self, x):
  6. return self.func(x)
  7. def preprocess(x):
  8. return x.view(-1, 1, 28, 28)

使用Sequential创建模型非常简单:

  1. model = nn.Sequential(
  2. Lambda(preprocess),
  3. nn.Conv2d(1, 16, kernel_size=3, stride=2, padding=1),
  4. nn.ReLU(),
  5. nn.Conv2d(16, 16, kernel_size=3, stride=2, padding=1),
  6. nn.ReLU(),
  7. nn.Conv2d(16, 10, kernel_size=3, stride=2, padding=1),
  8. nn.ReLU(),
  9. nn.AvgPool2d(4),
  10. Lambda(lambda x: x.view(x.size(0), -1)),
  11. )
  12. opt = optim.SGD(model.parameters(), lr=lr, momentum=0.9)
  13. fit(epochs, model, loss_func, opt, train_dl, valid_dl)

输出:

  1. 0 0.32739396529197695
  2. 1 0.25574398956298827

对DataLoader进行封装

我们的卷积神经网络已经非常简洁了,但是它只能运行在MNIST数据集上,原因如下:

  • 它假定输入是长度为28*28的向量
  • 它假定卷积神经网络最终输出是大小为4*4的网格(因为这是平均值池化操作时我们使用的核大小)

让我们摆脱这两种假定,这样我们的模型就可以运行在任意的2d单通道图像上。 首先,我们可以删除最初的Lambda层,并将数据预处理放在一个生成器中。

  1. def preprocess(x, y):
  2. return x.view(-1, 1, 28, 28), y
  3. class WrappedDataLoader:
  4. def __init__(self, dl, func):
  5. self.dl = dl
  6. self.func = func
  7. def __len__(self):
  8. return len(self.dl)
  9. def __iter__(self):
  10. batches = iter(self.dl)
  11. for b in batches:
  12. yield (self.func(*b))
  13. train_dl, valid_dl = get_data(train_ds, valid_ds, bs)
  14. train_dl = WrappedDataLoader(train_dl, preprocess)
  15. valid_dl = WrappedDataLoader(valid_dl, preprocess)

接下来,我们用nn.AdaptiveAvgPool2d替换nn.AvgPool2d,这个函数允许我们定义期望输出张量的大小,而不是定义已有输入的大小。 这样我们的模型就可以处理任意大小的输入了。

  1. model = nn.Sequential(
  2. nn.Conv2d(1, 16, kernel_size=3, stride=2, padding=1),
  3. nn.ReLU(),
  4. nn.Conv2d(16, 16, kernel_size=3, stride=2, padding=1),
  5. nn.ReLU(),
  6. nn.Conv2d(16, 10, kernel_size=3, stride=2, padding=1),
  7. nn.ReLU(),
  8. nn.AdaptiveAvgPool2d(1),
  9. Lambda(lambda x: x.view(x.size(0), -1)),
  10. )
  11. opt = optim.SGD(model.parameters(), lr=lr, momentum=0.9)

我们来试一下新的模型:

  1. fit(epochs, model, loss_func, opt, train_dl, valid_dl)

输出:

  1. 0 0.32888883714675904
  2. 1 0.31000419993400574

使用你的GPU

如果你有幸拥有支持CUDA的GPU(你可以租一个,大部分云服务提供商的价格使0.5$/每小时),那你可以用GPU来加速你的代码。 首先检查一下的GPU是否可以被PyTorch调用:

  1. print(torch.cuda.is_available())

输出:

  1. True

接下来,新建一个设备对象:

  1. dev = torch.device(
  2. "cuda") if torch.cuda.is_available() else torch.device("cpu")

然后更新一下preprocess函数将批运算移到GPU上计算

  1. def preprocess(x, y):
  2. return x.view(-1, 1, 28, 28).to(dev), y.to(dev)
  3. train_dl, valid_dl = get_data(train_ds, valid_ds, bs)
  4. train_dl = WrappedDataLoader(train_dl, preprocess)
  5. valid_dl = WrappedDataLoader(valid_dl, preprocess)

最后,我们可以把模型移动到GPU上。

  1. model.to(dev)
  2. opt = optim.SGD(model.parameters(), lr=lr, momentum=0.9)

你现在应该能发现模型运算变快了。

  1. fit(epochs, model, loss_func, opt, train_dl, valid_dl)

输出:

  1. 0 0.21190375366210937
  2. 1 0.18018000435829162

总结

现在我们有一个用PyTorch构建的通用数据管道和训练循环可以用来训练很多类型的模型。 想要知道训练一个模型有多么简单,可以参照mnist_sample这个例子。

当然了,你可能还想要在模型中加入很多其它的东西,比如数据扩充,超参数调整,训练监控,迁移学习等。这些特性可以在fastai库中获取到,该库使用与本教程中介绍的相同设计方法开发的,为想要扩展模型的学习者提供了合理的后续步骤。

在教程的开始部分,我们说了要通过例子对torch.nntorch.optimDatasetDataLoader进行说明。 现在我们来总结一下,我们都讲了些什么:

  • torch.nn
    • Module:创建一个可调用的,其表现类似于函数,但又可以包含状态(比如神经网络层的权重)的对象。该对象知道它包含的Parameter(s),并可以将梯度置为0,以及对梯度进行循环以更新权重等。
    • Parameter:是一个对张量的封装,它告诉Module在反向传播阶段更新权重。只有设置了requires_grad属性的张量会被更新。
    • functional:一个包含了梯度函数、损失函数等以及一些无状态的层,如卷积层和线性层的模块(通常使用F作为导入的别名)。
  • torch.optim:包含了优化器,比如在反向阶段更新Parameter中权重的SGD
  • Dataset:一个抽象接口,包含了__len____getitem__,还包含了PyTorch提供的类,如TensorDataset
  • DataLoader:接受任意的Dataset并生成一个迭代器可以批量返回数据。