2020年10月4号,国内已经5号凌晨了,依然在家学习。
今天是我写的第三个 Pytorch程序,从今天起也算是入门了。
这一次我想把之前自己手写的matlab实现的简易的传统的BP神经网络在Pytorch上重新实现一遍,想看看对比和效果差异。

这一次我设计的是一个四个层的全连接网络[784, 400, 200, 100, 10]的网络,输入层是400个节点,输出层是10个节点。其他的都是隐藏层,这里每一层节点的输出我用的是ReLu函数作为激活函数。测试数据依然是MNIST数据集。

核心步骤描述如下:

1:建立的网络模型如下:

结构十分清晰,完全就是按照上述的设计组建的模型,激活函数使用ReLu。

2:MNIST数据集采用运行阶段在网上下载的方式,如果指定目录已经存在该数据集,就会忽略掉download参数,跳过下载。

3:损失函数使用之前学习的交叉熵损失函数,梯度下降算则随机梯度下降。

4:为了我们方便观察整个训练的过程,我们在每一次迭代结束,都会现场使用模型去测试数据集上现场运行一把,看看实际的预测效果如何,分别记录下每次迭代过程中训练损失值,训练准确度,测试损失值,测试准确度,并且也是为了方便画图展示出来。

话不多说,我就直接上代码实例,代码的注释我都是用中文直接写的。

