如果不用激励函数(其实相当于激励函数是f(x) = x),在这种情况下你每一层输出实际上都是上层输入的线性函数。
这样就使得无论神经网络有多少层,输出都是输入的线性组合,与没有隐藏层效果相当,模型的表达力仍然不够。
我们决定引入非线性函数作为激励函数,这样深层神经网络才有意义(不再是输入的线性组合)。
本文将介绍深度学习中的4个常见的激活函数,从原函数公式、导数函数及二者的可视化来进行对比:
表示形式为tf.nn.sigmoid(x)
In [1]:
# import tensorflow as
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
%matplotlib inline
def sigmoid(x):
"""
返回sigmoid函数
"""
return 1 / (1 + np.exp(-x))
def plot_sigmoid():
# param:起点,终点,间距
x = np.arange(-10, 10, 0.2)
y = sigmoid(x)
plt.plot(x, y)
plt.grid()
plt.show()
if __name__ == '__main__':
plot_sigmoid()
该函数的导数为:
$$\begin{aligned} f^{\prime}(z) &=\left(\frac{1}{1+e{-z}}\right){\prime} \ &=\frac{e{-z}}{\left(1+e{-z}\right)^2} \ &=\frac{1+e{-z}-1}{\left(1+e{-z}\right)^2} \ &=\frac{1}{\left(1+e{-z}\right)}\left(1-\frac{1}{\left(1+e{-z}\right)}\right) \ &=f(z)(1-f(z)) \end{aligned} $$
另一种求解方法:
步骤1:
$$\begin{aligned} \frac{\mathrm{d} y}{\mathrm{~d} x} &=-\left(1+e{-x}\right){-2} \cdot\left(1+e{-x}\right){\prime} \ &=-\left(1+e{-x}\right){-2} \cdot 1 \cdot\left(e{-x}\right){\prime} \ &=-\left(1+e{-x}\right){-2} \cdot 1 \cdot\left(e^{-x}\right) \cdot(-x)^{\prime} \ &=-\left(1+e{-x}\right){-2} \cdot 1 \cdot\left(e^{-x}\right) \cdot(-1) \ &=\left(1+e{-x}\right){-2} \cdot\left(e^{-x}\right) \ &=\frac{e{-x}}{\left(1+e{-x}\right)^2} \end{aligned}$$
步骤2:
步骤3:
# import tensorflow as
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
%matplotlib inline
def der_sigmoid(x):
"""
返回sigmoid函数
"""
sig = 1 / (1 + np.exp(-x)) # sigmoid函数
return sig * (1 - sig) # 输出为导数函数
def plot_der_sigmoid():
# param:起点,终点,间距
x = np.arange(-10, 10, 0.2)
y = der_sigmoid(x) # 导数函数
plt.plot(x, y)
plt.grid()
plt.show()
if __name__ == '__main__':
plot_der_sigmoid()
Sigmoid函数是二分类算法,尤其是逻辑回归中算法中的常用激活函数;也可以作为较少层数的神经网络的激活函数。
它主要是有下面几个特点:
下面解释下Sigmoid最为致命的缺点(参考一位网友的回答,解释很全面)
对于Sigmoid激活函数最致命的缺点就是容易发生梯度弥散(Gradient vanishing)现象(当然也可能会发生梯度爆炸Exploding gradient,前面层的梯度通过模型训练变的很大,由于反向传播中链式法则的原因,导致后面层的梯度值会以指数级增大。 但是在Sigmoid激活函数中梯度保障发生的概率非常小),所谓梯度弥散故名思议就是梯度值越来越小。 在深度学习中,梯度更新是从后向前更新的,这也就是所谓的反向传播(Backpropagation algorithm),而反向传播的核心是链式法则。如果使用Sigmoid激活函数,训练的网络比较浅还比较好,但是一旦训练较深的神经网络,会导致每次传过来的梯度都会乘上小于1的值,多经过几层之后,梯度就会变得非常非常小(逐渐接近于0),因此梯度消失了,对应的参数得不到更新。 参考:https://zhuanlan.zhihu.com/p/104318223
现在神经网络中很少使用Sigmoid作为激活函数。
tanh 是一个双曲正切函数。tanh 函数和 sigmoid 函数的曲线相似。
该函数是一个奇函数,经过原点且严格单调递增。函数取值分布在(-1,1)之间。
下面是原函数的具体表达形式,我们还发现tanh函数和sigmoid函数的关系:
import math
import matplotlib.pyplot as plt
import numpy as np
import matplotlib as mpl
def tanh(x):
return (np.exp(x) - np.exp(-x)) / (np.exp(x) + np.exp(-x))
fig = plt.figure(figsize=(12, 8))
ax = fig.add_subplot(111)
x = np.linspace(-10, 10)
y = tanh(x)
ax.spines['top'].set_color('none')
ax.spines['right'].set_color('none')
ax.xaxis.set_ticks_position('bottom')
ax.spines['bottom'].set_position(('data', 0))
ax.set_xticks([-10, -5, 0, 5, 10])
ax.yaxis.set_ticks_position('left')
ax.spines['left'].set_position(('data', 0))
ax.set_yticks([-1, -0.5, 0.5, 1])
plt.plot(x, y, label="Tanh", color="red")
plt.legend()
plt.grid()
plt.show()
下面是具体的求导过程:
首先公布结果,导数函数为:
下面是具体过程:
步骤1、在除法的求导公式中有:
步骤2、定义两个中间变量:
步骤3、第一次求导
步骤4、其中有
步骤5、将a+b和a-b进行整理:
步骤6、化简上面的结果
$$\begin{array}{l} \frac{(a-b)^{\prime} \times(a+b)-(a-b) \times(a+b){\prime}}{(a+b)2} \ =\frac{(a+b)2-(a-b)2}{(a+b)^2} \ =1-\frac{(a-b)2}{(a+b)2} \ =1-\left(\frac{ex-e{-x}}{ex+e{-x}}\right)^2 \ =1-\operatorname{tanh}(x)^2 \end{array}$$
# 绘制图
import math
import matplotlib.pyplot as plt
import numpy as np
import matplotlib as mpl
def der_tanh(x):
"""
tanh(x)函数的导数求解
"""
tanh = (np.exp(x) - np.exp(-x)) / (np.exp(x) + np.exp(-x))
return 1-tanh ** 2
fig = plt.figure(figsize=(12, 8))
ax = fig.add_subplot(111)
x = np.linspace(-10, 10)
y = der_tanh(x)
ax.spines['top'].set_color('none')
ax.spines['right'].set_color('none')
ax.xaxis.set_ticks_position('bottom')
ax.spines['bottom'].set_position(('data', 0))
ax.set_xticks([-10, -5, 0, 5, 10])
ax.yaxis.set_ticks_position('left')
ax.spines['left'].set_position(('data', 0))
ax.set_yticks([-1, -0.5, 0.5, 1])
plt.plot(x, y, label="der_tanh", color="red")
plt.legend()
plt.grid()
plt.show()
# encoding:utf-8
import math
import matplotlib.pyplot as plt
import numpy as np
import matplotlib as mpl
mpl.rcParams['axes.unicode_minus']=False
def sigmoid(x):
return 1.0 / (1.0 + np.exp(-x))
fig = plt.figure(figsize=(10,6))
ax = fig.add_subplot(111)
x = np.linspace(-10, 10)
# sigmoid函数
y = sigmoid(x)
# tanh函数
tanh = 2*sigmoid(2*x) - 1
plt.xlim(-11,11)
plt.ylim(-1.1,1.1)
ax.spines['top'].set_color('none')
ax.spines['right'].set_color('none')
ax.xaxis.set_ticks_position('bottom')
ax.spines['bottom'].set_position(('data',0))
ax.set_xticks([-10,-5,0,5,10])
ax.yaxis.set_ticks_position('left')
ax.spines['left'].set_position(('data',0))
ax.set_yticks([-1,-0.5,0.5,1])
plt.plot(x,y,label="Sigmoid",color = "blue")
# tanh函数这里是2*x
plt.plot(2*x,tanh,label="Tanh", color = "red")
plt.legend()
plt.grid()
plt.show()
在一般的二分类问题中,Sigmoid函数一般用输出层,而tanh函数一般用于隐藏层(非固定,经验之谈)。
ReLu函数是目前深度学习中比较流行的一种激活函数。
ReLU函数, 也称之为线性整流函数(Rectified Linear Unit), 是神经网络结构中常用的非线性激活函数。下面是简写形式:
下面是完整的表达式:
$$\operatorname{ReLU}(x)=\left{\begin{array}{ll} 0, & x \leqslant 0 \ x, & x>0 \end{array}\right.$$
# import tensorflow as tf
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
%matplotlib inline
def relu(x):
"""
返回Relu函数
"""
# 不能直接使用max函数,报错
return np.maximum(0, x)
def plot_relu():
# param:起点,终点,间距
plt.figure(figsize=(12,6))
x = np.arange(-10, 10, 0.2)
y = relu(x) # 导数函数
plt.plot(x, y)
plt.grid()
plt.show()
if __name__ == '__main__':
plot_relu()
$$y{\prime}=f{\prime}(x)=\left{\begin{array}{ll} 1, & x>0 \ 0, & x \leq 0 \end{array}\right.$$
# import tensorflow as tf
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
%matplotlib inline
def der_relu(x):
"""
返回Relu导数函数
"""
# 重点:通过numpy的where函数进行判断
return np.where(x>0,1,0)
def plot_der_relu():
# param:起点,终点,间距
plt.figure(figsize=(12,6))
x = np.arange(-10, 10, 0.2)
y = der_relu(x) # 导数函数
plt.plot(x, y)
plt.grid()
plt.show()
if __name__ == '__main__':
plot_der_relu()
Leaky ReLu 函数是为了解决Relu函数负区间的取值为0而产生的。
在小于0的部分引入了一个斜率,使得小于0的部分取值不再全部是0(通常 a 的值为 0.01 左右)
$$y=f(x)=\left{\begin{array}{ll} x, & x>0 \ a \cdot x, & x \leq 0 \end{array}\right.$$
# import tensorflow as tf
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
%matplotlib inline
def leaky_relu(x):
"""
返回Leaky Relu函数图像
"""
a = 0.01 # 斜率a的取值
return np.maximum(a * x, x) # 通过np.maximum函数实现
def plot_Leaky_relu():
# param:起点,终点,间距
plt.figure(figsize=(12,6))
x = np.arange(-10, 10, 0.2)
y = leaky_relu(x) # 导数函数
plt.plot(x, y, color="red")
plt.grid()
plt.show()
if __name__ == '__main__':
plot_Leaky_relu()
$$y=f
(x)=\left{\begin{array}{ll}
1, & x>0 \
a, & x \leq 0
\end{array}\right.$$
# import tensorflow as tf
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
%matplotlib inline
def der_leaky_relu(x):
"""
返回Leaky Relu函数图像
"""
# 斜率a的取值
a = 0.1
# 通过where判断来取值
return np.where(x>0,1,a)
def plot_Leaky_relu():
# param:起点,终点,间距
plt.figure(figsize=(12,6))
x = np.arange(-10, 10, 0.2)
y = der_leaky_relu(x) # 导数函数
plt.plot(x, y, color = "red")
plt.grid()
plt.show()
if __name__ == '__main__':
plot_Leaky_relu()