Macierz Grama - Google

Macierz Grama

Z Wikipedii

(Przekierowano z Wyznacznik Grama)
Skocz do: nawigacji, szukaj

Macierz Grama dla zbioru wektorów v_1,\dots, v_n w przestrzeni unitarnejmacierz symetryczna, której wyrazy są iloczynami skalarnymi tych wektorów. Wartości macierzy dane są wzorem Gi,j = (vi,vj).

Spis treści

[edytuj] Postać

Macierz Grama ma postać:

\left[\begin{array}{c c c c c c}
(v_1, v_1)    & (v_1, v_2)    & (v_1, v_3)      & (v_1, v_4)     & \cdots  & (v_1, v_k) \\
(v_2, v_1)    & (v_2, v_2)    & (v_2, v_3)      & (v_2, v_4)     & \cdots  & (v_2, v_k) \\
(v_3, v_1)    & (v_3, v_2)    & (v_3, v_3)      & (v_3, v_4)     & \cdots  & (v_3, v_k) \\
\vdots  & \vdots  & \vdots  & \ddots & \ddots  & \vdots \\
(v_k, v_1)    & (v_k, v_2)    & (v_k, v_3)      & (v_k, v_4)     & \cdots  & (v_k, v_k) \\\end{array}\right]

Macierz ta jest stosowana do obliczania liniowej niezależności zbioru wektorów – zbiór wektorów jest liniowo niezależny wtedy i tylko wtedy, gdy wyznacznik jego macierzy Grama (wyznacznik Grama) jest niezerowy.

Nazwa macierzy pochodzi od duńskiego matematyka Jørgena Pedersena Grama.

[edytuj] Wyznacznik Grama

Wyznacznik Grama oznaczamy  G(v_1, v_2, \cdots , v_k)

Wartość wyznacznika Grama nie zależy od wyboru kolejności wektorów. Przestawienie dwóch wektorów w ciągu v_1, v_2, \cdots, v_k powoduje natchmiast prestawienie dwóch wierszy i dwóch kolumn w macierzy.

[edytuj] Zastosowanie – objętość

Wyznacznik Grama ma duże zastosowanie przy liczeniu objętości k-wymiarowych (objętość 1-dim to długość, 2-dim to powierzchnia)

Objętość równoległościanu rozpiętego przez wektory v_1, v_2, \cdots, v_k jest równa pierwiastkowi kwadratowemu z wyznacznika Grama

Wyznacznik Grama ma też zastosowanie przy obliczaniu k-objętości brył o bardziej skomplikowanym kształcie, np.: obliczanie długości krzywej

[edytuj] Przykład obliczania długości krzywej z zastosowaniem wyznacznika Grama

Dla przykładu policzymy długość okręgu o równaniu:

x2 + y2 = R2

okrąg sparametryzujemy przy pomocy t , które pochodzi z przedziału od 0 do . Stąd krzywą można opisać wektorem: (Rcos(t),Rsin(t)) a wektor styczny (obliczony jako pochodna) to: ( − Rsin(t),Rcos(t)) Wektory styczne rozpinają element objętości stąd: element objętości wynosi: dV = \sqrt{G(v_1, v_2, \cdots, v_k)} dt. W rozważanym przykładzie mamy jeden wektor styczny, a wyznacznik Grama jest równy długości tego wektora, czyli r. Aby obliczyć objętość liczymy:  \int\limits_0^{2 \pi} \sqrt{G(v_1, v_2, \cdots, v_k)} dt = \int\limits_0^{2 \pi} R dt =2 \pi R

Dokładnie analogicznie oblicza się powierzchnie, objętości brył, z tym że wykonywane jest całkowanie po większej liczbie zmiennych, a macierz Grama jest wyższego rzędu.

[edytuj] Sterowalność układu

Macierz Grama jest związana z liniowymi układami sterowania, która pozwala określić czy układ jest sterowalny.

[edytuj] Twierdzenie

Warunkiem koniecznym i wystarczającym sterowalności układu liniowego jest odwracalność macierzy Grama (dla pewnego tf > 0).

[edytuj] Wzory

Załóżmy, że układ jest niesterowalny, zatem będzie istniał taki niezerowy wektor v, że Gv = 0, a tym samym vTGv = 0. Po wprowadzeniu wektora vT oraz v "pod całkę" zauważamy, że zawartość całki można zapisać jako iloczyn wTw (w=B^Te^{-sA^T}v). Następnie zauważamy, że:

wTw = | | w | | 2.

Skoro G = 0, to | | w | | = 0, a tym samym w = 0. W ten sposób otrzymujemy (po dodatkowej transpozycji):

vTe sAB = 0.

Kolejne pochodne tego wzoru będą przedstawiały się (pomijając znak) jako:

\begin{matrix}v^TAe^{-sA}B=0 & v^TA^2e^{-sA}B=0 & ... &v^TA^{n-1}e^{-sA}B=0\end{matrix}, a po podstawieniu s = 0 otrzymamy ostatecznie iloczyn wektora vT z macierzą Kalmana.

[edytuj] Bibliografia

  1. Franciszek Leja: Rachunek różniczkowy i całkowy. Warszawa: PWN, 1976. 
  2. S. Zakrzewski: Algebra i Geometria. Warszawa. 

RMF FM: Lech Kaczyński zamówił czarter do Bruskeli
Kancelaria Prezydenta zamówiła lot czarterowy z Warszawy do Brukseli - dowiedziała się nieoficjalnie reporterka RMF FM Agnieszka Burzyńska. W razie braku porozumienia z premierem i braku zgody na wspólny lot, Lecha Kaczyńskiego zabierze na szczyt niewielki wynajęty samolot Embraer.
Jedna osoba ranna po eksplozji w zakładzie
Jedna osoba została ranna, a ponad 80 ewakuowano po tym jak w zakładzie przy ul. Technicznej w Łodzi doszło do eksplozji podczas ładowania akumulatorów - poinformował rzecznik straży pożarnej Artur Michalak.
Ćwiąkalski przekonuje: sądy 24-godzinne są za drogie
Minister sprawiedliwości Zbigniew Ćwiąkalski zapowiedział zmiany w działalności sądów 24 godzinnych.
70-lecie tragedii prawosławnych na Lubelszczyźnie
Dziś drugi dzień obchodów 70-ej rocznicy akcji burzenia cerkwi prawosławnych na Chełmszczyźnie i Południowym Podlasiu.
"Obecność prezydenta naruszyłaby dobry obyczaj"
Szef gabinetu politycznego premiera Sławomir Nowak powiedział że jednoczesna obecność premiera i prezydenta na szczycie w Brukseli naruszyłaby zarówno porządek konstytucyjny jak i dobry obyczaj.
Linki: Strona gwna