78553 img079 (18)

Wzór (4.4) na iteracje proste, w odniesieniu do równania (4.2), przekształconego do postaci (4.55) odpowiadającej metodzie Gaussa-Seidla, przyjmuje postać

(4.56)

*(o)^eR">

^x(*+i) = -{L + D)~^x -U ■ X(_k} + (L +D)b dla k = 0,1, 2, •••

W metodzie Gaussa-Seidla nie korzysta się jednak bezpośrednio z formuły (4.56), lecz przekształca się wzór na iteracje, otrzymując kolejno

i następnie

(L + D) •*(*+!)=-£/■*(*) +A ^D ■ ^x(k₊i) = ~^L ' ^x(k₊\) ~ U ' ^X(k) + b ,

(4.57)

(4.58)

skąd wynika następująca fonnuła obliczeniowa Gaussa-Seidla,

= -D ¹ • L • - D ^] ■ U ■ x+ D ¹ b dla k = 0,1, 2, • • • (

(4.59)

Po rozpisaniu we współrzędnych równania definiującego iteracje w metodzie Gaussa-Seidla otrzymuje się następujący układ równań

^x\,(k+i)

(^_fl12 '^x2,(k)

^aU '^X3,(k) ■■■ ^a,„-X,

ln ’^xn,(k)

^x2,(k+l)~~ ( ^a2\'^x\,(k+\) ^a23'^x3,(k) • ■ ■ ^_a2n ' ^Xn,(*))^H

-‘22

^a22

^c3,(*+l) ^:

(^—a31 '^xl,(/t+l) ^_fl32 ' ^x2,(i+l) ~ ^a2A ' ^x4,((t) ~ • ■ • ~ ^a3n '^xn.(k))⁺'

³33 ¹

^33

(4.60)

*"n,(£+l) ( ^an\'^x\,(k+\) ^an2'^x2,{k+\) ••• ^an,n-\ ' ^xn-\,(k+\)) '

dla k = 0,1, 2, .

Aby wyznaczyć współrzędne wektora je^i) z układu równań (4.60), należy wyznaczyć kolejno x_lj(jt+1) z pierwszego równania, następnie po podstawieniu już wyliczonej wartości x_]{k+\) do drugiego równania wyznaczyć X2,(t+i) i kontynuować w ten sposób obliczenia aż do wyznaczenia x_n^_k+V) z ostatniego równania, po uwzględnieniu w nim wyliczonych uprzednio wartości x_lj(i+1), x_uk+\), ■, ^x»-i#+i)-

Warunkiem koniecznym i dostatecznym zbieżności do wektora .v rozw iązania danego równania (4.2), dla dowolnego punktu początkowego jc₍₀₎ dla iteracji ciągu (x_(k))_k=Q j ₂...

kolejnych przybliżeń w metodzie Gaussa-Seidla, jest spełnienie przez promień spektralny macierzy G = ~(L + Z))^-1 • U warunku (4.30), to znaczy spełnienie nierówności

(4.61)

p({L + D)-^]-u)< 1 .

Wyszukiwarka

Podobne podstrony:
64785 img045 (37) 50 Na wstępie przedstawiono opis algorytmu iteracji prostej w zastosowaniu do znaj
DSC01524
lastscan3e 11. Stosując metodę iteracji prostej rozwiąż poniższe równanie, wymagan
18(1) 5 Ćwiczenie 7 Leżenie na plecach. Unieść głowę do góry, skręt głowy w prawo i w lewo,
OBJĘTOŚĆ BR VI. v = -^--g0l gdzie: Ogólny wzór na objętość całkowitej bryły obrotowej o równaniu
Pytania wzorcowe 4 11. Stosując metodę iteracji prostej rozwiąż poniższe równanie, wymagana dokładno
Rys. 93 c) Ponieważ płaszczyzna a jest rzutująca na rzutnię nA, rzut prostej prostopadłej do płaszcz

więcej podobnych podstron

78553 img079 (18)

78553 img079 (18)

(_fl12 'x2,(k)

(^_fl12 '^x2,(k)