Hessian matrix

From Wikipedia, the free encyclopedia

In mathematics, the Hessian matrix is the square matrix of second partial derivatives of a scalar-valued function. Given the real-valued function

f(x1,x2,...,xn),

if all second partial derivatives of f exist, then the Hessian matrix of f is the matrix

H(f)ij(x) = DiDjf(x)

where x = (x1,x2,...,xn). That is,

H(f) = \begin{bmatrix} \frac{\partial^2 f}{\partial x_1^2} & \frac{\partial^2 f}{\partial x_1\,\partial x_2} & \cdots & \frac{\partial^2 f}{\partial x_1\,\partial x_n} \\  \\ \frac{\partial^2 f}{\partial x_2\,\partial x_1} & \frac{\partial^2 f}{\partial x_2^2} & \cdots & \frac{\partial^2 f}{\partial x_2\,\partial x_n} \\  \\ \vdots & \vdots & \ddots & \vdots \\  \\ \frac{\partial^2 f}{\partial x_n\,\partial x_1} & \frac{\partial^2 f}{\partial x_n\,\partial x_2} & \cdots & \frac{\partial^2 f}{\partial x_n^2} \end{bmatrix}.

(some workers define the Hessian as the determinant of the above matrix). The term Hessian was coined by James Joseph Sylvester, named for German mathematician Ludwig Otto Hesse, who had used the term functional determinants. Hessian matrices are used in large-scale optimization problems within Newton-type methods; e.g., the BFGS method.

Contents

[edit] Mixed derivatives and symmetry of the Hessian

The mixed derivatives of f are the entries off the main diagonal in the Hessian. Assuming that they are continuous, the order of differentiation does not matter (Schwartz's theorem). For example,

\frac {\partial}{\partial x} \left( \frac { \partial f }{ \partial y} \right) =        \frac {\partial}{\partial y} \left( \frac { \partial f }{ \partial x} \right)

This can also be written (in reverse order) as:

f_{xy} = f_{yx} \,

In a formal statement: if the second derivatives of f are all continuous in a region D, then the Hessian of f is a symmetric matrix throughout D; see symmetry of second derivatives.

[edit] Critical points and discriminant

If the gradient of f (i.e. its derivative in the vector sense) is zero at some point x, then f has a critical point at x. The determinant of the Hessian at x is then called the discriminant. If this determinant is zero then x is called a degenerate critical point of f, this is also called a non-Morse critical point of f. Otherwise it is non-degenerate, this is called a Morse critical point of f.

[edit] Second derivative test

The following test can be applied at a non-degenerate critical point x. If the Hessian is positive definite at x, then f attains a local minimum at x. If the Hessian is negative definite at x, then f attains a local maximum at x. If the Hessian has both positive and negative eigenvalues then x is a saddle point for f (this is true even if x is degenerate). Otherwise the test is inconclusive.

Note that for positive semidefinite and negative semidefinite Hessians the test is inconclusive. However, more can be said from the point of view of Morse theory.

In view of what has just been said, the second derivative test for functions of one and two variables is simple. In one variable, the Hessian contains just one second derivative; if it is positive then x is a local minimum, if it is negative then x is a local maximum; if it is zero then the test is inconclusive. In two variables, the discriminant can be used, because the determinant is the product of the eigenvalues. If it is positive then the eigenvalues are both positive, or both negative. If it is negative then the two eigenvalues have different signs. If it is zero, then the second derivative test is inconclusive.

[edit] Bordered Hessian

A bordered Hessian is utilized for the second derivate test in certain constrained optimization problems. Given the function as before:

f(x1,x2,...,xn),

but adding a constraint function such that:

g(x1,x2,...,xn) = c,

the bordered Hessian appears as

H(f,g) = \begin{bmatrix} 0 & \frac{\partial g}{\partial x_1} & \frac{\partial g}{\partial x_2} & \cdots & \frac{\partial g}{\partial x_n} \\  \\ \frac{\partial g}{\partial x_1} & \frac{\partial^2 f}{\partial x_1^2} & \frac{\partial^2 f}{\partial x_1\,\partial x_2} & \cdots & \frac{\partial^2 f}{\partial x_1\,\partial x_n} \\  \\ \frac{\partial g}{\partial x_2} & \frac{\partial^2 f}{\partial x_2\,\partial x_1} & \frac{\partial^2 f}{\partial^2 x_2} & \cdots & \frac{\partial^2 f}{\partial x_2\,\partial x_n} \\  \\ \vdots & \vdots & \vdots & \ddots & \vdots \\  \\ \frac{\partial g}{\partial x_n} & \frac{\partial^2 f}{\partial x_n\,\partial x_1} & \frac{\partial^2 f}{\partial x_n\,\partial x_2} & \cdots & \frac{\partial^2 f}{\partial x_n^2} \end{bmatrix}.

Then the same rules of positive definite and negative definite apply as before.

[edit] Vector-valued functions

If f is instead vector-valued, i.e.

f = (f1,f2,...,fn),

then the array of second partial derivatives is not a matrix, but a tensor of rank 3.

[edit] See also