# -*- coding: utf-8 -*"""
Created on Fri Jul 27 17:47:03 2018@author: Administrator
"""
import numpy as np
import torch
from torchvision.datasets import mnist  # 导入 pytorch 内置的 mnist 数据from torch import nn
from torch.autograd import Variablefrom torch.utils.data import DataLoader
import matplotlib.pyplot as plt# Step 1:============================准备数据===================
# 定义一个对图像像素数据的标准化处理函数
# 变换到0~255的范围,在变换到0~1的范围
# 对数据进行标准化
# 对图像数据从矩阵形式变成一个 W*H的一维向量
def data_tf(img):img = np.array(img, dtype='float32') / 255img = (img - 0.5) / 0.5  # 标准化,img = img.reshape((-1,))  # 拉平img = torch.from_numpy(img)return img# 先来准备数据
# 使用内置函数下载 mnist 数据集,并且使用自定义的标准化函数对数据进行标准化
# download 参数是表明数据是要从网上下载么?如果该目录下已经存在数据集,就不会再下载了。
train_set = mnist.MNIST('./data', train = True, transform=data_tf, download = True)
test_set = mnist.MNIST('./data', train = False, transform=data_tf, download = True)
firstImg, firstImg_label = train_set[0]  # a为训练数据第一个的图像数据,a_label为训练数据第一个的标签
# 训练数据数量是60000
print(train_set)
# 测试数据数量是60000
print(test_set)
# 打印出第一个图像和其标签的值
print(firstImg.shape)
print(firstImg_label)# DataLoader本质上就是一个iterable(跟python的内置类型list等一样),并利用多进程来加速batch data的处理,使用yield来使用有限的内存
# 使用 pytorch 自带的 DataLoader 定义一个数据迭代器,也就是将数据进行排序标号,shuffle也就是打乱数据
# DataLoader是一个高效,简洁,直观的网络输入数据结构,便于使用和扩展
# 这种方式能加快数据计算速度,减少训练时间。
train_data = DataLoader(train_set, batch_size=64, shuffle=True)  # 训练数据
test_data = DataLoader(test_set, batch_size=128, shuffle=False)  # 测试数据
# 这里展示的是一个批量处理的数据,想象成之前学习的mini-batch,每次迭代处理一个小批量的数据。
# 训练数据是64个图像为一组数据,维度是[64, 784]
batch, batch_label = next(iter(train_data))
# 打印出一个批次的数图像和其标签,主要为了展示维度。
print(batch.shape)
print(batch_label.shape)# Step 2:============================定义模型===================
# 定义一个类,继承自 torch.nn.Module,torch.nn.Module是callable的类
# 在整个类里面重新定义一个标准的BP全连接神经网络,网络一共是四层,
# 层数定义:784, 400, 200, 100, 10
# 其中输入层784个节点,输出层是10个节点,分别代表10个数字,其他的层都是隐藏层。
# 我们使用了Relu的激活函数,而不是sigmoid激活函数
# 整个子类需要重写forward函数,
class BPNNModel(torch.nn.Module):def __init__(self):# 调用父类的初始化函数,必须要的super(BPNNModel, self).__init__()# 创建四个Sequential对象,Sequential是一个时序容器,将里面的小的模型按照序列建立网络self.layer1 = nn.Sequential(nn.Linear(784, 400), nn.ReLU())self.layer2 = nn.Sequential(nn.Linear(400, 200), nn.ReLU())self.layer3 = nn.Sequential(nn.Linear(200, 100), nn.ReLU())self.layer4 = nn.Sequential(nn.Linear(100, 10))def forward(self, img):# 每一个时序容器都是callable的,因此用法也是一样。img = self.layer1(img)img = self.layer2(img)img = self.layer3(img)img = self.layer4(img)return img# 创建和实例化一个整个模型类的对象
model = BPNNModel()
# 打印出整个模型
print(model)# Step 3:============================定义损失函数和优化器===================
# 定义 loss 函数,这里用的是交叉熵损失函数(Cross Entropy),这种损失函数之前博文也讲过的。
criterion = nn.CrossEntropyLoss()
# 我们优先使用随机梯度下降,lr是学习率: 0.1
optimizer = torch.optim.SGD(model.parameters(), 1e-1)# Step 4:============================开始训练网络===================
# 为了实时观测效果,我们每一次迭代完数据后都会,用模型在测试数据上跑一次,看看此时迭代中模型的效果。
# 用数组保存每一轮迭代中,训练的损失值和精确度,也是为了通过画图展示出来。
train_losses = []
train_acces = []
# 用数组保存每一轮迭代中,在测试数据上测试的损失值和精确度,也是为了通过画图展示出来。
eval_losses = []
eval_acces = []for e in range(20):# 4.1==========================训练模式==========================train_loss = 0train_acc = 0model.train()   # 将模型改为训练模式# 每次迭代都是处理一个小批量的数据,batch_size是64for im, label in train_data:im = Variable(im)label = Variable(label)# 计算前向传播,并且得到损失函数的值out = model(im)loss = criterion(out, label)# 反向传播,记得要把上一次的梯度清0,反向传播,并且step更新相应的参数。optimizer.zero_grad()loss.backward()optimizer.step()# 记录误差train_loss += loss.item()# 计算分类的准确率_, pred = out.max(1)num_correct = (pred == label).sum().item()acc = num_correct / im.shape[0]train_acc += acctrain_losses.append(train_loss / len(train_data))train_acces.append(train_acc / len(train_data))# 4.2==========================每次进行完一个训练迭代,就去测试一把看看此时的效果==========================# 在测试集上检验效果eval_loss = 0eval_acc = 0model.eval()  # 将模型改为预测模式# 每次迭代都是处理一个小批量的数据,batch_size是128for im, label in test_data:im = Variable(im)  # torch中训练需要将其封装即Variable,此处封装像素即784label = Variable(label)  # 此处为标签out = model(im)  # 经网络输出的结果loss = criterion(out, label)  # 得到误差# 记录误差eval_loss += loss.item()# 记录准确率_, pred = out.max(1)  # 得到出现最大值的位置,也就是预测得到的数即0—9num_correct = (pred == label).sum().item()  # 判断是否预测正确acc = num_correct / im.shape[0]  # 计算准确率eval_acc += acceval_losses.append(eval_loss / len(test_data))eval_acces.append(eval_acc / len(test_data))print('epoch: {}, Train Loss: {:.6f}, Train Acc: {:.6f}, Eval Loss: {:.6f}, Eval Acc: {:.6f}'.format(e, train_loss / len(train_data), train_acc / len(train_data),eval_loss / len(test_data), eval_acc / len(test_data)))plt.title('train loss')
plt.plot(np.arange(len(train_losses)), train_losses)
plt.plot(np.arange(len(train_acces)), train_acces)
plt.title('train acc')
plt.plot(np.arange(len(eval_losses)), eval_losses)
plt.title('test loss')
plt.plot(np.arange(len(eval_acces)), eval_acces)
plt.title('test acc')
plt.show()

这里有一些输出,我们解释下:

上图展示的是,原始图像数据中,训练数据有60000个,测试数据有10000个,其中第一个训练数据图像是784维度的向量,该图像代表的数字是5。经过DataLoader后,训练数据每一批量的数据是64个图像,是64*784维度的矩阵。

上图展示模型的设计样貌。都是参数模型,四层网络,一共有四层参数。

通过打印,以及结合最后看我们的画出来图像,可见此时模型的准确率已经达到了98%,比我用matlab当初设计的那个模型还要高出6个百分点,框架肯定是做了很多优化的。

《Pytorch - BP全连接神经网络模型》相关推荐

  1. ComeFuture英伽学院——2020年 全国大学生英语竞赛【C类初赛真题解析】(持续更新)

    视频:ComeFuture英伽学院--2019年 全国大学生英语竞赛[C类初赛真题解析]大小作文--详细解析 课件:[课件]2019年大学生英语竞赛C类初赛.pdf 视频:2020年全国大学生英语竞赛 ...

  2. ComeFuture英伽学院——2019年 全国大学生英语竞赛【C类初赛真题解析】大小作文——详细解析

    视频:ComeFuture英伽学院--2019年 全国大学生英语竞赛[C类初赛真题解析]大小作文--详细解析 课件:[课件]2019年大学生英语竞赛C类初赛.pdf 视频:2020年全国大学生英语竞赛 ...

  3. 信息学奥赛真题解析(玩具谜题)

    玩具谜题(2016年信息学奥赛提高组真题) 题目描述 小南有一套可爱的玩具小人, 它们各有不同的职业.有一天, 这些玩具小人把小南的眼镜藏了起来.小南发现玩具小人们围成了一个圈,它们有的面朝圈内,有的 ...

  4. 信息学奥赛之初赛 第1轮 讲解(01-08课)

    信息学奥赛之初赛讲解 01 计算机概述 系统基本结构 信息学奥赛之初赛讲解 01 计算机概述 系统基本结构_哔哩哔哩_bilibili 信息学奥赛之初赛讲解 02 软件系统 计算机语言 进制转换 信息 ...

  5. 信息学奥赛一本通习题答案(五)

    最近在给小学生做C++的入门培训,用的教程是信息学奥赛一本通,刷题网址 http://ybt.ssoier.cn:8088/index.php 现将部分习题的答案放在博客上,希望能给其他有需要的人带来 ...

  6. 信息学奥赛一本通习题答案(三)

    最近在给小学生做C++的入门培训,用的教程是信息学奥赛一本通,刷题网址 http://ybt.ssoier.cn:8088/index.php 现将部分习题的答案放在博客上,希望能给其他有需要的人带来 ...

  7. 信息学奥赛一本通 提高篇 第六部分 数学基础 相关的真题

    第1章   快速幂 1875:[13NOIP提高组]转圈游戏 信息学奥赛一本通(C++版)在线评测系统 第2 章  素数 第 3 章  约数 第 4 章  同余问题 第 5 章  矩阵乘法 第 6 章 ...

  8. 信息学奥赛一本通题目代码(非题库)

    为了完善自己学c++,很多人都去读相关文献,就比如<信息学奥赛一本通>,可又对题目无从下手,从今天开始,我将把书上的题目一 一的解析下来,可以做参考,如果有错,可以告诉我,将在下次解析里重 ...

  9. 信息学奥赛一本通(C++版) 刷题 记录

    总目录详见:https://blog.csdn.net/mrcrack/article/details/86501716 信息学奥赛一本通(C++版) 刷题 记录 http://ybt.ssoier. ...

  10. 最近公共祖先三种算法详解 + 模板题 建议新手收藏 例题: 信息学奥赛一本通 祖孙询问 距离

    首先什么是最近公共祖先?? 如图:红色节点的祖先为红色的1, 2, 3. 绿色节点的祖先为绿色的1, 2, 3, 4. 他们的最近公共祖先即他们最先相交的地方,如在上图中黄色的点就是他们的最近公共祖先 ...

最新文章

  1. mysql 对部分表binlog_MySQL抑制binlog日志中的BINLOG部分的方法
  2. 2017 Android 面试题 [ 基础与细节 ]
  3. 直线电机原理动画_最新的3D动画演示:同步电机原理及技术,一般人看不懂
  4. 手把手教你写高质量Android技术博客,画图工具,录像工具,Markdown写法
  5. JS If...Else
  6. 用viewpager实现图片轮播
  7. php封装app微信支付
  8. vofuria的开发(1)
  9. MyBatis(八)------MyBatis动态SQL
  10. 编程珠玑第一章习题解答
  11. 生物化学,材料化学必备!元素周期表
  12. 2021苹果AppleiOS开发证书申请详细图文流程
  13. mysql超卖问题处理_mysql 解决超卖问题的锁分析
  14. Airbnb暂停中国境内服务 中国民宿市场消化15万房源
  15. VMware连接U盘后无法显示U盘
  16. Forsage系统源码分析
  17. Access 的 AutoExec 宏
  18. python运维自动化脚本案例-python自动化运维脚本范例
  19. 搜狗都上市了,王小川还是单身,难道要找AI当女朋友?
  20. 阿里蚂蚁金服4面面经(已拿Offer)附答案!突如其来的意外之喜

热门文章

  1. 【转】JVM 架构解读
  2. Ado.Net读取Excel常见问题总结
  3. ext/iconv/.libs/iconv.o: In function `_php_iconv_strlen'
  4. 微软在动态语言支持上超越了Java?
  5. TypeScript入门教程 之 Let 关键字
  6. 面试精讲之面试考点及大厂真题 - 分布式专栏 23 分布式系统下分布式锁的实现
  7. 容器编排技术 -- Kubernetes kubectl patch 命令详解
  8. html段落前的空格,HTML空格:空格前后
  9. 【PHP】网站防止QQ拦截防红跳转代码
  10. 【Python】Pygame入门