Articulo de referencia

Matriz de bloques

En matemáticas , una matriz de bloques o una matriz particionada es una matriz que se interpreta como dividida en secciones llamadas bloques o submatrices . [ 1 ] [ 2 ] Intuitiv...

En matemáticas , una matriz de bloques o una matriz particionada es una matriz que se interpreta como dividida en secciones llamadas bloques o submatrices . [ 1 ] [ 2 ]

Intuitivamente, una matriz interpretada como una matriz de bloques puede visualizarse como la matriz original con una colección de líneas horizontales y verticales que la dividen o la particionan en una colección de matrices más pequeñas. [ 3 ] [ 2 ] Por ejemplo, la matriz de 3×4 que se presenta a continuación está dividida por líneas horizontales y verticales en cuatro bloques: el bloque superior izquierdo de 2×3, el bloque superior derecho de 2×1, el bloque inferior izquierdo de 1×3 y el bloque inferior derecho de 1×1.

[a11a12a13b1a21a22a23b2do1do2do3d]{\displaystyle \left[{\begin{array}{ccc|c}a_{11}&a_{12}&a_{13}&b_{1}\\a_{21}&a_{22}&a_{23}&b_{2}\\\hline c_{1}&c_{2}&c_{3}&d\end{array}}\right]}

Cualquier matriz puede interpretarse como una matriz de bloques de una o más maneras, y cada interpretación se define por cómo se dividen sus filas y columnas.

Esta noción se puede precisar más para unnorte{\displaystyle n}pormetro{\displaystyle m}matrizMETRO{\displaystyle M}mediante particiónnorte{\displaystyle n}en una coleccióngrupos de filas{\displaystyle {\text{grupos de filas}}}y luego particionandometro{\displaystyle m}en una coleccióncolgroups{\displaystyle {\text{colgroups}}}. La matriz original se considera entonces como el "total" de estos grupos, en el sentido de que la(i,j){\displaystyle (i,j)}una entrada de la matriz original se corresponde de forma uno a uno con alguna(s,t){\displaystyle (s,t)}entrada compensada de algunos(incógnita,y){\displaystyle (x,y)}, dóndeincógnitagrupos de filas{\displaystyle x\in {\text{rowgroups}}}yycolgroups{\displaystyle y\in {\text{colgroups}}}. [ 4 ]

El álgebra de matrices por bloques surge en general de biproductos en categorías de matrices. [ 5 ]

Una matriz de bloques de 168×168 elementos con submatrices de 12×12, 12×24, 24×12 y 24×24. Los elementos distintos de cero están en azul, los elementos cero están en gris.

Ejemplo

La matriz

PAG=[1227156233453367]{\displaystyle \mathbf {P} ={\begin{bmatrix}1&2&2&7\\1&5&6&2\\3&3&4&5\\3&3&6&7\end{bmatrix}}}

puede visualizarse como dividido en cuatro bloques, como

PAG=[1227156233453367].{\displaystyle \mathbf {P} =\left[{\begin{array}{cc|cc}1&2&2&7\\1&5&6&2\\\hline 3&3&4&5\\3&3&6&7\end{array}}\right].}

Las líneas horizontales y verticales no tienen un significado matemático especial, [ 6 ] [ 7 ] pero son una forma común de visualizar una partición. [ 6 ] [ 7 ] Mediante esta partición,PAG{\displaystyle P}se divide en cuatro bloques de 2×2, como

PAG11=[1215],PAG12=[2762],PAG21=[3333],PAG22=[4567].{\displaystyle {\begin{aligned}\mathbf {P} _{11}&={\begin{bmatrix}1&2\\1&5\end{bmatrix}},&\mathbf {P} _{12}&={\begin{bmatrix}2&7\\6&2\end{bmatrix}},\\[1ex]\mathbf {P} _{21}&={\begin{bmatrix}3&3\\3&3\end{bmatrix}},&\mathbf {P} _{22}&={\begin{bmatrix}4&5\\6&7\end{bmatrix}}.\end{aligned}}}

La matriz particionada se puede escribir entonces como [ 8 ]

PAG=[PAG11PAG12PAG21PAG22].{\displaystyle \mathbf {P} ={\begin{bmatrix}\mathbf {P} _{11}&\mathbf {P} _{12}\\\mathbf {P} _{21}&\mathbf {P} _{22}\end{bmatrix}}.}

Definición formal

DejarAdometro×norte{\displaystyle A\in \mathbb {C} ^{m\times n}}. Una partición deA{\displaystyle A}es una representación deA{\displaystyle A}en la forma

A=[A11A12A1qA21A22A2qApag1Apag2Apagq],{\displaystyle A={\begin{bmatrix}A_{11}&A_{12}&\cdots &A_{1q}\\A_{21}&A_{22}&\cdots &A_{2q}\\\vdots &\vdots &\ddots &\vdots \\A_{p1}&A_{p2}&\cdots &A_{pq}\end{bmatrix}},}

