Kalkulator Postaci Normalnej Jordana

O Kalkulator Postaci Normalnej Jordana

Kalkulator Postaci Normalnej Jordana generuje postać kanoniczną Jordana J macierzy kwadratowej A wraz z odwracalną macierzą przejścia P spełniającą relację podobieństwa P⁻¹AP = J. W przeciwieństwie do diagonalizacji, która nie udaje się dla macierzy wadliwych, postać Jordana istnieje dla każdej macierzy kwadratowej nad ciałem algebraicznie domkniętym — zastępuje ona reprezentację diagonalną sekwencją klatek Jordana, z których każda jest macierzą prawie diagonalną z wartością własną na przekątnej i jedynkami na nadprzekątnej. Narzędzie to oblicza wszystko przy użyciu dokładnej arytmetyki wymiernej, więc wynikowe J i P są dowodliwie poprawne — nie występują tu zaokrąglenia zmiennoprzecinkowe.

Co to jest postać normalna Jordana?

Dla macierzy A o rozmiarze n × n nad liczbami zespolonymi, postać normalna Jordana J jest macierzą blokowo-diagonalną

gdzie każda klatka Jordana J_k(λ) jest macierzą k × k z wartością λ na przekątnej, jedynkami na nadprzekątnej i zerami w pozostałych miejscach:

Wartości własne λ_i mogą się powtarzać w różnych klatkach; kluczowy jest wzór rozmiarów klatek, który jest kompletnym niezmiennikiem podobieństwa macierzy A.

Dlaczego potrzebujemy postaci Jordana, skoro mamy diagonalizację?

Nie każda macierz kwadratowa jest diagonalizowalna. Macierz nie jest diagonalizowalna, gdy pewna wartość własna ma mniej niezależnych wektorów własnych niż wynosi jej krotność algebraiczna — mówimy wtedy, że macierz jest wadliwa. Postać Jordana naprawia tę lukę poprzez wprowadzenie uogólnionych wektorów własnych, dając postać kanoniczną, która działa dla każdej macierzy.

Kluczowe pojęcia

Krotność algebraiczna vs geometryczna

Krotność algebraiczna wartości własnej λ to krotność λ jako pierwiastka wielomianu charakterystycznego p_A(λ) = det(λI − A). Krotność geometryczna to wymiar przestrzeni własnej, czyli dim ker(A − λI). Liczba klatek Jordana powiązanych z λ jest równa jej krotności geometrycznej, a całkowity rozmiar tych klatek jest równy jej krotności algebraicznej.

Uogólnione wektory własne i łańcuchy

Wektor v jest uogólnionym wektorem własnym rzędu k dla wartości własnej λ, jeśli (A − λI)^kv = 0, ale (A − λI)^k−1v ≠ 0. Zastosowanie operatora N = (A − λI) do uogólnionego wektora własnego rzędu k generuje wektor rzędu k−1, co pozwala uzyskać łańcuch Jordana:

Umieszczenie łańcucha w kolejności v₁, v₂, …, v_k jako kolumn macierzy P tworzy klatkę Jordana o rozmiarze k w odpowiednich wierszach/kolumnach macierzy J.

Drabina jąder i zliczanie klatek

Dla każdej wartości własnej λ definiujemy rosnący ciąg d_k = dim ker((A − λI)^k). Ciąg ten jest niemalejący i stabilizuje się na poziomie krotności algebraicznej λ. Liczba klatek Jordana poszczególnych rozmiarów jest wyznaczana z tej drabiny:

Jest to zliczanie oparte na diagramie Younga i jest ono dokładne — nie wymaga zgadywania. Kalkulator wyświetla tę drabinę dla każdej wartości własnej, aby umożliwić śledzenie dekompozycji krok po kroku.

Wielomian minimalny

Wielomian minimalny m_A(λ) to wielomian unormowany najniższego stopnia, który spełnia m_A(A) = 0. Po uzyskaniu postaci Jordana odczytanie go jest trywialne:

Macierz jest diagonalizowalna wtedy i tylko wtedy, gdy jej wielomian minimalny nie ma pierwiastków wielokrotnych, tzn. każda klatka Jordana ma rozmiar 1.

