视频链接:https://www.bilibili.com/video/BV1Y7411d7Ys?p=2&vd_source=b69e62de179bdd29b91104cb4fd157f8
代码:
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits import mplot3d
#y = 2 * x + 3
x_data = [1.0, 2.0, 3.0]
y_data = [5.0, 7.0, 9.0]
b = 3
def forward(x):
return x * w + b
def loss(x , y):
y_pred = forward(x)
return (y_pred - y)** 2
w_list = []
mse_list = []
b_list = []
for w in np.arange(0, 4.1, 0.1):
for b in np.arange(1, 5, 0.1):
print("w=", w)
print("b=", b)
l_sum = 0
for x_val, y_val in zip(x_data, y_data):
y_pre_val = forward(x_val)
loss_val = loss(x_val, y_val)
l_sum += loss_val
print('\t', x_val, y_val, y_pre_val, loss_val)
print('MSE:', l_sum / len(x_data))
b_list.append(b)
w_list.append(w)
mse_list.append(l_sum / len(x_data))
fig = plt.figure()
ax = plt.axes(projection='3d')
ax.plot_trisurf(w_list, b_list, mse_list)
plt.show()
效果图:
若有不对之处,还请指正!