Peatüki algus: Lineaarsete algebraliste võrrandisüsteemide lahendamise otsesed meetodid
Eelmine: Maatriksi LDMT-lahutus ja LDLT-lahutus
Järgmine: Positiivselt poolmääratud maatriksid


Positiivselt määratud süsteemid

Definitsioon 6.2.1. Maatriksit tex2html_wrap_inline1495 nimetatakse positiivselt määratuks, kui
displaymath1555
iga nullist erineva vektori tex2html_wrap_inline1631 korral.

Näide 6.2.1. Maatriks
displaymath1556
on positiivselt määratud, sest tex2html_wrap_inline1633 korral

displaymath1557

displaymath1558

Ülesanne 6.2.1.* Näidake, et maatriks
displaymath459
on positiivselt määratud.

Lause 6.2.1. Kui tex2html_wrap_inline1495 on positiivselt määratud maatriks ja maatriksi tex2html_wrap_inline1637 veeruvektorid on lineaarselt sõltumatud, siis maatriks
displaymath1559
on samuti positiivselt määratud.

Tõestus. Kui vektori tex2html_wrap_inline1639 korral peab paika seos
displaymath1560
siis
displaymath1561
ja maatriksi A positiivsest määratusest järeldub, et tex2html_wrap_inline1643 Kuna maatriksi X veeruvektorid on lineaarselt sõltumatud, siis tingimusest tex2html_wrap_inline1647 järeldub, et tex2html_wrap_inline1649 Seega, tingimustest tex2html_wrap_inline1639 ja tex2html_wrap_inline1653 järeldub, et tex2html_wrap_inline1655 st maatriks B on positiivselt määratud. tex2html_wrap_inline1659

Järeldus 6.2.1. Kui maatriks tex2html_wrap_inline1495 on positiivselt määratud, siis kõik maatriksi A alammaatriksid, mis on saadud maatriksi A samanumbriliste ridade ja veergude kustutamisel, on positiivselt määratud ja kõik maatriksi A peadiagonaali elemendid on positiivsed.

Tõestus. Kui tex2html_wrap_inline1669 on selliste naturaalarvuliste koordinaatidega vektor, et
displaymath1562
siis
displaymath1563
on maatriks, mis on saadud ühikmaatriksi In veergudest, indeksitega tex2html_wrap_inline1673 Seega, maatriksi X veeruvektorid on lineaarselt sõltumatud ja lause 6.2.1 põhjal on maatriks XTAX positiivselt määratud. Maatriks XTAX on maatriksi A alammaatriks, mis on saadud maatriksi A ridadest ja veergudest numbritega tex2html_wrap_inline1673 Seega, kõik maatriksi A alammaatriksid, mis on saadud maatriksi A samanumbriliste ridade j a veergude kustutamisel, on positiivselt määratud. Valides k=1, saame järelduse väite teise poole. tex2html_wrap_inline1519

Järeldus 6.2.2. Kui A tex2html_wrap_inline1697 on positiivselt määratud, siis maatriksil A leidub lahutus A=LDMT ja diagonaalmaatriksi tex2html_wrap_inline1703peadiagonaali kõik elemendid on positiivsed.

Tõestus. Järelduse 6.2.1 põhjal on kõik maatriksi A alammaatriksid A(1 : k, 1 : k) (1tex2html_wrap_inline1709ktex2html_wrap_inline1711 n) positiivselt määratud ja, seega, regulaarsed maatriksid ning lause 6.1.1 põhjal eksisteerib LDMT-lahutus. Võttes lauses 6.2.1 X=L-T, leiame, et maatriks
displaymath1564
on positiivselt määratud. Kuna maatriks MTL-T ülemine ühikdiagonaaliga kolmnurkmaatriks, siis maatriksitel B ja D on sama peadiagonaal ja seal peavad maatriksi B positiivse määratuse tõttu olema positiivsed elemendid. tex2html_wrap_inline1659

Lause 6.2.2 (Cholesky lahutus). Kui maatriks A tex2html_wrap_inline1697 on sümmeetriline ja positiivselt määratud, siis leidub täpselt üks alumine positiivse peadiagonaaliga kolmnurkmaatriks G, et
 
equation224