dóndeAijdometroi×nortej{\displaystyle A_{ij}\in \mathbb {C} ^{m_{i}\times n_{j}}}son submatrices contiguas,i=1pagmetroi=metro{\textstyle \sum _{i=1}^{p}m_{i}=m}, yj=1qnortej=norte{\textstyle \sum _{j=1}^{q}n_{j}=n}. [ 9 ] Los elementosAij{\displaystyle A_{ij}}de la partición se denominan bloques . [ 9 ]

Según esta definición, los bloques de cualquier columna deben tener todos el mismo número de columnas. [ 9 ] De manera similar, los bloques de cualquier fila deben tener el mismo número de filas. [ 9 ]

Métodos de particionamiento

Una matriz puede particionarse de muchas maneras. [ 9 ] Por ejemplo, una matrizA{\displaystyle A}Se dice que está particionado por columnas si se escribe como

A=(a1 a2  anorte),{\displaystyle A=(a_{1}\ a_{2}\ \cdots \ a_{n}),}

dóndeaj{\displaystyle a_{j}}es elj{\displaystyle j}columna deA{\displaystyle A}. [ 9 ] Una matriz también puede particionarse por filas :

A=[a1Ta2TametroT],{\displaystyle A={\begin{bmatrix}a_{1}^{T}\\a_{2}^{T}\\\vdots \\a_{m}^{T}\end{bmatrix}},}

dóndeaiT{\displaystyle a_{i}^{T}}es eli{\displaystyle i}-fila deA{\displaystyle A}. [ 9 ]

particiones comunes

A menudo, [ 9 ] nos encontramos con la partición 2×2

A=[A11A12A21A22],{\displaystyle A={\begin{bmatrix}A_{11}&A_{12}\\A_{21}&A_{22}\end{bmatrix}},}

particularmente en la forma dondeA11{\displaystyle A_{11}}es un escalar: [ 9 ]

A=[a11a12Ta21A22].{\displaystyle A={\begin{bmatrix}a_{11}&a_{12}^{T}\\a_{21}&A_{22}\end{bmatrix}}.}

Operaciones de matriz de bloques

Transponer

Dejar

A=[A11A12A1qA21A22A2qApag1Apag2Apagq]{\displaystyle A={\begin{bmatrix}A_{11}&A_{12}&\cdots &A_{1q}\\A_{21}&A_{22}&\cdots &A_{2q}\\\vdots &\vdots &\ddots &\vdots \\A_{p1}&A_{p2}&\cdots &A_{pq}\end{bmatrix}}}

dóndeAijdoki×j{\displaystyle A_{ij}\in \mathbb {C} ^{k_{i}\times \ell _{j}}}. (Esta matrizA{\displaystyle A}se reutilizará en §  Suma y §  Multiplicación .) Entonces su transpuesta es [ 9 ] [ 10 ]

AT=[A11TA21TApag1TA12TA22TApag2TA1qTA2qTApagqT],{\displaystyle A^{T}={\begin{bmatrix}A_{11}^{T}&A_{21}^{T}&\cdots &A_{p1}^{T}\\A_{12}^{T}&A_{22}^{T}&\cdots &A_{p2}^{T}\\\vdots &\vdots &\ddots &\vdots \\A_{1q}^{T}&A_{2q}^{T}&\cdots &A_{pq}^{T}\end{bmatrix}},}

y la misma ecuación se mantiene reemplazando la transpuesta por la transpuesta conjugada. [ 9 ]

Transposición de bloques

También se puede definir una forma especial de transposición de matrices para matrices de bloques, donde los bloques individuales se reordenan pero no se transponen.A=(Bij){\displaystyle A=(B_{ij})}ser unk×l{\displaystyle k\times l}matriz de bloques conmetro×norte{\displaystyle m\times n}bloquesBij{\displaystyle B_{ij}}, la transpuesta en bloque deA{\displaystyle A}es ell×k{\displaystyle l\times k}matriz de bloquesAB{\displaystyle A^{\mathcal {B}}}conmetro×norte{\displaystyle m\times n}bloques(AB)ij=Bji{\displaystyle \left(A^{\mathcal {B}}\right)_{ij}=B_{ji}}. [ 11 ] Al igual que con el operador de traza convencional, la transposición de bloques es una aplicación lineal tal que(A+do)B=AB+doB{\displaystyle (A+C)^{\mathcal {B}}=A^{\mathcal {B}}+C^{\mathcal {B}}}. [ 10 ] Sin embargo, en general la propiedad(Ado)B=doBAB{\displaystyle (AC)^{\mathcal {B}}=C^{\mathcal {B}}A^{\mathcal {B}}}no se sostiene a menos que los bloques deA{\displaystyle A}ydo{\displaystyle C}desplazarse.

Suma

Dejar

B=[B11B12B1sB21B22B2sBr1Br2Brs],{\displaystyle B={\begin{bmatrix}B_{11}&B_{12}&\cdots &B_{1s}\\B_{21}&B_{22}&\cdots &B_{2s}\\\vdots &\vdots &\ddots &\vdots \\B_{r1}&B_{r2}&\cdots &B_{rs}\end{bmatrix}},}

