Pytania dotyczące algorytmicznych / obliczeniowych aspektów algebry liniowej, w tym rozwiązania układów liniowych, problemów z najmniejszymi kwadratami, problemów własnych i innych podobnych kwestii.
Czy ktoś mógłby polecić metodę dla następującego problemu najmniejszych kwadratów: znajdź który minimalizuje: , gdzie jest jednostką (obrót) matryca.N ∑ i = 0 ( R x i - b i ) 2 → minR ∈ R3 × 3R∈R3×3R \in \mathbb{R}^{3 \times 3}∑i = 0N.( R xja- bja)2)→ min∑i=0N(Rxi−bi)2→min\sum\limits_{i=0}^N (Rx_i - …
Rozważmy symetryczny dodatnio określona tridiagonal system liniowy , gdzie i . Biorąc pod uwagę trzy wskaźniki , jeśli przyjmiemy tylko rzędy równań ściśle między i hold, możemy wyeliminować zmienne pośrednie, aby uzyskać równanie w postaci gdzie . To równanie odnosi wartość do niezależnie od wpływu „zewnętrznego” (powiedzmy, jeśli wprowadzono ograniczenie …
Jak na przykład MATLAB oblicza SVD danej macierzy? Zakładam, że odpowiedź prawdopodobnie obejmuje obliczenie wektorów własnych i wartości własnych A*A'. Jeśli tak jest, chciałbym również wiedzieć, jak to oblicza?
Dano mi macierz , która jest symetryczna, odwracalna, dodatnia określona i gęsta. Muszę sprawdzić, czy gdzie J jest macierzą wszystkich jedynek.Q det ( Q ) = det ( 12 I - Q - J )12×1212×1212 \times 12QQQJdet ( Q ) = det ( 12 I- Q - J)( 1 )det(Q)=det(12I−Q−J)(1)\det(Q) …
Mam niejednorodny układ liniowy A x = bAx=b Ax=b gdzie jest rzeczywistą macierzą n × n przy n ≤ 4 . Nullspace z A jest gwarancją wymiaru zerowym więc równanie ma unikalną odwrotną x = A - 1 b . Ponieważ wynik wchodzi w prawą stronę ODE, którą zamierzam rozwiązać …
Twierdzenie o rozkładzie Cholesky'ego stwierdza zatem, że każda prawdziwa symetryczna macierz dodatnio określona ma rozkład Choleskiego gdzie jest dolną macierzą trójkątną.M = L L ⊤ L.MMMM=LL⊤M=LL⊤M= LL^\topL.LL Biorąc pod uwagę , wiemy już, że istnieją szybkie algorytmy do obliczania jego współczynnika Cholesky'ego .LM.MMLLL Załóżmy teraz, że otrzymałem prostokątną macierz , …
Rozwiązuję dla ogromnej rzadkiej dodatniej określonej macierzy za pomocą metody gradientu sprzężonego (CG). Czy można obliczyć wyznacznik podstawie informacji uzyskanych podczas rozwiązania?A AA x = bZAx=bAx=bZAZAAZAZAA
Na przykład nVidia ma CUBLAS, który obiecuje 7-14-krotne przyspieszenie. Naiwnie nie jest to teoretyczna przepustowość żadnej z kart GPU nVidii. Jakie są wyzwania związane z przyspieszaniem algebry liniowej na GPU i czy są już dostępne szybsze trasy algebry liniowej?
Załóżmy, że jest prawdziwą macierzą symetryczną i podano jej rozkład wartości własnej V Λ V T. Łatwo jest zobaczyć, co dzieje się z wartościami własnymi sumy A + c I, gdzie c jest stałą skalarną (patrz to pytanie ). Czy możemy wyciągnąć wnioski w ogólnym przypadku A + D, gdzie …
Załóżmy, że A jest ogólną macierzą rzadką i chcę obliczyć wartości własne. Nie wiem, jak wykryć wielokrotność wartości własnych. O ile mi wiadomo, w specjalnym przypadku znalezienia wielomianowych pierwiastków metodą macierzy towarzyszącej, możemy zastosować RRQR w celu wykrycia wielokrotności pierwiastków.
Biorąc pod uwagę system gdzie , przeczytałem, że w przypadku gdy iteracja Jacobiego jest używana jako solver, metoda nie zbiegnie się, jeśli ma wartość niezerową składnik zerowa przestrzeni . Jak zatem można formalnie stwierdzić, że pod warunkiem, że ma niezerowy składnik obejmujący zerową przestrzeń , metoda Jacobiego jest niespójna? Zastanawiam …
Robię diagonalizację Lanczosa dużej rzadkiej macierzy (~ 2 miliony elementów). Prawie wszystkie kroki w algorytmie Lanzcosa są wykonywane równolegle na GPU, z wyjątkiem diagonalizacji macierzy Lanczosa w celu sprawdzenia zbieżności. W tym celu korzystałem z algorytmu TQLI z receptur numerycznych. Czy istnieją metody znalezienia układu macierzystego macierzy tridiagonalnej, które są …
Czy istnieje szybszy sposób obliczenia błędów standardowych dla problemów z regresją liniową niż odwrócenie ? Tutaj zakładam, że mamy regresję:X′XX′XX'X y=Xβ+ε,y=Xβ+ε,y=X\beta+\varepsilon, gdzie jest macierzą n × k , ay jest wektorem n × 1 .XXXn×kn×kn\times kyyyn×1n×1n\times 1 Dla znalezienia najmniejszych kwadratów rozwiązania problemu jest niepraktyczne nic zrobić z , można …
W informatyce często spotykamy duże układy liniowe, które musimy rozwiązać za pomocą niektórych (skutecznych) środków, np. Metod bezpośrednich lub iteracyjnych. Jeśli skupimy się na tym drugim, jak możemy ustalić, że iteracyjna metoda rozwiązywania dużych układów liniowych jest w praktyce zbieżna? Oczywiste jest, że możemy przeprowadzić analizę prób i błędów (por. …
Używamy plików cookie i innych technologii śledzenia w celu poprawy komfortu przeglądania naszej witryny, aby wyświetlać spersonalizowane treści i ukierunkowane reklamy, analizować ruch w naszej witrynie, i zrozumieć, skąd pochodzą nasi goście.
Kontynuując, wyrażasz zgodę na korzystanie z plików cookie i innych technologii śledzenia oraz potwierdzasz, że masz co najmniej 16 lat lub zgodę rodzica lub opiekuna.