找回密码
 立即注册
首页 业界区 安全 1 线性模型

1 线性模型

赘暨逢 6 小时前
所学习的B站视频链接 1 线性模型

目录

  • 1 模型原理

    • 1.1数据集的预处理
    • 1.2 模型的引入
    • 1.3 使用的方法

  • 2 代码编写

1 模型原理

这里我们讲解的是一个线性模型的使用方法Linear Model
在使用模型的过程中主要有四个步骤

  • dataset
  • model
  • Training
  • inferring
1.1数据集的预处理

1.png

就像上面说到的,在期末考试中已经得到了很多个数据,这些数据是学习时间和最后取得分的关系,求解并推理当\(x=4\)的时候这个分数大概是多少\(y=?\)。
2.png

将数据集交付给我们的数据模型,然后我们训练完毕后求推理的目标。
同时在学习过程中我们将已经有答案的数据称为训练集合,在训练集合中还需要分开使用测试集合。
将还未得到输出的数据称为训练模型。
因此我们将自己的数据集合在拿到的时候需要将整个数据分成两份,一份是训练使用,另一份是我们的测试模型
3.png

在处理采样数据集合的过程中我们基本上是无法求出数据的真是的分布的哈。
如果在训练的过程中在训练集合的时候,误差特别的小,甚至于将误差都已经学习进来了,这个时候这个训练可以说是过拟合了。因此我们很有可能再一次将训练集进行分开,一部分是训练集,另一部分是开发集,用来泛化模型,让我们的模型更加适合各个误差和输出
4.png

1.2 模型的引入

在这里所谓的模型就是想要找到一个比较好的函数\(y =  f(x)\),能够比较完美的适合这个输入于输出。
5.png

现在我们就使用一个线性模型来进行举例。
y上面有一个小帽子表示的是我们自己的预测模型
6.png

7.png

我们要求的就是两个参数\(w , b\)
现在我们给出咱们的输入参数和输出的训练集模型
8.png

利用上面这个玩意儿来确定\(w,b\)并尝试对各个散点进行拟合。
1.3 使用的方法

首先先随便找一组数据看看,这个假设的数据可以引出一组\(\hat{y}\),这个时候将\(\hat{y}\)和真实数值\(y\)之间的误差也就是下文的\(loss\)(单个样本)。

\[loss=(\hat{y}-y)^{2}=(x*\omega-y)^{2}\]
这个就是咱们的损失函数,要让这个拟合的越好,那么必定有损失函数的数据越小。
如果有多个样本的话咱们需要计算所有样本的平均损失。

\[$cost=\frac{1}{N}\sum_{n=1}^N(\hat{y}_n-y_n)^2$\]
因此我们真正需要的就是求出一组数据\(w,b\),使得该模型的平均损失达到最小即可。
9.png

MSE平均平方误差。
这个时候可以根据这样的公式进行损失上面的一个计算。
10.png

那我们怎么样求出这样的好数据呢,如果数据比较少,那么方法非常简单,我们可以使用暴力枚举法进行求解,也就是穷举法,求出在某一段部分中损失和w之间的关系曲线,并给出曲线最低点的数据即可。
2 代码编写