dóndeBijdometroi×nortej{\displaystyle B_{ij}\in \mathbb {C} ^{m_{i}\times n_{j}}}y dejarA{\displaystyle A}sea ​​la matriz definida en §  Transpuesta . (Esta matrizB{\displaystyle B}se reutilizará en §  Multiplicación .) Entonces sipag=r{\displaystyle p=r},q=s{\displaystyle q=s},ki=metroi{\displaystyle k_{i}=m_{i}}, yj=nortej{\displaystyle \ell _{j}=n_{j}}, entonces [ 9 ]

A+B=[A11+B11A12+B12A1q+B1qA21+B21A22+B22A2q+B2qApag1+Bpag1Apag2+Bpag2Apagq+Bpagq].{\displaystyle A+B={\begin{bmatrix}A_{11}+B_{11}&A_{12}+B_{12}&\cdots &A_{1q}+B_{1q}\\A_{21}+B_{21}&A_{22}+B_{22}&\cdots &A_{2q}+B_{2q}\\\vdots &\vdots &\ddots &\vdots \\A_{p1}+B_{p1}&A_{p2}+B_{p2}&\cdots &A_{pq}+B_{pq}\end{bmatrix}}.}

Multiplicación

Es posible utilizar un producto matricial particionado por bloques que solo implica álgebra sobre submatrices de los factores. Sin embargo, la partición de los factores no es arbitraria y requiere " particiones conformes " [ 12 ] entre dos matrices.A{\displaystyle A}yB{\displaystyle B}de tal manera que se definan todos los productos de submatrices que se utilizarán. [ 13 ]

Dos matricesA{\displaystyle A}yB{\displaystyle B}Se dice que están particionados conformemente para el producto.AB{\displaystyle AB}, cuandoA{\displaystyle A}yB{\displaystyle B}se dividen en submatrices y si la multiplicaciónAB{\displaystyle AB}se lleva a cabo tratando las submatrices como si fueran escalares, pero manteniendo el orden, y cuando se definen todos los productos y sumas de las submatrices involucradas.

Arak M. Mathai y Hans J. Haubold, Álgebra lineal: un curso para físicos e ingenieros [ 14 ]

DejarA{\displaystyle A}Sea la matriz definida en §  Transpuesta , y seaB{\displaystyle B}Sea la matriz definida en §  Adición . Entonces el producto de la matriz

do=AB{\displaystyle C=AB}

se puede realizar por bloques, lo que producedo{\displaystyle C}como un(pag×s){\displaystyle (p\times s)}matriz. Las matrices en la matriz resultantedo{\displaystyle C}se calculan multiplicando: [ 6 ]

doij=k=1qAikBkj.{\displaystyle C_{ij}=\sum _{k=1}^{q}A_{ik}B_{kj}.}

O bien, utilizando la notación de Einstein que implícitamente suma sobre índices repetidos:

doij=AikBkj.{\displaystyle C_{ij}=A_{ik}B_{kj}.}

Representandodo{\displaystyle C}como matriz, tenemos [ 9 ]

do=AB=[i=1qA1iBi1i=1qA1iBi2i=1qA1iBisi=1qA2iBi1i=1qA2iBi2i=1qA2iBisi=1qApagiBi1i=1qApagiBi2i=1qApagiBis].{\displaystyle C=AB={\begin{bmatrix}\sum _{i=1}^{q}A_{1i}B_{i1}&\sum _{i=1}^{q}A_{1i}B_{i2}&\cdots &\sum _{i=1}^{q}A_{1i}B_{is}\\\sum _{i=1}^{q}A_{2i}B_{i1}&\sum _{i=1}^{q}A_{2i}B_{i2}&\cdots &\sum _{i=1}^{q}A_{2i}B_{is}\\\vdots &\vdots &\ddots &\vdots \\\sum _{i=1}^{q}A_{pi}B_{i1}&\sum _{i=1}^{q}A_{pi}B_{i2}&\cdots &\sum _{i=1}^{q}A_{pi}B_{is}\end{bmatrix}}.}

Inversión

Si una matriz se divide en cuatro bloques, se puede invertir bloque a bloque de la siguiente manera:

PAG=[ABdoD]1=[A1+A1B(DdoA1B)1doA1A1B(DdoA1B)1(DdoA1B)1doA1(DdoA1B)1],{\displaystyle {\begin{aligned}P&={\begin{bmatrix}A&B\\C&D\end{bmatrix}}^{-1}\\[1ex]&={\begin{bmatrix}A^{-1}+A^{-1}B\left(D-CA^{-1}B\right)^{-1}CA^{-1}&-A^{-1}B\left(D-CA^{-1}B\right)^{-1}\\-\left(D-CA^{-1}B\right)^{-1}CA^{-1}&\left(D-CA^{-1}B\right)^{-1}\end{bmatrix}},\end{aligned}}}

donde A y D son bloques cuadrados de tamaño arbitrario, y B y C son compatibles con ellos para la partición. Además, A y el complemento de Schur de A en P : P / A = DCA −1 B deben ser invertibles. [ 15 ]

