Hessian matrix

This article describes an analogue for functions of multiple variables of the following term/fact/notion for functions of one variable: second derivative

Definition

Definition in terms of Jacobian matrix and gradient vector

Suppose $f$ is a real-valued function of $n$ variables $x_{1},x_{2},\dots ,x_{n}$. The 'Hessian matrix of $f$ is a $n\times n$-matrix-valued function with domain a subset of the domain of $f$, defined as follows: the Hessian matrix at any point in the domain is the Jacobian matrix of the gradient vector of $f$ at the point. In point-free notation, we denote by $H(f)$ the Hessian matrix function, and we define it as:

$H(f)=J(\nabla f)$

Interpretation as second derivative

The Hessian matrix function is the correct notion of second derivative for a real-valued function of $n$ variables. Here's why:

• The correct notion of first derivative for a scalar-valued function of multiple variables is the gradient vector, so the correct notion of first derivative for $f$ is $\nabla f$.
• The gradient vector $\nabla f$ is itself a vector-valued function with $n$-dimensional inputs and $n$-dimensional outputs. The correct notion of derivative for that is the Jacobian matrix, with $n$-dimensional inputs and outputs valued in $n\times n$-matrices.

Thus, the Hessian matrix is the correct notion of second derivative.

Definition in terms of second-order partial derivatives

For further information, refer: Relation between Hessian matrix and second-order partial derivatives

Wherever the Hessian matrix for a function exists, its entries can be described as second-order partial derivatives of the function. Explicitly, for a function $f$ is a real-valued function of $n$ variables $x_{1},x_{2},\dots ,x_{n}$, the Hessian matrix $H(f)$ is a $n\times n$-matrix-valued function whose $(ij)^{th}$ entry is the second-order partial derivative $\partial ^{2}f/(\partial x_{j}\partial x_{i})$, which is the same as $f_{x_{i}x_{j}}$. Note that the diagonal entries give second-order pure partial derivatives whereas the off-diagonal entries give second-order mixed partial derivatives.

Some people choose to define the Hessian matrix as the matrix whose entries are the second-order partial derivatives as indicated here. However, that is not quite the correct definition of Hessian matrix because it is possible for all the second-order partial derivatives to exist but for the function to not be twice differentiable at the point. The main disadvantage of defining the Hessian matrix in the more expansive sense (i.e., in terms of second-order partial derivatives) is that all the important results about the Hessian matrix crucially rely on the function being twice differentiable, so we don't actually gain anything by using the more expansive definition.

Computationally useful definition at a point

For a function of two variables at a point

Suppose $f$ is a real-valued function of two variables $x,y$ and $(x_{0},y_{0})$ is a point in the domain of $f$ at which $f$ is twice differentiable. In particular, this means that all the four second-order partial derivatives exist at $(x_{0},y_{0})$, i.e., the two pure second-order partials $f_{xx}(x_{0},y_{0}),f_{yy}(x_{0},y_{0})$ exist, and so do the two second-order mixed partial derivatives $f_{xy}(x_{0},y_{0})$ and $f_{yx}(x_{0},y_{0})$. Then, the Hessian matrix of $f$ at $(x_{0},y_{0})$, denoted $H(f)(x_{0},y_{0})$, can be expressed explicitly as a $2\times 2$ matrix of real numbers defined as follows:

${\begin{pmatrix}f_{xx}(x_{0},y_{0})&f_{xy}(x_{0},y_{0})\\f_{yx}(x_{0},y_{0})&f_{yy}(x_{0},y_{0})\\\end{pmatrix}}$

