Proiettando lo spazio nullo di


11

Dato il sistema dove A R n × n , leggo che, nel caso in cui l'iterazione Jacobi sia utilizzata come solutore, il metodo non converge se b ha un componente diverso da zero nello spazio nullo di A . Quindi, come si può affermare formalmente che, a condizione che b abbia un componente diverso da zero che copra lo spazio nullo di A , il metodo Jacobi non è convergente? Mi chiedo come possa essere formalizzato matematicamente, poiché parte della soluzione ortogonale allo spazio nullo converge.

Ax=b,
ARn×nbAbA

Pertanto, proiettando lo spazio nullo di su ciascuna iterazione, converge (o?).A

.........

Sono particolarmente interessato al caso di cui L è una matrice laplaciana simmetrica con lo spazio nullo attraversato da un vettore 1 n = [ 1 1 ] TR n , e b ha una componente zero in lo spazio nullo di L , J b = b , dove J = I - 1

Lx=b,
L1n=[11]TRnbL
Jb=b,
è la matrice di centraggio. Ciò implica che ogni iterato di Jacobi avrà lo spazio nullo diLproiettato fuori, cioè ogni iterato saràcentrato? Lo sto chiedendo da allora non sarebbe necessario proiettare lo spazio nullo diLdalle iterate di Jacobi (o, in altre parole,centrarele iterate).J=I1n1n1nTLL

Questa domanda potrebbe essere rilevante anche per te: scicomp.stackexchange.com/questions/1505/…
shuhalo,

Grazie. In realtà ho fatto un estratto dei miei commenti lì, poiché la domanda merita attenzione da sola. Tuttavia, quanto sopra non è stato affrontato (almeno non formalizzato).
usero,

Oh, peccato per me, non ho verificato che fosse una tua domanda.
shuhalo,

@JedBrown La tua risposta su scicomp.stackexchange.com/questions/1505/… ha ispirato questa domanda. Penso che meriti una considerazione indipendente. Immagino che sarai in grado di considerare le domande precedenti.
usero,

Risposte:


7

AAATATu=0Ax=buTb=uTAx=0bAT

Ma se questo è il caso, esiste una soluzione e nel caso quadrato ce ne sono infiniti.

AATA


bAAAA

AATA

Ab

1
AA=IBx0=bxn+1=b+BxnAu=0uTb=0uTB=uTuTxnè costante per induzione, quindi zero. - Ma perché ti importa del metodo Jacobi? È molto lento!
Arnold Neumaier,

BAdiag(A)cIcR
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.