Proyectando el espacio nulo de

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,$$ $A\in\mathbb{R}^{n\times n}$$b$$A$$b$$A$

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 ] T ∈ R n , y b tiene un componente cero en el espacio nulo de L , J b = b , donde J = I -

$$Lx=b,$$ $L$$1_n=[1\dots 1]^T\in\mathbb{R}^n$$b$$L$ $$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=I-\frac{1}{n}1_n1_n^T$$L$$L$

$A$$A$$A^T$$A^Tu=0$$Ax=b$$u^Tb=u^TAx=0$$b$$A^T$

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

$A$$A^T$$A$