Wyk lad 2
Uk lady r´ owna´ n liniowych
1 Podstawowe poj ecia zwi
,azane z uk ladami r´
,owna´ n liniowych
Uk ladem m r´ owna´ n liniowych o niewiadomych x
1, x
2, . . ., x
nnazywamy uk lad r´ owna´ n postaci:
a
11x
1+ a
12x
2+ . . . + a
1nx
n= b
1a
21x
1+ a
22x
2+ . . . + a
2nx
n= b
2. . . . . . . . . . . . . . . . . . . . . a
m1x
1+ a
m2x
2+ . . . + a
mnx
n= b
m, (1)
gdzie wsp´ o lczynniki a
ij(dla i = 1, . . . , m; j = 1, . . . , n) oraz elementy b
i(dla i = 1, . . . , m) nale˙z a do cia la R. Uk lad ten nazywamy jednorodnym, gdy b
, 1= b
2= . . . = b
m= 0.
Macierz a wsp´
,o lczynnik´ ow uk ladu (1) nazywamy macierz:
A =
a
11a
12. . . a
1na
21a
22. . . a
2n.. . .. . . .. .. . a
m1a
m2. . . a
mn
, (2)
za´ s macierz a uzupe lnion
,a uk ladu (1) nazywamy macierz:
,A
u=
a
11a
12. . . a
1nb
1a
21a
22. . . a
2nb
2.. . .. . . .. .. . .. . a
m1a
m2. . . a
mnb
m
. (3)
Rozwi azaniem uk ladu (1) nazywamy taki ci
,ag liczbowy (c
, 1, c
2, . . . , c
n), ˙ze po zast apieniu
,w r´ ownaniach tego uk ladu niewiadomych x
ielementami c
idla i = 1, 2, . . . , n otrzymujemy r´ owno´ sci prawdziwe w ciele R, tj. gdy A · C = B, gdzie B jest kolumn a wyraz´
,ow wolnych tzn.
B =
b
1b
2.. . b
m
oraz C =
c
1c
2.. . c
n
. (4)
Wynika st ad, ˙ze przy tych oznaczeniach uk lad (1) mo˙zna zapisa´
,c w postaci macierzowej
A · X = B. (5)
Definicja 2.1. Powiemy, ˙ze r´ ownanie
a
1x
1+ a
2x
2+ . . . + a
nx
n= b (6)
jest kombinacj a liniow
,a r´
,owna´ n uk ladu (1), je˙zeli istniej a c
, 1, c
2, . . . , c
m∈ R (zwane wsp´ o l- czynnikami tej kombinacji) takie, ˙ze po pomno˙zeniu stronami i-tego r´ ownania przez c
idla i = 1, . . . , m i dodaniu stronami otrzymanych r´ owna´ n uzyskamy r´ ownanie (1), tzn.
b =
m
X
i=1
c
ib
ioraz a
j=
m
X
i=1
c
ia
ijdla j = 1, 2, . . . , n. (7) Uwaga 2.2. Zauwa˙zmy, ˙ze je´ sli r´ ownanie (6) jest kombinacj a liniow
,a pewnych r´
,owna´ n uk ladu (1), to jest ono tak˙ze kombinacj a liniow
,a wszystkich r´
,owna´ n tego uk ladu (brakuj ace
,wsp´ o lczynniki s a r´
,owne 0).
Twierdzenie 2.3. Ka˙zde rozwi azanie uk ladu (1) jest rozwi
,azaniem ka˙zdego r´
,ownania b ed
,a-
,cego kombinacj a liniow
,a r´
,owna´ n uk ladu (1).
Dow´ od. Za l´ o˙zmy, ˙ze r´ ownanie (6) jest kombinacj a liniow
,a o wsp´
,o lczynnikach c
1, c
2, . . . , c
mr´ owna´ n uk ladu (1) i niech (p
1, p
2, . . . , p
n) b edzie rozwi
,azaniem uk ladu (1). Wtedy
,
a
11p
1+ a
12p
2+ . . . + a
1np
n= b
1a
21p
1+ a
22p
2+ . . . + a
2np
n= b
2. . . . . . . . . . . . . . . . . . . . . . . a
m1p
1+ a
m2p
2+ . . . + a
mnp
n= b
m,
wi ec po pomno˙zeniu obu stron i-tej r´
,owno´ sci przez c
idla i = 1, . . . , m i dodaniu stronami otrzymanych r´ owno´ sci uzyskamy na mocy wzor´ ow (7), ˙ze a
1p
1+ a
2p
2+ . . . + a
np
n= b, czyli (p
1, p
2, . . . , p
n) jest rozwi azaniem r´
,ownania (6). 2
Definicja 2.4. Dwa uk lady r´ owna´ n liniowych (U
1) i (U
2) z n niewiadomymi x
1, . . . , x
nnazywamy r´ ownowa ˙znymi, gdy ka˙zde r´ ownanie uk ladu (U
1) jest kombinacj a liniow
,a r´
,owna´ n uk ladu (U
2) i odwrotnie. Piszemy wtedy (U
1) ≡ (U
2).
Z twierdzenia 2.3 wynika od razu nast epuj
,ace
,Twierdzenie 2.5. R´ ownowa˙zne uk lady r´ owna´ n liniowych posiadaj a identyczne zbiory
,rozwi aza´
,n.
Definicja 2.6. Uk lad r´ owna´ n liniowych z n niewiadomymi x
1, x
2, . . . , x
nnazywamy sprzecz- nym, gdy r´ ownanie 0 · x
1+ 0 · x
2+ . . . + 0 · x
n= 1 jest kombinacj a liniow
,a r´
,owna´ n tego uk ladu.
Poniewa˙z r´ ownanie 0 · x
1+ 0 · x
2+ . . . + 0 · x
n= 1 nie posiada rozwi azania, wi
,ec z powy˙zszej
,definicji oraz z twierdzenia 2.3 wynika od razu nast epuj
,ace
,Twierdzenie 2.7. Sprzeczny uk lad r´ owna´ n liniowych nie posiada rozwi azania.
,Lemat 2.8. Za l´ o˙zmy, ˙ze i-te r´ ownanie dla i = 1, . . . , k uk ladu r´ owna´ n
a
011x
1+ a
012x
2+ . . . + a
01nx
n= b
01a
021x
1+ a
022x
2+ . . . + a
02nx
n= b
02. . . . . . . . . . . . . . . . . . . . . . . a
0k1x
1+ a
0k2x
2+ . . . + a
0knx
n= b
0k, (8)
jest kombinacj a liniow
,a r´
,owna´ n uk ladu (1) o wsp´ o lczynnikach c
i1, c
i2, . . . , c
im. Je˙zeli r´ ownanie (6) jest kombinacj a liniow
,a r´
,owna´ n uk ladu (8) o wsp´ o lczynnikach c
1, c
2, . . . , c
k, to r´ ownanie (6) jest kombinacj a liniow
,a r´
,owna´ n uk ladu (1) o wsp´ o lczynnikach
k
X
i=1
c
ic
i1,
k
X
i=1
c
ic
i2, . . . ,
k
X
i=1
c
ic
im.
Dow´ od. Na mocy (7) mamy, ˙ze dla i = 1, 2, . . . , k: a
0ij=
m
X
t=1
c
ita
tjdla j = 1, 2, . . . , n
oraz b
0i=
m
X
t=1
c
itb
t. Ponadto a
j=
k
X
i=1
c
ia
0ijdla j = 1, 2, . . . , n oraz b =
k
X
i=1
c
ib
0i. Zatem dla j = 1, 2, . . . , n:
a
j=
k
X
i=1
c
i·
m
X
t=1
c
ita
tj!
=
k
X
i=1 m
X
t=1
c
ic
ita
tj=
m
X
t=1 k
X
i=1
c
ic
ita
tj=
m
X
t=1 k
X
i=1
c
ic
it!
a
tjoraz
b =
k
X
i=1
c
i·
m
X
t=1
c
itb
t!
=
k
X
i=1 m
X
t=1
c
ic
itb
t=
m
X
t=1 k
X
i=1
c
ic
itb
t=
m
X
t=1 k
X
i=1
c
ic
it!
b
t, sk ad na mocy
,(7) mamy tez e.
,2
Z lematu wynika od razu, ˙ze je´ sli r´ ownanie 0 · x
1+ 0 · x
2+ . . . + 0 · x
n= a, gdzie a 6= 0 jest kombinacj a liniow
,a r´
,owna´ n uk ladu (1), to uk lad ten jest sprzeczny. Ponadto z lematu wynikaj a
,od razu nast epuj
,ace twierdzenia:
,Twierdzenie 2.9. Za l´ o˙zmy, ˙ze uk lady r´ owna´ n liniowych (U
1) i (U
2) z n niewiadomymi x
1, x
2, . . . , x
ns a r´
,ownowa˙zne. W´ owczas uk lad (U
1) jest sprzeczny wtedy i tylko wtedy, gdy uk lad (U
2) jest sprzeczny.
Twierdzenie 2.10. Niech (U
1), (U
2), (U
3) b ed
,a uk ladami r´
,owna´ n liniowych z n niewiado- mymi x
1, x
2, . . . , x
n. W´ owczas:
(i) (U
1) ≡ (U
1);
(ii) je˙zeli (U
1) ≡ (U
2), to (U
2) ≡ (U
1);
(iii) je˙zeli (U
1) ≡ (U
2) i (U
2) ≡ (U
3), to (U
1) ≡ (U
3).
Problem rozwi azania uk ladu r´
,owna´ n liniowych polega na znalezieniu wszystkich rozwi aza´
,n tego uk ladu. Bardzo u˙zyteczne przy rozwi azywaniu tego problemu s
,a tzw. operacje elemen-
,tarne.
2 Operacje elementarne nad uk ladem r´ owna´ n liniowych
(i). Zamiana miejscami r´ ownania i-tego z r´ ownaniem j-tym (i 6= j) oznaczana przez r
i↔ r
j. Operacja odwrotna: r
i↔ r
j.
(ii). Pomno˙zenie i-tego r´ ownania przez niezerowy skalar a oznaczana przez a · r
i. Operacja odwrotna:
1a· r
i.
(iii). Dodanie do i-tego r´ ownania r´ ownania j-tego (i 6= j) pomno˙zonego przez dowolny skalar a
oznaczana przez r
i+a·r
j. Przy tej operacji zmieniamy tylko r´ ownanie i-te! Operacja odwrotna:
r
i+ (−a) · r
j.
(iv). Wykre´ slenie powtarzaj acych si
,e kopii pewnego r´
,ownania.
(v). Wykre´ slenie r´ owna´ n postaci 0 · x
1+ 0 · x
2+ . . . + 0 · x
n= 0 (gdy liczba r´ owna´ n jest wi eksza
,od 1).
(vi). Zamiana kolejno´ sci niewiadomych x
ioraz x
j(i 6= j) w ka˙zdym r´ ownaniu oznaczana przez x
i↔ x
j. W wyniku zastosowania tej operacji r´ ownanie
a
1x
1+ . . . + a
ix
i+ . . . + a
jx
j+ . . . + a
nx
n= b przechodzi na r´ ownanie
a
1x
1+ . . . + a
jx
j+ . . . + a
ix
i+ . . . + a
nx
n= b.
Z definicji uk lad´ ow r´ ownowa˙znych mamy zatem, ˙ze je´ sli uk lad (U
2) powstaje z uk ladu (U
1) przez wykonanie operacji elementarnej, to uk lady (U
1) i (U
2) s a r´
,ownowa˙zne. Zatem z twierdze´ n 2.5, 2.9 i 2.10 przez prost a indukcj
,e otrzymujemy st
,ad od razu nast
,epuj
,ace
,Twierdzenie 2.11. Za l´ o˙zmy, ˙ze uk lad (U
0) r´ owna´ n liniowych powstaje z uk ladu (U ) przez kolejne wykonanie sko´ nczonej liczby operacji elementarnych. W´ owczas uk lady te s a r´
,ownowa˙zne, maj a te same zbiory rozwi
,aza´
,n i uk lad (U ) jest sprzeczny wtedy i tylko wtedy, gdy uk lad (U
0) jest sprzeczny.
Na twierdzeniu 2.11 opiera si e metoda rozwi
,azywania uk lad´
,ow r´ owna´ n liniowych zwana me- tod a eliminacji Gaussa.
,3 Metoda eliminacji Gaussa
Metoda eliminacji Gaussa polega na znalezieniu dla danego uk ladu (1) (w kt´ orym a
ij6= 0 dla pewnych i, j) przy pomocy operacji elementarnych r´ ownowa˙znego mu (czyli posiadaj acego
,taki sam zbi´ or rozwi aza´
,n) uk ladu (9), kt´ ory po ewentualnej permutacji niewiadomych x
1, . . . , x
nma posta´ c:
x
1+ c
1, k+1x
k+1+ . . . + c
1nx
n= d
1x
2+ c
2, k+1x
k+1+ . . . + c
2nx
n= d
2x
3+ c
3, k+1x
k+1+ . . . + c
3nx
n= d
3. .. .. . .. . .. .
x
k+ c
k, k+1x
k+1+ . . . + c
knx
n= d
k0 = d
k+1. (9)
Je˙zeli d
k+16= 0, to uk lad (9) nie ma rozwi azania, a wi
,ec te˙z uk lad (1) nie ma rozwi
,azania i jest
,sprzeczny.
Je˙zeli d
k+1= 0 i k = n, to uk lad (1) posiada dok ladnie jedno rozwi azanie:
,x
1= d
1, x
2= d
2, . . . , x
n= d
n. (10)
Je˙zeli d
k+1= 0 oraz k < n, to x
k+1, x
k+2, . . . , x
ns a dowolnymi skalarami (nazywamy je
,parametrami), za´ s pozosta le niewiadome wyliczamy z r´ owna´ n uk ladu (9), tzn.
x
i= d
i− c
i k+1x
k+1− . . . − c
inx
ndla i = 1, 2, . . . , k. (11) Aby sprowadzi´ c uk lad (1) do postaci (9) nale˙zy najpierw przy pomocy operacji elementarnych przekszta lci´ c go do uk ladu postaci:
x
1+ a
012x
2+ . . . + a
01nx
n= b
01a
021x
1+ a
022x
2+ . . . + a
02nx
n= b
02.. . .. . . .. .. . a
0m1x
1+ a
0m2x
2+ . . . + a
0mnx
n= b
0m. (12)
Robimy to np. w ten spos´ ob, ˙ze najpierw znajdujemy element a
ij6= 0, a nast epnie przez
,operacje: x
1↔ x
j, r
1↔ r
i,
a1ij
· r
1doprowadzamy uk lad (1) do postaci (12).
Nast epnie przy pomocy r´
,ownania pierwszego eliminujemy zmienn a x
, 1z pozosta lych r´ owna´ n uk ladu (12) przez wykonanie operacji: r
2− a
021· r
1, r
3− a
031· r
1,...,r
m− a
0m1· r
1. Otrzymamy w´ owczas uk lad postaci:
x
1+ a”
12x
2+ . . . + a”
1nx
n= b”
1a”
22x
2+ . . . + a”
2nx
n= b”
2.. . . .. .. . a”
m2x
2+ . . . + a”
mnx
n= b”
m. (13)
Z kolei stosujemy nasz algorytm do uk ladu:
a”
12x
2+ . . . + a”
1nx
n= b”
1a”
22x
2+ . . . + a”
2nx
n= b”
2.. . . .. .. . a”
m2x
2+ . . . + a”
mnx
n= b”
m(14)
nie ruszaj ac pierwszego r´
,ownania uk ladu (13). Po sko´ nczonej liczbie krok´ ow uzyskamy uk lad postaci:
x
1+ e
12x
2+ e
13x
3+ . . . + e
1kx
k+ e
1 k+1x
k+1+ . . . + e
1nx
n= f
1x
2+ e
23x
3+ . . . + e
2kx
k+ e
2 k+1x
k+1+ . . . + e
2nx
n= f
2x
3+ . . . + e
3kx
k+ e
3 k+1x
k+1+ . . . + e
3nx
n= f
3. .. .. .
x
k+ e
k k+1x
k+1+ . . . + e
knx
n= f
k0 = f
k+1.
Je˙zeli f
k+16= 0, to otrzymany uk lad jest sprzeczny, a wi ec te˙z uk lad (1) jest sprzeczny.
,Je˙zeli za´ s f
k+1= 0, to przy pomocy operacji: r
1− e
1k· r
k, r
2− e
2k· r
k,...,r
k−1− e
k−1 k· r
keliminujemy zmienn a x
, kz pocz atkowych k − 1 r´
,owna´ n. P´ o´ zniej eliminujemy zmienn a x
, k−1z wcze´ sniejszych r´ owna´ n przy pomocy k − 1-szego r´ ownania, itd. W ko´ ncu, po sko´ nczonej liczbie krok´ ow, uzyskamy w ten spos´ ob uk lad (9).
Om´ owiony wy˙zej spos´ ob rozwi azywania uk ladu r´
,owna´ n metod a Gaussa zawiera du˙zo ele-
,ment´ ow dowolnych. Dowolno´ s´ c zachodzi na ka˙zdym etapie rozwa˙za´ n, poniewa˙z mo˙zemy eli- minowa´ c dowoln a niewiadom
,a (pod warunkiem, ˙ze odpowiedni wsp´
,o lczynnik nie r´ owna si e 0).
,Opr´ ocz tego dowolna jest r´ ownie˙z kolejno´ s´ c r´ owna´ n w danym uk ladzie. Je˙zeli np. w jakikolwiek spos´ ob zmienimy kolejno´ s´ c r´ owna´ n w wyj´ sciowym uk ladzie, to proces stopniowego eliminowa- nia niewiadomych przebiega´ c b edzie inaczej. Jednak zawsze musimy otrzyma´
,c t e sam
,a liczb
,e
,parametr´ ow!
W praktyce proces rozwi azywania uk ladu (1) mo˙zemy znacznie upro´
,sci´ c, je˙zeli zamiast prze- kszta lce´ n uk ladu r´ owna´ n b edziemy przekszta lca´
,c jego macierz uzupe lnion a A
, u. Oczywiste jest,
˙ze ka˙zdej operacji elementarnej uk ladu (1) odpowiada odpowiednia operacja elementarna ma- cierzy A
u, a mianowicie:
operacji r
i↔ r
jodpowiada operacja w
i↔ w
j, operacji a · r
iodpowiada operacja a · w
i,
operacji r
i+ a · r
jodpowiada operacja w
i+ a · w
j,
wykre´ slaniu i-tego r´ ownania odpowiada wykre´ slanie i-tego wiersza,
operacji x
i↔ x
jodpowiada operacja k
i↔ k
j(nale˙zy przy tym pami eta´
,c, ˙ze nie wolno rusza´ c ostatniej kolumny i na koniec nale˙zy jeszcze uwzgl edni´
,c wszystkie przenumerowania niewiado- mych!).
Przyk lad 2.12. Stosuj ac metod
,e eliminacji Gaussa rozwi
,a˙zemy nad cia lem R uk lad r´owna´
,n:
x
1− x
2− 9x
3+ 6x
4+ 7x
5+ 10x
6= 3
− 6x
3+ 4x
4+ 2x
5+ 3x
6= 2
− 3x
3+ 2x
4− 11x
5− 15x
6= 1 .
B edziemy wykonywali rachunki na macierzy uzupe lnionej naszego uk ladu:
,
1 1 −9 6 7 10 3
0 0 −6 4 2 3 2
0 0 −3 2 −11 −15 1
w1−3w3, w2−2w3
≡
1 1 0 0 40 55 0
0 0 0 0 24 33 0
0 0 −3 2 −11 −15 1
x2↔x4
≡
x1
1
x0
4 x0
3 x1
2x5
40 55
x60
0 0 0 0 24 33 0
0 2 −3 0 −11 −15 1
w2↔w3
≡
x1
1
x0
4 x0
3 x1
2x5
40 55
x60 0 2 −3 0 −11 −15 1
0 0 0 0 24 33 0
x3↔x5
≡
x1
1
x0
4x5
40
x1
2 x0
3x6
55 0 0 2 −11 0 −3 −15 1
0 0 24 0 0 33 0
1
2w2,241w3
≡
x1
1
x0
4x5
40
x1
2 x0
3x6
55 0 0 1 −
1120 −
32−
152 120 0 1 0 0
1180
w1−40w3, w2+112w3
≡
x1
1
x0
4 x0
5 x1
2 x0
3 x0
60 0 1 0 0 −
32 161 120 0 1 0 0
1180
. Zatem zmiennymi bazowymi s a x
, 2, x
3, x
6oraz x
1= −x
2, x
4=
12+
32x
3−
161x
6, x
5= −
118x
6. St ad uk lad posiada niesko´
,nczenie wiele
rozwi aza´
,n danych wzorami:
x
1= −a, x
2= a, x
3= b, x
4=
12+
32b −
161c, x
5= −
118c, x
6= c, gdzie a, b, c s a dowolnymi
,liczbami rzeczywistymi. 2
Przyk lad 2.13. Stosuj ac metod
,e eliminacji Gaussa rozwi
,a˙zemy uk lad r´
,owna´ n:
2x
1+ x
2− x
3+ x
4= 1 3x
1− 2x
2+ 2x
3− 3x
4= 2 5x
1+ x
2− x
3+ 2x
4= −1 2x
1− x
2+ x
3− 3x
4= 4
.
Rachunki b edziemy wykonywali na macierzy uzupe lnionej A
, unaszego uk ladu. Mamy, ˙ze
A
u=
2 1 −1 1 1
3 −2 2 −3 2
5 1 −1 2 −1
2 −1 1 −3 4
x1↔x2
≡
x2
1
x2
1x3
−1
x1
41
−2 3 2 −3 2
1 5 −1 2 −1
−1 2 1 −3 4
w2+2w1, w3−w1, w4+w1
≡
x2
1
x2
1x3
−1
x1
41
0 7 0 −1 4
0 3 0 1 −2
0 4 0 −2 5
x1↔x4
≡
x2
1
x1
4x3
−1
x2
11
0 −1 0 7 4
0 1 0 3 −2
0 −2 0 4 5
w3+w2, w4−2w2
≡
x2
1
x1
4x3
−1
x2
11
0 −1 0 7 4
0 0 0 −10 2
0 0 0 10 −3
w4+w3
≡
x2
1
x1
4x3
−1
x2
11
0 −1 0 7 4
0 0 0 10 2
0 0 0 0 −1
. Zatem nasz uk lad jest sprzeczny
i nie posiada rozwi azania, bo ostatnie r´
,ownanie ma posta´ c:
0 · x
2+ 0 · x
4+ 0 · x
3+ 0 · x
1= −1. 2
Przyk lad 2.14. Stosuj ac metod
,e eliminacji Gaussa rozwi
,a˙zemy uk lad r´
,owna´ n:
2x
1− x
2+ x
3− x
4= 1
2x
1− x
2− 3x
4= 2
3x
1− x
3+ x
4= −3
2x
1+ 2x
2− 2x
3+ 5x
4= −6 .
Rachunki b edziemy wykonywali na macierzy uzupe lnionej A
, unaszego uk ladu. Mamy, ˙ze
A
u=
2 −1 1 −1 1
2 −1 0 −3 2
3 0 −1 1 −3
2 2 −2 5 −6
x1↔x3
≡
x3
1
x2
−1
x2
1x4
−1 1
0 −1 2 −3 2
−1 0 3 1 −3
−2 2 2 5 −6
w3+w1, w4+2w1
≡
x3
1
x2
−1
x2
1x4
−1 1
0 −1 2 −3 2
0 −1 5 0 −2
0 0 6 3 −4
(−1)w2
≡
x3
1
x2
−1
x2
1x4
−1 1
0 1 −2 3 −2
0 −1 5 0 −2
0 0 6 3 −4
.
Wykonujemy operacj e w
, 3+ w
2:
x3
1
x2
−1
x2
1x4
−1 1
0 1 −2 3 −2
0 0 3 3 −4
0 0 6 3 −4
w4−2w3
≡
x3
1
x2
−1
x2
1x4
−1 1
0 1 −2 3 −2
0 0 3 3 −4
0 0 0 −3 4
.
Wykonujemy operacje
13w
3i (−
13)w
4:
x3
1
x2
−1
x2
1x4
−1 1
0 1 −2 3 −2
0 0 1 1 −
430 0 0 1 −
43
w1+w4, w2−3w4, w3−w4
≡
x3
1
x2
−1
x2
1 x0
4−
130 1 −2 0 2
0 0 1 0 0
0 0 0 1 −
43
.
Wykonujemy operacje w
1− 2w
3i w
2+ 2w
3:
x3
1
x2
−1
x0
1 x0
4−
130 1 0 0 2
0 0 1 0 0
0 0 0 1 −
43
w1+w2
≡
x3