Python 梯度下降法（二）：RMSProp Optimize

Python 梯度下降法（二）：RMSProp Optimize

结合第一篇文章一起看：Python 梯度下降法（一）：Gradient Descent-CSDN博客

一、数学原理

1.1 介绍

RMSProp（Root Mean Square Propagation）是一种自适应学习率优化算法，广泛用于深度学习中的梯度下降优化。它通过调整每个参数的学习率来解决传统梯度下降法中学习率固定的问题，从而加速收敛并提高性能。

RMSProp 的核心思想是对每个参数的学习率进行自适应调整。它通过维护一个指数加权移动平均（Exponential Moving Average, EMA）的梯度平方值来调整学习率：

• 对于梯度较大的参数，降低其学习率。
• 对于梯度较小的参数，增加其学习率。
  这种方法可以有效缓解梯度下降中的震荡问题，尤其是在非凸优化问题中。

1.2 公式

符号说明：

$$\begin{array}{lc}\theta & :需要优化的参数向量\\ J(\theta )& :损失函数\\ g_t& :在第t次迭代时损失函数关于\theta 的梯度，{\nabla }_{\theta }J({\theta }_t)\\ \rho & :衰减率，常用值为0.9\\ \eta & :学习率，需要手动设置\\ ϵ& :一个及小的参数，无限趋近于零，避免不会出现零(10^{-8})\\ s_t& :指数加权移动平均\end{array}$$

1. 初始化参数为：${\theta }_0,s_0=0$

2. 迭代更新，每次迭代$t$中，更新指数加权移动平均：
   $$s_t=\rho s_{t-1}+(1-\rho )g_t\odot g_t$$ $s_t$可以理解为对梯度平方的一个平滑估计，它更关注近期的梯度信息，$\rho$控制了历史信息的衰减程度。

3. 计算自适应学习率，公式为$\frac{\eta }{\sqrt{s_t+ϵ}}$，其中，分母$\sqrt{s_t+ϵ}$起到了归一化梯度的作用，使得学习率可以更具梯度的尺寸进行自适应的调整。

4. 更新参数：${\theta }_{t+1}={\theta }_t-\frac{\eta }{\sqrt{s_t+ϵ}}{g}_t$

二、代码实现

2.1 函数代码

RMSProp优化算法实现：

# 定义RMSProp优化算法
def rmsprop_optimizer(X, y, eta, num_iter=1000, epsilon=1e-5, rho=0.9):"""X: 数据 x  mxny: 数据 y  nx1eta: 学习率  num_iter: 迭代次数epsilon： 无穷小rho: 衰减率"""m, n = X.shape theta, s = np.zeros((n, 1)), np.zeros((n, 1))  # 初始参数 nx1，以及指数加权移动平均 nx1loss_ = []  # 存储损失率的变化，便于绘图for _ in range(num_iter):# 计算预测值h = np.dot(X, theta)# 计算误差error = h - yloss_.append(np.mean(error**2) / 2)# 计算梯度gradient = (1/m) * np.dot(X.T, error)s = rho * s + (1 - rho) * np.pow(gradient, 2)  # 利用广播机制来进行运算每个维度上的平滑估计theta = theta - np.multiply(eta / np.sqrt(s + epsilon), gradient)  # Hadamar productreturn theta, loss_

2.2 总代码

