传统机器学习算法依赖人工设计特征、提取特征,而深度学习依赖算法自动提取特征。深度学习模仿人类大脑的运行方式,从大量数据中学习特征,这也是深度学习被看做黑盒子、可解释性差的原因。
随着算力的提升,深度学习可以处理图像,文本,音频,视频等各种内容,主要应用领域有:
我们要学习的深度学习(Deep Learning)是神经网络的一个子领域,主要关注更深层次的神经网络结构,也就是深层神经网络(Deep Neural Networks,DNNs)。所以,我们需要先搞清楚什么是神经网络!
神经网络(Neural Networks)是一种模拟人脑神经元网络结构的计算模型,用于处理复杂的模式识别、分类和预测等任务。生物神经元如下图:
生物学:
人脑可以看做是一个生物神经网络,由众多的神经元连接而成
人脑神经元处理信息的过程:
神经网络由多个互相连接的节点(即人工神经元)组成。
人工神经元(Artificial Neuron)是神经网络的基本构建单元,模仿了生物神经元的工作原理。其核心功能是接收输入信号,经过加权求和和非线性激活函数处理后,输出结果。
人工神经元接受多个输入信息,对它们进行加权求和,再经过激活函数处理,最后将这个结果输出。
如果有 n 个输入 x 1 , x 2 , … , x n x_1, x_2, \ldots, x_n x1,x2,…,xn,权重分别为 w 1 , w 2 , … , w n w_1, w_2, \ldots, w_n w1,w2,…,wn,偏置为 b b b,则神经元的输出 y y y 表示为:
z = ∑ i = 1 n w i ⋅ x i + b y = σ ( z ) z=\sum_{i=1}^nw_i\cdot x_i+b \\ y=\sigma(z) z=i=1∑nwi⋅xi+by=σ(z)
其中, σ ( z ) \sigma(z) σ(z) 是激活函数。
例如:
线性回归:
y = ∑ i = 1 n w i ⋅ x i + b y=\sum_{i=1}^nw_i\cdot x_i+b \\ y=i=1∑nwi⋅xi+b
线性回归不需要激活函数
逻辑回归:
z = ∑ i = 1 n w i ⋅ x i + b y = σ ( z ) = s i g m o i d ( z ) = 1 1 + e − z z=\sum_{i=1}^nw_i\cdot x_i+b \\ y=\sigma(z)=sigmoid(z)=\frac{1}{1+e^{-z}} z=i=1∑nwi⋅xi+by=σ(z)=sigmoid(z)=1+e−z1
人工神经元和生物神经元对比如下表:
生物神经元 | 人工神经元 |
---|---|
细胞核 | 节点 (加权求和 + 激活函数) |
树突 | 输入 |
轴突 | 带权重的连接 |
突触 | 输出 |
神经网络是由大量人工神经元按层次结构连接而成的计算模型。每一层神经元的输出作为下一层的输入,最终得到网络的输出。
神经网络有下面三个基础层(Layer)构建而成:
输入层(Input): 神经网络的第一层,负责接收外部数据,不进行计算。
隐藏层(Hidden): 位于输入层和输出层之间,进行特征提取和转换。隐藏层一般有多层,每一层有多个神经元。
输出层(Output): 网络的最后一层,产生最终的预测结果或分类结果
我们使用多个神经元来构建神经网络,相邻层之间的神经元相互连接,并给每一个连接分配一个权重,经典如下:
注意:同一层的各个神经元之间是没有连接的。
前馈神经网络(Feedforward Neural Network,FNN)是一种最基本的神经网络结构,其特点是信息从输入层经过隐藏层单向传递到输出层,没有反馈或循环连接。
全连接神经网络(Fully Connected Neural Network,FCNN)是前馈神经网络的一种,每一层的神经元与上一层的所有神经元全连接,常用于图像分类、文本分类等任务。
外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传
如上图,网络中每个神经元:
z 1 = x 1 ∗ w 1 + x 2 ∗ w 2 + b 1 z 2 = x 1 ∗ w 1 + x 2 ∗ w 2 + b 2 z 3 = x 1 ∗ w 1 + x 2 ∗ w 2 + b 3 z_1 = x_1*w_1 + x_2*w_2+b_1 \\ z_2 = x_1*w_1 + x_2*w_2+b_2 \\ z_3 = x_1*w_1 + x_2*w_2+b_3 z1=x1∗w1+x2∗w2+b1z2=x1∗w1+x2∗w2+b2z3=x1∗w1+x2∗w2+b3
说明:三个等式中的w1和w2在这里只是为了方便表示对应x1和x2的权重,实际三个等式中的w值是不同的。
向量x为: [ x 1 , x 2 ] [x_1,x_2] [x1,x2]
向量w: ( w 1 , w 2 w 1 , w 2 w 1 , w 2 ) \begin{pmatrix}w_1,w_2\\w_1,w_2\\w_1,w_2 \end{pmatrix} w1,w2w1,w2w1,w2 ,其形状为(3,2),3是神经元节点个数,2是向量x的个数
向量z: [ z 1 , z 2 , z 3 ] [z_1,z_2,z_3] [z1,z2,z3]
向量b: [ b 1 , b 2 , b 3 ] [b_1,b_2,b_3] [b1,b2,b3]
所以用向量表示为:
z = ( z 1 , z 2 , z 3 ) = ( x 1 , x 2 ) ( w 1 , w 1 , w 1 w 2 , w 2 , w 2 ) + ( b 1 , b 2 , b 3 ) = ( x 1 , x 2 ) ( w 1 , w 2 w 1 , w 2 w 1 , w 2 ) T + ( b 1 , b 2 , b 3 ) = x w T + b z = \begin{pmatrix}z_1,z_2,z_3 \end{pmatrix}=\begin{pmatrix}x_1,x_2 \end{pmatrix}\begin{pmatrix}w_1,w_1,w_1\\w_2,w_2,w_2\end{pmatrix}+\begin{pmatrix}b_1,b_2,b_3 \end{pmatrix}=\begin{pmatrix}x_1,x_2 \end{pmatrix}\begin{pmatrix}w_1,w_2\\w_1,w_2\\w_1,w_2 \end{pmatrix}^T + \begin{pmatrix}b_1,b_2,b_3 \end{pmatrix}=xw^T+b z=(z1,z2,z3)=(x1,x2)(w1,w1,w1w2,w2,w2)+(b1,b2,b3)=(x1,x2) w1,w2w1,w2w1,w2 T+(b1,b2,b3)=xwT+b
创建一个最基本的全连接神经网络(也称为多层感知机,MLP)通常需要以下步骤和方法:
示例:创建一个全连接神经网络,主要步骤包括:
你可以根据实际任务调整网络结构、损失函数和优化器等。
import torch
from torch import nn
from torch import optim
from torch.nn import functional as F
torch.manual_seed(42)
# 定义全连接神经网络模型
class MyFcnn(nn.Module):
def __init__(self, input_size):
# 父类初始化
super(MyFcnn, self).__init__()
# 定义线性层1
self.fc1 = nn.Linear(input_size, 64)
# 初始化w权重
nn.init.kaiming_uniform_(self.fc1.weight)
# 初始化b偏置
nn.init.zeros_(self.fc1.bias)
# 定义线性层2,输入要和第一层的输出一致
self.fc2 = nn.Linear(64, 32)
# 初始化w权重
nn.init.kaiming_uniform_(self.fc2.weight)
# 初始化b偏置
nn.init.zeros_(self.fc2.bias)
# 定义线性层3,输入要和第二层的输出一致
self.fc3 = nn.Linear(32, 1)
# 初始化w权重
nn.init.xavier_uniform_(self.fc3.weight)
# 初始化b偏置
nn.init.zeros_(self.fc3.bias)
def forward(self, x):
x = F.relu(self.fc1(x))
x = F.relu(self.fc2(x))
x = F.sigmoid(self.fc3(x))
return x
# 初始化模型
input_size = 10
model = MyFcnn(input_size)
print(model)
# 定义损失函数
criterion = nn.BCELoss()
# 定义优化器
optimizer = optim.Adam(model.parameters(), lr=0.01)
# 训练模型
def train():
model.train()
# 数据准备
# 100个样本,每个样本10个特征
x = torch.randn(100, input_size)
# 随机生成二分类标签
y = torch.randint(0, 2, (100, 1)).float()
# 迭代次数
epochs = 10
for epoch in range(epochs):
# 前向传播
y_pred = model(x)
# 计算损失
loss = criterion(y_pred, y)
# 梯度清零
optimizer.zero_grad()
# 反向传播
loss.backward()
# 更新梯度
optimizer.step()
# 打印训练信息
print(f"Epoch [{epoch + 1}/{epochs}], Loss: {loss.item():.4f}")
# 模型验证
X_test = torch.randn(20, input_size) # 20个测试样本
y_test = torch.randint(0, 2, (20, 1)).float() # 测试标签
def eval():
model.eval()
with torch.no_grad():
y_pred = model(X_test)
y_pred = (y_pred > 0.5).float()
accuracy = (y_pred == y_test).float().mean().item()
print(f"Test Accuracy: {accuracy:.4f}")
if __name__ == '__main__':
train()
eval()
激活函数的作用是在隐藏层引入非线性,使得神经网络能够学习和表示复杂的函数关系,使网络具备非线性能力,增强其表达能力。
通过认识线性和非线性的基础概念,深刻理解激活函数存在的价值。
如果在隐藏层不使用激活函数,那么整个神经网络会表现为一个线性模型。我们可以通过数学推导来展示这一点。
假设:
一层网络的情况
对于单层网络(输入层到输出层),如果没有激活函数,输出 a ( 1 ) \mathbf{a}^{(1)} a(1) 可以表示为:
a ( 1 ) = W ( 1 ) x + b ( 1 ) \mathbf{a}^{(1)} = \mathbf{W}^{(1)} \mathbf{x} + \mathbf{b}^{(1)} a(1)=W(1)x+b(1)
两层网络的情况
假设我们有两层网络,且每层都没有激活函数,则:
将 a ( 1 ) \mathbf{a}^{(1)} a(1)代入到 a ( 2 ) \mathbf{a}^{(2)} a(2)中,可以得到:
a ( 2 ) = W ( 2 ) ( W ( 1 ) x + b ( 1 ) ) + b ( 2 ) \mathbf{a}^{(2)} = \mathbf{W}^{(2)} (\mathbf{W}^{(1)} \mathbf{x} + \mathbf{b}^{(1)}) + \mathbf{b}^{(2)} a(2)=W(2)(W(1)x+b(1))+b(2)
a ( 2 ) = W ( 2 ) W ( 1 ) x + W ( 2 ) b ( 1 ) + b ( 2 ) \mathbf{a}^{(2)} = \mathbf{W}^{(2)} \mathbf{W}^{(1)} \mathbf{x} + \mathbf{W}^{(2)} \mathbf{b}^{(1)} + \mathbf{b}^{(2)} a(2)=W(2)W(1)x+W(2)b(1)+b(2)
我们可以看到,输出 a ( 2 ) \mathbf{a}^{(2)} a(2)是输入 x \mathbf{x} x的线性变换,因为: a ( 2 ) = W ′ x + b ′ \mathbf{a}^{(2)} = \mathbf{W}' \mathbf{x} + \mathbf{b}' a(2)=W′x+b′
其中 W ′ = W ( 2 ) W ( 1 ) \mathbf{W}' = \mathbf{W}^{(2)} \mathbf{W}^{(1)} W′=W(2)W(1), b ′ = W ( 2 ) b ( 1 ) + b ( 2 ) \mathbf{b}' = \mathbf{W}^{(2)} \mathbf{b}^{(1)} + \mathbf{b}^{(2)} b′=W(2)b(1)+b(2)。
多层网络的情况
如果有 L L L层,每层都没有激活函数,则第 l l l层的输出为: a ( l ) = W ( l ) a ( l − 1 ) + b ( l ) \mathbf{a}^{(l)} = \mathbf{W}^{(l)} \mathbf{a}^{(l-1)} + \mathbf{b}^{(l)} a(l)=W(l)a(l−1)+b(l)
通过递归代入,可以得到:
a ( L ) = W ( L ) W ( L − 1 ) ⋯ W ( 1 ) x + W ( L ) W ( L − 1 ) ⋯ W ( 2 ) b ( 1 ) + W ( L ) W ( L − 1 ) ⋯ W ( 3 ) b ( 2 ) + ⋯ + b ( L ) \mathbf{a}^{(L)} = \mathbf{W}^{(L)} \mathbf{W}^{(L-1)} \cdots \mathbf{W}^{(1)} \mathbf{x} + \mathbf{W}^{(L)} \mathbf{W}^{(L-1)} \cdots \mathbf{W}^{(2)} \mathbf{b}^{(1)} + \mathbf{W}^{(L)} \mathbf{W}^{(L-1)} \cdots \mathbf{W}^{(3)} \mathbf{b}^{(2)} + \cdots + \mathbf{b}^{(L)} a(L)=W(L)W(L−1)⋯W(1)x+W(L)W(L−1)⋯W(2)b(1)+W(L)W(L−1)⋯W(3)b(2)+⋯+b(L)
表达式可简化为:
a ( L ) = W ′ ′ x + b ′ ′ \mathbf{a}^{(L)} = \mathbf{W}'' \mathbf{x} + \mathbf{b}'' a(L)=W′′x+b′′
其中, W ′ ′ \mathbf{W}'' W′′ 是所有权重矩阵的乘积, b ′ ′ \mathbf{b}'' b′′是所有偏置项的线性组合。
如此可以看得出来,无论网络多少层,意味着:
整个网络就是线性模型,无法捕捉数据中的非线性关系。
激活函数是引入非线性特性、使神经网络能够处理复杂问题的关键。
我们可以通过可视化的方式去理解非线性的拟合能力:https://playground.tensorflow.org/
激活函数通过引入非线性来增强神经网络的表达能力,对于解决线性模型的局限性至关重要。由于反向传播算法(BP)用于更新网络参数,因此激活函数必须是可微的,也就是说能够求导的。
Sigmoid激活函数是一种常见的非线性激活函数,特别是在早期神经网络中应用广泛。它将输入映射到0到1之间的值,因此非常适合处理概率问题。
Sigmoid函数的数学表达式为:
f ( x ) = σ ( x ) = 1 1 + e − x f(x) = \sigma(x) = \frac{1}{1 + e^{-x}} f(x)=σ(x)=1+e−x1
其中, e e e 是自然常数(约等于2.718), x x x 是输入。
将任意实数输入映射到 (0, 1)之间,因此非常适合处理概率场景。
sigmoid函数一般只用于二分类的输出层。
微分性质: 导数计算比较方便,可以用自身表达式来表示:
σ ′ ( x ) = σ ( x ) ⋅ ( 1 − σ ( x ) ) \sigma'(x)=\sigma(x)\cdot(1-\sigma(x)) σ′(x)=σ(x)⋅(1−σ(x))
通过代码实现函数和导函数绘制:
import torch
import matplotlib.pyplot as plt
# plt支持中文
plt.rcParams["font.sans-serif"] = ["SimHei"]
plt.rcParams["axes.unicode_minus"] = False
def test001():
# 一行两列绘制图像
_, ax = plt.subplots(1, 2)
# 绘制函数图像
x = torch.linspace(-10, 10, 100)
y = torch.sigmoid(x)
# 网格
ax[0].grid(True)
ax[0].set_title("sigmoid 函数曲线图")
ax[0].set_xlabel("x")
ax[0].set_ylabel("y")
# 在第一行第一列绘制sigmoid函数曲线图
ax[0].plot(x, y)
# 绘制sigmoid导数曲线图
x = torch.linspace(-10, 10, 100, requires_grad=True)
# y = torch.sigmoid(x) * (1 - torch.sigmoid(x))
# 自动求导
torch.sigmoid(x).sum().backward()
ax[1].grid(True)
ax[1].set_title("sigmoid 函数导数曲线图", color="red")
ax[1].set_xlabel("x")
ax[1].set_ylabel("y")
# ax[1].plot(x.detach().numpy(), y.detach())
# 用自动求导的结果绘制曲线图
ax[1].plot(x.detach().numpy(), x.grad.detach().numpy())
# 设置曲线颜色
ax[1].lines[0].set_color("red")
plt.show()
if __name__ == "__main__":
test001()
tanh(双曲正切)是一种常见的非线性激活函数,常用于神经网络的隐藏层。tanh 函数也是一种S形曲线,输出范围为 ( − 1 , 1 ) (−1,1) (−1,1)。
tanh数学表达式为:
t a n h ( x ) = e x − e − x e x + e − x {tanh}(x) = \frac{e^x - e^{-x}}{e^x + e^{-x}} tanh(x)=ex+e−xex−e−x
输出范围: 将输入映射到 ( − 1 , 1 ) (-1, 1) (−1,1)之间,因此输出是零中心的。相比于Sigmoid函数,这种零中心化的输出有助于加速收敛。
对称性: Tanh函数是关于原点对称的奇函数,因此在输入为0时,输出也为0。这种对称性有助于在训练神经网络时使数据更平衡。
平滑性: Tanh函数在整个输入范围内都是连续且可微的,这使其非常适合于使用梯度下降法进行优化。
d d x tanh ( x ) = 1 − tanh 2 ( x ) \frac{d}{dx} \text{tanh}(x) = 1 - \text{tanh}^2(x) dxdtanh(x)=1−tanh2(x)
绘制代码:
import torch
import matplotlib.pyplot as plt
# plt支持中文
plt.rcParams["font.sans-serif"] = ["SimHei"]
plt.rcParams["axes.unicode_minus"] = False
def test001():
# 一行两列绘制图像
_, ax = plt.subplots(1, 2)
# 绘制函数图像
x = torch.linspace(-10, 10, 100)
y = torch.tanh(x)
# 网格
ax[0].grid(True)
ax[0].set_title("tanh 函数曲线图")
ax[0].set_xlabel("x")
ax[0].set_ylabel("y")
# 在第一行第一列绘制tanh函数曲线图
ax[0].plot(x, y)
# 绘制tanh导数曲线图
x = torch.linspace(-10, 10, 100, requires_grad=True)
# y = torch.tanh(x) * (1 - torch.tanh(x))
# 自动求导:需要标量才能反向传播
torch.tanh(x).sum().backward()
ax[1].grid(True)
ax[1].set_title("tanh 函数导数曲线图", color="red")
ax[1].set_xlabel("x")
ax[1].set_ylabel("x.grad")
# ax[1].plot(x.detach().numpy(), y.detach())
# 用自动求导的结果绘制曲线图
ax[1].plot(x.detach().numpy(), x.grad.detach().numpy())
# 设置曲线颜色
ax[1].lines[0].set_color("red")
plt.show()
if __name__ == "__main__":
test001()
ReLU(Rectified Linear Unit)是深度学习中最常用的激活函数之一,它的全称是修正线性单元。ReLU 激活函数的定义非常简单,但在实践中效果非常好。
ReLU 函数定义如下:
ReLU ( x ) = max ( 0 , x ) \text{ReLU}(x) = \max(0, x) ReLU(x)=max(0,x)
即 R e L U ReLU ReLU对输入 x x x进行非线性变换:
∙ 当 x > 0 时,ReLU ( x ) = x ∙ 当 x ≤ 0 时,ReLU ( x ) = 0 \bullet\quad\text{当 }x>0\text{ 时,ReLU}(x)=x\text{}\\\bullet\quad\text{当 }x\leq0\text{ 时,ReLU}(x)=0\text{} ∙当 x>0 时,ReLU(x)=x∙当 x≤0 时,ReLU(x)=0
计算简单:ReLU 的计算非常简单,只需要对输入进行一次比较运算,这在实际应用中大大加速了神经网络的训练。
ReLU 函数的导数是分段函数:
ReLU ′ ( x ) = { 1 , if x > 0 0 , if x ≤ 0 \text{ReLU}'(x)=\begin{cases}1,&\text{if } x>0\\0,&\text{if }x\leq0\end{cases} ReLU′(x)={1,0,if x>0if x≤0
缓解梯度消失问题:相比于 Sigmoid 和 Tanh 激活函数,ReLU 在正半区的导数恒为 1,这使得深度神经网络在训练过程中可以更好地传播梯度,不存在饱和问题。
稀疏激活:ReLU在输入小于等于 0 时输出为 0,这使得 ReLU 可以在神经网络中引入稀疏性(即一些神经元不被激活),这种稀疏性可以减少网络中的冗余信息,提高网络的效率和泛化能力。
神经元死亡:由于 R e L U ReLU ReLU在 x ≤ 0 x≤0 x≤0时输出为 0 0 0,如果某个神经元输入值是负,那么该神经元将永远不再激活,成为“死亡”神经元。随着训练的进行,网络中可能会出现大量死亡神经元,从而会降低模型的表达能力。
参考代码如下:
import torch
import torch.nn.functional as F
import matplotlib.pyplot as plt
# 中文问题
plt.rcParams["font.sans-serif"] = ["SimHei"]
plt.rcParams["axes.unicode_minus"] = False
def test006():
# 输入数据x
x = torch.linspace(-20, 20, 1000)
y = F.relu(x)
# 绘制一行2列
_, ax = plt.subplots(1, 2)
ax[0].plot(x.numpy(), y.numpy())
# 显示坐标格子
ax[0].grid()
ax[0].set_title("relu 激活函数")
ax[0].set_xlabel("x")
ax[0].set_ylabel("y")
# 绘制导数函数
x = torch.linspace(-20, 20, 1000, requires_grad=True)
F.relu(x).sum().backward()
ax[1].plot(x.detach().numpy(), x.grad.numpy())
ax[1].grid()
ax[1].set_title("relu 激活函数导数", color="red")
# 设置绘制线色颜色
ax[1].lines[0].set_color("red")
ax[1].set_xlabel("x")
ax[1].set_ylabel("x.grad")
plt.show()
if __name__ == "__main__":
test006()
Leaky ReLU是一种对 ReLU 函数的改进,旨在解决 ReLU 的一些缺点,特别是Dying ReLU 问题。Leaky ReLU 通过在输入为负时引入一个小的负斜率来改善这一问题。
Leaky ReLU 函数的定义如下:
Leaky ReLU ( x ) = { x , if x > 0 α x , if x ≤ 0 \text{Leaky ReLU}(x)=\begin{cases}x,&\text{if } x>0\\\alpha x,&\text{if } x\leq0\end{cases} Leaky ReLU(x)={x,αx,if x>0if x≤0
其中, α \alpha α 是一个非常小的常数(如 0.01),它控制负半轴的斜率。这个常数 α \alpha α是一个超参数,可以在训练过程中可自行进行调整。
参考代码:
import torch
import torch.nn.functional as F
import matplotlib.pyplot as plt
# 中文设置
plt.rcParams["font.sans-serif"] = ["SimHei"]
plt.rcParams["axes.unicode_minus"] = False
def test006():
x = torch.linspace(-5, 5, 200)
# 设置leaky_relu的负斜率超参数
slope = 0.03
y = F.leaky_relu(x, slope)
# 一行两列
_, ax = plt.subplots(1, 2)
# 开始绘制函数曲线图
ax[0].plot(x, y)
ax[0].set_title("Leaky ReLU 函数曲线图")
ax[0].set_xlabel("x")
ax[0].set_ylabel("y")
ax[0].grid(True)
# 绘制leaky_relu的梯度曲线图
x = torch.linspace(-5, 5, 200, requires_grad=True)
F.leaky_relu(x, slope).sum().backward()
ax[1].plot(x.detach().numpy(), x.grad)
ax[1].set_title("Leaky ReLU 梯度曲线图", color="red")
ax[1].set_xlabel("x")
ax[1].set_ylabel("x.grad")
ax[1].grid(True)
# 设置线的颜色
ax[1].lines[0].set_color("red")
plt.show()
if __name__ == "__main__":
test006()
Softmax激活函数通常用于分类问题的输出层,它能够将网络的输出转换为概率分布,使得输出的各个类别的概率之和为 1。Softmax 特别适合用于多分类问题。
假设神经网络的输出层有 n n n个节点,每个节点的输入为 z i z_i zi,则 Softmax 函数的定义如下:
S o f t m a x ( z i ) = e z i ∑ j = 1 n e z j \mathrm{Softmax}(z_i)=\frac{e^{z_i}}{\sum_{j=1}^ne^{z_j}} Softmax(zi)=∑j=1nezjezi
给定输入向量 z = [ z 1 , z 2 , … , z n ] z=[z_1,z_2,…,z_n] z=[z1,z2,…,zn]
1.指数变换:对每个 z i z_i zi进行指数变换,得到 t = [ e z 1 , e z 2 , . . . , e z n ] t = [e^{z_1},e^{z_2},...,e^{z_n}] t=[ez1,ez2,...,ezn],使z的取值区间从 ( − ∞ , + ∞ ) (-\infty,+\infty) (−∞,+∞)变为 ( 0 , + ∞ ) (0,+\infty) (0,+∞)
2.将所有指数变换后的值求和,得到 s = e z 1 + e z 2 + . . . + e z n = Σ j = 1 n e z j s = e^{z_1} + e^{z_2} + ... + e^{z_n} = \Sigma_{j=1}^ne^{z_j} s=ez1+ez2+...+ezn=Σj=1nezj
3.将t中每个 e z i e^{z_i} ezi除以归一化因子s,得到概率分布:
s o f t m a x ( z ) = [ e z 1 s , e z 2 s , . . . , e z n s ] = [ e z 1 Σ j = 1 n e z j , e z 2 Σ j = 1 n e z j , . . . , e z n Σ j = 1 n e z j ] softmax(z) =[\frac{e^{z_1}}{s},\frac{e^{z_2}}{s},...,\frac{e^{z_n}}{s}]=[\frac{e^{z_1}}{\Sigma_{j=1}^ne^{z_j}},\frac{e^{z_2}}{\Sigma_{j=1}^ne^{z_j}},...,\frac{e^{z_n}}{\Sigma_{j=1}^ne^{z_j}}] softmax(z)=[sez1,sez2,...,sezn]=[Σj=1nezjez1,Σj=1nezjez2,...,Σj=1nezjezn]
即:
S o f t m a x ( z i ) = e z i ∑ j = 1 n e z j \mathrm{Softmax}(z_i)=\frac{e^{z_i}}{\sum_{j=1}^ne^{z_j}} Softmax(zi)=∑j=1nezjezi
从上述公式可以看出:
每个输出值在 (0,1)之间
Softmax()对向量的值做了改变,但其位置不变
所有输出值之和为1,即
s u m ( s o f t m a x ( z ) ) = e z 1 s + e z 2 s + . . . + e z n s = s s = 1 sum(softmax(z)) =\frac{e^{z_1}}{s}+\frac{e^{z_2}}{s}+...+\frac{e^{z_n}}{s}=\frac{s}{s}=1 sum(softmax(z))=sez1+sez2+...+sezn=ss=1
将输出转化为概率:通过 S o f t m a x Softmax Softmax,可以将网络的原始输出转化为各个类别的概率,从而可以根据这些概率进行分类决策。
概率分布: S o f t m a x Softmax Softmax的输出是一个概率分布,即每个输出值 Softmax ( z i ) \text{Softmax}(z_i) Softmax(zi)都是一个介于 0 0 0和 1 1 1之间的数,并且所有输出值的和为 1:
∑ i = 1 n Softmax ( z i ) = 1 \sum_{i=1}^n\text{Softmax}(z_i)=1 i=1∑nSoftmax(zi)=1
突出差异: S o f t m a x Softmax Softmax会放大差异,使得概率最大的类别的输出值更接近 1 1 1,而其他类别更接近 0 0 0。
在实际应用中, S o f t m a x Softmax Softmax常与交叉熵损失函数Cross-Entropy Loss结合使用,用于多分类问题。在反向传播中, S o f t m a x Softmax Softmax的导数计算是必需的。
设 p i = S o f t m a x ( z i ) ,则对于 z i 的导数为: ∙ 当 i = j 时: ∂ p i ∂ z i = e z i ( Σ j = 1 n e z j ) − e z i e z i ( Σ j = 1 n e z j ) 2 = p i ( 1 − p i ) ∙ 当 i ≠ j 时 : ∂ p i ∂ z j = 0 ( Σ j = 1 n e z j ) − e z i e z j ( Σ j = 1 n e z j ) 2 = − p i p j \begin{aligned} &\text{设 }p_i=\mathrm{Softmax}(z_i)\text{,则对于 }z_i\text{ 的导数为:} \\ &\bullet\text{ 当 }i=j\text{ 时:} \\ &&&\frac{\partial p_i}{\partial z_i}=\frac{e^{z_i}(\Sigma_{j=1}^ne^{z_j})-e^{z_i}e^{z_i}}{(\Sigma_{j=1}^ne^{z_j})^2}=p_i(1-p_i) \\ & \bullet\text{ 当 }i\neq j\text{ 时}: \\ &&&\frac{\partial p_i}{\partial z_j}=\frac{0(\Sigma_{j=1}^ne^{z_j})-e^{z_i}e^{z_j}}{(\Sigma_{j=1}^ne^{z_j})^2} =-p_{i}p_{j} \end{aligned} 设 pi=Softmax(zi),则对于 zi 的导数为:∙ 当 i=j 时:∙ 当 i=j 时:∂zi∂pi=(Σj=1nezj)2ezi(Σj=1nezj)−eziezi=pi(1−pi)∂zj∂pi=(Σj=1nezj)20(Σj=1nezj)−eziezj=−pipj
S o f t m a x ( z i ) = e z i − max ( z ) ∑ j = 1 n e z j − max ( z ) \mathrm{Softmax}(z_i)=\frac{e^{z_i-\max(z)}}{\sum_{j=1}^ne^{z_j-\max(z)}} Softmax(zi)=∑j=1nezj−max(z)ezi−max(z)
解释:
z i − max ( z ) z_i-\max(z) zi−max(z)是一个非正数,由于 e z i − m a x ( z ) e^{z_i−max(z)} ezi−max(z) 的形式,当 z i z_i zi 接近 max(z) 时, e z i − m a x ( z ) e^{z_i−max(z)} ezi−max(z) 的值会接近 1,而当 z i z_i zi 远小于 max(z) 时, e z i − m a x ( z ) e^{z_i−max(z)} ezi−max(z) 的值会接近 0。这使得 Softmax 函数的输出中,最大值对应的概率会相对较大,而其他值对应的概率会相对较小,从而提高数值稳定性。
这种调整不会改变 S o f t m a x Softmax Softmax的概率分布结果,因为从数学的角度讲相当于分子、分母都除以了 e max ( z ) e^{\max(z)} emax(z)。
在 PyTorch 中,torch.nn.functional.softmax
函数就自动处理了数值稳定性问题。
代码参考如下:
import torch
import torch.nn as nn
# 表示4分类,每个样本全连接后得到4个得分,下面示例模拟的是两个样本的得分
input_tensor = torch.tensor([[-1.0, 2.0, -3.0, 4.0], [-2, 3, -3, 9]])
softmax = nn.Softmax()
output_tensor = softmax(input_tensor)
# 关闭科学计数法
torch.set_printoptions(sci_mode=False)
print("输入张量:", input_tensor)
print("输出张量:", output_tensor)
输出结果:
输入张量: tensor([[-1., 2., -3., 4.],
[-2., 3., -3., 9.]])
输出张量: tensor([[ 0.0059, 0.1184, 0.0008, 0.8749],
[ 0.0000, 0.0025, 0.0000, 0.9975]])
更多激活函数可以查看官方文档:https://pytorch.org/docs/stable/nn.html#non-linear-activations-weighted-sum-nonlinearity
那这么多激活函数应该如何选择呢?实际没那么纠结