Metoda redukcji Gaussa – Jordana, Metody numeryczne Scilab


Metoda redukcji Gaussa - Jordana.

     

    Redukcja Gaussa - Jordana  jest innym sposobem eliminacji niewiadomych, niż w opisanej wyżej metodzie Gaussa. Różnica  polega  na tym, że redukcję robimy w całej k-tej kolumnie, a nie tylko poniżej elementu głównego. Układ równań postaci

                                   0x01 graphic
                            (1.36)

przekształcamy w następujący sposób: Pierwsze równanie dzielimy obustronnie przez 0x01 graphic
, a następnie od i-tego wiersza, i=2, 3, ..., n, odejmujemy wiersz pierwszy pomnożony przez 0x01 graphic
. Otrzymujemy układ równań postaci

                                      0x01 graphic
                                (1.37)

Następnie drugie równanie dzielimy obustronnie przez 0x01 graphic
 i od  i-tego wiersza,

i = 1, 3, 4, ..., n,  odejmujemy wiersz drugi, pomnożony przez 0x01 graphic
, otrzymując w ten sposób układ równań

                                      0x01 graphic
                               (1.38)

Po  wykonaniu (n-1) eliminacji otrzymujemy układ równań postaci

                                                0x01 graphic
                                          (1.39)

zatem  układ   (1.39)  jest rozwiązaniem układu równań  (1.36).

Metoda ta umożliwia rozwiązywanie układu równań przy oszczędniejszym wykorzystaniu pamięci operacyjnej, niż w metodzie Gaussa - mniejsza jest liczba współczynników układu, które w czasie obliczeń muszą jednocześnie znajdować się w pamięci komputera. Można zatem metodą redukcji Gaussa - Jordana rozwiązać większe układy równań [4].

Przykład 3.     Rozwiązywanie układów równań metodą redukcji Gaussa - Jordana

Rozwiązać układ równań:

                                              0x01 graphic
                                             (2.8)

Pierwsze równanie układu (2.8) dzielimy obustronnie przez współczynnik przy x1 w pierwszym równaniu ( a11=1).

W wyniku dzielenia otrzymujemy równanie

                                                    0x01 graphic
                                                    (2.9)

Następnie mnożymy równanie (2.9) przez  współczynnik przy x1 w drugim równaniu (a21=4). W wyniku tego mnożenia otrzymujemy równanie

                                                    0x01 graphic
                                              (2.10)

Od drugiego równania układu (2.8) odejmujemy stronami otrzymane w wyniku mnożenia równanie (2.10). Otrzymujemy równanie z wyeliminowaną zmienną x1.

                                                       0x01 graphic
                                               (2.11)

Następnie przystępujemy do eliminacji zmiennej x1 z trzeciego równania układu (2.8). Podobnie jak poprzednio, mnożymy równanie (2.9), ale tym razem przez współczynnik przy x1 w trzecim równaniu układu (2.8)  ( przez a31=1). Otrzymujemy równanie

                                                      0x01 graphic
                                               (2.12)

Od trzeciego równania układu (2.8) odejmujemy stronami otrzymane w wyniku mnożenia równanie (2.12). Otrzymujemy równanie

                                                     0x01 graphic
                                                  (2.13)

Równania  (2.9), (2.11) i (2.13) tworzą pierwszy układ zredukowany z wyeliminowaną zmienną  x1 z drugiego i trzeciego równania.

                                            0x01 graphic
                                         (2.14)

Przystępujemy teraz do eliminacji  zmiennej x2 z pierwszego i trzeciego równania. Dzielimy w tym celu stronami drugie równanie układu równań (2.14) przez współczynnik przy x2 czyli a22= -5 . Otrzymujemy w wyniku tego działania równanie

                                                       0x01 graphic
                                                    (2.15)

Następnie mnożymy równanie (2.15) przez współczynnik przy x2 w pierwszym równaniu układu (2.14) czyli a12=2. Otrzymujemy w wyniku mnożenia równanie

                                                 0x01 graphic
                         (2.16)

