Matriz y determinante jacobianos

En cálculo vectorial, la matriz jacobiana de una función vectorial de varias variables es la matriz cuyos elementos son las derivadas parciales de primer orden de dicha función. Si esta matriz es cuadrada, su determinante se llama el determinante jacobiano.

Tanto la matriz jacobiana como el determinante jacobiano reciben su nombre en honor al matemático Carl Gustav Jacobi.

Definición[editar]

Sea $f:\mathbb {R} ^{n}\longrightarrow \mathbb {R} ^{m}$ una función cuyas derivadas parciales de primer orden existen en todo $\mathbb {R} ^{n}$ y denotemos $f_{1},f_{2},...,f_{m}$ a sus componentes escalares. Se define la matriz jacobiana de $f$ en un punto $x\in \mathbb {R} ^{n}$ como:
${\text{J}}_{f}(x)={\begin{pmatrix}\displaystyle {\frac {\partial f_{1}}{\partial x_{1}}}(x)&\displaystyle {\frac {\partial f_{1}}{\partial x_{2}}}(x)&\dots &\displaystyle {\frac {\partial f_{1}}{\partial x_{n}}}(x)\\\displaystyle {\frac {\partial f_{2}}{\partial x_{1}}}(x)&\displaystyle {\frac {\partial f_{2}}{\partial x_{2}}}(x)&\dots &\displaystyle {\frac {\partial f_{2}}{\partial x_{n}}}(x)\\\vdots &\vdots &\ddots &\vdots \\\displaystyle {\frac {\partial f_{m}}{\partial x_{1}}}(x)&\displaystyle {\frac {\partial f_{m}}{\partial x_{2}}}(x)&\dots &\displaystyle {\frac {\partial f_{m}}{\partial x_{n}}}(x)\\\end{pmatrix}}={\begin{pmatrix}\nabla f_{1}(x)\\\vdots \\\nabla f_{m}(x)\end{pmatrix}}$
donde $\nabla f_{i}$ es el gradiente de la $i$ -ésima componente escalar.

Cuando $m=n$ , la matriz Jacobiana es cuadrada. Su determinante es conocido como el determinante Jacobiano de $f$ .
Cuando $m=1$ , esto es, cuando $f:\mathbb {R} ^{n}\longrightarrow \mathbb {R}$ es un campo escalar, la matriz Jacobiana se reduce a un vector fila. Este vector fila con todas las derivadas parciales de primer orden es el gradiente de $f$ , es decir, ${\text{J}}_{f}(x)=\nabla f(x)$ .
Cuando $m=n=1$ , esto es, cuando $f:\mathbb {R} \longrightarrow \mathbb {R}$ es una función una variable, la matriz Jacobiana sólo tiene una entrada: la derivada de la función $f$ en el punto.

Propiamente deberíamos hablar más que de matriz jacobiana, de diferencial jacobiana o aplicación lineal jacobiana ya que la forma de la matriz dependerá de la base o coordenadas elegidas. Es decir, dadas dos bases diferentes la aplicación lineal jacobiana tendrá componentes diferentes aún tratándose del mismo objeto matemático. La propiedad básica de la "matriz" jacobiana es la siguiente, dada una aplicación cualquiera $\mathbf {F} :\mathbb {R} ^{n}\to \mathbb {R} ^{m}$ continua, es decir $\mathbf {F} \in {\mathcal {C}}^{(k)}(\mathbb {R} ^{n},\mathbb {R} ^{m})$ se dirá que es diferenciable si existe una aplicación lineal ${\boldsymbol {\lambda }}\in {\mathcal {L}}(\mathbb {R} ^{n},\mathbb {R} ^{m})$ tal que:

\lim _{\|\mathbf {x} -\mathbf {y} \|\to 0}{\frac {\|(\mathbf {F} (\mathbf {x} )-\mathbf {F} (\mathbf {y} ))-{\boldsymbol {\lambda }}(\mathbf {x} -\mathbf {y} )\|}{\|\mathbf {x} -\mathbf {y} \|}}=0

Si $\mathbf {p}$ es un punto en $\mathbb {R} ^{n}$ y $\mathbf {F}$ es diferenciable en $\mathbf {p}$ entonces su diferencial está dada por J_F(p). En este caso, la aplicación lineal descrita por J_F(p) es la mejor aproximación lineal de $\mathbf {F}$ cerca del punto $\mathbf {p}$ , de esta manera:

\mathbf {F} (\mathbf {x} )\approx \mathbf {F} (\mathbf {p} )+J_{\mathbf {F} }(\mathbf {p} )(\mathbf {x} -\mathbf {p} )

para $\mathbf {x}$ cerca de $\mathbf {p}$ . O con mayor precisión:

