Articulo de referencia

Lema del determinante de matrices

En matemáticas , en particular en álgebra lineal , el lema del determinante de una matriz calcula el determinante de la suma de una matriz invertible A y el producto diádico , u...

En matemáticas , en particular en álgebra lineal , el lema del determinante de una matriz calcula el determinante de la suma de una matriz invertible A y el producto diádico , u v T , de un vector columna u y un vector fila v T . [ 1 ] [ 2 ]

Declaración

Supongamos que A es una matriz cuadrada invertible y u y v son vectores columna. Entonces, el lema del determinante de matrices establece que

det(A+vT)=(1+vTA1)det(A).{\displaystyle \det(\mathbf {A} +\mathbf {uv} ^{\textsf {T}})=(1+\mathbf {v} ^{\textsf {T}}\mathbf {A} ^{-1}\mathbf {u} )\,\det(\mathbf {A} )\,.}

Aquí, uv T es el producto exterior de dos vectores u y v .

El teorema también puede enunciarse en términos de la matriz adjunta de A :

det(A+vT)=det(A)+vTadj(A),{\displaystyle \det(\mathbf {A} +\mathbf {uv} ^{\textsf {T}})=\det(\mathbf {A} )+\mathbf {v} ^{\textsf {T}}\mathrm {adj} (\mathbf {A} )\mathbf {u} \,,}

en cuyo caso se aplica independientemente de si la matriz A es invertible o no.

Prueba

Primero, la demostración del caso especial A = I se deduce de la igualdad: [ 3 ]

(I0vT1)(I+vT01)(I0vT1)=(I01+vT).{\displaystyle {\begin{pmatrix}\mathbf {I} &0\\\mathbf {v} ^{\textsf {T}}&1\end{pmatrix}}{\begin{pmatrix}\mathbf {I} +\mathbf {uv} ^{\textsf {T}}&\mathbf {u} \\0&1\end{pmatrix}}{\begin{pmatrix}\mathbf {I} &0\\-\mathbf {v} ^{\textsf {T}}&1\end{pmatrix}}={\begin{pmatrix}\mathbf {I} &\mathbf {u} \\0&1+\mathbf {v} ^{\textsf {T}}\mathbf {u} \end{pmatrix}}.}

El determinante del lado izquierdo es el producto de los determinantes de las tres matrices. Dado que la primera y la tercera matriz son triangulares con diagonal unitaria, sus determinantes son simplemente 1. El determinante de la matriz central es el valor que buscamos. El determinante del lado derecho es simplemente (1 + v T u ). Por lo tanto, tenemos el resultado:

det(I+vT)=1+vT.{\displaystyle \det(\mathbf {I} +\mathbf {uv} ^{\textsf {T}})=1+\mathbf {v} ^{\textsf {T}}\mathbf {u} .}

Entonces, el caso general se puede encontrar estableciendo u como A −1 u :

det(A+vT)=det(A)det(I+(A1)vT)=det(A)(1+vT(A1)).{\displaystyle {\begin{aligned}\det(\mathbf {A} +\mathbf {uv} ^{\textsf {T}})&=\det(\mathbf {A} )\det(\mathbf {I} +(\mathbf {A} ^{-1}\mathbf {u} )\mathbf {v} ^{\textsf {T}})\\&=\det(\mathbf {A} )\left(1+\mathbf {v} ^{\textsf {T}}(\mathbf {A} ^{-1}\mathbf {u} )\right).\end{aligned}}}

Solicitud

Si el determinante y la inversa de A ya se conocen, la fórmula proporciona una forma numéricamente económica de calcular el determinante de A corregido por la matriz uv T. El cálculo es relativamente económico porque no es necesario calcular el determinante de A + uv T desde cero (lo cual, en general, es costoso). Utilizando vectores unitarios para u y/o v , se pueden manipular columnas, filas o elementos individuales [ 4 ] de A y calcular un determinante actualizado correspondiente de forma relativamente económica.

Cuando se utiliza el lema del determinante de matrices junto con la fórmula de Sherman-Morrison , tanto la inversa como el determinante pueden actualizarse convenientemente de forma conjunta.

Generalización

Supongamos que A es una matriz invertible de n x n y U y V son matrices de n x m . Entonces

det(A+UVT)=det(Imetro+VTA1U)det(A).{\displaystyle \det(\mathbf {A} +\mathbf {UV} ^{\textsf {T}})=\det(\mathbf {I_{m}} +\mathbf {V} ^{\textsf {T}}\mathbf {A} ^{-1}\mathbf {U} )\det(\mathbf {A} ).}

En el caso especialA=Inorte{\displaystyle \mathbf {A} =\mathbf {I_ {n}}}Esta es la identidad Weinstein-Aronszajn .

Dada además una matriz invertible m x m W , la relación también puede expresarse como

det(A+UWVT)=det(W1+VTA1U)det(W)det(A).{\displaystyle \det(\mathbf {A} +\mathbf {UWV} ^{\textsf {T}})=\det(\mathbf {W} ^{-1}+\mathbf {V} ^{\textsf {T}}\mathbf {A} ^{-1}\mathbf {U} )\det(\mathbf {W} )\det(\mathbf {A} ).}

Véase también

Referencias

  1. Harville, DA (1997). Álgebra matricial desde la perspectiva de un estadístico . Nueva York: Springer-Verlag. ISBN 0-387-94978-X.
  2. Brookes, M. (2005). "The Matrix Reference Manual (online)" .
  3. Ding, J.; Zhou, A. (2007). "Autovalores de matrices actualizadas de rango uno con algunas aplicaciones" . Applied Mathematics Letters . 20 (12): 1223– 1226. doi : 10.1016/j.aml.2006.11.016 . ISSN 0893-9659 . 
  4. ↑ William H. Press; Brian P. Flannery ; Saul A. Teukolsky; William T. Vetterling (1992). Numerical Recipes in C: The Art of Scientific Computing . Cambridge University Press. pp. 73. ISBN  0-521-43108-5.