Odejmujemy równanie (2.16) od pierwszego równania układu równań (2.14). Otrzymujemy równanie

                                                 0x01 graphic
                                                        (2.17)

Eliminujemy teraz  x2  z trzeciego równania. Mnożymy w tym celu stronami równanie (2.15) przez współczynnik przy x2 w trzecim równaniu, czyli a32=-3. W wyniku otrzymujemy równanie

                                               0x01 graphic
                                                   (2.18)

odejmujemy teraz stronami równanie  (2. 18) od trzeciego równania układu (2.14). Otrzymujemy równanie

                                                       0x01 graphic
                                                          (2.19)

Równania (2.17),  (2.15) oraz  (2.19) tworzą drugi układ zredukowany z wyeliminowaną zmienną x2 z pierwszego i trzeciego równania układu równań.

                                         0x01 graphic
                                                 (2.20)

Należy teraz wyeliminować zmienną  xz pierwszego i drugiego równania układu równań (2.20). Trzecie równanie układu (2.20) dzielimy obustronnie przez współczynnik przy x3, czyli  a33= 0x01 graphic
. Otrzymujemy równanie             

                                                               0x01 graphic
                                                              (2.21)

Równanie  (2.21) mnożymy przez współczynnik przy x3 w pierwszym równaniu układu równań (2.20) czyli 0x01 graphic
. Następnie otrzymane w ten sposób równanie odejmujemy od równania pierwszego układu (2.20). W wyniku odejmowania otrzymujemy równanie

                                                                    0x01 graphic
                                                               (2.22)

Mnożymy następnie równanie (2.21) przez współczynnik przy x3 w drugim równaniu układu równań (2.20), czyli 0x01 graphic
. Otrzymane równanie odejmujemy stronami od równania drugiego układu (2.20).  Otrzymujemy w wyniku równanie

                                                                 0x01 graphic
                                                               (2.23)

Równania (2.21), (2.22) oraz (2.23) tworzą trzeci układ zredukowany z wyeliminowaną zmienną x3 z pierwszego i drugiego równania

                                              0x01 graphic
                                                       (2.24)

Jak widać układ równań (2.24) jest rozwiązaniem układu równań  (2.8) czyli rozwiązaniem naszego zadania.



Wyszukiwarka

Podobne podstrony:
Metoda eleminacji Gaussa, Metody numeryczne Scilab
Metoda Jacobiego, Metody numeryczne Scilab
Metoda Eulera, Metody numeryczne Scilab
Metoda Rungego – Kutty, Metody numeryczne Scilab
metoda siecznych, Elektrotechnika, SEM3, Metody numeryczne, egzamin metody numeryczn
METODA BAIRSTOWA, Politechnika, Lab. Metody numeryczne
Cw 9 DUO, Politechnika Lubelska, Studia, Studia, sem III, pen, METODY NUMERYCZNE, Scilab
metoda siecznych, Elektrotechnika, SEM3, Metody numeryczne, egzamin metody numeryczn
METODA BAIRSTOWA, Politechnika, Lab. Metody numeryczne
Metoda redukcyjna Gaussa
metoda el gaussa, sprawozdania PWR, metody numeryczne w5
Metody numeryczne Metoda węzłowa
metoda regula falsi, Elektrotechnika, SEM3, Metody numeryczne, egzamin metody numeryczn
barcz,metody numeryczne, metoda iteracji prostych
Metody numeryczne, Metoda Eulera, LABORATORIUM Z
Metody numeryczne, newton 1, Metoda ta służy do obliczenia przybliżonej wartości pierwiastka równani
metoda grupowa, gik, gik, I sem, zz przodki, II sem, numerki, od chłopaków, metody numeryczne, metod
Metoda RK sprawko, Automatyka i robotyka air pwr, VI SEMESTR, Metody numeryczne
Metody numeryczne, metoda Rungego-Kutte grzesiek kucharczyk, Akademia Górniczo-Hutnicza

więcej podobnych podstron