Soluzione di sistemi lineari con metodi iterativi

Messaggioda alexmazz » 11/09/2019, 09:38

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
alexmazz
Starting Member
Starting Member
 
Messaggio: 18 di 44
Iscritto il: 05/07/2017, 21:16

Re: Soluzione di sistemi lineari con metodi iterativi

Messaggioda feddy » 11/09/2019, 09:46

Ragiona sul raggio spettrale
Avatar utente
feddy
Moderatore
Moderatore
 
Messaggio: 2571 di 5934
Iscritto il: 26/06/2016, 00:25
Località: SISSA

Re: Soluzione di sistemi lineari con metodi iterativi

Messaggioda alexmazz » 11/09/2019, 10:59

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 ?
alexmazz
Starting Member
Starting Member
 
Messaggio: 19 di 44
Iscritto il: 05/07/2017, 21:16

Re: Soluzione di sistemi lineari con metodi iterativi

Messaggioda feddy » 12/09/2019, 09:52

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$
Avatar utente
feddy
Moderatore
Moderatore
 
Messaggio: 2572 di 5934
Iscritto il: 26/06/2016, 00:25
Località: SISSA

Re: Soluzione di sistemi lineari con metodi iterativi

Messaggioda alexmazz » 12/09/2019, 13:29

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?
alexmazz
Starting Member
Starting Member
 
Messaggio: 20 di 44
Iscritto il: 05/07/2017, 21:16

Re: Soluzione di sistemi lineari con metodi iterativi

Messaggioda feddy » 12/09/2019, 14:28

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:
Avatar utente
feddy
Moderatore
Moderatore
 
Messaggio: 2574 di 5934
Iscritto il: 26/06/2016, 00:25
Località: SISSA

Re: Soluzione di sistemi lineari con metodi iterativi

Messaggioda alexmazz » 12/09/2019, 16:11

Si scusami, hai ragione.
Grazie mille :D
alexmazz
Starting Member
Starting Member
 
Messaggio: 21 di 44
Iscritto il: 05/07/2017, 21:16

Re: Soluzione di sistemi lineari con metodi iterativi

Messaggioda feddy » 12/09/2019, 16:36

Di nulla ! :)
Avatar utente
feddy
Moderatore
Moderatore
 
Messaggio: 2575 di 5934
Iscritto il: 26/06/2016, 00:25
Località: SISSA


Torna a Analisi Numerica e Ricerca Operativa

Chi c’è in linea

Visitano il forum: Nessuno e 1 ospite