Activation function

From Infogalactic: the planetary knowledge core
Jump to: navigation, search

<templatestyles src="Module:Hatnote/styles.css"></templatestyles>

In computational networks, the activation function of a node defines the output of that node given an input or set of inputs. A standard computer chip circuit can be seen as a digital network of activation functions that can be "ON" (1) or "OFF" (0), depending on input. This is similar to the behavior of the linear perceptron in neural networks. However, it is the nonlinear activation function that allows such networks to compute nontrivial problems using only a small number of nodes. In artificial neural networks this function is also called transfer function (not to be confused with a linear system’s transfer function).

Functions

In biologically inspired neural networks, the activation function is usually an abstraction representing the rate of action potential firing in the cell. In its simplest form, this function is binary—that is, either the neuron is firing or not. The function looks like \phi(v_i)=U(v_i), where U is the Heaviside step function. In this case a large number of neurons must be used in computation beyond linear separation of categories.

A line of positive slope may also be used to reflect the increase in firing rate that occurs as input current increases. The function would then be of the form \phi(v_i)=\mu v_i, where \mu is the slope. This activation function is linear, and therefore has the same problems as the binary function. In addition, networks constructed using this model have unstable convergence because neuron inputs along favored paths tend to increase without bound, as this function is not normalizable.

All problems mentioned above can be handled by using a normalizable sigmoid activation function. One realistic model stays at zero until input current is received, at which point the firing frequency increases quickly at first, but gradually approaches an asymptote at 100% firing rate. Mathematically, this looks like \phi(v_i)=U(v_i)\tanh(v_i), where the hyperbolic tangent function can also be replaced by any sigmoid function. This behavior is realistically reflected in the neuron, as neurons cannot physically fire faster than a certain rate. This model runs into problems, however, in computational networks as it is not differentiable, a requirement in order to calculate backpropagation.

The final model, then, that is used in multilayer perceptrons is a sigmoidal activation function in the form of a hyperbolic tangent. Two forms of this function are commonly used: \phi(v_i)=\tanh(v_i) whose range is normalized from -1 to 1, and \phi(v_i) = (1+\exp(-v_i))^{-1} is vertically translated to normalize from 0 to 1. The latter model is often considered more biologically realistic, but it runs into theoretical and experimental difficulties with certain types of computational problems.

Alternative structures

A special class of activation functions known as radial basis functions (RBFs) are used in RBF networks, which are extremely efficient as universal function approximators. These activation functions can take many forms, but they are usually found as one of three functions:

  • Gaussian: \,\phi(v_i)=\exp\left(-\frac{\|v_i-c_i\|^2}{2\sigma^2}\right)
  • Multiquadratics: \,\phi(v_i)=\sqrt{\|v_i-c_i\|^2 + a^2}
  • Inverse multiquadratics: \,\phi(v_i)=(\|v_i-c_i\|^2 + a^2)^{-1/2}

where c_i is the vector representing the function center and a and \sigma are parameters affecting the spread of the radius.

Support vector machines (SVMs) can effectively utilize a class of activation functions that includes both sigmoids and RBFs. In this case, the input is transformed to reflect a decision boundary hyperplane based on a few training inputs called support vectors x. The activation function for the hidden layer of these machines is referred to as the inner product kernel, K(v_i,x)=\phi(v_i). The support vectors are represented as the centers in RBFs with the kernel equal to the activation function, but they take a unique form in the perceptron as

\,\phi(v_i)=\tanh\left(\beta_1 + \beta_0 \sum_j v_{i,j} x_j \right),

where \beta_0 and \beta_1 must satisfy certain conditions for convergence. These machines can also accept arbitrary-order polynomial activation functions where

\,\phi(v_i)=\left(1+\sum_j v_{i,j} x_j \right)^p.[1]

Activation function having types:

  1. Identity function.
  2. Binary step function.
  3. Bipolar step function.
  4. Sigmoidal function.
    1. Binary sigmoidal function.
    2. Bipolar sigmoidal function.
  5. Ramp function.

Comparison of activation functions

