Mam trochę trudności z próbą zrozumienia artykułu. W pracy zastosowano metodę spektralną do obliczenia wartości własnej pochodzącej z układu sprzężonych ODE. Wypiszę teraz tylko jedno równanie, ponieważ wystarczy przejść do sedna moich pytań. Równanie to
Rozwijam jakiś większy kod do wykonywania obliczeń wartości własnych ogromnych rzadkich macierzy w kontekście fizyki obliczeniowej. Moje procedury sprawdzam na prostym oscylatorze harmonicznym w jednym wymiarze, ponieważ wartości własne są dobrze znane analitycznie. Robiąc to i porównując własne...
Mam dwa wykresy z prawie n ~ 100000 węzłów każdy. Na obu wykresach każdy węzeł jest podłączony dokładnie do 3 innych węzłów, więc macierz przylegania jest symetryczna i bardzo rzadka. Najtrudniejsze jest to, że potrzebuję wszystkich wartości własnych macierzy przylegania, ale nie wektorów...
Mam listę macierzy symetrycznych, które muszę sprawdzić pod kątem dodatniej półokreśloności (tzn. Ich wartości własne są nieujemne).LL{\cal L} Powyższy komentarz sugeruje, że można to zrobić, obliczając odpowiednie wartości własne i sprawdzając, czy nie są one ujemne (być może trzeba zająć się...
Załóżmy, że jest symetryczną, dodatnią określoną macierzą. A jest na tyle duże, że rozwiązanie A x = b jest drogie .A∈Rn×nA∈Rn×nA\in\mathbb{R}^{n\times n}AAAAx=bAx=bAx=b Czy istnieje iteracyjny algorytm znajdowania najmniejszej wartości własnej , który nie obejmuje odwracania A w każdej...
Mam zestaw danych, który powoli się zmienia i muszę śledzić wektory własne / wartości własne macierzy kowariancji. Używałem scipy.linalg.eigh, ale jest zbyt drogi i nie korzysta z faktu, że mam już rozkład, który jest tylko nieznacznie niepoprawny. Czy ktoś może zaproponować lepsze podejście do...
Każdy prawdziwy matryca można ograniczyć do rzeczywistego postaci Schur T = U T A U stosując prostopadłe similiary przekształcenie U . Tutaj macierz T ma postać quasi-trójkątną z 1 na 1 lub 2 na 2 blokami na głównej przekątnej. Każdy 1 o 1 odpowiada blok do rzeczywistej wartości własnej A i co 2 2...
Mam duży problem wartości własnej sześciennej: (A0+λA1+λ2A2+λ3A3)x=0.(A0+λA1+λ2A2+λ3A3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Mógłbym to rozwiązać, przechodząc na liniowy problem wartości własnych, ale spowodowałoby to,...