刘二大人Pytorch实战系列

课程内容

课程代码

import numpy as np
import matplotlib.pyplot as plt

x_data = [1.0, 2.0, 3.0]
y_data = [2.0, 4.0, 6.0]


def forward(x):
    return x * w


def loss(x, y):
    y_pred = forward(x)
    return (y_pred - y) ** 2


w_list = []
mse_list = []
for w in np.arange(0.0, 4.1, 0.1):
    print('w=', w)
    loss_sum = 0
    for x_val, y_val in zip(x_data, y_data):
        y_pred_val = forward(x_val)
        loss_val = loss(x_val, y_val)
        loss_sum += loss_val
        print('\t', x_val, y_val, y_pred_val, loss_val)
    print('MSE=', loss_sum / 3)
    w_list.append(w)
    mse_list.append(loss_sum / 3)
plt.plot(w_list, mse_list)
plt.ylabel('Loss')
plt.xlabel('w')
plt.show()

课后作业

作业代码

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
from matplotlib import cm

x_data = [1.0, 2.0, 4.0]
y_data = [7.0, 10.0, 16.0]


def model(x):
    y_pred = w * x + b
    return y_pred


def loss(x, y):
    y_pred = model(x)
    return (y_pred - y) ** 2


w_list = []
b_list = []
mse_list = []
for w in np.arange(1.0, 4.1, 0.1):
    print('w=', w)
    loss_sum = 0
    for b in np.arange(2.0, 5.1, 0.1):
        print('b=', b)
        for x_val, y_val in zip(x_data, y_data):
            y_pred_val = model(x_val)
            loss_val = loss(x_val, y_val)
            loss_sum += loss_val
            print(x_val, y_val, y_pred_val, loss_val)
        print('MSE=', loss_sum / 3)
        mse_list.append(loss_sum / 3)
        b_list.append(b)
    w_list.append(w)
fig = plt.figure()
ax = fig.gca(projection='3d')

w = np.arange(1.0, 4.1, 0.1)  # w的值
b = np.arange(2.0, 5.1, 0.1)  # b的值
X, Y = np.meshgrid(w, b)

Z = (7-(X+Y))**2+(10-(2*X+b))**2+(16-(4*X+b))**2
surf = ax.plot_surface(X, Y, Z, cmap=cm.gist_rainbow)
fig.colorbar(surf, shrink=0.5, aspect=5)

plt.show()

说明：

1. meshgrid()函数输入是两个一维数组，输出是一个二维数组，就类似于生成底部x-y平面。
2. z的值我是通过粗暴的手段进行计算的
   输出的图为：