主頁 > 後端開發 > 使用cnn,bpnn,lstm實作mnist資料集的分類

使用cnn,bpnn,lstm實作mnist資料集的分類

2023-02-14 07:18:53 後端開發

1.cnn

import torch
import torch.nn as nn
import torch.nn.functional as F
import torch.optim as optim
from torchvision import datasets, transforms

# 設定亂數種子
torch.manual_seed(0)

# 超引數
EPOCH = 1  # 訓練整批資料的次數
BATCH_SIZE = 50
DOWNLOAD_MNIST = False  # 表示還沒有下載資料集,如果資料集下載好了就寫False

# 加載 MNIST 資料集
train_dataset = datasets.MNIST(
    root="./mnist",
    train=True,#True表示是訓練集
    transform=transforms.ToTensor(),
    download=False)
test_dataset = datasets.MNIST(
    root="./mnist",
    train=False,#Flase表示測驗集
    transform=transforms.ToTensor(),
    download=False)

# 將資料集放入 DataLoader 中
train_loader = torch.utils.data.DataLoader(
    dataset=train_dataset,
    batch_size=100,#每個批次讀取的資料樣本數
    shuffle=True)#是否將資料打亂,在這種情況下為True,表示每次讀取的資料是隨機的
test_loader = torch.utils.data.DataLoader(dataset=test_dataset, batch_size=100, shuffle=False)

# 為了節約時間, 我們測驗時只測驗前2000個
test_x = torch.unsqueeze(test_dataset.test_data, dim=1).type(torch.FloatTensor)[
         :2000] / 255.  # shape from (2000, 28, 28) to (2000, 1, 28, 28), value in range(0,1)
test_y = test_dataset.test_labels[:2000]

# 定義卷積神經網路模型
class CNN(nn.Module):
    def __init__(self):
        super(CNN, self).__init__()
        self.conv1 = nn.Conv2d(#輸入影像的大小為(28,28,1)
            in_channels=1,#當前輸入特征圖的個數
            out_channels=32,#輸出特征圖的個數
            kernel_size=3,#卷積核大小,在一個3*3空間里對當前輸入的特征影像進行特征提取
            stride=1,#步長:卷積視窗每隔一個單位滑動一次
            padding=1)#如果希望卷積后大小跟原來一樣,需要設定padding=(kernel_size-1)/2
        #第一層結束后影像大小為(28,28,32)32是輸出影像個數,28計算方法為(h-k+2p)/s+1=(28-3+2*1)/1 +1=28
        self.pool = nn.MaxPool2d(kernel_size=2, stride=2)#可以縮小輸入影像的尺寸,同時也可以防止過擬合
        #通過池化層之后影像大小變為(14,14,32)
        self.conv2 = nn.Conv2d(#輸入影像大小為(14,14,32)
            in_channels=32,#第一層的輸出特征圖的個數當做第二層的輸入特征圖的個數
            out_channels=64,
            kernel_size=3,
            stride=1,
            padding=1)#二層卷積之后影像大小為(14,14,64)

        self.fc = nn.Linear(64 * 7 * 7, 10)#10表示最終輸出的

    # 下面定義x的傳播路線
    def forward(self, x):
        x = self.pool(F.relu(self.conv1(x)))# x先通過conv1
        x = self.pool(F.relu(self.conv2(x)))# 再通過conv2
        x = x.view(-1, 64 * 7 * 7)
        x = self.fc(x)
        return x

# 實體化卷積神經網路模型
model = CNN()

# 定義損失函式和優化器
criterion = nn.CrossEntropyLoss()
#lr(學習率)是控制每次更新的引數的大小的超引數
optimizer = torch.optim.Adam(model.parameters(), lr=0.01)