Jak działa ten kalkulator

1. Analiza macierzy — akceptowane są liczby całkowite, ułamki (np. 1/2) oraz liczby dziesiętne, które są konwertowane na dokładne liczby wymierne (fractions.Fraction).
2. Obliczanie wielomianu charakterystycznego przy użyciu algorytmu Faddieewa–Leverriera, który unika symbolicznego rozwinięcia wyznacznika i działa w czasie O(n⁴) przy użyciu dokładnej arytmetyki.
3. Znajdowanie wymiernych wartości własnych poprzez twierdzenie o pierwiastkach wymiernych — każdy wymierny pierwiastek p/q wielomianu o współczynnikach całkowitych spełnia p ∣ wyraz wolny oraz q ∣ współczynnik przy najwyższej potędze. Każdy znaleziony pierwiastek jest dzielony, a poszukiwanie jest powtarzane.
4. Budowanie drabiny jąder dla każdej wartości własnej λ poprzez obliczanie dim ker((A − λI)^k) przy użyciu wymiernej postaci schodkowej (RREF), aż ciąg ustabilizuje się na krotności algebraicznej.
5. Wybór wektorów startowych łańcucha od największego jądra w dół do najmniejszego, rozszerzając bazę zawsze, gdy wymagana jest nowa klatka Jordana. Każdy wektor startowy jest następnie wielokrotnie mnożony przez (A − λI), aby uzyskać wektory łańcucha.
6. Składanie J i P poprzez grupowanie łańcuchów dla każdej wartości własnej (najpierw klatki o największym rozmiarze), umieszczając wektory łańcucha jako kolumny P i wypełniając J wartościami własnymi oraz jedynkami nad przekątną.
7. Weryfikacja dokładna faktu, że P⁻¹ A P = J przy użyciu arytmetyki liczb całkowitych — wynik jest gwarantowany, ponieważ wszystkie obliczenia pośrednie są wymierne.

Przykład krok po kroku

Rozważmy wadliwą macierz 3 × 3

• Wielomian charakterystyczny: \(p_A(\lambda) = (\lambda - 5)^3\). Pojedyncza wartość własna λ = 5 o krotności algebraicznej 3.
• Drabina jądra dla λ = 5: \(d_1 = 1\), \(d_2 = 2\), \(d_3 = 3\). Przyrosty wynoszą 1, 1, 1 → jedna klatka Jordana o rozmiarze 3.
• Postać Jordana: \(J = \begin{bmatrix} 5 & 1 & 0 \\ 0 & 5 & 1 \\ 0 & 0 & 5 \end{bmatrix}\), z krotnością geometryczną 1 i indeksem 3.
• Wielomian minimalny: \(m_A(\lambda) = (\lambda - 5)^3\) — taki sam jak wielomian charakterystyczny, ponieważ istnieje tylko jedna klatka Jordana.

Zastosowania Postaci Normalnej Jordana

• Eksponenty macierzy i liniowe równania różniczkowe zwyczajne (ODE) — dla układu o stałych współczynnikach x′ = Ax, rozwiązaniem w postaci zamkniętej jest \(e^{tA}x_0\), a \(e^{tA}\) jest proste do obliczenia, gdy A jest zapisana w postaci Jordana.
• Potęgi macierzy — \(A^k = P J^k P^{-1}\), a klatki Jordana posiadają jawne wzory na ich potęgi.
• Rachunek funkcyjny — \(f(A) = P f(J) P^{-1}\) uogólnia się na dowolną funkcję analityczną f, o ile f jest określona w otoczeniu widma.
• Teoria sterowania — stabilność układów liniowych zależy od wartości własnych oraz rozmiarów klatek Jordana (przypadki graniczne wymagają sprawdzenia największej klatki dla granicznej wartości własnej).
• Klasyfikacja operatorów liniowych — dwie macierze są podobne wtedy i tylko wtedy, gdy mają tę samą postać Jordana, więc postać ta jest kompletnym niezmiennikiem.

Często zadawane pytania

Co to jest postać normalna Jordana macierzy?

