Articulo de referencia

Descomposición matricial

Diagrama que resume las relaciones entre las clases de matrices y las factorizaciones de matrices comunes. En la disciplina matemática del álgebra lineal , una descomposición ma...

Diagrama que resume las relaciones entre las clases de matrices y las factorizaciones de matrices comunes.

En la disciplina matemática del álgebra lineal , una descomposición matricial o factorización matricial consiste en factorizar una matriz en un producto de matrices. Existen diversas descomposiciones matriciales; cada una se utiliza en una clase particular de problemas.

Ejemplo

En el análisis numérico , se utilizan diferentes descomposiciones para implementar algoritmos matriciales eficientes .

Por ejemplo, al resolver un sistema de ecuaciones linealesAincógnita=b{\displaystyle A\mathbf {x} =\mathbf {b} }La matriz A se puede descomponer mediante la descomposición LU . La descomposición LU factoriza una matriz en una matriz triangular inferior L y una matriz triangular superior U. Los sistemasL(Uincógnita)=b{\displaystyle L(U\mathbf {x} )=\mathbf {b} }yUincógnita=L1b{\displaystyle U\mathbf {x} =L^{-1}\mathbf {b} }requieren menos sumas y multiplicaciones para resolverse, en comparación con el sistema original.Aincógnita=b{\displaystyle A\mathbf {x} =\mathbf {b} }, aunque en aritmética inexacta como la de punto flotante se podrían requerir muchos más dígitos .

De manera similar, la descomposición QR expresa A como QR con Q una matriz ortogonal y R una matriz triangular superior. El sistema Q ( R x ) = b se resuelve mediante R x = Q T b = c , y el sistema R x = c se resuelve mediante ' sustitución hacia atrás '. El número de sumas y multiplicaciones requeridas es aproximadamente el doble que el de usar el solucionador LU, pero no se requieren más dígitos en aritmética inexacta porque la descomposición QR es numéricamente estable .

descomposición LU

reducción de LU

Descomposición LU en bloques

factorización de rangos

Descomposición colérica

  • Aplicable a: matriz cuadrada , hermitiana y definida positivaA{\displaystyle A}
  • Descomposición:A=UU{\displaystyle A=U^{*}U}, dóndeU{\displaystyle U}es triangular superior con entradas diagonales reales positivas
  • Comentario: si la matrizA{\displaystyle A}es hermitiana y semidefinida positiva, entonces tiene una descomposición de la formaA=UU{\displaystyle A=U^{*}U}si las entradas diagonales deU{\displaystyle U}se permite que sean cero
  • Unicidad: para matrices definidas positivas, la descomposición de Cholesky es única. Sin embargo, no lo es en el caso de matrices semidefinidas positivas.
  • Comentario: siA{\displaystyle A}es real y simétrico,U{\displaystyle U}tiene todos los elementos reales
  • Comentario: Una alternativa es la descomposición LDL , que puede evitar la extracción de raíces cuadradas.

descomposición QR

  • Aplicable a: matriz A de m × n con columnas linealmente independientes.
  • Descomposición:A=QR{\displaystyle A=QR}dóndeQ{\displaystyle Q}es una matriz unitaria de tamaño m por m , yR{\displaystyle R}es una matriz triangular superior de tamaño m por n
  • Singularidad: En general no es único, pero siA{\displaystyle A}es de rango completo , entonces existe un únicoR{\displaystyle R}que tiene todos los elementos diagonales positivos. SiA{\displaystyle A}es cuadrado, tambiénQ{\displaystyle Q}es único.
  • Comentario: La descomposición QR proporciona una forma eficaz de resolver el sistema de ecuaciones.Aincógnita=b{\displaystyle A\mathbf {x} =\mathbf {b} }. El hecho de queQ{\displaystyle Q}es ortogonal significa queQTQ=I{\displaystyle Q^{\mathrm {T} }Q=I}, de modo queAincógnita=b{\displaystyle A\mathbf {x} =\mathbf {b} }es equivalente aRincógnita=QTb{\displaystyle R\mathbf {x} =Q^{\mathsf {T}}\mathbf {b} }, lo cual es muy fácil de resolver ya queR{\displaystyle R}es triangular .

Factorización RRQR

descomposición interpolativa

