Peatüki algus: LU-lahutus
Eelmine: Kolmnurksete võrrandisüsteemide lahendamine
Järgmine: QR-lahutus


Gaussi teisendus ja LU-lahutus

Teatud tingimustel on võrrandisüsteemi tex2html_wrap_inline709 süsteemimaatriks A esitatav alumise ühikdiagonaaliga, st peadiagonaali elementideks on ühed, kolmnurkmaatriksi L ja ülemise kolmnurkmaatriksi U korrutisena ja lahendada tuleb kaks kolmnurkse maatriksiga võrrandisüsteemi.

Lause 1.2.1 (LU-meetod). Kui tex2html_wrap_inline719 kus L on alumine ühikdiagonaaliga kolmnurkmaatriks ja U on ülemine regulaarne kolmnurkmaatriks, ning tex2html_wrap_inline725 siis tex2html_wrap_inline727 ning süsteemi lahendamiseks tuleb alguses lahendada süsteem tex2html_wrap_inline729 ja seejärel lahendada süsteem tex2html_wrap_inline731

Näide 1.2.1. Lahendame süsteemi
displaymath665
LU-meetodil. Kuna
displaymath666
siis lause 1.2.1 põhjal tuleb, esiteks, lahendada süsteem
displaymath667
Selle süsteemi lahendiks on tex2html_wrap_inline735 ja tex2html_wrap_inline737 Teiseks, lahendades süsteemi
displaymath668
leiame, et tex2html_wrap_inline739 ja tex2html_wrap_inline741 Seega, tex2html_wrap_inline743

Lineaaralgebra põhikursuses võrrandisüsteemide lahendamiseks esitatud Gaussi võte on rakendatav ka LU-lahutuse korral. Olgu tex2html_wrap_inline747 kusjuures tex2html_wrap_inline749 Kui
displaymath669
ja
 equation274
siis

displaymath670

Definitsioon 1.2.1. Maatriksit Mk kujul (2) nimetatakse Gaussi maatriksiks, komponente tex2html_wrap_inline753 nimetatakse Gaussi kordajateks ja vektorit tex2html_wrap_inline755 Gaussi vektoriks ning Gaussi maatriksiga tex2html_wrap_inline757määratud teisendust nimetatakse Gaussi teisenduseks.

Definitsioon 1.2.2. Maatriksi tex2html_wrap_inline759 k-ndaks juhtelemendiks nimetatakse suurust
displaymath671
kusjuures tex2html_wrap_inline763 ja tex2html_wrap_inline765

Kui tex2html_wrap_inline767 siis maatriksi A nullist erinevate juhtelementide korral on leitavad Gaussi maatriksid tex2html_wrap_inline771 nii, et tex2html_wrap_inline773 on ülemine kolmnurkmaatriks.

Näide 1.2.2. Vaatleme Gaussi maatriksite M1 ja tex2html_wrap_inline777ning ülemise kolmnurkmaatriksi U leidmist maatriksi
displaymath672
korral. Seose (2) abil leiame, et
displaymath673

displaymath674
Seega,
displaymath675
ja
displaymath676

displaymath677
ning

displaymath678

Märgime, et maatriks tex2html_wrap_inline781 on ülemine kolmnurkne veergudes ühest kuni (k-1)-ni ja Gaussi maatriksi Mk elementide arvutamiseks kasutatakse maatriksi vektorit A(k-1)(k:m,k), kusjuures maatriksi Mk leidmine on võimalik, kui tex2html_wrap_inline791. Lisaks, tex2html_wrap_inline793 Kui valida
displaymath679
siis
displaymath680
Rõhutame, et meie käsitluses alumine kolmnurkmaatriks L on ühikdiagonaaliga maatriks.

Lause 1.2.2. Kui peamiinorid, so tex2html_wrap_inline7971: k , 1: tex2html_wrap_inline8011: n-1), siis maatriksil tex2html_wrap_inline807 on olemas LU-lahutus. Kui regulaarsel maatriksil A eksisteerib LU-lahutus, siis see on ühene ja tex2html_wrap_inline815

Tõestus. Oletame, et k-1 sammu on sooritatud ja leitud on maatriks tex2html_wrap_inline819 Element akk(k-1) on maatriksi A k-s juhtelement ja tex2html_wrap_inline7971:k, 1:tex2html_wrap_inline831 Seega, kui A(1:k, 1:k) on regulaarne, siis tex2html_wrap_inline839 ja eksisteerib maatriksi A jaoks LU-lahutus. Oletades, et regulaarsel maatriksil A on kaks LU-lahutust A=L1U1 ja A=L2U2, saame, et L1U1=L2U2 ehk L2-1L1=U2U1-1. Kuna L2-1L1 o n alumine kolmnurkmaatriks, mille peadiagonaalil on ühed ja U2U1-1on ülemine kolmnukmaatriks, siis L2-1L1=I ja U2U1-1=I ning L2=L1 ja tex2html_wrap_inline867

