Gramova–Schmidtova ortogonalizace: Porovnání verzí

Smazaný obsah Přidaný obsah
pridani prikladu s lauchliho matici
→‎Ztráta ortogonality: ztrata ortogonality obecne
Řádek 76:
=== Ztráta ortogonality ===
 
Nechť <math>\hat{Q}_n</math> (<math>\hat{Q}_n=Q_n+E_n</math>) je matice vektorů spočtených pomocí některé varianty Gramova-Schmidtova algoritmu v počítači se standardní konečnou aritmetikou s plovoucí řádovou čárkou, tj. <math>\hat{Q}_n=Q_n+E_n\approx Q_n</math> a <math>\hat{Q}_n^T\hat{Q}_n\approx I_n</math>. Veličina
 
:<math>\|\hat{Q}_n^T\hat{Q}_n-I_n\|_2</math>
Řádek 87:
 
kde <math>\kappa_2(A)</math> je [[Podmíněnost matice|podmíněnost]] matice <math>A</math>.
Uvažujeme-li standardní aritmetiku se [[Strojová přesnost|strojovou přesností]] <math>\epsilon_M\approx2.22\times10^{-16}</math> (double), pak ztráta ortogonality odpovídající jednotlivým výše zmíněným algoritmům aplikovaným na danou Lauchcliho matici, je vve druhém sloupci následující tabulce:tabulky. Ve třetím sloupci je obecný vztah platný pro libovolnou matici <math>A</math>
 
{| class="wikitable"
|+
! Algoritmus
! Ztráta ortogonality (Lauchliho matice)
! Ztráta ortogonality (obecně)
|-
|CGS
|<math>2.2\times10^{-2}</math>
|<math>\kappa_2^2(A)\epsilon</math>
|-
|MGS
|<math>2.2\times10^{-9}</math>
|<math>\kappa_2(A)\epsilon</math>
|-
|ICGS
|<math>2.4\times10^{-16}</math>
|<math>\epsilon</math>
|-
|}
 
== Vztah Gramova-Schmidtova algoritmu a QR rozkladu ==