\lim _{\|\mathbf {x} -\mathbf {p} \|\to 0}{\frac {\|\mathbf {F} (\mathbf {x} )-\mathbf {F} (\mathbf {p} )-J_{\mathbf {F} }(\mathbf {p} )(\mathbf {x} -\mathbf {p} )\|}{\|\mathbf {x} -\mathbf {p} \|}}=0

En ciertos espacios vectoriales de dimensión no finita, formados por funciones, puede generalizarse el concepto de matriz jacobiana definiendo una aplicación lineal jacobiana.

Determinante jacobiano[editar]

Si $m=n$ entonces $\mathbf {F}$ es una función que va de $\mathbb {R} ^{n}$ a $\mathbb {R} ^{n}$ y en este caso la matriz jacobiana es una matriz cuadrada, por lo que podemos calcular su determinante, este es conocido como el determinante jacobiano. El determinante jacobiano en ocasiones es conocido simplemente como “el Jacobiano”.

El determinante jacobiano en un punto dado nos da información importante sobre el comportamiento de $\mathbf {F}$ cerca de ese punto. Una función continuamente diferenciable $\mathbf {F}$ es invertible cerca del punto $\mathbf {p} \in \mathbb {R} ^{n}$ si el determinante jacobiano en $\mathbf {p}$ es no nulo. Este es el teorema de la función inversa. Más aún, el valor absoluto del determinante en $\mathbf {p}$ nos da el factor con el cual $\mathbf {F}$ expande o contrae su volumen cerca de $\mathbf {p}$ .

Inversa[editar]

De acuerdo al teorema de la función inversa, la matriz inversa de la matriz Jacobiana de una función invertible es la matriz Jacobiana de la función inversa. Esto es, si el Jacobiano de una función $\mathbf {F} :\mathbb {R} ^{n}\to \mathbb {R} ^{n}$ es continua y no singular en el punto $\mathbf {p} \in \mathbb {R} ^{n}$ entonces $\mathbf {F}$ es invertible cuando se restringe a un entorno de $\mathbf {p}$ y

\mathbf {J} _{\mathbf {F} ^{-1}}\circ \mathbf {F} =\mathbf {J} _{\mathbf {F} }^{-1}

Si el determinante jacobiano es diferente de cero en un punto entonces la función es localmente invertible cerca de este punto, esto es, existe un entorno de este punto en el que la función es invertible.

Puntos críticos[editar]

Si $\mathbf {F} :\mathbb {R} ^{n}\to \mathbb {R} ^{m}$ es una función diferenciable, un punto crítico de $\mathbf {F}$ es un punto en el que el rango de la matriz jacobiana es no maximal.

En el caso en que $m=n=k$ , un punto es crítico si el determinante jacobiano es cero.

Aplicaciones de la matriz jacobiana[editar]

La matriz Jacobiana establece una relación entre dos sistemas de referencia (espacios vectoriales), de manera que mediante una transformación (el jacobiano o matriz jacobiana), se pueda expresar un vector de un primer sistema, a otro sistema. Algo así como un traductor de un espacio vectorial a otro.

En Robótica de Manipuladores, se usa para deducir posibles puntos de una configuración que traerán dificultades para el movimiento, llamados puntos singulares o singularidades.

En el análisis dinámico, los “jacobianos del eslabón”, que permiten relacionar posición de cada eslabón con su movimiento de translación y de rotación.

Ejemplos[editar]

Ejemplo 1[editar]

La matriz jacobiana de la función $\mathbf {F} :\mathbb {R} ^{3}\to \mathbb {R} ^{3}$ dada por

\mathbf {F} (x_{1},x_{2},x_{3})=(x_{1},5x_{3},4x_{2}^{2}-2x_{3})

cuyas funciones componentes son

{\begin{aligned}f_{1}&=x_{1}\\f_{2}&=5x_{3}\\f_{3}&=4x_{2}^{2}-2x_{3}\end{aligned}}

es

\mathbf {J} _{\mathbf {F} }(x_{1},x_{2},x_{3})={\begin{bmatrix}{\cfrac {\partial f_{1}}{\partial x_{1}}}&{\cfrac {\partial f_{1}}{\partial x_{2}}}&{\cfrac {\partial f_{1}}{\partial x_{3}}}\\{\cfrac {\partial f_{2}}{\partial x_{1}}}&{\cfrac {\partial f_{2}}{\partial x_{2}}}&{\cfrac {\partial f_{2}}{\partial x_{3}}}\\{\cfrac {\partial f_{3}}{\partial x_{1}}}&{\cfrac {\partial f_{3}}{\partial x_{2}}}&{\cfrac {\partial f_{3}}{\partial x_{3}}}\end{bmatrix}}={\begin{bmatrix}1&0&0\\0&0&5\\0&8x_{2}&-2\end{bmatrix}}

