In mathematics, the Hessian matrix or Hessian is a square matrix of second-order partial derivatives of a scalar-valued function, or scalar field. It describes the local curvature of a function of many variables. The Hessian matrix was developed in the 19th century by the German mathematician Ludwig Otto Hesse and later named after him. Hesse originally used the term "functional determinants".
Specifically, suppose f : ℝn → ℝ is a function taking as input a vector x ∈ ℝn and outputting a scalar f(x) ∈ ℝ; if all second partial derivatives of f exist and are continuous over the domain of the function, then the Hessian matrix H of f is a square n×n matrix, usually defined and arranged as follows:
The Hessian matrix can be considered related to the Jacobian matrix by H(f)(x) = J(∇f)(x).
where J is the Jacobian matrix. The full Hessian matrix can be difficult to compute in practice; in such situations, quasi-Newton algorithms have been developed that use approximations to the Hessian. One of the most popular quasi-Newton algorithms is BFGS.
Mixed derivatives and symmetry of the Hessian
|This section does not cite any references or sources. (April 2013)|
If the gradient (the vector of the partial derivatives) of a function f is zero at some point x, then f has a critical point (or stationary 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, or a non-Morse critical point of f. Otherwise it is non-degenerate, and called a Morse critical point of f.
Second derivative test
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. Otherwise the test is inconclusive. This implies that, at a local minimum (resp. a local maximum), the Hessian is positive-semi-definite (resp. negative semi-definite).
Note that for positive semidefinite and negative semidefinite Hessians the test is inconclusive (yet a conclusion can be made that f is locally convex or concave respectively). However, more can be said from the point of view of Morse theory.
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, and if it is negative then x is a local maximum; if it is zero then the test is inconclusive. In two variables, the determinant 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.
More generally, the second-order conditions that are sufficient for a local minimum or maximum can be expressed in terms of the sequence of principal (upper-leftmost) minors (determinants of sub-matrices) of the Hessian; these conditions are a special case of those given in the next section for bordered Hessians for constrained optimization—the case in which the number of constraints is zero.
A bordered Hessian is used for the second-derivative test in certain constrained optimization problems. Given the function f considered previously, but adding a constraint function g such that g(x) = c, the bordered Hessian appears as
If there are, say, m constraints then the zero in the north-west corner is an m × m block of zeroes, and there are m border rows at the top and m border columns at the left.
The above rules stating that extrema are characterized (among critical points with a nonsingular Hessian) by a positive-definite or negative-definite Hessian cannot apply here since a bordered Hessian can neither be negative-definite nor positive-definite, as zT H z = 0 if z is any vector whose sole non-zero entry is its first.
The second derivative test consists here of sign restrictions of the determinants of a certain set of n - m submatrices of the bordered Hessian. Intuitively, one can think of the m constraints as reducing the problem to one with n - m free variables. (For example, the maximization of f(x1,x2,x3) subject to the constraint x1+x2+x3 = 1 can be reduced to the maximization of f(x1,x2,1-x1-x2) without constraint.)
Specifically, sign conditions are imposed on the sequence of principal minors (determinants of upper-left-justified sub-matrices) of the bordered Hessian, the smallest minor consisting of the truncated first 2m+1 rows and columns, the next consisting of the truncated first 2m+2 rows and columns, and so on, with the last being the entire bordered Hessian. There are thus n–m minors to consider. A sufficient condition for a local maximum is that these minors alternate in sign with the smallest one having the sign of (–1)m+1. A sufficient condition for a local minimum is that all of these minors have the sign of (–1)m (in the unconstrained case of m=0 these conditions coincide with the conditions for the unbordered Hessian to be negative definite or positive definite respectively).
|This section does not cite any references or sources. (April 2013)|
If f is instead a vector field f : ℝn → ℝm, i.e.
then the collection of second partial derivatives is not a n×n matrix, but rather a third order tensor. This can be thought of as an array of m Hessian matrices, one for each component of f:
This tensor promptly degenerates to the usual Hessian matrix when m = 1.
Generalizations to Riemannian manifolds
- by ,
where we have taken advantage of the first covariant derivative of a function being the same as its ordinary derivative. Choosing local coordinates we obtain the local expression for the Hessian as
where are the Christoffel symbols of the connection. Other equivalent forms for the Hessian are given by
- and .
- The determinant of the Hessian matrix is a covariant; see Invariant of a binary form
- Polarization identity, useful for rapid calculations involving Hessians.
- Jacobian matrix
- Hessian equations
- The Hessian matrix is commonly used for expressing image processing operators in image processing and computer vision (see the Laplacian of Gaussian (LoG) blob detector, the determinant of Hessian (DoH) blob detector and scale space).
- Binmore, Ken; Davies, Joan (2007). Calculus Concepts and Methods. Cambridge University Press. p. 190. ISBN 9780521775410. OCLC 717598615.
- Nocedal, Jorge; Wright, Stephen (2000). Numerical Optimization. Springer Verlag. ISBN 978-0387987934.
- Neudecker, Heinz; Magnus, Jan R. (1988). Matrix Differential Calculus with Applications in Statistics and Econometrics. New York: John Wiley & Sons. p. 136. ISBN 0-471-91516-5.
- Chiang, Alpha C. (1984). Fundamental Methods of Mathematical Economics (Third ed.). McGraw-Hill. p. 386.