这里采用的穷举法进行实现的
具体的文件路径[穷举法的线性模型](E:\code\4.DeepLearning_basis\2 WithLiu2Da\1 线性模型.py)
在深度学习过程之中对可视化的操作和处理是十分的重要的,一定要学会进行绘图的绘制。
让模型在训练的过程之中时实绘制可视化图形,毕竟这玩意儿一开始就是一两天,一周之类的。
这个模型是最简单的$$y = wx$$
  1. import numpy as np
  2. import matplotlib.pyplot as plt
  3. x_data = [1.0 , 2.0 , 3.0 , 4.0 , 5.0]
  4. y_data = [2.0 , 4.0 , 6.0 , 8.0 , 10.0]
  5. def forward(x):
  6.     return x * w
  7. def MSE_loss(x , y):
  8.     y_hat = forward(x)
  9.     res = (y_hat -  y)**(2)
  10.     return res
  11. w_list = []
  12. mse_loss = []
  13. for w in np.arange(0.0 , 4.1 , 0.1):
  14.     w_list.append(w)
  15.     sum = 0
  16.     for x , y in zip(x_data , y_data):
  17.         sum += MSE_loss(x , y)
  18.     sum = sum/(len(x_data))
  19.     mse_loss.append(sum)
  20. plt.plot(w_list , mse_loss)
  21. plt.scatter(w_list , mse_loss , color = 'r' ,s = 20 )
  22. plt.xlabel('W')
  23. plt.ylabel('MSE')
  24. plt.show()
  25. print(w_list)
  26. print(mse_loss)
复制代码
11.png

稍微难一点点的$$y = wx+b$$
  1. import numpy as np
  2. import matplotlib.pyplot as plt
  3. from mpl_toolkits.mplot3d import Axes3D
  4. x_data = [1.0, 2.0, 3.0, 4.0, 5.0]
  5. y_data = [4.0, 6.0, 8.0, 10.0, 12.0]
  6. def forward(x, w, b):
  7.     return x * w + b
  8. def loss(x, y, w, b):
  9.     return ((forward(x, w, b) - y) * (forward(x, w, b) - y))
  10. # 生成w和b的网格
  11. w_range = np.arange(0, 4.1, 0.1)
  12. b_range = np.arange(0, 4.1, 0.1)
  13. W, B = np.meshgrid(w_range, b_range)
  14. # 计算每个(w,b)组合的MSE
  15. MSE = np.zeros_like(W)
  16. for i in range(len(w_range)):
  17.     for j in range(len(b_range)):
  18.         w_val = w_range[i]
  19.         b_val = b_range[j]
  20.         total_loss = 0
  21.         for x, y in zip(x_data, y_data):
  22.             total_loss += loss(x, y, w_val, b_val)
  23.         MSE[j, i] = total_loss / len(x_data)  # 注意索引顺序
  24. min_idx = np.unravel_index(np.argmin(MSE, axis=None), MSE.shape)
  25. min_w = w_range[min_idx[1]]
  26. min_b = b_range[min_idx[0]]
  27. min_mse = MSE[min_idx]
  28. fig = plt.figure(figsize=(12,9))
  29. ax = fig.add_subplot(111 , projection='3d')
  30. ax.scatter(min_w, min_b, min_mse, color='red', s=100,
  31.           label=f'Min MSE: {min_mse:.2f}\nw={min_w:.1f}, b={min_b:.1f}')
  32. surf = ax.plot_surface(W , B , MSE , cmap='viridis',alpha = 0.8 , linewidth = 0 , antialiased = True)
  33. ax.set_xlabel('Weight (w)', fontsize=12)
  34. ax.set_ylabel('Bias (b)', fontsize=12)
  35. ax.set_zlabel('MSE Loss', fontsize=12)
  36. ax.set_title('MSE Loss Surface for Linear Regression y = w*x + b', fontsize=14)
  37. fig.colorbar(surf, ax=ax, shrink=0.5, aspect=20, label='MSE Value')
  38. ax.legend()
  39. ax.view_init(elev=25, azim=40)            
  40. plt.tight_layout()
  41. plt.show()
  42. print(f"Minimum MSE: {min_mse:.4f} at w={min_w:.1f}, b={min_b:.1f}")
复制代码
12.png

这个是用穷举法的使用,来模拟所谓的线性模型。
以后也是画图,但是要绘制的是误差和轮数之间的关系。这一点就非常的重要了。
13.png


来源:程序园用户自行投稿发布,如果侵权,请联系站长删除
免责声明:如果侵犯了您的权益,请联系站长,我们会及时删除侵权内容,谢谢合作!
您需要登录后才可以回帖 登录 | 立即注册