De forma equivalente, permutando los bloques: [ 16 ]

PAG=[ABdoD]1=[(ABD1do)1(ABD1do)1BD1D1do(ABD1do)1D1+D1do(ABD1do)1BD1].{\displaystyle {\begin{aligned}P&={\begin{bmatrix}A&B\\C&D\end{bmatrix}}^{-1}\\[1ex]&={\begin{bmatrix}\left(A-BD^{-1}C\right)^{-1}&-\left(A-BD^{-1}C\right)^{-1}BD^{-1}\\-D^{-1}C\left(A-BD^{-1}C\right)^{-1}&D^{-1}+D^{-1}C\left(A-BD^{-1}C\right)^{-1}BD^{-1}\end{bmatrix}}.\end{aligned}}}

Aquí, D y el complemento de Schur de D en P : P / D = ABD −1 C deben ser invertibles.

Si A y D son ambos invertibles, entonces:

[ABdoD]1=[(ABD1do)100(DdoA1B)1][IBD1doA1I].{\displaystyle {\begin{bmatrix}A&B\\C&D\end{bmatrix}}^{-1}={\begin{bmatrix}\left(A-BD^{-1}C\right)^{-1}&0\\0&\left(D-CA^{-1}B\right)^{-1}\end{bmatrix}}{\begin{bmatrix}I&-BD^{-1}\\-CA^{-1}&I\end{bmatrix}}.}

Según la identidad de Weinstein-Aronszajn , una de las dos matrices en la matriz diagonal por bloques es invertible exactamente cuando la otra lo es. La inversión de la matriz por bloques también permite aprovechar la eficiencia de los algoritmos de multiplicación rápida de matrices , lo que permite realizar la inversión en tiempoO(norteω){\displaystyle O({n^{\omega }})}para 2.37ω<3{\displaystyle ~2.37\leq \omega <3}[ 17 ] , Sec. 11, págs. 413-414.

Cálculo de inversas de submatrices a partir de la inversa completa

Debido a la simetría entre una matriz y su inversa en la fórmula de inversión por bloques, si una matriz P y su inversa P −1 se particionan conformemente:

PAG=[ABdoD],PAG1=[miFGRAMOH]{\displaystyle P={\begin{bmatrix}{A}&{B}\\{C}&{D}\end{bmatrix}},\quad P^{-1}={\begin{bmatrix}{E}&{F}\\{G}&{H}\end{bmatrix}}}

Entonces, la inversa de cualquier submatriz principal se puede calcular a partir de los bloques correspondientes de P −1 :

A1=miFH1GRAMO{\displaystyle {A}^{-1}={E}-{FH}^{-1}{G}}D1=HGRAMOmi1F{\displaystyle {D}^{-1}={H}-{GE}^{-1}{F}}

Esta relación se deriva del reconocimiento de que E −1 = ABD −1 C (el complemento de Schur) y de la aplicación de la misma fórmula de inversión de bloques con los roles de P y P −1 invertidos. [ 18 ] [ 19 ]

Determinante

La fórmula para el determinante de un2×2{\displaystyle 2\times 2}La matriz anterior sigue siendo válida, bajo supuestos adicionales apropiados, para una matriz compuesta por cuatro submatrices.A,B,do,D{\displaystyle A,B,C,D}conA{\displaystyle A}yD{\displaystyle D}cuadrado. La fórmula más sencilla de este tipo, que puede demostrarse utilizando la fórmula de Leibniz o una factorización que involucre el complemento de Schur , es [ 16 ].det[A0doD]=det(A)det(D)=det[AB0D].{\displaystyle \det {\begin{bmatrix}A&0\\C&D\end{bmatrix}}=\det(A)\det(D)=\det {\begin{bmatrix}A&B\\0&D\end{bmatrix}}.}

Utilizando esta fórmula, podemos derivar que los polinomios característicos de[A0doD]{\displaystyle {\begin{bmatrix}A&0\\C&D\end{bmatrix}}}y[AB0D]{\displaystyle {\begin{bmatrix}A&B\\0&D\end{bmatrix}}}son iguales y equivalentes al producto de los polinomios característicos deA{\displaystyle A}yD{\displaystyle D}. Además, si[A0doD]{\displaystyle {\begin{bmatrix}A&0\\C&D\end{bmatrix}}}o[AB0D]{\displaystyle {\begin{bmatrix}A&B\\0&D\end{bmatrix}}}es diagonalizable , entoncesA{\displaystyle A}yD{\displaystyle D}También son diagonalizables. Lo contrario es falso; simplemente compruébelo.[1101]{\displaystyle {\begin{bmatrix}1&1\\0&1\end{bmatrix}}}.

SiA{\displaystyle A}es invertible , uno tiene [ 16 ]

det[ABdoD]=det(A)det(DdoA1B),{\displaystyle \det {\begin{bmatrix}A&B\\C&D\end{bmatrix}}=\det(A)\det \left(D-CA^{-1}B\right),}

y siD{\displaystyle D}es invertible, uno tiene [ 20 ] [ 16 ]