import numpy as np
import matplotlib.pyplot as plt
plt.rcParams['font.sans-serif'] = ['Microsoft YaHei'] # 定义RMSProp优化算法
def rmsprop_optimizer(X, y, eta, num_iter=1000, epsilon=1e-5, rho=0.9):"""X: 数据 x  mxny: 数据 y  nx1eta: 学习率  num_iter: 迭代次数epsilon： 无穷小rho: 衰减率"""m, n = X.shape theta, s = np.zeros((n, 1)), np.zeros((n, 1))  # 初始参数 nx1，以及指数加权移动平均 nx1loss_ = []  # 存储损失率的变化，便于绘图for _ in range(num_iter):# 计算预测值h = np.dot(X, theta)# 计算误差error = h - yloss_.append(np.mean(error**2) / 2)# 计算梯度gradient = (1/m) * np.dot(X.T, error)s = rho * s + (1 - rho) * np.pow(gradient, 2)  # 利用广播机制来进行运算每个维度上的平滑估计theta = theta - np.multiply(eta / np.sqrt(s + epsilon), gradient)  # Hadamar productreturn theta, loss_# 生成一些示例数据
X = 2 * np.random.rand(100, 1)
y = 4 + 3 * X + np.random.randn(100, 1)
# 添加偏置项
X_b = np.c_[np.ones((100, 1)), X]# 设置超参数
eta = 0.1# RMSProp优化算法
theta, loss_ = rmsprop_optimizer(X_b, y, eta)print("最优参数 theta:")
print(theta)
plt.plot(range(len(loss_)), loss_, label="损失函数图像")
plt.title("损失函数图像")
plt.xlabel("迭代次数")
plt.ylabel("损失值")
plt.show()

可以发现，其对于损失值的下降性能也较好，损失率也较为稳定。

三、代码优化

3.1 存在问题

• 未使用小批量数据：该代码在每次迭代时使用了全部的训练数据 X 和 y 来计算梯度，这相当于批量梯度下降的方式。在处理大规模数据集时，这种方式可能会导致计算效率低下，并且可能会陷入局部最优解。可以参考之前小批量梯度下降的代码，引入小批量数据的处理，以提高算法的效率和泛化能力。
• 缺乏数据预处理：在实际应用中，输入数据 X 可能需要进行预处理，例如归一化或标准化，以确保不同特征具有相似的尺度，从而加快算法的收敛速度。(这里不进行解决，参考特征缩放：数据归一化-CSDN博客)
• 缺乏收敛判断：代码只是简单地进行了固定次数的迭代，没有设置收敛条件。在实际应用中，可以添加收敛判断，例如当损失值的变化小于某个阈值时提前停止迭代，以节省计算资源。

这里引入Mini-batch Gradient Descent，以及收敛判断，减少计算资源

3.2 收敛判断

# 收敛判断，设定阈值，进行收敛判断
# 满足条件即停止，减少系统资源的使用
if len(loss_) > 1 and abs(loss_[-1] - loss_[-2]) < threshold:print(f"Converged at iteration {_ + 1}")break  # 注意，这里不能使用return

3.3 函数代码

# 定义RMSProp优化算法
def rmsprop_optimizer(X, y, eta, batch_size=32, num_iter=1000, epsilon=1e-5, rho=0.9, threshold=1e-3):"""X: 数据 x  mxn，可以在传入数据之前进行数据的归一化y: 数据 y  nx1eta: 学习率  batch_size: 批量数据的大小num_iter: 迭代次数epsilon： 无穷小rho: 衰减率threshold: 收敛阈值"""m, n = X.shape theta, s = np.zeros((n, 1)), np.zeros((n, 1))  # 初始参数 nx1，以及指数加权移动平均 nx1loss_ = []  # 存储损失率的变化，便于绘图num_batchs = m // batch_sizefor _ in range(num_iter):# 打乱数据集shuffled_indices = np.random.permutation(m)X_shuffled = X[shuffled_indices]y_shuffled = y[shuffled_indices]loss_temp = []  # 存储每次小批量样本生成的值for batch in range(num_batchs):# 选取小批量样本start_index = batch * batch_sizeend_index = start_index + batch_sizexi = X_shuffled[start_index:end_index]yi = y_shuffled[start_index:end_index]# 计算预测值h = np.dot(xi, theta)# 计算误差error = h - yiloss_temp.append(np.mean(error**2) / 2)# 计算梯度gradient = (1/m) * np.dot(xi.T, error)s = rho * s + (1 - rho) * np.pow(gradient, 2)  # 利用广播机制来进行运算每个维度上的平滑估计theta = theta - np.multiply(eta / np.sqrt(s + epsilon), gradient)  # Hadamar productloss_.append(np.mean(loss_temp))# 收敛判断if len(loss_) > 1 and abs(loss_[-1] - loss_[-2]) < threshold:print(f"Converged at iteration {_ + 1}")breakreturn theta, loss_