Tõestus. Lause 6.1.2 põhjal leiduvad ja on üheselt määratud ühikdigonaaliga alumine kolmnurkmaatriks L ja diagonaalmaatriks tex2html_wrap_inline1501 nii, et kehtib lahutus (3), so A=LDLT. Järelduse 6.2.2 kohaselt on maatriksi D elemendid dk positiivsed. Seetõttu maatriks
displaymath1565
on positiivse peadiagonaaliga alumine kolmnurkmaatriks ja paika peab seos (5). Lahutuse (5) ühesus järeldub lahutuse (3) ühesusest. tex2html_wrap_inline1519

Lahutus (5) on tuntud Cholesky lahutuse nime all. Maatriksit G nimetatakse maatriksi A Cholesky kolmnurkmaatriksiks. Sümmeetrilise ja positiivselt määratud maatriksiga A võrrandisüsteemi
displaymath1566
lahendamiseks tuleb, esiteks, leida maatriksi A Cholesky kolmnurkmaatriksiks G. Teiseks, tuleb lahendada kolmnurkmaatriksiga süsteem
displaymath1567
Kolmandaks, tuleb lahendada süsteem

displaymath1568

Cholesky lahutust on võimalik leida samm-sammult.

Lause 6.2.3. Kui maatriks A tex2html_wrap_inline1697 on sümmeetriline ja positiivselt määratud, siis tähistuse
displaymath1569
korral on maatriks A esitatav kujul
 equation269
kusjuures tex2html_wrap_inline1765 Maatriks tex2html_wrap_inline1767 on positiivselt määratud. Kui
displaymath1570
siis A=GGT, kusjuures

displaymath1571

Tõestus. Kontrollime lahutuse (6) õigsust:
displaymath1572

displaymath1573

displaymath1574
Kui
displaymath1575
siis
displaymath1576

displaymath1577
Kuna maatriks A on positiivselt määratud ja maatriksi X veeruvektorite süsteem on lineaarselt sõltumatu, siis lause 6.2.1 põhjal on positiivselt määratud ka maatriks
displaymath1578
ja järelduse 6.2.1 väitel on maatriks tex2html_wrap_inline1767 samuti positiivselt määratud ning võime analoogiliselt maatriksi A lahutusega blokkideks jaotada ka maatriks tex2html_wrap_inline1767 blokkideks, jne.

Näide 6.2.1. Leiame maatriksi
displaymath1579
LU-lahutuse, LDMT-lahutuse, LDLT-lahutuse ja Cholesky lahutuse.

Nendime, et kui maatriksi A peamiinorid on nullist erinevad, siis maatriksi A viimisel Gaussi teisenduse abil kolmnurkkujule, leiame samaaegselt nii maatriksi L kui ka maatriksi U. Nimelt, alumise kolmnurkmaatriksi L element lij tex2html_wrap_inline1799 võrdub teguriga, millega tuleb j-ndat rida korrutada mahalahutamisel i-ndast reast i-ndas reas oleva elemendi nullistamisel. Leiame, et
displaymath1580
ja
displaymath1581
ning
displaymath1582
Teades maatriksi A jaoks LU-lahutust, leiame maatriksi A jaoks LDMT-lahutuse, LDLT-lahutuse ja Cholesky lahutuse:
displaymath1583
ja
displaymath1584
Leiame maatriksi A Cholesky lahutuse ka samm-sammult, kasutades lauses 6.2.3 esitatud algoritmi. Kuna esimesel sammul
displaymath1585
siis
displaymath1586
Järgneval sammul
displaymath1587
ja
displaymath1588
Seega
displaymath1589
ja
displaymath1590
ning
displaymath1591

Ülesanne 6.2.2.* Leidke positiivselt määratud maatriksi
displaymath460
Cholesky lahutus.

Ülesanne 6.2.3.* Lahendage võrrandisüsteem tex2html_wrap_inline575 kus
displaymath461
kui on teada maatriksi A Cholesky lahutus
displaymath462


Peatüki algus: Lineaarsete algebraliste võrrandisüsteemide lahendamise otsesed meetodid
Eelmine: Maatriksi LDMT-lahutus ja LDLT-
Järgmine: Positiivselt poolmääratud maatriksid