det[ABdoD]=det(D)det(ABD1do).{\displaystyle \det {\begin{bmatrix}A&B\\C&D\end{bmatrix}}=\det(D)\det \left(A-BD^{-1}C\right).}

Si los bloques son matrices cuadradas del mismo tamaño, se cumplen otras fórmulas. Por ejemplo, sido{\displaystyle C}yD{\displaystyle D}desplazamiento (es decir,doD=Ddo{\displaystyle CD=DC}), entonces [ 21 ]det[ABdoD]=det(ADBdo).{\displaystyle \det {\begin{bmatrix}A&B\\C&D\end{bmatrix}}=\det(AD-BC).} Afirmaciones similares son válidas cuandoAB=BA{\displaystyle AB=BA},Ado=doA{\displaystyle AC=CA}, oBD=DB{\displaystyle BD=DB} . Es decir, siAdo=doA{\displaystyle AC=CA}, entonces det[ABdoD]=det(ADdoB).{\displaystyle \det {\begin{bmatrix}A&B\\C&D\end{bmatrix}}=\det(AD-CB).} Observe el cambio en el orden dedo{\displaystyle C}yB{\displaystyle B}(tenemosdoB{\displaystyle CB}en lugar deBdo{\displaystyle BC}). De manera similar, siBD=DB{\displaystyle BD=DB}, entoncesAD{\displaystyle AD}debe ser reemplazado porDA{\displaystyle DA}(es decir, obtenemosdet(DABdo){\displaystyle \det(DA-BC)}) y siAB=BA{\displaystyle AB=BA}, entonces deberíamos tenerdet(DAdoB){\displaystyle \det(DA-CB)}. Nota: para los dos últimos resultados, hay que usar la conmutatividad del anillo subyacente, pero no para los dos primeros.

Esta fórmula se ha generalizado a matrices compuestas por más de2×2{\displaystyle 2\times 2}bloques, nuevamente bajo condiciones de conmutatividad apropiadas entre los bloques individuales. [ 22 ]

ParaA=D{\displaystyle A=D}yB=do{\displaystyle B=C}, la siguiente fórmula se cumple (incluso siA{\displaystyle A}yB{\displaystyle B}no se desplazan diariamente) [ 16 ]det[ABBA]=det(AB)det(A+B).{\displaystyle \det {\begin{bmatrix}A&B\\B&A\end{bmatrix}}=\det(A-B)\det(A+B).}

Tipos especiales de matrices de bloques

Sumas directas y matrices diagonales por bloques

Suma directa

Para cualesquiera matrices arbitrarias A (de tamaño m  × n ) y B (de tamaño p × q ), tenemos la suma directa de A y B , denotada por A B y definida como [ 10 ].     

AB=[a11a1norte00ametro1ametronorte0000b11b1q00bpag1bpagq].{\displaystyle {A}\oplus {B}={\begin{bmatrix}a_{11}&\cdots &a_{1n}&0&\cdots &0\\\vdots &\ddots &\vdots &\vdots &\ddots &\vdots \\a_{m1}&\cdots &a_{mn}&0&\cdots &0\\0&\cdots &0&b_{11}&\cdots &b_{1q}\\\vdots &\ddots &\vdots &\vdots &\ddots &\vdots \\0&\cdots &0&b_{p1}&\cdots &b_{pq}\end{bmatrix}}.}

Por ejemplo,

[132231][1601]=[13200231000001600001].{\displaystyle {\begin{bmatrix}1&3&2\\2&3&1\end{bmatrix}}\oplus {\begin{bmatrix}1&6\\0&1\end{bmatrix}}={\begin{bmatrix}1&3&2&0&0\\2&3&1&0&0\\0&0&0&1&6\\0&0&0&0&1\end{bmatrix}}.}

Esta operación se generaliza de forma natural a matrices de dimensiones arbitrarias (siempre que A y B tengan el mismo número de dimensiones).

Tenga en cuenta que cualquier elemento de la suma directa de dos espacios vectoriales de matrices puede representarse como una suma directa de dos matrices.

Matrices diagonales por bloques

Una matriz diagonal por bloques es una matriz por bloques que es una matriz cuadrada tal que los bloques de la diagonal principal son matrices cuadradas y todos los bloques fuera de la diagonal son matrices nulas. [ 16 ] Es decir, una matriz diagonal por bloques A tiene la forma

A=[A1000A2000Anorte]{\displaystyle {A}={\begin{bmatrix}A_{1}&0&\cdots &0\\0&A_{2}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &A_{n}\end{bmatrix}}}

donde A k es una matriz cuadrada para todo k = 1, ..., n . En otras palabras, la matriz A es la suma directa de A 1 , ..., A n . [ 16 ] También se puede indicar como A 1 A 2 ⊕ ... ⊕ A n [ 10 ] o diag( A 1 , A 2 , ..., A n ) [ 10 ] (este último es el mismo formalismo utilizado para una matriz diagonal ). Cualquier matriz cuadrada puede considerarse trivialmente una matriz diagonal por bloques con un solo bloque.      