3.4 总代码

import numpy as np
import matplotlib.pyplot as plt
plt.rcParams['font.sans-serif'] = ['Microsoft YaHei'] # 定义RMSProp优化算法
def rmsprop_optimizer(X, y, eta, batch_size=32, num_iter=1000, epsilon=1e-5, rho=0.9, threshold=1e-3):"""X: 数据 x  mxn，可以在传入数据之前进行数据的归一化y: 数据 y  nx1eta: 学习率  batch_size: 批量数据的大小num_iter: 迭代次数epsilon： 无穷小rho: 衰减率threshold: 收敛阈值"""m, n = X.shape theta, s = np.zeros((n, 1)), np.zeros((n, 1))  # 初始参数 nx1，以及指数加权移动平均 nx1loss_ = []  # 存储损失率的变化，便于绘图num_batchs = m // batch_sizefor _ in range(num_iter):# 打乱数据集shuffled_indices = np.random.permutation(m)X_shuffled = X[shuffled_indices]y_shuffled = y[shuffled_indices]loss_temp = []  # 存储每次小批量样本生成的值for batch in range(num_batchs):# 选取小批量样本start_index = batch * batch_sizeend_index = start_index + batch_sizexi = X_shuffled[start_index:end_index]yi = y_shuffled[start_index:end_index]# 计算预测值h = np.dot(xi, theta)# 计算误差error = h - yiloss_temp.append(np.mean(error**2) / 2)# 计算梯度gradient = (1/m) * np.dot(xi.T, error)s = rho * s + (1 - rho) * np.pow(gradient, 2)  # 利用广播机制来进行运算每个维度上的平滑估计theta = theta - np.multiply(eta / np.sqrt(s + epsilon), gradient)  # Hadamar productloss_.append(np.mean(loss_temp))  # 使用平均值作为参考# 收敛判断if len(loss_) > 1 and abs(loss_[-1] - loss_[-2]) < threshold:print(f"Converged at iteration {_ + 1}")breakreturn theta, loss_# 生成一些示例数据
X = 2 * np.random.rand(100, 1)
y = 4 + 3 * X + np.random.randn(100, 1)
# 添加偏置项
X_b = np.c_[np.ones((100, 1)), X]# 设置超参数
eta = 0.1# RMSProp优化算法
theta, loss_ = rmsprop_optimizer(X_b, y, eta)print("最优参数 theta:")
print(theta)
plt.plot(range(len(loss_)), loss_, label="损失函数图像")
plt.title("损失函数图像")
plt.xlabel("迭代次数")
plt.ylabel("损失值")
plt.show()

四、优缺点

4.1 优点

• 自适应学习率：RMSProp 能够根据参数的梯度变化情况自适应地调整学习率。对于梯度较大的参数，学习率会自动减小；对于梯度较小的参数，学习率会相对增大。这使得算法在处理不同尺度的梯度时更加稳定，有助于加快收敛速度。
• 缓解 Adagrad 学习率衰减过快问题：与 Adagrad 算法不同，RMSProp 使用指数加权移动平均来计算梯度平方的累积值，避免了 Adagrad 中学习率单调递减且后期学习率过小的问题，使得算法在训练后期仍然能够继续更新参数。

4.2 缺点

• 对超参数敏感：RMSProp 的性能依赖于超参数$\eta$和$\rho$的选择。如果超参数设置不当，可能会导致算法收敛速度慢或者无法收敛到最优解。
• 可能陷入局部最优：和其他基于梯度的优化算法一样，RMSProp 仍然有可能陷入局部最优解，尤其是在损失函数具有复杂的地形时。