# 訓練模型
for epoch in range(1):
    for i, (images, labels) in enumerate(train_loader):
        outputs = model(images)  # 先將資料放到cnn中計算output
        loss = criterion(outputs, labels)# 輸出和真實標簽的loss,二者位置不可顛倒
        optimizer.zero_grad()# 清除之前學到的梯度的引數
        loss.backward()  # 反向傳播,計算梯度
        optimizer.step()#應用梯度
        if i % 50 == 0:
            data_all = model(test_x)#不分開寫就會出現ValueError: too many values to unpack (expected 2)
            last_layer = data_all
            test_output = data_all
            pred_y = torch.max(test_output, 1)[1].data.numpy()
            accuracy = float((pred_y == test_y.data.numpy()).astype(int).sum()) / float(test_y.size(0))
            print('Epoch: ', epoch, '| train loss: %.4f' % loss.data.numpy(), '| test accuracy: %.4f' % accuracy)

# print 10 predictions from test data
data_all1 = model(test_x[:10])
test_output = data_all1
_ = data_all1
pred_y = torch.max(test_output, 1)[1].data.numpy()
print(pred_y, 'prediction number')
print(test_y[:10].numpy(), 'real number')

2.bpnn

import torch
import torch.nn as nn
import torch.nn.functional as F
import torch.optim as optim
from torchvision import datasets, transforms
import torchvision
DOWNLOAD_MNIST = False  # 表示還沒有下載資料集,如果資料集下載好了就寫False
BATCH_SIZE = 50
LR = 0.01  # 學習率
# 下載mnist手寫資料集
train_loader = torchvision.datasets.MNIST(
    root='./mnist/',  # 保存或提取的位置  會放在當前檔案夾中
    train=True,  # true說明是用于訓練的資料,false說明是用于測驗的資料
    transform=torchvision.transforms.ToTensor(),  # 轉換PIL.Image or numpy.ndarray

    download=DOWNLOAD_MNIST,  # 已經下載了就不需要下載了
)

test_loader = torchvision.datasets.MNIST(
    root='./mnist/',
    train=False  # 表明是測驗集
)
train_data = https://www.cnblogs.com/twq46/p/torch.utils.data.DataLoader(dataset=train_loader, batch_size=BATCH_SIZE, shuffle=True)

# 為了節約時間, 我們測驗時只測驗前2000個
test_x = torch.unsqueeze(test_loader.test_data, dim=1).type(torch.FloatTensor)[
         :2000] / 255.  # shape from (2000, 28, 28) to (2000, 1, 28, 28), value in range(0,1)
test_y = test_loader.test_labels[:2000]

# 定義模型
class BPNN(nn.Module):
    def __init__(self):
        super(BPNN, self).__init__()
        self.fc1 = nn.Linear(28 * 28, 512)#定義了一個全連接層fc1,該層的輸入是28 * 28個數字,輸出是512個數字
        self.fc2 = nn.Linear(512, 512)
        self.fc3 = nn.Linear(512, 10)

    def forward(self, x):#x是輸入的影像
        x = x.view(-1, 28 * 28)#將輸入x的形狀轉換為二維,分別是batch_size和28 * 28
        x = F.relu(self.fc1(x))#將x通過第1個全連接層fc1進行計算,并將結果通過ReLU激活函式處理
        x = F.relu(self.fc2(x))
        x = self.fc3(x)
        #Softmax函式是一種分類模型中常用的激活函式,它能將輸入資料映射到(0,1)范圍內,并且滿足所有元素的和為1
        return F.log_softmax(x, dim=1)#dim=1表示對每一行的資料進行運算

# 初始化模型
bpnn = BPNN()
print(bpnn)
# 定義損失函式和優化器
optimizer = torch.optim.Adam(bpnn.parameters(), lr=LR)  # optimize all parameters
loss_func = nn.CrossEntropyLoss()  # the target label is not one-hotted
#
# criterion = nn.NLLLoss()
# optimizer = optim.SGD(bpnn.parameters(), lr=0.01, momentum=0.5)

# 訓練模型
for epoch in range(1):
    for step, (b_x,b_y) in enumerate(train_data):
        b_x = b_x.view(-1, 28, 28)  # reshape x to (batch, time_step, input_size)

        output = bpnn(b_x)
        loss = loss_func(output, b_y)
        optimizer.zero_grad()
        loss.backward()
        optimizer.step()
        if step % 50 == 0:
            test_x = test_x.view(-1, 28, 28)
            test_output = bpnn(test_x)
            pred_y = torch.max(test_output, 1)[1].data.squeeze()
            acc = (pred_y == test_y).sum().float() / test_y.size(0)
            print('Epoch: ', epoch, '| train loss: %.4f' % loss.data.float(), 'test acc: ', acc.numpy())