Para el determinante y la traza , se cumplen las siguientes propiedades:

detA=detA1××detAnorte,{\displaystyle {\begin{aligned}\det {A}&=\det {A}_{1}\times \cdots \times \det {A}_{n},\end{aligned}}}[ 23 ] [ 24 ] y
trA=trA1++trAnorte.{\displaystyle {\begin{aligned}\operatorname {tr} {A}&=\operatorname {tr} {A}_{1}+\cdots +\operatorname {tr} {A}_{n}.\end{aligned}}}[ 16 ] [ 24 ]

Una matriz diagonal por bloques es invertible si y solo si cada uno de sus bloques diagonales principales es invertible, y en este caso su inversa es otra matriz diagonal por bloques dada por [ 25 ].[A1000A2000Anorte]1=[A11000A21000Anorte1].{\displaystyle {\begin{bmatrix}{A}_{1}&{0}&\cdots &{0}\\{0}&{A}_{2}&\cdots &{0}\\\vdots &\vdots &\ddots &\vdots \\{0}&{0}&\cdots &{A}_{n}\end{bmatrix}}^{-1}={\begin{bmatrix}{A}_{1}^{-1}&{0}&\cdots &{0}\\{0}&{A}_{2}^{-1}&\cdots &{0}\\\vdots &\vdots &\ddots &\vdots \\{0}&{0}&\cdots &{A}_{n}^{-1}\end{bmatrix}}.}

Los valores propios [ 26 ] y los vectores propios deA{\displaystyle {A}}son simplemente los de laAk{\displaystyle {A}_{k}}s combinados. [ 24 ]

Matrices tridiagonales por bloques

Una matriz tridiagonal por bloques es otra matriz por bloques especial, que al igual que la matriz diagonal por bloques es una matriz cuadrada , con matrices cuadradas (bloques) en la diagonal inferior, la diagonal principal y la diagonal superior, siendo todos los demás bloques matrices nulas. Esencialmente es una matriz tridiagonal , pero tiene submatrices en lugar de escalares.A{\displaystyle A}tiene la forma

A=[B1do10A2B2do2AkBkdokAnorte1Bnorte1donorte10AnorteBnorte]{\displaystyle {A}={\begin{bmatrix}B_{1}&C_{1}&&&\cdots &&0\\A_{2}&B_{2}&C_{2}&&&&\\&\ddots &\ddots &\ddots &&&\vdots \\&&A_{k}&B_{k}&C_{k}&&\\\vdots &&&\ddots &\ddots &\ddots &\\&&&&A_{n-1}&B_{n-1}&C_{n-1}\\0&&\cdots &&&A_{n}&B_{n}\end{bmatrix}}}

dóndeAk{\displaystyle {A}_{k}},Bk{\displaystyle {B}_{k}}ydok{\displaystyle {C}_{k}}son submatrices cuadradas de la diagonal inferior, principal y superior respectivamente. [ 27 ] [ 28 ]

Las matrices tridiagonales por bloques se encuentran frecuentemente en soluciones numéricas de problemas de ingeniería (por ejemplo, dinámica de fluidos computacional ). Existen métodos numéricos optimizados para la factorización LU [ 29 ] y, por lo tanto, algoritmos de solución eficientes para sistemas de ecuaciones con una matriz tridiagonal por bloques como matriz de coeficientes. El algoritmo de Thomas , utilizado para la solución eficiente de sistemas de ecuaciones que involucran una matriz tridiagonal, también puede aplicarse mediante operaciones matriciales a matrices tridiagonales por bloques (véase también Descomposición LU por bloques ).

Matrices triangulares de bloques

Unnorte×norte{\displaystyle n\times n}matrizA{\displaystyle A}es triangular superior de bloque (o triangular superior de bloque [ 30 ] ) si hay enteros positivosnorte1,,nortek{\displaystyle n_{1},\ldots ,n_{k}}de tal manera quenorte=norte1+norte2++nortek{\displaystyle n=n_{1}+n_{2}+\ldots +n_{k}}y A=[A11A12A1k0A22A2k00Akk],{\displaystyle A={\begin{bmatrix}A_{11}&A_{12}&\cdots &A_{1k}\\0&A_{22}&\cdots &A_{2k}\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &A_{kk}\end{bmatrix}},} donde la matrizAij{\displaystyle A_{ij}}esnortei×nortej{\displaystyle n_{i}\times n_{j}}a pesar dei,j=1,,k{\displaystyle i,j=1,\ldots ,k}. [ 26 ] [ 30 ] De manera similar,A{\displaystyle A}es el bloque inferior triangular si A=[A1100A21A220Ak1Ak2Akk],{\displaystyle A={\begin{bmatrix}A_{11}&0&\cdots &0\\A_{21}&A_{22}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\A_{k1}&A_{k2}&\cdots &A_{kk}\end{bmatrix}},} dóndeAij{\displaystyle A_{ij}}esnortei×nortej{\displaystyle n_{i}\times n_{j}}a pesar dei,j=1,,k{\displaystyle i,j=1,\ldots ,k}. [ 26 ]

