Метод релаксации

Метод релаксации - приближенный метод решения систем линейных уравнений.

Система линейных уравнений

${\displaystyle \left\{{\begin{matrix}a_{11}x_{1}+\ldots +a_{1n}x_{n}&=&b_{1}\\a_{21}x_{1}+\ldots +a_{2n}x_{n}&=&b_{2}\\&\ldots &\\a_{n1}x_{1}+\ldots +a_{nn}x_{n}&=&b_{n}\\\end{matrix}}\right.}$

приводится к виду

${\displaystyle \left\{{\begin{matrix}-x_{1}+b_{12}x_{2}+\ldots +b_{1n}x_{n}+c_{1}&=&0\\&\ldots &\\b_{n1}x_{1}+b_{n2}x_{2}+\ldots -x_{n}+c_{n}&=&0\\\end{matrix}}\right.}$

где ${\displaystyle b_{ij}=-{\frac {a_{ij}}{a_{ii}}}}$, ${\displaystyle (i\neq j)}$, ${\displaystyle c_{i}={\frac {b_{i}}{a_{ii}}}}$

Находятся невязки ${\displaystyle R_{j}}$:

${\displaystyle \left\{{\begin{matrix}R_{1}^{(0)}&=&c_{1}-x_{1}^{(0)}+\sum _{j=2}^{n}b_{1j}x_{j}^{(0)}\\R_{2}^{(0)}&=&c_{2}-x_{2}^{(0)}+\sum _{j=1,j\neq 2}^{n}b_{2j}x_{j}^{(0)}\\&\ldots &\\R_{n}^{(0)}&=&c_{n}-x_{n}^{(0)}+\sum _{j=1}^{n-1}b_{nj}x_{j}^{(0)}\\\end{matrix}}\right.}$

Выбирается начальное приближениее ${\displaystyle X^{(0)}=0}$. На каждом шаге необходимо обратить в ноль максимальную невязку: ${\displaystyle R_{s}^{(k)}=\delta x_{s}^{(k)}\Rightarrow R_{s}^{(k+1)}=0,R_{i}^{(k+1)}=R_{i}^{(k)}+b_{is}\delta x_{s}^{(k)}}$.

Условие остановки: ${\displaystyle R_{j}^{(k)}|<\varepsilon ,\forall j={\overline {1,n}}}$.

Ответ находится по формуле: ${\displaystyle x_{i}\approx x_{i}^{(0)}+\sum _{j}\delta x_{i}^{(j)}}$.