plik


WykBad trzeci 1 Metoda Gaussa-Seidla EiT, sem. 2, 2014/2015 2 Ax = b x = Gx + c A = AL + AD + AU a11 0 . 0 0 0 . 0 0 a12 . a1n a 0 0 . a2n 0 a22 . 0 0 . 0 21 AD = AL = AU = . . . . . . . . . . . . a an2 . 0 0 0 . ann n1 0 0 . 0 (AL + AD)x = - AU x + b (AL + AD + AU)x = b (AL + AD)x(k+1) = - AU x(k) + b k = 0, 1, & AD x(k+1) = -AD-1 AL x(k+1) - AD-1 AU x(k ) + AD-1 b k = 0,1,... 3 Macierz G a12 a13 a1n b1 x1,(k+1) = - x2,(k) - x3,(k) -K- xn,(k) + a11 a11 a11 a11 a21 a23 a2n b2 x2,(k+1) = - x3,(k) -K- xn,(k) + - a22 x1,(k+1) a22 a22 a22 k = 0,1,... a31 a32 a34 a3n b3 x3,(k+1) = - x4,(k) -K- xn,(k) + - a33 x1,(k+1) - a33 x2,(k+1) a33 a33 a33 .. ....... .. . .. .. ......... . .. .. ..... .. .. .. . .. .. ..... .. .. .. . .. .. ..... .. .. .. . .. .. ..... an1 an2 ann-1 bn xn,(k+1) = - xn-1,(k+1) + - ann x1,(k+1) - ann x2,(k+1) -K ann ann bi c = i = 1, 2,..., n wektor c 4 aii x(0) = c Zaczynamy obliczenia od przyjcia, |e KoDczymy obliczenia, gdy x(k+1) - x(k ) < e x(k +1) = x* Rozwizanie 5 g11 g12 g1n g g22 g2n 21 G = gn2 gnn n1 g Warunki zbie|no[ci n gij < 1 j = 1,2,...,n i=1 n gij < 1 i = 1,2,...,n j=1 6

Wyszukiwarka

Podobne podstrony:
Pierwszy wyklad 14?z tła
Drugi wykład 14?z tła
Czwarty wykład 14?z tła
Czwarty wykład 14?z tła
Czwarty wykład? 2014?z tła
Czwarty wykład? 2014?z tła
Szósty wykład 14 bez tła
Wyklad 2?kultet?rmatozy tla autoimmunologicznego
Sieci komputerowe wyklady dr Furtak
Wykład 05 Opadanie i fluidyzacja
WYKŁAD 1 Wprowadzenie do biotechnologii farmaceutycznej

więcej podobnych podstron