线性模型:y=wx+b
并绘制3D图
不太理解绘图那里np.array(MSE_list[0])是什么意思
程序总报错:Argument Z must be 2-dimensional. 运行不出来
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, 3.0]
y_data = [2.0, 4.0, 6.0]
#模型函数
def forward(x, w, b):
return x * w + b
#损失函数
def loss(x, y):
y_pre = forward(x, w, b)
return (y_pre - y)*(y_pre - y)
W, B = np.arange(0.0, 4.1, 0.1), np.arange(-2.0, 2.1, 0.1) # 规定 W,B 的区间
w, b = np.meshgrid(W, B, indexing='ij')
w_list = []
b_list = []
MSE_list = []
for w in np.arange(0.0, 4.1, 0.1):
for b in np.arange(-2.1, 2.1, 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, w, b)
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/3)
w_list.append(w)
b_list.append(b)
MSE_list.append(l_sum/3)
#绘制3D曲面
fig = plt.figure(figsize=(10,10), dpi=300)
ax = Axes3D(fig)
surf = ax.plot_surface(w, b, np.array(MSE_list[0]), rstride=1, cstride=1, cmap=cm.rainbow, linewidth=0, antialiased=False)
fig.colorbar(surf, shrink=0.5, aspect=5)
ax.set_xlabel("weight")
ax.set_ylabel("bias")
ax.set_zlabel("loss")
plt.show()