descomposición en valores propios

  • También llamada descomposición espectral .
  • Aplicable a: matriz cuadrada A con vectores propios linealmente independientes (no necesariamente valores propios distintos).
  • Descomposición:A=VDV1{\displaystyle A=VDV^{-1}}, donde D es una matriz diagonal formada a partir de los valores propios de A , y las columnas de V son los vectores propios correspondientes de A.
  • Existencia: Una matriz A de n × n siempre tiene n valores propios (complejos), que pueden ordenarse (de más de una manera) para formar una matriz diagonal D de n × n y una matriz correspondiente de columnas no nulas V que satisface la ecuación de valores propios.AV=VD{\displaystyle AV=VD}. V{\displaystyle V}Es invertible si y solo si los n autovectores son linealmente independientes (es decir, cada autovalor tiene una multiplicidad geométrica igual a su multiplicidad algebraica ). Una condición suficiente (pero no necesaria) para que esto ocurra es que todos los autovalores sean diferentes (en este caso, la multiplicidad geométrica y algebraica son iguales a 1).
  • Comentario: Siempre se pueden normalizar los autovectores para que tengan longitud uno (véase la definición de la ecuación de autovalores).
  • Comentario: Toda matriz normal A (es decir, matriz para la cualAA=AA{\displaystyle AA^{*}=A^{*}A}, dóndeA{\displaystyle A^{*}}es una transpuesta conjugada ) se puede descomponer en valores propios. Para una matriz normal A (y solo para una matriz normal), los vectores propios también se pueden hacer ortonormales (VV=I{\displaystyle VV^{*}=I}) y la descomposición en valores propios se lee comoA=VDV{\displaystyle A=VDV^{*}}En particular, todas las matrices unitarias , hermíticas o antihermíticas (en el caso de valores reales, todas las matrices ortogonales , simétricas o antisimétricas , respectivamente) son normales y, por lo tanto, poseen esta propiedad.
  • Comentario: Para cualquier matriz simétrica real A , la descomposición en valores propios siempre existe y se puede escribir comoA=VDVT{\displaystyle A=VDV^{\mathsf {T}}}donde tanto D como V son variables de valor real.
  • Comentario: La descomposición en valores propios es útil para comprender la solución de un sistema de ecuaciones diferenciales ordinarias lineales o ecuaciones en diferencias lineales. Por ejemplo, la ecuación en diferenciasincógnitat+1=Aincógnitat{\displaystyle x_{t+1}=Ax_{t}}partiendo de la condición inicialincógnita0=do{\displaystyle x_{0}=c}se resuelve medianteincógnitat=Atdo{\displaystyle x_{t}=A^{t}c}, lo cual es equivalente aincógnitat=VDtV1do{\displaystyle x_{t}=VD^{t}V^{-1}c}donde V y D son las matrices formadas a partir de los autovectores y autovalores de A. Dado que D es diagonal, al elevarla a la potenciaDt{\displaystyle D^{t}}, simplemente implica elevar cada elemento de la diagonal a la potencia t . Esto es mucho más fácil de hacer y comprender que elevar A a la potencia t , ya que A generalmente no está en la diagonal.

descomposición de Jordan

La forma normal de Jordan y la descomposición de Jordan-Chevalley

  • Aplicable a: matriz cuadrada A
  • Comentario: la forma normal de Jordan generaliza la descomposición en valores propios a casos donde hay valores propios repetidos y no se pueden diagonalizar; la descomposición de Jordan-Chevalley hace esto sin elegir una base.

Descomposición de Schur

Descomposición de Schur real

  • Aplicable a: matriz cuadrada A
  • Descomposición: Esta es una versión de la descomposición de Schur dondeV{\displaystyle V}yS{\displaystyle S}solo contienen números reales. Siempre se puede escribirA=VSVT{\displaystyle A=VSV^{\mathsf {T}}}donde V es una matriz ortogonal real ,VT{\displaystyle V^{\mathsf {T}}}es la transpuesta de V , y S es una matriz triangular superior por bloques llamada forma de Schur real . Los bloques en la diagonal de S son de tamaño 1×1 (en cuyo caso representan valores propios reales) o 2×2 (en cuyo caso se derivan de pares de valores propios conjugados complejos ).

