Ortogonalizacja Grama-Schmidta

Z Wikipedii

Masz nowe wiadomości (różnica z poprzednią wersją).

Ortogonalizacja Grama-Schmidta to metoda za pomocą której można przekształcić zbiór liniowo niezależnych wektorów przestrzeni unitarnej w zbiór wektorów ortogonalnych. Przestrzenie liniowe, rozpinane przez zbiory przed i po ortogonalizacji są tożsame, tak więc proces może służyć do ortogonalizowania bazy.

Proces został nazwany na cześć Jørgena Grama, matematyka duńskiego, oraz Erharda Schmidta, matematyka niemieckiego.

[edytuj] Proces ortogonalizacji

Operator rzutowania ortogonalnego wektora $v$ na wektor $u$ definiujemy jako:

$\mathrm{proj}_{\mathbf{u}}\,\mathbf{v} = {\langle \mathbf{v}, \mathbf{u}\rangle\over\langle \mathbf{u}, \mathbf{u}\rangle}\mathbf{u}.$

Wówczas dla układu k wektorów{v₁, …, v_k} proces przebiega następująco:

$\mathbf{u}_1 = \mathbf{v}_1,$

Dwa pierwsze kroki procesu ortogonalizacji

$\mathbf{u}_2 = \mathbf{v}_2-\mathrm{proj}_{\mathbf{u}_1}\,\mathbf{v}_2,$

$\mathbf{u}_3 = \mathbf{v}_3-\mathrm{proj}_{\mathbf{u}_1}\,\mathbf{v}_3-\mathrm{proj}_{\mathbf{u}_2}\,\mathbf{v}_3,$

$\vdots$

$\mathbf{u}_k = \mathbf{v}_k-\sum_{j=1}^{k-1}\mathrm{proj}_{\mathbf{u}_j}\,\mathbf{v}_k,$

Otrzymany zbiór {u₁, …, u_k} jest zbiorem wektorów ortogonalnych.

Aby zbudować w ten sposób zbiór ortonormalny, każdy wektor należy podzielić przez jego normę:

$\mathbf{e}_n = {\mathbf{u}_n\over||\mathbf{u}_n||}, n=1, 2, ..., k$

Dowód ortogonalności tak otrzymanego układu opiera się na indukcji.

Proces ortogonalizacji pozwala na wskazanie bazy ortogonalnej w dowolnej n-wymiarowej przestrzeni unitarnej. Jej istnienie można wykazać na bazie lematu Kuratowskiego-Zorna.

[edytuj] Funkcje ciągłe

Jeżeli iloczyn skalarny funkcji ciągłych jest określony wzorem:

$\langle f,g\rangle _w = \int\limits_a^b w(x) f(x) g(x) dx.$

gdzie $w (x)$ jest funkcją wagową, to dla zbioru funkcji liniowo niezależnych przekształcenie w zbiór funkcji ortogonalnych przebiega następująco:

g 0 (x) = f 0 (x)

$g_i(x) = f_i(x) - \sum_{j=0}^{i-1} g_j(t)\frac{\int\limits_a^b w(t) f_i(t)g_j(t) dt}{\int\limits_a^b w(t) g_j^2(t) dt}$

Iloczyn skalarny funkcji $g i (x)$ i $g j (x)$ dla różnych $i$ , $j$ wynosi (bez straty ogólności przyjmijmy, że $i > j$ ):

$\int\limits_a^b w(t) g_j(t) g_i(t) dt = \int\limits_a^b w(t) g_j(t) f_i(t) dt - \int\limits_a^b w(s) g_j(s) g_j(s)\frac{\int\limits_a^b w(t) f_i(t)g_j(t) dt}{\int\limits_a^b w(t) g_j^2(t) dt} ds +$

$- \sum_{k=0 \and k\ne j}^{i-1} \frac{\int\limits_a^b w(t) f_i(t)g_k(t) dt}{\int\limits_a^b w(t) g_k^2(t) dt} \int\limits_a^b w(s) g_j(s) g_k(s) ds$

Jeśli dla wszystkich różnych par $j, k$ mniejszych od $i$ iloczyn skalarny wynosi 0, to:

$\int\limits_a^b w(t) g_j(t) g_i(t) dt = \int\limits_a^b w(t) g_j(t) f_i(t) dt - \int\limits_a^b w(s) g_j(s) g_j(s)\frac{\int\limits_a^b w(t) f_i(t)g_k(t) dt}{\int\limits_a^b w(t) g_k^2(t) dt} ds$

$\int\limits_a^b w(t) g_j(t) g_i(t) dt = \int\limits_a^b w(t) f_i(t) g_j(t) dt - \frac{\int\limits_a^b w(t) f_i(t)g_j(t) dt}{\int\limits_a^b w(t) g_j^2(t) dt} \int\limits_a^b w(s) g_j^2(s) ds$

$\int\limits_a^b w(t) g_j(t) g_i(t) dt = \int\limits_a^b w(t) f_i(t) g_j(t) dt - \int\limits_a^b w(t) f_i(t)g_j(t) dt = 0$

W wielu przypadkach wynikowy układ $g_0,g_1,\cdots ,g_n$ obliczony numerycznie może nie być ortogonalny. Jednym ze sposób poprawienia wyniku jest reortogonizacja - zastosowanie algorytmu ortogonizacyjnego dla elementów $g i$ zamiast $f i$ . Prostszym sposobem jest jednak zastosowanie zmodyfikowanego algorytmu Grama-Schmidta:

   for(i=0;i<=n;i++)
   {
     g[i]=f[i]; 
     for(j=i+1;j<=n;j++)
        f[j]=f[j] - <math>\frac{(g[i],f[i])}{(g[i],g[i])}\cdot g[i]</math>
   }

Gdzie:

$(g[i],f[i])=\sum_{j=1}^{n}f(x_i)\cdot g(x_i) \cdot w(x_i)$

[edytuj] Zobacz też

przegląd zagadnień z zakresu matematyki

Źródło: "http://pl.wikipedia.org../../../o/r/t/Ortogonalizacja_Grama-Schmidta_10a3.html"

Kategoria: Algebra liniowa

See also ebooksgratis.com: no banners, no cookies, totally FREE.

Ortogonalizacja Grama-Schmidta

Z Wikipedii

[edytuj] Proces ortogonalizacji

[edytuj] Funkcje ciągłe

[edytuj] Zobacz też

Views

nawigacja

techniczne

zmiany

Szukaj

W innych językach