Pytorch-神经网络

神经网络可以通过torch.nn工具包创建,在前面我们已经简单介绍了autogradnn依靠autograd定义模型并实现微分。一个nn.Module包含layers和一个forward(input)方法。
神经网络一个常见的训练过程如下:
(1)定义含有可训练参数的神经网络;
(2)遍历输入数据集;
(3)输入数据经过网络进行的相关运算(正向传播);
(4)计算损失函数(loss);
(5)梯度反向传播;
(6)更新网络权重,最简单的一个更新方法:

weight = weight - learning_rate * gradient

Define network

首先定义一个简单的神经网络:

import torch
import torch.nn as nn
import torch.nn.functional as F

class Net(nn.Module):
    def __init__(self):
        super(Net, self).__init__()
        self.conv1 = nn.Conv2d(1, 6, 3)  # input channel:1; output channel:6; kernel size:3*3
        self.conv2 = nn.Conv2d(6, 16, 3)  # input channel:6; output channel:16; kernel size:3*3
        self.fc1 = nn.Linear(16 * 6 * 6, 120)  # full connect
        self.fc2 = nn.Linear(120, 84)  # full connect
        self.fc3 = nn.Linear(84, 10)  # full connect

    def forward(self, x):
        x = F.max_pool2d(F.relu(self.conv1(x)), (2, 2))
        x = F.max_pool2d(F.relu(self.conv2(x)), (2, 2))
        x = x.view(-1, self.num_flat_features(x))
        x = F.relu(self.fc1(x))
        x = F.relu(self.fc2(x))
        x = self.fc3(x)
        return x

    def num_flat_features(self, x):
        size = x.size()[1:]
        num_features = 1
        for s in size:
            num_features *= s
        return num_features

net = Net()
print(net)

输出结果:

Net(
  (conv1): Conv2d(1, 6, kernel_size=(3, 3), stride=(1, 1))
  (conv2): Conv2d(6, 16, kernel_size=(3, 3), stride=(1, 1))
  (fc1): Linear(in_features=576, out_features=120, bias=True)
  (fc2): Linear(in_features=120, out_features=84, bias=True)
  (fc3): Linear(in_features=84, out_features=10, bias=True)
)

上面的网络结构示意图:


模型可训练参数可以通过net.parameters()查看:

params = list(net.parameters())
for i in range(len(params)):
    print(params[i].size())  # weight & bias

输出结果:

torch.Size([6, 1, 3, 3])    #conv1 weight
torch.Size([6])             #conv1 bias
torch.Size([16, 6, 3, 3])
torch.Size([16])
torch.Size([120, 576])
torch.Size([120])
torch.Size([84, 120])
torch.Size([84])
torch.Size([10, 84])
torch.Size([10])

现在向网络中输入一个随机的样本:

input = torch.randn(1, 1, 32, 32)
out = net(input)
print(out)

输出结果:

tensor([[ 0.0688, -0.0847,  0.0974, -0.1122,  0.0627,  0.0300, -0.0888, -0.0528,
          0.0811,  0.0312]], grad_fn=<AddmmBackward>)

torch.nn只支持mini-batches,比如nn.Conv2d输入的是4D Tensor (nSamples * nchannels * height * width),如果只有一个输入样本,可以使用input.unsqueeze(0)添加一个batch维度。

Loss Function

损失函数将(output, target)作为输入,计算结果用来估计output到target的距离。在nn包中有多种不同的损失函数,一个简单的损失函数是nn.MSELoss,用来计算output和target之间的均方差。
举一个例子:

target = torch.randn(10)    # random create target
target = target.view(1, -1) # reshape target
criterion = nn.MSELoss()
loss = criterion(out, target)
print(loss)

得到结果:

tensor(1.5616, grad_fn=<MseLossBackward>)

当我们调用loss.backward(),整个图会求loss的偏导,图中所有设置requires_grad=True的tensor,它们的.grad tensor都会累积梯度。

Backprop

为了反向传播误差,我们需要调用loss.backward()。我们需要清理已经存在的梯度,不然的话梯度会累积到之前的值上。
现在我们调用loss.backward(),并查看conv1的bias梯度在反向传播前后的变化:

net.zero_grad()
print(net.conv1.bias.grad)  # before backward
loss.backward()
print(net.conv1.bias.grad)  # after backward

输出结果:

None
tensor([-0.0133,  0.0076,  0.0047, -0.0117, -0.0077,  0.0037])

这里的输出结果很奇怪,官网给出的结果是中在backward之前的梯度是0,而我输出的结果是None,不知道问题出在哪里。

Update weights

最后就剩下更新权重了,最简单的一个方法就是随机梯度下降(Stochastic Gradient Descent, SGD):

weight = weight - learning_rate * gradient

一个简单的python实现:

learning_rate = 0.01
for f in net.parameters():
    f.data.sub_(f.grad.data * learning_rate)

除了上面介绍的SGD,还有Nesterov-SGD, Adam, RMSProp等等,Pytorch为了方便用户使用,提供了torch.optim包,里面有大部分常用的优化方法。
使用torch.optim的方法如下:

import torch.optim as optim

# create optimizer
optimizer = optim.SGD(net.parameters(), lr=0.01)

# in your training loop:
optimizer.zero_grad()   # zero the gradient buffers
output = net(input)
loss = criterion(output, target)
loss.backward()
optimizer.step()    # Does the update

结语

本节代码可以从github上下载,github上提供了cpp和python两个版本的入门程序,仅供大家参考。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 221,820评论 6 515
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 94,648评论 3 399
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 168,324评论 0 360
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 59,714评论 1 297
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 68,724评论 6 397
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 52,328评论 1 310
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,897评论 3 421
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,804评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 46,345评论 1 318
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 38,431评论 3 340
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,561评论 1 352
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 36,238评论 5 350
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,928评论 3 334
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,417评论 0 24
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,528评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,983评论 3 376
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,573评论 2 359