Yakobi üsulu

Yakobi üsulu — rəqəmsal xətti cəbrdə diaqonal dominant xətti bərabərliklərin həllinin tapılması alqoritmi. Hər bir diaqonal element həll edilir və təxmini dəyər daxil edilir. Proses həllə yaxınlaşana kimi davam etdirilir. Bu üsula Karl Qustav Yakob Yakobinin adı verilib.

Təsviri

Fərz edək ki,

A\mathbf {x} =\mathbf {b}

n dərəcəli xətti bərabərliklərdir, burada:

$A={\begin{bmatrix}a_{11}&a_{12}&\cdots &a_{1n}\\a_{21}&a_{22}&\cdots &a_{2n}\\\vdots &\vdots &\ddots &\vdots \\a_{n1}&a_{n2}&\cdots &a_{nn}\end{bmatrix}},\qquad \mathbf {x} ={\begin{bmatrix}x_{1}\\x_{2}\\\vdots \\x_{n}\end{bmatrix}},\qquad \mathbf {b} ={\begin{bmatrix}b_{1}\\b_{2}\\\vdots \\b_{n}\end{bmatrix}}.$

Sonra A matrisi diaqonal D komponentinə və onun qalığı R matrisinə bölünür:

A=D R\qquad {\text{where}}\qquad D={\begin{bmatrix}a_{11}&0&\cdots &0\\0&a_{22}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &a_{nn}\end{bmatrix}}{\text{ and }}R={\begin{bmatrix}0&a_{12}&\cdots &a_{1n}\\a_{21}&0&\cdots &a_{2n}\\\vdots &\vdots &\ddots &\vdots \\a_{n1}&a_{n2}&\cdots &0\end{bmatrix}}.

Bunun həlli təkrarlanmaqla belə tapılır

\mathbf {x} ^{(k 1)}=D^{-1}(\mathbf {b} -R\mathbf {x} ^{(k)}),

burada $\mathbf {x} ^{(k)}$ , $\mathbf {x}$ -nin k dərəcəli approksimasiyası yaxud təkrarlanması və $\mathbf {x} ^{(k 1)}$ , $\mathbf {x}$ -nin növbəti yaxud k 1 dərəcəli təkrarlanmasıdır. Element əsaslı formula beləcə aşağıdakı kimidir:

x_{i}^{(k 1)}={\frac {1}{a_{ii}}}\left(b_{i}-\sum _{j\neq i}a_{ij}x_{j}^{(k)}\right),\quad i=1,2,\ldots ,n.

x_i^{(k 1)} hesablanması x^(k)-də özündən başqa hər bir elementin olmasını tələb edir.

Nümunə

Xətti bərabərlik sistemi $Ax=b$ formasında və onun ilkin fərz edilən həlli $x^{(0)}$ verilib

A={\begin{bmatrix}2&1\\5&7\\\end{bmatrix}},\ b={\begin{bmatrix}11\\13\\\end{bmatrix}}\quad {\text{and}}\quad x^{(0)}={\begin{bmatrix}1\\1\\\end{bmatrix}}.

Biz $x$ hesablamaq üçün yuxarıda verilən $x^{(k 1)}=D^{-1}(b-Rx^{(k)})$ bərabərliyindən istifadə edirik. Əvvəlcə biz bərabərliyi daha rahat olan $D^{-1}(b-Rx^{(k)})=Tx^{(k)} C$ formasında yazırıq, burada $T=-D^{-1}R$ və $C=D^{-1}b$ . Nəzərə alın ki, $R=L U$ , burada $L$ və $U$ , $A$ matrisinin aşağı və yuxarı hissələridir. Verilən dəyərlərə əsasən