Matriz jacobiana y determinante

Parte de una serie de artículos sobre
Cálculo
${\displaystyle \int _{a}^{b}f'(t)\,dt=f(b)-f(a)}$
Teorema fundamental Límites Continuidad Teorema de Rolle Teorema del valor medio Teorema de la función inversa
Diferencial Definiciones Derivada  ( generalizaciones ) Diferencial infinitesimal de una función total Conceptos Notación de diferenciación Segunda derivada Diferenciación implícita Diferenciación logarítmica Tarifas relacionadas Teorema de Taylor Reglas e identidades Suma Producto Cadena Fuerza Cociente La regla del Hospital Inverso El general Leibniz La fórmula de Faà di Bruno Reynolds
Integral Listas de integrales Transformada integral Regla integral de Leibniz Definiciones Antiderivada Integral  ( impropia ) Integral de Riemann Integración de Lebesgue Integración de contornos Integral de funciones inversas Integración por Regiones Discos Conchas cilíndricas Sustitución  ( trigonométrica , tangente de medio ángulo , Euler ) Fórmula de Euler Fracciones parciales ( método de Heaviside ) Cambiar orden Fórmulas de reducción Derivando bajo el signo integral Algoritmo de riesgo
Serie Geométrico  ( aritmético-geométrico ) Armónico Alterno Fuerza Binomio Taylor Pruebas de convergencia Límite de sumando (prueba de término) Relación Raíz Integral Comparación directa Comparación de límites Serie alternada Condensación de Cauchy Dirichlet Abel
Vector Gradiente Divergencia Rizo Laplaciano Derivada direccional Identidades Teoremas Gradiente Verduras De Stokes Divergencia Stokes generalizado Descomposición de Helmholtz
Multivariable Formalismos Matriz Tensor Exterior Geométrico Definiciones Derivada parcial Integral múltiple Integral de línea Integral de superficie Integral de volumen Jacobiano arpillera
Avanzado Cálculo sobre el espacio euclidiano Funciones generalizadas Límite de distribuciones
Especializado Fraccionario Malliavin Estocástico Variaciones
Miscelánea Precálculo Historia Glosario Lista de temas Abeja de integración Análisis matemático Análisis no estándar
en a mi

En cálculo vectorial , la matriz jacobiana ( / dʒ ə ˈ k oʊ b i ə n / , [ ^{1] [2] [3]} / dʒ ɪ -, j ɪ -/ ) de una función vectorial de varias variables es la matriz de todas sus derivadas parciales de primer orden . Cuando esta matriz es cuadrada , es decir, cuando la función toma como entrada el mismo número de variables que el número de componentes vectoriales de su salida, su determinante se denomina determinante jacobiano . Tanto la matriz como (si corresponde) el determinante se denominan a menudo simplemente jacobiano en la literatura. ^[4] Se denominan así en honor a Carl Gustav Jacob Jacobi .

Supóngase que f  : R ⁿ → R ^m es una función tal que cada una de sus derivadas parciales de primer orden existe en R ⁿ . Esta función toma un punto x ∈ R ⁿ como entrada y produce el vector f ( x ) ∈ R ^m como salida. Entonces la matriz jacobiana de f , denotada J _f ∈ R ^{m × n} , se define de manera que su ( i , j ) ^ésima entrada es , o explícitamente donde es la transpuesta (vector fila) del gradiente del -ésimo componente.${\textstyle {\frac {\partial f_{i}}{\partial x_{j}}}}$$${\displaystyle \mathbf {J_{f}} ={\begin{bmatrix}{\dfrac {\partial \mathbf {f} }{\partial x_{1}}}&\cdots &{\dfrac {\partial \mathbf {f} }{\partial x_{n}}}\end{bmatrix}}={\begin{bmatrix}\nabla ^{\mathrm {T} }f_{1}\\\vdots \\\nabla ^{\mathrm {T} }f_{m}\end{bmatrix}}={\begin{bmatrix}{\dfrac {\partial f_{1}}{\partial x_{1}}}&\cdots &{\dfrac {\partial f_{1}}{\partial x_{n}}}\\\vdots &\ddots &\vdots \\{\dfrac {\partial f_{m}}{\partial x_{1}}}&\cdots &{\dfrac {\partial f_{m}}{\partial x_{n}}}\end{bmatrix}}}$$${\displaystyle \nabla ^{\mathrm {T} }f_{i}}$${\displaystyle i}$

La matriz jacobiana, cuyas entradas son funciones de x , se denota de varias maneras; otras notaciones comunes incluyen D f , , y . ^{[5] [6]} Algunos autores definen la jacobiana como la transpuesta de la forma dada anteriormente.${\displaystyle \nabla \mathbf {f} }$${\displaystyle {\frac {\partial (f_{1},\ldots ,f_{m})}{\partial (x_{1},\ldots ,x_{n})}}}$

La matriz jacobiana representa la diferencial de f en cada punto donde f es diferenciable. En detalle, si h es un vector de desplazamiento representado por una matriz columna , el producto matricial J ( x ) ⋅ h es otro vector de desplazamiento, que es la mejor aproximación lineal del cambio de f en un entorno de x , si f ( x ) es diferenciable en x . ^[a] Esto significa que la función que asigna y a f ( x ) + J ( x ) ⋅ ( y – x ) es la mejor aproximación lineal de f ( y ) para todos los puntos y cercanos a x . La función lineal h → J ( x ) ⋅ h se conoce como la derivada o la diferencial de f en x .

Cuando m = n , la matriz jacobiana es cuadrada, por lo que su determinante es una función bien definida de x , conocida como determinante jacobiano de f . Contiene información importante sobre el comportamiento local de f . En particular, la función f tiene una función inversa diferenciable en un entorno de un punto x si y solo si el determinante jacobiano es distinto de cero en x (véase la conjetura jacobiana para un problema relacionado de invertibilidad global ). El determinante jacobiano también aparece cuando se cambian las variables en integrales múltiples (véase la regla de sustitución para variables múltiples ).

Cuando m = 1 , es decir, cuando f  : R ⁿ → R es una función escalar , la matriz jacobiana se reduce al vector fila ; este vector fila de todas las derivadas parciales de primer orden de f es la transpuesta del gradiente de f , es decir . Especializándonos aún más, cuando m = n = 1 , es decir, cuando f  : R → R es una función escalar de una sola variable, la matriz jacobiana tiene una sola entrada; esta entrada es la derivada de la función f .${\displaystyle \nabla ^{\mathrm {T} }f}$${\displaystyle \mathbf {J} _{f}=\nabla ^{T}f}$

Estos conceptos llevan el nombre del matemático Carl Gustav Jacob Jacobi (1804-1851).

La matriz jacobiana de una función vectorial de varias variables generaliza el gradiente de una función escalar de varias variables, que a su vez generaliza la derivada de una función escalar de una sola variable. En otras palabras, la matriz jacobiana de una función escalar de varias variables es (la transpuesta de) su gradiente y el gradiente de una función escalar de una sola variable es su derivada.

En cada punto en el que una función es diferenciable, su matriz jacobiana también puede considerarse como la descripción de la cantidad de "estiramiento", "rotación" o "transformación" que la función impone localmente cerca de ese punto. Por ejemplo, si ( x ′, y ′) = f ( x , y ) se utiliza para transformar suavemente una imagen, la matriz jacobiana J _f ( x , y ) , describe cómo se transforma la imagen en la vecindad de ( x , y ) .

Si una función es diferenciable en un punto, su diferencial está dada en coordenadas por la matriz jacobiana. Sin embargo, una función no necesita ser diferenciable para que su matriz jacobiana esté definida, ya que solo se requiere que existan sus derivadas parciales de primer orden .

Si f es diferenciable en un punto p en R ⁿ , entonces su diferencial está representada por J _f ( p ) . En este caso, la transformación lineal representada por J _f ( p ) es la mejor aproximación lineal de f cerca del punto p , en el sentido de que

$${\displaystyle \mathbf {f} (\mathbf {x} )-\mathbf {f} (\mathbf {p} )=\mathbf {J} _{\mathbf {f} }(\mathbf {p} )(\mathbf {x} -\mathbf {p} )+o(\|\mathbf {x} -\mathbf {p} \|)\quad ({\text{as }}\mathbf {x} \to \mathbf {p} ),}$$

donde o (‖ x − p ‖) es una cantidad que se acerca a cero mucho más rápido que la distancia entre x y p a medida que x se acerca a p . Esta aproximación se especializa en la aproximación de una función escalar de una sola variable por su polinomio de Taylor de grado uno, es decir

$${\displaystyle f(x)-f(p)=f'(p)(x-p)+o(x-p)\quad ({\text{as }}x\to p).}$$

En este sentido, el jacobiano puede considerarse como una especie de " derivada de primer orden " de una función vectorial de varias variables. En particular, esto significa que el gradiente de una función escalar de varias variables también puede considerarse como su "derivada de primer orden".

Las funciones diferenciables componibles f  : R ⁿ → R ^m y g  : R ^m → R ^k satisfacen la regla de la cadena , es decir, para x en R ⁿ .${\displaystyle \mathbf {J} _{\mathbf {g} \circ \mathbf {f} }(\mathbf {x} )=\mathbf {J} _{\mathbf {g} }(\mathbf {f} (\mathbf {x} ))\mathbf {J} _{\mathbf {f} }(\mathbf {x} )}$

El jacobiano del gradiente de una función escalar de varias variables tiene un nombre especial: la matriz hessiana , que en cierto sentido es la " segunda derivada " de la función en cuestión.

Si m = n , entonces f es una función de R ⁿ sobre sí misma y la matriz jacobiana es una matriz cuadrada . Podemos entonces formar su determinante , conocido como determinante jacobiano . El determinante jacobiano a veces se denomina simplemente "el jacobiano".

El determinante jacobiano en un punto dado proporciona información importante sobre el comportamiento de f cerca de ese punto. Por ejemplo, la función continuamente diferenciable f es invertible cerca de un punto p ∈ R ⁿ si el determinante jacobiano en p no es cero. Este es el teorema de la función inversa . Además, si el determinante jacobiano en p es positivo , entonces f conserva la orientación cerca de p ; si es negativo , f invierte la orientación. El valor absoluto del determinante jacobiano en p nos da el factor por el cual la función f expande o contrae volúmenes cerca de p ; por eso aparece en la regla general de sustitución .

El determinante jacobiano se utiliza cuando se realiza un cambio de variables al evaluar una integral múltiple de una función sobre una región dentro de su dominio. Para tener en cuenta el cambio de coordenadas, la magnitud del determinante jacobiano surge como un factor multiplicativo dentro de la integral. Esto se debe a que el elemento dV de n dimensiones es, en general, un paralelepípedo en el nuevo sistema de coordenadas, y el volumen n de un paralelepípedo es el determinante de sus vectores de arista.

El jacobiano también se puede utilizar para determinar la estabilidad de los equilibrios de sistemas de ecuaciones diferenciales aproximando el comportamiento cerca de un punto de equilibrio.

De acuerdo con el teorema de la función inversa , la matriz inversa de la matriz jacobiana de una función invertible f  : R ⁿ → R ⁿ es la matriz jacobiana de la función inversa . Es decir, la matriz jacobiana de la función inversa en un punto p es y el determinante jacobiano es$${\displaystyle \mathbf {J} _{\mathbf {f} ^{-1}}(\mathbf {p} )={\mathbf {J} _{\mathbf {f} }^{-1}(\mathbf {f} ^{-1}(\mathbf {p} ))},}$$

${\displaystyle \det(\mathbf {J} _{\mathbf {f} ^{-1}}(\mathbf {p} ))={\frac {1}{\det(\mathbf {J} _{\mathbf {f} }(\mathbf {f} ^{-1}(\mathbf {p} )))}}}$.

Si el jacobiano es continuo y no singular en el punto p en R ⁿ , entonces f es invertible cuando se restringe a algún entorno de p . En otras palabras, si el determinante jacobiano no es cero en un punto, entonces la función es localmente invertible cerca de este punto.

La conjetura jacobiana (no demostrada) está relacionada con la invertibilidad global en el caso de una función polinómica, es decir, una función definida por n polinomios en n variables. Afirma que, si el determinante jacobiano es una constante distinta de cero (o, equivalentemente, que no tiene ningún cero complejo), entonces la función es invertible y su inversa es una función polinómica.

Si f  : R ⁿ → R ^m es una función diferenciable , un punto crítico de f es un punto donde el rango de la matriz jacobiana no es máximo. Esto significa que el rango en el punto crítico es menor que el rango en algún punto vecino. En otras palabras, sea k la dimensión máxima de las bolas abiertas contenidas en la imagen de f ; entonces un punto es crítico si todos los menores de rango k de f son cero.

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

Consideremos una función f  : R ² → R ² , con ( x , y ) ↦ ( f ₁ ( x , y ), f ₂ ( x , y )), dada por$${\displaystyle \mathbf {f} \left({\begin{bmatrix}x\\y\end{bmatrix}}\right)={\begin{bmatrix}f_{1}(x,y)\\f_{2}(x,y)\end{bmatrix}}={\begin{bmatrix}x^{2}y\\5x+\sin y\end{bmatrix}}.}$$

Entonces tenemos y y la matriz jacobiana de f es y el determinante jacobiano es$${\displaystyle f_{1}(x,y)=x^{2}y}$$$${\displaystyle f_{2}(x,y)=5x+\sin y}$$$${\displaystyle \mathbf {J} _{\mathbf {f} }(x,y)={\begin{bmatrix}{\dfrac {\partial f_{1}}{\partial x}}&{\dfrac {\partial f_{1}}{\partial y}}\\[1em]{\dfrac {\partial f_{2}}{\partial x}}&{\dfrac {\partial f_{2}}{\partial y}}\end{bmatrix}}={\begin{bmatrix}2xy&x^{2}\\5&\cos y\end{bmatrix}}}$$$${\displaystyle \det(\mathbf {J} _{\mathbf {f} }(x,y))=2xy\cos y-5x^{2}.}$$

La transformación de coordenadas polares ( r , φ ) a coordenadas cartesianas ( x , y ), viene dada por la función F : R ⁺ × [0, 2 π ) → R ² con componentes:$${\displaystyle {\begin{aligned}x&=r\cos \varphi ;\\y&=r\sin \varphi .\end{aligned}}}$$$${\displaystyle \mathbf {J} _{\mathbf {F} }(r,\varphi )={\begin{bmatrix}{\frac {\partial x}{\partial r}}&{\frac {\partial x}{\partial \varphi }}\\[0.5ex]{\frac {\partial y}{\partial r}}&{\frac {\partial y}{\partial \varphi }}\end{bmatrix}}={\begin{bmatrix}\cos \varphi &-r\sin \varphi \\\sin \varphi &r\cos \varphi \end{bmatrix}}}$$

El determinante jacobiano es igual a r . Esto se puede utilizar para transformar integrales entre los dos sistemas de coordenadas:$${\displaystyle \iint _{\mathbf {F} (A)}f(x,y)\,dx\,dy=\iint _{A}f(r\cos \varphi ,r\sin \varphi )\,r\,dr\,d\varphi .}$$

La transformación de coordenadas esféricas ( ρ , φ , θ ) ^[7] a coordenadas cartesianas ( x , y , z ), viene dada por la función F : R ⁺ × [0, π ) × [0, 2 π ) → R ³ con componentes:$${\displaystyle {\begin{aligned}x&=\rho \sin \varphi \cos \theta ;\\y&=\rho \sin \varphi \sin \theta ;\\z&=\rho \cos \varphi .\end{aligned}}}$$

La matriz jacobiana para este cambio de coordenadas es El determinante es ρ ² sen φ . Dado que dV = dx dy dz es el volumen de un elemento de volumen diferencial rectangular (porque el volumen de un prisma rectangular es el producto de sus lados), podemos interpretar dV = ρ ² sen φ dρ dφ dθ como el volumen del elemento de volumen diferencial esférico . A diferencia del volumen del elemento de volumen diferencial rectangular, el volumen de este elemento de volumen diferencial no es una constante y varía con las coordenadas ( ρ y φ ). Puede utilizarse para transformar integrales entre los dos sistemas de coordenadas:$${\displaystyle \mathbf {J} _{\mathbf {F} }(\rho ,\varphi ,\theta )={\begin{bmatrix}{\dfrac {\partial x}{\partial \rho }}&{\dfrac {\partial x}{\partial \varphi }}&{\dfrac {\partial x}{\partial \theta }}\\[1em]{\dfrac {\partial y}{\partial \rho }}&{\dfrac {\partial y}{\partial \varphi }}&{\dfrac {\partial y}{\partial \theta }}\\[1em]{\dfrac {\partial z}{\partial \rho }}&{\dfrac {\partial z}{\partial \varphi }}&{\dfrac {\partial z}{\partial \theta }}\end{bmatrix}}={\begin{bmatrix}\sin \varphi \cos \theta &\rho \cos \varphi \cos \theta &-\rho \sin \varphi \sin \theta \\\sin \varphi \sin \theta &\rho \cos \varphi \sin \theta &\rho \sin \varphi \cos \theta \\\cos \varphi &-\rho \sin \varphi &0\end{bmatrix}}.}$$ $${\displaystyle \iiint _{\mathbf {F} (U)}f(x,y,z)\,dx\,dy\,dz=\iiint _{U}f(\rho \sin \varphi \cos \theta ,\rho \sin \varphi \sin \theta ,\rho \cos \varphi )\,\rho ^{2}\sin \varphi \,d\rho \,d\varphi \,d\theta .}$$

La matriz jacobiana de la función F  : R ³ → R ⁴ con componentes es$${\displaystyle {\begin{aligned}y_{1}&=x_{1}\\y_{2}&=5x_{3}\\y_{3}&=4x_{2}^{2}-2x_{3}\\y_{4}&=x_{3}\sin x_{1}\end{aligned}}}$$$${\displaystyle \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}}}\\[1em]{\dfrac {\partial y_{2}}{\partial x_{1}}}&{\dfrac {\partial y_{2}}{\partial x_{2}}}&{\dfrac {\partial y_{2}}{\partial x_{3}}}\\[1em]{\dfrac {\partial y_{3}}{\partial x_{1}}}&{\dfrac {\partial y_{3}}{\partial x_{2}}}&{\dfrac {\partial y_{3}}{\partial x_{3}}}\\[1em]{\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&\sin x_{1}\end{bmatrix}}.}$$

