Metodo di Gauss-Seidel

Da Wikipedia, l'enciclopedia libera.
Jump to navigation Jump to search

In analisi numerica il metodo di Gauss-Seidel è un metodo iterativo, simile al metodo di Jacobi, per la risoluzione di un sistema lineare, scritto nella forma matriciale

Algoritmo[modifica | modifica wikitesto]

Come per il metodo di Jacobi, la matrice viene scritta come differenza di due matrici: Nel metodo di Gauss-Seidel si prendono matrice triangolare inferiore, e matrice triangolare superiore con diagonale nulla.

Preso un qualunque vettore si costruisce quindi una successione di vettori per iterazione (come per il metodo di Jacobi)

Come sistema lineare[modifica | modifica wikitesto]

Come per il metodo di Jacobi, la ricorrenza si può ottenere riscrivendo il sistema di equazioni lineari in modo da isolare una variabile per ogni riga. In questo caso, tuttavia, ogni nuova componente di viene immediatamente sostituita per il calcolo delle successive

Ad esempio, se con il metodo di Jacobi in un sistema 3x3 si utilizza la ricorrenza

con il metodo di Gauss-Seidel si utilizza la ricorrenza

Convergenza[modifica | modifica wikitesto]

Come per il metodo di Jacobi, la successione converge indipendentemente dalla scelta di se e solo se la matrice ha tutti autovalori con valore assoluto minore di 1.

Una condizione sufficiente perché questo avvenga è che sia una matrice a diagonale dominante per righe in senso stretto. (Questa condizione implica la non singolarità di quindi l'unicità della soluzione.)

Per risolvere il sistema lineare si utilizza quindi una successione che converge verso la soluzione del sistema lineare. Siano e , la successione è costruita come segue:

Aspetti computazionali[modifica | modifica wikitesto]

A differenza del metodo di Jacobi, il metodo di Gauss-Seidel richiede di tenere in memoria un solo vettore, purché i calcoli non vengano svolti in parallelo.

Considerando il fatto che per calcolare la componente -esima del nuovo vettore della successione si utilizzano le componenti già calcolate, si osserva che una iterazione del metodo di Gauss-Seidel equivale a due iterazioni del metodo di Jacobi (che utilizza soltanto le componenti del vecchio vettore ), quindi nel metodo di Gauss-Seidel si eseguano meno iterazioni e, quindi, converge più velocemente. Tuttavia la complessità computazionale non varia ed entrambi i due metodi hanno una complessità di iterazioni.

Voci correlate[modifica | modifica wikitesto]

Controllo di autoritàGND (DE4156115-6
Matematica Portale Matematica: accedi alle voci di Wikipedia che trattano di matematica