219 votos

Inverso de la suma de matrices

Tengo dos matrices cuadradas: A y B. Se conoce A1 y quiero calcular (A+B)1. ¿Existen teoremas que ayuden a calcular la inversa de la suma de matrices? En general, B1 no se conoce, pero si es necesario, se puede asumir que B1 también se conoce.

4 votos

Este artículo puede ser de ayuda, aunque en este momento no tengo acceso a él.

0 votos

@Arturo: Sé que pueden no ser invertibles, pero vamos a suponer que lo son. @Adrian: Desafortunadamente, no tengo acceso directo a jstor.

0 votos

192voto

Lorin Hochstein Puntos 11816

En general, A+B no siempre tiene que ser invertible, incluso cuando A y B lo son. Pero uno podría preguntarse si se puede tener una fórmula bajo la suposición adicional de que A+B es invertible.

Como señaló Adrián Barquero, existe un artículo de Ken Miller publicado en la revista Mathematics Magazine en 1981 que aborda esto.

Él demuestra lo siguiente:

Lema. Si A y A+B son invertibles, y B tiene rango 1, entonces sea g=trace(BA1). Entonces g1 y (A+B)1=A111+gA1BA1.

A partir de este lema, podemos tomar un A+B general que sea invertible y escribirlo como A+B=A+B1+B2++Br, donde cada Bi tiene rango 1 y tal que cada A+B1++Bk es invertible (dicha descomposición siempre existe si A+B es invertible y rank(B)=r). Entonces se obtiene:

Teorema. Sean A y A+B matrices no singulares, y sea B de rango r>0. Sea B=B1++Br, donde cada Bi tiene rango 1, y cada Ck+1=A+B1++Bk es no singular. Tomando C1=A, entonces Ck+11=Ck1gkCk1BkCk1 donde gk=11+trace(Ck1Bk). En particular, (A+B)1=Cr1grCr1BrCr1.

(Si el rango de B es 0, entonces B=0, por lo que (A+B)1=A1).

2 votos

Gracias, estaba buscando algo así.

30 votos

El lema es la fórmula de Sherman-Morrison, ¿verdad?

2 votos

¿Se puede utilizar este teorema para encontrar el inverso de [gμν+χkμkνk2] donde g es el tensor métrico de Minkowski y los k son cuatrivectores? Por favor, consulte esta pregunta en Physics.SE: physics.stackexchange.com/q/141613/31965 Gracias.

72voto

osama Puntos 16

Se muestra en Sobre la derivación del inverso de la suma de matrices que

(A+B)1=A1A1B(A+B)1.

Esta ecuación no puede ser utilizada para calcular (A+B)1, pero es útil para el análisis de perturbaciones donde B es una perturbación de A. Hay varias otras variaciones de la forma anterior (ver ecuaciones (22)-(26) en este documento).

Este resultado es bueno porque solo requiere que A y A+B sean no singulares. En comparación, la identidad SMW o el artículo de Ken Miller (como se menciona en las otras respuestas) requiere algunas condiciones de no singularidad o rango de B.

14 votos

¿Qué pasa con el caso de (A+λI)1? Supongamos que A es PSD.

5 votos

Esto sigue directamente de la Identidad de la Matriz de Woodbury. Deja C=I, V=I. en.m.wikipedia.org/wiki/Woodbury_matrix_identity

1 votos

46voto

Chris Puntos 133

(A+B)1=A1A1BA1+A1BA1BA1A1BA1BA1BA1+

si A1B<1 o BA1<1 (aquí significa norma). Esta es simplemente la expansión de Taylor de la función de inversión junto con información básica sobre convergencia.

(publicado prácticamente al mismo tiempo que mjqxxx)

46voto

leeand00 Puntos 6410

Esto lo encontré accidentalmente.

Supongamos que se dan A y B, donde A y A+B son invertibles. Ahora queremos conocer la expresión de (A+B)1 sin imponer el inverso completo. Ahora seguimos la intuición de esta manera. Supongamos que podemos expresar (A+B)1=A1+X, a continuación presentaremos un método simple y directo para calcular X (A+B)1=A1+X (A1+X)(A+B)=I A1A+XA+A1B+XB=I X(A+B)=A1B \begin{equation} X = - A^{-1} B ( A + B)^{-1} \end{equation> \begin{equation> X = - A^{-1} B (A^{-1} + X) \begin{equation} (I + A^{-1}B) X = - A^{-1} B A^{-1} \end{equation> \begin{equation> X = - (I + A^{-1}B)^{-1} A^{-1} B A^{-1}

Este lema es una simplificación del lema presentado por Ken Miller, 1981

5 votos

¿Dónde encontraste esto? ¿Puedes proporcionar una cita?

4 votos

¿Cómo es esto una simplificación del lema mostrado en Ken Miller 1981? ¿Estamos hablando de "Sobre el Inverso de la Suma de Matrices" u otro trabajo? (En cualquier caso, encuentro esta propiedad bastante útil, solo necesito citarla correctamente).

4 votos

Interesante notar que la línea 3 es una ecuación de Sylvester.

41voto

Wangyan Li Puntos 132

Me sorprende que nadie se dé cuenta de que es un caso especial de la bien conocida lema de inversa de matriz o [identidad de matriz de Woodbury], dice,

(A+UCV)1=A1A1U(C1+VA1U)1VA1 ,

solo configura U=V=I, inmediatamente obtiene

(A+C)1=A1A1(C1+A1)1A1 .

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X