Este ejemplo muestra que la matriz jacobiana no necesita ser una matriz cuadrada.

El determinante jacobiano de la función F  : R ³ → R ³ con componentes es$${\displaystyle {\begin{aligned}y_{1}&=5x_{2}\\y_{2}&=4x_{1}^{2}-2\sin(x_{2}x_{3})\\y_{3}&=x_{2}x_{3}\end{aligned}}}$$$${\displaystyle {\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}}=-8x_{1}{\begin{vmatrix}5&0\\x_{3}&x_{2}\end{vmatrix}}=-40x_{1}x_{2}.}$$

De esto vemos que F invierte la orientación cerca de aquellos puntos donde x ₁ y x ₂ tienen el mismo signo; la función es localmente invertible en todas partes excepto cerca de los puntos donde x ₁ = 0 o x ₂ = 0. Intuitivamente, si uno comienza con un objeto diminuto alrededor del punto (1, 2, 3) y aplica F a ese objeto, obtendrá un objeto resultante con aproximadamente 40 × 1 × 2 = 80 veces el volumen del original, con la orientación invertida.

Consideremos un sistema dinámico de la forma , donde es la derivada (componente a componente) de con respecto al parámetro de evolución (tiempo), y es diferenciable. Si , entonces es un punto estacionario (también llamado estado estable ). Por el teorema de Hartman-Grobman , el comportamiento del sistema cerca de un punto estacionario está relacionado con los valores propios de , el jacobiano de en el punto estacionario. ^[8] Específicamente, si todos los valores propios tienen partes reales que son negativas, entonces el sistema es estable cerca del punto estacionario. Si cualquier valor propio tiene una parte real que es positiva, entonces el punto es inestable. Si la parte real más grande de los valores propios es cero, la matriz jacobiana no permite una evaluación de la estabilidad. ^[9]${\displaystyle {\dot {\mathbf {x} }}=F(\mathbf {x} )}$${\displaystyle {\dot {\mathbf {x} }}}$${\displaystyle \mathbf {x} }$ ${\displaystyle t}$${\displaystyle F\colon \mathbb {R} ^{n}\to \mathbb {R} ^{n}}$${\displaystyle F(\mathbf {x} _{0})=0}$${\displaystyle \mathbf {x} _{0}}$${\displaystyle \mathbf {J} _{F}\left(\mathbf {x} _{0}\right)}$${\displaystyle F}$

Un sistema cuadrado de ecuaciones no lineales acopladas se puede resolver iterativamente mediante el método de Newton . Este método utiliza la matriz jacobiana del sistema de ecuaciones.

El jacobiano sirve como matriz de diseño linealizada en regresión estadística y ajuste de curvas ; véase mínimos cuadrados no lineales . El jacobiano también se utiliza en matrices aleatorias, momentos, sensibilidad local y diagnósticos estadísticos. ^{[10] [11]}

• Colector central
• Matriz de Hesse
• Empuje hacia adelante (diferencial)

• Gandolfo, Giancarlo (1996). "Estática comparativa y el principio de correspondencia". Dinámica económica (tercera edición). Berlín: Springer. pp. 305–330. ISBN 3-540-60988-1.
• Protter, Murray H. ; Morrey, Charles B. Jr. (1985). "Transformaciones y jacobianos". Cálculo intermedio (segunda edición). Nueva York: Springer. pp. 412–420. ISBN 0-387-96058-9.

• "Jacobiano", Enciclopedia de Matemáticas , EMS Press , 2001 [1994]
• Mathworld Una explicación más técnica de los jacobianos