{{#widget:YouTube|id=47WX0VfWS8k}}

For a function of multiple variables at a point

Suppose $f$ is a real-valued function of multiple variables $(x_{1},x_{2},\dots ,x_{n})$. Suppose $(a_{1},a_{2},\dots ,a_{n})$ is a point in the domain of $f$ at which $f$ is twice differentiable. In other words, $a_{1},a_{2},\dots ,a_{n}$ are real numbers and the point has coordinates $x_{1}=a_{1},x_{2}=a_{2},\dots ,x_{n}=a_{n}$. Suppose, further, that all the second-order partials (pure and mixed) of $f$ with respect to these variables exist at the point $(a_{1},a_{2},\dots ,a_{n})$. Then, the Hessian matrix of $f$ at $(a_{1},a_{2},\dots ,a_{n})$, denoted $H(f)(a_{1},a_{2},\dots ,a_{n})$, is a $n\times n$ matrix of real numbers that can be expressed explicitly as follows:

The $(ij)^{th}$ entry (i.e., the entry in the $i^{th}$ row and $j^{th}$ column) is $f_{x_{i}x_{j}}(a_{1},a_{2},\dots ,a_{n})$. This is the same as ${\frac {\partial ^{2}}{\partial x_{j}\partial x_{i}}}f(x_{1},x_{2},\dots ,x_{n})|_{(x_{1},x_{2},\dots ,x_{n})=(a_{1},a_{2},\dots ,a_{n})}$. Note that in the two notations, the order in which we write the partials differs because the convention differs (left-to-right versus right-to-left).

The matrix looks like this:

${\begin{pmatrix}f_{x_{1}x_{1}}(a_{1},a_{2},\dots ,a_{n})&f_{x_{1}x_{2}}(a_{1},a_{2},\dots ,a_{n})&\dots &f_{x_{1}x_{n}}(a_{1},a_{2},\dots ,a_{n})\\f_{x_{2}x_{1}}(a_{1},a_{2},\dots ,a_{n})&f_{x_{2}x_{2}}(a_{1},a_{2},\dots ,a_{n})&\dots &f_{x_{2}x_{n}}(a_{1},a_{2},\dots ,a_{n})\\\cdot &\cdot &\cdot &\cdot \\\cdot &\cdot &\cdot &\cdot \\\cdot &\cdot &\cdot &\cdot \\f_{x_{n}x_{1}}(a_{1},a_{2},\dots ,a_{n})&f_{x_{n}x_{2}}(a_{1},a_{2},\dots ,a_{n})&\dots &f_{x_{n}x_{n}}(a_{1},a_{2},\dots ,a_{n})\\\end{pmatrix}}$

{{#widget:YouTube|id=FIRMeFAeYqc}}

Definition as a function

For a function of two variables

Suppose $f$ is a real-valued function of two variables $x,y$. The Hessian matrix of $f$, denoted $H(f)$, is a $2\times 2$ matrix-valued function that sends each point to the Hessian matrix at that point, if that matrix is defined. It is defined as:

$(x_{0},y_{0})\mapsto H(f)(x_{0},y_{0})={\begin{pmatrix}f_{xx}(x_{0},y_{0})&f_{xy}(x_{0},y_{0})\\f_{yx}(x_{0},y_{0})&f_{yy}(x_{0},y_{0})\\\end{pmatrix}}$

In the point-free notation, we can write this as:

$H(f)={\begin{pmatrix}f_{xx}&f_{xy}\\f_{yx}&f_{yy}\\\end{pmatrix}}$

{{#widget:YouTube|id=39VO16DieuQ}}

For a function of multiple variables

Suppose $f$ is a function of variables $x_{1},x_{2},\dots ,x_{n}$. The Hessian matrix of $f$, denoted $H(f)$, is a $n\times n$ matrix-valued function that sends each point to the Hessian matrix at that point, if the matrix is defined. It is defined as:

$(a_{1},a_{2},\dots ,a_{n})\mapsto H(f)(a_{1},a_{2},\dots ,a_{n})$

In the point-free notation, we can write it as:

${\begin{pmatrix}f_{x_{1}x_{1}}&f_{x_{1}x_{2}}&\dots &f_{x_{1}x_{n}}\\f_{x_{2}x_{1}}&f_{x_{2}x_{2}}&\dots &f_{x_{2}x_{n}}\\\cdot &\cdot &\cdot &\cdot \\\cdot &\cdot &\cdot &\cdot \\\cdot &\cdot &\cdot &\cdot \\f_{x_{n}x_{1}}&f_{x_{n}x_{2}}&\dots &f_{x_{n}x_{n}}\\\end{pmatrix}}$

{{#widget:YouTube|id=DeFoV-NfjQQ}}

Under continuity assumptions

If we assume that all the second-order partials of $f$ are continuous functions everywhere, then the following happens:

• The Hessian matrix of $f$ at any point is a symmetric matrix, i.e., its $(ij)^{th}$ entry equals its $(ji)^{th}$ entry. This follows from Clairaut's theorem on equality of mixed partials.
• We can think of the Hessian matrix as the second derivative of the function, i.e., it is a matrix describing the second derivative.
• $f$ is twice differentiable as a function. Hence, the Hessian matrix of $f$ is the same as the Jacobian matrix of the gradient vector $\nabla f$, where the latter is viewed as a vector-valued function.

Note that the final conclusion actually only requires the existence of the gradient vector, hence it holds even if the second-order partials are not continuous.