Some desirable properties in an activation function include:

  • Nonlinear: When the activation function is non-linear, then a two-layer neural network can be proven to be a universal function approximator.[2] The identity activation function does not satisfy this property. When multiple layers use the identity activation function, the entire network is equivalent to a single-layer model.
  • Continuously differentiable: This property is necessary for enabling gradient-based optimization methods. The binary step activation function is not differentiable at 0, and it differentiates to 0 for all other values, so gradient-based methods can make no progress with it.[3]
  • Range: When the range of the activation function is finite, gradient-based training methods tend to be more stable, because pattern presentations significantly affect only limited weights. When the range is infinite, training is generally more efficient because pattern presentations significantly affect most of the weights. In the latter case, smaller learning rates are typically necessary.[citation needed]
  • Monotonic: When the activation function is monotonic, the error surface associated with a single-layer model is guaranteed to be convex.[4]
  • Smooth Functions with a Monotonic derivative have been shown to generalize better in some cases. The argument for these properties suggests that such activation functions are more consistent with Occam's razor.[5]
  • f(x)\approx x when x \approx 0:[clarification needed] This property enables the neural network to train efficiently when its weights are initialized with small random values. When the activation function does not satisfy this property, special care must be used when initializing the weights [6]

The following table compares the properties of several activation functions:

Name Plot Equation Derivative Range Order of continuity Monotonic Derivative Monotonic f(x)\approx x when x \approx 0
Identity File:Activation identity.svg f(x)=x f'(x)=1 (-\infty,\infty) C^\infty Yes Yes Yes
Binary step File:Activation binary step.svg f(x) = \left \{	\begin{array}{rcl}
	0 & \mbox{for} & x < 0\\
	1 & \mbox{for} & x \ge 0\end{array} \right. f'(x) = \left \{	\begin{array}{rcl}
	0 & \mbox{for} & x \ne 0\\
	? & \mbox{for} & x = 0\end{array} \right. \{0,1\} C^{-1} Yes No No
Logistic (a.k.a Soft step) File:Activation logistic.svg f(x)=\frac{1}{1+e^{-x}} f'(x)=f(x)(1-f(x)) (0,1) C^\infty Yes No No
TanH File:Activation tanh.svg f(x)=\tanh(x)=\frac{2}{1+e^{-2x}}-1 f'(x)=1-f(x)^2 (-1,1) C^\infty Yes No Yes
ArcTan File:Activation arctan.svg f(x)=\tan^{-1}(x) f'(x)=\frac{1}{x^2+1} (-\frac{\pi}{2},\frac{\pi}{2}) C^\infty Yes No Yes
Softsign [7] File:Activation softsign.png f(x)=\frac{x}{1+|x|} f'(x)=\frac{1}{(1+|x|)^2} (-1,1) C^1 Yes No Yes
Rectified Linear Unit (ReLU)[8] File:Activation rectified linear.svg f(x) = \left \{	\begin{array}{rcl}
	0 & \mbox{for} & x < 0\\
	x & \mbox{for} & x \ge 0\end{array} \right. f'(x) = \left \{	\begin{array}{rcl}
	0 & \mbox{for} & x < 0\\
	1 & \mbox{for} & x \ge 0\end{array} \right. [0,\infty) C^0 Yes Yes Yes
Parameteric Rectified Linear Unit (PReLU)[9] File:Activation prelu.svg f(x) = \left \{	\begin{array}{rcl}
	\alpha x & \mbox{for} & x < 0\\
	x & \mbox{for} & x \ge 0\end{array} \right. f'(x) = \left \{	\begin{array}{rcl}
	\alpha & \mbox{for} & x < 0\\
	1 & \mbox{for} & x \ge 0\end{array} \right. (-\infty,\infty) C^0 Yes Yes Yes
Exponential Linear Unit (ELU)[10] File:Activation elu.svg f(x) = \left \{	\begin{array}{rcl}
	\alpha(e^x - 1) & \mbox{for} & x < 0\\
	x & \mbox{for} & x \ge 0\end{array} \right. f'(x) = \left \{	\begin{array}{rcl}
	f(x) + \alpha & \mbox{for} & x < 0\\
	1 & \mbox{for} & x \ge 0\end{array} \right. (-\alpha,\infty) C^1 when
\alpha=1
Yes Yes Yes iff \alpha\approx 1
SoftPlus[11] File:Activation softplus.svg f(x)=\log_e(1+e^x) f'(x)=\frac{1}{1+e^{-x}} (0,\infty) C^\infty Yes Yes No
Bent identity File:Activation bent identity.svg f(x)=\frac{\sqrt{x^2 + 1} - 1}{2} + x f'(x)=\frac{x}{2\sqrt{x^2 + 1}} + 1 (-\infty,\infty) C^\infty Yes Yes Yes
SoftExponential [12] File:Activation soft exponential.svg f(\alpha,x) = \left \{	\begin{array}{rcl}
	-\frac{\log_e(1-\alpha (x + \alpha))}{\alpha} & \mbox{for} & \alpha < 0\\
	x & \mbox{for} & \alpha = 0\\
        \frac{e^{\alpha x} - 1}{\alpha} + \alpha & \mbox{for} & \alpha > 0\end{array}\right. f'(\alpha,x) = \left \{	\begin{array}{rcl}
	\frac{1}{1-\alpha (\alpha + x)} & \mbox{for} & \alpha < 0\\
	e^{\alpha x} & \mbox{for} & \alpha \ge 0\end{array} \right. (-\infty,\infty) C^\infty Yes Yes Yes iff \alpha\approx 0
Sinusoid File:Activation sinusoid.svg f(x)=\sin(x) f'(x)=\cos(x) [-1,1] C^\infty No No Yes
Sinc File:Activation sinc.svg f(x)=\left \{	\begin{array}{rcl}
        1 & \mbox{for} & x = 0\\
	\frac{\sin(x)}{x} & \mbox{for} & x \ne 0\end{array} \right. f'(x)=\left \{	\begin{array}{rcl}
	0 & \mbox{for} & x = 0\\
	\frac{\cos(x)}{x} - \frac{\sin(x)}{x^2} & \mbox{for} & x \ne 0\end{array} \right. [\approx-.217234,1] C^\infty No No No
Gaussian File:Activation gaussian.svg f(x)=e^{-x^2} f'(x)=-2xe^{-x^2} (0,1] C^\infty No No No

See also

References

  1. Lua error in package.lua at line 80: module 'strict' not found.
  2. Cybenko, George. "Approximation by superpositions of a sigmoidal function." Mathematics of control, signals and systems 2.4 (1989): 303-314.
  3. Snyman, Jan. Practical mathematical optimization: an introduction to basic optimization theory and classical and new gradient-based algorithms. Vol. 97. Springer Science & Business Media, 2005.
  4. Wu, Huaiqin. "Global stability analysis of a general class of discontinuous neural networks with linear growth activation functions." Information Sciences 179.19 (2009): 3432-3441.
  5. Gashler, Michael S., and Stephen C. Ashmore. "Training Deep Fourier Neural Networks to Fit Time-Series Data." Intelligent Computing in Bioinformatics. Springer International Publishing, 2014. 48-55, arXiv:1405.2262
  6. Sussillo, David, and L. F. Abbott. "Random walks: Training very deep nonlinear feed-forward networks with smart initialization." CoRR, arXiv:1412.6558 (2014): 286.
  7. Glorot, Xavier and Bengio, Yoshua."Understanding the difficulty of training deep feedforward neural networks". In Proceedings of the International Conference on Artificial Intelligence and Statistics (AISTATS’10). Society for Artificial Intelligence and Statistics. 2010. http://jmlr.org/proceedings/papers/v9/glorot10a/glorot10a.pdf
  8. Nair, Vinod, and Geoffrey E. Hinton. "Rectified linear units improve restricted boltzmann machines." Proceedings of the 27th International Conference on Machine Learning (ICML-10). 2010.
  9. arXiv:1502.01852
  10. arXiv:1511.07289v3
  11. Glorot, Xavier, Antoine Bordes, and Yoshua Bengio. "Deep sparse rectifier neural networks." International Conference on Artificial Intelligence and Statistics. 2011.
  12. Luke B. Godfrey and Gashler, Michael S. A Continuum among Logarithmic, Linear, and Exponential Functions, and Its Potential to Improve Generalization in Neural Networks. In Proceedings of the 7th International Joint Conference on Knowledge Discovery, Knowledge Engineering and Knowledge Management: KDIR, pages 481-486. Lisbon, Portugal, November, 2015, arXiv:1602.01321.