Jacobian matrix and determinant
Part of a series of articles about  
Calculus  





Specialized


In vector calculus, the Jacobian matrix (/dʒəˈkoʊbiən/,^{[1]}^{[2]}^{[3]} /dʒɪ,
Suppose f : ℝ^{n} → ℝ^{m} is a function each of whose firstorder partial derivatives exist on ℝ^{n}. This function takes a point x ∈ ℝ^{n} as input and produces the vector f(x) ∈ ℝ^{m} as output. Then the Jacobian matrix of f is defined to be an m×n matrix, denoted by J, whose (i,j)th entry is , or explicitly
This matrix, whose entries are functions of x, is also denoted variously by Df, J_{f}, and ∂(f_{1},...,f_{m})/∂(x_{1},...,x_{n}). (However, some literature defines the Jacobian as the transpose of the matrix given above.)
The Jacobian matrix represents the differential of f at every point where f is differentiable. In detail, with respect to a given point x∈ ℝ^{n}, the linear transformation represented by J(x) takes a position vector in ℝ^{n} from x as reference point as input and produces the position vector in ℝ^{m} from f(x) as reference point obtained by multiplying by J(x) as output. If f is differentiable at some point x^{[a]}, then this is the linear transformation that best approximates f for points close x, and is known as the derivative or the differential of f at x.
When m = n, the Jacobian matrix is square, so its determinant is a welldefined function of x, known as the Jacobian determinant of f. It carries important information about the local behavior of f. In particular, the function f has locally in the neighborhood of a point x an inverse function that is differentiable if and only if the Jacobian determinant is nonzero at x (see Jacobian conjecture). The Jacobian determinant also appears when changing the variables in multiple integrals (see substitution rule for multiple variables).
When m = 1, that is when f : ℝ^{n} → ℝ is a scalarvalued function, the Jacobian matrix reduces to a row vector. This row vector of all firstorder partial derivatives of f is the transpose of the gradient of f, i.e. . Here we are adopting the convention that the gradient vector is a column vector. Specialising further, when m = n = 1, that is when f : ℝ → ℝ is a scalarvalued function of a single variable, the Jacobian matrix has a single entry. This entry is the derivative of the function f.
These concepts are named after the mathematician Carl Gustav Jacob Jacobi (1804–1851).
Contents
Jacobian matrix
The Jacobian of a vectorvalued function in several variables generalizes the gradient of a scalarvalued function in several variables, which in turn generalizes the derivative of a scalarvalued function of a single variable. In other words, the Jacobian matrix of a scalarvalued function in several variables is (the transpose of) its gradient and the gradient of a scalarvalued function of a single variable is its derivative.
At each point where a function is differentiable, its Jacobian matrix can also be thought of as describing the amount of "stretching", "rotating" or "transforming" that the function imposes locally near that point. For example, if (x′, y′) = f(x, y) is used to smoothly transform an image, the Jacobian matrix J_{f}(x, y), describes how the image in the neighborhood of (x, y) is transformed.
If a function is differentiable at a point, its differential is given in coordinates by the Jacobian matrix. However a function does not need to be differentiable for its Jacobian matrix to be defined, since only its firstorder partial derivatives are required to exist.
If f is differentiable at a point p in ℝ^{n}, then its differential is represented by J_{f}(p). In this case, the linear transformation represented by J_{f}(p) is the best linear approximation of f near the point p, in the sense that
where o(‖x − p‖) is a quantity that approaches zero much faster than the distance between x and p does as x approaches p. This approximation specializes to the approximation of a scalar function of a single variable by its Taylor polynomial of degree one, namely
 .
In this sense, the Jacobian may be regarded as a kind of "firstorder derivative" of a vectorvalued function of several variables. In particular, this means that the gradient of a scalarvalued function of several variables may too be regarded as its "firstorder derivative".
Composable differentiable functions f : ℝ^{n} → ℝ^{m} and g : ℝ^{m} → ℝ^{k} satisfy the chain rule, namely for x in ℝ^{n}.
The Jacobian of the gradient of a scalar function of several variables has a special name: the Hessian matrix, which in a sense is the "second derivative" of the function in question.
Jacobian determinant
If m = n, then f is a function from ℝ^{n} to itself and the Jacobian matrix is a square matrix. We can then form its determinant, known as the Jacobian determinant. The Jacobian determinant is sometimes referred to as "the Jacobian".
The Jacobian determinant at a given point gives important information about the behavior of f near that point. For instance, the continuously differentiable function f is invertible near a point p ∈ ℝ^{n} if the Jacobian determinant at p is nonzero. This is the inverse function theorem. Furthermore, if the Jacobian determinant at p is positive, then f preserves orientation near p; if it is negative, f reverses orientation. The absolute value of the Jacobian determinant at p gives us the factor by which the function f expands or shrinks volumes near p; this is why it occurs in the general substitution rule.
The Jacobian determinant is used when making a change of variables when evaluating a multiple integral of a function over a region within its domain. To accommodate for the change of coordinates the magnitude of the Jacobian determinant arises as a multiplicative factor within the integral. This is because the ndimensional dV element is in general a parallelepiped in the new coordinate system, and the nvolume of a parallelepiped is the determinant of its edge vectors.
The Jacobian can also be used to solve systems of differential equations at an equilibrium point or approximate solutions near an equilibrium point. Its applications include determining the stability of the diseasefree equilibrium in disease modelling.^{[5]}
Inverse
According to the inverse function theorem, the matrix inverse of the Jacobian matrix of an invertible function is the Jacobian matrix of the inverse function. That is, if the Jacobian of the function f : ℝ^{n} → ℝ^{n} is continuous and nonsingular at the point p in ℝ^{n}, then f is invertible when restricted to some neighborhood of p and
Conversely, if the Jacobian determinant is not zero at a point, then the function is locally invertible near this point, that is, there is a neighbourhood of this point in which the function is invertible.
The (unproved) Jacobian conjecture is related to global invertibility in the case of a polynomial function, that is a function defined by n polynomials in n variables. It asserts that, if the Jacobian determinant is a nonzero constant (or, equivalently, that it does not have any complex zero), then the function is invertible and its inverse is a polynomial function.
Critical points
If f : ℝ^{n} → ℝ^{m} is a differentiable function, a critical point of f is a point where the rank of the Jacobian matrix is not maximal. This means that the rank at the critical point is lower than the rank at some neighbour point. In other words, let k be the maximal dimension of the open balls contained in the image of f; then a point is critical if all minors of rank k of f are zero.
In the case where 1 = m = n = k, a point is critical if the Jacobian determinant is zero.
Examples
Example 1
Consider the function f : ℝ^{2} → ℝ^{2}, with (x, y) ↦ (f_{1}(x, y), f_{2}(x, y)), given by
Then we have
and
and the Jacobian matrix of f is
and the Jacobian determinant is
Example 2: polarCartesian transformation
The transformation from polar coordinates (r, φ) to Cartesian coordinates (x, y), is given by the function F: ℝ^{+} × [0, 2π) → ℝ^{2} with components:
The Jacobian determinant is equal to r. This can be used to transform integrals between the two coordinate systems:
Example 3: sphericalCartesian transformation
The transformation from spherical coordinates (r, θ, φ) to Cartesian coordinates (x, y, z), is given by the function F: ℝ^{+} × [0, π] × [0, 2 π) → ℝ^{3} with components:
The Jacobian matrix for this coordinate change is
The determinant is r^{2} sin θ. As an example, since dV = dx dy dz this determinant implies that the differential volume element dV = r^{2} sin θ dr dθ dφ. Unlike for a change of Cartesian coordinates, this determinant is not a constant, and varies with coordinates (r and θ).
Example 4
The Jacobian matrix of the function F : ℝ^{3} → ℝ^{4} with components
is
This example shows that the Jacobian need not be a square matrix.
Example 5
The Jacobian determinant of the function F : ℝ^{3} → ℝ^{3} with components
is
From this we see that F reverses orientation near those points where x_{1} and x_{2} have the same sign; the function is locally invertible everywhere except near points where x_{1} = 0 or x_{2} = 0. Intuitively, if one starts with a tiny object around the point (1, 2, 3) and apply F to that object, one will get a resulting object with approximately 40 × 1 × 2 = 80 times the volume of the original one, with orientation reversed.
Other uses
The Jacobian serves as a linearized design matrix in statistical regression and curve fitting; see nonlinear least squares.
Dynamical systems
Consider a dynamical system of the form , where is the (componentwise) derivative of with respect to the evolution parameter (time), and is differentiable. If , then is a stationary point (also called a steady state). By the Hartman–Grobman theorem, the behavior of the system near a stationary point is related to the eigenvalues of , the Jacobian of at the stationary point.^{[6]} Specifically, if the eigenvalues all have real parts that are negative, then the system is stable near the stationary point, if any eigenvalue has a real part that is positive, then the point is unstable. If the largest real part of the eigenvalues is zero, the Jacobian matrix does not allow for an evaluation of the stability.^{[7]}
Newton's method
A square system of coupled nonlinear equations can be solved iteratively by Newton's method. This method uses the Jacobian matrix of the system of equations.
Surface analysis
Let n = 2 so the Jacobian is a 2 × 2 real matrix. Suppose a surface diffeomorphism f: U → V in the neighborhood of p in U is written The matrix can be interpreted as a complex number: ordinary, split, or dual. Furthermore, since is invertible, the complex number has a polar decomposition or an alternative planar decomposition.
And again, each such complex number represents a group action on the tangent plane at p. The action is dilation by the norm of the complex number, and rotation respecting angle, hyperbolic angle, or slope, according to the case of Such action corresponds to a conformal mapping.
See also
References
 ^ "Jacobian  Definition of Jacobian in English by Oxford Dictionaries". Oxford Dictionaries  English. Archived from the original on 1 December 2017. Retrieved 2 May 2018.
 ^ "the definition of jacobian". Dictionary.com. Archived from the original on 1 December 2017. Retrieved 2 May 2018.
 ^ Team, Forvo. "Jacobian pronunciation: How to pronounce Jacobian in English". forvo.com. Retrieved 2 May 2018.
 ^ W., Weisstein, Eric. "Jacobian". mathworld.wolfram.com. Archived from the original on 3 November 2017. Retrieved 2 May 2018.
 ^ Smith? RJ (2015). "The Joys of the Jacobian". Chalkdust. 2: 10–17.
 ^ Arrowsmith, D. K.; Place, C. M. (1992). "The Linearization Theorem". Dynamical Systems: Differential Equations, Maps, and Chaotic Behaviour. London: Chapman & Hall. pp. 77–81. ISBN 0412390809.
 ^ Hirsch, Morris; Smale, Stephen (1974). Differential equations, dynamical systems and linear algebra.
Further reading
 Gandolfo, Giancarlo (1996). Economic Dynamics (Third ed.). Berlin: Springer. pp. 305–330. ISBN 3540609881.
External links
 Hazewinkel, Michiel, ed. (2001) [1994], "Jacobian", Encyclopedia of Mathematics, Springer Science+Business Media B.V. / Kluwer Academic Publishers, ISBN 9781556080104
 Mathworld A more technical explanation of Jacobians
Cite error: There are <ref group=loweralpha>
tags or {{efn}}
templates on this page, but the references will not show without a {{reflistgroup=loweralpha}}
template or {{notelist}}
template (see the help page).