La matriz jacobiana NO siempre es cuadrada. Véase el siguiente ejemplo.

Ejemplo 2[editar]

Supóngase la función $\mathbf {F} :\mathbb {R} ^{3}\to \mathbb {R} ^{4}$ , cuyas funciones componentes son:

y_{1}=x_{1}

y_{2}=5x_{3}\,

y_{3}=4x_{2}^{2}-2x_{3}\,

y_{4}=x_{3}\operatorname {sen}(x_{1})

tiene asociada como matriz jacobiana

\mathbf {J} _{\mathbf {F} }(x_{1},x_{2},x_{3})={\begin{bmatrix}{\dfrac {\partial y_{1}}{\partial x_{1}}}&{\dfrac {\partial y_{1}}{\partial x_{2}}}&{\dfrac {\partial y_{1}}{\partial x_{3}}}\\[3pt]{\dfrac {\partial y_{2}}{\partial x_{1}}}&{\dfrac {\partial y_{2}}{\partial x_{2}}}&{\dfrac {\partial y_{2}}{\partial x_{3}}}\\[3pt]{\dfrac {\partial y_{3}}{\partial x_{1}}}&{\dfrac {\partial y_{3}}{\partial x_{2}}}&{\dfrac {\partial y_{3}}{\partial x_{3}}}\\[3pt]{\dfrac {\partial y_{4}}{\partial x_{1}}}&{\dfrac {\partial y_{4}}{\partial x_{2}}}&{\dfrac {\partial y_{4}}{\partial x_{3}}}\\\end{bmatrix}}={\begin{bmatrix}1&0&0\\0&0&5\\0&8x_{2}&-2\\x_{3}\cos x_{1}&0&\operatorname {sen} x_{1}\end{bmatrix}}.

Ejemplo 3: Transformación a coordenadas polares[editar]

La transformación de coordenadas polares $(r,\theta )$ a coordenadas cartesianas $(x,y)$ está dada por la función $F:\mathbb {R} ^{+}\times [0,2\pi )\longrightarrow \mathbb {R} ^{2},\quad F(r,\theta )=(r\cos(\theta ),r\operatorname {sen}(\theta ))$ .

La matriz jacobiana de $F$ está dada por:

{\text{J}}_{F}(r,\theta )={\begin{pmatrix}\cos(\theta )&-r\operatorname {sen}(\theta )\\\operatorname {sen}(\theta )&r\cos(\theta )\end{pmatrix}}

y su determinante, o el jacobiano, es pues:

{\begin{aligned}\det({\text{J}}_{F}(r,\theta ))&=r\cos ^{2}(\theta )+r\operatorname {sen} ^{2}(\theta )\\&=r[\cos ^{2}(\theta )+\operatorname {sen} ^{2}(\theta )]\\&=r\end{aligned}}

y esto puede ser utilizado para transformar integrales entre dos sistemas de coordenadas:

\iint _{F(A)}f(x,y)dxdy=\iint _{A}f(r\cos \theta ,r\operatorname {sen} \theta )|J_{F}|=\iint _{A}f(r\cos \theta ,r\operatorname {sen} \theta )|r|

Pero al cogerse siempre $r>0$ con las coordenadas cartesianas, tenemos que

$\iint _{F(A)}f(x,y)dxdy=\iint _{A}f(r\cos \theta ,r\operatorname {sen} \theta )rdrd\theta$

Ejemplo 4[editar]

El determinante jacobiano de la función $\mathbf {F} :\mathbb {R} ^{3}\to \mathbb {R} ^{3}$ dada por:

\mathbf {F} (x_{1},x_{2},x_{3})=(5x_{2},4x_{1}^{2}-2\operatorname {sen}(x_{2}x_{3}),x_{2}x_{3})

es:

{\begin{aligned}\det(\mathbf {J} _{\mathbf {F} }(x_{1},x_{2},x_{3}))&={\begin{vmatrix}0&5&0\\8x_{1}&-2x_{3}\cos(x_{2}x_{3})&-2x_{2}\cos(x_{2}x_{3})\\0&x_{3}&x_{2}\end{vmatrix}}\\&=-5\cdot {\begin{vmatrix}8x_{1}&-2x_{2}\cos(x_{2}&x_{3})\\0&x_{2}\end{vmatrix}}\\&=-40x_{1}x_{2}\end{aligned}}

El teorema de la función inversa garantiza que la función es localmente invertible en todo el dominio excepto quizá donde $x_{1}=0$ o $x_{2}=0$ (es decir, los valores para los que el determinante se hace cero). Si imaginamos un objeto pequeño centrado en el punto (1,1,1) y le aplicamos F, tendremos un objeto aproximadamente 40 veces más voluminoso que el original.