Näide 1.2.3.* Leiame maatriksi
displaymath301
LU-lahutuse.
Leiame Gaussi maatriksi M1 maatriksi A korral:
displaymath302

displaymath303
Seega
displaymath304
ja

displaymath305

Näide 1.2.4.* Leiame maatriksi
displaymath306
LU-lahutuse.
Leiame Gaussi maatriksi M1 maatriksi A korral:
displaymath307

displaymath308
Seega
displaymath309
ja kuna maatriks M1A on ülemine kolmnurkmaatriks, siis M2=I ning

displaymath310

displaymath311

Näide 1.2.5.* Lahendame süsteemi tex2html_wrap_inline381, kus
displaymath312
kasutades LU-lahutust.
Näites 1.2.4 on leitud maatriksi A jaoks LU -lahutus:
displaymath311
Süsteemi tex2html_wrap_inline389 so
displaymath314
lahendamisel saame, et
displaymath315
Süsteemi tex2html_wrap_inline391 so
displaymath316
lahendamisel leiame, et

displaymath317

Ülesanne 1.2.1.* Leidke maatriksi A jaoks LU-lahutus, kui
displaymath318

Ülesanne 1.2.2.* Lahendage süsteem tex2html_wrap_inline381, kus
displaymath319
kasutades LU-lahutust.

Ka ristkülikmaatriksi tex2html_wrap_inline759, mille peamiinorid on nullist erinevad, so
displaymath681
jaoks eksisteerib LU-lahutus.

Näide 1.2.6. Järgnevad seosed peavad paika
displaymath682

displaymath683
Algebra põhikursusest on teada, et Gaussi elimineerimismeetodi vahetu rakendamine , seega ka LU-lahutuse vahetu teostamine, nurjub, kui vähemalt üks maatriksi peamiinoreist on singulaarne. Osutub, et regulaarse maatriksi korral on peale maatriksi ridade järjekorra sobivat muutmist võimalik teisendatud maatriksi LU-lahutus. Ridade (ka veergude) järjekorra muutmiseks kasutatakse permutatsioonimaatrikseid.

Definitsioon 1.2.3. Permutatsioonimaatriksiks tex2html_wrap_inline877 nimetatakse maatriksit, mis on saadud ühikmaatriksist I ridade järjekorra muutmise teel.

Näide 1.2.7 . Vaatleme, kuidas mõjub tex2html_wrap_inline881 maatriksile A selle korrutamine ühe konkreetse permutatsioonimaatriksiga P.
displaymath684
Permutatsioonimaatriksiga P vasakult korrutamisel saame uue maatriksi, mis on lähtemaatriksist saadud ridade samal viisil ümberpaigutamisel, kui on saadud maatriks P maatriksist I. Korrutades paremalt,
displaymath685
saame uue maatriksi, mis on lähtemaatriksist saadud veergudede samal viisil ümberpaigutamisel, kui on saadud maatriks P maatriksist I veergude ümberpaigutamise teel. Peab paika

Lause 1.2.3. Kui tex2html_wrap_inline897 ja tex2html_wrap_inline899 siis leidub selline permutatsioonimaatriks tex2html_wrap_inline901 et maatriksi PA kõik peamiinorid on nullist erinevad ja seega eksisteerib LU-lahutus

Näide 1.2.8.* Olgu
displaymath320
Leiame sobiva permutatsioonimaatriksi tex2html_wrap_inline407 korral maatriksi PA jaoks LU-lahutuse.
Vahetame ära maatriksi A esimese ja teise rea, st võtame
displaymath321
ning leiame maatriksi
displaymath322
korral Gaussi maatriksi
displaymath323
Seega
displaymath324
ja
displaymath325
ning
displaymath326
Järelikult,
displaymath327
ja

displaymath328

Ülesanne 1.2.3.* Leiame sobiva permutatsioonimaatriksi P korral maatriksi PA jaoks LU-lahutuse, kui
displaymath329
Peatüki algus: LU-lahutus
Eelmine: Kolmnurksete võrrandisüsteemide lahendamine
Järgmine: QR-lahutus