Descomposición de QZ

  • También llamada: descomposición de Schur generalizada
  • Aplicable a: matrices cuadradas A y B
  • Comentario: existen dos versiones de esta descomposición: compleja y real.
  • Descomposición (versión compleja):A=QSZ{\displaystyle A=QSZ^{*}}yB=QTZ{\displaystyle B=QTZ^{*}}donde Q y Z son matrices unitarias , el superíndice * representa la transpuesta conjugada , y S y T son matrices triangulares superiores .
  • Comentario: en la descomposición QZ compleja, las razones de los elementos diagonales de S a los elementos diagonales correspondientes de T ,λi=Sii/Tii{\displaystyle \lambda _{i}=S_{ii}/T_{ii}}, son los autovalores generalizados que resuelven el problema generalizado de autovalores.Av=λBv{\displaystyle A\mathbf {v} =\lambda B\mathbf {v} }(dóndeλ{\displaystyle \lambda }es un escalar desconocido y v es un vector no nulo desconocido).
  • Descomposición (versión real):A=QSZT{\displaystyle A=QSZ^{\mathsf {T}}}yB=QTZT{\displaystyle B=QTZ^{\mathsf {T}}}donde A , B , Q , Z , S y T son matrices que contienen solo números reales. En este caso, Q y Z son matrices ortogonales , el superíndice T representa la transposición , y S y T son matrices triangulares superiores por bloques . Los bloques en la diagonal de S y T son de tamaño 1×1 o 2×2.

Factorización de Takagi

  • Aplicable a: matriz cuadrada, compleja y simétrica A.
  • Descomposición:A=VDVT{\displaystyle A=VDV^{\mathsf {T}}}donde D es una matriz diagonal real no negativa y V es unitaria .VT{\displaystyle V^{\mathsf {T}}}denota la transpuesta de la matriz V.
  • Comentario: Los elementos diagonales de D son las raíces cuadradas no negativas de los valores propios deAA=VD2V1{\displaystyle AA^{*}=VD^{2}V^{-1}}.
  • Comentario: V puede ser complejo incluso si A es real.
  • Comentario: Este no es un caso especial de la descomposición en valores propios (ver arriba), que utilizaV1{\displaystyle V^{-1}}en lugar deVT{\displaystyle V^{\mathsf {T}}}Además, si A no es real, no es hermitiano y la forma que utilizaV{\displaystyle V^{*}}Tampoco aplica.

Descomposición en valores singulares

  • Aplicable a: matriz A de m por n .
  • Descomposición:A=UDV{\displaystyle A=UDV^{*}}donde D es una matriz diagonal no negativa y U y V satisfacenUU=I,VV=I{\displaystyle U^{*}U=I,V^{*}V=I}. AquíV{\displaystyle V^{*}}es la transpuesta conjugada de V (o simplemente la transpuesta , si V contiene solo números reales), e I denota la matriz identidad (de alguna dimensión).
  • Comentario: Los elementos diagonales de D se denominan valores singulares de A.
  • Comentario: Al igual que la descomposición en valores propios descrita anteriormente, la descomposición en valores singulares implica encontrar direcciones base a lo largo de las cuales la multiplicación de matrices es equivalente a la multiplicación escalar, pero tiene mayor generalidad ya que la matriz en consideración no tiene por qué ser cuadrada.
  • Unicidad: los valores singulares deA{\displaystyle A}Siempre están determinados de forma única.U{\displaystyle U}yV{\displaystyle V}No es necesario que sea único en general.

Descomposiciones invariantes a escala

Se refiere a variantes de descomposiciones matriciales existentes, como la SVD, que son invariantes con respecto al escalado diagonal.

  • Aplicable a: matriz A de m por n .
  • Descomposición en valores singulares invariante a la escala unitaria:A=DUSVmi{\displaystyle A=DUSV^{*}E}, donde S es una matriz diagonal no negativa única de valores singulares invariantes a escala, U y V son matrices unitarias ,V{\displaystyle V^{*}}es la transpuesta conjugada de V y de las matrices diagonales positivas D y E.
  • Comentario: Es análogo a la SVD excepto que los elementos diagonales de S son invariantes con respecto a la multiplicación izquierda y/o derecha de A por matrices diagonales no singulares arbitrarias, a diferencia de la SVD estándar para la cual los valores singulares son invariantes con respecto a la multiplicación izquierda y/o derecha de A por matrices unitarias arbitrarias.
  • Comentario: Es una alternativa a la SVD estándar cuando se requiere invariancia con respecto a transformaciones diagonales en lugar de unitarias de A.
  • Unicidad: Los valores singulares invariantes a escala deA{\displaystyle A}(dadas por los elementos diagonales de S ) siempre están determinadas de forma única. Las matrices diagonales D y E , y las matrices unitarias U y V , no son necesariamente únicas en general.
  • Comentario: Las matrices U y V no son las mismas que las obtenidas mediante la descomposición en valores singulares (SVD).

