Ostatnio porównywałem różne nieliniowe solwery z scipy i byłem pod szczególnym wrażeniem przykładu Newtona-Kryłowa w książce kucharskiej Scipy, w którym rozwiązują równanie różniczkowe drugiego rzędu z nieliniowym terminem reakcji w około 20 wierszach kodu. Zmodyfikowałem przykładowy kod, aby...
Zaimplementowałem solver z Eulerem wstecznym w Pythonie 3 (używając numpy). Dla własnej wygody i jako ćwiczenie napisałem również małą funkcję, która oblicza przybliżoną różnicę skończoną różnicy gradientu, aby nie zawsze musiałem określać analitycznie jakobian (jeśli to w ogóle...
Próbuję dowiedzieć się, czy istnieje szybszy sposób obliczenia wszystkich wartości własnych i wektorów własnych bardzo dużej i rzadkiej macierzy przylegania niż przy użyciu scipy.sparse.linalg.eigsh O ile mi wiadomo, metody te wykorzystują tylko rzadkość i atrybuty symetrii macierzy. Macierz...
Pisałem przybornik systemu sterowania od zera i wyłącznie w Python3 (bezwstydna wtyczka:) harold. Na podstawie moich wcześniejszych badań zawsze narzekałem na solver Riccati care.mz powodów technicznych / nieistotnych. Dlatego piszę własny zestaw procedur. Jednej rzeczy, której nie mogę znaleźć,...