2

我正在制作一个简单的 PyTorch 神经网络来逼近 x = [0, 2pi] 上的正弦函数。这是一个简单的架构，我使用不同的深度学习库来测试我是否了解如何使用它。神经网络在未经训练时总是产生一条水平直线，而在训练时，在 y = 0 处产生一条直线。一般来说，它总是在 y =（函数的平均值）处产生一条直线。这让我相信它的前支柱部分有问题，因为边界不应该只是未经训练的直线。这是网络的代码：

class Net(nn.Module):
    def __init__(self):
      super(Net, self).__init__()
      self.model = nn.Sequential(
      nn.Linear(1, 20),
      nn.Sigmoid(),
      nn.Linear(20, 50),
      nn.Sigmoid(),
      nn.Linear(50, 50),
      nn.Sigmoid(),
      nn.Linear(50, 1)
      )

    def forward(self, x):
        x = self.model(x)
        return x

这是训练循环

def train(net, trainloader, valloader, learningrate, n_epochs):
    net = net.train()
    loss = nn.MSELoss()
    optimizer = torch.optim.SGD(net.parameters(), lr = learningrate)

    for epoch in range(n_epochs):

        for X, y in trainloader:
            X = X.reshape(-1, 1)
            y = y.view(-1, 1)
            optimizer.zero_grad()

            outputs = net(X)

            error   = loss(outputs, y)
            error.backward()
            #net.parameters()  net.parameters() * learningrate
            optimizer.step()

        total_loss = 0
        for X, y in valloader:
            X = X.reshape(-1, 1).float()
            y = y.view(-1, 1)
            outputs = net(X)
            error   = loss(outputs, y)
            total_loss += error.data

        print('Val loss for epoch', epoch, 'is', total_loss / len(valloader) )

它被称为：

net = Net()
losslist = train(net, trainloader, valloader, .0001, n_epochs = 4)

其中 trainloader 和 valloader 是训练和验证加载器。谁能帮我看看这有什么问题？我知道它不是学习率，因为它是我在其他框架中使用的，而且我知道这不是我使用 SGD 或 sigmoid 激活函数的事实，尽管我怀疑错误出在某处的激活函数中。

有谁知道如何解决这一问题？谢谢。

4

2 回答 2

1

在玩了一些超参数、修改网络和更改优化器（遵循这个优秀的配方）之后，我最终将行更改optimizer = torch.optim.SGD(net.parameters(), lr = learningrate)为optimizer = torch.optim.Adam(net.parameters())（使用了默认的优化器参数），运行 100 个 epoch 和批量大小等于 1。

使用了以下代码（仅在 CPU 上测试）：

import torch
import torch.nn as nn
from torch.utils import data
import numpy as np
import matplotlib.pyplot as plt

# for reproducibility
torch.manual_seed(0)
np.random.seed(0)

class Dataset(data.Dataset):

    def __init__(self, init, end, n):

        self.n = n
        self.x = np.random.rand(self.n, 1) * (end - init) + init
        self.y = np.sin(self.x)

    def __len__(self):

        return self.n

    def __getitem__(self, idx):

        x = self.x[idx, np.newaxis]
        y = self.y[idx, np.newaxis]

        return torch.Tensor(x), torch.Tensor(y)


class Net(nn.Module):
    def __init__(self):
        super(Net, self).__init__()
        self.model = nn.Sequential(
        nn.Linear(1, 20),
        nn.Sigmoid(),
        nn.Linear(20, 50),
        nn.Sigmoid(),
        nn.Linear(50, 50),
        nn.Sigmoid(),
        nn.Linear(50, 1)
        )

    def forward(self, x):
        x = self.model(x)
        return x

def train(net, trainloader, valloader, n_epochs):

    loss = nn.MSELoss()
    # Switch the two following lines and run the code
    # optimizer = torch.optim.SGD(net.parameters(), lr = 0.0001)
    optimizer = torch.optim.Adam(net.parameters())

    for epoch in range(n_epochs):

        net.train()
        for x, y in trainloader:
            optimizer.zero_grad()
            outputs = net(x).view(-1)
            error   = loss(outputs, y)
            error.backward()
            optimizer.step()

        net.eval()
        total_loss = 0
        for x, y in valloader:
            outputs = net(x)
            error   = loss(outputs, y)
            total_loss += error.data

        print('Val loss for epoch', epoch, 'is', total_loss / len(valloader) )    

    net.eval()

    f, (ax1, ax2) = plt.subplots(1, 2, sharey=True)

    def plot_result(ax, dataloader):
        out, xx, yy = [], [], []
        for x, y in dataloader:
            out.append(net(x))
            xx.append(x)
            yy.append(y)
        out = torch.cat(out, dim=0).detach().numpy().reshape(-1)
        xx = torch.cat(xx, dim=0).numpy().reshape(-1)
        yy = torch.cat(yy, dim=0).numpy().reshape(-1)
        ax.scatter(xx, yy, facecolor='green')
        ax.scatter(xx, out, facecolor='red')
        xx = np.linspace(0.0, 3.14159*2, 1000)
        ax.plot(xx, np.sin(xx), color='green')

    plot_result(ax1, trainloader)
    plot_result(ax2, valloader)
    plt.show()


train_dataset = Dataset(0.0, 3.14159*2, 100)
val_dataset = Dataset(0.0, 3.14159*2, 30)

params = {'batch_size': 1,
          'shuffle': True,
          'num_workers': 4}

trainloader = data.DataLoader(train_dataset, **params)
valloader = data.DataLoader(val_dataset, **params)

net = Net()
losslist = train(net, trainloader, valloader, n_epochs = 100)

Adam 优化器的结果： 在此处输入图像描述

SGD 优化器的结果： 在此处输入图像描述

于 2019-09-22T02:29:57.663 回答
1

一般来说，它总是在 y =（函数的平均值）处产生一条直线。

通常，这意味着 NN 到目前为止只成功训练了最后一层。正如 ViniciusArruda 所示，您需要对其进行更长时间的训练或进行更好的优化。

编辑：进一步解释.. 当只训练最后一层时，NN 有效地尝试在不知道输入 X 的情况下猜测输出 y。在这种情况下，它可以做出的最佳猜测是平均值。这样，它可以最大限度地减少其 MSE 损失。

于 2019-09-22T02:37:14.967 回答