QR rozklad

QR rozklad dané matice je způsob, jak zapsat tuto matici jako součin dvou matic, z nichž jedna je ortogonální (tj. její sloupce tvoří ortonormální systém) a druhá je v horním trojúhelníkovém tvaru. (Pozor, nezaměňovat QR rozklad s QR algoritmem, který slouží k výpočtu vlastních čísel čtvercové matice.)

Definice

Nechť , QR rozkladem nazýváme vztah

,

kde má vzájemně ortonormální sloupce (tj. ) a je v horním trojúhelníkovém tvaru (tj. pro všechna ).

Lineárně nezávislé sloupce

Pokud má matice lineárně nezávislé sloupce, pak

,

kde je unitární (v reálném případě ortogonální) matice, a je horní trojúhelníková regulární matice.

Označme , sloupce matic , platí

,

přičemž značí lineární obal. Tedy a obsahuje ortonormální bázi prostoru generovaného sloupci matice .

Pokud navíc volíme diagonální prvky matice kladné, je QR pak rozklad

jednoznačný. Je-li , tedy je-li regulární, pak a nulový blok v matici neexistují, , a tedy i QR rozklad lze volit jednoznačný.

Lineárně závislé sloupce

Pokud má rozkládaná matice lineárně závislé sloupce, QR rozklad zpravidla uvažujeme tak, aby i nadále platilo . Nechť , pak

,

kde oproti předchozímu případu a je v horním schodovitém tvaru (pokud je pak blok a nulový blok v matici neexistují).

Vždy existuje permutace sloupců matice realizovaná permutační maticí tak, že

,

kde je horní trojúhelníková regulární matice, kterou lze volit tak, že její diagonální prvky jsou kladné.


Výpočet QR rozkladu

QR rozklad lze provést pomocí klasického nebo modifikovaného Gramova-Schmidtova algoritmu (případně s iteračním zpřesněním), nebo pomocí Householderových nebo Givensových transformačních matic. Při reálném výpočtu (tj. v aritmetice s konečnou přesností) se všechny zmíněné postupy výrazně liší v přesnosti a rychlosti výpočtu. Přesnost je klíčovým faktorem zejména v případě, že matice obsahuje lineárně závislé sloupce.

LQ rozklad

LQ rozkladem matice nazveme transponovaný a komplexně sdružený (tzv. hermitovsky sdružený) QR rozklad matice . Tedy, je-li

,

kde je v dolním trojúhelníkovém tvaru, představuje LQ rozklad matice .

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.