Mètode de Gauss-Seidel

En àlgebra lineal numèrica, el mètode de Gauss-Seidel, també conegut com a mètode de Liebmann o mètode de desplaçament successiu, és un mètode iteratiu utilitzat per resoldre un sistema d'equacions lineals. Porta el nom dels matemàtics alemanys Carl Friedrich Gauss i Philipp Ludwig von Seidel, i és similar al mètode de Jacobi. Encara que es pot aplicar a qualsevol matriu amb elements diferents de zero a les diagonals, la convergència només està garantida si la matriu és estrictament diagonal dominant, o simètrica i definida positiva. Només es va esmentar en una carta privada de Gauss al seu alumne Christian Ludwig Gerling el 1823. Una publicació no va ser lliurada abans de 1874 per Seidel.

El mètode de Gauss-Seidel és una tècnica iterativa per resoldre un sistema quadrat de n equacions lineals amb x desconegut:$A𝐱=𝐛.$Es defineix per la iteració

$L_{*}{𝐱}^{(k+1)}=𝐛-U{𝐱}^{(k)},$

on ${𝐱}^{(k)}$ és la k- èsima aproximació o iteració de $𝐱,{𝐱}^{(k+1)}$ és la següent o k + 1 iteració de $𝐱$, i la matriu A es descompon en un component triangular inferior $L_{*}$, i un component triangular estrictament superior $U$ és a dir, $A=L_{*}+U$.

Amb més detall, escriu A, x i b en els seus components:$A=\left[\begin{array}{cccc}{a}_{11}& a_{12}& \cdots & a_{1n}\\ a_{21}& a_{22}& \cdots & a_{2n}\\ ⋮& ⋮& \ddots & ⋮\\ a_{n1}& a_{n2}& \cdots & a_{nn}\end{array}\right],𝐱=\left[\begin{array}{c}{x}_1\\ x_2\\ ⋮\\ x_n\end{array}\right],𝐛=\left[\begin{array}{c}{b}_1\\ b_2\\ ⋮\\ b_n\end{array}\right].$Aleshores, la descomposició de A en la seva component triangular inferior i la seva component triangular estrictament superior ve donada per:

$A=L_{*}+U\text{where}{L}_{*}=\left[\begin{array}{cccc}{a}_{11}& 0& \cdots & 0\\ a_{21}& a_{22}& \cdots & 0\\ ⋮& ⋮& \ddots & ⋮\\ a_{n1}& a_{n2}& \cdots & a_{nn}\end{array}\right],U=\left[\begin{array}{cccc}0& a_{12}& \cdots & a_{1n}\\ 0& 0& \cdots & a_{2n}\\ ⋮& ⋮& \ddots & ⋮\\ 0& 0& \cdots & 0\end{array}\right].$

El sistema d'equacions lineals es pot reescriure com:

$L_{*}𝐱=𝐛-U𝐱$

El mètode Gauss-Seidel ara resol el costat esquerre d'aquesta expressió per a x, utilitzant el valor anterior per a x al costat dret. Analíticament, això es pot escriure com:

${𝐱}^{(k+1)}=L_{*}^{-1}(𝐛-U{𝐱}^{(k)}).$

Tanmateix, aprofitant la forma triangular de $L_{*}$, els elements de x ^{(k +1)} es poden calcular seqüencialment mitjançant la substitució directa:

$x_i^{(k+1)}=\frac{1}{a_{ii}}(b_i-\sum _{j=1}^{i-1}{a}_{ij}{x}_j^{(k+1)}-\sum _{j=i+1}^n{a}_{ij}{x}_j^{(k)}),i=1,2,\dots ,n.$