Postać normalna Jordana (zwana również postacią kanoniczną Jordana) to macierz prawie diagonalna J podobna do oryginalnej macierzy A, co oznacza, że istnieje odwracalna macierz P taka, że P⁻¹AP = J. Przekątna J zawiera wartości własne A, a tuż nad przekątną znajdują się jedynki, które pojawiają się wewnątrz klatek Jordana, gdy macierz A nie jest diagonalizowalna. Każda macierz kwadratowa nad liczbami zespolonymi posiada postać normalną Jordana, unikalną z dokładnością do kolejności klatek.

Kiedy macierz nie jest diagonalizowalna?

Macierz nie jest diagonalizowalna, gdy co najmniej jedna wartość własna ma mniej liniowo niezależnych wektorów własnych niż wynosi jej krotność algebraiczna — luka ta jest wypełniana przez klatki Jordana o rozmiarze 2 lub większym. Równoważnie, macierz nie jest diagonalizowalna, gdy jej wielomian minimalny ma pierwiastek wielokrotny. Takie macierze nazywamy wadliwymi.

Jak definiuje się uogólnione wektory własne?

Uogólniony wektor własny rzędu k dla wartości własnej λ to niezerowy wektor v taki, że (A − λI)^kv = 0, ale (A − λI)^k−1v jest niezerowy. Zastosowanie (A − λI) do uogólnionego wektora własnego rzędu k daje wektor rzędu k−1, tworząc łańcuch. Łańcuchy te tworzą kolumny macierzy przejścia P w dekompozycji Jordana.

Jaka jest różnica między krotnością algebraiczną a geometryczną?

Krotność algebraiczna wartości własnej λ to liczba określająca, ile razy występuje ona jako pierwiastek wielomianu charakterystycznego. Krotność geometryczna to wymiar jej przestrzeni własnej — liczba liniowo niezależnych wektorów własnych. Krotność geometryczna jest równa liczbie klatek Jordana dla danej λ, podczas gdy krotność algebraiczna jest równa całkowitemu rozmiarowi wszystkich tych klatek. Równe krotności oznaczają, że wartość własna generuje tylko klatki o rozmiarze 1.

Jak ten kalkulator znajduje rozmiary klatek Jordana?

Dla każdej wartości własnej λ kalkulator oblicza wymiary d_k = dim ker((A − λI)^k) dla k = 1, 2, … aż sekwencja ustabilizuje się na krotności algebraicznej. Liczba klatek Jordana o rozmiarze co najmniej k jest równa d_k − d_k−1. Odejmując kolejne wyrazy, otrzymujemy dokładną liczbę klatek każdego rozmiaru. Te obliczenia oparte na diagramach Younga są dokładne i wykorzystują arytmetykę wymierną.

Czy kalkulator obsługuje macierze z niewymiernymi lub zespolonymi wartościami własnymi?

Kalkulator używa dokładnej arytmetyki wymiernej, co wymaga, aby wartości własne były liczbami wymiernymi. Gdy wielomian charakterystyczny ma czynniki, które nie rozkładają się nad liczbami wymiernymi, narzędzie pokazuje przybliżone numerycznie zespolone wartości własne dla pozostałego czynnika, ale nie generuje pełnej postaci Jordana, ponieważ dokładna arytmetyka jest niezbędna do prawidłowego określenia rozmiarów klatek. Przeskaluj lub zmodyfikuj macierz tak, aby wszystkie wartości własne były wymierne, aby uzyskać pełną dekompozycję Jordana.

Co to jest wielomian minimalny i jak jest tutaj obliczany?

Wielomian minimalny m(λ) to wielomian unormowany najniższego stopnia, który zeruje macierz A, czyli m(A) = 0. Jest on równy iloczynowi czynników (λ − λ_i)^indeks_i po wszystkich różnych wartościach własnych, gdzie indeks jest rozmiarem największej klatki Jordana dla wartości własnej λ_i. Kalkulator odczytuje indeks bezpośrednio z obliczonej struktury blokowej, więc wielomian minimalny jest darmowym produktem ubocznym dekompozycji Jordana.

Dalsza lektura

• Postać Jordana — Wikipedia
• Uogólniony wektor własny — Wikipedia (EN)
• Wielomian minimalny — Wikipedia
• Algorytm Faddieewa–Leverriera — Wikipedia (EN)