Generalized Hebbian algorithm

The generalized Hebbian algorithm, also known in the literature as Sanger's rule, is a linear feedforward neural network for unsupervised learning with applications primarily in principal components analysis.

First defined in 1989,[1] it is similar to Oja's rule in its formulation and stability, except it can be applied to networks with multiple outputs.

The name originates because of the similarity between the algorithm and a hypothesis made by Donald Hebb[2] about the way in which synaptic strengths in the brain are modified in response to experience, i.e., that changes are proportional to the correlation between the firing of pre- and post-synaptic neurons.

[3] Consider a problem of learning a linear code for some data.

should be the highest principal component vectors.

The generalized Hebbian algorithm is an iterative algorithm to find the highest principal component vectors, in an algorithmic form that resembles unsupervised Hebbian learning in neural networks.

Consider a one-layered neural network with

The linear code vectors are the connection strengths, that is,

is the synaptic weight or connection strength between the

The generalized Hebbian algorithm learning rule is of the form where

is the learning rate parameter.

[4] In matrix form, Oja's rule can be written and the Gram-Schmidt algorithm is where w(t) is any matrix, in this case representing synaptic weights, Q = η x xT is the autocorrelation matrix, simply the outer product of inputs, diag is the function that diagonalizes a matrix, and lower is the function that sets all matrix elements on or above the diagonal equal to 0.

We can combine these equations to get our original rule in matrix form, where the function LT sets all matrix elements above the diagonal equal to 0, and note that our output y(t) = w(t) x(t) is a linear neuron.

[1] [5] Oja's rule is the special case where

[6] One can think of the generalized Hebbian algorithm as iterating Oja's rule.

is learned, and it has the same direction as the largest principal component vector is learned, with length determined by

has converged, the second neuron is effectively running another Oja's rule on the modified input vectors, defined by

, which we know is the input vector with the first principal component removed.

Therefore, the second neuron learns to code for the second principal component.

By induction, this results in finding the top-

principal components for arbitrary

The generalized Hebbian algorithm is used in applications where a self-organizing map is necessary, or where a feature or principal components analysis can be used.

Examples of such cases include artificial intelligence and speech and image processing.

Its importance comes from the fact that learning is a single-layer process—that is, a synaptic weight changes only depending on the response of the inputs and outputs of that layer, thus avoiding the multi-layer dependence associated with the backpropagation algorithm.

It also has a simple and predictable trade-off between learning speed and accuracy of convergence as set by the learning rate parameter η.

[5] As an example, (Olshausen and Field, 1996)[7] performed the generalized Hebbian algorithm on 8-by-8 patches of photos of natural scenes, and found that it results in Fourier-like features.

The features are the same as the principal components found by principal components analysis, as expected, and that, the features are determined by the

variance matrix of the samples of 8-by-8 patches.

In other words, it is determined by the second-order statistics of the pixels in images.

They criticized this as insufficient to capture higher-order statistics which are necessary to explain the Gabor-like features of simple cells in the primary visual cortex.

Features learned by generalized Hebbian algorithm running on 8-by-8 patches of Caltech 101 .
Features found by Principal Component Analysis on the same Caltech 101 dataset.