Pytorch_Gradient Descent

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

# 准备训练集的数据
x_data = [1.0, 2.0, 3.0]
y_data = [2.0, 4.0, 6.0]
w = 1.0 # 设置初始权重为1

# 计算yhat ,yhat = x * w
def forward(x):
    return x * w

# 计算cost
def cost(xs, ys):
    cost = 0
    for x, y in zip(xs, ys): #zip 拼装
        y_pred = forward(x)
        cost += (y_pred - y) ** 2
    return cost / len(xs)

# 计算梯度
def gradient(xs, ys):
    grad = 0
    for x, y in zip(xs, ys):
        grad += 2 * x * (x * w - y)
    return grad / len(xs)

print('Predict (before training)', 4, forward(4))

mse_list = []
for epoch in range(100):
    cost_val = cost(x_data, y_data)
    grad_val = gradient(x_data, y_data)
    w -= 0.01 * grad_val
    mse_list.append(cost_val)
    print('epoch:', epoch, 'w = ', w, 'loss = ', cost_val)
print('predict (after training)', 4, forward(4))

# 绘图
w_list = np.arange(0,100,1)
plt.plot(w_list,mse_list)
plt.ylabel('cost')
plt.xlabel('epoch')
plt.show()

# 随机梯度下降法
# 计算yhat ,yhat = x * w
def forward(x):
    return x * w

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

# 计算梯度
def gradient(x, y):
    return 2 * x * (x * w - y)

print('Predict (before training)', 4, forward(4))

mse_list = []
for epoch in range(100):
    for x, y, in zip(x_data, y_data):
        grad = gradient(x, y)
        w = w - 0.01 * grad
        print("\tgrad:", x, y, grad)
        l = loss(x, y)
    mse_list.append(l)
    print('epoch:', epoch, 'w = ', w, 'loss = ', l)
print('predict (after training)', 4, forward(4))

# 绘图
w_list = np.arange(0,100,1)
plt.plot(w_list,mse_list)
plt.ylabel('cost')
plt.xlabel('epoch')
plt.show()