test_output = bpnn(test_x[:10].view(-1, 28, 28))
pred_y = torch.max(test_output, 1)[1].data.numpy().squeeze()
print(pred_y, 'prediction number')
print(test_y[:10], 'real number')
# # 評估模型
# bpnn.eval()
# correct = 0
# with torch.no_grad():
#     for data, target in test_loader:
#         output = bpnn(data)
#         pred = output.argmax(dim=1, keepdim=True)
#         correct += pred.eq(target.view_as(pred)).sum().item()
#
# print('Test accuracy:', correct / len(test_loader.dataset))


3.lstm

import torch
from torch import nn
import torchvision.datasets as dsets
import torchvision.transforms as transforms
import matplotlib.pyplot as plt
import numpy as np

torch.manual_seed(1)  # reproducible

# Hyper Parameters
EPOCH = 1  # 訓練整批資料多少次, 為了節約時間, 我們只訓練一次
BATCH_SIZE = 64
TIME_STEP = 28  # rnn 時間步數 / 圖片高度
INPUT_SIZE = 28  # rnn 每步輸入值 / 圖片每行像素
LR = 0.01  # learning rate
DOWNLOAD_MNIST = False  # 如果你已經下載好了mnist資料就寫上 Fasle

# Mnist 手寫數字
train_data = https://www.cnblogs.com/twq46/p/dsets.MNIST(
    root='./mnist/',  # 保存或者提取位置
    train=True,  # this is training data
    transform=transforms.ToTensor(),  # 轉換 PIL.Image or numpy.ndarray 成
    # torch.FloatTensor (C x H x W), 訓練的時候 normalize 成 [0.0, 1.0] 區間
    download=DOWNLOAD_MNIST,  # 沒下載就下載, 下載了就不用再下了
)

test_data = https://www.cnblogs.com/twq46/p/dsets.MNIST(root='./mnist/', train=False)

# 批訓練 50samples, 1 channel, 28x28 (50, 1, 28, 28)
train_loader = torch.utils.data.DataLoader(dataset=train_data, batch_size=BATCH_SIZE, shuffle=True)

# 為了節約時間, 我們測驗時只測驗前2000個
test_x = torch.unsqueeze(test_data.test_data, dim=1).type(torch.FloatTensor)[
         :2000] / 255.  # shape from (2000, 28, 28) to (2000, 1, 28, 28), value in range(0,1)
test_y = test_data.test_labels[:2000]

#LSTM默認input(seq_len,batch,feature)
class Lstm(nn.Module):
    def __init__(self):
        super(Lstm, self).__init__()

        self.Lstm = nn.LSTM(  # LSTM 效果要比 nn.RNN() 好多了
            input_size=28,  # 圖片每行的資料像素點,輸入特征的大小
            hidden_size=64,  # lstm模塊的數量相當于bp網路影藏層神經元的個數
            num_layers=1,  # 隱藏層的層數
            batch_first=True,  # input & output 會是以 batch size 為第一維度的特征集 e.g. (batch, time_step, input_size)
        )

        self.out = nn.Linear(64, 10)  # 輸出層,接入線性層

    def forward(self, x):  # 必須有這個方法
        # x shape (batch, time_step, input_size)
        # r_out shape (batch, time_step, output_size)包含每個序列的輸出結果
        # h_n shape (n_layers, batch, hidden_size)只包含最后一個序列的輸出結果,LSTM 有兩個 hidden states, h_n 是分線, h_c 是主線
        # h_c shape (n_layers, batch, hidden_size)只包含最后一個序列的輸出結果
        r_out, (h_n, h_c) = self.Lstm(x, None)  # None 表示 hidden state 會用全0的 state
        # 當RNN運行結束時刻,(h_n, h_c)表示最后的一組hidden states,這里用不到

        # 選取最后一個時間點的 r_out 輸出
        # 這里 r_out[:, -1, :] 的值也是 h_n 的值
        out = self.out(r_out[:, -1, :])  # (batch_size, time step, input),這里time step選擇最后一個時刻
        # output_np = out.detach().numpy()  # 可以使用numpy的sciview監視每次結果
        return out


