Proyectando el espacio nulo de

11

Dado el sistema donde A R n × n , leí que, en caso de que la iteración de Jacobi se use como solucionador, el método no convergerá si b tiene un componente distinto de cero en el espacio nulo de A . Entonces, ¿cómo podría uno declarar formalmente eso, siempre que b tenga un componente distinto de cero que abarque el espacio nulo de A , el método de Jacobi no es convergente? Me pregunto cómo podría formalizarse matemáticamente, ya que parte de la solución ortogonal al espacio nulo converge.

Ax=b,
ARn×nbAbA

Por lo tanto, al proyectar el espacio nulo de fuera de cada iteración, converge (o?).A

.........

Estoy particularmente interesado en el caso de donde L es una matriz laplaciana simétrica con el espacio nulo atravesado por un vector 1 n = [ 1 ... 1 ] TR n , y b tiene un componente cero en el espacio nulo de L , J b = b , donde J = I - 1

Lx=b,
L1n=[11]TRnbL
Jb=b,
es la matriz de centrado. ¿Eso implica que cada iteración de Jacobi tendrá el espacio nulo deLproyectado, es decir, cada iteración estarácentrada? Pregunto esto desde entonces, no habría necesidad de proyectar el espacio nulo deL apartir de iteraciones de Jacobi (o, en otras palabras,centrarlas iteraciones).J=I1n1n1nTLL
usero
fuente
Esta pregunta también puede ser relevante para usted: scicomp.stackexchange.com/questions/1505/…
shuhalo
Gracias. Realmente hice un extracto de mis comentarios allí, ya que la pregunta merece atención por sí misma. Sin embargo, lo anterior no fue abordado (al menos no formalizado).
usero
Oh, qué vergüenza, no comprobé que era tu propia pregunta.
shuhalo
@JedBrown Su respuesta en scicomp.stackexchange.com/questions/1505/… inspiró esta pregunta. Creo que merece una consideración independiente. Supongo que podrás considerar las preguntas anteriores.
usero

Respuestas:

7

AAATATu=0Ax=buTb=uTAx=0bAT

Pero si este es el caso, existe una solución, y en el caso cuadrado hay infinitos.

AATA

Arnold Neumaier
fuente
bAAAA
AATA
Ab
1
AA=IBx0=bxn+1=b+BxnAu=0uTb=0uTB=uTuTxnes constante por inducción, por lo tanto cero. - Pero, ¿por qué te importa el método Jacobi? Es muy lento!
Arnold Neumaier
BAdiag(A)cIcR