Matrices de Toeplitz de bloques

Una matriz de Toeplitz por bloques es otro tipo de matriz por bloques especial, que contiene bloques que se repiten a lo largo de las diagonales de la matriz, al igual que una matriz de Toeplitz tiene elementos repetidos a lo largo de la diagonal.

Una matrizA{\displaystyle A}es bloquear Toeplitz siA(i,j)=A(k,l){\displaystyle A_{(i,j)}=A_{(k,l)}}a pesar deki=lj{\displaystyle k-i=l-j}, eso es,

A=[A1A2A3A4A1A2A5A4A1],{\displaystyle A={\begin{bmatrix}A_{1}&A_{2}&A_{3}&\cdots \\A_{4}&A_{1}&A_{2}&\cdots \\A_{5}&A_{4}&A_{1}&\cdots \\\vdots &\vdots &\vdots &\ddots \end{bmatrix}},}

dóndeAiFnortei×metroi{\displaystyle A_{i}\in \mathbb {F} ^{n_{i}\times m_{i}}}. [ 26 ]

Matrices de Hankel por bloques

Una matrizA{\displaystyle A}es bloquear Hankel siA(i,j)=A(k,l){\displaystyle A_{(i,j)}=A_{(k,l)}}a pesar dei+j=k+l{\displaystyle i+j=k+l}, eso es,

A=[A1A2A3A2A3A4A3A4A5],{\displaystyle A={\begin{bmatrix}A_{1}&A_{2}&A_{3}&\cdots \\A_{2}&A_{3}&A_{4}&\cdots \\A_{3}&A_{4}&A_{5}&\cdots \\\vdots &\vdots &\vdots &\ddots \end{bmatrix}},}

dóndeAiFnortei×metroi{\displaystyle A_{i}\in \mathbb {F} ^{n_{i}\times m_{i}}}. [ 26 ]

Véase también

  • Producto de Kronecker (producto directo de matrices que da como resultado una matriz de bloques)
  • Forma normal de Jordan (forma canónica de un operador lineal en un espacio vectorial complejo de dimensión finita)
  • Algoritmo de Strassen (algoritmo para la multiplicación de matrices más rápido que el algoritmo convencional de multiplicación de matrices)

