Pagina 1 di 1

Soluzione di sistemi lineari con metodi iterativi

MessaggioInviato: 11/09/2019, 10:38
da alexmazz
Salve a tutti, vorrei sapere se sto ragionando bene. L'esercizio è il seguente:
dato il seguente metodo di iterazione $X^(k+1) = D^-1*(D+A)X^(k) +B$; con $D$ matrice diagonale $D_(ii) = A_(ii)$
Verificare se può convergere alla soluzione del sistema lineare $AX = B$ per ogni scelta dell'approssimazione iniziale $X^(0)$

Io procederei così: verificherei se la soluzione del sistema $X=A^-1B$ è punto unito della funzione di iterazione, sostituendo $X=A^-1B$ ad ogni $X$ della funzione. In teoria se ottengo un'identità il metodo dovrebbe convergere, altrimenti no. E' corretto?
Vi ringrazio in anticipo

Re: Soluzione di sistemi lineari con metodi iterativi

MessaggioInviato: 11/09/2019, 10:46
da feddy
Ragiona sul raggio spettrale

Re: Soluzione di sistemi lineari con metodi iterativi

MessaggioInviato: 11/09/2019, 11:59
da alexmazz
Però $A$ è una matrice generica, non ho gli effettivi valori degli elementi che la compongono. Come posso quindi ragionare sul raggio spettrale della matrice di iterazione $D^−1(D+A)$ , per verificare che sia minore di 1 ?

Re: Soluzione di sistemi lineari con metodi iterativi

MessaggioInviato: 12/09/2019, 10:52
da feddy
Hai ragione, ho risposto di fretta.
Se prendi un punto fisso di tale iterazione, questo soddisfa $x=D^{-1}(D+A)x +b$, da cui segue subito $Ax = -Db$. Dunque, se $x$ converge a qualcosa, non è soluzione del sistema $Ax=b$

Re: Soluzione di sistemi lineari con metodi iterativi

MessaggioInviato: 12/09/2019, 14:29
da alexmazz
feddy ha scritto:Hai ragione, ho risposto di fretta.
Se prendi un punto fisso di tale iterazione, questo soddisfa $x=D^{-1}(D+A)x +b$, da cui segue subito $Ax = -Db$. Dunque, se $x$ converge a qualcosa, non è soluzione del sistema $Ax=b$



Sì ok, mi trovo con quello che dici, ma aggiungerei che se la matrice $A$ ha elementi diagonali tutti pari a $-1$, avrei $D=-I$ con $I$ matrice identità. Quindi il teorema in questo caso può convergere alla soluzione del sistema $AX=b$, giusto?

Re: Soluzione di sistemi lineari con metodi iterativi

MessaggioInviato: 12/09/2019, 15:28
da feddy
Certo, solo nel caso in cui $D=-I_n$ si ha convergenza alla soluzione.

Oltretutto, nota che se $A$ avesse sulla diagonale un elemento uguale a zero, allora l'inversa di $D$ non esiste nemmeno.

P.S. E' la $x$ che converge, non il teorema :lol:

Re: Soluzione di sistemi lineari con metodi iterativi

MessaggioInviato: 12/09/2019, 17:11
da alexmazz
Si scusami, hai ragione.
Grazie mille :D

Re: Soluzione di sistemi lineari con metodi iterativi

MessaggioInviato: 12/09/2019, 17:36
da feddy
Di nulla ! :)