Lstm = Lstm()
print(Lstm)

optimizer = torch.optim.Adam(Lstm.parameters(), lr=LR)  # optimize all parameters
loss_func = nn.CrossEntropyLoss()  # the target label is not one-hotted

# training and testing
for epoch in range(EPOCH):
    for step, (x, b_y) in enumerate(train_loader):  # gives batch data
        b_x = x.view(-1, 28, 28)  # reshape x to (batch, time_step, input_size)

        output = Lstm(b_x)  # rnn output
        loss = loss_func(output, b_y)  # cross entropy loss
        optimizer.zero_grad()  # clear gradients for this training step
        loss.backward()  # backpropagation, compute gradients
        optimizer.step()  # apply gradients

        # output_np = output.detach().numpy()

        if step % 50 == 0:
            test_x = test_x.view(-1, 28, 28)
            test_output = Lstm(test_x)
            pred_y = torch.max(test_output, 1)[1].data.squeeze()
            acc = (pred_y == test_y).sum().float() / test_y.size(0)
            print('Epoch: ', epoch, '| train loss: %.4f' % loss.data.float(), 'test acc: ', acc.numpy())

test_output = Lstm(test_x[:10].view(-1, 28, 28))
pred_y = torch.max(test_output, 1)[1].data.numpy().squeeze()
print(pred_y, 'prediction number')
print(test_y[:10], 'real number')


轉載請註明出處,本文鏈接:https://www.uj5u.com/houduan/543725.html

標籤:Python

上一篇:《Terraform 101 從入門到實踐》 Terraform在公有云Azure上的應用

下一篇:多變數兩兩相互關系聯合分布圖的Python繪制

標籤雲
其他(157675) Python(38076) JavaScript(25376) Java(17977) C(15215) 區塊鏈(8255) C#(7972) AI(7469) 爪哇(7425) MySQL(7132) html(6777) 基礎類(6313) sql(6102) 熊猫(6058) PHP(5869) 数组(5741) R(5409) Linux(5327) 反应(5209) 腳本語言(PerlPython)(5129) 非技術區(4971) Android(4554) 数据框(4311) css(4259) 节点.js(4032) C語言(3288) json(3245) 列表(3129) 扑(3119) C++語言(3117) 安卓(2998) 打字稿(2995) VBA(2789) Java相關(2746) 疑難問題(2699) 细绳(2522) 單片機工控(2479) iOS(2429) ASP.NET(2402) MongoDB(2323) 麻木的(2285) 正则表达式(2254) 字典(2211) 循环(2198) 迅速(2185) 擅长(2169) 镖(2155) 功能(1967) .NET技术(1958) Web開發(1951) python-3.x(1918) HtmlCss(1915) 弹簧靴(1913) C++(1909) xml(1889) PostgreSQL(1872) .NETCore(1853) 谷歌表格(1846) Unity3D(1843) for循环(1842)