Se pueden derivar descomposiciones análogas invariantes de escala a partir de otras descomposiciones de matrices; por ejemplo, para obtener valores propios invariantes de escala. [ 3 ] [ 4 ]

descomposición de Hessenberg

  • Aplicable a: matriz cuadrada A.
  • Descomposición:A=PAGHPAG{\displaystyle A=PHP^{*}}dóndeH{\displaystyle H}es la matriz de Hessenberg yPAG{\displaystyle P}es una matriz unitaria .
  • Comentario: a menudo es el primer paso en la descomposición de Schur.

descomposición ortogonal completa

  • También conocida como: descomposición UTV , descomposición ULV , descomposición URV .
  • Aplicable a: matriz A de m por n .
  • Descomposición:A=UTV{\displaystyle A=UTV^{*}}donde T es una matriz triangular y U y V son matrices unitarias .
  • Comentario: Similar a la descomposición en valores singulares y a la descomposición de Schur.

Otras descomposiciones

descomposición polar

  • Aplicable a : cualquier matriz cuadrada compleja A.
  • Descomposición:A=UPAG{\displaystyle A=UP}(descomposición polar derecha) oA=PAGU{\displaystyle A=P'U}(descomposición polar izquierda), donde U es una matriz unitaria y P y P' son matrices hermíticas semidefinidas positivas .
  • Unicidad:PAG{\displaystyle P}siempre es único e igual aAA{\displaystyle {\sqrt {A^{*}A}}}(que siempre es hermitiana y semidefinida positiva). SiA{\displaystyle A}es invertible, entoncesU{\displaystyle U}es único.
  • Comentario: Dado que cualquier matriz hermitiana admite una descomposición espectral con una matriz unitaria,PAG{\displaystyle P}se puede escribir comoPAG=VDV{\displaystyle P=VDV^{*}}. DesdePAG{\displaystyle P}es semidefinida positiva, todos los elementos enD{\displaystyle D}son no negativos. Dado que el producto de dos matrices unitarias es unitario, tomandoW=UV{\displaystyle W=UV}uno puede escribirA=U(VDV)=WDV{\displaystyle A=U(VDV^{*})=WDV^{*}}que es la descomposición en valores singulares. Por lo tanto, la existencia de la descomposición polar es equivalente a la existencia de la descomposición en valores singulares.

Descomposición polar algebraica

  • Aplicable a: matriz cuadrada, compleja y no singular A. [ 5 ]
  • Descomposición:A=QS{\displaystyle A=QS}donde Q es una matriz ortogonal compleja y S es una matriz simétrica compleja.
  • Singularidad: SiATA{\displaystyle A^{\mathsf {T}}A}Si no tiene valores propios reales negativos, entonces la descomposición es única. [ 6 ]
  • Comentario: La existencia de esta descomposición es equivalente aAAT{\displaystyle AA^{\mathsf {T}}}ser similar aATA{\displaystyle A^{\mathsf {T}}A}. [ 7 ]
  • Comentario: Una variante de esta descomposición esA=Rdo{\displaystyle A=RC}, donde R es una matriz real y C es una matriz circular . [ 6 ]

Descomposición de Mostow

  • Aplicable a: matriz cuadrada, compleja y no singular A. [ 8 ] [ 9 ]
  • Descomposición:A=UmiiMETROmiS{\displaystyle A=Ue^{iM}e^{S}}donde U es unitaria, M es real antisimétrica y S es real simétrica.
  • Comentario: La matriz A también se puede descomponer comoA=U2miS2miiMETRO2{\displaystyle A=U_{2}e^{S_{2}}e^{iM_{2}}}, donde U 2 es unitaria, M 2 es real antisimétrica y S 2 es real simétrica. [ 6 ]

Forma normal de Sinkhorn

  • Aplicable a: matriz cuadrada real A con elementos estrictamente positivos.
  • Descomposición:A=D1SD2{\displaystyle A=D_{1}SD_{2}}donde S es doblemente estocástica y D 1 y D 2 son matrices diagonales reales con elementos estrictamente positivos.

Descomposición sectorial

  • Aplicable a: matriz cuadrada compleja A con rango numérico contenido en el sectorSα={rmiiθdor>0,|θ|α<π2}{\displaystyle S_{\alpha }=\left\{re^{i\theta }\in \mathbb {C} \mid r>0,|\theta |\leq \alpha <{\frac {\pi }{2}}\right\}}.
  • Descomposición:A=doZdo{\displaystyle A=CZC^{*}}, donde C es una matriz compleja invertible yZ=diagnóstico(miiθ1,,miiθnorte){\displaystyle Z=\operatorname {diag} \left(e^{i\theta _{1}},\ldots ,e^{i\theta _{n}}\right)}con todo|θj|α{\displaystyle \left|\theta _{j}\right|\leq \alpha }. [ 10 ] [ 11 ]

La forma normal de Williamson

  • Aplicable a: matriz real cuadrada definida positiva A de orden 2 n ×2 n .
  • Descomposición:A=STdiagnóstico(D,D)S{\displaystyle A=S^{\mathsf {T}}\operatorname {diag} (D,D)S}, dóndeSSp(2norte){\displaystyle S\in {\text{Sp}}(2n)}es una matriz simpléctica y D es una matriz diagonal no negativa de n por n . [ 12 ]

raíz cuadrada de matriz

  • Descomposición:A=BB{\displaystyle A=BB}, no es único en general.
  • En el caso de semidefinido positivoA{\displaystyle A}, hay un semidefinido positivo únicoB{\displaystyle B}de tal manera queA=BB=BB{\displaystyle A=B^{*}B=BB}.

Generalizaciones

Existen análogos de las factorizaciones SVD, QR, LU y Cholesky para cuasimatrices y cmatrices o matrices continuas . [ 13 ] Una "cuasimatrice" es, como una matriz, un esquema rectangular cuyos elementos están indexados, pero un índice discreto se reemplaza por un índice continuo. De igual modo, una "cmatrix" es continua en ambos índices. Como ejemplo de una cmatrix, se puede pensar en el núcleo de un operador integral .

Estas factorizaciones se basan en trabajos iniciales de Fredholm (1903) , Hilbert (1904) y Schmidt (1907) . Para una descripción y traducción al inglés de los trabajos fundamentales, véase Stewart (2011) .

Véase también

Referencias

Notas

  1. Sin embargo, si se utiliza una matriz no cuadrada, la matriz U tendrá la misma forma rectangular que la matriz original A. Por lo tanto, sería incorrectodenominar a la matriz U triangular superior, ya que el término correcto sería que U es la «forma escalonada por filas» de A. Aparte de esto, no existen diferencias en la factorización LU para matrices cuadradas y no cuadradas.

Citas

  1. Lay, David C. (2016). Álgebra lineal y sus aplicaciones . Steven R. Lay, Judith McDonald (Quinta  edición global). Harlow. pág.  142. ISBN 978-1-292-09223-2OCLC 920463015 {{cite book}}: CS1 mantenimiento: falta el editor de ubicación ( enlace )
  2. Piziak, R.; Odell, PL (1 de junio de 1999). "Factorización de rango completo de matrices". Mathematics Magazine . 72 (3): 193. doi : 10.2307/2690882 . JSTOR 2690882 . 
  3. Uhlmann, JK (2018), "Una inversa matricial generalizada que es consistente con respecto a las transformaciones diagonales", SIAM Journal on Matrix Analysis and Applications , 239 (2): 781–800 , doi : 10.1137/17M113890X
  4. Uhlmann, JK (2018), "Una inversa matricial generalizada que preserva el rango para la consistencia con respecto a la similitud", IEEE Control Systems Letters , 3 : 91–95 , arXiv : 1804.07334 , doi : 10.1109/LCSYS.2018.2854240 , ISSN 2475-1456 , S2CID 5031440  
  5. Choudhury y Horn 1987 , págs. 219–225 
  6. 1 2 3 Bhatia, Rajendra (2013-11-15). "La descomposición bipolar". Álgebra lineal y sus aplicaciones . 439 (10): 3031– 3037. doi : 10.1016/j.laa.2013.09.006 .
  7. Horn y Merino 1995 , págs. 43–92 
  8. Mostow, GD (1955), Algunos nuevos teoremas de descomposición para grupos semisimples , Mem. Amer. Math. Soc., vol. 14, American Mathematical Society, pp . 31–54  
  9. Nielsen, Frank; Bhatia, Rajendra (2012). Geometría de la información matricial . Springer. pág. 224. arXiv : 1007.4402 . doi : 10.1007/978-3-642-30232-9 . ISBN  978-3-642-30232-9. S2CID 118466496 . 
  10. Zhang, Fuzhen (30 de junio de 2014). "Una descomposición matricial y sus aplicaciones" . Álgebra lineal y multilineal . 63 (10): 2033–2042 . doi : 10.1080/03081087.2014.933219 . S2CID 19437967 . 
  11. Drury, SW (noviembre de 2013). "Desigualdades determinantes de Fischer y la conjetura de Higham" . Álgebra lineal y sus aplicaciones . 439 (10): 3129– 3133. doi : 10.1016/j.laa.2013.08.031 .
  12. Idel, Martin; Soto Gaona, Sebastián; Wolf, Michael M. (2017-07-15). "Límites de perturbación para la forma normal simpléctica de Williamson". Álgebra lineal y sus aplicaciones . 525 : 45–58 . arXiv : 1609.01338 . doi : 10.1016/j.laa.2017.03.013 . S2CID 119578994 . 
  13. Townsend y Trefethen 2015

Bibliografía

  • Choudhury, Dipa; Horn, Roger A. (abril de 1987). "Un análogo ortogonal-simétrico complejo de la descomposición polar". SIAM Journal on Algebraic and Discrete Methods . 8 (2): 219– 225. doi : 10.1137/0608019 .
  • Fredholm, I. (1903), "Sur une classe d''equations fonctionnelles", Acta Mathematica (en francés), 27 : 365– 390, doi : 10.1007/bf02421317
  • Hilbert, D. (1904), "Grundzüge einer allgemeinen Theorie der linearen Integralgleichungen", Nachr. Königl. Ges. Gött (en alemán), 1904 : 49– 91
  • Horn, Roger A.; Merino, Dennis I. (enero de 1995). "Equivalencia contragrediente: una forma canónica y algunas aplicaciones" . Álgebra lineal y sus aplicaciones . 214 : 43–92 . doi : 10.1016/0024-3795(93)00056-6 .
  • Meyer, CD (2000), Análisis matricial y álgebra lineal aplicada , SIAM , ISBN 978-0-89871-454-8
  • Schmidt, E. (1907), "Zur Theorie der linearen und nichtlinearen Integralgleichungen. I Teil. Entwicklung willkürlichen Funktionen nach System vorgeschriebener" , Mathematische Annalen (en alemán), 63 (4): 433– 476, doi : 10.1007/bf01449770
  • Simon, C.; Blume, L. (1994). Matemáticas para economistas . Norton. ISBN 978-0-393-95733-4.
  • Stewart, GW (2011), Fredholm, Hilbert, Schmidt: tres artículos fundamentales sobre ecuaciones integrales (PDF) , consultado el 6 de enero de 2015.
  • Townsend, A.; Trefethen, LN (2015), "Análogos continuos de factorizaciones matriciales", Proc. R. Soc. A , 471 (2173) 20140585, Bibcode : 2014RSPSA.47140585T , doi : 10.1098/rspa.2014.0585 , PMC 4277194 , PMID 25568618  
  • Jun, Lu (2021), Descomposición numérica de matrices y sus aplicaciones modernas: Un primer curso riguroso , arXiv : 2107.02579

  • Calculadora de matrices en línea archivada el 12/12/2008 en Wayback Machine.
  • Cálculo de la descomposición matricial de Wolfram Alpha  » Descomposición LU y QR
  • Enciclopedia Springer de Matemáticas  » Factorización matricial
  • GraphLab es una biblioteca de filtrado colaborativo que permite la implementación paralela a gran escala de métodos de descomposición de matrices (en C++) para sistemas multinúcleo.