„Konjugált gradiens módszer” változatai közötti eltérés

[nem ellenőrzött változat][nem ellenőrzött változat]
Tartalom törölve Tartalom hozzáadva
Nummatek (vitalap | szerkesztései)
Nincs szerkesztési összefoglaló
Nummat2 (vitalap | szerkesztései)
51. sor:
:<math>\mathbf{p}_0 := \mathbf{r}_0 \,</math>
:<math>k := 0 \, </math>
:'''repeatismétlés'''
::<math>\alpha_k := \frac{\mathbf{r}_k^{\mathrm{T}} \mathbf{r}_k}{\mathbf{p}_k^{\mathrm{T}} \mathbf{A p}_k} \, </math>
::<math>\mathbf{x}_{k+1} := \mathbf{x}_k + \alpha_k \mathbf{p}_k \, </math>
::<math>\mathbf{r}_{k+1} := \mathbf{r}_k - \alpha_k \mathbf{A p}_k \, </math>
::'''ifha''' '''r'''<sub>''k''+1</sub> iselegendően sufficiently smallkicsiny '''thenakkor''' exitfejezze be a loopciklust '''endvége ifha'''
::<math>\beta_k := \frac{\mathbf{r}_{k+1}^{\mathrm{T}} \mathbf{r}_{k+1}}{\mathbf{r}_k^{\mathrm{T}} \mathbf{r}_k} \, </math>
::<math>\mathbf{p}_{k+1} := \mathbf{r}_{k+1} + \beta_k \mathbf{p}_k \, </math>
::<math>k := k + 1 \, </math>
:'''endvége repeatismétlés'''
:TheAz result iseredmény: '''x'''<sub>''k''+1</sub>
 
===Példa kód a konjugált gradiens módszerre [[GNU Octave|Octave]] programnyelvben===