刘二大人PyTorch深度学习实践第2讲线性模型作业

np.meshgrid函数Numpy Meshgrid函数https://geek-docs.com/numpy/numpy-tutorials/numpy-meshgrid-function.html

zip函数python中zip()函数的用法https://blog.csdn.net/csdn15698845876/article/details/73411541

三维绘图Python三维绘图--Matplotlibhttps://blog.csdn.net/u014636245/article/details/82799573

作业

实现代码

# 实现线性模型（y=wx+b）并输出损失函数的3D图像
import numpy as np
import matplotlib.pyplot as plt

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) ** 2


# 设定w，b的范围
W, B = np.arange(0.0, 4.1, 0.1).round(1), np.arange(-2.0, 2.1, 0.1).round(1)  # round(1)保留小数点后1位
print("W=", W)
print("B=", B)
# 将w，b组合成二维矩阵
w, b = np.meshgrid(W, B)
print("w=", w)
print("b=", b)

# 计算损失函数
l_sum = 0
# zip函数
print("zip", list(zip(x_data, y_data)))
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
mse = np.round(l_sum / len(x_data), 1)
print('MSE=', mse)

# 绘制三维曲面图
ax = plt.axes(projection='3d')  # 创建三维坐标轴对象
# 还可以通过从mpl_toolkits.mplot3d导入对象Axes3D来实现
# from mpl_toolkits.mplot3d import Axes3D
# fig=plt.figure()
# ax = Axes3D(fig)
ax.plot_surface(w, b, mse)  # 绘制三维曲面
ax.set_xlabel('w')
ax.set_ylabel('b')
ax.set_zlabel('MSE')
plt.show()

输出图像