Notas

  1. Eves, Howard (1980). Teoría elemental de matrices (  edición reimpresa). Nueva York: Dover. pág . 37. ISBN  0-486-63946-0. Recuperado el 24 de abril de 2013 . Veremos que a veces es conveniente subdividir una matriz en bloques rectangulares de elementos. Esto nos lleva a considerar las llamadas matrices particionadas o de bloques .
  2. 1 2 Dobrushkin, Vladimir. "Matrices de partición" . Álgebra lineal con Mathematica . Consultado el 24 de marzo de 2024 .
  3. Anton, Howard (1994). Álgebra lineal elemental (7.ª ed.). Nueva York: John Wiley. pág. 30. ISBN   0-471-58742-7Una matriz puede subdividirse o particionarse en matrices más pequeñas insertando líneas horizontales y verticales entre filas y columnas seleccionadas .
  4. Indhumathi, D.; Sarala, S. (2014-05-16). "Análisis de fragmentos y generación de casos de prueba utilizando la medida F para pruebas aleatorias adaptativas y pruebas aleatorias adaptativas basadas en bloques particionados" (PDF) . Revista Internacional de Aplicaciones Informáticas . 93 (6): 13. Bibcode : 2014IJCA...93f..11I . doi : 10.5120/16218-5662 .
  5. Macedo, HD; Oliveira, JN (2013). "Escriba álgebra lineal: un enfoque orientado al biproducto". Science of Computer Programming . 78 (11): 2160– 2191. arXiv : 1312.4818 . doi : 10.1016/j.scico.2012.07.012 .
  6. 1 2 3 Johnston, Nathaniel (2021). Introducción al álgebra lineal y matricial . Cham, Suiza: Springer Nature. pp. 30, 425. ISBN  978-3-030-52811-9.
  7. 1 2 Johnston, Nathaniel (2021). Álgebra lineal y matricial avanzada . Cham, Suiza: Springer Nature. pág. 298. ISBN  978-3-030-52814-0.
  8. Jeffrey, Alan (2010). Operaciones matriciales para ingenieros y científicos: una guía esencial en álgebra lineal . Dordrecht [Países Bajos]; Nueva York: Springer. pág. 54. ISBN  978-90-481-9273-1OCLC 639165077 
  9. 1 2 3 4 5 6 7 8 9 10 11 12 13 Stewart, Gilbert W. (1998). Algoritmos matriciales. 1: Descomposiciones básicas . Filadelfia, PA: Soc. de Matemáticas Industriales y Aplicadas. págs. 18–20 . ISBN  978-0-89871-414-2.
  10. 1 2 3 4 5 Gentle, James E. (2007). Álgebra matricial: teoría, cálculos y aplicaciones en estadística . Textos de Springer en estadística. Nueva York, NY: Springer New York Springer e-books. págs. 47, 487. ISBN  978-0-387-70873-7.
  11. Mackey, D. Steven (2006). Linealizaciones estructuradas para polinomios matriciales (PDF) (Tesis). Universidad de Manchester. ISSN 1749-9097 . OCLC 930686781 .  
  12. Eves, Howard (1980). Teoría elemental de matrices ( edición reimpresa). Nueva York: Dover. pág . 37. ISBN   0-486-63946-0. Recuperado el 24 de abril de 2013 . Una partición como en el Teorema 1.9.4 se llama partición conforme de A y B .
  13. Anton, Howard (1994). Álgebra lineal elemental (7.ª ed.). Nueva York: John Wiley. pág. 36. ISBN   0-471-58742-7... siempre que los tamaños de las submatrices de A y B sean tales que se puedan realizar las operaciones indicadas.
  14. Mathai, Arakaparampil M.; Haubold, Hans J. (2017). Álgebra lineal: un curso para físicos e ingenieros . Libro de texto de De Gruyter. Berlín Boston: De Gruyter. pág. 162. ISBN  978-3-11-056259-0.
  15. Bernstein, Dennis (2005). Matemáticas matriciales . Princeton University Press. pág. 44. ISBN  0-691-11802-7.
  16. 1 2 3 4 5 6 7 8 Abadir, Karim M.; Magnus, enero R. (2005). Álgebra matricial . Prensa de la Universidad de Cambridge. págs.97 , 100, 106, 111, 114, 118. ISBN  9781139443647.
  17. Pan, V. (1984). "¿Cómo podemos acelerar la multiplicación de matrices?". SIAM Review . 26 (3): 393– 415. doi : 10.1137/1026076 .
  18. "¿Se conoce esta fórmula para la inversa de un bloque de matriz en términos de la inversa de la matriz completa?" . MathOverflow .
  19. Escalante-B., Alberto N.; Wiskott, Laurenz (2016). "SFA mejorado basado en grafos: la preservación de la información complementa el principio de lentitud" . Machine Learning . arXiv : 1412.4679 . doi : 10.1007/s10994-016-5563-y .
  20. Taboga, Marco (2021). "Determinante de una matriz de bloques", Lecciones de álgebra matricial.
  21. Silvester, JR (2000). "Determinantes de matrices de bloques" (PDF) . Math. Gaz . 84 (501): 460– 467. doi : 10.2307/3620776 . JSTOR 3620776. Archivado del original (PDF) el 18 de marzo de 2015. Recuperado el 25 de junio de 2021 . 
  22. Sothanaphan, Nat (enero de 2017). "Determinantes de matrices de bloques con bloques no conmutativos". Álgebra lineal y sus aplicaciones . 512 : 202–218 . arXiv : 1805.06027 . doi : 10.1016/j.laa.2016.10.004 . S2CID 119272194 . 
  23. Quarteroni, Alfio; Sacco, Riccardo; Saleri, Fausto (2000). Matemáticas numéricas . Textos de matemáticas aplicadas. Nueva York: Springer. pp. 10, 13. ISBN  978-0-387-98959-4.
  24. ^ George, Raju K .; Ajayakumar, Abhijith (2024). "Un curso de álgebra lineal" . Textos Universitarios en Ciencias Matemáticas : 35, 407. doi : 10.1007/978-981-99-8680-4 . ISBN 978-981-99-8679-8ISSN 2731-9318 
  25. Prince, Simon JD (2012). Visión por computadora: modelos, aprendizaje e inferencia . Nueva York: Cambridge University Press. pág. 531. ISBN  978-1-107-01179-3.
  26. 1 2 3 4 5 Bernstein, Dennis S. (2009). Matemáticas matriciales: teoría, hechos y fórmulas (2.ª ed.). Princeton, NJ: Princeton University Press. pp. 168, 298. ISBN   978-0-691-14039-1.
  27. Dietl, Guido KE (2007). Estimación y detección lineal en subespacios de Krylov . Fundamentos en procesamiento de señales, comunicaciones y redes. Berlín; Nueva York: Springer. pp. 85, 87. ISBN  978-3-540-68478-7OCLC 85898525 
  28. Horn, Roger A.; Johnson, Charles R. (2017). Análisis matricial (Segunda edición, reimpresión corregida ). Nueva York, NY: Cambridge University Press. pág. 36. ISBN   978-0-521-83940-2.
  29. Datta, Biswa Nath (2010). Álgebra lineal numérica y aplicaciones (2.ª ed.). Filadelfia, Pa: SIAM. pág. 168. ISBN   978-0-89871-685-6.
  30. 1 2 Stewart, Gilbert W. (2001). Algoritmos matriciales. 2: Sistemas propios . Filadelfia, Pa: Soc. de Matemáticas Industriales y Aplicadas. pág. 5. ISBN  978-0-89871-503-3.

Referencias

  • Strang, Gilbert (1999). "Clase 3: Multiplicación y matrices inversas" . Material didáctico del MIT Open Course. 18:30–21:10.