5. Rozwiązywanie układów równań liniowych Wprowadzenie (5.1) Układ równań z niewiadomymi: Wprowadzając oznaczenia macierzowe układ równań można zapisać w postaci n n AX = A 0 . a 11 x 1 +a 12 x 2 +... +a 1n x n =a 10 , a 21 x 1 + a 22 x 2 +... +a 2n x n =a 20 , .............................................. , a n1 x 1 + a n2 x 2 +...+ a nn x n = a n0 .
32
Embed
5. Rozwiązywanie układów równań liniowych · 2014-10-02 · Metody skończone rozwiązywania układów równań liniowych (5.2) Metoda Eliminacji Gaussa (5.2.2) Rozwiązywanie
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Zakłada się, że elementy są liczbami rzeczywistymi. Macierz nazywa się macierzą główną układu równań, wektor kolumną wyrazów wolnych a wektor wektorem niewiadomych. Rozwiązaniem rozważanego układu równań nazywany jest każdy wektor
o składowych rzeczywistych, spełniający równanie
Twierdzenie (Cramer): Układ równań liniowych ma dokładnie jedno rozwiązanie gdy
Metody rozwiązywania równań liniowych dzieli się na dwie klasy:
- metody skończone - pozwalają uzyskać rozwiązanie po wykonaniu skończonej liczby działań arytmetycznych a otrzymane rozwiązania obarczone są tylko błędami zaokrągleń.
- metody iteracyjne - polegają na wyznaczaniu ciągu wektorów zbieżnego do rozwiązania układu. Uzyskane rozwiązania są obarczone błędem metody oraz zaokrągleń. Metody te pozwalają jednak wyznaczyć rozwiązania z dowolną, z góry zadaną dokładniością.
A X̄=A0 .
X 0 , X 1 , , X k ,
a ij Aa ijA0
X
x j , j=1 n ,
AX=A0det A≠0.
X=[ x1 x2 xn ]T
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda wyznaczników Cramera (5.2.1)
Układ równań liniowych można rozwiązać stosując wzory Cramera
Def. Podwyznacznikiem danego wyznacznika nazywamy każdy wyznacznik, który otrzymamy usuwając z macierzy danego wyznacznika pewną liczbę wierszy i taką samą liczbę kolumn, zachowując kolejność elementów pozostałych.
Def. Minorem wyznacznika przynależnym do elementu macierzy tego wyznacznika nazywamy ten podwyznacznik danego wyznacznika, który otrzymamy usuwając z macierzy danego wyznacznika wiersz oraz kolumnę, na przecięciu których znajduje się ten element (usuwamy i-ty wiersz oraz k-tą kolumnę).
AX=A0
a ik
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda wyznaczników Cramera (5.2.1)
Uwaga! Minor wyznacznika przynależny do elementu jest minorem głównym wyznacznika
Def. Dopełnieniem algebraicznym elementu wyznacznika nazywamy iloczyn minora tego wyznacznika przynależnego do elementu oraz czynnika
Tw. (Laplace'a) Wyznacznik jest równy sumie iloczynów każdego elementu dowolnego wiersza (kolumny) i odpowiadającego temu elementowi dopełnienia algebraicznego.
Przykład 1: Stosując wzory Cramera wyprowadzić wzory na rozwiązanie układu dwóch równań liniowych z dwiema niewiadomymi
a11 x1+a12 x2=b1 ,a21 x1+a22 x2=b2 .
x1, x2
Aik
a ii
a ika ik
−1ik .
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda wyznaczników Cramera (5.2.1)
Rozwiązanie:
Przykład 2: Stosując wzory Cramera wyprowadzić wzory na rozwiązanie układu trzech równań liniowych z trzema niewiadomymi x1, x2, x3
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda wyznaczników Cramera (5.2.1) Rozwiązanie:
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda wyznaczników Cramera (5.2.1)
Złożoność obliczeniowa metody wyznaczników Cramera:
Stosowanie wzorów Cramera nakłada konieczność obliczenia wyznaczników stopnia . Stosując opisaną wcześniej klasyczną metodę obliczania wyznacznika należy wykonać dodawań oraz mnożeń. Łącznie należy wykonać działań arytmetycznych zmiennoprzecinko-wych. W tabeli poniżej zamieszczono czasy obliczania wyznacznika przez komputer wykonujący operacji zmiennoprzecinkowych na sekundę.
Aby rozwiązać układ równań należy obliczyć wyznaczni- ków stopnia oraz wykonać dzieleń. Całkowita liczba działań arytmetycz- nych jest zatem równa Z powodu dużej złożoności obliczeniowej wzory Cramera stosuje się tylko dla
n1n
n! n⋅n !n !n1!
108
AX=A0 n1n n
n1⋅n1 !n .n4.
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Rozwiązywanie układu równań przebiega w dwóch etapach:
- etap postępowania prostego - układ równań zostaje przekształcony do postaci równoważnej (mającej takie samo rozwiązanie, co układ wyjściowy)
z trójkątną górną macierzą główną układu,
- etap postępowania odwrotnego - układ w postaci równoważnej zostaje rozwiązany.
Metoda Eliminacji Gaussa - etap postępowania prostego
Wyjściowy układ równań:
(1)
Przekształcanie macierzy głównej układu w macierz trójkątną odbywa sie w krokach.
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Etap postępowania prostego (cd)
Krok 1: eliminacja niewiadomej z równań nr
Zakładając, że element dzieli się pierwsze równanie przez i otrzymuje
(2)
gdzie
Następnie od i-tego ( ) równania układu (1) odejmuje się równanie (2) pomnożone przez otrzymując
(3)
gdzie
a110 ≠0
x1
a i10
2,3, , n .
a110
i=2,3, , n
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Etap postępowania prostego (cd)
Krok 2: eliminacja niewiadomej z równań nr Schemat obliczeń kroku 1 powtarza się w odniesieniu do ukłdu równań (3). Tak więc zakładając, że element dzieli się pierwsze równanie układu (3) przez i otrzymuje
(4)
gdzie
Następnie od i-tego ( ) równania układu (3) odejmuje się równanie (4) pomnożone przez otrzymując
(5)
gdzie
a221≠0
x2
a i21
3,4,, n.
a221
i=3,4, , n
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Etap postępowania prostego (cd)
Krok 3,4, ..., n-1: eliminacja kolejnych niewiadomych w sposób analogiczny do algorytmu postępowania z kroków 1 i 2. Po wykonaniu obliczeń w kroku, otrzymuje się równania
(6) gdzie
oraz (7)
Krok n: zakłada się, że element Równanie (7) dzieli się przez i otrzymuje
(8)
n−1
x3 , , xn−1
an nn−1≠0 .
an nn−1
xn=an0n .
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Etap postępowania odwrotnego:
Rozwiącuje się równoważny układ równań uzyskany w wyniku postępowania prostego
Algorytm rozwiązania tego układu równań jest oczywisty. Z ostatniego równania odczytuje się bezpośrednio Wartość tę podstawia się do równań Z równania można wyznaczyć Obliczoną wartość podstawia sie do równań Teraz z równania można wyznaczyć itd.
xn .n−1,n−2, ,1 . xn−1 .n−1
n−2, n−3, ,1 .n−1,n−2, ,1 .
n−2xn−2
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Wnioski: Układ równań liniowych mozna rozwiązać metodą eliminacji Gaussa, gdy elementy są różne od zera. Ponieważ pomiędzy tymi elementami a minorami macierzy głównej układu zachodzą związki wyrażone w twierdzeniu Twierdzenie. Zachodzą równości
Stąd prawdziwe jest twierdzenie Twierdzenie. Układ równań liniowych można rozwiązać metodą eliminacji Gaussa, jeśli wszystkie minory główne macierzy głównej układu są różne od zera
ak k(k−1 ) (k=1 , 2 ,… , n)
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Przykład : Stosując metodę eliminacji Gaussa rozwiązać układ równań liniowych
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Z pełnym wyborem elementu podstawowego
Przytoczone wczesniej twierdzenie określające warunki stosowalności metody eliminacji Gaussa pokazało istotne jej ograniczenia (aby układ równań liniowych można było rozwiązać metodą eliminacji Gaussa, wszystkie minory główne macierzy głównej układu muszą być różne od zera).
Algorytm ten można jednak zmodyfikować w taki sposób, aby można go było stosować do układu Cramera, tzn do układu z
Przyjmijmy, że układ równań liniowych jest rozwiązywany metodą eliminacji Gaussa i że zostało już wykonanych kroków etapu postępowania prostego. Wyjściowy układ równań został przekształcony do postaci
(9)
det A≠0.
k−1AX=A0
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Algorytm z pełnym wyborem elementu podstawowego
1. wyszukuje się element spełniający warunek
2. przestawia się w układzie (9) równanie z równaniem oraz kolumę kolumną ,
3. eliminuje się niewiadomą z równań zgodnie z algorytmem k-tego kroku prostej eliminacji Gaussa.
Jeśli to żaden element podstawowy w tej metodzie nie jest równy zero.
Przykład : Rozwiązać układ równań liniowych stosując metodę eliminacji Gaussa z pełnym wyborem elementu podstawowego
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Algorytm z pełnym wyborem elementu podstawowego
Rozwiązanie: Postępowanie proste
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Algorytm z pełnym wyborem elementu podstawowego
Rozwiązanie (cd): Postępowanie proste
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Algorytm z pełnym wyborem elementu podstawowego
Rozwiązanie (cd): Postępowanie proste
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Algorytm z pełnym wyborem elementu podstawowego
Proces zamiany elementów w macierzy oraz wektorach i najprościej można dokonać możąc powyższe zmienne przez odpowiednio dobrane macierze elementarnej permutacji .
Macierz taka różni się od macierzy jednostkowej tylko elementami i , stąd pełna informacja o macierzy dana jest tylko dwiema liczbami i . Można pokazać, że:
czyli macierz elementarnej permutacji jest symetryczna i ortogonalna.
Lewostronne przemnożenie danej macierzy przez (t.j. ) odpowiada przestawieniu jej wierszy o numerach i , a przemnożenie prawo- stronne (t.j. ) równoważne jest przestawieniu odpowiednich kolumn.
A0A X
Pk ,ik
Pk ,ik=[1..0 . . 1. .. .1 . . 0
..1
]pk ,k=p ik , ik=0 pk ,ik=p ik , k=1k ik
Pk ,ik=Pk ,i kT =Pk ,i k
−1
A Pk ,ik Pk ,ik A
k ikAPk ,i k
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Algorytm z pełnym wyborem elementu podstawowego
W omawianym przykładzie mamy:
W kroku 1 postępowania prostego należy zamienić miejscami wiersze
oraz . Macierz elementarnej
permutacji ma wówczas postać:
Macierz oraz wektor z zamienionymi wierszami otrzymamy w wyniku operacji
Metody skończone rozwiązywania układów równań liniowych (5.2)
Metoda Eliminacji Gaussa (5.2.2)
Algorytm z pełnym wyborem elementu podstawowego
W kroku 4 (ostatnim) postępowania prostego należy jedynie sprowadzić do współczynnik macierzy . Ostatecznie po zakończeniu postępowania prostego otrzymujemy układ równań z macierzą główną trójkątną górną i wektorem
Teraz można rozwiązać powyższy układ równań - rozpocząć etap postępowa-nia odwrotnego algorytmu eliminacji Gaussa. Ostatecznie otrzymuje się wektor , którego elementy należy odpowiednio przepisać do wektora .