熱門瀏覽
  • 【C++】Microsoft C++、C 和匯編程式檔案

    ......

    uj5u.com 2020-09-10 00:57:23 more
  • 例外宣告

    相比于斷言適用于排除邏輯上不可能存在的狀態,例外通常是用于邏輯上可能發生的錯誤。 例外宣告 Item 1:當函式不可能拋出例外或不能接受拋出例外時,使用noexcept 理由 如果不打算拋出例外的話,程式就會認為無法處理這種錯誤,并且應當盡早終止,如此可以有效地阻止例外的傳播與擴散。 示例 //不可 ......

    uj5u.com 2020-09-10 00:57:27 more
  • Codeforces 1400E Clear the Multiset(貪心 + 分治)

    鏈接:https://codeforces.com/problemset/problem/1400/E 來源:Codeforces 思路:給你一個陣列,現在你可以進行兩種操作,操作1:將一段沒有 0 的區間進行減一的操作,操作2:將 i 位置上的元素歸零。最終問:將這個陣列的全部元素歸零后操作的最少 ......

    uj5u.com 2020-09-10 00:57:30 more
  • UVA11610 【Reverse Prime】

    本人看到此題沒有翻譯,就附帶了一個自己的翻譯版本 思考 這一題,它的第一個要求是找出所有 $7$ 位反向質數及其質因數的個數。 我們應該需要質數篩篩選1~$10^{7}$的所有數,這里就不慢慢介紹了。但是,重讀題,我們突然發現反向質數都是 $7$ 位,而將它反過來后的數字卻是 $6$ 位數,這就說明 ......

    uj5u.com 2020-09-10 00:57:36 more
  • 統計區間素數數量

    1 #pragma GCC optimize(2) 2 #include <bits/stdc++.h> 3 using namespace std; 4 bool isprime[1000000010]; 5 vector<int> prime; 6 inline int getlist(int ......

    uj5u.com 2020-09-10 00:57:47 more
  • C/C++編程筆記:C++中的 const 變數詳解,教你正確認識const用法

    1、C中的const 1、區域const變數存放在堆疊區中,會分配記憶體(也就是說可以通過地址間接修改變數的值)。測驗代碼如下: 運行結果: 2、全域const變數存放在只讀資料段(不能通過地址修改,會發生寫入錯誤), 默認為外部聯編,可以給其他源檔案使用(需要用extern關鍵字修飾) 運行結果: ......

    uj5u.com 2020-09-10 00:58:04 more
  • 【C++犯錯記錄】VS2019 MFC添加資源不懂如何修改資源宏ID

    1. 首先在資源視圖中,添加資源 2. 點擊新添加的資源,復制自動生成的ID 3. 在解決方案資源管理器中找到Resource.h檔案,編輯,使用整個專案搜索和替換的方式快速替換 宏宣告 4. Ctrl+Shift+F 全域搜索,點擊查找全部,然后逐個替換 5. 為什么使用搜索替換而不使用屬性視窗直 ......

    uj5u.com 2020-09-10 00:59:11 more
  • 【C++犯錯記錄】VS2019 MFC不懂的批量添加資源

    1. 打開資源頭檔案Resource.h,在其中預先定義好宏 ID(不清楚其實ID值應該設定多少,可以先新建一個相同的資源項,再在這個資源的ID值的基礎上遞增即可) 2. 在資源視圖中選中專案資源,按F7編輯資源檔案,按 ID 型別 相對路徑的形式添加 資源。(別忘了先把檔案拷貝到專案中的res檔案 ......

    uj5u.com 2020-09-10 01:00:19 more
  • C/C++編程筆記:關于C++的參考型別,專供新手入門使用

    今天要講的是C++中我最喜歡的一個用法——參考,也叫別名。 參考就是給一個變數名取一個變數名,方便我們間接地使用這個變數。我們可以給一個變數創建N個參考,這N + 1個變數共享了同一塊記憶體區域。(參考型別的變數會占用記憶體空間,占用的記憶體空間的大小和指標型別的大小是相同的。雖然參考是一個物件的別名,但 ......

    uj5u.com 2020-09-10 01:00:22 more
  • 【C/C++編程筆記】從頭開始學習C ++:初學者完整指南

    眾所周知,C ++的學習曲線陡峭,但是花時間學習這種語言將為您的職業帶來奇跡,并使您與其他開發人員區分開。您會更輕松地學習新語言,形成真正的解決問題的技能,并在編程的基礎上打下堅實的基礎。 C ++將幫助您養成良好的編程習慣(即清晰一致的編碼風格,在撰寫代碼時注釋代碼,并限制類內部的可見性),并且由 ......

    uj5u.com 2020-09-10 01:00:41 more
最新发布
  • Rust中的智能指標:Box<T> Rc<T> Arc<T> Cell<T> RefCell<T> Weak

    Rust中的智能指標是什么 智能指標(smart pointers)是一類資料結構,是擁有資料所有權和額外功能的指標。是指標的進一步發展 指標(pointer)是一個包含記憶體地址的變數的通用概念。這個地址參考,或 ” 指向”(points at)一些其 他資料 。參考以 & 符號為標志并借用了他們所 ......

    uj5u.com 2023-04-20 07:24:10 more
  • Java的值傳遞和參考傳遞

    值傳遞不會改變本身,參考傳遞(如果傳遞的值需要實體化到堆里)如果發生修改了會改變本身。 1.基本資料型別都是值傳遞 package com.example.basic; public class Test { public static void main(String[] args) { int ......

    uj5u.com 2023-04-20 07:24:04 more
  • [2]SpinalHDL教程——Scala簡單入門

    第一個 Scala 程式 shell里面輸入 $ scala scala> 1 + 1 res0: Int = 2 scala> println("Hello World!") Hello World! 檔案形式 object HelloWorld { /* 這是我的第一個 Scala 程式 * 以 ......

    uj5u.com 2023-04-20 07:23:58 more
  • 理解函式指標和回呼函式

    理解 函式指標 指向函式的指標。比如: 理解函式指標的偽代碼 void (*p)(int type, char *data); // 定義一個函式指標p void func(int type, char *data); // 宣告一個函式func p = func; // 將指標p指向函式func ......

    uj5u.com 2023-04-20 07:23:52 more
  • Django筆記二十五之資料庫函式之日期函式

    本文首發于公眾號:Hunter后端 原文鏈接:Django筆記二十五之資料庫函式之日期函式 日期函式主要介紹兩個大類,Extract() 和 Trunc() Extract() 函式作用是提取日期,比如我們可以提取一個日期欄位的年份,月份,日等資料 Trunc() 的作用則是截取,比如 2022-0 ......

    uj5u.com 2023-04-20 07:23:45 more
  • 一天吃透JVM面試八股文

    什么是JVM? JVM,全稱Java Virtual Machine(Java虛擬機),是通過在實際的計算機上仿真模擬各種計算機功能來實作的。由一套位元組碼指令集、一組暫存器、一個堆疊、一個垃圾回收堆和一個存盤方法域等組成。JVM屏蔽了與作業系統平臺相關的資訊,使得Java程式只需要生成在Java虛擬機 ......

    uj5u.com 2023-04-20 07:23:31 more
  • 使用Java接入小程式訂閱訊息!

    更新完微信服務號的模板訊息之后,我又趕緊把微信小程式的訂閱訊息給實作了!之前我一直以為微信小程式也是要企業才能申請,沒想到小程式個人就能申請。 訊息推送平臺🔥推送下發【郵件】【短信】【微信服務號】【微信小程式】【企業微信】【釘釘】等訊息型別。 https://gitee.com/zhongfuch ......

    uj5u.com 2023-04-20 07:22:59 more
  • java -- 緩沖流、轉換流、序列化流

    緩沖流 緩沖流, 也叫高效流, 按照資料型別分類: 位元組緩沖流:BufferedInputStream,BufferedOutputStream 字符緩沖流:BufferedReader,BufferedWriter 緩沖流的基本原理,是在創建流物件時,會創建一個內置的默認大小的緩沖區陣列,通過緩沖 ......

    uj5u.com 2023-04-20 07:22:49 more
  • Java-SpringBoot-Range請求頭設定實作視頻分段傳輸

    老實說,人太懶了,現在基本都不喜歡寫筆記了,但是網上有關Range請求頭的文章都太水了 下面是抄的一段StackOverflow的代碼...自己大修改過的,寫的注釋挺全的,應該直接看得懂,就不解釋了 寫的不好...只是希望能給視頻網站開發的新手一點點幫助吧. 業務場景:視頻分段傳輸、視頻多段傳輸(理 ......

    uj5u.com 2023-04-20 07:22:42 more
  • Windows 10開發教程_編程入門自學教程_菜鳥教程-免費教程分享

    教程簡介 Windows 10開發入門教程 - 從簡單的步驟了解Windows 10開發,從基本到高級概念,包括簡介,UWP,第一個應用程式,商店,XAML控制元件,資料系結,XAML性能,自適應設計,自適應UI,自適應代碼,檔案管理,SQLite資料庫,應用程式到應用程式通信,應用程式本地化,應用程式 ......

    uj5u.com 2023-04-20 07:22:35 more