Kalkulator rozkładu LU macierzy

Witaj w kalkulatorze rozkładu LU macierzy, kompleksowym narzędziu do algebry liniowej, które rozkłada dowolną macierz kwadratową na iloczyn macierzy dolnotrójkątnej (L) i górnotrójkątnej (U) przy użyciu eliminacji Gaussa z częściowym wyborem elementu podstawowego. Otrzymasz szczegółową eliminację krok po kroku, interaktywną animację rozkładu oraz automatyczną weryfikację. Idealne dla studentów, inżynierów i każdego, kto pracuje z układami równań liniowych.

Co to jest rozkład LU?

Rozkład LU (zwany również faktoryzacją LU) wyraża macierz kwadratową \(A\) jako iloczyn dwóch macierzy trójkątnych:

Gdzie:

• L (dolnotrójkątna): posiada 1 na przekątnej oraz niezerowe elementy tylko poniżej przekątnej. Elementy te są mnożnikami użytymi podczas eliminacji Gaussa.
• U (górnotrójkątna): posiada niezerowe elementy tylko na przekątnej i powyżej niej. Jest to postać schodkowa wierszowa macierzy.

Gdy stosowany jest częściowy wybór elementu podstawowego (aby uniknąć zerowych elementów głównych i poprawić stabilność numeryczną), faktoryzacja przyjmuje postać:

Gdzie \(P\) to macierz permutacji, która rejestruje zamiany wierszy wykonane podczas eliminacji.

Jak korzystać z tego kalkulatora

1. Wprowadź macierz: Wpisz macierz kwadratową, umieszczając wiersze w oddzielnych liniach lub oddzielając je średnikami. Elementy mogą być oddzielone spacjami, przecinkami lub tabulatorami. Obsługuje rozmiary do 8×8.
2. Ustaw precyzję dziesiętną: Wybierz, ile miejsc po przecinku (2-10) ma być wyświetlanych w wynikach.
3. Kliknij Rozłóż: Kalkulator wykona faktoryzację LU z częściowym wyborem elementu podstawowego i pokaże wyniki.
4. Przejrzyj wyniki: Zbadaj macierze L, U i P, animowany rozkład oraz szczegóły eliminacji krok po kroku.

Rozwiązywanie układów liniowych za pomocą rozkładu LU

Rozkład LU jest szczególnie przydatny przy rozwiązywaniu układów równań liniowych \(Ax = b\). Po uzyskaniu \(PA = LU\), rozwiązanie staje się procesem dwuetapowym:

Krok 1: Podstawienie w przód

Rozwiąż \(Ly = Pb\) dla \(y\). Ponieważ \(L\) jest dolnotrójkątna, jest to proste — zacznij od górnego równania i idź w dół:

Krok 2: Podstawienie wstecz

Rozwiąż \(Ux = y\) dla \(x\). Ponieważ \(U\) jest górnotrójkątna, zacznij od dolnego równania i idź w górę:

Obliczanie wyznacznika

Wyznacznik macierzy \(A\) można wydajnie obliczyć z czynników LU:

Gdzie \(s\) to liczba zamian wierszy (pivotów), a \(U_{ii}\) to elementy na przekątnej macierzy \(U\). Ponieważ \(\det(L) = 1\) (wszystkie elementy na przekątnej to 1) oraz \(\det(P) = (-1)^s\), wzór wynika z zależności \(\det(P)\det(A) = \det(L)\det(U)\).

Dlaczego częściowy wybór elementu podstawowego?

Bez zamiany wierszy, rozkład LU zawodzi, jeśli jakikolwiek element główny wynosi zero. Nawet gdy elementy główne są niezerowe, ale małe, wynik obliczeń może być obarczony poważnymi błędami numerycznymi. Częściowy wybór elementu podstawowego wybiera największy dostępny element główny w każdej kolumnie, co:

• Zapobiega dzieleniu przez zero
• Minimalizuje wzrost błędów zaokrągleń
• Gwarantuje, że mnożniki w L spełniają warunek \(|L_{ij}| \leq 1\)
• Zapewnia, że każdą nieosobliwą macierz można rozłożyć

Zastosowania rozkładu LU

LU vs inne rozkłady

• LU vs QR: LU jest szybszy (\(O(\frac{2}{3}n^3)\) vs \(O(\frac{4}{3}n^3)\)), ale mniej stabilny numerycznie. QR jest preferowany w problemach najmniejszych kwadratów.
• LU vs Cholesky: Rozkład Choleskiego (\(A = LL^T\)) działa tylko dla macierzy symetrycznych dodatnio określonych, ale jest dwukrotnie szybszy i bardziej stabilny niż ogólny rozkład LU.
• LU vs eliminacja Gaussa: Rozkład LU jest eliminacją Gaussa, ale postać sfaktoryzowana L i U może być ponownie użyta do wydajnego rozwiązywania wielu wektorów prawych stron.

Często zadawane pytania

Co to jest rozkład LU?

Rozkład LU (zwany również faktoryzacją LU) to metoda polegająca na rozkładzie kwadratowej macierzy A na iloczyn macierzy dolnotrójkątnej L i górnotrójkątnej U, tak aby A = LU (lub PA = LU z częściowym wyborem elementu podstawowego). Macierz L ma jedynki na przekątnej i przechowuje mnożniki eliminacji, podczas gdy U jest wynikiem eliminacji Gaussa.

Dlaczego w rozkładzie LU potrzebny jest częściowy wybór elementu podstawowego?

Częściowy wybór elementu podstawowego (pivoting) polega na zamianie wierszy w celu umieszczenia największej wartości bezwzględnej na pozycji głównej. Zapobiega to dzieleniu przez zero, gdy element główny wynosi zero, oraz zmniejsza błędy numeryczne spowodowane dzieleniem przez bardzo małe liczby. Przy zastosowaniu tej metody faktoryzacja przyjmuje postać PA = LU, gdzie P jest macierz permutacji rejestrującą zamiany wierszy.

Jakie są zastosowania rozkładu LU?

Rozkład LU jest używany do wydajnego rozwiązywania układów równań liniowych (Ax = b), obliczania wyznaczników macierzy, znajdowania macierzy odwrotnych oraz analizy stabilności numerycznej. Jest szczególnie wydajny przy rozwiązywaniu wielu układów z tą samą macierzą współczynników, ale różnymi wektorami prawych stron, ponieważ faktoryzację wykonuje się tylko raz.

Jak rozwiązać Ax = b za pomocą rozkładu LU?

Po obliczeniu PA = LU, rozwiązanie Ax = b odbywa się w dwóch krokach: najpierw rozwiązuje się Ly = Pb stosując podstawienie w przód (łatwe, bo L jest dolnotrójkątna), a następnie rozwiązuje się Ux = y stosując podstawienie wstecz (łatwe, bo U jest górnotrójkątna). Ten proces jest znacznie szybszy niż eliminacja Gaussa przy rozwiązywaniu wielu układów.

Czy każdą macierz kwadratową można rozłożyć metodą LU?

Nie każda macierz kwadratowa posiada rozkład LU bez zastosowania zamiany wierszy. Macierz posiada faktoryzację LU wtedy i tylko wtedy, gdy wszystkie jej główne minory wiodące są niezerowe. Jednakże, stosując częściowy wybór elementu podstawowego (PA = LU), każdą nieosobliwą macierz kwadratową można rozłożyć. Macierze osobliwe mogą zawieść, jeśli napotkany zostanie zerowy element główny.

Dodatkowe zasoby

• Rozkład LU - Wikipedia
• Macierz trójkątna - Wikipedia
• Eliminacja Gaussa - Wikipedia