Peatüki algus: Singulaarlahutus
Eelmine: Singulaarlahutuse omadused
Järgmine: Pseudo-pöördmaatriks


Singulaarlahutuse algoritm

Algoritm 3.3.1. Maatriksi tex2html_wrap_inline841 singulaarlahutuse (2) leidmiseks tuleb:

I Leida maatriksi ATA omaväärtused ja paigutada need kahanemise järjekorras.
II Leida maatriksi ATA nullist erinevate omaväärtuste arv r.
III Leida saadud omaväärtustele vastavad maatriksi ATA ortonormeeritud omavektorid ja paigutada need samas järjekorras maatriksi tex2html_wrap_inline927 veeruvektoreiks.
IV Moodustada diagonaalmaatriks tex2html_wrap_inline1001 mille peadiagonaalile paigutada ruutjuured tex2html_wrap_inline1003 punktis I saadud maatriksi tex2html_wrap_inline1005 tex2html_wrap_inline837 esimesest omaväärtusest kahanevas järjestuses.
V Leida maatriksi tex2html_wrap_inline1009esimesed veeruvektorid:
 
equation427
VI Lisada maatriksisse U ülejäänud m-r veeruvektorit Gram-Schmidt'i ortogonaliseerimisprotsessi abil. tex2html_wrap_inline765

Näide 3.3.1. Leiame maatriksi
displaymath957
singulaarlahutuse.

I Leiame tex2html_wrap_inline1017 omaväärtused: tex2html_wrap_inline1019
II Leiame maatriksi ATA nullist erinevate omaväärtuste arvu: r=2 .
III Leiame maatriksi ATA omaväärtustele tex2html_wrap_inline1029 ja tex2html_wrap_inline1031 v astavad ortonormeeritud omavektorid tex2html_wrap_inline1033 ja tex2html_wrap_inline1035 ning moodustame maatriksi
displaymath958
IV Leiame singulaarmaatriksi tex2html_wrap_inline1037
displaymath959
mille peadiagonaalil on ruutjuured maatriksi ATA omaväärtustest (kahanevas järjestuses) ja ülejäänud maatriksi tex2html_wrap_inline833 elemendid on nullid.
V Leiame valemite (9) abil maatriksi tex2html_wrap_inline1043 kaks esimest veeruvektorit,
displaymath960
ja

displaymath961
VI Vektori tex2html_wrap_inline1045 arvutamiseks leiame Gram-Schmidt'i protsessiga algul vektori tex2html_wrap_inline1047, mis on risti vektoritega tex2html_wrap_inline1049 ja tex2html_wrap_inline1051:
displaymath962
Normeerides vektori tex2html_wrap_inline1053 saame
displaymath963
Seega
displaymath964
ja maatriksi A singulaarlahutuseks on:

displaymath965

Näide 3.3.2. Leiame maatriksi tex2html_wrap_inline1057 singulaarlahutuse.

I Leiame maatriksi ATA omaväärtused:
displaymath966
II Leiame maatriksi ATA nullist erinevate omaväärtuste arvu: r=1.
III Leiame maatriksi ATA omavektorid:
displaymath967

displaymath968
Kuna omaväärtus 0 on kordne, siis vektori tex2html_wrap_inline1067 leidmiseks on kasutatud Gram-Schmidt'i ortogonaliseerimisprotsessi. ''Kleebime'' kokku ortogonaalmaatriksi V:

displaymath969
IV Koostame singulaarmaatriksi:
displaymath970
V Arvutame valemite (9) abil maatriksi U ainsa veeruvektori:
displaymath971
Seega, maatriksi A singulaarlahutus on

displaymath972

Näide 3.3.3.* Leiame maatriksi
displaymath197
singulaarlahutuse.
Antud tex2html_wrap_inline237-maatriksil A on ülimalt kolm nullist erinevat singulaarväärtust. Järelikult on otstarbekas leida maatriksi A nullist erinevad singulaarväärtused tex2html_wrap_inline243-maatriksi AA T abil (mitte tex2html_wrap_inline247-maatriksi ATA abil). Kuna
displaymath198
siis maatriksi AAT karakteristlik võrrand on
displaymath199
ehk
displaymath200
ja selle võrrandi lahendeiks on tex2html_wrap_inline253 ning tex2html_wrap_inline255 Kuna tex2html_wrap_inline257 ja maatriksi tex2html_wrap_inline259 on tex2html_wrap_inline237-maatriks, siis maatriksi tex2html_wrap_inline259 peadiagonaalil on maatriksi A singulaarväärtused kahanevas järjekorras ning ülejäänud maatriksi tex2html_wrap_inline259 elemendid on n
displaymath201
Maatriksi U veeruvektoreiks on maatriksi AAT ortonormeeritud omavektorid:
displaymath202

displaymath203

displaymath204
Vektorite tex2html_wrap_inline273 ja tex2html_wrap_inline275 ''kleepimisel'' saame maatriksi
displaymath205
Vastavalt seosele (6) leiame maatriksi V kolm esimest veeruvektorit (maatriksi tex2html_wrap_inline259 peadiagonaalil on kolm nullist erinevat elementi) valemi
displaymath206
abil. Seega
displaymath207
Vektori tex2html_wrap_inline281 arvutamiseks leiame Gram-Schmitd'i ortogonaliseerimisprotsessiga algul vektori tex2html_wrap_inline283, mis on risti vektoritega tex2html_wrap_inline285 tex2html_wrap_inline287 ja tex2html_wrap_inline289:
displaymath208

displaymath209
Kuna tex2html_wrap_inline291 siis
displaymath210
ja
displaymath211
Teostame kontrolli
displaymath212

displaymath213
ja

displaymath214

displaymath215

Ülesanne 3.3.1.* Kasutades näites 3.3.3 saadud maatriksi A singulaarlahutust leida maatriksi A veeruvektorite alamruumi tex2html_wrap_inline299 parempoolse nullruumi tex2html_wrap_inline301 reavektorite alamruumi tex2html_wrap_inline303 ja vasakpoolse nullruumi tex2html_wrap_inline305 baasid.

Ülesanne 3.3.2.* Leida maatriksi tex2html_wrap_inline309 singulaarlahutus ja QR-lahutus.

Ülesanne 3.3.3.* Leida maatriksi tex2html_wrap_inline315 singulaarlahutus.