社区首页 >专栏 >Softmax和Cross-entropy是什么关系？

Softmax和Cross-entropy是什么关系？

公众号机器学习与AI生成创作

发布于 2020-09-14 09:27:02

1.5K00

代码可运行

文章被收录于专栏：机器学习与生成对抗网络机器学习与生成对抗网络

运行总次数：0

代码可运行

来自 | 知乎作者 | 董鑫

https://www.zhihu.com/question/294679135/answer/885285177

文仅分享，著作权归作者，侵删

softmax 虽然简单，但是其实这里面有非常的多细节值得一说。

我们挨个捋一捋。

1. 什么是 Softmax？

首先，softmax 的作用是把一个序列，变成概率。

他能够保证：

所有的值都是 [0, 1] 之间的（因为概率必须是 [0, 1]）
所有的值加起来等于 1

从概率的角度解释 softmax 的话，就是

2. 文档里面跟 Softmax 有关的坑

这里穿插一个“小坑”，很多deep learning frameworks的文档里面（PyTorch，TensorFlow）是这样描述 softmax 的，

take logits and produce probabilities

很明显，这里面的 logits 就是全连接层（经过或者不经过 activation都可以）的输出， probability 就是 softmax 的输出结果。这里 logits 有些地方还称之为 unscaled log probabilities。这个就很意思了，unscaled probability可以理解，那又为什么全连接层直接出来结果会和 log 有关系呢？

原因有两个：

因为全连接层出来的结果，其实是无界的（有正有负），这个跟概率的定义不一致，但是你如果他看成概率的 log，就可以理解了。
softmax 的作用，我们都知道是 normalize probability。在 softmax 里面，输入

都是在指数上的

，所有把

想成 log of probability 也就顺理成章了。

3. Softmax 就是 Soft 版本的 ArgMax

好的，我们把话题拉回到 softmax。

softmax，顾名思义就是 soft 版本的 argmax。我们来看一下为什么？

举个栗子，假如 softmax 的输入是：

softmax 的结果是：

我们稍微改变一下输入，把 3 改大一点，变成 5，输入是

softmax 的结果是：

可见 softmax 是一种非常明显的 “马太效应”：强（大）的更强（大），弱（小）的更弱（小）。假如你要选一个最大的数出来，这个其实就是叫 hardmax。那么 softmax 呢，其实真的就是 soft 版本的 max，以一定的概率选一个最大值出来。在hardmax中，真正最大的那个数，一定是以1(100%) 的概率被选出来，其他的值根本一点机会没有。但是在 softmax 中，所有的值都有机会被作为最大值选出来。只不过，由于 softmax 的 “马太效应”，次大的数，即使跟真正最大的那个数差别非常少，在概率上跟真正最大的数相比也小了很多。

所以，前面说，“softmax 的作用是把一个序列，变成概率。” 这个概率不是别的，而是被选为 max 的概率。

这种 soft 版本的 max 在很多地方有用的上。因为 hard 版本的 max 好是好，但是有很严重的梯度问题，求最大值这个函数本身的梯度是非常非常稀疏的（比如神经网络中的 max pooling），经过hardmax之后，只有被选中的那个变量上面才有梯度，其他都是没有梯度。这对于一些任务（比如文本生成等）来说几乎是不可接受的。所以要么用 hard max 的变种，比如Gumbel，

Categorical Reparameterization with Gumbel-Softmax

链接：https://arxiv.org/abs/1611.01144

亦或是 ARSM

ARSM: Augment-REINFORCE-Swap-Merge Estimator for Gradient Backpropagation Through Categorical Variable

链接：http://proceedings.mlr.press/v97/yin19c.html

，要么就直接 softmax。

4. Softmax 的实现以及数值稳定性

softmax 的代码实现看似是比较简单的，直接套上面的公式就好

def softmax(x):
    """Compute the softmax of vector x."""
    exps = np.exp(x)
    return exps / np.sum(exps)

但是这种方法非常的不稳定。因为这种方法要算指数，只要你的输入稍微大一点，比如：

分母上就是

很明显，在计算上一定会溢出。解决方法也比较简单，就是我们在分子分母上都乘上一个系数，减小数值大小，同时保证整体还是对的

把常数 C 吸收进指数里面

这里的D是可以随便选的，一般可以选成

具体实现可以写成这样


def stablesoftmax(x):
    """Compute the softmax of vector x in a numerically stable way."""
    shiftx = x - np.max(x)
    exps = np.exp(shiftx)
    return exps / np.sum(exps)