What is it?
For , the matrix of all partial derivatives : the best linear approximation of near a point. Its determinant measures how stretches volumes.
Why does it exist?
When both input and output are vectors, "the derivative" has to say how each output responds to each input: an table of rates. Arranged as a matrix it composes by multiplication (chain rule), inverts when is locally invertible, and its determinant gives the change-of-variables factor in integrals.
Intuition
Zoom in on a smooth map until it looks linear: the Jacobian is that linear map. A tiny square at is mapped to a tiny parallelogram spanned by the columns of ; its area is times the square's. For a robot arm, column is how fast the hand moves when joint turns.
Formal definition
Inverse function theorem: if and , is invertible near with .
Formulas
- robot velocity kinematics
- densities under an invertible map (normalizing flows)
Example
Polar coordinates : , . That is the familiar factor in .
Why does it matter?
Robotics lives on Jacobians (velocity kinematics, inverse kinematics, singularities where ). Newton's method for systems inverts one per step. Normalizing flows need cheap . Autodiff frameworks never build Jacobians explicitly: they compute Jacobian–vector (JVP) and vector–Jacobian (VJP) products.
Where it shows up in computing
Maps joint velocities to end-effector velocities; its singularities are configurations where the arm loses mobility.
IK solvers iterate (pseudo-inverse) or damped least squares.
The extended Kalman filter propagates covariances with the Jacobians of the dynamics and measurement models.
Where it shows up in AI
Normalizing flows compute exact likelihoods via , designing layers with cheap determinants.
AD computes JVPs (forward mode) and VJPs (reverse mode) without materializing .
Where is it used?
Computing topics reachable from here, through the chain of ideas that leads to them:
3D Computer graphics
- Multiple integrals and change of variables→Monte Carlo methods→The rendering equation★★★★★
- Equilibria and stability→Attractors→Chaos and sensitivity to initial conditions→Fractals→Procedural generation and noise★★★★★
- Multiple integrals and change of variables→Green's theorem→Mesh processing (discrete differential geometry)★★★★★
⚛ Physics and simulation
- Multiple integrals and change of variables→Surface integrals and flux→Divergence theorem (Gauss)→Electromagnetism (Maxwell's equations)★★★★★
- Multiple integrals and change of variables→Surface integrals and flux→Divergence theorem (Gauss)→Fluid dynamics and CFD★★★★★
- Equilibria and stability→Attractors→Chaos and sensitivity to initial conditions→Weather and climate modelling★★★★★
- Equilibria and stability→Stiffness and implicit methods→Physics engines★★★★★
- Equilibria and stability→Stiffness and implicit methods→Physics engines→N-body gravitational simulation★★★★★
- Equilibria and stability→Bifurcations→Population and epidemic models★★★★★
What depends on it
Exercises
Compute the Jacobian of and its determinant. Where is not locally invertible?
Solution
, . Only at the origin ( is in complex form).
A 2-link planar arm with lengths has hand position . When is ?
Solution
: zero when or — arm fully stretched or folded. There the hand cannot move radially.