pytorch_homework_#1

 代码：

import numpy as np
import matplotlib 
import matplotlib.pyplot as plt
from scipy.interpolate import griddata

x_data = [1.0, 2.0, 3,0]
y_data = [3.0, 5.0, 7.0]

def forward(x):
  return x * w + b

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


w_list = []
b_list = []
mse_list = []

for w in np.arange(0.0, 4.1, 0.1):
  for b in np.arange(0.0, 2.0, 0.5):
    print('w=', w, 'b=', b)
    l_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)
      l_sum += loss_val
      print('\t', x_val, y_val, y_pred_val, loss_val)
    print('MSE=', l_sum/3)
    w_list.append(w)
    b_list.append(b)
    mse_list.append(l_sum/3)

X, Y= np.meshgrid(w_list, b_list)
Z = griddata((w_list, b_list), mse_list, (X, Y), method='cubic')

ax3 = plt.axes(projection='3d')
ax3.plot_surface(X, Y, Z, cmap='rainbow')
ax3.set_xlabel('w')
ax3.set_ylabel('b')
ax3.set_zlabel('MSE')
plt.show()  

结果：